Generalized Linear Covariance Analysis
NASA Technical Reports Server (NTRS)
Carpenter, James R.; Markley, F. Landis
2014-01-01
This talk presents a comprehensive approach to filter modeling for generalized covariance analysis of both batch least-squares and sequential estimators. We review and extend in two directions the results of prior work that allowed for partitioning of the state space into solve-for'' and consider'' parameters, accounted for differences between the formal values and the true values of the measurement noise, process noise, and textita priori solve-for and consider covariances, and explicitly partitioned the errors into subspaces containing only the influence of the measurement noise, process noise, and solve-for and consider covariances. In this work, we explicitly add sensitivity analysis to this prior work, and relax an implicit assumption that the batch estimator's epoch time occurs prior to the definitive span. We also apply the method to an integrated orbit and attitude problem, in which gyro and accelerometer errors, though not estimated, influence the orbit determination performance. We illustrate our results using two graphical presentations, which we call the variance sandpile'' and the sensitivity mosaic,'' and we compare the linear covariance results to confidence intervals associated with ensemble statistics from a Monte Carlo analysis.
Generalized Linear Covariance Analysis
NASA Technical Reports Server (NTRS)
Carpenter, J. Russell; Markley, F. Landis
2008-01-01
We review and extend in two directions the results of prior work on generalized covariance analysis methods. This prior work allowed for partitioning of the state space into "solve-for" and "consider" parameters, allowed for differences between the formal values and the true values of the measurement noise, process noise, and a priori solve-for and consider covariances, and explicitly partitioned the errors into subspaces containing only the influence of the measurement noise, process noise, and a priori solve-for and consider covariances. In this work, we explicitly add sensitivity analysis to this prior work, and relax an implicit assumption that the batch estimator s anchor time occurs prior to the definitive span. We also apply the method to an integrated orbit and attitude problem, in which gyro and accelerometer errors, though not estimated, influence the orbit determination performance. We illustrate our results using two graphical presentations, which we call the "variance sandpile" and the "sensitivity mosaic," and we compare the linear covariance results to confidence intervals associated with ensemble statistics from a Monte Carlo analysis.
Linear Algebraic Method for Non-Linear Map Analysis
Yu,L.; Nash, B.
2009-05-04
We present a newly developed method to analyze some non-linear dynamics problems such as the Henon map using a matrix analysis method from linear algebra. Choosing the Henon map as an example, we analyze the spectral structure, the tune-amplitude dependence, the variation of tune and amplitude during the particle motion, etc., using the method of Jordan decomposition which is widely used in conventional linear algebra.
NASA Technical Reports Server (NTRS)
Redding, David C.; Breckenridge, William G.
1990-01-01
A new, coordinate-free version of the exact ray-trace equations for optical systems consisting of conic reflecting, refracting and reference surfaces is presented. These equations are differentiated to obtain closed-form optical sensitivity dyadics. For computation, the sensitivities are evaluated in a single global coordinate frame and combined in linearized ray-trace matrix difference equations that propagate the rays and the sensitivities from element to element. One purpose of this analysis is to create optical models that can be directly integrated with models of the instrument structure and control systems for dynamic simulation.
Cluster Analysis by Linear Contrasts.
ERIC Educational Resources Information Center
Shafto, Michael
The purpose of this paper is to suggest a technique of cluster analysis which is similar in aim to the Interactive Intercolumnar Correlation Analysis (IICA), though different in detail. Two methods are proposed for extracting a single bipolar factor (a "contrast compenent") directly from the initial similarities matrix. The advantages of this…
Linear covariance analysis for gimbaled pointing systems
NASA Astrophysics Data System (ADS)
Christensen, Randall S.
Linear covariance analysis has been utilized in a wide variety of applications. Historically, the theory has made significant contributions to navigation system design and analysis. More recently, the theory has been extended to capture the combined effect of navigation errors and closed-loop control on the performance of the system. These advancements have made possible rapid analysis and comprehensive trade studies of complicated systems ranging from autonomous rendezvous to vehicle ascent trajectory analysis. Comprehensive trade studies are also needed in the area of gimbaled pointing systems where the information needs are different from previous applications. It is therefore the objective of this research to extend the capabilities of linear covariance theory to analyze the closed-loop navigation and control of a gimbaled pointing system. The extensions developed in this research include modifying the linear covariance equations to accommodate a wider variety of controllers. This enables the analysis of controllers common to gimbaled pointing systems, with internal states and associated dynamics as well as actuator command filtering and auxiliary controller measurements. The second extension is the extraction of power spectral density estimates from information available in linear covariance analysis. This information is especially important to gimbaled pointing systems where not just the variance but also the spectrum of the pointing error impacts the performance. The extended theory is applied to a model of a gimbaled pointing system which includes both flexible and rigid body elements as well as input disturbances, sensor errors, and actuator errors. The results of the analysis are validated by direct comparison to a Monte Carlo-based analysis approach. Once the developed linear covariance theory is validated, analysis techniques that are often prohibitory with Monte Carlo analysis are used to gain further insight into the system. These include the creation
Reachability analysis of rational eigenvalue linear systems
NASA Astrophysics Data System (ADS)
Xu, Ming; Chen, Liangyu; Zeng, Zhenbing; Li, Zhi-bin
2010-12-01
One of the key problems in the safety analysis of control systems is the exact computation of reachable state spaces for continuous-time systems. Issues related to the controllability and observability of these systems are well-studied in systems theory. However, there are not many results on reachability, even for general linear systems. In this study, we present a large class of linear systems with decidable reachable state spaces. This is approached by reducing the reachability analysis to real root isolation of exponential polynomials. Furthermore, we have implemented this method in a Maple package based on symbolic computation and applied to several examples successfully.
Thermal analysis of a linear infrared lamp
Nakos, J.T.
1982-01-01
A theoretical and experimental analysis of an infrared lamp is presented based on radiant heat transfer theory. The analysis is performed on a specific type of linear lamp which has a coiled tungsten filament surrounded by a fused quartz envelope. The purpose of the study was to model the lamp thermally, not electrically, to arrive at a better understanding of the operation of the lamp.
Numerical analysis method for linear induction machines.
NASA Technical Reports Server (NTRS)
Elliott, D. G.
1972-01-01
A numerical analysis method has been developed for linear induction machines such as liquid metal MHD pumps and generators and linear motors. Arbitrary phase currents or voltages can be specified and the moving conductor can have arbitrary velocity and conductivity variations from point to point. The moving conductor is divided into a mesh and coefficients are calculated for the voltage induced at each mesh point by unit current at every other mesh point. Combining the coefficients with the mesh resistances yields a set of simultaneous equations which are solved for the unknown currents.
Linear digital imaging system fidelity analysis
NASA Technical Reports Server (NTRS)
Park, Stephen K.
1989-01-01
The combined effects of imaging gathering, sampling and reconstruction are analyzed in terms of image fidelity. The analysis is based upon a standard end-to-end linear system model which is sufficiently general so that the results apply to most line-scan and sensor-array imaging systems. Shift-variant sampling effects are accounted for with an expected value analysis based upon the use of a fixed deterministic input scene which is randomly shifted (mathematically) relative to the sampling grid. This random sample-scene phase approach has been used successfully by the author and associates in several previous related papers.
Visual Analysis and Dissemination of Scientific Literature Collections with SurVis.
Beck, Fabian; Koch, Sebastian; Weiskopf, Daniel
2016-01-01
Bibliographic data such as collections of scientific articles and citation networks have been studied extensively in information visualization and visual analytics research. Powerful systems have been built to support various types of bibliographic analysis, but they require some training and cannot be used to disseminate the insights gained. In contrast, we focused on developing a more accessible visual analytics system, called SurVis, that is ready to disseminate a carefully surveyed literature collection. The authors of a survey may use our Web-based system to structure and analyze their literature database. Later, readers of the survey can obtain an overview, quickly retrieve specific publications, and reproduce or extend the original bibliographic analysis. Our system employs a set of selectors that enable users to filter and browse the literature collection as well as to control interactive visualizations. The versatile selector concept includes selectors for textual search, filtering by keywords and meta-information, selection and clustering of similar publications, and following citation links. Agreement to the selector is represented by word-sized sparkline visualizations seamlessly integrated into the user interface. Based on an analysis of the analytical reasoning process, we derived requirements for the system. We developed the system in a formative way involving other researchers writing literature surveys. A questionnaire study with 14 visual analytics experts confirms that SurVis meets the initially formulated requirements. PMID:26529698
Linear Covariance Analysis and Epoch State Estimators
NASA Technical Reports Server (NTRS)
Markley, F. Landis; Carpenter, J. Russell
2012-01-01
This paper extends in two directions the results of prior work on generalized linear covariance analysis of both batch least-squares and sequential estimators. The first is an improved treatment of process noise in the batch, or epoch state, estimator with an epoch time that may be later than some or all of the measurements in the batch. The second is to account for process noise in specifying the gains in the epoch state estimator. We establish the conditions under which the latter estimator is equivalent to the Kalman filter.
Linear Covariance Analysis and Epoch State Estimators
NASA Technical Reports Server (NTRS)
Markley, F. Landis; Carpenter, J. Russell
2014-01-01
This paper extends in two directions the results of prior work on generalized linear covariance analysis of both batch least-squares and sequential estimators. The first is an improved treatment of process noise in the batch, or epoch state, estimator with an epoch time that may be later than some or all of the measurements in the batch. The second is to account for process noise in specifying the gains in the epoch state estimator. We establish the conditions under which the latter estimator is equivalent to the Kalman filter.
Runtime Analysis of Linear Temporal Logic Specifications
NASA Technical Reports Server (NTRS)
Giannakopoulou, Dimitra; Havelund, Klaus
2001-01-01
This report presents an approach to checking a running program against its Linear Temporal Logic (LTL) specifications. LTL is a widely used logic for expressing properties of programs viewed as sets of executions. Our approach consists of translating LTL formulae to finite-state automata, which are used as observers of the program behavior. The translation algorithm we propose modifies standard LTL to B chi automata conversion techniques to generate automata that check finite program traces. The algorithm has been implemented in a tool, which has been integrated with the generic JPaX framework for runtime analysis of Java programs.
A Navigational Analysis of Linear and Non-Linear Hypermedia Interfaces.
ERIC Educational Resources Information Center
Hall, Richard H.; Balestra, Joel; Davis, Miles
The purpose of this experiment was to assess the effectiveness of a comprehensive model for the analysis of hypermap navigation patterns through a comparison of navigation patterns associated with a traditional linear interface versus a non-linear "hypermap" interface. Twenty-six general psychology university students studied material on bipolar…
Double Linear Damage Rule for Fatigue Analysis
NASA Technical Reports Server (NTRS)
Halford, G.; Manson, S.
1985-01-01
Double Linear Damage Rule (DLDR) method for use by structural designers to determine fatigue-crack-initiation life when structure subjected to unsteady, variable-amplitude cyclic loadings. Method calculates in advance of service how many loading cycles imposed on structural component before macroscopic crack initiates. Approach eventually used in design of high performance systems and incorporated into design handbooks and codes.
Advanced analysis technique for the evaluation of linear alternators and linear motors
NASA Astrophysics Data System (ADS)
Holliday, Jeffrey C.
1995-12-01
A method for the mathematical analysis of linear alternator and linear motor devices and designs is described, and an example of its use is included. The technique seeks to surpass other methods of analysis by including more rigorous treatment of phenomena normally omitted or coarsely approximated such as eddy braking, non-linear material properties, and power losses generated within structures surrounding the device. The technique is broadly applicable to linear alternators and linear motors involving iron yoke structures and moving permanent magnets. The technique involves the application of Amperian current equivalents to the modeling of the moving permanent magnet components within a finite element formulation. The resulting steady state and transient mode field solutions can simultaneously account for the moving and static field sources within and around the device.
Advanced analysis technique for the evaluation of linear alternators and linear motors
NASA Technical Reports Server (NTRS)
Holliday, Jeffrey C.
1995-01-01
A method for the mathematical analysis of linear alternator and linear motor devices and designs is described, and an example of its use is included. The technique seeks to surpass other methods of analysis by including more rigorous treatment of phenomena normally omitted or coarsely approximated such as eddy braking, non-linear material properties, and power losses generated within structures surrounding the device. The technique is broadly applicable to linear alternators and linear motors involving iron yoke structures and moving permanent magnets. The technique involves the application of Amperian current equivalents to the modeling of the moving permanent magnet components within a finite element formulation. The resulting steady state and transient mode field solutions can simultaneously account for the moving and static field sources within and around the device.
Employment of CB models for non-linear dynamic analysis
NASA Technical Reports Server (NTRS)
Klein, M. R. M.; Deloo, P.; Fournier-Sicre, A.
1990-01-01
The non-linear dynamic analysis of large structures is always very time, effort and CPU consuming. Whenever possible the reduction of the size of the mathematical model involved is of main importance to speed up the computational procedures. Such reduction can be performed for the part of the structure which perform linearly. Most of the time, the classical Guyan reduction process is used. For non-linear dynamic process where the non-linearity is present at interfaces between different structures, Craig-Bampton models can provide a very rich information, and allow easy selection of the relevant modes with respect to the phenomenon driving the non-linearity. The paper presents the employment of Craig-Bampton models combined with Newmark direct integration for solving non-linear friction problems appearing at the interface between the Hubble Space Telescope and its solar arrays during in-orbit maneuvers. Theory, implementation in the FEM code ASKA, and practical results are shown.
Bounded Linear Stability Margin Analysis of Nonlinear Hybrid Adaptive Control
NASA Technical Reports Server (NTRS)
Nguyen, Nhan T.; Boskovic, Jovan D.
2008-01-01
This paper presents a bounded linear stability analysis for a hybrid adaptive control that blends both direct and indirect adaptive control. Stability and convergence of nonlinear adaptive control are analyzed using an approximate linear equivalent system. A stability margin analysis shows that a large adaptive gain can lead to a reduced phase margin. This method can enable metrics-driven adaptive control whereby the adaptive gain is adjusted to meet stability margin requirements.
A linearized Euler analysis of unsteady transonic flows in turbomachinery
Hall, K.C.; Clark, W.S.; Lorence, C.B. . Dept. of Mechanical Engineering and Materials Science)
1994-07-01
A computational method for efficiently predicting unsteady transonic flows in two- and three-dimensional cascades is presented. The unsteady flow is modeled using a linearized Euler analysis whereby the unsteady flow field is decomposed into a nonlinear mean flow plus a linear harmonically varying unsteady flow. The equations that govern the perturbation flow, the linearized Euler equations, are linear variable coefficient equations. For transonic flows containing shocks, shock capturing is used to model the shock impulse (the unsteady load due to the harmonic motion of the shock). A conservative Lax-Wendroff scheme is used to obtain a set of linearized finite volume equations that describe the harmonic small disturbance behavior of the flow. Conditions under which such a discretization will correctly predict the shock impulse are investigated. Computational results are presented that demonstrate the accuracy and efficiency of the present method as well as the essential role of unsteady shock impulse loads on the flutter stability of fans.
Virtual Estimator for Piecewise Linear Systems Based on Observability Analysis
Morales-Morales, Cornelio; Adam-Medina, Manuel; Cervantes, Ilse; Vela-Valdés and, Luis G.; García Beltrán, Carlos Daniel
2013-01-01
This article proposes a virtual sensor for piecewise linear systems based on observability analysis that is in function of a commutation law related with the system's outpu. This virtual sensor is also known as a state estimator. Besides, it presents a detector of active mode when the commutation sequences of each linear subsystem are arbitrary and unknown. For the previous, this article proposes a set of virtual estimators that discern the commutation paths of the system and allow estimating their output. In this work a methodology in order to test the observability for piecewise linear systems with discrete time is proposed. An academic example is presented to show the obtained results. PMID:23447007
A failsafe analysis using NASTRAN's piecewise linear analysis and a nine node linear crack element
NASA Technical Reports Server (NTRS)
Wilkinson, R. F.; Kelley, J. W.
1975-01-01
A two-dimensional crack element was implemented into NASTRAN as a user dummy element and used to study failsafe characteristics of the C5A fuselage. The element is formulated from Reitsner's functional requiring that it satisfy compatability with the linear boundary displacement elements in NASTRAN. Its accuracy is demonstrated by analyzing for the stress intensity factors of two simple crack configurations for which there are classic solutions.
Mathematical modelling and linear stability analysis of laser fusion cutting
NASA Astrophysics Data System (ADS)
Hermanns, Torsten; Schulz, Wolfgang; Vossen, Georg; Thombansen, Ulrich
2016-06-01
A model for laser fusion cutting is presented and investigated by linear stability analysis in order to study the tendency for dynamic behavior and subsequent ripple formation. The result is a so called stability function that describes the correlation of the setting values of the process and the process' amount of dynamic behavior.
Linear regression analysis of survival data with missing censoring indicators.
Wang, Qihua; Dinse, Gregg E
2011-04-01
Linear regression analysis has been studied extensively in a random censorship setting, but typically all of the censoring indicators are assumed to be observed. In this paper, we develop synthetic data methods for estimating regression parameters in a linear model when some censoring indicators are missing. We define estimators based on regression calibration, imputation, and inverse probability weighting techniques, and we prove all three estimators are asymptotically normal. The finite-sample performance of each estimator is evaluated via simulation. We illustrate our methods by assessing the effects of sex and age on the time to non-ambulatory progression for patients in a brain cancer clinical trial. PMID:20559722
Real-time linear predictive analysis of speech using multimicroprocessors
Seethardman, S.; Radhakrishnan, T.; Suen, C.Y.
1982-01-01
Many applications of linear predictive coding (often known as LPC) of speech signals require a system capable of performing the complete LPC analysis in real time. This paper describes a pipeline network consisting of several general purpose microprocessors, primarily suitable for complete LPC analysis of a 10-pole model with a sampling frequency of 10 khz and a frame rate of 100 hz in real time. The proposed system is different from the previous systems, which either employed special purpose hardware or produced an analysis at a lower frame rate. 27 references.
Functional linear models for association analysis of quantitative traits.
Fan, Ruzong; Wang, Yifan; Mills, James L; Wilson, Alexander F; Bailey-Wilson, Joan E; Xiong, Momiao
2013-11-01
Functional linear models are developed in this paper for testing associations between quantitative traits and genetic variants, which can be rare variants or common variants or the combination of the two. By treating multiple genetic variants of an individual in a human population as a realization of a stochastic process, the genome of an individual in a chromosome region is a continuum of sequence data rather than discrete observations. The genome of an individual is viewed as a stochastic function that contains both linkage and linkage disequilibrium (LD) information of the genetic markers. By using techniques of functional data analysis, both fixed and mixed effect functional linear models are built to test the association between quantitative traits and genetic variants adjusting for covariates. After extensive simulation analysis, it is shown that the F-distributed tests of the proposed fixed effect functional linear models have higher power than that of sequence kernel association test (SKAT) and its optimal unified test (SKAT-O) for three scenarios in most cases: (1) the causal variants are all rare, (2) the causal variants are both rare and common, and (3) the causal variants are common. The superior performance of the fixed effect functional linear models is most likely due to its optimal utilization of both genetic linkage and LD information of multiple genetic variants in a genome and similarity among different individuals, while SKAT and SKAT-O only model the similarities and pairwise LD but do not model linkage and higher order LD information sufficiently. In addition, the proposed fixed effect models generate accurate type I error rates in simulation studies. We also show that the functional kernel score tests of the proposed mixed effect functional linear models are preferable in candidate gene analysis and small sample problems. The methods are applied to analyze three biochemical traits in data from the Trinity Students Study. PMID:24130119
Linear Analysis of Drift Ballooning Modes in Tokamak Edge Plasmas
NASA Astrophysics Data System (ADS)
Tangri, Varun; Kritz, Arnold; Rafiq, Tariq; Pankin, Alexei
2012-10-01
The H-mode pedestal structure depends on the linear stability of drift ballooning modes (DBMs) in many H-mode pedestal models. Integrated modeling that uses these pedestal models requires fast evaluation of linear stability of DBMs. Linear analysis of DBMs is also needed in the computations of effective diffusivities associated with anomalous transport that is driven by the DBMs in tokamak edge plasmas. In this study several numerical techniques of linear analysis of the DBMs are investigated. Differentiation matrix based spectral methods are used to compute the physical eigenvalues of the DBMs. The model for DBMs used here consists of six differential equations [T. Rafiq et al. Phys. Plasmas, 17, 082511, (2010)]. It is important to differentiate among non-physical (numerical) modes and physical modes. The determination of the number of eigenvalues is solved by a computation of the `nearest' and `ordinal' distances. The Finite Difference, Hermite and Sinc based differentiation matrices are used. It is shown that spectral collocation methods are more accurate than finite difference methods. The technique that has been developed for calculating eigenvalues is quite general and is relevant in the computation of other modes that utilize the ballooning mode formalism.
Linear stability analysis of swirling turbulent flows with turbulence models
NASA Astrophysics Data System (ADS)
Gupta, Vikrant; Juniper, Matthew
2013-11-01
In this paper, we consider the growth of large scale coherent structures in turbulent flows by performing linear stability analysis around a mean flow. Turbulent flows are characterized by fine-scale stochastic perturbations. The momentum transfer caused by these perturbations affects the development of larger structures. Therefore, in a linear stability analysis, it is important to include the perturbations' influence. One way to do this is to include a turbulence model in the stability analysis. This is done in the literature by using eddy viscosity models (EVMs), which are first order turbulence models. We extend this approach by using second order turbulence models, in this case explicit algebraic Reynolds stress models (EARSMs). EARSMs are more versatile than EVMs, in that they can be applied to a wider range of flows, and could also be more accurate. We verify our EARSM-based analysis by applying it to a channel flow and then comparing the results with those from an EVM-based analysis. We then apply the EARSM-based stability analysis to swirling pipe flows and Taylor-Couette flows, which demonstrates the main benefit of EARSM-based analysis. This project is supported by EPSRC and Rolls-Royce through a Dorothy Hodgkin Research Fellowship.
Linear and non-linear deflection analysis of thick rectangular plates. 2: Numerical applications
NASA Astrophysics Data System (ADS)
Bencharif, N.; Ng, S. F.
1994-03-01
Variational methods are widely used for the solution of complex differential equations in mechanics for which exact solutions are not possible. The finite difference method, although well known as an efficient numerical method, was applied in the past only for the analysis of linear and non-linear thin plates. In this paper the suitability of the method for the analysis of non-linear deflection of thick plates is studied for the first time. While there are major differences between small deflection and large deflection plate theories, the former can be treated as a particular case of the latter, when the centre deflection of the plate is less than or equal to 0.2-0.25 of the thickness of the plate. The finite difference method as applied here is a modified finite difference approach to the ordinary finite difference method generally used for the solution of thin plate problems. In this analysis thin plates are treated as a particular case of the corresponding thick plate when the boundary conditions of the plates are taken into account. The method is first applied to investigate the deflection behaviour of clamped and simply supported square isotropic thick plates. After the validity of the method is established, it is then extended to the solution of rectangular thick plates of various aspect ratios and thicknesses. Generally, beginning with the use of a limited number of mesh sizes for a given plate aspect ratio and boundary conditions, a general solution of the problem including the investigation of accuracy and convergence was extended to rectangular thick plates by providing more detailed functions satisfying the rectangular mesh sizes generated automatically by the program. Whenever possible results obtained by the present method are compared with existing solutions in the technical literature obtained by much more laborious methods and close agreements are found. The significant number of results presented here are not currently available in the technical
[Relations between biomedical variables: mathematical analysis or linear algebra?].
Hucher, M; Berlie, J; Brunet, M
1977-01-01
The authors, after a short reminder of one pattern's structure, stress on the possible double approach of relations uniting the variables of this pattern: use of fonctions, what is within the mathematical analysis sphere, use of linear algebra profiting by matricial calculation's development and automatiosation. They precise the respective interests on these methods, their bounds and the imperatives for utilization, according to the kind of variables, of data, and the objective for work, understanding phenomenons or helping towards decision. PMID:922174
Credibility analysis of risk classes by generalized linear model
NASA Astrophysics Data System (ADS)
Erdemir, Ovgucan Karadag; Sucu, Meral
2016-06-01
In this paper generalized linear model (GLM) and credibility theory which are frequently used in nonlife insurance pricing are combined for reliability analysis. Using full credibility standard, GLM is associated with limited fluctuation credibility approach. Comparison criteria such as asymptotic variance and credibility probability are used to analyze the credibility of risk classes. An application is performed by using one-year claim frequency data of a Turkish insurance company and results of credible risk classes are interpreted.
NASA Astrophysics Data System (ADS)
Kozubek, M.; Rozanov, E.; Krizan, P.
2014-09-01
The stratosphere is influenced by many external forcings (natural or anthropogenic). There are many studies which are focused on this problem and that is why we can compare our results with them. This study is focused on the variability and trends of temperature and circulation characteristics (zonal and meridional wind component) in connection with different phenomena variation in the stratosphere and lower mesosphere. We consider the interactions between the troposphere-stratosphere-lower mesosphere system and external and internal phenomena, e.g. solar cycle, QBO, NAO or ENSO using multiple linear techniques. The analysis was applied to the period 1979-2012 based on the current reanalysis data, mainly the MERRA reanalysis dataset (Modern Era Retrospective-analysis for Research and Applications) for pressure levels: 1000-0.1 hPa. We do not find a strong temperature signal for solar flux over the tropics about 30 hPa (ERA-40 results) but the strong positive signal has been observed near stratopause almost in the whole analyzed area. This could indicate that solar forcing is not represented well in the higher pressure levels in MERRA. The analysis of ENSO and ENSO Modoki shows that we should take into account more than one ENSO index for similar analysis. Previous studies show that the volcanic activity is important parameter. The signal of volcanic activity in MERRA is very weak and insignificant.
Testing the assumptions of linear prediction analysis in normal vowels
NASA Astrophysics Data System (ADS)
Little, M. A.
This paper develops an improved surrogate data test to show experimental evidence, for all the simple vowels of US English, for both male and female speakers, that Gaussian linear prediction analysis, a ubiquitous technique in current speech technologies, cannot be used to extract all the dynamical structure of real speech time series. The test provides robust evidence undermining the validity of these linear techniques, supporting the assumptions of either dynamical nonlinearity and/or non-Gaussianity common to more recent, complex, efforts at dynamical modelling speech time series. However, an additional finding is that the classical assumptions cannot be ruled out entirely, and plausible evidence is given to explain the success of the linear Gaussian theory as a weak approximation to the true, nonlinear/non-Gaussian dynamics. This supports the use of appropriate hybrid linear/nonlinear/non-Gaussian modelling. With a calibrated calculation of statistic and particular choice of experimental protocol, some of the known systematic problems of the method of surrogate data testing are circumvented to obtain results to support the conclusions to a high level of significance.
A linearized Euler analysis of unsteady flows in turbomachinery
NASA Technical Reports Server (NTRS)
Hall, Kenneth C.; Crawley, Edward F.
1987-01-01
A method for calculating unsteady flows in cascades is presented. The model, which is based on the linearized unsteady Euler equations, accounts for blade loading shock motion, wake motion, and blade geometry. The mean flow through the cascade is determined by solving the full nonlinear Euler equations. Assuming the unsteadiness in the flow is small, then the Euler equations are linearized about the mean flow to obtain a set of linear variable coefficient equations which describe the small amplitude, harmonic motion of the flow. These equations are discretized on a computational grid via a finite volume operator and solved directly subject to an appropriate set of linearized boundary conditions. The steady flow, which is calculated prior to the unsteady flow, is found via a Newton iteration procedure. An important feature of the analysis is the use of shock fitting to model steady and unsteady shocks. Use of the Euler equations with the unsteady Rankine-Hugoniot shock jump conditions correctly models the generation of steady and unsteady entropy and vorticity at shocks. In particular, the low frequency shock displacement is correctly predicted. Results of this method are presented for a variety of test cases. Predicted unsteady transonic flows in channels are compared to full nonlinear Euler solutions obtained using time-accurate, time-marching methods. The agreement between the two methods is excellent for small to moderate levels of flow unsteadiness. The method is also used to predict unsteady flows in cascades due to blade motion (flutter problem) and incoming disturbances (gust response problem).
Density functional steric analysis of linear and branched alkanes.
Ess, Daniel H; Liu, Shubin; De Proft, Frank
2010-12-16
Branched alkane hydrocarbons are thermodynamically more stable than straight-chain linear alkanes. This thermodynamic stability is also manifest in alkane bond separation energies. To understand the physical differences between branched and linear alkanes, we have utilized a novel density functional theory (DFT) definition of steric energy based on the Weizäcker kinetic energy. Using the M06-2X functional, the total DFT energy was partitioned into a steric energy term (E(s)[ρ]), an electrostatic energy term (E(e)[ρ]), and a fermionic quantum energy term (E(q)[ρ]). This analysis revealed that branched alkanes have less (destabilizing) DFT steric energy than linear alkanes. The lower steric energy of branched alkanes is mitigated by an equal and opposite quantum energy term that contains the Pauli component of the kinetic energy and exchange-correlation energy. Because the steric and quantum energy terms cancel, this leaves the electrostatic energy term that favors alkane branching. Electrostatic effects, combined with correlation energy, explains why branched alkanes are more stable than linear alkanes. PMID:21086970
A linearized Euler analysis of unsteady flows in turbomachinery
NASA Astrophysics Data System (ADS)
Hall, Kenneth C.; Crawley, Edward F.
1987-06-01
A method for calculating unsteady flows in cascades is presented. The model, which is based on the linearized unsteady Euler equations, accounts for blade loading shock motion, wake motion, and blade geometry. The mean flow through the cascade is determined by solving the full nonlinear Euler equations. Assuming the unsteadiness in the flow is small, then the Euler equations are linearized about the mean flow to obtain a set of linear variable coefficient equations which describe the small amplitude, harmonic motion of the flow. These equations are discretized on a computational grid via a finite volume operator and solved directly subject to an appropriate set of linearized boundary conditions. The steady flow, which is calculated prior to the unsteady flow, is found via a Newton iteration procedure. An important feature of the analysis is the use of shock fitting to model steady and unsteady shocks. Use of the Euler equations with the unsteady Rankine-Hugoniot shock jump conditions correctly models the generation of steady and unsteady entropy and vorticity at shocks. In particular, the low frequency shock displacement is correctly predicted. Results of this method are presented for a variety of test cases. Predicted unsteady transonic flows in channels are compared to full nonlinear Euler solutions obtained using time-accurate, time-marching methods. The agreement between the two methods is excellent for small to moderate levels of flow unsteadiness. The method is also used to predict unsteady flows in cascades due to blade motion (flutter problem) and incoming disturbances (gust response problem).
Density Functional Steric Analysis of Linear and Branched Alkanes
Ess, Daniel H.; Liu, Shubin; De Proft, Frank
2010-11-18
Branched alkane hydrocarbons are thermodynamically more stable than straight-chain linear alkanes. This thermodynamic stability is also manifest in alkane bond separation energies. To understand the physical differences between branched and linear alkanes, we have utilized a novel density functional theory (DFT) definition of steric energy based on the Weizäcker kinetic energy. Using the M06-2X functional, the total DFT energy was partitioned into a steric energy term (E_{e}[[ρ]), an electrostatic energy term (E_{e}[ρ]), and a fermionic quantum energy term (E_{q}[[ρ]). This analysis revealed that branched alkanes have less (destabilizing) DFT steric energy than linear alkanes. The lower steric energy of branched alkanes is mitigated by an equal and opposite quantum energy term that contains the Pauli component of the kinetic energy and exchange-correlation energy. Because the steric and quantum energy terms cancel, this leaves the electrostatic energy term that favors alkane branching. Electrostatic effects, combined with correlation energy, explains why branched alkanes are more stable than linear alkanes.
Feasible logic Bell-state analysis with linear optics.
Zhou, Lan; Sheng, Yu-Bo
2016-01-01
We describe a feasible logic Bell-state analysis protocol by employing the logic entanglement to be the robust concatenated Greenberger-Horne-Zeilinger (C-GHZ) state. This protocol only uses polarization beam splitters and half-wave plates, which are available in current experimental technology. We can conveniently identify two of the logic Bell states. This protocol can be easily generalized to the arbitrary C-GHZ state analysis. We can also distinguish two N-logic-qubit C-GHZ states. As the previous theory and experiment both showed that the C-GHZ state has the robustness feature, this logic Bell-state analysis and C-GHZ state analysis may be essential for linear-optical quantum computation protocols whose building blocks are logic-qubit entangled state. PMID:26877208
Feasible logic Bell-state analysis with linear optics
NASA Astrophysics Data System (ADS)
Zhou, Lan; Sheng, Yu-Bo
2016-02-01
We describe a feasible logic Bell-state analysis protocol by employing the logic entanglement to be the robust concatenated Greenberger-Horne-Zeilinger (C-GHZ) state. This protocol only uses polarization beam splitters and half-wave plates, which are available in current experimental technology. We can conveniently identify two of the logic Bell states. This protocol can be easily generalized to the arbitrary C-GHZ state analysis. We can also distinguish two N-logic-qubit C-GHZ states. As the previous theory and experiment both showed that the C-GHZ state has the robustness feature, this logic Bell-state analysis and C-GHZ state analysis may be essential for linear-optical quantum computation protocols whose building blocks are logic-qubit entangled state.
Aerodynamic preliminary analysis system. Part 1: Theory. [linearized potential theory
NASA Technical Reports Server (NTRS)
Bonner, E.; Clever, W.; Dunn, K.
1978-01-01
A comprehensive aerodynamic analysis program based on linearized potential theory is described. The solution treats thickness and attitude problems at subsonic and supersonic speeds. Three dimensional configurations with or without jet flaps having multiple non-planar surfaces of arbitrary planform and open or closed slender bodies of non-circular contour may be analyzed. Longitudinal and lateral-directional static and rotary derivative solutions may be generated. The analysis was implemented on a time sharing system in conjunction with an input tablet digitizer and an interactive graphics input/output display and editing terminal to maximize its responsiveness to the preliminary analysis problem. Nominal case computation time of 45 CPU seconds on the CDC 175 for a 200 panel simulation indicates the program provides an efficient analysis for systematically performing various aerodynamic configuration tradeoff and evaluation studies.
Feasible logic Bell-state analysis with linear optics
Zhou, Lan; Sheng, Yu-Bo
2016-01-01
We describe a feasible logic Bell-state analysis protocol by employing the logic entanglement to be the robust concatenated Greenberger-Horne-Zeilinger (C-GHZ) state. This protocol only uses polarization beam splitters and half-wave plates, which are available in current experimental technology. We can conveniently identify two of the logic Bell states. This protocol can be easily generalized to the arbitrary C-GHZ state analysis. We can also distinguish two N-logic-qubit C-GHZ states. As the previous theory and experiment both showed that the C-GHZ state has the robustness feature, this logic Bell-state analysis and C-GHZ state analysis may be essential for linear-optical quantum computation protocols whose building blocks are logic-qubit entangled state. PMID:26877208
Linear discriminant analysis with misallocation in training samples
NASA Technical Reports Server (NTRS)
Chhikara, R. (Principal Investigator); Mckeon, J.
1982-01-01
Linear discriminant analysis for a two-class case is studied in the presence of misallocation in training samples. A general appraoch to modeling of mislocation is formulated, and the mean vectors and covariance matrices of the mixture distributions are derived. The asymptotic distribution of the discriminant boundary is obtained and the asymptotic first two moments of the two types of error rate given. Certain numerical results for the error rates are presented by considering the random and two non-random misallocation models. It is shown that when the allocation procedure for training samples is objectively formulated, the effect of misallocation on the error rates of the Bayes linear discriminant rule can almost be eliminated. If, however, this is not possible, the use of Fisher rule may be preferred over the Bayes rule.
Meshless analysis of shear deformable shells: the linear model
NASA Astrophysics Data System (ADS)
Costa, Jorge C.; Tiago, Carlos M.; Pimenta, Paulo M.
2013-10-01
This work develops a kinematically linear shell model departing from a consistent nonlinear theory. The initial geometry is mapped from a flat reference configuration by a stress-free finite deformation, after which, the actual shell motion takes place. The model maintains the features of a complete stress-resultant theory with Reissner-Mindlin kinematics based on an inextensible director. A hybrid displacement variational formulation is presented, where the domain displacements and kinematic boundary reactions are independently approximated. The resort to a flat reference configuration allows the discretization using 2-D Multiple Fixed Least-Squares (MFLS) on the domain. The consistent definition of stress resultants and consequent plane stress assumption led to a neat formulation for the analysis of shells. The consistent linear approximation, combined with MFLS, made possible efficient computations with a desired continuity degree, leading to smooth results for the displacement, strain and stress fields, as shown by several numerical examples.
Near infrared reflectance analysis by Gauss-Jordan linear algebra
NASA Astrophysics Data System (ADS)
Honigs, D. E.; Freelin, J. M.; Hieftje, G. M.
1983-02-01
Near-infrared reflectance analysis (NIRA) is an analytical technique that uses the near-infrared diffuse reflectance of a sample at several discrete wavelengths to predict the concentration of one or more of the chemical species in that sample. However, because near-infrared bands from solid samples are both abundant and broad, the reflectance at a given wavelength usually contains contributions from several sample components, requiring extensive calculations on overlapped bands. In the present study, these calculations have been performed using an approach similar to that employed in multi-component spectrophotometry, but with Gauss-Jordan linear algebra serving as the computational vehicle. Using this approach, correlations for percent protein in wheat flour and percent benzene in hydrocarbons have been obtained and are evaluated. The advantages of a linear-algebra approach over the common one employing stepwise regression are explored.
ERIC Educational Resources Information Center
Gonzalez-Vega, Laureano
1999-01-01
Using a Computer Algebra System (CAS) to help with the teaching of an elementary course in linear algebra can be one way to introduce computer algebra, numerical analysis, data structures, and algorithms. Highlights the advantages and disadvantages of this approach to the teaching of linear algebra. (Author/MM)
Analysis of the transverse electromagnetic mode linearly tapered slot antenna
NASA Technical Reports Server (NTRS)
Janaswamy, R.; Schaubert, D. H.; Pozar, D. M.
1986-01-01
A theoretical analysis of the radiation characteristics of the transverse electromagnetic mode linearly tapered slot antenna is presented. The theory presented is valid for antennas with air dielectric and forms the basis for analysis of the more popular dielectric-supported antennas. The method of analysis involves two steps. In the first step, the aperture distribution in the flared slot is determined. In the second step, the equivalent magnetic current in the slot is treated as radiating in the presence of a conducting half-plane, and the far-field components are obtained. Design curves for the variation of the 3-dB and 10-dB beamwidths as a function of the antenna length, with the flare angle as a parameter, are presented.
Analysis of the Tem Mode Linearly Tapered Slot Antenna
NASA Technical Reports Server (NTRS)
Janaswamy, R.; Schaubert, D. H.; Pozar, D. M.
1985-01-01
The theoretical analysis of the radiation characteristics of the TEM mode Linearly Tapered Slot Antenna (LTSA) is presented. The theory presented is valid for antennas with air dielectric and forms the basis for analysis of the more popular dielectric-supported antennas. The method of analysis involves two steps. In the first step, the aperture distribution in the flared slot is determined. In the second step, the equivalent magnetic current in the slot is treated as radiating in the presence of a conducting half-plane and the far-field components are obtained. Detailed comparison with experiment is made and excellent agreement is obtained. Design curves for the variation of the 3 dB and 10 dB beamwidths as a function of the antenna length, with the flare angle as a parameter, are presented.
NASA Astrophysics Data System (ADS)
Kbaier Ben Ismail, Dhouha; Lazure, Pascal; Puillat, Ingrid
2016-10-01
In marine sciences, many fields display high variability over a large range of spatial and temporal scales, from seconds to thousands of years. The longer recorded time series, with an increasing sampling frequency, in this field are often nonlinear, nonstationary, multiscale and noisy. Their analysis faces new challenges and thus requires the implementation of adequate and specific methods. The objective of this paper is to highlight time series analysis methods already applied in econometrics, signal processing, health, etc. to the environmental marine domain, assess advantages and inconvenients and compare classical techniques with more recent ones. Temperature, turbidity and salinity are important quantities for ecosystem studies. The authors here consider the fluctuations of sea level, salinity, turbidity and temperature recorded from the MAREL Carnot system of Boulogne-sur-Mer (France), which is a moored buoy equipped with physico-chemical measuring devices, working in continuous and autonomous conditions. In order to perform adequate statistical and spectral analyses, it is necessary to know the nature of the considered time series. For this purpose, the stationarity of the series and the occurrence of unit-root are addressed with the Augmented-Dickey Fuller tests. As an example, the harmonic analysis is not relevant for temperature, turbidity and salinity due to the nonstationary condition, except for the nearly stationary sea level datasets. In order to consider the dominant frequencies associated to the dynamics, the large number of data provided by the sensors should enable the estimation of Fourier spectral analysis. Different power spectra show a complex variability and reveal an influence of environmental factors such as tides. However, the previous classical spectral analysis, namely the Blackman-Tukey method, requires not only linear and stationary data but also evenly-spaced data. Interpolating the time series introduces numerous artifacts to the
Linear stability analysis for hydrothermal alteration of kimberlitic rocks
NASA Astrophysics Data System (ADS)
Afanasyev, Andrey; Belyaeva, Ekaterina
2016-04-01
The influx of groundwater into hot kimberlite deposits results in the reaction of water with olivine-rich rocks. The products of the reaction are serpentine and release of latent heat. The rise of temperature due to the heat release increases the rate of the reaction. Under certain conditions, this self-speeding up of the reaction can result in instabilities associated with a significantly higher final serpentinisation in slightly warmer regions of the kimberlite deposit. We conduct linear stability analysis of serpentinisation in an isolated volume of porous kimberlitic rocks saturated with water and an inert gas. There is a counteracting interplay between the heat release tending to destabilise the uniform distribution of parameters and the heat conduction tending to stabilise it by smoothing out temperature perturbations. We determine the critical spatial scale separating the parameters where one phenomenon dominates over another. The perturbations of longer-than-critical length grow, whereas the perturbations of shorter-than-critical length fade. The analytical results of the linear stability analysis are supported by direct numerical simulations using a full nonlinear model. Keywords: Hydrothermal systems, volcaniclastic deposits, phase transitions, instability analysis, numerical solutions
Modal analysis to accommodate slap in linear structures.
Segalman, Daniel Joseph; Roy, Anthony M.; Starr, Michael James
2004-10-01
The generalized momentum balance (GMB) methods, explored chiefly by Shabana and his co-workers, treat slap or collision in linear structures as sequences of impulses, thereby maintaining the linearity of the structures throughout. Further, such linear analysis is facilitated by modal representation of the structures. These methods are discussed here and extended. Simulations on a simple two-rod problem demonstrate how this modal impulse approximation affects the system both directly after each impulse as well as over the entire collision. Furthermore, these simulations illustrate how the GMB results differ from the exact solution and how mitigation of these artifacts is achieved. Another modal method discussed in this paper is the idea of imposing piecewise constant forces over short, yet finite, time intervals during contact. The derivation of this method is substantially different than that of the GMB method, yet the numerical results show similar behavior, adding credence to both models. Finally, a novel method combining these two approaches is introduced. The new method produces physically reasonable results that are numerically very close to the exact solution of the collision of two rods. This approach avoids most of the non physical, numerical artifacts of interpenetration or chatter present in the first two methods.
Coupling analysis of linear vibration energy harvesting systems
NASA Astrophysics Data System (ADS)
Wang, Xu; Liang, Xingyu; Shu, Gequn; Watkins, Simon
2016-03-01
This paper has disclosed the relationship of vibration energy harvester performance with dimensionless force factor. Numerical ranges of the dimensionless force factor have been defined for cases of weak, moderate and strong coupling. The relationships of coupling loss factor, dimensionless force factor, critical coupling strength, coupling quotient, electro-mechanical coupling factor, damping loss factor and modal densities have been established in linear vibration energy harvester systems. The new contribution of this paper is to determine a frequency range where the vibration energy harvesting systems are in a weak coupling and the statistical energy analysis is applicable.
Analysis of integral controls in linear quadratic regulator design
NASA Technical Reports Server (NTRS)
Slater, G. L.
1979-01-01
The application of linear optimal control to the design of systems with integral control action on specified outputs is considered. Using integral terms in a quadratic performance index, an asymptotic analysis is used to determine the effect of variable quadratic weights on the eigenvalues and eigenvectors of the closed loop system. It is shown that for small integral terms the placement of integrator poles and gain calculation can be effectively decoupled from placement of the primary system eigenvalues. This technique is applied to the design of integral controls for a STOL aircraft outer loop guidance system.
Dynamic analysis of space-related linear and non-linear structures
Bosela, P.A.; Shaker, F.J.; Fertis, D.G.
1990-01-01
In order to be cost effective, space structures must be extremely light weight, and subsequently, very flexible structures. The power system for Space Station Freedom is such a structure. Each array consists of a deployable truss mast and a split blanket of photovoltaic solar collectors. The solar arrays are deployed in orbit, and the blanket is stretched into position as the mast is extended. Geometric stiffness due to the preload make this an interesting non-linear problem. The space station will be subjected to various dynamic loads, during shuttle docking, solar tracking, attitude adjustment, etc. Accurate prediction of the natural frequencies and mode shapes of the space station components, including the solar arrays, is critical for determining the structural adequacy of the components, and for designing a dynamic controls system. The process used in developing and verifying the finite element dynamic model of the photo-voltaic arrays is documented. Various problems were identified, such as grounding effects due to geometric stiffness, large displacement effects, and pseudo-stiffness (grounding) due to lack of required rigid body modes. Analysis techniques, such as development of rigorous solutions using continuum mechanics, finite element solution sequence altering, equivalent systems using a curvature basis, Craig-Bampton superelement approach, and modal ordering schemes were utilized. The grounding problems associated with the geometric stiffness are emphasized.
Dynamic analysis of space-related linear and non-linear structures
NASA Technical Reports Server (NTRS)
Bosela, Paul A.; Shaker, Francis J.; Fertis, Demeter G.
1990-01-01
In order to be cost effective, space structures must be extremely light weight, and subsequently, very flexible structures. The power system for Space Station Freedom is such a structure. Each array consists of a deployable truss mast and a split blanket of photo-voltaic solar collectors. The solar arrays are deployed in orbit, and the blanket is stretched into position as the mast is extended. Geometric stiffness due to the preload make this an interesting non-linear problem. The space station will be subjected to various dynamic loads, during shuttle docking, solar tracking, attitude adjustment, etc. Accurate prediction of the natural frequencies and mode shapes of the space station components, including the solar arrays, is critical for determining the structural adequacy of the components, and for designing a dynamic control system. The process used in developing and verifying the finite element dynamic model of the photo-voltaic arrays is documented. Various problems were identified, such as grounding effects due to geometric stiffness, large displacement effects, and pseudo-stiffness (grounding) due to lack of required rigid body modes. Analysis techniques, such as development of rigorous solutions using continuum mechanics, finite element solution sequence altering, equivalent systems using a curvature basis, Craig-Bampton superelement approach, and modal ordering schemes were utilized. The grounding problems associated with the geometric stiffness are emphasized.
Linear stability analysis for hydrothermal alteration of kimberlitic rocks
NASA Astrophysics Data System (ADS)
Afanasyev, Andrey; Belyaeva, Ekaterina
2016-06-01
The influx of groundwater into hot kimberlite deposits results in the reaction of water with olivine-rich rocks. The products of the reaction are serpentine and release of latent heat. The rise of temperature due to the heat release increases the rate of the reaction. Under certain conditions, this self-speeding up of the reaction can result in instabilities associated with a significantly higher final serpentinization in slightly warmer regions of the kimberlite deposit. We conduct linear stability analysis of serpentinization in an isolated volume of porous kimberlitic rocks saturated with water and an inert gas. There is a counteracting interplay between the heat release tending to destabilize the uniform distribution of parameters and the heat conduction tending to stabilize it by smoothing out temperature perturbations. We determine the critical spatial scale separating the parameters where one phenomenon dominates over another. The perturbations of longer-than-critical length grow, whereas the perturbations of shorter-than-critical length fade. The analytical results of the linear stability analysis are supported by direct numerical simulations using a full nonlinear model.
Linear Stability Analysis of an Acoustically Vaporized Droplet
NASA Astrophysics Data System (ADS)
Siddiqui, Junaid; Qamar, Adnan; Samtaney, Ravi
2015-11-01
Acoustic droplet vaporization (ADV) is a phase transition phenomena of a superheat liquid (Dodecafluoropentane, C5F12) droplet to a gaseous bubble, instigated by a high-intensity acoustic pulse. This approach was first studied in imaging applications, and applicable in several therapeutic areas such as gas embolotherapy, thrombus dissolution, and drug delivery. High-speed imaging and theoretical modeling of ADV has elucidated several physical aspects, ranging from bubble nucleation to its subsequent growth. Surface instabilities are known to exist and considered responsible for evolving bubble shapes (non-spherical growth, bubble splitting and bubble droplet encapsulation). We present a linear stability analysis of the dynamically evolving interfaces of an acoustically vaporized micro-droplet (liquid A) in an infinite pool of a second liquid (liquid B). We propose a thermal ADV model for the base state. The linear analysis utilizes spherical harmonics (Ynm, of degree m and order n) and under various physical assumptions results in a time-dependent ODE of the perturbed interface amplitudes (one at the vapor/liquid A interface and the other at the liquid A/liquid B interface). The perturbation amplitudes are found to grow exponentially and do not depend on m. Supported by KAUST Baseline Research Funds.
DYNAMIC NON LINEAR IMPACT ANALYSIS OF FUEL CASK CONTAINMENT VESSELS
Leduc, D
2008-06-10
Large fuel casks present challenges when evaluating their performance in the accident sequence specified in 10CFR 71. Testing is often limited because of cost, difficulty in preparing test units and the limited availability of facilities which can carry out such tests. In the past, many casks were evaluated without testing using simplified analytical methods. This paper details the use of dynamic non-linear analysis of large fuel casks using advanced computational techniques. Results from the dynamic analysis of two casks, the T-3 Spent Fuel Cask and the Hanford Un-irradiated Fuel Package are examined in detail. These analyses are used to fully evaluate containment vessel stresses and strains resulting from complex loads experienced by cask components during impacts. Importantly, these advanced analytical analyses are capable of examining stresses in key regions of the cask including the cask closure. This paper compares these advanced analytical results with the results of simplified cask analyses like those detailed in NUREG 3966.
NASA Technical Reports Server (NTRS)
Watson, Brian C.; Kamat, Manohar P.
1992-01-01
Concurrent computing environments provide the means to achieve very high performance for finite element analysis of systems, provided the algorithms take advantage of multiple processors. The authors have examined several algorithms for both linear and nonlinear finite element analysis. The performance of these algorithms on an Alliant FX/80 parallel supercomputer has been studied. For single load case linear analysis, the optimal solution algorithm is strongly problem dependent. For multiple load cases or nonlinear analysis through a modified Newton-Raphson method, decomposition algorithms are shown to have a decided advantage over element-by-element preconditioned conjugate gradient algorithms.
Interactive Analysis of Hyperspectral Data under Linearity Constraints
NASA Astrophysics Data System (ADS)
Schmidt, A.; Treguier, E.; Schmidt, F.; Moussaoui, S.; Pelloquin, C.
2010-12-01
Large data sets delivered by imaging spectrometers are interesting in many ways in the Planetary Sciences. Due to the size of the data and lack of ground truth, which often prohibit conventional exploratory data analysis methods, interactive but unsupervised analysis methods could be a way of discovering relevant information about the sources that make up the data. In this work, we investigate some of the opportunities and limitations of such analyses based on non-negative matrix approximation in planetary settings. Since typically there often is no ground truth to compare to, the degrees of freedom inherent in the aforementioned approximation techniques often has to be constrained by users to discover physically valid sources and patterns. One way of going about this is to present users with different valid solutions have them choose the one or ones that fit their knowledge of the environment best. Recent developments have made it possible to exploit linear mixing constraints and present results to users in real or near-real time; thus, the approach has become practicable. The general setting of the problem is as follows: By considering P pixels of an hyperspectral image acquired at L frequency bands, the observed spectra are gathered in a PxL data matrix X. Each row of this matrix contains a measured spectrum at a pixel with spatial index p=1..P. According to the linear mixing model, the p-th spectrum, 1<=p<=P, can be expressed as a linear combination of r, 1<=r<=R, pure spectra of the surface components. Thus, X=AxS+E, E being an error matrix, should be minimised, where X, A, and S have only non-negative entries. The rows of matrix S now contain the pure surface spectra of the R components, and each entry of A corresponds to the abundance of the r-th component in pixel with spatial index p. For a qualitative and quantitative description of the observed scene composition, the estimation problem consists of finding matrices S and A which allow to explain the data
Detection of linear impermeable barriers by transient pressure analysis
Martinez-Romero, N.; Cinco-Ley, H.
1983-05-01
This work presents a new simple method to detect a linear impermeable barrier by analysis of transient pressure data. This technique is based on the desuperposition method (negative superposition) discussed by some authors and considers the calculation the pressure change caused by the presence of the barrier. The pressure change is analyzed to estimate the distance between the well and the barrier by using the type curve matching technique. Type curves are provided for both drawdown and build tests. The advantage of the technique presented in this work is that pressure data can be analyzed even if the second semilog straight line, whose slope is twice the slope of the first semilog straight line, is not reached. Examples of application are discussed for illustration.
Analysis of Power Model for Linear Plasma Device
NASA Astrophysics Data System (ADS)
Zhang, Weiwei; Deng, Baiquan; Zuo, Haoyi; Zheng, Xianjun; Cao, Xiaogang; Xue, Xiaoyan; Ou, Wei; Cao, Zhi; Gou, Fujun
2016-08-01
A single cathode linear plasma device has been designed and constructed to investigate the interactions between plasma and materials at the Sichuan University. In order to further investigate the Ohmic power of the device, the output heat load on the specimen and electric potential difference (between cathode and anode) have been tested under different discharge currents. This special power distribution in the radial direction of the plasma discharge channel has also been discussed and described by some improved integral equations in this paper; it can be further simplified as P ∝ α‑2 in one-parameter. Besides, we have measured the power loss of the channel under different discharge currents by the calorimetric method, calculated the effective power of the device and evaluated the performances of the plasma device through the power efficiency analysis. supported by International Thermonuclear Experimental Reactor (ITER) Program (No. 2013GB114003) and National Natural Science Foundation of China (Nos. 11275135 and 11475122)
Linear analysis of incompressible Rayleigh-Taylor instability in solids.
Piriz, A R; Cela, J J López; Tahir, N A
2009-10-01
The study of the linear stage of the incompressible Rayleigh-Taylor instability in elastic-plastic solids is performed by considering thick plates under a constant acceleration that is also uniform except for a small sinusoidal ripple in the horizontal plane. The analysis is carried out by using an analytical model based on the Newton second law and it is complemented with extensive two-dimensional numerical simulations. The conditions for marginal stability that determine the instability threshold are derived. Besides, the boundary for the transition from the elastic to the plastic regime is obtained and it is demonstrated that such a transition is not a sufficient condition for instability. The model yields complete analytical solutions for the perturbation amplitude evolution and reveals the main physical process that governs the instability. The theory is in general agreement with the numerical simulations and provides useful quantitative results. Implications for high-energy-density-physics experiments are also discussed. PMID:19905434
Adjusting for matching and covariates in linear discriminant analysis
Asafu-Adjei, Josephine K.; Sampson, Allan R.; Sweet, Robert A.; Lewis, David A.
2013-01-01
In studies that compare several diagnostic or treatment groups, subjects may not only be measured on a certain set of feature variables, but also be matched on a number of demographic characteristics and measured on additional covariates. Linear discriminant analysis (LDA) is sometimes used to identify which feature variables best discriminate among groups, while accounting for the dependencies among the feature variables. We present a new approach to LDA for multivariate normal data that accounts for the subject matching used in a particular study design, as well as covariates not used in the matching. Applications are given for post-mortem tissue data with the aim of comparing neurobiological characteristics of subjects with schizophrenia with those of normal controls, and for a post-mortem tissue primate study comparing brain biomarker measurements across three treatment groups. We also investigate the performance of our approach using a simulation study. PMID:23640791
Spectrum Analysis of the Linearized Relativistic Landau Equation
NASA Astrophysics Data System (ADS)
Luo, Lan; Yu, Hongjun
2016-05-01
In this work we prove the complete spectrum structure of the linearized relativistic Landau equation in L^2 by using the semigroup theory and the linear operator perturbation theory. Our results include the physical interesting Coulombic interaction.
Linear harmonic analysis of free-piston Stirling engines
Chen, N.C.J.; Griffin, F.P.
1986-06-01
The equations that govern the behavior of free-piston Stirling engines are nonlinear differential equations. Traditional solution methods have been time-stepping integrations that can be plagued by numerical instabilities and can use large amounts of computer time. Closed-form analytical solutions are possible if the working gas behaves isothermally or if the nonlinear terms in the governing equations are replaced with accurate approximations. An almost closed-form solution method, called the linear harmonic analysis (LHA), has been developed for free-piston Stirling engine applications by representing all of the periodic variables with harmonic functions. The solution method accounts for the important thermodynamic losses that are coupled together in free-piston engines, yet it is efficient enough for optimization studies. The LHA method was compared with a standard numerical integration method to verify its mathematical accuracy. The LHA and numerical predictions for a sample free-piston Stirling engine configuration differed by <5% for all important parameters. Sensitivity studies using the LHA method have also shown that the thermodynamic loss assumptions used in an analysis can have a significant impact on the predicted dynamic behavior of a free-piston Stirling engine.
Moment method analysis of linearly tapered slot antennas
NASA Technical Reports Server (NTRS)
Koeksal, Adnan
1993-01-01
A method of moments (MOM) model for the analysis of the Linearly Tapered Slot Antenna (LTSA) is developed and implemented. The model employs an unequal size rectangular sectioning for conducting parts of the antenna. Piecewise sinusoidal basis functions are used for the expansion of conductor current. The effect of the dielectric is incorporated in the model by using equivalent volume polarization current density and solving the equivalent problem in free-space. The feed section of the antenna including the microstripline is handled rigorously in the MOM model by including slotline short circuit and microstripline currents among the unknowns. Comparison with measurements is made to demonstrate the validity of the model for both the air case and the dielectric case. Validity of the model is also verified by extending the model to handle the analysis of the skew-plate antenna and comparing the results to those of a skew-segmentation modeling results of the same structure and to available data in the literature. Variation of the radiation pattern for the air LTSA with length, height, and taper angle is investigated, and the results are tabulated. Numerical results for the effect of the dielectric thickness and permittivity are presented.
Stratospheric ozone time series analysis using dynamical linear models
NASA Astrophysics Data System (ADS)
Laine, Marko; Kyrölä, Erkki
2013-04-01
We describe a hierarchical statistical state space model for ozone profile time series. The time series are from satellite measurements by the SAGE II and GOMOS instruments spanning years 1984-2012. The original data sets are combined and gridded monthly using 10 degree latitude bands, and covering 20-60 km with 1 km vertical spacing. Model components include level, trend, seasonal effect with solar activity, and quasi biennial oscillations as proxy variables. A typical feature of an atmospheric time series is that they are not stationary but exhibit both slowly varying and abrupt changes in the distributional properties. These are caused by external forcing such as changes in the solar activity or volcanic eruptions. Further, the data sampling is often nonuniform, there are data gaps, and the uncertainty of the observations can vary. When observations are combined from various sources there will be instrument and retrieval method related biases. The differences in sampling lead also to uncertainties. Standard classical ARIMA type of statistical time series methods are mostly useless for atmospheric data. A more general approach makes use of dynamical linear models and Kalman filter type of sequential algorithms. These state space models assume a linear relationship between the unknown state of the system and the observations and for the process evolution of the hidden states. They are still flexible enough to model both smooth trends and sudden changes. The above mentioned methodological challenges are discussed, together with analysis of change points in trends related to recovery of stratospheric ozone. This work is part of the ESA SPIN and ozone CCI projects.
State-variable analysis of non-linear circuits with a desk computer
NASA Technical Reports Server (NTRS)
Cohen, E.
1981-01-01
State variable analysis was used to analyze the transient performance of non-linear circuits on a desk top computer. The non-linearities considered were not restricted to any circuit element. All that is required for analysis is the relationship defining each non-linearity be known in terms of points on a curve.
Prediction of sumoylation sites in proteins using linear discriminant analysis.
Xu, Yan; Ding, Ya-Xin; Deng, Nai-Yang; Liu, Li-Ming
2016-01-15
Sumoylation is a multifunctional post-translation modification (PTM) in proteins by the small ubiquitin-related modifiers (SUMOs), which have relations to ubiquitin in molecular structure. Sumoylation has been found to be involved in some cellular processes. It is very significant to identify the exact sumoylation sites in proteins for not only basic researches but also drug developments. Comparing with time exhausting experiment methods, it is highly desired to develop computational methods for prediction of sumoylation sites as a complement to experiment in the post-genomic age. In this work, three feature constructions (AAIndex, position-specific amino acid propensity and modification of composition of k-space amino acid pairs) and five different combinations of them were used to construct features. At last, 178 features were selected as the optimal features according to the Mathew's correlation coefficient values in 10-fold cross validation based on linear discriminant analysis. In 10-fold cross-validation on the benchmark dataset, the accuracy and Mathew's correlation coefficient were 86.92% and 0.6845. Comparing with those existing predictors, SUMO_LDA showed its better performance. PMID:26432000
Linear stability analysis of inclined two-layer stratified flows
NASA Astrophysics Data System (ADS)
Negretti, M. Eletta; Socolofsky, Scott A.; Jirka, Gerhard H.
2008-09-01
Two-layer stratified flows are commonly observed in geophysical and environmental contexts. At the interface between the two layers, both velocity shear and buoyancy interplay, resulting in various modes of instability. Results from a temporal linear stability analysis of a two-layer stratified exchange flow under the action of a mean advection are presented, investigating the effect of a mild bottom slope on the stability of the interface. The spatial acceleration is directly included in the governing stability equations. The results demonstrate that increasing the bottom slope has a similar effect on the stability of the flow as does increasing the ratio R of the thickness of the velocity mixing layer δν to that of the density layer δρ as it causes the flow to be more unstable to the Kelvin-Helmholtz instabilities. The transition from Kelvin-Helmholtz modes to stable flow occurs at lower Richardson numbers and wavenumbers compared to the horizontal two-layer flow. Kelvin-Helmholtz modes are decreasingly amplified for 1
An analysis of Seasat SAR for detecting geologic linears
NASA Technical Reports Server (NTRS)
Yan, S.-Y.; Philipson, W. R.; Teng, W. L.
1982-01-01
The value of Seasat synthetic aperture radar (SAR) imagery for detecting geologic linears was assessed in a study of an 89,000 sq km section of New York's Adirondack Mountains. A photographic print of optically processed, 1:500,000 scale SAR imagery (one look direction) was analyzed visually, and the detected linears were compared to those recorded on a 1:250,000 scale geologic map. Eighty percent of the 4,170 km of mapped, geologic linears were detected with the SAR imagery. Moreover, nearly 6,900 km of unmapped linears were also detected. Of these, an estimated 90 percent could be observed on high altitude aerial photographs. The relationship between SAR image detection of linears and the different types of indicators (e.g., straight valleys or shorelines) is reported.
Non-linear dynamic analysis of anisotropic cylindrical shells
Lakis, A.A.; Selmane, A.; Toledano, A.
1996-12-01
A theory to predict the influence of geometric non-linearities on the natural frequencies of an empty anisotropic cylindrical shell is presented in this paper. It is a hybrid of finite element and classical thin shell theories. Sanders-Koiter non-linear and strain-displacement relations are used. Displacement functions are evaluated using linearized equations of motion. Modal coefficients are then obtained for these displacement functions. Expressions for the mass, linear and non-linear stiffness matrices are derived through the finite element method. The uncoupled equations are solved with the help of elliptic functions. The period and frequency variations are first determined as a function of shell amplitudes and then compared with the results in the literature.
Energy and environmental analysis of a linear concentrating photovoltaic system
NASA Astrophysics Data System (ADS)
Kerzmann, Tony
The world is facing an imminent energy supply crisis. In order to sustain and increase our energy supply in an environmentally-conscious manner, it is necessary to advance renewable technologies. Despite this urgency, however, it is paramount to consider the larger environmental effects associated with using renewable energy resources. This research is meant to better understand linear concentrating photovoltaics (LCPVs) from an engineering and environmental standpoint. In order to analyze the LCPV system, a simulation and life cycle assessment (LCA) were developed. The LCPV system serves two major purposes: it produces electricity, and waste heat is collected for heating use. There are three parts to the LCPV simulation. The first part simulates the multijunction cell output so as to calculate the temperature-dependent electricity generation. The second part simulates the cell cooling and waste heat recovery system using a model consisting of heat transfer and fluid flow equations. The waste heat recovery in the LCPV system was linked to a hot water storage system, which was also modeled. Coupling the waste heat recovery simulation and the hot water storage system gives an overall integrated system that is useful for system design, optimization, and acts as a stepping stone for future multijunction cell Photovoltaic/Thermal (PV/T) systems. Finally, all of the LCPV system components were coded in Engineering Equation Solver (EES) and were used in an energy analysis under actual weather and solar conditions for the Phoenix, AZ, region. The life cycle assessment for the LCPV system allowed for an environmental analysis of the system where areas of the highest environmental impact were pinpointed. While conducting the LCA research, each component of the system was analyzed from a resource extraction, production, and use standpoint. The collective production processes of each LCPV system component were gathered into a single inventory of materials and energy flows
Linearization of the full activated sludge model No 1 for interaction analysis.
Benhalla, Abdelhay; Houssou, Mohamed; Charif, Moussa
2010-08-01
This paper deals with the linearization of the full activated sludge model No 1 (ASM1) in the scope of interaction analysis. For consistency, the linearization procedure is developed and validated within the BSM1 simulation benchmark framework. It is based on reaction rate approximation by linear combinations of states. The linear rate models are identified and incorporated in the mass balance equations, yielding a linear locally equivalent to the ASM1 model. Linear models for anoxic and aerated compartments are proposed. It is observed that the presented models track very closely the nonlinear ASM1 responses to various influent data. The key feature of this linearization strategy is that the gotten linear version of the ASM1 model is linear time invariant (LTI) and that it conserves the states biological interpretation and the original ASM1 dimension. It allows, therefore, application of interaction analysis methods and makes it possible to determine motivated control configurations for the ASM1 model. PMID:20131068
Analysis of non-linearity in differential wavefront sensing technique.
Duan, Hui-Zong; Liang, Yu-Rong; Yeh, Hsien-Chi
2016-03-01
An analytical model of a differential wavefront sensing (DWS) technique based on Gaussian Beam propagation has been derived. Compared with the result of the interference signals detected by quadrant photodiode, which is calculated by using the numerical method, the analytical model has been verified. Both the analytical model and numerical simulation show milli-radians level non-linearity effect of DWS detection. In addition, the beam clipping has strong influence on the non-linearity of DWS. The larger the beam clipping is, the smaller the non-linearity is. However, the beam walking effect hardly has influence on DWS. Thus, it can be ignored in laser interferometer. PMID:26974079
LCFIPlus: A framework for jet analysis in linear collider studies
NASA Astrophysics Data System (ADS)
Suehara, Taikan; Tanabe, Tomohiko
2016-02-01
We report on the progress in flavor identification tools developed for a future e+e- linear collider such as the International Linear Collider (ILC) and Compact Linear Collider (CLIC). Building on the work carried out by the LCFIVertex collaboration, we employ new strategies in vertex finding and jet finding, and introduce new discriminating variables for jet flavor identification. We present the performance of the new algorithms in the conditions simulated using a detector concept designed for the ILC. The algorithms have been successfully used in ILC physics simulation studies, such as those presented in the ILC Technical Design Report.
Analysis of ERTS-1 linear features in New York State
NASA Technical Reports Server (NTRS)
Isachsen, Y. W. (Principal Investigator); Fakundiny, R. H.; Forster, S. W.
1974-01-01
The author has identified the following significant results. All ERTS-1 linears confirmed to date have topographic expression although they may appear as featureless tonal linears on the imagery. A bias is unavoidably introduced against any linears which may parallel raster lines, lithological trends, or the azimuth of solar illumination. Ground study of ERTS-1 topographic lineaments in the Adirondacks indicates: outcrops along linears are even more rare than expected, fault breccias are found along some NNE lineaments, chloritization and slickensiding without brecciation characterize one EW lineament whereas closely-spaced jointing plus a zone of plastic shear define another. Field work in the Catskills suggests that the prominent new NNE lineaments may be surface manifestations of normal faulting in the basement, and that it may become possible to map major joint sets over extensive plateau regions directly on the imagery. Fall and winter images each display some unique linears, and long linears on the fall image commonly appear as aligned segments on the winter scene. A computer-processed color composite image permitted the extraction or additional information on the shaded side of mountains.
NASA Technical Reports Server (NTRS)
Jamison, J. W.
1994-01-01
CFORM was developed by the Kennedy Space Center Robotics Lab to assist in linear control system design and analysis using closed form and transient response mechanisms. The program computes the closed form solution and transient response of a linear (constant coefficient) differential equation. CFORM allows a choice of three input functions: the Unit Step (a unit change in displacement); the Ramp function (step velocity); and the Parabolic function (step acceleration). It is only accurate in cases where the differential equation has distinct roots, and does not handle the case for roots at the origin (s=0). Initial conditions must be zero. Differential equations may be input to CFORM in two forms - polynomial and product of factors. In some linear control analyses, it may be more appropriate to use a related program, Linear Control System Design and Analysis (KSC-11376), which uses root locus and frequency response methods. CFORM was written in VAX FORTRAN for a VAX 11/780 under VAX VMS 4.7. It has a central memory requirement of 30K. CFORM was developed in 1987.
Linearized analysis of Richtmyer-Meshkov flow for elastic materials
NASA Astrophysics Data System (ADS)
Plohr, Jeeyeon N.; Plohr, Bradley J.
2005-08-01
We present a study of Richtmyer Meshkov flow for elastic materials. This flow, in which a material interface is struck by a shock wave, was originally investigated for gases, where growth of perturbations of the interface is observed. Here we consider two elastic materials in frictionless contact. The governing system of equations comprises conservation laws supplemented by constitutive equations. To analyse it, we linearize the equations around a one-dimensional background solution under the assumption that the perturbation is small. The background problem defines a Riemann problem that is solved numerically; its solution contains transmitted and reflected shock waves in the longitudinal modes. The linearized Rankine Hugoniot condition provides the interface conditions at the longitudinal and shear waves; the frictionless material interface conditions are also linearized. The resulting equations, a linear system of partial differential equations, is solved numerically using a finite-difference method supplemented by front tracking. In verifying the numerical code, we reproduce growth of the interface in the gas case. For the elastic case, in contrast, we find that the material interface remains bounded: the non-zero shear stiffness stabilizes the flow. In particular, the linear theory remains valid at late time. Moreover, we identify the principal mechanism for the stability of Richtmyer Meshkov flow for elastic materials: the vorticity deposited on the material interface during shock passage is propagated away by the shear waves, whereas for gas dynamics it stays on the interface.
Non-linear dynamic analysis of geared systems, part 2
NASA Technical Reports Server (NTRS)
Singh, Rajendra; Houser, Donald R.; Kahraman, Ahmet
1990-01-01
A good understanding of the steady state dynamic behavior of a geared system is required in order to design reliable and quiet transmissions. This study focuses on a system containing a spur gear pair with backlash and periodically time-varying mesh stiffness, and rolling element bearings with clearance type non-linearities. A dynamic finite element model of the linear time-invariant (LTI) system is developed. Effects of several system parameters, such as torsional and transverse flexibilities of the shafts and prime mover/load inertias, on free and force vibration characteristics are investigated. Several reduced order LTI models are developed and validated by comparing their eigen solution with the finite element model results. Several key system parameters such as mean load and damping ratio are identified and their effects on the non-linear frequency response are evaluated quantitatively. Other fundamental issues such as the dynamic coupling between non-linear modes, dynamic interactions between component non-linearities and time-varying mesh stiffness, and the existence of subharmonic and chaotic solutions including routes to chaos have also been examined in depth.
An implementation analysis of the linear discontinuous finite element method
Becker, T. L.
2013-07-01
This paper provides an implementation analysis of the linear discontinuous finite element method (LD-FEM) that spans the space of (l, x, y, z). A practical implementation of LD includes 1) selecting a computationally efficient algorithm to solve the 4 x 4 matrix system Ax = b that describes the angular flux in a mesh element, and 2) choosing how to store the data used to construct the matrix A and the vector b to either reduce memory consumption or increase computational speed. To analyze the first of these, three algorithms were selected to solve the 4 x 4 matrix equation: Cramer's rule, a streamlined implementation of Gaussian elimination, and LAPACK's Gaussian elimination subroutine dgesv. The results indicate that Cramer's rule and the streamlined Gaussian elimination algorithm perform nearly equivalently and outperform LAPACK's implementation of Gaussian elimination by a factor of 2. To analyze the second implementation detail, three formulations of the discretized LD-FEM equations were provided for implementation in a transport solver: 1) a low-memory formulation, which relies heavily on 'on-the-fly' calculations and less on the storage of pre-computed data, 2) a high-memory formulation, which pre-computes much of the data used to construct A and b, and 3) a reduced-memory formulation, which lies between the low - and high-memory formulations. These three formulations were assessed in the Jaguar transport solver based on relative memory footprint and computational speed for increasing mesh size and quadrature order. The results indicated that the memory savings of the low-memory formulation were not sufficient to warrant its implementation. The high-memory formulation resulted in a significant speed advantage over the reduced-memory option (10-50%), but also resulted in a proportional increase in memory consumption (5-45%) for increasing quadrature order and mesh count; therefore, the practitioner should weigh the system memory constraints against any
Linear Analysis and Verification Suite for Edge Turbulence
Myra, J R; Umansky, M
2008-04-24
The edge and scrape-off-layer region of a tokamak plasma is subject to well known resistive and ideal instabilities that are driven by various curvature- and sheath-related mechanisms. While the boundary plasma is typically strongly turbulent in experiments, it is useful to have computational tools that can analyze the linear eigenmode structure, predict quantitative trends in growth rates and elucidate and the underlying drive mechanisms. Furthermore, measurement of the linear growth rate of unstable modes emerging from a known, established equilibrium configuration provides one of the few quantitative ways of rigorously benchmarking large-scale plasma turbulence codes with each other and with a universal standard. In this report, a suite of codes that can describe linearized, nonlocal (e.g. separatrix-spanning) modes in axisymmetric (realistic divertor), toroidal geometry is discussed. Examples of several benchmark comparisons are given, and future development plans for a new eigenvalue edge code are presented.
NASA Technical Reports Server (NTRS)
Nguyen, D. T.; Al-Nasra, M.; Zhang, Y.; Baddourah, M. A.; Agarwal, T. K.; Storaasli, O. O.; Carmona, E. A.
1991-01-01
Several parallel-vector computational improvements to the unconstrained optimization procedure are described which speed up the structural analysis-synthesis process. A fast parallel-vector Choleski-based equation solver, pvsolve, is incorporated into the well-known SAP-4 general-purpose finite-element code. The new code, denoted PV-SAP, is tested for static structural analysis. Initial results on a four processor CRAY 2 show that using pvsolve reduces the equation solution time by a factor of 14-16 over the original SAP-4 code. In addition, parallel-vector procedures for the Golden Block Search technique and the BFGS method are developed and tested for nonlinear unconstrained optimization. A parallel version of an iterative solver and the pvsolve direct solver are incorporated into the BFGS method. Preliminary results on nonlinear unconstrained optimization test problems, using pvsolve in the analysis, show excellent parallel-vector performance indicating that these parallel-vector algorithms can be used in a new generation of finite-element based structural design/analysis-synthesis codes.
Powerful tool for design analysis of linear control systems
Maddux, Jr, A S
1982-05-10
The methods for designing linear controls for electronic or mechanical systems have been understood and put to practice. What has not been readily available to engineers, however, is a practical, quick and inexpensive method for analyzing these linear control (feedback) systems once they have been designed into the electronic or mechanical hardware. Now, the PET, manufactured by Commodore Business Machines (CBM), operating with several peripherals via the IEEE 488 Bus, brings to the engineer for about $4000 a complete set of office tools for analyzing these system designs.
Linearized analysis for compressible laminar boundary layers by ECVL technique
NASA Astrophysics Data System (ADS)
Soong, Chyi-Yeou
1987-04-01
The effective convective velocity linearization (ECVL) technique of Soong (1984) is applied to simplify the equations for steady two-dimensional compressible laminar boundary-layer flow of a perfect gas with constant Prandtl number and specific heat. The derivation of the linearized momentum equation is explained, and numerical results for zero-pressure-gradient and pressure-gradient flows are presented graphically and compared with experimental measurements and analytical results obtained with other approaches. The accuracy and simplicity of the ECVL method are demonstrated, and its practical applicability is indicated.
Simulated Analysis of Linear Reversible Enzyme Inhibition with SCILAB
ERIC Educational Resources Information Center
Antuch, Manuel; Ramos, Yaquelin; Álvarez, Rubén
2014-01-01
SCILAB is a lesser-known program (than MATLAB) for numeric simulations and has the advantage of being free software. A challenging software-based activity to analyze the most common linear reversible inhibition types with SCILAB is described. Students establish typical values for the concentration of enzyme, substrate, and inhibitor to simulate…
Linear matrix inequalities for analysis and control of linear vector second-order systems
Adegas, Fabiano D.; Stoustrup, Jakob
2014-10-06
Many dynamical systems are modeled as vector second-order differential equations. This paper presents analysis and synthesis conditions in terms of LMI with explicit dependence in the coefficient matrices of vector second-order systems. These conditions benefit from the separation between the Lyapunov matrix and the system matrices by introducing matrix multipliers, which potentially reduce conservativeness in hard control problems. Multipliers facilitate the usage of parameter-dependent Lyapunov functions as certificates of stability of uncertain and time-varying vector second-order systems. The conditions introduced in this work have the potential to increase the practice of analyzing and controlling systems directly in vector second-order form.
NASA Technical Reports Server (NTRS)
Rizzi, Stephen A.
2003-01-01
The use of stress predictions from equivalent linearization analyses in the computation of high-cycle fatigue life is examined. Stresses so obtained differ in behavior from the fully nonlinear analysis in both spectral shape and amplitude. Consequently, fatigue life predictions made using this data will be affected. Comparisons of fatigue life predictions based upon the stress response obtained from equivalent linear and numerical simulation analyses are made to determine the range over which the equivalent linear analysis is applicable.
Linearized stability analysis of gravastars in noncommutative geometry
NASA Astrophysics Data System (ADS)
Lobo, Francisco S. N.; Garattini, Remo
2013-12-01
In this work, we find exact gravastar solutions in the context of noncommutative geometry, and explore their physical properties and characteristics. The energy density of these geometries is a smeared and particle-like gravitational source, where the mass is diffused throughout a region of linear dimension due to the intrinsic uncertainty encoded in the coordinate commutator. These solutions are then matched to an exterior Schwarzschild spacetime. We further explore the dynamical stability of the transition layer of these gravastars, for the specific case of β = M 2/ α < 1.9, where M is the black hole mass, to linearized spherically symmetric radial perturbations about static equilibrium solutions. It is found that large stability regions exist and, in particular, located sufficiently close to where the event horizon is expected to form.
Dielectric loss analysis using linear resonators with different impedances
NASA Astrophysics Data System (ADS)
Sarabi, Bahman; Khalil, M. S.; Khasawneh, M. A.; Stoutimore, M. J. A.; Gladchenko, Sergiy; Wellstood, F. C.; Lobb, C. J.; Osborn, K. D.
2012-02-01
It is known that amorphous dielectrics are a major source of decoherence in superconducting qubits due to energy absorption by two-level systems coupled to the electric fields. Linear resonators have been applied extensively to study loss in different dielectrics used in qubit circuits due to their versatility and relative simplicity in design, fabrication and measurement. We have designed linear resonators with multi-turn inductors and parallel-plate capacitors with resonance frequencies of 4.8-6.4 GHz. We achieve substantially different L/C values and capacitor volumes by varying the number of inductance turns and the area of the capacitors. We will present results of continuous wave measurements with SiNx capacitors and show how loss tangent and phase noise are related to impedance and capacitor volume.
HIGH RESOLUTION FOURIER ANALYSIS WITH AUTO-REGRESSIVE LINEAR PREDICTION
Barton, J.; Shirley, D.A.
1984-04-01
Auto-regressive linear prediction is adapted to double the resolution of Angle-Resolved Photoemission Extended Fine Structure (ARPEFS) Fourier transforms. Even with the optimal taper (weighting function), the commonly used taper-and-transform Fourier method has limited resolution: it assumes the signal is zero beyond the limits of the measurement. By seeking the Fourier spectrum of an infinite extent oscillation consistent with the measurements but otherwise having maximum entropy, the errors caused by finite data range can be reduced. Our procedure developed to implement this concept adapts auto-regressive linear prediction to extrapolate the signal in an effective and controllable manner. Difficulties encountered when processing actual ARPEFS data are discussed. A key feature of this approach is the ability to convert improved measurements (signal-to-noise or point density) into improved Fourier resolution.
Analysis of photo linear elements, Laramie Mountains, Wyoming
NASA Technical Reports Server (NTRS)
Blackstone, D. L., Jr.
1973-01-01
The author has identified the following significant results. Photo linear features in the Precambrian rocks of the Laramie Mountains are delineated, and the azimuths plotted on rose diagrams. Three strike directions are dominant, two of which are in the northeast quadrant. Laramide folds in the Laramie basin to the west of the mountains appear to have the same trend, and apparently have been controlled by response of the basement along fractures such as have been measured from the imagery.
Liapunov functions for non-linear difference equation stability analysis.
NASA Technical Reports Server (NTRS)
Park, K. E.; Kinnen, E.
1972-01-01
Liapunov functions to determine the stability of non-linear autonomous difference equations can be developed through the use of auxiliary exact difference equations. For this purpose definitions are introduced for the gradient of an implicit function of a discrete variable, a principal sum, a definite sum and an exact difference equation, and a theorem for exactness of a difference form is proved. Examples illustrate the procedure.
Linear and Nonlinear Analysis of Brain Dynamics in Children with Cerebral Palsy
ERIC Educational Resources Information Center
Sajedi, Firoozeh; Ahmadlou, Mehran; Vameghi, Roshanak; Gharib, Masoud; Hemmati, Sahel
2013-01-01
This study was carried out to determine linear and nonlinear changes of brain dynamics and their relationships with the motor dysfunctions in CP children. For this purpose power of EEG frequency bands (as a linear analysis) and EEG fractality (as a nonlinear analysis) were computed in eyes-closed resting state and statistically compared between 26…
A linear mixture analysis-based compression for hyperspectral image analysis
C. I. Chang; I. W. Ginsberg
2000-06-30
In this paper, the authors present a fully constrained least squares linear spectral mixture analysis-based compression technique for hyperspectral image analysis, particularly, target detection and classification. Unlike most compression techniques that directly deal with image gray levels, the proposed compression approach generates the abundance fractional images of potential targets present in an image scene and then encodes these fractional images so as to achieve data compression. Since the vital information used for image analysis is generally preserved and retained in the abundance fractional images, the loss of information may have very little impact on image analysis. In some occasions, it even improves analysis performance. Airborne visible infrared imaging spectrometer (AVIRIS) data experiments demonstrate that it can effectively detect and classify targets while achieving very high compression ratios.
Tangent linear analysis of the Mosaic land surface model
NASA Astrophysics Data System (ADS)
Yang, Runhua; Cohn, Stephen E.; da Silva, Arlindo; Joiner, Joanna; Houser, Paul R.
2003-01-01
In this study, a tangent linear eigenanalysis is applied to the Mosaic land surface model (LSM) [, 1992] to examine the impacts of the model internal dynamics and physics on the land surface state variability. The tangent linear model (TLM) of the Mosaic LSM is derived numerically for two sets of basic states and two tile types of land condition, grass and bare soil. An additional TLM, for the soil moisture subsystem of this LSM, is derived analytically for the same cases to obtain explicit expressions for the eigenvalues. An eigenvalue of the TLM determines a characteristic timescale, and the corresponding eigenvector, or mode, describes a particular coupling among the perturbed states. The results show that (1) errors in initial conditions tend to decay with e-folding times given by the characteristic timescales; (2) the LSM exhibits a wide range of internal variability, modes mainly representing surface temperature and surface moisture perturbations exhibit short timescales, whereas modes mainly representing deep soil temperature perturbations and moisture transfer throughout the entire soil column exhibit much longer timescales; (3) the modes of soil moisture tend to be weakly coupled with other perturbed variables, and the mode representing the deep soil temperature perturbation has a consistent e-folding time across the experiments; (4) the key parameters include soil moisture, soil layer depth, and soil hydraulic parameters. The results agree qualitatively with previous findings. However, tangent linear eigenanalysis provides a new approach to the quantitative substantiation of those findings. Also, it reveals the evolution and the coupling of the perturbed land states that are useful for the development of land surface data assimilation schemes. One must be careful when generalizing the quantitative results since they are obtained with respect to two specific basic states and two simple land conditions. Also, the methodology employed here does not apply
Analysis of integrity and microstructure of linear friction welded Waspaloy
Chamanfar, A.; Jahazi, M.; Gholipour, J.; Wanjara, P.; Yue, S.
2015-06-15
Nickel-base superalloy, Waspaloy, was linear friction welded (LFWed) under different axial shortening conditions of 2.0, 3.4, and 4.6 mm. The tensile properties and microhardness of the weldments were investigated in the as-LFWed condition and compared with those in the post-weld heat treated (PWHTed) condition. Mechanical properties were related to microstructures following examination by optical microscopy, high resolution scanning electron microscopy, and electron backscatter diffraction (EBSD). Analyses of the EBSD results in terms of the misorientation angle distribution, which represents the stored energy, were performed. In the as-LFWed condition, the yield strength (YS) and ultimate tensile strength (UTS) increased with axial shortening due to greater expulsion of the softened interfacial material toward the periphery as flash. In contrast, with increasing axial shortening the total elongation initially remained constant and then decreased. This was also related to the expulsion of the softened interfacial material into the bifurcated flash. Extensive dissolution of the strengthening phase (γ′) in the weld area during linear friction welding (LFW) contributed to the lower YS and UTS in the as-welded condition compared to the PWHTed condition where the γ′ particles were recovered. After performing post-weld heat treatment (PWHT), the total elongation improved due to the relaxation of stored energy and grain growth in the thermomechanically affected zone (TMAZ). - Highlights: • Tensile property and microstructure in Waspaloy linear friction welds were studied. • Yield strength and ultimate tensile strength increased with axial shortening. • Elongation initially remained constant and then decreased with axial shortening. • Post-weld heat treat recovered dissolved γ′ particles and increased weld strength. • Stored energy relaxation during post-weld heat treatment improved weld elongation.
Unitary discrete linear canonical transform: analysis and application.
Zhao, Liang; Healy, John J; Sheridan, John T
2013-03-01
The numerical approximation of the linear canonical transforms (LCTs) is important in modeling coherent wave field propagation through first-order optical systems and in many digital signal processing applications. The continuous LCTs are unitary, but discretization can destroy this property. We present a sufficient condition on the sampling rates chosen in the discretization to ensure unitarity. We discuss the various subsets of the unitary matrices examined in this paper that have been proposed elsewhere. We offer a proof of the existence of all of the unitary matrices we discuss. We examine the consequences of these results, particularly in relation to the use of discrete transforms in iterative phase retrieval applications. PMID:23458814
Analysis and application of minimum variance discrete linear system identification
NASA Technical Reports Server (NTRS)
Kotob, S.; Kaufman, H.
1977-01-01
An on-line minimum variance (MV) parameter identifier is developed which embodies both accuracy and computational efficiency. The formulation results in a linear estimation problem with both additive and multiplicative noise (AMN). The resulting filter which utilizes both the covariance of the parameter vector itself and the covariance of the error in identification is proven to be mean-square convergent and mean-square consistent. The MV parameter identification scheme is then used to construct a stable state and parameter estimation algorithm.
Structure/load dependent vectors for linear structural dynamic analysis
NASA Technical Reports Server (NTRS)
Qin, Jiangning; Nguyen, Duc T.
1992-01-01
The dynamic solution vectors yielded by the present structure/load dependent-vectors method for large-scale linear structural dynamic analyses involving complex loadings can be used as starting vectors, so that both structure and load characteristics are encompassed by the basis vectors. The method is shown to entail fewer vectors than current alternatives for a given level of accuracy, especially in the cases of structures that have external concentrated masses. Numerical results are presented which illustrate the advantages of this dependent-vectors method relative to other reduction methods.
Discriminative analysis of non-linear brain connectivity in schizophrenia: an fMRI Study
Su, Longfei; Wang, Lubin; Shen, Hui; Feng, Guiyu; Hu, Dewen
2013-01-01
Background: Dysfunctional integration of distributed brain networks is believed to be the cause of schizophrenia, and resting-state functional connectivity analyses of schizophrenia have attracted considerable attention in recent years. Unfortunately, existing functional connectivity analyses of schizophrenia have been mostly limited to linear associations. Objective: The objective of the present study is to evaluate the discriminative power of non-linear functional connectivity and identify its changes in schizophrenia. Method: A novel measure utilizing the extended maximal information coefficient was introduced to construct non-linear functional connectivity. In conjunction with multivariate pattern analysis, the new functional connectivity successfully discriminated schizophrenic patients from healthy controls with relative higher accuracy rate than the linear measure. Result: We found that the strength of the identified non-linear functional connections involved in the classification increased in patients with schizophrenia, which was opposed to its linear counterpart. Further functional network analysis revealed that the changes of the non-linear and linear connectivity have similar but not completely the same spatial distribution in human brain. Conclusion: The classification results suggest that the non-linear functional connectivity provided useful discriminative power in diagnosis of schizophrenia, and the inverse but similar spatial distributed changes between the non-linear and linear measure may indicate the underlying compensatory mechanism and the complex neuronal synchronization underlying the symptom of schizophrenia. PMID:24155713
NASA Technical Reports Server (NTRS)
Snow, L. S.; Kuhn, A. E.
1975-01-01
Previous error analyses conducted by the Guidance and Dynamics Branch of NASA have used the Guidance Analysis Program (GAP) as the trajectory simulation tool. Plans are made to conduct all future error analyses using the Space Vehicle Dynamics Simulation (SVDS) program. A study was conducted to compare the inertial measurement unit (IMU) error simulations of the two programs. Results of the GAP/SVDS comparison are presented and problem areas encountered while attempting to simulate IMU errors, vehicle performance uncertainties and environmental uncertainties using SVDS are defined. An evaluation of the SVDS linear error analysis capability is also included.
Linear analysis of ion cyclotron interaction in a multicomponent plasma
NASA Technical Reports Server (NTRS)
Gendrin, R.; Ashour-Abdalla, M.; Omura, Y.; Quest, K.
1984-01-01
The mechanism by which hot anisotropic protons generate electromagnetic ion cyclotron waves in a plasma containing cold H(+) and He(+) ions is quantitatively studied. Linear growth rates (both temporal and spatial) are computed for different plasma parameters: concentration, temperature,and anisotropy of cold He(+) ions and of hot protons. It is shown that: (1) for parameters typical of the geostationary altitude the maximum growth rates are not drastically changed when a small proportion (about 1 to 20 percent) of cold He(+) ions is present; (2) because of the important cyclotron absorption by thermal He(+) ions in the vicinity of the He(+) gyrofrequency, waves which could resonate with the bulk of the He(+) distribution cannot be generated. Therefore quasi-linear effects, in a homogeneous medium at least, cannot be responsible for the heating of He(+) ions which is often observed in conjunction with ion cyclotron waves. The variation of growth rate versus wave number is also studied for its importance in selecting suitable parameters in numerical simulation experiments.
A quasi-linear control theory analysis of timesharing skills
NASA Technical Reports Server (NTRS)
Agarwal, G. C.; Gottlieb, G. L.
1977-01-01
The compliance of the human ankle joint is measured by applying 0 to 50 Hz band-limited gaussian random torques to the foot of a seated human subject. These torques rotate the foot in a plantar-dorsal direction about a horizontal axis at a medial moleolus of the ankle. The applied torques and the resulting angular rotation of the foot are measured, digitized and recorded for off-line processing. Using such a best-fit, second-order model, the effective moment of inertia of the ankle joint, the angular viscosity and the stiffness are calculated. The ankle joint stiffness is shown to be a linear function of the level of tonic muscle contraction, increasing at a rate of 20 to 40 Nm/rad/Kg.m. of active torque. In terms of the muscle physiology, the more muscle fibers that are active, the greater the muscle stiffness. Joint viscosity also increases with activation. Joint stiffness is also a linear function of the joint angle, increasing at a rate of about 0.7 to 1.1 Nm/rad/deg from plantar flexion to dorsiflexion rotation.
Design and Analysis of Tubular Permanent Magnet Linear Wave Generator
Si, Jikai; Feng, Haichao; Su, Peng; Zhang, Lufeng
2014-01-01
Due to the lack of mature design program for the tubular permanent magnet linear wave generator (TPMLWG) and poor sinusoidal characteristics of the air gap flux density for the traditional surface-mounted TPMLWG, a design method and a new secondary structure of TPMLWG are proposed. An equivalent mathematical model of TPMLWG is established to adopt the transformation relationship between the linear velocity of permanent magnet rotary generator and the operating speed of TPMLWG, to determine the structure parameters of the TPMLWG. The new secondary structure of the TPMLWG contains surface-mounted permanent magnets and the interior permanent magnets, which form a series-parallel hybrid magnetic circuit, and their reasonable structure parameters are designed to get the optimum pole-arc coefficient. The electromagnetic field and temperature field of TPMLWG are analyzed using finite element method. It can be included that the sinusoidal characteristics of air gap flux density of the new secondary structure TPMLWG are improved, the cogging force as well as mechanical vibration is reduced in the process of operation, and the stable temperature rise of generator meets the design requirements when adopting the new secondary structure of the TPMLWG. PMID:25050388
A MODFLOW package to linearize stream depletion analysis
NASA Astrophysics Data System (ADS)
Ou, Gengxin; Li, Ruopu; Pun, Mahesh; Osborn, Colby; Bradley, Jesse; Schneider, Jim; Chen, Xun-Hong
2016-01-01
The conventional numerical method is computationally intensive and prone to numerical noises for stream depletion analyses using MODFLOW. In this study, a new MODFLOW package has been developed to improve the computational efficiency and reduce the noises for each simulation. Using the assumption of unchanged flow coefficients between the baseline and scenario runs, the nonlinear groundwater flow system is linearized for solving the flow equations. The new package has been successfully applied to a regional groundwater model in Nebraska. The results show that the numerical noises, commonly identified in conventional approach, have been significantly reduced and a twenty-fold speedup has been achieved for a regional groundwater model in Nebraska. The results suggest this package can be adapted to be a component of optimization tools for water management scenario analyses especially when a large number of scenario model runs are involved.
A triangular thin shell finite element: Linear analysis
NASA Technical Reports Server (NTRS)
Thomas, G. R.; Gallagher, R. H.
1975-01-01
The formulation of the linear stiffness matrix for a doubly-curved triangular thin shell element, using a modified potential energy principle, is described. The strain energy component of the potential energy is expressed in terms of displacements and displacement gradients by use of consistent Koiter strain-displacement equations. The element inplane and normal displacement fields are approximated by complete cubic polynomials. The interelement displacement admissibility conditions are met in the global representation by imposition of constraint conditions on the interelement boundaries; the constraints represent the modification of the potential energy. Errors due to the nonzero strains under rigid body motion are shown to be of small importance for practical grid refinements through performance of extensive comparison analyses.
Computational aspects of sensitivity calculations in linear transient structural analysis
NASA Technical Reports Server (NTRS)
Greene, W. H.; Haftka, R. T.
1991-01-01
The calculation of sensitivities of displacements, velocities, accelerations, and stresses in linear, structural, and transient response problems is studied. Several existing sensitivity calculation methods and two new methods are compared for three example problems. Approximation vectors such as vibration mode shapes are used to reduce the dimensionality of the finite model. To be computationally practical in large-order problems, the overall finite difference methods must use the approximation vectors from the original design in the analyses of the perturbed models. This was found to result in poor convergence of stress sensitivities in several cases. Two semianalytical techniques are developed to overcome this poor convergence. Both new methods result in very good convergence of the stress sensitivities; the computational cost is much less than would result if the vibration modes were recalculated and then used in an overall finite difference method.
Linear Stability Regime Analysis of the Compressible Reacting Mixing Layer
NASA Technical Reports Server (NTRS)
Day, M. J.; Reynolds, William C.; Mansour, N. N.; Rai, Man Mohan (Technical Monitor)
1995-01-01
Previous investigations have shown that a compressible reacting mixing layer can develop two peaks in the mean density weighted vorticity profile. Linear stability analyses show that at these peaks two distinct 'outer' instability modes appear in addition to the more common central mode, which exists unaccompanied in incompressible nonreacting flows. The present study parametrically analyzes the effects of compressibility, heat release, stoichiometry, and density ratio on the amplification rate and obliquity of each stability mode. The mean profiles used in the spatial stability calculation are generated by self-similar solutions of the compressible boundary layer equations combined with the assumption of infinitely fast chemistry. It is shown that the influence of stoichiometry and density ratio on the peaks of the density weighted vorticity profile determines which modes will dominate. Of particular interest are the conditions where two modes are equally amplified, causing the mixing layer to develop into a 'colayer' structure.
Linear MHD Stability Analysis of the SSPX Spheromak
NASA Astrophysics Data System (ADS)
Jayakumar, R.; Cohen, B. I.; Hooper, E. B.; Lodestro, L. L.; McLean, H. S.; Pearlstein, L. D.; Wood, R.; Turnbull, A. D.; Sovinec, C.
2007-11-01
Good correlation between the toroidal mode numbers of measured magnetic fluctuations in high temperature SSPX plasmas and presence of low-order rational surfaces in the reconstructed q profiles, suggests that the quality of magnetic surfaces in SSPX is sufficiently good for applying standard linear MHD stability analyses. Previously we have reported on benchmarking the code NIMROD against GATO, with good agreement in growth rates for ideal-MHD internal kinks and an external kinks with no current on open field lines (for equilibria imported from the code Corsica). Recent stability analyses also show that presence of low order rational surfaces causes internal modes to become unstable. We will report on the progress in applying these tools for assessing beta limits in SSPX, using NIMROD analyses including current on open field lines and for comparison with experiments.
NASA Astrophysics Data System (ADS)
Warmenbol, Eugène
2004-11-01
Gold remains a rare material during the whole Bronze Age. It is a metal that circulates amongst the elites exclusively. The exchange of gold artefacts over very long distances is fundamental in the build-up of socio-political relations. Fifty artefacts found during subaquatic excavations in the river Lesse in the cave of Han-sur-Lesse have been submitted to quantitative elemental analysis using PIXE, a method never used before on archaeological material found in Belgium. The results show that there are clearly three groups of artefacts. The first one (discs, "baskets", pearls) proves of local manufacture, as also suggested by the discovery of two ingots. The second one (gilded rings) could have been locally made too, but is a lot less homogeneous. The third one (filigreed and granulated beads or pin's heads) must be imported and most probably originates in the Mediterranean.
Sensitivity Analysis of Parameters in Linear-Quadratic Radiobiologic Modeling
Fowler, Jack F.
2009-04-01
Purpose: Radiobiologic modeling is increasingly used to estimate the effects of altered treatment plans, especially for dose escalation. The present article shows how much the linear-quadratic (LQ) (calculated biologically equivalent dose [BED] varies when individual parameters of the LQ formula are varied by {+-}20% and by 1%. Methods: Equivalent total doses (EQD2 = normalized total doses (NTD) in 2-Gy fractions for tumor control, acute mucosal reactions, and late complications were calculated using the linear- quadratic formula with overall time: BED = nd (1 + d/ [{alpha}/{beta}]) - log{sub e}2 (T - Tk) / {alpha}Tp, where BED is BED = total dose x relative effectiveness (RE = nd (1 + d/ [{alpha}/{beta}]). Each of the five biologic parameters in turn was altered by {+-}10%, and the altered EQD2s tabulated; the difference was finally divided by 20. EQD2 or NTD is obtained by dividing BED by the RE for 2-Gy fractions, using the appropriate {alpha}/{beta} ratio. Results: Variations in tumor and acute mucosal EQD ranged from 0.1% to 0.45% per 1% change in each parameter for conventional schedules, the largest variation being caused by overall time. Variations in 'late' EQD were 0.4% to 0.6% per 1% change in the only biologic parameter, the {alpha}/{beta} ratio. For stereotactic body radiotherapy schedules, variations were larger, up to 0.6 to 0.9 for tumor and 1.6% to 1.9% for late, per 1% change in parameter. Conclusions: Robustness occurs similar to that of equivalent uniform dose (EUD), for the same reasons. Total dose, dose per fraction, and dose-rate cause their major effects, as well known.
Three-dimensional linear system analysis for breast tomosynthesis.
Zhao, Bo; Zhao, Wei
2008-12-01
The optimization of digital breast tomosynthesis (DBT) geometry and reconstruction is crucial for the clinical translation of this exciting new imaging technique. In the present work, the authors developed a three-dimensional (3D) cascaded linear system model for DBT to investigate the effects of detector performance, imaging geometry, and image reconstruction algorithm on the reconstructed image quality. The characteristics of a prototype DBT system equipped with an amorphous selenium flat-panel detector and filtered backprojection reconstruction were used as an example in the implementation of the linear system model. The propagation of signal and noise in the frequency domain was divided into six cascaded stages incorporating the detector performance, imaging geometry, and reconstruction filters. The reconstructed tomosynthesis imaging quality was characterized by spatial frequency dependent presampling modulation transfer function (MTF), noise power spectrum (NPS), and detective quantum efficiency (DQE) in 3D. The results showed that both MTF and NPS were affected by the angular range of the tomosynthesis scan and the reconstruction filters. For image planes parallel to the detector (in-plane), MTF at low frequencies was improved with increase in angular range. The shape of the NPS was affected by the reconstruction filters. Noise aliasing in 3D could be introduced by insufficient voxel sampling, especially in the z (slice-thickness) direction where the sampling distance (slice thickness) could be more than ten times that for in-plane images. Aliasing increases the noise at high frequencies, which causes degradation in DQE. Application of a reconstruction filter that limits the frequency components beyond the Nyquist frequency in the z direction, referred to as the slice thickness filter, eliminates noise aliasing and improves 3D DQE. The focal spot blur, which arises from continuous tube travel during tomosynthesis acquisition, could degrade DQE significantly
Microstructural Analysis of Linear Friction-Welded 718 Plus Superalloy
NASA Astrophysics Data System (ADS)
Vishwakarma, K. R.; Ojo, O. A.; Wanjara, P.; Chaturvedi, M. C.
2014-12-01
The microstructure of Allvac 718 Plus (718 Plus; ATI, Pittsburgh, PA, USA) superalloy was examined after linear friction welding (LFW) and after standard postweld heat treatment (PWHT). The liquid phase reaction of second-phase precipitates, which are known to constitutionally liquate during conventional fusion welding, was observed in the thermomechanically affected zone (TMAZ) of the welded material. These phases included MC-type carbides, Ti-rich carbonitrides, and δ phase precipitates. This observation is contrary to the general assumption that LFW is a completely solid-state joining process. However, unlike conventional fusion welding processes that cause heat-affected zone liquation cracking in 718 Plus and many other superalloys, the LFW process did not cause cracking in 718 Plus superalloy despite the liquation of precipitates. This absence of cracking during joining is attributed to the applied compressive stress during the forging stage of the LFW process. Also, no cracking was observed after PWHT, although PWHT resulted in a microstructure that had a nonhomogeneous distribution of precipitates in the weld and the TMAZ.
Hamiltonian analysis for linearly acceleration-dependent Lagrangians
NASA Astrophysics Data System (ADS)
Cruz, Miguel; Gómez-Cortés, Rosario; Molgado, Alberto; Rojas, Efraín
2016-06-01
We study the constrained Ostrogradski-Hamilton framework for the equations of motion provided by mechanical systems described by second-order derivative actions with a linear dependence in the accelerations. We stress out the peculiar features provided by the surface terms arising for this type of theories and we discuss some important properties for this kind of actions in order to pave the way for the construction of a well defined quantum counterpart by means of canonical methods. In particular, we analyse in detail the constraint structure for these theories and its relation to the inherent conserved quantities where the associated energies together with a Noether charge may be identified. The constraint structure is fully analyzed without the introduction of auxiliary variables, as proposed in recent works involving higher order Lagrangians. Finally, we also provide some examples where our approach is explicitly applied and emphasize the way in which our original arrangement results in propitious for the Hamiltonian formulation of covariant field theories.
ERIC Educational Resources Information Center
Preacher, Kristopher J.; Curran, Patrick J.; Bauer, Daniel J.
2006-01-01
Simple slopes, regions of significance, and confidence bands are commonly used to evaluate interactions in multiple linear regression (MLR) models, and the use of these techniques has recently been extended to multilevel or hierarchical linear modeling (HLM) and latent curve analysis (LCA). However, conducting these tests and plotting the…
Application of Local Linear Embedding to Nonlinear Exploratory Latent Structure Analysis
ERIC Educational Resources Information Center
Wang, Haonan; Iyer, Hari
2007-01-01
In this paper we discuss the use of a recent dimension reduction technique called Locally Linear Embedding, introduced by Roweis and Saul, for performing an exploratory latent structure analysis. The coordinate variables from the locally linear embedding describing the manifold on which the data reside serve as the latent variable scores. We…
Simplified seawater alkalinity analysis: Use of linear array spectrometers
NASA Astrophysics Data System (ADS)
Yao, Wensheng; Byrne, Robert H.
1998-08-01
Modified spectrophotometric procedures are presented for the determination of seawater total alkalinity using rapid scan linear array spectrometers. Continuous monitoring of solution pH allows titrations to be terminated at relatively high pH, whereby excess acid terms are very small. Excess acid concentrations are quantified using the sulfonephthalein indicators, bromocresol green and bromocresol purple. The outlined spectrophotometric procedures require no thermal equilibration of samples. Using bromocresol green, solution pH T ([H +] T in moles per kg of solution) is given as: pHT=4.2699+0.002578(35- S)+ log((R(25)-0.00131)/(2.3148-0.1299 R(25))) - log(1-0.001005S) and R(25)= R( t){1+0.00909(25- t)}, where 29⩽S⩽37, 13° C⩽t⩽32° C, and R( t) is the absorbance ratio ( A616/ A444) at temperature t and salinity S. Using bromocresol purple, the solution pH T is given as pH T=5.8182+0.00129(35- S)+log(( R(25)-0.00381)/(2.8729-0.05104 R(25))) and R(25)= R( t){1+0.01869(25- t)}, where 29⩽S⩽37, 13° C⩽t⩽32° C, and R( t)= A589/ A432. Alkalinity measurements using bromocresol purple had a precision on the order of 0.3 μmol kg -1 and were within 0.3-0.9 μmol kg -1 of the alkalinities of certified seawater reference materials.
The Effect of Data Scaling on Dual Prices and Sensitivity Analysis in Linear Programs
ERIC Educational Resources Information Center
Adlakha, V. G.; Vemuganti, R. R.
2007-01-01
In many practical situations scaling the data is necessary to solve linear programs. This note explores the relationships in translating the sensitivity analysis between the original and the scaled problems.
NASA Astrophysics Data System (ADS)
Goyal, Deepak
Textile composites have a wide variety of applications in the aerospace, sports, automobile, marine and medical industries. Due to the availability of a variety of textile architectures and numerous parameters associated with each, optimal design through extensive experimental testing is not practical. Predictive tools are needed to perform virtual experiments of various options. The focus of this research is to develop a better understanding of linear elastic response, plasticity and material damage induced nonlinear behavior and mechanics of load flow in textile composites. Textile composites exhibit multiple scales of complexity. The various textile behaviors are analyzed using a two-scale finite element modeling. A framework to allow use of a wide variety of damage initiation and growth models is proposed. Plasticity induced non-linear behavior of 2x2 braided composites is investigated using a modeling approach based on Hill's yield function for orthotropic materials. The mechanics of load flow in textile composites is demonstrated using special non-standard postprocessing techniques that not only highlight the important details, but also transform the extensive amount of output data into comprehensible modes of behavior. The investigations show that the damage models differ from each other in terms of amount of degradation as well as the properties to be degraded under a particular failure mode. When compared with experimental data, predictions of some models match well for glass/epoxy composite whereas other's match well for carbon/epoxy composites. However, all the models predicted very similar response when damage factors were made similar, which shows that the magnitude of damage factors are very important. Full 3D as well as equivalent tape laminate predictions lie within the range of the experimental data for a wide variety of braided composites with different material systems, which validated the plasticity analysis. Conclusions about the effect of
Computer analysis of general linear networks using digraphs.
NASA Technical Reports Server (NTRS)
Mcclenahan, J. O.; Chan, S.-P.
1972-01-01
Investigation of the application of digraphs in analyzing general electronic networks, and development of a computer program based on a particular digraph method developed by Chen. The Chen digraph method is a topological method for solution of networks and serves as a shortcut when hand calculations are required. The advantage offered by this method of analysis is that the results are in symbolic form. It is limited, however, by the size of network that may be handled. Usually hand calculations become too tedious for networks larger than about five nodes, depending on how many elements the network contains. Direct determinant expansion for a five-node network is a very tedious process also.
Consistent linearization method for finite-element analysis of viscoelastic materials
Smith, P.D.; Pelessone, D.
1983-05-01
A method of formulating material models for viscoelastic analysis using the finite-element method is presented. The method, named consistent linearization, includes the influence of creep in the material stiffness in a theoretically ideal manner. This method has been applied to the linear viscoelastic analysis of graphite subject to irradiation. Previously, using the initial strain method, short time steps had been required to avoid a numerical instability associated with the rapid transient creep. Using the consistent linearization method a factor of 15 reduction in computer time was achieved for the same accuracy.
Linear error analysis of slope-area discharge determinations
Kirby, W.H.
1987-01-01
The slope-area method can be used to calculate peak flood discharges when current-meter measurements are not possible. This calculation depends on several quantities, such as water-surface fall, that are subject to large measurement errors. Other critical quantities, such as Manning's n, are not even amenable to direct measurement but can only be estimated. Finally, scour and fill may cause gross discrepancies between the observed condition of the channel and the hydraulic conditions during the flood peak. The effects of these potential errors on the accuracy of the computed discharge have been estimated by statistical error analysis using a Taylor-series approximation of the discharge formula and the well-known formula for the variance of a sum of correlated random variates. The resultant error variance of the computed discharge is a weighted sum of covariances of the various observational errors. The weights depend on the hydraulic and geometric configuration of the channel. The mathematical analysis confirms the rule of thumb that relative errors in computed discharge increase rapidly when velocity heads exceed the water-surface fall, when the flow field is expanding and when lateral velocity variation (alpha) is large. It also confirms the extreme importance of accurately assessing the presence of scour or fill. ?? 1987.
Linear stability analysis of three-dimensional compressible boundary layers
NASA Technical Reports Server (NTRS)
Malik, Mujeeb R.; Orszag, Steven A.
1987-01-01
A compressible stability analysis computer code is developed. The code uses a matrix finite-difference method for local eigenvale solution when a good guess for the eigenvalue is available and is significantly more computationally efficient than the commonly used inital-value approach. The local eigenvalue search procedure also results in eigenfunctions and, at little extra work, group velocities. A globally convergent eigenvalue procedure is also developed that may be used when no guess for the eigenvalue is available. The global problem is formulated in such a way that no unstable spurious modes appear so that the method is suitable for use in a black-box stability code. Sample stability calculations are presented for the boundary layer profiles of an LFC swept wing.
An analysis method for control reconfigurability of linear systems
NASA Astrophysics Data System (ADS)
Wang, Dayi; Duan, Wenjie; Liu, Chengrui
2016-01-01
The reconfigurability of control systems is further researched based on the function-objective model (FOM). The establishment of the FOM has been published in the authors' former paper, solving the problem whether the system is reconfigurable without losing the desired control objective. Based on the FOM, the importance factor, the risk factor and the k th reconfigurability factor are proposed to evaluate the fault risks of all components and the system reconfigurability with k faults. These factors show which components should be improved and which faults cannot be tolerated. The analysis results are very useful for enhancing the fault-tolerance performances of the control systems by improving system designs. A satellite model is utilized to illustrate the proposed method.
The Stability of Radiatively Cooling Jets I. Linear Analysis
NASA Technical Reports Server (NTRS)
Hardee, Philip E.; Stone, James M.
1997-01-01
The results of a spatial stability analysis of a two-dimensional slab jet, in which optically thin radiative cooling is dynamically important, are presented. We study both magnetized and unmagnetized jets at external Mach numbers of 5 and 20. We model the cooling rate by using two different cooling curves: one appropriate to interstellar gas, and the other to photoionized gas of reduced metallicity. Thus, our results will be applicable to both protostellar (Herbig-Haro) jets and optical jets from active galactic nuclei. We present analytical solutions to the dispersion relations in useful limits and solve the dispersion relations numerically over a broad range of perturbation frequencies. We find that the growth rates and wavelengths of the unstable Kelvin-Helmholtz (K-H) modes are significantly different from the adiabatic limit, and that the form of the cooling function strongly affects the results. In particular, if the cooling curve is a steep function of temperature in the neighborhood of the equilibrium state, then the growth of K-H modes is reduced relative to the adiabatic jet. On the other hand, if the cooling curve is a shallow function of temperature, then the growth of K-H modes can be enhanced relative to the adiabatic jet by the increase in cooling relative to heating in overdense regions. Inclusion of a dynamically important magnetic field does not strongly modify the important differences between an adiabatic jet and a cooling jet, provided the jet is highly supermagnetosonic and not magnetic pressure-dominated. In the latter case, the unstable modes behave more like the transmagnetosonic magnetic pressure-dominated adiabatic limit. We also plot fluid displacement surfaces associated with the various waves in a cooling jet in order to predict the structures that might arise in the nonlinear regime. This analysis predicts that low-frequency surface waves and the lowest order body modes will be the most effective at producing observable features in
Biomedical image analysis using Markov random fields & efficient linear programing.
Komodakis, Nikos; Besbes, Ahmed; Glocker, Ben; Paragios, Nikos
2009-01-01
Computer-aided diagnosis through biomedical image analysis is increasingly considered in health sciences. This is due to the progress made on the acquisition side, as well as on the processing one. In vivo visualization of human tissues where one can determine both anatomical and functional information is now possible. The use of these images with efficient intelligent mathematical and processing tools allows the interpretation of the tissues state and facilitates the task of the physicians. Segmentation and registration are the two most fundamental tools in bioimaging. The first aims to provide automatic tools for organ delineation from images, while the second focuses on establishing correspondences between observations inter and intra subject and modalities. In this paper, we present some recent results towards a common formulation addressing these problems, called the Markov Random Fields. Such an approach is modular with respect to the application context, can be easily extended to deal with various modalities, provides guarantees on the optimality properties of the obtained solution and is computationally efficient. PMID:19963682
GraTeLPy: graph-theoretic linear stability analysis
2014-01-01
Background A biochemical mechanism with mass action kinetics can be represented as a directed bipartite graph (bipartite digraph), and modeled by a system of differential equations. If the differential equations (DE) model can give rise to some instability such as multistability or Turing instability, then the bipartite digraph contains a structure referred to as a critical fragment. In some cases the existence of a critical fragment indicates that the DE model can display oscillations for some parameter values. We have implemented a graph-theoretic method that identifies the critical fragments of the bipartite digraph of a biochemical mechanism. Results GraTeLPy lists all critical fragments of the bipartite digraph of a given biochemical mechanism, thus enabling a preliminary analysis on the potential of a biochemical mechanism for some instability based on its topological structure. The correctness of the implementation is supported by multiple examples. The code is implemented in Python, relies on open software, and is available under the GNU General Public License. Conclusions GraTeLPy can be used by researchers to test large biochemical mechanisms with mass action kinetics for their capacity for multistability, oscillations and Turing instability. PMID:24572152
Modal cost analysis for linear matrix-second-order systems
NASA Technical Reports Server (NTRS)
Skelton, R. E.; Hughes, P. C.
1980-01-01
Reduced models and reduced controllers for systems governed by matrix-second-order differential equations are obtained by retaining those modes which make the largest contributions to quadratic control objectives. Such contributions, expressed in terms of modal data, used as mode truncation criteria, allow the statement of the specific control objectives to influence the early model reduction from very high order models which are available, for example, from finite element methods. The relative importance of damping, frequency, and eigenvector in the mode truncation decisions are made explicit for each of these control objectives: attitude control, vibration suppression and figure control. The paper also shows that using modal cost analysis (MCA) on the closed loop modes of the optimally controlled system allows the construction of reduced control policies which feedback only those closed loop modal coordinates which are most critical to the quadratic control performance criterion. In this way, the modes which should be controlled (and hence the modes which must be observable by choice of measurements), are deduced from truncations of the optimal controller.
Bounded Linear Stability Analysis - A Time Delay Margin Estimation Approach for Adaptive Control
NASA Technical Reports Server (NTRS)
Nguyen, Nhan T.; Ishihara, Abraham K.; Krishnakumar, Kalmanje Srinlvas; Bakhtiari-Nejad, Maryam
2009-01-01
This paper presents a method for estimating time delay margin for model-reference adaptive control of systems with almost linear structured uncertainty. The bounded linear stability analysis method seeks to represent the conventional model-reference adaptive law by a locally bounded linear approximation within a small time window using the comparison lemma. The locally bounded linear approximation of the combined adaptive system is cast in a form of an input-time-delay differential equation over a small time window. The time delay margin of this system represents a local stability measure and is computed analytically by a matrix measure method, which provides a simple analytical technique for estimating an upper bound of time delay margin. Based on simulation results for a scalar model-reference adaptive control system, both the bounded linear stability method and the matrix measure method are seen to provide a reasonably accurate and yet not too conservative time delay margin estimation.
NASA Technical Reports Server (NTRS)
Cooke, C. H.
1975-01-01
STICAP (Stiff Circuit Analysis Program) is a FORTRAN 4 computer program written for the CDC-6400-6600 computer series and SCOPE 3.0 operating system. It provides the circuit analyst a tool for automatically computing the transient responses and frequency responses of large linear time invariant networks, both stiff and nonstiff (algorithms and numerical integration techniques are described). The circuit description and user's program input language is engineer-oriented, making simple the task of using the program. Engineering theories underlying STICAP are examined. A user's manual is included which explains user interaction with the program and gives results of typical circuit design applications. Also, the program structure from a systems programmer's viewpoint is depicted and flow charts and other software documentation are given.
Linear analysis of a coaxial-waveguide gyrotron traveling-wave tube
Hung, C.L.
2006-03-15
Linear theory provides an efficient analysis model for the preliminary design of a gyrotron traveling-wave tube (gyro-TWT). This study presents a linear theory, which is applicable to amplifications or self-excited oscillations induced by absolute instabilities in a coaxial waveguide of finite length. The effects of wall losses are incorporated in the theoretical formalism. The validity of the linear theory is verified by comparison with calculation results obtained using an existing self-consistent nonlinear theory. The linear theory is applied to analyze a TE{sub 01} mode coaxial gyro-TWT at the fundamental cyclotron harmonic. Numerical analysis of coupling between the beam cyclotron mode and the waveguide mode provides physical insight into the wave-growing mechanisms of various oscillations. The critical parameters for the onset of threatening oscillation modes are analyzed to determine the stable operating conditions. Finally, the dependencies of small-signal amplifications on system parameters are studied in great detail.
Numerical analysis of linear buckling of wind turbine blade with different trailing bonding models
NASA Astrophysics Data System (ADS)
Zhang, J. D.; Xu, Y.
2013-12-01
The work focus on the linear buckling analysis of wind turbine blade with different trailing bonding models. Based on finite element model, it has been demonstrated that there are some differences for buckling load factor between different models. Several different models are valid for buckling analysis.
Using Linear versus Quadratic Rules in Predictive and Descriptive Discriminant Analysis.
ERIC Educational Resources Information Center
McGee, Jennifer
Both predictive discriminant analysis (PDA) and descriptive discriminant analysis (DDA) require a decision to pool group covariance matrices, or alternatively, to retain separate group covariance matrices when the group covariance matrices are too dissimilar to pool. Pooling the group covariance matrices involves the so-called "linear" rule,…
Augmenting Visual Analysis in Single-Case Research with Hierarchical Linear Modeling
ERIC Educational Resources Information Center
Davis, Dawn H.; Gagne, Phill; Fredrick, Laura D.; Alberto, Paul A.; Waugh, Rebecca E.; Haardorfer, Regine
2013-01-01
The purpose of this article is to demonstrate how hierarchical linear modeling (HLM) can be used to enhance visual analysis of single-case research (SCR) designs. First, the authors demonstrated the use of growth modeling via HLM to augment visual analysis of a sophisticated single-case study. Data were used from a delayed multiple baseline…
Meta-Analysis in Higher Education: An Illustrative Example Using Hierarchical Linear Modeling
ERIC Educational Resources Information Center
Denson, Nida; Seltzer, Michael H.
2011-01-01
The purpose of this article is to provide higher education researchers with an illustrative example of meta-analysis utilizing hierarchical linear modeling (HLM). This article demonstrates the step-by-step process of meta-analysis using a recently-published study examining the effects of curricular and co-curricular diversity activities on racial…
Application of variational and Galerkin equations to linear and nonlinear finite element analysis
NASA Technical Reports Server (NTRS)
Yu, Y.-Y.
1974-01-01
The paper discusses the application of the variational equation to nonlinear finite element analysis. The problem of beam vibration with large deflection is considered. The variational equation is shown to be flexible in both the solution of a general problem and in the finite element formulation. Difficulties are shown to arise when Galerkin's equations are used in the consideration of the finite element formulation of two-dimensional linear elasticity and of the linear classical beam.
Polynomial elimination theory and non-linear stability analysis for the Euler equations
NASA Technical Reports Server (NTRS)
Kennon, S. R.; Dulikravich, G. S.; Jespersen, D. C.
1986-01-01
Numerical methods are presented that exploit the polynomial properties of discretizations of the Euler equations. It is noted that most finite difference or finite volume discretizations of the steady-state Euler equations produce a polynomial system of equations to be solved. These equations are solved using classical polynomial elimination theory, with some innovative modifications. This paper also presents some preliminary results of a new non-linear stability analysis technique. This technique is applicable to determining the stability of polynomial iterative schemes. Results are presented for applying the elimination technique to a one-dimensional test case. For this test case, the exact solution is computed in three iterations. The non-linear stability analysis is applied to determine the optimal time step for solving Burgers' equation using the MacCormack scheme. The estimated optimal time step is very close to the time step that arises from a linear stability analysis.
Garfinkle, David; Pretorius, Frans; Yunes, Nicolas
2010-08-15
We perform a linear stability analysis of dynamical Chern-Simons modified gravity in the geometric optics approximation and find that it is linearly stable on the backgrounds considered. Our analysis also reveals that gravitational waves in the modified theory travel at the speed of light in Minkowski spacetime. However, on a Schwarzschild background the characteristic speed of propagation along a given direction splits into two modes, one subluminal and one superluminal. The width of the splitting depends on the azimuthal components of the propagation vector, is linearly proportional to the mass of the black hole, and decreases with the third inverse power of the distance from the black hole. Radial propagation is unaffected, implying that as probed by gravitational waves the location of the event horizon of the spacetime is unaltered. The analysis further reveals that when a high frequency, pure gravitational wave is scattered from a black hole, a scalar wave of comparable amplitude is excited, and vice versa.
NASA Technical Reports Server (NTRS)
Reddy, T. S. R.; Srivastava, R.; Mehmed, Oral
2002-01-01
An aeroelastic analysis system for flutter and forced response analysis of turbomachines based on a two-dimensional linearized unsteady Euler solver has been developed. The ASTROP2 code, an aeroelastic stability analysis program for turbomachinery, was used as a basis for this development. The ASTROP2 code uses strip theory to couple a two dimensional aerodynamic model with a three dimensional structural model. The code was modified to include forced response capability. The formulation was also modified to include aeroelastic analysis with mistuning. A linearized unsteady Euler solver, LINFLX2D is added to model the unsteady aerodynamics in ASTROP2. By calculating the unsteady aerodynamic loads using LINFLX2D, it is possible to include the effects of transonic flow on flutter and forced response in the analysis. The stability is inferred from an eigenvalue analysis. The revised code, ASTROP2-LE for ASTROP2 code using Linearized Euler aerodynamics, is validated by comparing the predictions with those obtained using linear unsteady aerodynamic solutions.
Analysis of an inventory model for both linearly decreasing demand and holding cost
NASA Astrophysics Data System (ADS)
Malik, A. K.; Singh, Parth Raj; Tomar, Ajay; Kumar, Satish; Yadav, S. K.
2016-03-01
This study proposes the analysis of an inventory model for linearly decreasing demand and holding cost for non-instantaneous deteriorating items. The inventory model focuses on commodities having linearly decreasing demand without shortages. The holding cost doesn't remain uniform with time due to any form of variation in the time value of money. Here we consider that the holding cost decreases with respect to time. The optimal time interval for the total profit and the optimal order quantity are determined. The developed inventory model is pointed up through a numerical example. It also includes the sensitivity analysis.
NASA Technical Reports Server (NTRS)
Rankin, C. C.
1988-01-01
A consistent linearization is provided for the element-dependent corotational formulation, providing the proper first and second variation of the strain energy. As a result, the warping problem that has plagued flat elements has been overcome, with beneficial effects carried over to linear solutions. True Newton quadratic convergence has been restored to the Structural Analysis of General Shells (STAGS) code for conservative loading using the full corotational implementation. Some implications for general finite element analysis are discussed, including what effect the automatic frame invariance provided by this work might have on the development of new, improved elements.
Stability analysis, non-linear pulsations and mass loss of models for 55 Cygni (HD 198478)
NASA Astrophysics Data System (ADS)
Yadav, Abhay Pratap; Glatzel, Wolfgang
2016-04-01
55 Cygni is a variable supergiant. Recent observational studies revealed that this star pulsates in pressure, gravity and strange modes. The pulsations seem to be associated with episodes of mass loss. In this paper we present a theoretical study of stellar models with parameters close to that of 55 Cygni. A linear non-adiabatic stability analysis with respect to radial perturbations is performed and the evolution of instabilities into the non-linear regime is followed by numerical simulation. Our study indicates that the mass of 55 Cygni lies below 28 M⊙. As the final consequence of the instabilities the non-linear simulations revealed finite amplitude pulsations with periods consistent with the observations. The non-linear results also indicate a connection between pulsations and mass loss and allow for an estimate of the mean mass-loss rate. It is consistent with the observed values.
Tuning of power system controllers using symbolic eigensensitivity analysis and linear programming
Xu, L.; Ahmed-Zaid, S.
1995-02-01
In this paper, a systematic method for the tuning of multiple power system controllers using symbolic eigensensitivity analysis and linear programming is presented. The concept of eigenvalue sensitivity is used here to formulate the linear first-order variation of the real part of the dominant system eigenvalue as a symbolic function of controller parameters. This eigenvalue variation is minimized iteratively subject to linear equality and inequality constraints. All controller parameters are tuned correctly when the objective function meets the desired performance criteria. In this method, no special assumptions are made on the type of power devices and linear feedback controllers present in the system. The proposed method is illustrated with three examples of a single-machine system with a power system stabilizer and a controller for a thyristor-controlled series capacitor. The tested configurations are supplemented with nonlinear time-domain simulations validating the small-signal stability results.
Lyra, Wladimir
2014-07-01
Recently, Klahr and Hubbard claimed that a hydrodynamical linear overstability exists in protoplanetary disks, powered by buoyancy in the presence of thermal relaxation. We analyze this claim, confirming it through rigorous compressible linear analysis. We model the system numerically, reproducing the linear growth rate for all cases studied. We also study the saturated properties of the overstability in the shearing box, finding that the saturated state produces finite amplitude fluctuations strong enough to trigger the subcritical baroclinic instability (SBI). Saturation leads to a fast burst of enstrophy in the box, and a large-scale vortex develops in the course of the next ≈100 orbits. The amount of angular momentum transport achieved is of the order of α ≈ 10{sup –3}, as in compressible SBI models. For the first time, a self-sustained three-dimensional vortex is produced from linear amplitude perturbation of a quiescent base state.
Development of a linearized unsteady Euler analysis for turbomachinery blade rows
NASA Technical Reports Server (NTRS)
Verdon, Joseph M.; Montgomery, Matthew D.; Kousen, Kenneth A.
1995-01-01
A linearized unsteady aerodynamic analysis for axial-flow turbomachinery blading is described in this report. The linearization is based on the Euler equations of fluid motion and is motivated by the need for an efficient aerodynamic analysis that can be used in predicting the aeroelastic and aeroacoustic responses of blade rows. The field equations and surface conditions required for inviscid, nonlinear and linearized, unsteady aerodynamic analyses of three-dimensional flow through a single, blade row operating within a cylindrical duct, are derived. An existing numerical algorithm for determining time-accurate solutions of the nonlinear unsteady flow problem is described, and a numerical model, based upon this nonlinear flow solver, is formulated for the first-harmonic linear unsteady problem. The linearized aerodynamic and numerical models have been implemented into a first-harmonic unsteady flow code, called LINFLUX. At present this code applies only to two-dimensional flows, but an extension to three-dimensions is planned as future work. The three-dimensional aerodynamic and numerical formulations are described in this report. Numerical results for two-dimensional unsteady cascade flows, excited by prescribed blade motions and prescribed aerodynamic disturbances at inlet and exit, are also provided to illustrate the present capabilities of the LINFLUX analysis.
A Three-Dimensional Linearized Unsteady Euler Analysis for Turbomachinery Blade Rows
NASA Technical Reports Server (NTRS)
Montgomery, Matthew D.; Verdon, Joseph M.
1997-01-01
A three-dimensional, linearized, Euler analysis is being developed to provide an efficient unsteady aerodynamic analysis that can be used to predict the aeroelastic and aeroacoustic responses of axial-flow turbo-machinery blading.The field equations and boundary conditions needed to describe nonlinear and linearized inviscid unsteady flows through a blade row operating within a cylindrical annular duct are presented. A numerical model for linearized inviscid unsteady flows, which couples a near-field, implicit, wave-split, finite volume analysis to a far-field eigenanalysis, is also described. The linearized aerodynamic and numerical models have been implemented into a three-dimensional linearized unsteady flow code, called LINFLUX. This code has been applied to selected, benchmark, unsteady, subsonic flows to establish its accuracy and to demonstrate its current capabilities. The unsteady flows considered, have been chosen to allow convenient comparisons between the LINFLUX results and those of well-known, two-dimensional, unsteady flow codes. Detailed numerical results for a helical fan and a three-dimensional version of the 10th Standard Cascade indicate that important progress has been made towards the development of a reliable and useful, three-dimensional, prediction capability that can be used in aeroelastic and aeroacoustic design studies.
Development of a Linearized Unsteady Euler Analysis with Application to Wake/Blade-Row Interactions
NASA Technical Reports Server (NTRS)
Verdon, Joseph M.; Montgomery, Matthew D.; Chuang, H. Andrew
1999-01-01
A three-dimensional, linearized, Euler analysis is being developed to provide a comprehensive and efficient unsteady aerodynamic analysis for predicting the aeroacoustic and aeroelastic responses of axial-flow turbomachinery blading. The mathematical models needed to describe nonlinear and linearized, inviscid, unsteady flows through a blade row operating within a cylindrical annular duct are presented in this report. A numerical model for linearized inviscid unsteady flows, which couples a near-field, implicit, wave-split, finite volume analysis to far-field eigen analyses, is also described. The linearized aerodynamic and numerical models have been implemented into the three-dimensional unsteady flow code, LINFLUX. This code is applied herein to predict unsteady subsonic flows driven by wake or vortical excitations. The intent is to validate the LINFLUX analysis via numerical results for simple benchmark unsteady flows and to demonstrate this analysis via application to a realistic wake/blade-row interaction. Detailed numerical results for a three-dimensional version of the 10th Standard Cascade and a fan exit guide vane indicate that LINFLUX is becoming a reliable and useful unsteady aerodynamic prediction capability that can be applied, in the future, to assess the three-dimensional flow physics important to blade-row, aeroacoustic and aeroelastic responses.
NASA Astrophysics Data System (ADS)
Wang, Xu
2016-02-01
This paper establishes coupling loss factor of linear vibration energy harvesting systems in a framework of statistical energy analysis under parameter variations and random excitations. The new contributions of this paper are to define the numerical ranges of the dimensionless force factor for the weak, moderate and strong coupling and to study the connections of dimensionless force factor, coupling loss factor, coupling quotient, critical coupling strength, electro-mechanical coupling factor, damping loss factor and modal densities in linear vibration energy harvesting systems. The motivation of this paper is to enable statistical energy analysis of linear vibration energy harvesting systems for reliable performance predictions and design optimisation under parameter variations of materials and manufacturing processes and random ambient environmental excitations.
NASA Astrophysics Data System (ADS)
Rahamim, Rony; Mahlab, Uri; Dahan, David
2009-10-01
We propose and demonstrate a novel approach to identify linear and nonlinear propagation regimes of an optical signal in an optical fiber link by using chaos analysis. We show that the chaotic characteristics of a propagating optical signal are affected by both the chromatic dispersion and the nonlinear effects in the optical fiber. Linear or nonlinear behavior is detected by determining the maximum Lyapunov exponent of the signal and the use of the recurrence plot technique. An experimental demonstration is performed using 10-Gbps signal propagation in a 100-km fiber link with different launched optical powers. Chaos analysis shows a clear identification of the linear and nonlinear optical propagation regimes by using a classification scheme based on a multilayer neural network. Numerical simulations confirm the experimental results.
A Signal Transmission Technique for Stability Analysis of Multivariable Non-Linear Control Systems
NASA Technical Reports Server (NTRS)
Jackson, Mark; Zimpfer, Doug; Adams, Neil; Lindsey, K. L. (Technical Monitor)
2000-01-01
Among the difficulties associated with multivariable, non-linear control systems is the problem of assessing closed-loop stability. Of particular interest is the class of non-linear systems controlled with on/off actuators, such as spacecraft thrusters or electrical relays. With such systems, standard describing function techniques are typically too conservative, and time-domain simulation analysis is prohibitively extensive, This paper presents an open-loop analysis technique for this class of non-linear systems. The technique is centered around an innovative use of multivariable signal transmission theory to quantify the plant response to worst case control commands. The technique has been applied to assess stability of thruster controlled flexible space structures. Examples are provided for Space Shuttle attitude control with attached flexible payloads.
Non-linear canonical correlation for joint analysis of MEG signals from two subjects
Campi, Cristina; Parkkonen, Lauri; Hari, Riitta; Hyvärinen, Aapo
2013-01-01
Traditional stimulus-based analysis methods of magnetoencephalography (MEG) data are often dissatisfactory when applied to naturalistic experiments where two or more subjects are measured either simultaneously or sequentially. To uncover the commonalities in the brain activity of the two subjects, we propose a method that searches for linear transformations that output maximally correlated signals between the two brains. Our method is based on canonical correlation analysis (CCA), which provides linear transformations, one for each subject, such that the temporal correlation between the transformed MEG signals is maximized. Here, we present a non-linear version of CCA which measures the correlation of energies and allows for a variable delay between the time series to accommodate, e.g., leader–follower changes. We test the method with simulations and with MEG data from subjects who received the same naturalistic stimulus sequence. The method may help analyse future experiments where the two subjects are measured simultaneously while engaged in social interaction. PMID:23785311
[The linearity analysis of ultrahigh temperature FTIR spectral emissivity measurement system].
Wang, Zong-wei; Dai, Jing-min; He, Xiao-wa; Yang, Chun-ling
2012-02-01
To study thermal radiation properties of special materials at high temperature in aerospace fields, the ultrahigh temperature spectral emissivity measurement system with Fourier spectrometer has been established. The linearity of system is the guarantee of emissivity measurement precision. Through measuring spectral radiation signals of a blackbody source at different temperatures, the function relations between spectral signal values and blackbody spectral radiation luminance of every spectrum points were calculated with the method of multi-temperature and multi-spectrum linear fitting. The spectral radiation signals of blackbody were measured between 1 000 degrees C and 2 000 degrees C in the spectral region from 3 to 20 microm. The linear relations between spectral signal and theory line at wavelength of 4 microm were calculated and introduced. The spectral response is well good between 4 and 18 microm, the spectral linearity are less than 1% except CO2 strong absorption spectrum regions. The results show that when the errors of measured spectrum radiation and linear fitting theory lines are certain, the higher the temperature, the smaller the spectral errors on emissivity. The linearity analysis of spectrum response is good at eliminating errors caused by individual temperature' disturbance to the spectra. PMID:22512159
Analysis and improvement for a linearized seafloor elastic parameter inversion method
NASA Astrophysics Data System (ADS)
Liu, Yangting; Liu, Xuewei; Ning, Hongxiao
2016-05-01
AVO inversion is an effective seismic exploration method to predict elastic parameters. In this paper, we review and analyze the linearized AVO inversion method previously published for seafloor elastic parameters, and present a modification strategy. Before the linearized inversion is performed, a proper near-angle range in which the relationship between the reflection coefficient and sine-squared incidence angle is linear needs to be provided. However, the near-angle range is determined by the elastic parameters which are to be estimated by inversion. Therefore, only an approximated value of the near-angle range can be provided for the linearized inversion. Model tests show that a too large near-angle range may cause inversion fault, and a too small near-angle range may cause unreliable estimation. Further analysis shows that the estimation stability can be further improved even though the linearized inversion is performed under an exact near-angle range. To mitigate the strong dependence on the near-angle range, we use the seafloor elastic parameters estimated from the linearized method as the initial model for an unconstrained optimization method. Compared with the previously published method, the modified method is more robust to noisy data and shows less dependence on the near-angle range.
Barnaby, H.; Tausch, H.J.; Turfler, R.; Cole, P.; Baker, P.; Pease, R.L.
1996-12-01
A methodology is presented for the identification of circuit total dose response mechanisms in bipolar linear microcircuits irradiated at high and low dose rates. This methodology includes manual circuit analysis, circuit simulations with SPICE using extracted device parameters, and selective irradiations of portions of the circuit using a scanning electron microscope.
Micosoft Excel Sensitivity Analysis for Linear and Stochastic Program Feed Formulation
Technology Transfer Automated Retrieval System (TEKTRAN)
Sensitivity analysis is a part of mathematical programming solutions and is used in making nutritional and economic decisions for a given feed formulation problem. The terms, shadow price and reduced cost, are familiar linear program (LP) terms to feed formulators. Because of the nonlinear nature of...
Nonmodal linear stability analysis of miscible viscous fingering in porous media
NASA Astrophysics Data System (ADS)
Hota, Tapan Kumar; Pramanik, Satyajit; Mishra, Manoranjan
2015-11-01
The nonmodal linear stability of miscible viscous fingering in a two-dimensional homogeneous porous medium has been investigated. The linearized perturbed equations for Darcy's law coupled with a convection-diffusion equation is discretized using a finite difference method. The resultant initial value problem is solved by a fourth-order Runge-Kutta method, followed by a singular value decomposition of the propagator matrix. Particular attention is given to the transient behavior rather than the long-time behavior of eigenmodes predicted by the traditional modal analysis. The transient behaviors of the response to external excitations and the response to initial conditions are studied by examining the ɛ -pseudospectra structures and the largest energy growth function, respectively. With the help of nonmodal stability analysis we demonstrate that at early times the displacement flow is dominated by diffusion and the perturbations decay. At later times, when convection dominates diffusion, perturbations grow. Furthermore, we show that the dominant perturbation that experiences the maximum amplification within the linear regime lead to the transient growth. These two important features were previously unattainable in the existing linear stability methods for miscible viscous fingering. To explore the relevance of the optimal perturbation obtained from nonmodal analysis, we performed direct numerical simulations using a highly accurate pseudospectral method. Furthermore, a comparison of the present stability analysis with existing modal and initial value approach is also presented. It is shown that the nonmodal stability results are in better agreement than the other existing stability analyses, with those obtained from direct numerical simulations.
Principal Component Analysis: Resources for an Essential Application of Linear Algebra
ERIC Educational Resources Information Center
Pankavich, Stephen; Swanson, Rebecca
2015-01-01
Principal Component Analysis (PCA) is a highly useful topic within an introductory Linear Algebra course, especially since it can be used to incorporate a number of applied projects. This method represents an essential application and extension of the Spectral Theorem and is commonly used within a variety of fields, including statistics,…
A Hierarchical Linear Model with Factor Analysis Structure at Level 2
ERIC Educational Resources Information Center
Miyazaki, Yasuo; Frank, Kenneth A.
2006-01-01
In this article the authors develop a model that employs a factor analysis structure at Level 2 of a two-level hierarchical linear model (HLM). The model (HLM2F) imposes a structure on a deficient rank Level 2 covariance matrix [tau], and facilitates estimation of a relatively large [tau] matrix. Maximum likelihood estimators are derived via the…
Li, M.; Sun, Z.; Zeng, D.
1996-12-31
Linear stability theory was applied in the present paper to analyze the stability of the basic state solution of thermocapillary convection in a liquid bridge with liquid encapsulation. Discretizing the linearized disturbance equations by using finite-difference approximation, stability analysis is evolved to a complex generalized eigenvalue problem with a complicated band structure of matrix. The influence of dimensionless parameters on the stability of the system is revealed through solving the complex generalized eigenvalue problem by inverse iteration. The results provide a theoretical and numerical foundation for crystal growth by float-zone method and for other engineering applications.
NASA Technical Reports Server (NTRS)
Oden, J. Tinsley; Fly, Gerald W.; Mahadevan, L.
1987-01-01
A hybrid stress finite element method is developed for accurate stress and vibration analysis of problems in linear anisotropic elasticity. A modified form of the Hellinger-Reissner principle is formulated for dynamic analysis and an algorithm for the determination of the anisotropic elastic and compliance constants from experimental data is developed. These schemes were implemented in a finite element program for static and dynamic analysis of linear anisotropic two dimensional elasticity problems. Specific numerical examples are considered to verify the accuracy of the hybrid stress approach and compare it with that of the standard displacement method, especially for highly anisotropic materials. It is that the hybrid stress approach gives much better results than the displacement method. Preliminary work on extensions of this method to three dimensional elasticity is discussed, and the stress shape functions necessary for this extension are included.
NASA Astrophysics Data System (ADS)
Qiao, Lingbo; Wang, Yingxin; Zhao, Ziran; Chen, Zhiqiang
2014-07-01
Terahertz (THz) time-domain spectroscopy is considered as an attractive tool for the analysis of chemical composition. The traditional methods for identification and quantitative analysis of chemical compounds by THz spectroscopy are all based on full-spectrum data. However, intrinsic features of the THz spectrum only lie in absorption peaks due to existence of disturbances, such as unexpected components, scattering effects, and barrier materials. We propose a strategy that utilizes Lorentzian parameters of THz absorption peaks, extracted by a multiscale linear fitting method, for both identification of pure chemicals and quantitative analysis of mixtures. The multiscale linear fitting method can automatically remove background content and accurately determine Lorentzian parameters of the absorption peaks. The high recognition rate for 16 pure chemical compounds and the accurate predicted concentrations for theophylline-lactose mixtures demonstrate the practicability of our approach.
NASA Astrophysics Data System (ADS)
Peeling, S. M.; Ponting, K. M.
1991-12-01
Linear discriminant analysis is used to generate speech data transformations. This transformed data is then used within the Airborne Reconnaissance Mission (ARM) continuous speech recognition system. The aim of the ARM project is accurate recognition of continuously spoken airborne reconnaissance reports using a speech recognition system based on phoneme level hidden Markov models. A fuller description of a linear discriminant analysis, which is applied to speaker dependent data in the ARM system, is given. Preliminary results are presented from experiments using transformed data alone and also in conjunction with one, or both, of the word transition penalties and variable frame rate analysis. Speaker dependent results are reported which are significantly better then the best obtained previously.
A Three-Dimensional Linearized Unsteady Euler Analysis for Turbomachinery Blade Rows
NASA Technical Reports Server (NTRS)
Montgomery, Matthew D.; Verdon, Joseph M.
1996-01-01
A three-dimensional, linearized, Euler analysis is being developed to provide an efficient unsteady aerodynamic analysis that can be used to predict the aeroelastic and aeroacoustic response characteristics of axial-flow turbomachinery blading. The field equations and boundary conditions needed to describe nonlinear and linearized inviscid unsteady flows through a blade row operating within a cylindrical annular duct are presented. In addition, a numerical model for linearized inviscid unsteady flow, which is based upon an existing nonlinear, implicit, wave-split, finite volume analysis, is described. These aerodynamic and numerical models have been implemented into an unsteady flow code, called LINFLUX. A preliminary version of the LINFLUX code is applied herein to selected, benchmark three-dimensional, subsonic, unsteady flows, to illustrate its current capabilities and to uncover existing problems and deficiencies. The numerical results indicate that good progress has been made toward developing a reliable and useful three-dimensional prediction capability. However, some problems, associated with the implementation of an unsteady displacement field and numerical errors near solid boundaries, still exist. Also, accurate far-field conditions must be incorporated into the FINFLUX analysis, so that this analysis can be applied to unsteady flows driven be external aerodynamic excitations.
A spectral analysis of the domain decomposed Monte Carlo method for linear systems
Slattery, S. R.; Wilson, P. P. H.; Evans, T. M.
2013-07-01
The domain decomposed behavior of the adjoint Neumann-Ulam Monte Carlo method for solving linear systems is analyzed using the spectral properties of the linear operator. Relationships for the average length of the adjoint random walks, a measure of convergence speed and serial performance, are made with respect to the eigenvalues of the linear operator. In addition, relationships for the effective optical thickness of a domain in the decomposition are presented based on the spectral analysis and diffusion theory. Using the effective optical thickness, the Wigner rational approximation and the mean chord approximation are applied to estimate the leakage fraction of stochastic histories from a domain in the decomposition as a measure of parallel performance and potential communication costs. The one-speed, two-dimensional neutron diffusion equation is used as a model problem to test the models for symmetric operators. In general, the derived approximations show good agreement with measured computational results. (authors)
Radiative-neutron-capture gamma-ray analysis by a linear combination technique
Tanner, A.B.; Bhargava, R.C.; Senftle, F.E.; Brinkerhoff, J.M.
1972-01-01
The linear combination technique, when applied to a gamma-ray spectrum, gives a single number indicative of the extent to which the spectral lines of a sought element are present in a complex spectrum. Spectra are taken of the sought element and of various other substances whose spectra interfere with that of the sought element. A weighting function is then computed for application to spectra of unknown materials. The technique was used to determine calcium by radiative-neutron-capture gamma-ray analysis in the presence of interfering elements, notably titanium, and the results were compared with those for two popular methods of peak area integration. Although linearity of response was similar for the methods, the linear combination technique was much better at rejecting interferences. For analyses involving mixtures of unknown composition the technique consequently offers improved sensitivity. ?? 1972.
Analysis of linear and mildly nonlinear relationships using pooled subject data.
Poon, C S
1988-02-01
The detection of mild nonlinearities and/or state-dependent variability in otherwise linear physiological relationships is generally difficult in the presence of significant measurement errors. Conventional approaches using pooled subject data to increase the degree of freedom for statistical inference are enervated by the resultant introduction of intersubject variability. This paper proposes a new, simple method of pooling multiple subject data for linearity analysis. With the use of a special standardization procedure for the individual response curves, this method allows sensitive detection of occult nonlinearities as well as any state-dependent variability in the underlying relationship. Application of this analytic approach to reported hypercapnic exercise-response data in eight healthy subjects showed that 1) the hypercapnic ventilation-CO2 output relationship is nonlinear with a downward concavity; and 2) the ventilation-tidal volume relationship, which is linear at low tidal volume values, is similar in hypercapnic exercise as in resting hypercapnia or eucapnic exercise. PMID:3372442
Sammer, G
1998-05-01
In the investigation of heart rate and heart rate variability, the discrimination between mental workload, physical activity and respiration is known to be methodologically difficult. At most, heart rate variability measures are more likely to be coarse-grained measures with variability confounded by heart rate. Moreover, the spectral analysis of heart rate variability shows broad-band frequency characteristics, pointing towards non-stationarity or non-linearity. From this it is suggested to focus on non-linear dynamic analyses that are variance-insensitive. The experimental section of the paper focuses on the estimation of two non-linear measures for both heartbeat dynamics and respiration, the correlation dimension indicating complexity and the Lyapunov exponents indicating predictability. The results indicate that the complexity of heart dynamics is related to the type of task and that the predictability of heart dynamics is related to the amount of load. PMID:9613233
Analysis and synthesis of phase shifting algorithms based on linear systems theory
NASA Astrophysics Data System (ADS)
Servin, M.; Estrada, J. C.
2012-08-01
We review and update a recently published formalism for the theory of linear Phase Shifting Algorithms (PSAs) based on linear filtering (systems) theory, mainly using the Frequency Transfer Function (FTF). The FTF has been for decades the standard tool in Electrical Engineering to analyze and synthesize their linear systems. Given the well defined FTF approach (matured over the last century), it clarifies, in our view, many not fully understood properties of PSAs. We present easy formulae for the spectra of the PSAs (the FTF magnitude), their Signal to Noise (S/N) power-ratio gain, their detuning robustness, and their harmonic rejection in terms of the FTF. This paper has more practical appeal than previous publications by the same authors, hoping to enrich the understanding of this PSA's theory as applied to the analysis and synthesis of temporal interferometry algorithms in Optical Metrology.
On the application of a linearized unsteady potential-flow analysis to fan-tip cascades
NASA Technical Reports Server (NTRS)
Usab, W. J., Jr.; Verdon, J. M.
1986-01-01
A linearized potential flow analysis, which accounts for the effects of nonuniform steady flow phenomena on the unsteady response to prescribed blade motions, has been applied to five two-dimensional cascade configurations. These include a flat-plate cascade and three cascades which are representative of the tip sections of current fan designs. Here the blades are closely spaced, highly staggered, and operate at low mean incidence. The fifth configuration is a NASA Lewis cascade of symmetric biconvex airfoils for which experimental measurements are available. Numerical solutions are presented that clearly illustrate the effects and importance of blade geometry and mean blade loading on the linearized unsteady response at high subsonic inlet Mach number and high blade-vibrational frequency. In addition, a good qualitative agreement is shown between the analytical predictions and experimental measurements for the cascade of symmetric biconvex airfoils. Finally, recommendations on the research needed to extend the range of application of linearized unsteady aerodynamic analyses are provided.
Spherically symmetric analysis on open FLRW solution in non-linear massive gravity
Chiang, Chien-I; Izumi, Keisuke; Chen, Pisin E-mail: izumi@phys.ntu.edu.tw
2012-12-01
We study non-linear massive gravity in the spherically symmetric context. Our main motivation is to investigate the effect of helicity-0 mode which remains elusive after analysis of cosmological perturbation around an open Friedmann-Lemaitre-Robertson-Walker (FLRW) universe. The non-linear form of the effective energy-momentum tensor stemming from the mass term is derived for the spherically symmetric case. Only in the special case where the area of the two sphere is not deviated away from the FLRW universe, the effective energy momentum tensor becomes completely the same as that of cosmological constant. This opens a window for discriminating the non-linear massive gravity from general relativity (GR). Indeed, by further solving these spherically symmetric gravitational equations of motion in vacuum to the linear order, we obtain a solution which has an arbitrary time-dependent parameter. In GR, this parameter is a constant and corresponds to the mass of a star. Our result means that Birkhoff's theorem no longer holds in the non-linear massive gravity and suggests that energy can probably be emitted superluminously (with infinite speed) on the self-accelerating background by the helicity-0 mode, which could be a potential plague of this theory.
Dynamic Analysis of Flexible Slider-Crank Mechanisms with Non-Linear Finite Element Method
NASA Astrophysics Data System (ADS)
CHEN, J.-S.; HUANG, C.-L.
2001-09-01
Previous research in finite element formulation of flexible mechanisms usually neglected high order terms in the strain-energy function. In particular, the quartic term of the displacement gradient is always neglected due to the common belief that it is not important in the dynamic analysis. In this paper, we show that this physical intuition is not always valid. By retaining all the high order terms in the strain-energy function the equations of motion naturally become non-linear, which can then be solved by the Newmark method. In the low-speed range it is found that the dynamic responses predicted by non-linear and linear approaches indeed make no significant difference. However, when the rotation speed increases up to about one-fifth of the fundamental bending natural frequency of the connecting rod, simplified approaches begin to incur noticeable error. Specifically, for a connecting rod with a slenderness ratio of 0·01 the conventional simplified approaches overestimate the vibration amplitude almost 10-fold when the rotation speed is comparable to the fundamental natural frequency of the connecting rod. Therefore, non-linear finite element formulation taking into account the complete non-linear strain is needed in analyzing high-speed flexible mechnisms with slender links.
NASA Astrophysics Data System (ADS)
Preuss, R. D.
1985-12-01
The theory of multipulse linear predictive coding (LPC) analysis is extended to include the possible presence of acoustic noise, as for a telephone near a busy road. Models are developed assuming two signals are provided: the primary signal is the output of a microphone which samples the combined acoustic fields of the noise and the speech, while the secondary signal is the output of a microphone which samples the acoustic field of the noise alone. Analysis techniques to extract the multipulse LPC parameters from these two signals are developed; these techniques are developed as approximations to maximum likelihood analysis for the given model.
Improvement of the narrowband linear predictive coder. Part 1: Analysis improvements
NASA Astrophysics Data System (ADS)
Kang, G. S.; Everett, S. S.
1982-12-01
The narrowband linear predictive coder (LPC) is widely used in both civilian and military applications. Yet, in spite of the many improvements over the years, it is still not universally acceptable to general users. This report presents improvements to five aspects of the LPC analysis which have drawn little attention in the past. These improvements are frequency spreading of sibilant sound spectra, adaptive placement of the analysis window at onsets, modified LPC analysis for sustained vowels, enhancement of LPC spectra in noise, and an automatic gain control. These improvements may be incorporated in the existing narrowband LPC without interfering with the speech sampling rate, the frame rate, or the parameter coding.
Non-linear analysis of PESA-Lo electrostatic analyzer data and solar wind temperature anisotropies
NASA Astrophysics Data System (ADS)
Djordjevic, B. Z.; Maruca, B.; Bale, S. D.; Wilson, L. B., III; Larson, D. E.
2015-12-01
In this study, non-linear fitting techniques are applied to ion measurements from the Wind spacecraft's PESA-Lo electrostatic analyzer. Previous studies have relied primarily on moments-analyses, which, although satisfactory for simple distributions and density calculations, often return unreasonable values for higher order moments (e.g., temperature) and fail to account for non-thermal effects (e.g., temperature anisotropy and beams) and multiple ion-species. A Levenberg-Marquadt non-linear algorithm is applied to the PESA-Lo data in order to calculate the characteristic parameters of the proton, alpha-particle, and beam distributions. This analysis is augmented with calibration data from the WIND Faraday cups and magnetic-field data from WIND/MFI. Preliminary results from this non-linear analysis indicate that it indeed provides higher-quality ion parameters than the existing moments-analysis. When this analysis is complete, the set of bulk-parameter values will be suitable for studies of microinstabilities in the solar wind and of possible correlation between magnetic field fluctuations and non-thermal properties of the ion distributions. Applications of thermodynamic principles, such as the Boltzmann H-theorem, will allow for further characterization of the non-thermal properties of the solar wind.
ERIC Educational Resources Information Center
Phillips, Gary W.
The usefulness of path analysis as a means of better understanding various linear models is demonstrated. First, two linear models are presented in matrix form using linear structural relations (LISREL) notation. The two models, regression and factor analysis, are shown to be identical although the research question and data matrix to which these…
Failure mode analysis of a post-tension anchored dam using linear finite element analysis
NASA Astrophysics Data System (ADS)
Corn, Aimee
There are currently over 84,000 dams in the United States, and the average age of those dams is 52 years. Concrete gravity dams are the second most common dam type, with more than 3,000 in the United States. Current engineering technology and technical understanding of hydrologic and seismic events has resulted in significant increases to the required design loads for most dams; therefore, many older dams do not have adequate safety for extreme loading events. Concrete gravity dams designed and constructed in the early 20th century did not consider uplift pressures beneath the dam, which reduces the effective weight of the structure. One method that has been used to enhance the stability of older concrete gravity dams includes the post-tension anchor (PTA) system. Post-tensioning infers modifying cured concrete and using self-equilibrating elements to increase the weight of the section, which provides added stability. There is a lack of historical evidence regarding the potential failure mechanisms for PTA concrete gravity dams. Of particular interest, is how these systems behave during large seismic events. The objective of this thesis is to develop a method by which the potential failure modes during a seismic event for a PTA dam can be evaluated using the linear elastic finite element method of analysis. The most likely potential failure modes (PFM) for PTA designs are due to tensile failure and shear failure. A numerical model of a hypothetical project was developed to simulate PTAs in the dam. The model was subjected to acceleration time-history motions that simulated the seismic loads. The results were used to evaluate the likelihood of tendon failure due to both tension and shear. The results from the analysis indicated that the PTA load increased during the seismic event; however, the peak load in the tendons was less than the gross ultimate tensile strength (GUTS) and would not be expected to result in tensile failure at the assumed project. The analysis
Linear analysis of the heave dynamics of a bag and finger air cushion vehicle skirt
NASA Astrophysics Data System (ADS)
Ma, T.; Sullivan, P. A.
1986-09-01
A linear analysis of the heave dynamics of an air-cushion vehicle bag and finger skirt is presented. A simplified geometry is considered; this is a two-dimensional section of the skirt without interior compartmentation. The bag is modeled as a membrane having distributed mass and viscoelasticity, and the fingers are modeled as rigid bodies having both mass and moment of inertia. A finite-element technique is used to discretize the equations of motion of the bag, but otherwise standard linear analysis techniques are used to obtain predictions of frequency response and stability characteristics. The stability results confirm the experimental observation that the dominant factor controlling the onset of skirt bounce is the bag-to-cushion pressure ratio.
Chi, Xianglin; Zhou, Jianhua; Shi, Ping; Liu, Chengyu
2016-02-01
The linear analysis for heart rate variability (HRV), including time domain method, frequency domain method and time-frequency analysis, has reached a lot of consensus. The non-linear analysis has also been widely applied in biomedical and clinical researches. However, for non-linear HRV analysis, especially for short-term non-linear HRV analysis, controversy still exists, and a unified standard and conclusion has not been formed. This paper reviews and discusses three short-term non-linear HRV analysis methods (fractal dimension, entropy and complexity) and their principles, progresses and problems in clinical application in detail, in order to provide a reference for accurate application in clinical medicine. PMID:27382764
Linear and nonlinear dynamic analysis by boundary element method. Ph.D. Thesis, 1986 Final Report
NASA Technical Reports Server (NTRS)
Ahmad, Shahid
1991-01-01
An advanced implementation of the direct boundary element method (BEM) applicable to free-vibration, periodic (steady-state) vibration and linear and nonlinear transient dynamic problems involving two and three-dimensional isotropic solids of arbitrary shape is presented. Interior, exterior, and half-space problems can all be solved by the present formulation. For the free-vibration analysis, a new real variable BEM formulation is presented which solves the free-vibration problem in the form of algebraic equations (formed from the static kernels) and needs only surface discretization. In the area of time-domain transient analysis, the BEM is well suited because it gives an implicit formulation. Although the integral formulations are elegant, because of the complexity of the formulation it has never been implemented in exact form. In the present work, linear and nonlinear time domain transient analysis for three-dimensional solids has been implemented in a general and complete manner. The formulation and implementation of the nonlinear, transient, dynamic analysis presented here is the first ever in the field of boundary element analysis. Almost all the existing formulation of BEM in dynamics use the constant variation of the variables in space and time which is very unrealistic for engineering problems and, in some cases, it leads to unacceptably inaccurate results. In the present work, linear and quadratic isoparametric boundary elements are used for discretization of geometry and functional variations in space. In addition, higher order variations in time are used. These methods of analysis are applicable to piecewise-homogeneous materials, such that not only problems of the layered media and the soil-structure interaction can be analyzed but also a large problem can be solved by the usual sub-structuring technique. The analyses have been incorporated in a versatile, general-purpose computer program. Some numerical problems are solved and, through comparisons
Non-linear dynamic analysis of ancient masonry structures by 3D rigid block models
NASA Astrophysics Data System (ADS)
Orduña, Agustin; Ayala, A. Gustavo
2015-12-01
This work presents a formulation for non-linear dynamic analysis of unreinforced masonry structures using rigid block models. This procedure is akin to the distinct element family of methods, nevertheless, we assume that small displacements occur and, therefore, the formulation does not involve the search for new contacts between blocks. This proposal is also related to the rigid element method, although, in this case we use full three-dimensional models and a more robust interface formulation.
Linearized stability analysis and design of a flyback dc-dc boost regulator.
NASA Technical Reports Server (NTRS)
Wester, G. W.
1973-01-01
Analytic expressions for the small-signal power-stage describing functions of a switched dc-dc boost regulator are derived from an approximate continuous circuit model which is developed by a time-averaging technique. Closed-loop stability is attained through the design of frequency compensation of the loop gain. Open- and closed-loop regulator output impedances are derived from the linearized models for the given configuration. The analysis and design are compared with and confirmed by breadboard measurements.
Nonmodal linear stability analysis of miscible viscous fingering in porous media.
Hota, Tapan Kumar; Pramanik, Satyajit; Mishra, Manoranjan
2015-11-01
The nonmodal linear stability of miscible viscous fingering in a two-dimensional homogeneous porous medium has been investigated. The linearized perturbed equations for Darcy's law coupled with a convection-diffusion equation is discretized using a finite difference method. The resultant initial value problem is solved by a fourth-order Runge-Kutta method, followed by a singular value decomposition of the propagator matrix. Particular attention is given to the transient behavior rather than the long-time behavior of eigenmodes predicted by the traditional modal analysis. The transient behaviors of the response to external excitations and the response to initial conditions are studied by examining the ε-pseudospectra structures and the largest energy growth function, respectively. With the help of nonmodal stability analysis we demonstrate that at early times the displacement flow is dominated by diffusion and the perturbations decay. At later times, when convection dominates diffusion, perturbations grow. Furthermore, we show that the dominant perturbation that experiences the maximum amplification within the linear regime lead to the transient growth. These two important features were previously unattainable in the existing linear stability methods for miscible viscous fingering. To explore the relevance of the optimal perturbation obtained from nonmodal analysis, we performed direct numerical simulations using a highly accurate pseudospectral method. Furthermore, a comparison of the present stability analysis with existing modal and initial value approach is also presented. It is shown that the nonmodal stability results are in better agreement than the other existing stability analyses, with those obtained from direct numerical simulations. PMID:26651779
LETTER TO THE EDITOR: Analysis of `gauge modes' in linearized relativity
NASA Astrophysics Data System (ADS)
Matzner, Richard A.; Huq, Mijan; Botero, Alonso; Choi, Dae Il; Kask, Ullar; Lara, Juan; Liebling, Steven; Neilsen, David; Premadi, Premana; Shoemaker, Deirdre
1997-01-01
By writing the complete set of 3 + 1 (ADM) equations for linearized waves, we are able to demonstrate the properties of the initial data and of the evolution of a wave problem set by Alcubierre and Schutz. We show that the gauge modes and constraint error modes arise in a straightforward way in the analysis, and are of a form which will be controlled in any well specified convergent computational discretization of the differential equations.
NASA Technical Reports Server (NTRS)
Jobson, Daniel J.; Rahman, Zia-Ur; Woodell, Glenn A.; Hines, Glenn D.
2004-01-01
Noise is the primary visibility limit in the process of non-linear image enhancement, and is no longer a statistically stable additive noise in the post-enhancement image. Therefore novel approaches are needed to both assess and reduce spatially variable noise at this stage in overall image processing. Here we will examine the use of edge pattern analysis both for automatic assessment of spatially variable noise and as a foundation for new noise reduction methods.
NASA Astrophysics Data System (ADS)
Polster, J.
2001-01-01
Linear reaction systems consist by definition of first-order reaction steps. Linearly independent reactions are independent of reaction order. Each reaction mechanism consists of a distinct number ( s) of linearly independent reaction steps. Thus, the mechanism A→B→C→D can be described by three linearly independent reactions as it is also true for the multiple equilibria A ⇋ B, C ⇋ D, E ⇋ F. A general method is developed for the spectroscopic-kinetic analysis of linear reactions ( s=3) on the basis of three-dimensional absorbance ( A) diagrams (A λ 1vs. A λ 2vs. A λ 3) . A distorted 'absorbance tetrahedron' can be constructed from the curve running in the absorbance space (called Mauser space). The tetrahedron is generated by tangents and osculating planes belonging to the initial point and endpoint of curve (measured). Planes being parallel to the tetrahedral surfaces and running through the points of curve, can be constructed and brought to intersection with the corresponding sides of tetrahedron. The quantities z i are introduced with the help of distance relationships on the sides of tetrahedron. The differentiation of z i with respect to time ( z i˙) leads to equations which are linearly dependent on z i. The solution of these differential equations provides the eigenvalues ( r i) searched. The results obtained are in accordance with Theorem 2 of kinetics (two strictly linear reaction systems having the same number of linearly independent reaction steps cannot be distinguished from each other by purely spectroscopic means). The procedure of evaluation is demonstrated by the investigation of the spontaneous hydrolyses of 4-methylumbelliferyl p-trimethylammonium cinnamate chloride, o-nitrophenylacetate and cinnamoylimidazole in borax buffer (0.1 M; pH=8.7; 10% acetonitrile; temperature 25.0°C).
Linear analysis of transient growth in stably-stratified, turbulent channel flow
NASA Astrophysics Data System (ADS)
Del Alamo, Juan Carlos; Yanez, Carlos; Garcia-Villalba, Manuel
2011-11-01
We studied stably-stratified, fully-developed, turbulent channel flow using linear stability analysis. The analysis considered the mean velocity and density profiles extracted from DNS calculations as base flow, and included a linear model to represent the energy dissipation and scalar diffusion felt at the large scales as a consequence of the small scales. The flow was found to be asymptotically stable in all cases but transient growth of initial perturbations was observed. The perturbations showing maximal transient growth corresponded well with spanwise waves in the center of the channel and with streaks in the near wall region, which were both observed in the DNS. In particular, their sizes and convection velocities were reasonably well predicted by the linear model. Component-wise analysis revealed that, while the streaks were formed by the vertical stirring of mean shear, the transient amplification of the spanwise waves was reminiscent of the Orr mechanism. Supporte by the German Research Foundation, Project GA 1360/2-1.
The non-linear analysis of multi-support rotor-bearing systems
Kicinski, J.; Drozdowski, R.
1995-12-31
This paper contains selected parts of the simulation research of large rotor machines (200 MW power turbine-sets). These investigations were based on a non-linear theoretical model and the NLDW computer program, and were carried out in the Institute of Fluid-Flow Machinery of PAS. A trial has been performed of the optimization of system-dynamic properties, through the suitable selection of thermally deformed bearing-bush centers line -- the so called ``hot`` line -- (due to a rotor`s geodesic line), as well as the selection of the external fixing stiffness of bearing supports. Examples are also included of the orbits of selected system nodes for two differently powered turbine-sets. On this basis, an analysis of the stability of those turbines was achieved. A significant objective of this paper is also to point out some possibilities of applying the simulation research, based on a non-linear description of the system, to the diagnostics of rotor-machinery. Non-linear analysis facilitates the possibility of easily generating vibration spectra, as well as creating simulation waterfall graphs. These properties of nonlinear analysis create convenient conditions for gaining specific diagnostic information.
Finite Element Analysis and Optimization of Flexure Bearing for Linear Motor Compressor
NASA Astrophysics Data System (ADS)
Khot, Maruti; Gawali, Bajirao
Nowadays linear motor compressors are commonly used in miniature cryocoolers instead of rotary compressors because rotary compressors apply large radial forces to the piston, which provide no useful work, cause large amount of wear and usually require lubrication. Recent trends favour flexure supported configurations for long life. The present work aims at designing and geometrical optimization of flexure bearings using finite element analysis and the development of design charts for selection purposes. The work also covers the manufacturing of flexures using different materials and the validation of the experimental finite element analysis results.
Design and dynamic analysis of permanent magnet linear synchronous machine for servoapplication
NASA Astrophysics Data System (ADS)
Jang, Seok-Myeong; You, Dae-Joon; Jang, Won-Bum; Han, Sang-Chul
2006-04-01
In servosystem demanding precision dynamic characteristics, application of the permanent magnet linear synchronous machines (PMLSMs) has an advantage of analysis convenience by simple geometry and thrust ripple reduction from the sinusoidal back electromotive force and excited stator. Therefore, this paper presents the design of surface-mounted PMLSM with a slotless iron-cored stator according to coil turns to satisfy the rate thrust. Also, from the dynamic analysis for servoapplication of manufactured motor with heavy mass, we offer an accurate range of the dc link voltage and acceleration in rate speed. This is applied to speed reference profile considering the system characteristics in total length of moving position.
An improved multiple linear regression and data analysis computer program package
NASA Technical Reports Server (NTRS)
Sidik, S. M.
1972-01-01
NEWRAP, an improved version of a previous multiple linear regression program called RAPIER, CREDUC, and CRSPLT, allows for a complete regression analysis including cross plots of the independent and dependent variables, correlation coefficients, regression coefficients, analysis of variance tables, t-statistics and their probability levels, rejection of independent variables, plots of residuals against the independent and dependent variables, and a canonical reduction of quadratic response functions useful in optimum seeking experimentation. A major improvement over RAPIER is that all regression calculations are done in double precision arithmetic.
Classical linear-control analysis applied to business-cycle dynamics and stability
NASA Technical Reports Server (NTRS)
Wingrove, R. C.
1983-01-01
Linear control analysis is applied as an aid in understanding the fluctuations of business cycles in the past, and to examine monetary policies that might improve stabilization. The analysis shows how different policies change the frequency and damping of the economic system dynamics, and how they modify the amplitude of the fluctuations that are caused by random disturbances. Examples are used to show how policy feedbacks and policy lags can be incorporated, and how different monetary strategies for stabilization can be analytically compared. Representative numerical results are used to illustrate the main points.
Non-linear Creep Analysis of Ceramic Specimen Using Finite Element Method
NASA Astrophysics Data System (ADS)
Saini, Jaswinder Singh; Khera, Saurabh
2016-07-01
In the present work the stress analysis of a ceramic tensile specimen is obtained. The effects of specimen geometry along with the pin loading are considered in the stress distribution calculations. Thereafter, the optimization based on a set of constraints is performed on the specimen with pinhole location, pinhole diameter, head width, neck radius and gauge length as its design variables. The work is then extended for the non-linear analysis for creep. A mathematical model is developed which is implemented using C++ code.
Non-linear Creep Analysis of Ceramic Specimen Using Finite Element Method
NASA Astrophysics Data System (ADS)
Saini, Jaswinder Singh; Khera, Saurabh
2016-03-01
In the present work the stress analysis of a ceramic tensile specimen is obtained. The effects of specimen geometry along with the pin loading are considered in the stress distribution calculations. Thereafter, the optimization based on a set of constraints is performed on the specimen with pinhole location, pinhole diameter, head width, neck radius and gauge length as its design variables. The work is then extended for the non-linear analysis for creep. A mathematical model is developed which is implemented using C++ code.
Flexible analysis of digital PCR experiments using generalized linear mixed models.
Vynck, Matthijs; Vandesompele, Jo; Nijs, Nele; Menten, Björn; De Ganck, Ariane; Thas, Olivier
2016-09-01
The use of digital PCR for quantification of nucleic acids is rapidly growing. A major drawback remains the lack of flexible data analysis tools. Published analysis approaches are either tailored to specific problem settings or fail to take into account sources of variability. We propose the generalized linear mixed models framework as a flexible tool for analyzing a wide range of experiments. We also introduce a method for estimating reference gene stability to improve accuracy and precision of copy number and relative expression estimates. We demonstrate the usefulness of the methodology on a complex experimental setup. PMID:27551671
A FORTRAN program for the analysis of linear continuous and sample-data systems
NASA Technical Reports Server (NTRS)
Edwards, J. W.
1976-01-01
A FORTRAN digital computer program which performs the general analysis of linearized control systems is described. State variable techniques are used to analyze continuous, discrete, and sampled data systems. Analysis options include the calculation of system eigenvalues, transfer functions, root loci, root contours, frequency responses, power spectra, and transient responses for open- and closed-loop systems. A flexible data input format allows the user to define systems in a variety of representations. Data may be entered by inputing explicit data matrices or matrices constructed in user written subroutines, by specifying transfer function block diagrams, or by using a combination of these methods.
Pleiotropy analysis of quantitative traits at gene level by multivariate functional linear models.
Wang, Yifan; Liu, Aiyi; Mills, James L; Boehnke, Michael; Wilson, Alexander F; Bailey-Wilson, Joan E; Xiong, Momiao; Wu, Colin O; Fan, Ruzong
2015-05-01
In genetics, pleiotropy describes the genetic effect of a single gene on multiple phenotypic traits. A common approach is to analyze the phenotypic traits separately using univariate analyses and combine the test results through multiple comparisons. This approach may lead to low power. Multivariate functional linear models are developed to connect genetic variant data to multiple quantitative traits adjusting for covariates for a unified analysis. Three types of approximate F-distribution tests based on Pillai-Bartlett trace, Hotelling-Lawley trace, and Wilks's Lambda are introduced to test for association between multiple quantitative traits and multiple genetic variants in one genetic region. The approximate F-distribution tests provide much more significant results than those of F-tests of univariate analysis and optimal sequence kernel association test (SKAT-O). Extensive simulations were performed to evaluate the false positive rates and power performance of the proposed models and tests. We show that the approximate F-distribution tests control the type I error rates very well. Overall, simultaneous analysis of multiple traits can increase power performance compared to an individual test of each trait. The proposed methods were applied to analyze (1) four lipid traits in eight European cohorts, and (2) three biochemical traits in the Trinity Students Study. The approximate F-distribution tests provide much more significant results than those of F-tests of univariate analysis and SKAT-O for the three biochemical traits. The approximate F-distribution tests of the proposed functional linear models are more sensitive than those of the traditional multivariate linear models that in turn are more sensitive than SKAT-O in the univariate case. The analysis of the four lipid traits and the three biochemical traits detects more association than SKAT-O in the univariate case. PMID:25809955
Dynamic Linear Model Analysis of Optical Imaging Data Acquired from the Human Neocortex
Lavine, Michael; Haglund, Michael M.; Hochman, Daryl W.
2011-01-01
The amount of light absorbed and scattered by neocortical tissue is altered by neuronal activity. Imaging of intrinsic optical signals (ImIOS), a technique for mapping these activity-evoked optical changes with an imaging detector, has the potential to be useful for both clinical and experimental investigations of the human neocortex. However, its usefulness for human studies is currently limited because intraoperatively acquired ImIOS data is noisy. To improve the reliability and usefulness of ImIOS for human studies, it is desirable to find appropriate methods for the removal of noise artifacts and its statistical analysis. Here we develop a Bayesian, dynamic linear modeling approach that appears to address these problems. A dynamic linear model (DLM) was constructed that included cyclic components to model the heartbeat and respiration artifacts, and a local linear component to model the activity-evoked response. The robustness of the model was tested on a set of ImIOS data acquired from the exposed cortices of six human subjects illuminated with either 535 nm or 660 nm light. The DLM adequately reduced noise artifacts in these data while reliably preserving their activity-evoked optical responses. To demonstrate how these methods might be used for intraoperative neurosurgical mapping, optical data acquired from a single human subject during direct electrical stimulation of the cortex were quantitatively analyzed. This example showed that the DLM can be used to provide quantitative information about human ImIOS data that is not available through qualitative analysis alone. PMID:21640137
Schultze, J.; Hemez, F.
2000-11-01
This research presents a new method to improve analytical model fidelity for non-linear systems. The approach investigates several mechanisms to assist the analyst in updating an analytical model based on experimental data and statistical analysis of parameter effects. The first is a new approach at data reduction called feature extraction. This is an expansion of the classic update metrics to include specific phenomena or character of the response that is critical to model application. This is an extension of the familiar linear updating paradigm of utilizing the eigen-parameters or FRF's to include such devices as peak acceleration, time of arrival or standard deviation of model error. The next expansion of the updating process is the inclusion of statistical based parameter analysis to quantify the effects of uncertain or significant effect parameters in the construction of a meta-model. This provides indicators of the statistical variation associated with parameters as well as confidence intervals on the coefficients of the resulting meta-model. Also included in this method is the investigation of linear parameter effect screening using a partial factorial variable array for simulation. This is intended to aid the analyst in eliminating from the investigation the parameters that do not have a significant variation effect on the feature metric. Finally an investigation of the model to replicate the measured response variation is examined.
Worst-case analysis and linear parameter-varying gain-scheduled control of aerospace systems
NASA Astrophysics Data System (ADS)
Shin, Jong-Yeob
In this thesis, two main subjects are discussed. The first is a worst-case performance analysis, the second is a linear parameter varying (LPV) synthesis using a blending approach. On the first subject, a linear fractional transformation (LFT) model of the linearized X-38 Crew Return Vehicle (CRV) has been developed to facilitate the analysis of its flight control system. The LFT model represents uncertainty in nine aerodynamic stability derivatives at a given flight condition. The X-38 LFT model, combined with a controller at specific flight conditions, is used to determine the aerodynamic coefficients within a predefined set that result in the worst-case performance and worst-case gain/phase margins of the closed-loop system. LPV and mu controllers are synthesized for the X-38 CRV lateral-directional axes over the candidate flight envelope and compared with the baseline gain-scheduled classical control design. Worst-case analysis of the LPV and mu controllers are compared with the baseline gain-scheduled classical control design. Analysis and time simulations show that the LPV controller achieves significant performance and robustness improvements when compared to a linear mu controller and the baseline gain-scheduled controller. On the second subject, a quasi-LPV model of the F-16 longitudinal axes was developed using three methods: Jacobian linearization, state transformation and function substitution. Time simulations of quasi-LPV models show that the quasi-LPV models developed using state transformation and function substitution accurately represent the nonlinear dynamics of the F-16 longitudinal axes. In designing an LPV controller for the F-16 longitudinal axes, the function substitution quasi-LPV models are used since these quasi-LPV models can represent the nonlinear dynamics at non-trim points. Two LPV controllers are synthesized for the F-16 longitudinal axes for two separated flight envelopes: low and high altitude regions. Blending these controllers
Refining and end use study of coal liquids II - linear programming analysis
Lowe, C.; Tam, S.
1995-12-31
A DOE-funded study is underway to determine the optimum refinery processing schemes for producing transportation fuels that will meet CAAA regulations from direct and indirect coal liquids. The study consists of three major parts: pilot plant testing of critical upgrading processes, linear programming analysis of different processing schemes, and engine emission testing of final products. Currently, fractions of a direct coal liquid produced form bituminous coal are being tested in sequence of pilot plant upgrading processes. This work is discussed in a separate paper. The linear programming model, which is the subject of this paper, has been completed for the petroleum refinery and is being modified to handle coal liquids based on the pilot plant test results. Preliminary coal liquid evaluation studies indicate that, if a refinery expansion scenario is adopted, then the marginal value of the coal liquid (over the base petroleum crude) is $3-4/bbl.
A 3-D Magnetic Analysis of a Linear Alternator For a Stirling Power System
NASA Technical Reports Server (NTRS)
Geng, Steven M.; Schwarze, Gene E.; Niedra, Janis M.
2000-01-01
The NASA Glenn Research Center and the Department of Energy (DOE) are developing advanced radioisotope Stirling convertors, under contract with Stirling Technology Company (STC), for space applications. Of critical importance to the successful development of the Stirling convertor for space power applications is the development of a lightweight and highly efficient linear alternator. This paper presents a 3-D finite element method (FEM) approach for evaluating Stirling convertor linear alternators. Preliminary correlations with open-circuit voltage measurements provide an encouraging level of confidence in the model. Spatial plots of magnetic field strength (H) are presented in the region of the exciting permanent magnets. These plots identify regions of high H, where at elevated temperature and under electrical load, the potential to alter the magnetic moment of the magnets exists. This implies the need for further testing and analysis.
NASA Astrophysics Data System (ADS)
Singleton, V. L.; Gantzer, P.; Little, J. C.
2007-02-01
An existing linear bubble plume model was improved, and data collected from a full-scale diffuser installed in Spring Hollow Reservoir, Virginia, were used to validate the model. The depth of maximum plume rise was simulated well for two of the three diffuser tests. Temperature predictions deviated from measured profiles near the maximum plume rise height, but predicted dissolved oxygen profiles compared very well with observations. A sensitivity analysis was performed. The gas flow rate had the greatest effect on predicted plume rise height and induced water flow rate, both of which were directly proportional to gas flow rate. Oxygen transfer within the hypolimnion was independent of all parameters except initial bubble radius and was inversely proportional for radii greater than approximately 1 mm. The results of this work suggest that plume dynamics and oxygen transfer can successfully be predicted for linear bubble plumes using the discrete-bubble approach.
Non-Linear Analysis of Mode II Fracture in the end Notched Flexure Beam
NASA Astrophysics Data System (ADS)
Rizov, V.
2016-03-01
Analysis is carried-out of fracture in the End Notched Flex- ure (ENF) beam configuration, taking into account the material nonlin- earity. For this purpose, the J-integral approach is applied. A non-linear model, based on the Classical beam theory is used. The mechanical be- haviour of the ENF configuration is described by the Ramberg-Osgood stress-strain curve. It is assumed that the material possesses the same properties in tension and compression. The influence is evaluated of the material constants in the Ramberg-Osgood stress-strain equation on the fracture behaviour. The effect of the crack length on the J-integral value is investigated, too. The analytical approach, developed in the present paper, is very useful for parametric analyses, since the simple formulae obtained capture the essentials of the non-linear fracture in the ENF con- figuration.
Analysis of a Linear System for Variable-Thrust Control in the Terminal Phase of Rendezvous
NASA Technical Reports Server (NTRS)
Hord, Richard A.; Durling, Barbara J.
1961-01-01
A linear system for applying thrust to a ferry vehicle in the 3 terminal phase of rendezvous with a satellite is analyzed. This system requires that the ferry thrust vector per unit mass be variable and equal to a suitable linear combination of the measured position and velocity vectors of the ferry relative to the satellite. The variations of the ferry position, speed, acceleration, and mass ratio are examined for several combinations of the initial conditions and two basic control parameters analogous to the undamped natural frequency and the fraction of critical damping. Upon making a desirable selection of one control parameter and requiring minimum fuel expenditure for given terminal-phase initial conditions, a simplified analysis in one dimension practically fixes the choice of the remaining control parameter. The system can be implemented by an automatic controller or by a pilot.
Optical analysis and thermal management of 2-cell strings linear concentrating photovoltaic system
NASA Astrophysics Data System (ADS)
Reddy, K. S.; Kamnapure, Nikhilesh R.
2015-09-01
This paper presents the optical and thermal analyses for a linear concentrating photovoltaic/thermal collector under different operating conditions. Linear concentrating photovoltaic system (CPV) consists of a highly reflective mirror, a receiver and semi-dual axis tracking mechanism. The CPV receiver embodies two strings of triple-junction cells (100 cells in each string) adhered to a mild steel circular tube mounted at the focal length of trough. This system provides 560 W of electricity and 1580 W of heat which needs to be dissipated by active cooling. The Al2O3/Water nanofluid is used as heat transfer fluid (HTF) flowing through circular receiver for CPV cells cooling. Optical analysis of linear CPV system with 3.35 m2 aperture and geometric concentration ratio (CR) of 35 is carried out using Advanced System Analysis Program (ASAP) an optical simulation tool. Non-uniform intensity distribution model of solar disk is used to model the sun in ASAP. The impact of random errors including slope error (σslope), tracking error (σtrack) and apparent change in sun's width (σsun) on optical performance of collector is shown. The result from the optical simulations shows the optical efficiency (ηo) of 88.32% for 2-cell string CPV concentrator. Thermal analysis of CPV receiver is carried out with conjugate heat transfer modeling in ANSYS FLUENT-14. Numerical simulations of Al2O3/Water nanofluid turbulent forced convection are performed for various parameters such as nanoparticle volume fraction (φ), Reynolds number (Re). The addition of the nanoparticle in water enhances the heat transfer in the ranges of 3.28% - 35.6% for φ = 1% - 6%. Numerical results are compared with literature data which shows the reasonable agreement.
Kazmi, Syed Zaki Hassan; Zhang, Henggui; Aziz, Wajid; Monfredi, Oliver; Abbas, Syed Ali; Shah, Saeed Arif; Kazmi, Syeda Sobia Hassan; Butt, Wasi Haider
2016-01-01
The dynamical fluctuations in the rhythms of biological systems provide valuable information about the underlying functioning of these systems. During the past few decades analysis of cardiac function based on the heart rate variability (HRV; variation in R wave to R wave intervals) has attracted great attention, resulting in more than 17000-publications (PubMed list). However, it is still controversial about the underling mechanisms of HRV. In this study, we performed both linear (time domain and frequency domain) and nonlinear analysis of HRV data acquired from humans and animals to identify the relationship between HRV and heart rate (HR). The HRV data consists of the following groups: (a) human normal sinus rhythm (n = 72); (b) human congestive heart failure (n = 44); (c) rabbit sinoatrial node cells (SANC; n = 67); (d) conscious rat (n = 11). In both human and animal data at variant pathological conditions, both linear and nonlinear analysis techniques showed an inverse correlation between HRV and HR, supporting the concept that HRV is dependent on HR, and therefore, HRV cannot be used in an ordinary manner to analyse autonomic nerve activity of a heart. PMID:27336907
Linear Covariance Analysis For Proximity Operations Around Asteroid 2008 EV5
NASA Technical Reports Server (NTRS)
Wright, Cinnamon A.; Bhatt, Sagar; Woffinden, David; Strube, Matthew; D'Souza, Chris
2015-01-01
The NASA initiative to collect an asteroid, the Asteroid Robotic Redirect Mission (ARRM), is currently investigating the option of retrieving a boulder from an asteroid, demonstrating planetary defense with an enhanced gravity tractor technique, and returning it to a lunar orbit. Techniques for accomplishing this are being investigated by the Satellite Servicing Capabilities Office (SSCO) at NASA GSFC in collaboration with JPL, NASA JSC, LaRC, and Draper Laboratory, Inc. Two critical phases of the mission are the descent to the boulder and the Enhanced Gravity Tractor demonstration. A linear covariance analysis is done for these phases to assess the feasibility of these concepts with the proposed design of the sensor and actuator suite of the Asteroid Redirect Vehicle (ARV). The sensor suite for this analysis includes a wide field of view camera, LiDAR, and an IMU. The proposed asteroid of interest is currently the C-type asteroid 2008 EV5, a carbonaceous chondrite that is of high interest to the scientific community. This paper presents an overview of the linear covariance analysis techniques and simulation tool, provides sensor and actuator models, and addresses the feasibility of descending to the surface of the asteroid within allocated requirements as well as the possibility of maintaining a halo orbit to demonstrate the Enhanced Gravity Tractor technique.
Zhang, Henggui; Aziz, Wajid; Monfredi, Oliver; Abbas, Syed Ali; Shah, Saeed Arif; Kazmi, Syeda Sobia Hassan; Butt, Wasi Haider
2016-01-01
The dynamical fluctuations in the rhythms of biological systems provide valuable information about the underlying functioning of these systems. During the past few decades analysis of cardiac function based on the heart rate variability (HRV; variation in R wave to R wave intervals) has attracted great attention, resulting in more than 17000-publications (PubMed list). However, it is still controversial about the underling mechanisms of HRV. In this study, we performed both linear (time domain and frequency domain) and nonlinear analysis of HRV data acquired from humans and animals to identify the relationship between HRV and heart rate (HR). The HRV data consists of the following groups: (a) human normal sinus rhythm (n = 72); (b) human congestive heart failure (n = 44); (c) rabbit sinoatrial node cells (SANC; n = 67); (d) conscious rat (n = 11). In both human and animal data at variant pathological conditions, both linear and nonlinear analysis techniques showed an inverse correlation between HRV and HR, supporting the concept that HRV is dependent on HR, and therefore, HRV cannot be used in an ordinary manner to analyse autonomic nerve activity of a heart. PMID:27336907
NASA Astrophysics Data System (ADS)
El-Hanbaly, A. M.; El-Shewy, E. K.; Sallah, M.; Darweesh, H. F.
2015-05-01
The propagation of linear and nonlinear dust acoustic waves in a homogeneous unmagnetized, collisionless and dissipative dusty plasma consisted of extremely massive, micron-sized, negative dust grains has been investigated. The Boltzmann distribution is suggested for electrons whereas vortex-like distribution for ions. In the linear analysis, the dispersion relation is obtained, and the dependence of damping rate of the waves on the carrier wave number , the dust kinematic viscosity coefficient and the ratio of the ions to the electrons temperatures is discussed. In the nonlinear analysis, the modified Korteweg-de Vries-Burgers (mKdV-Burgers) equation is derived via the reductive perturbation method. Bifurcation analysis is discussed for non-dissipative system in the absence of Burgers term. In the case of dissipative system, the tangent hyperbolic method is used to solve mKdV-Burgers equation, and yield the shock wave solution. The obtained results may be helpful in better understanding of waves propagation in the astrophysical plasmas as well as in inertial confinement fusion laboratory plasmas.
Non-linear dynamic analysis of a flexible rotor supported on porous oil journal bearings
NASA Astrophysics Data System (ADS)
Laha, S. K.; Kakoty, S. K.
2011-03-01
In the present paper, the non-linear dynamic analysis of a flexible rotor with a rigid disk under unbalance excitation mounted on porous oil journal bearings at the two ends is carried out. The system equation of motion is obtained by finite element formulation of Timoshenko beam and the disk. The non-linear oil-film forces are calculated from the solution of the modified Reynolds equation simultaneously with Darcy's equation. The system equation of motion is then solved by the Wilson- θ method. Bifurcation diagrams, Poincaré maps, time response, journal trajectories, FFT-spectrum, etc. are obtained to study the non-linear dynamics of the rotor-bearing system. The effect of various non-dimensional rotor-bearing parameters on the bifurcation characteristics of the system is studied. It is shown that the system undergoes Hopf bifurcation as the speed increases. Further, slenderness ratio, material properties of the rotor, ratio of disk mass to shaft mass and permeability of the porous bush are shown to have profound effect on the bifurcation characteristics of the rotor-bearing system.
Hadjiloucas, Sillas; Galvão, Roberto K H; Bowen, John W
2002-12-01
We provide a unified framework for a range of linear transforms that can be used for the analysis of terahertz spectroscopic data, with particular emphasis on their application to the measurement of leaf water content. The use of linear transforms for filtering, regression, and classification is discussed. For illustration, a classification problem involving leaves at three stages of drought and a prediction problem involving simulated spectra are presented. Issues resulting from scaling the data set are discussed. Using Lagrange multipliers, we arrive at the transform that yields the maximum separation between the spectra and show that this optimal transform is equivalent to computing the Euclidean distance between the samples. The optimal linear transform is compared with the average for all the spectra as well as with the Karhunen-Loève transform to discriminate a wet leaf from a dry leaf. We show that taking several principal components into account is equivalent to defining new axes in which data are to be analyzed. The procedure shows that the coefficients of the Karhunen-Loève transform are well suited to the process of classification of spectra. This is in line with expectations, as these coefficients are built from the statistical properties of the data set analyzed. PMID:12469746
Non-linear analysis and calculation of the performance of a shelving protection system by FEM
NASA Astrophysics Data System (ADS)
García Nieto, P. J.; del Coz Díaz, J. J.; Vilán Vilán, J. A.; Suárez Sierra, J. L.
2012-12-01
The aim of this paper consists on the study, analysis and calculation of the efficiency of a shelving protection system by means of the finite element method (FEM). These shelving protection systems are intended to prevent the eventual damage due to the impacts of transport elements in motion, such as: forklifts, dumpers, hand pallet trucks, and so on. The impact loads may threaten the structural integrity of the shelving system. The present structural problem is highly non-linear, due to the simultaneous presence of the following nonlinearities: material non-linearity (plasticity in this case), geometrical non-linearity (large displacements) and contact-type boundary conditions (between the rigid body and the protection system). A total of forty eight different FEM models are built varying the thickness of the steel plate (4, 5 and 6 mm), the impact height (0.1, 0.2, 0.3 and 0.4 meters) and the impact direction (head-on collision and side impact). Once the models are solved, the stress distribution, the overall displacements and the absorbed impact energy were calculated. In order to determine the best shelving protection's candidate, some constraints must be taken into account: the maximum allowable stress (235 MPa), the maximum displacement (0.05 m) and the absorbed impact energy (400 J according to the European Standard Rule PREN-15512). Finally, the most important results are shown and conclusions of this study are exposed.
A 3-D Magnetic Analysis of a Stirling Convertor Linear Alternator Under Load
NASA Technical Reports Server (NTRS)
Geng, Steven M.; Schwarze, Gene E.; Niedra, Janis M.; Regan, Timothy F.
2001-01-01
The NASA Glenn Research Center (GRC), the Department of Energy (DOE), and the Stirling Technology Company (STC) are developing Stirling convertors for Stirling Radioisotope Power Systems (SRPS) to provide electrical power for future NASA deep space missions. STC is developing the 55-We Technology Demonstration Convertor (TDC) under contract to DOE. Of critical importance to the successful development of the Stirling convertor for space power applications is the development of a lightweight and highly efficient linear alternator. This paper presents a 3-dimensional finite element method (FEM) approach for evaluating Stirling convertor linear alternators. The model extends a magnetostatic analysis previously reported at the 35th Intersociety Energy Conversion Engineering Conference (IECEC) to include the effects of the load current. STC's 55-We linear alternator design was selected to validate the model. Spatial plots of magnetic field strength (H) are presented in the region of the exciting permanent magnets. The margin for permanent magnet demagnetization is calculated at the expected magnet operating temperature for the near earth environment and for various average magnet temperatures. These thermal conditions were selected to represent a worst-case condition for the planned deep space missions. This paper presents plots that identify regions of high H where the potential to alter the magnetic moment of the magnets exists.
Dhanya, S; Kumari Roshni, V S
2016-01-01
Textures play an important role in image classification. This paper proposes a high performance texture classification method using a combination of multiresolution analysis tool and linear regression modelling by channel elimination. The correlation between different frequency regions has been validated as a sort of effective texture characteristic. This method is motivated by the observation that there exists a distinctive correlation between the image samples belonging to the same kind of texture, at different frequency regions obtained by a wavelet transform. Experimentally, it is observed that this correlation differs across textures. The linear regression modelling is employed to analyze this correlation and extract texture features that characterize the samples. Our method considers not only the frequency regions but also the correlation between these regions. This paper primarily focuses on applying the Dual Tree Complex Wavelet Packet Transform and the Linear Regression model for classification of the obtained texture features. Additionally the paper also presents a comparative assessment of the classification results obtained from the above method with two more types of wavelet transform methods namely the Discrete Wavelet Transform and the Discrete Wavelet Packet Transform. PMID:26835234
Development of a linearized unsteady aerodynamic analysis for cascade gust response predictions
NASA Technical Reports Server (NTRS)
Verdon, Joseph M.; Hall, Kenneth C.
1990-01-01
A method for predicting the unsteady aerodynamic response of a cascade of airfoils to entropic, vortical, and acoustic gust excitations is being developed. Here, the unsteady flow is regarded as a small perturbation of a nonuniform isentropic and irrotational steady background flow. A splitting technique is used to decompose the linearized unsteady velocity into rotational and irrotational parts leading to equations for the complex amplitudes of the linearized unsteady entropy, rotational velocity, and velocity potential that are coupled only sequentially. The entropic and rotational velocity fluctuations are described by transport equations for which closed-form solutions in terms of the mean-flow drift and stream functions can be determined. The potential fluctuation is described by an inhomogeneous convected wave equation in which the source term depends on the rotational velocity field, and is determined using finite-difference procedures. The analytical and numerical techniques used to determine the linearized unsteady flow are outlined. Results are presented to indicate the status of the solution procedure and to demonstrate the impact of blade geometry and mean blade loading on the aerodynamic response of cascades to vortical gust excitations. The analysis described herein leads to very efficient predictions of cascade unsteady aerodynamic response phenomena making it useful for turbomachinery aeroelastic and aeroacoustic design applications.
A spectral analysis of the domain decomposed Monte Carlo method for linear systems
Slattery, Stuart R.; Evans, Thomas M.; Wilson, Paul P. H.
2015-09-08
The domain decomposed behavior of the adjoint Neumann-Ulam Monte Carlo method for solving linear systems is analyzed using the spectral properties of the linear oper- ator. Relationships for the average length of the adjoint random walks, a measure of convergence speed and serial performance, are made with respect to the eigenvalues of the linear operator. In addition, relationships for the effective optical thickness of a domain in the decomposition are presented based on the spectral analysis and diffusion theory. Using the effective optical thickness, the Wigner rational approxi- mation and the mean chord approximation are applied to estimate the leakage frac- tion of random walks from a domain in the decomposition as a measure of parallel performance and potential communication costs. The one-speed, two-dimensional neutron diffusion equation is used as a model problem in numerical experiments to test the models for symmetric operators with spectral qualities similar to light water reactor problems. We find, in general, the derived approximations show good agreement with random walk lengths and leakage fractions computed by the numerical experiments.
A spectral analysis of the domain decomposed Monte Carlo method for linear systems
Slattery, Stuart R.; Evans, Thomas M.; Wilson, Paul P. H.
2015-09-08
The domain decomposed behavior of the adjoint Neumann-Ulam Monte Carlo method for solving linear systems is analyzed using the spectral properties of the linear oper- ator. Relationships for the average length of the adjoint random walks, a measure of convergence speed and serial performance, are made with respect to the eigenvalues of the linear operator. In addition, relationships for the effective optical thickness of a domain in the decomposition are presented based on the spectral analysis and diffusion theory. Using the effective optical thickness, the Wigner rational approxi- mation and the mean chord approximation are applied to estimate the leakagemore » frac- tion of random walks from a domain in the decomposition as a measure of parallel performance and potential communication costs. The one-speed, two-dimensional neutron diffusion equation is used as a model problem in numerical experiments to test the models for symmetric operators with spectral qualities similar to light water reactor problems. We find, in general, the derived approximations show good agreement with random walk lengths and leakage fractions computed by the numerical experiments.« less
Linear Analysis of the m=0 Instability for a Visco-resistive Hall MHD Plasma
Oliver, B.V.; Genoni, T.C.; Mehlhorn, T.A.
2006-01-05
We present a comprehensive analysis of the linear dispersion relation for the axisymmetric (m=0) compressible interchange instability of Bennett equilibria in a visco-resistive, Hall MHD plasma. The full anisotropic stress tensor with Braginskii viscous coefficients is considered. The eigenvalues are obtained numerically. For small axial mode number, Hall currents enhance the growth rates, whereas in the limit of high mode number the growth rates are suppressed, eventually resulting in wave cutoff, even in the ideal limit. For the visco-resistive plasma the unstable spectra are weakly dependent on the off-diagonal elements of the stress tensor.
NASA Astrophysics Data System (ADS)
Brake, M. R.
2011-06-01
The analysis of continuous systems with piecewise-linear constraints in their domains have previously been limited to either numerical approaches, or analytical methods that are constrained in the parameter space, boundary conditions, or order of the system. The present analysis develops a robust method for studying continuous systems with arbitrary boundary conditions and discrete piecewise-linear constraints. A superposition method is used to generate homogeneous boundary conditions, and modal analysis is used to find the displacement of the system in each state of the piecewise-linear constraint. In order to develop a mapping across each slope discontinuity in the piecewise-linear force-deflection profile, a variational calculus approach is taken that minimizes the L 2 energy norm between the previous and current states. An approach for calculating the finite-time Lyapunov exponents is presented in order to determine chaotic regimes. To illustrate this method, two examples are presented: a pinned-pinned beam with a deadband constraint, and a leaf spring coupled with a connector pin immersed in a viscous fluid. The pinned-pinned beam example illustrates the method for a non-operator based analysis. Results are used to show that the present method does not necessitate the need of a large number of basis functions to adequately map the displacement and velocity of the system across states. In the second example, the leaf spring is modeled as a clamped-free beam. The interaction between the beam and the connector pin is modeled with a preload and a penalty stiffness. Several experiments are conducted in order to validate aspects of the leaf spring model. From the results of the convergence and parameter studies, a high correlation between the finite-time Lyapunov exponents and the contact time per period of the excitation is observed. The parameter studies also indicate that when the system's parameters are changed in order to reduce the magnitude of the impact
Some non-linear interactions in polytropic gas cosmology: phase space analysis
NASA Astrophysics Data System (ADS)
Khurshudyan, Martiros
2015-11-01
There are various cosmological models with polytropic equation of state associated to dark energy. Polytropic EoS has important applications in astrophysics, therefore a study of it on cosmological framework continues to be interesting. From the other hand, there are various forms of interactions phenomenologically involved into the darkness of the universe able to solve important cosmological problems. This is a motivation for us to perform a phase space analysis of various cosmological scenarios where non-linear interacting polytropic gas models are involved. Dark matter is taken to be a pressureless fluid.
Statistical mechanical analysis of the linear vector channel in digital communication
NASA Astrophysics Data System (ADS)
Takeda, Koujin; Hatabu, Atsushi; Kabashima, Yoshiyuki
2007-11-01
A statistical mechanical framework to analyze linear vector channel models in digital wireless communication is proposed for a large system. The framework is a generalization of that proposed for code-division multiple-access systems in Takeda et al (2006 Europhys. Lett. 76 1193) and enables the analysis of the system in which the elements of the channel transfer matrix are statistically correlated with each other. The significance of the proposed scheme is demonstrated by assessing the performance of an existing model of multi-input multi-output communication systems.
NASA Technical Reports Server (NTRS)
Chun, E.; Rosner, R.
1993-01-01
We study the linear stability of an optically thin uniform radiating plasma subject to nonlocal heat transport. We derive the dispersion relation appropriate to this problem, and the marginal wavenumbers for instability. Our analysis indicates that nonlocal heat transport acts to reduce the stabilizing influence of thermal conduction, and that there are critical values for the electron mean free path such that the plasma is always unstable. Our results may be applied to a number of astrophysical plasmas, one such example being the halos of clusters of galaxies.
Finite element analysis of direct thrust-controlled linear induction motor
Kwon, B.I.; Woo, K.I.; Kim, S. . Dept. of Electrical Engineering)
1999-05-01
This paper describes the finite element analysis of a direct thrust-controlled linear induction motor (LIM). The time-stepping finite element method and the moving mesh technique are used to calculate the dynamic characteristics of LIM during the direct thrust control. Because LIM has the end effect, thrust correction coefficient is introduced to predict an actual thrust in control. The simulation results, the thrust and the stator flux linkage are shown below and the stator current is compared with an experimental one.
Jiang, Zhaolin; Shen, Nuo; Zhou, Jianwei
2013-01-01
We first give the style spectral decomposition of a special skew circulant matrix C and then get the style decomposition of arbitrary skew circulant matrix by making use of the Kronecker products between the elements of first row in skew circulant and the special skew circulant C. Besides that, we obtain the singular value of skew circulant matrix as well. Finally, we deal with the optimal backward perturbation analysis for the linear system with skew circulant coefficient matrix on the base of its style spectral decomposition. PMID:24369488
Modeling and analysis of aircraft non-linear components for harmonics analysis
Karimi, K.J.; Voss, J.
1995-12-31
Modern commercial aircraft Electric Power Systems (EPS) include many nonlinear components which produce harmonics. The addition of all the current harmonics could result in a power system with unacceptable levels of voltage distortion. It is important to be able to predict the levels of voltage distortion at early program stages to correct any potential problems and avoid costly redesigns. In this paper the nature and sources of harmonic producing equipment are described. These sources of harmonics and their effect on aircraft power system operation are described. Models for various aircraft non-linear components are developed in this paper. These component models are used in a model of the Boeing 777 EPS which is used to calculate voltage harmonics for various airplane configurations and flight conditions. A description of this model and the models used for various components are given. Tests performed to validate these models are described. Comparison of experimental results with analytical model predictions are given.
The linear statistical d.c. model of GaAs MESFET using factor analysis
NASA Astrophysics Data System (ADS)
Dobrzanski, Lech
1995-02-01
The linear statistical model of the GaAs MESFET's current generator is obtained by means of factor analysis. Three different MESFET deterministic models are taken into account in the analysis: the Statz model (ST), the Materka-type model (MT) and a new proprietary model of MESFET with implanted channel (PLD). It is shown that statistical models obtained using factor analysis provide excellent generation of the multidimensional random variable representing the drain current of MESFET. The method of implementation of the statistical model into the SPICE program is presented. It is proved that for a strongly limited number of Monte Carlo analysis runs in that program, the statistical models considered in each case (ST, MT and PLD) enable good reconstruction of the empirical factor structure. The empirical correlation matrix of model parameters is not reconstructed exactly by statistical modelling, but values of correlation matrix elements obtained from simulated data are within the confidence intervals for the small sample. This paper proves that a formal approach to statistical modelling using factor analysis is the right path to follow, in spite of the fact, that CAD systems (PSpice[MicroSim Corp.], Microwave Harmonica[Compact Software]) are not designed properly for generation of the multidimensional random variable. It is obvious that further progress in implementation of statistical methods in CAD software is required. Furthermore, a new approach to the MESFET's d.c. model is presented. The separate functions, describing the linear as well as the saturated region of MESFET output characteristics, are combined in the single equation. This way of modelling is particularly suitable for transistors with an implanted channel.
Ma, Rongfei
2015-01-01
In this paper, ammonia quantitative analysis based on miniaturized Al ionization gas sensor and non-linear bistable dynamic model was proposed. Al plate anodic gas-ionization sensor was used to obtain the current-voltage (I-V) data. Measurement data was processed by non-linear bistable dynamics model. Results showed that the proposed method quantitatively determined ammonia concentrations. PMID:25975362
PLATSIM: An efficient linear simulation and analysis package for large-order flexible systems
NASA Technical Reports Server (NTRS)
Maghami, Periman; Kenny, Sean P.; Giesy, Daniel P.
1995-01-01
PLATSIM is a software package designed to provide efficient time and frequency domain analysis of large-order generic space platforms implemented with any linear time-invariant control system. Time domain analysis provides simulations of the overall spacecraft response levels due to either onboard or external disturbances. The time domain results can then be processed by the jitter analysis module to assess the spacecraft's pointing performance in a computationally efficient manner. The resulting jitter analysis algorithms have produced an increase in speed of several orders of magnitude over the brute force approach of sweeping minima and maxima. Frequency domain analysis produces frequency response functions for uncontrolled and controlled platform configurations. The latter represents an enabling technology for large-order flexible systems. PLATSIM uses a sparse matrix formulation for the spacecraft dynamics model which makes both the time and frequency domain operations quite efficient, particularly when a large number of modes are required to capture the true dynamics of the spacecraft. The package is written in MATLAB script language. A graphical user interface (GUI) is included in the PLATSIM software package. This GUI uses MATLAB's Handle graphics to provide a convenient way for setting simulation and analysis parameters.
Hopfgartner, Gérard; Varesio, Emmanuel; Tschäppät, Viviane; Grivet, Chantal; Bourgogne, Emmanuel; Leuthold, Luc Alexis
2004-08-01
Recently, linear ion traps (LITs) have been combined with quadrupole (Q), time-of-flight (TOF) and Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry (MS). LITs can be used either as ion accumulation devices or as commercially available, stand-alone mass spectrometers with MSn capabilities. The combination of triple quadrupole MS with LIT technology in the form of an instrument of configuration QqLIT, using axial ejection, is particularly interesting, because this instrument retains the classical triple quadrupole scan functions such as selected reaction monitoring (SRM), product ion (PI), neutral loss (NL) and precursor ion (PC) while also providing access to sensitive ion trap experiments. For small molecules, quantitative and qualitative analysis can be performed using the same instrument. In addition, for peptide analysis, the enhanced multiply charged (EMC) scan allows an increase in selectivity, while the time-delayed fragmentation (TDF) scan provides additional structural information. Various methods of operating the hybrid instrument are described for the case of the commercial Q TRAP (AB/MDS Sciex) and applications to drug metabolism analysis, quantitative confirmatory analysis, peptides analysis and automated nanoelectrospray (ESI-chip-MS) analysis are discussed. PMID:15329837
Analysis of non linear partially standing waves from 3D velocity measurements
NASA Astrophysics Data System (ADS)
Drevard, D.; Rey, V.; Svendsen, Ib; Fraunie, P.
2003-04-01
Surface gravity waves in the ocean exhibit an energy spectrum distributed in both frequency and direction of propagation. Wave data collection is of great importance in coastal zones for engineering and scientific studies. In particular, partially standing waves measurements near coastal structures and steep or barred beaches may be a requirement, for instance for morphodynamic studies. The aim of the present study is the analysis of partially standing surface waves icluding non-linear effects. According to 1st order Stokes theory, synchronous measurements of horizontal and vertical velocity components allow calculation of rate of standing waves (Drevard et al, 2003). In the present study, it is demonstrated that for deep water conditions, partially standing 2nd order Stokes waves induced velocity field is still represented by the 1st order solution for the velocity potential contrary to the surface elevation which exhibits harmonic components. For intermediate water depth, harmonic components appear not only in the surface elevation but also in the velocity fields, but their weight remains much smaller, because of the vertical decreasing wave induced motion. For irregular waves, the influence of the spectrum width on the non-linear effects in the analysis is discussed. Keywords: Wave measurements ; reflection ; non-linear effects Acknowledgements: This work was initiated during the stay of Prof. Ib Svendsen, as invited Professor, at LSEET in autumn 2002. This study is carried out in the framework of the Scientific French National Programmes PNEC ART7 and PATOM. Their financial supports are acknowledged References: Drevard, D., Meuret, A., Rey, V. Piazzola, J. And Dolle, A.. (2002). "Partially reflected waves measurements using Acoustic Doppler Velocimeter (ADV)", Submitted to ISOPE 03, Honolulu, Hawaii, May 2003.
NASA Astrophysics Data System (ADS)
Choudhury, Prakriti Pal; Sharma, Prateek
2016-04-01
We perform global linear stability analysis and idealized numerical simulations in global thermal balance to understand the condensation of cold gas from hot/virial atmospheres (coronae), in particular the intracluster medium (ICM). We pay particular attention to geometry (e.g. spherical versus plane-parallel) and the nature of the gravitational potential. Global linear analysis gives a similar value for the fastest growing thermal instability modes in spherical and Cartesian geometries. Simulations and observations suggest that cooling in haloes critically depends on the ratio of the cooling time to the free-fall time (tcool/tff). Extended cold gas condenses out of the ICM only if this ratio is smaller than a threshold value close to 10. Previous works highlighted the difference between the nature of cold gas condensation in spherical and plane-parallel atmospheres; namely, cold gas condensation appeared easier in spherical atmospheres. This apparent difference due to geometry arises because the previous plane-parallel simulations focused on in situ condensation of multiphase gas but spherical simulations studied condensation anywhere in the box. Unlike previous claims, our non-linear simulations show that there are only minor differences in cold gas condensation, either in situ or anywhere, for different geometries. The amount of cold gas depends on the shape of tcool/tff; gas has more time to condense if gravitational acceleration decreases towards the centre. In our idealized plane-parallel simulations with heating balancing cooling in each layer, there can be significant mass/energy/momentum transfer across layers that can trigger condensation and drive tcool/tff far beyond the critical value close to 10.
NASA Technical Reports Server (NTRS)
Slater, G. L.; Shelley, Stuart; Jacobson, Mark
1993-01-01
In this paper, the design, analysis, and test of a low cost, linear proof mass actuator for vibration control is presented. The actuator is based on a linear induction coil from a large computer disk drive. Such disk drives are readily available and provide the linear actuator, current feedback amplifier, and power supply for a highly effective, yet inexpensive, experimental laboratory actuator. The device is implemented as a force command input system, and the performance is virtually the same as other, more sophisticated, linear proof mass systems.
Non-linear analysis and the design of Pumpkin Balloons: stress, stability and viscoelasticity
NASA Astrophysics Data System (ADS)
Rand, J. L.; Wakefield, D. S.
Tensys have a long-established background in the shape generation and load analysis of architectural stressed membrane structures Founded upon their inTENS finite element analysis suite these activities have broadened to encompass lighter than air structures such as aerostats hybrid air-vehicles and stratospheric balloons Winzen Engineering couple many years of practical balloon design and fabrication experience with both academic and practical knowledge of the characterisation of the non-linear viscoelastic response of the polymeric films typically used for high-altitude scientific balloons Both companies have provided consulting services to the NASA Ultra Long Duration Balloon ULDB Program Early implementations of pumpkin balloons have shown problems of geometric instability characterised by improper deployment and these difficulties have been reproduced numerically using inTENS The solution lies in both the shapes of the membrane lobes and also the need to generate a biaxial stress field in order to mobilise in-plane shear stiffness Balloons undergo significant temperature and pressure variations in flight The different thermal characteristics between tendons and film can lead to significant meridional stress Fabrication tolerances can lead to significant local hoop stress concentrations particularly adjacent to the base and apex end fittings The non-linear viscoelastic response of the envelope film acts positively to help dissipate stress concentrations However creep over time may produce lobe geometry variations that may
Factor Analysis of Linear Type Traits and Their Relation with Longevity in Brazilian Holstein Cattle
Kern, Elisandra Lurdes; Cobuci, Jaime Araújo; Costa, Cláudio Napolis; Pimentel, Concepta Margaret McManus
2014-01-01
In this study we aimed to evaluate the reduction in dimensionality of 20 linear type traits and more final score in 14,943 Holstein cows in Brazil using factor analysis, and indicate their relationship with longevity and 305 d first lactation milk production. Low partial correlations (−0.19 to 0.38), the medium to high Kaiser sampling mean (0.79) and the significance of the Bartlett sphericity test (p<0.001), indicated correlations between type traits and the suitability of these data for a factor analysis, after the elimination of seven traits. Two factors had autovalues greater than one. The first included width and height of posterior udder, udder texture, udder cleft, loin strength, bone quality and final score. The second included stature, top line, chest width, body depth, fore udder attachment, angularity and final score. The linear regression of the factors on several measures of longevity and 305 d milk production showed that selection considering only the first factor should lead to improvements in longevity and 305 milk production. PMID:25050015
A linear stability analysis for nonlinear, grey, thermal radiative transfer problems
Wollaber, Allan B.; Larsen, Edward W.
2011-02-20
We present a new linear stability analysis of three time discretizations and Monte Carlo interpretations of the nonlinear, grey thermal radiative transfer (TRT) equations: the widely used 'Implicit Monte Carlo' (IMC) equations, the Carter Forest (CF) equations, and the Ahrens-Larsen or 'Semi-Analog Monte Carlo' (SMC) equations. Using a spatial Fourier analysis of the 1-D Implicit Monte Carlo (IMC) equations that are linearized about an equilibrium solution, we show that the IMC equations are unconditionally stable (undamped perturbations do not exist) if {alpha}, the IMC time-discretization parameter, satisfies 0.5 < {alpha} {<=} 1. This is consistent with conventional wisdom. However, we also show that for sufficiently large time steps, unphysical damped oscillations can exist that correspond to the lowest-frequency Fourier modes. After numerically confirming this result, we develop a method to assess the stability of any time discretization of the 0-D, nonlinear, grey, thermal radiative transfer problem. Subsequent analyses of the CF and SMC methods then demonstrate that the CF method is unconditionally stable and monotonic, but the SMC method is conditionally stable and permits unphysical oscillatory solutions that can prevent it from reaching equilibrium. This stability theory provides new conditions on the time step to guarantee monotonicity of the IMC solution, although they are likely too conservative to be used in practice. Theoretical predictions are tested and confirmed with numerical experiments.
NASA Astrophysics Data System (ADS)
Dunstan, Jocelyn; Lee, Kyoung Jin; Park, Simon; Goldstein, Raymond E.
A novel form of convection was observed in a suspension of non-motile Photobacterium phosphoreum bacteria. The pattern resembles classical bioconvection, however this strain has limited if any motility, which excludes this possible explanation. After performing a series of control experiments we found that the convection was actually driven by the evaporation of the salty bacterial medium, and the same kind of plumes were observed using polystyrene beads suspended in water with salt added. A mathematical model was formulated for the process and studied using a linear stability analysis and finite element method simulations, reproducing most of the observed experimental features. From the linear stability analysis, a threshold in salt concentration to observe convective motion was obtained, as well as the wavelength of the pattern at the onset of the instability. This was complemented by finite element simulations, which produced plume dynamics remarkably similar to the experimental observations. Evaporation-driven convection on the millimeter scale has not been studied extensively, and its effect may have been underestimated in other experiments.
The Analysis and Construction of Perfectly Matched Layers for the Linearized Euler Equations
NASA Technical Reports Server (NTRS)
Hesthaven, J. S.
1997-01-01
We present a detailed analysis of a recently proposed perfectly matched layer (PML) method for the absorption of acoustic waves. The split set of equations is shown to be only weakly well-posed, and ill-posed under small low order perturbations. This analysis provides the explanation for the stability problems associated with the split field formulation and illustrates why applying a filter has a stabilizing effect. Utilizing recent results obtained within the context of electromagnetics, we develop strongly well-posed absorbing layers for the linearized Euler equations. The schemes are shown to be perfectly absorbing independent of frequency and angle of incidence of the wave in the case of a non-convecting mean flow. In the general case of a convecting mean flow, a number of techniques is combined to obtain a absorbing layers exhibiting PML-like behavior. The efficacy of the proposed absorbing layers is illustrated though computation of benchmark problems in aero-acoustics.
Lee, S. H.; van der Werf, J. H. J.
2016-01-01
Summary: We have developed an algorithm for genetic analysis of complex traits using genome-wide SNPs in a linear mixed model framework. Compared to current standard REML software based on the mixed model equation, our method is substantially faster. The advantage is largest when there is only a single genetic covariance structure. The method is particularly useful for multivariate analysis, including multi-trait models and random regression models for studying reaction norms. We applied our proposed method to publicly available mice and human data and discuss the advantages and limitations. Availability and implementation: MTG2 is available in https://sites.google.com/site/honglee0707/mtg2. Contact: hong.lee@une.edu.au Supplementary information: Supplementary data are available at Bioinformatics online. PMID:26755623
Linear Closed-form Solution and Finite-element Analysis of an Active Tensegrity Unit
NASA Astrophysics Data System (ADS)
Kmeť, Stanislav; Platko, Peter
2012-11-01
Results of the linear closed form solution of an active or adaptive tensegrity unit, as well as its numerical analysis using finite element method are presented in the paper. The shape of the unit is an octahedral cell with a square base and it is formed by thirteen members (four bottom and four top cables, four edge struts and one central strut). The central strut is designed as an actuator that allows for an adjustment of the shape of the unit which leads to changes of tensile forces in the cables. Due to the diagonal symmetry of the 3D tensegrity unit the closed-form analysis is based on the 2D solution of the equivalent planar biconvex cable system with one central strut under a vertical point load.
NASA Astrophysics Data System (ADS)
Iwaoka, Nobuyuki; Hagita, Katsumi; Takano, Hiroshi
2014-03-01
On the basis of relaxation mode analysis (RMA), we present an efficient method to estimate the linear viscoelasticity of polymer melts in a molecular dynamics (MD) simulation. Slow relaxation phenomena appeared in polymer melts cause a problem that a calculation of the stress relaxation function in MD simulations, especially in the terminal time region, requires large computational efforts. Relaxation mode analysis is a method that systematically extracts slow relaxation modes and rates of the polymer chain from the time correlation of its conformations. We show the computational cost may be drastically reduced by combining a direct calculation of the stress relaxation function based on the Green-Kubo formula with the relaxation rates spectra estimated by RMA. N. I. acknowledges the Graduate School Doctoral Student Aid Program from Keio University.
Analysis of a tubular linear motor with soft magnetic composites for reciprocating compressors
NASA Astrophysics Data System (ADS)
Hsieh, Min-Fu; Hu, Kai-Hsiang
2008-04-01
This paper presents the analysis of a single-phase tubular linear motor formed with soft magnetic composites for direct driving reciprocating compressors. This direct-drive design can solve the potential cylinder wearing problem caused by the crank side force in conventional reciprocating compressors. The overall size may also be reduced without transmission mechanism. To produce sufficient thrust for the cooling cycles, the designed motor features a structure of three-dimensional flux distributions that allows an increase of slot-fill factor to enhance the thrust density. Moreover, the motor makes use of the inherent cogging force which contributes to the overall thrust. Finite element analysis is employed to verify the performance. The results show that the motor is capable of producing high thrust with a compact size compared to other types of motors.
NASA Astrophysics Data System (ADS)
Eguizabal, Alma; Laughney, Ashley M.; Garcia Allende, Pilar Beatriz; Krishnaswamy, Venkataramanan; Wells, Wendy A.; Paulsen, Keith D.; Pogue, Brian W.; López-Higuera, José M.; Conde, Olga M.
2013-02-01
Texture analysis of light scattering in tissue is proposed to obtain diagnostic information from breast cancer specimens. Light scattering measurements are minimally invasive, and allow the estimation of tissue morphology to guide the surgeon in resection surgeries. The usability of scatter signatures acquired with a micro-sampling reflectance spectral imaging system was improved utilizing an empirical approximation to the Mie theory to estimate the scattering power on a per-pixel basis. Co-occurrence analysis is then applied to the scattering power images to extract the textural features. A statistical analysis of the features demonstrated the suitability of the autocorrelation for the classification of notmalignant (normal epithelia and stroma, benign epithelia and stroma, inflammation), malignant (DCIS, IDC, ILC) and adipose tissue, since it reveals morphological information of tissue. Non-malignant tissue shows higher autocorrelation values while adipose tissue presents a very low autocorrelation on its scatter texture, being malignant the middle ground. Consequently, a fast linear classifier based on the consideration of just one straightforward feature is enough for providing relevant diagnostic information. A leave-one-out validation of the linear classifier on 29 samples with 48 regions of interest showed classification accuracies of 98.74% on adipose tissue, 82.67% on non-malignant tissue and 72.37% on malignant tissue, in comparison with the biopsy H and E gold standard. This demonstrates that autocorrelation analysis of scatter signatures is a very computationally efficient and automated approach to provide pathological information in real-time to guide surgeon during tissue resection.
Design and analysis of an unconventional permanent magnet linear machine for energy harvesting
NASA Astrophysics Data System (ADS)
Zeng, Peng
This Ph.D. dissertation proposes an unconventional high power density linear electromagnetic kinetic energy harvester, and a high-performance two-stage interface power electronics to maintain maximum power abstraction from the energy source and charge the Li-ion battery load with constant current. The proposed machine architecture is composed of a double-sided flat type silicon steel stator with winding slots, a permanent magnet mover, coil windings, a linear motion guide and an adjustable spring bearing. The unconventional design of the machine is that NdFeB magnet bars in the mover are placed with magnetic fields in horizontal direction instead of vertical direction and the same magnetic poles are facing each other. The derived magnetic equivalent circuit model proves the average air-gap flux density of the novel topology is as high as 0.73 T with 17.7% improvement over that of the conventional topology at the given geometric dimensions of the proof-of-concept machine. Subsequently, the improved output voltage and power are achieved. The dynamic model of the linear generator is also developed, and the analytical equations of output maximum power are derived for the case of driving vibration with amplitude that is equal, smaller and larger than the relative displacement between the mover and the stator of the machine respectively. Furthermore, the finite element analysis (FEA) model has been simulated to prove the derived analytical results and the improved power generation capability. Also, an optimization framework is explored to extend to the multi-Degree-of-Freedom (n-DOF) vibration based linear energy harvesting devices. Moreover, a boost-buck cascaded switch mode converter with current controller is designed to extract the maximum power from the harvester and charge the Li-ion battery with trickle current. Meanwhile, a maximum power point tracking (MPPT) algorithm is proposed and optimized for low frequency driving vibrations. Finally, a proof
Discriminative analysis of non-linear brain connectivity for leukoaraiosis with resting-state fMRI
NASA Astrophysics Data System (ADS)
Lai, Youzhi; Xu, Lele; Yao, Li; Wu, Xia
2015-03-01
Leukoaraiosis (LA) describes diffuse white matter abnormalities on CT or MR brain scans, often seen in the normal elderly and in association with vascular risk factors such as hypertension, or in the context of cognitive impairment. The mechanism of cognitive dysfunction is still unclear. The recent clinical studies have revealed that the severity of LA was not corresponding to the cognitive level, and functional connectivity analysis is an appropriate method to detect the relation between LA and cognitive decline. However, existing functional connectivity analyses of LA have been mostly limited to linear associations. In this investigation, a novel measure utilizing the extended maximal information coefficient (eMIC) was applied to construct non-linear functional connectivity in 44 LA subjects (9 dementia, 25 mild cognitive impairment (MCI) and 10 cognitively normal (CN)). The strength of non-linear functional connections for the first 1% of discriminative power increased in MCI compared with CN and dementia, which was opposed to its linear counterpart. Further functional network analysis revealed that the changes of the non-linear and linear connectivity have similar but not completely the same spatial distribution in human brain. In the multivariate pattern analysis with multiple classifiers, the non-linear functional connectivity mostly identified dementia, MCI and CN from LA with a relatively higher accuracy rate than the linear measure. Our findings revealed the non-linear functional connectivity provided useful discriminative power in classification of LA, and the spatial distributed changes between the non-linear and linear measure may indicate the underlying mechanism of cognitive dysfunction in LA.
ANALYSIS OF QUIET-SUN INTERNETWORK MAGNETIC FIELDS BASED ON LINEAR POLARIZATION SIGNALS
Orozco Suarez, D.; Bellot Rubio, L. R.
2012-05-20
We present results from the analysis of Fe I 630 nm measurements of the quiet Sun taken with the spectropolarimeter of the Hinode satellite. Two data sets with noise levels of 1.2 Multiplication-Sign 10{sup -3} and 3 Multiplication-Sign 10{sup -4} are employed. We determine the distribution of field strengths and inclinations by inverting the two observations with a Milne-Eddington model atmosphere. The inversions show a predominance of weak, highly inclined fields. By means of several tests we conclude that these properties cannot be attributed to photon noise effects. To obtain the most accurate results, we focus on the 27.4% of the pixels in the second data set that have linear polarization amplitudes larger than 4.5 times the noise level. The vector magnetic field derived for these pixels is very precise because both circular and linear polarization signals are used simultaneously. The inferred field strength, inclination, and filling factor distributions agree with previous results, supporting the idea that internetwork (IN) fields are weak and very inclined, at least in about one quarter of the area occupied by the IN. These properties differ from those of network fields. The average magnetic flux density and the mean field strength derived from the 27.4% of the field of view with clear linear polarization signals are 16.3 Mx cm{sup -2} and 220 G, respectively. The ratio between the average horizontal and vertical components of the field is approximately 3.1. The IN fields do not follow an isotropic distribution of orientations.
Linear Spectral Analysis of Plume Emissions Using an Optical Matrix Processor
NASA Technical Reports Server (NTRS)
Gary, C. K.
1992-01-01
Plume spectrometry provides a means to monitor the health of a burning rocket engine, and optical matrix processors provide a means to analyze the plume spectra in real time. By observing the spectrum of the exhaust plume of a rocket engine, researchers have detected anomalous behavior of the engine and have even determined the failure of some equipment before it would normally have been noticed. The spectrum of the plume is analyzed by isolating information in the spectrum about the various materials present to estimate what materials are being burned in the engine. Scientists at the Marshall Space Flight Center (MSFC) have implemented a high resolution spectrometer to discriminate the spectral peaks of the many species present in the plume. Researchers at the Stennis Space Center Demonstration Testbed Facility (DTF) have implemented a high resolution spectrometer observing a 1200-lb. thrust engine. At this facility, known concentrations of contaminants can be introduced into the burn, allowing for the confirmation of diagnostic algorithms. While the high resolution of the measured spectra has allowed greatly increased insight into the functioning of the engine, the large data flows generated limit the ability to perform real-time processing. The use of an optical matrix processor and the linear analysis technique described below may allow for the detailed real-time analysis of the engine's health. A small optical matrix processor can perform the required mathematical analysis both quicker and with less energy than a large electronic computer dedicated to the same spectral analysis routine.
Indoor calibration of Sky Quality Meters: Linearity, spectral responsivity and uncertainty analysis
NASA Astrophysics Data System (ADS)
Pravettoni, M.; Strepparava, D.; Cereghetti, N.; Klett, S.; Andretta, M.; Steiger, M.
2016-09-01
The indoor calibration of brightness sensors requires extremely low values of irradiance in the most accurate and reproducible way. In this work the testing equipment of an ISO 17025 accredited laboratory for electrical testing, qualification and type approval of solar photovoltaic modules was modified in order to test the linearity of the instruments from few mW/cm2 down to fractions of nW/cm2, corresponding to levels of simulated brightness from 6 to 19 mag/arcsec2. Sixteen Sky Quality Meter (SQM) produced by Unihedron, a Canadian manufacturer, were tested, also assessing the impact of the ageing of their protective glasses on the calibration coefficients and the drift of the instruments. The instruments are in operation on measurement points and observatories at different sites and altitudes in Southern Switzerland, within the framework of OASI, the Environmental Observatory of Southern Switzerland. The authors present the results of the calibration campaign: linearity; brightness calibration, with and without protective glasses; transmittance measurement of the glasses; and spectral responsivity of the devices. A detailed uncertainty analysis is also provided, according to the ISO 17025 standard.
Spatial Analysis Based Health and Safety Risk Assessment for Linear Construction Projects
NASA Astrophysics Data System (ADS)
Atay, H.; Toz, G.
2012-07-01
This paper describes an on-going study that aims to develop a web-based spatial decision support system model for proactive health and safety management in linear construction projects. Currently, health and safety management is usually performed reactively instead of proactive management since hazard identification and risk assessment is mostly performed on paper based documents that are not effectively used at site. This leads to accidents and fatalities at construction sites. The proposed system automatically identifies the spatial risks according to the topographic and layout map of the site, project specification and health and safety regulations by means of spatial analysis. It enables the workers and management personnel to access the possible hazards and thematic risk map of any portion of the construction site for linear projects. Finally, the described approach provides the proposed mitigation measures for the identified hazards. The developed system is expected to raise awareness in H&S among workers and engineers, and increase participation of workers to health and safety management.