Science.gov

Sample records for nonlinear space decomposition

  1. Some nonlinear space decomposition algorithms

    SciTech Connect

    Tai, Xue-Cheng; Espedal, M.

    1996-12-31

    Convergence of a space decomposition method is proved for a general convex programming problem. The space decomposition refers to methods that decompose a space into sums of subspaces, which could be a domain decomposition or a multigrid method for partial differential equations. Two algorithms are proposed. Both can be used for linear as well as nonlinear elliptic problems and they reduce to the standard additive and multiplicative Schwarz methods for linear elliptic problems. Two {open_quotes}hybrid{close_quotes} algorithms are also presented. They converge faster than the additive one and have better parallelism than the multiplicative method. Numerical tests with a two level domain decomposition for linear, nonlinear and interface elliptic problems are presented for the proposed algorithms.

  2. Generalized spectral decomposition for stochastic nonlinear problems

    SciTech Connect

    Nouy, Anthony Le Maitre, Olivier P.

    2009-01-10

    We present an extension of the generalized spectral decomposition method for the resolution of nonlinear stochastic problems. The method consists in the construction of a reduced basis approximation of the Galerkin solution and is independent of the stochastic discretization selected (polynomial chaos, stochastic multi-element or multi-wavelets). Two algorithms are proposed for the sequential construction of the successive generalized spectral modes. They involve decoupled resolutions of a series of deterministic and low-dimensional stochastic problems. Compared to the classical Galerkin method, the algorithms allow for significant computational savings and require minor adaptations of the deterministic codes. The methodology is detailed and tested on two model problems, the one-dimensional steady viscous Burgers equation and a two-dimensional nonlinear diffusion problem. These examples demonstrate the effectiveness of the proposed algorithms which exhibit convergence rates with the number of modes essentially dependent on the spectrum of the stochastic solution but independent of the dimension of the stochastic approximation space.

  3. Nonlinear mode decomposition: A noise-robust, adaptive decomposition method

    NASA Astrophysics Data System (ADS)

    Iatsenko, Dmytro; McClintock, Peter V. E.; Stefanovska, Aneta

    2015-09-01

    The signals emanating from complex systems are usually composed of a mixture of different oscillations which, for a reliable analysis, should be separated from each other and from the inevitable background of noise. Here we introduce an adaptive decomposition tool—nonlinear mode decomposition (NMD)—which decomposes a given signal into a set of physically meaningful oscillations for any wave form, simultaneously removing the noise. NMD is based on the powerful combination of time-frequency analysis techniques—which, together with the adaptive choice of their parameters, make it extremely noise robust—and surrogate data tests used to identify interdependent oscillations and to distinguish deterministic from random activity. We illustrate the application of NMD to both simulated and real signals and demonstrate its qualitative and quantitative superiority over other approaches, such as (ensemble) empirical mode decomposition, Karhunen-Loève expansion, and independent component analysis. We point out that NMD is likely to be applicable and useful in many different areas of research, such as geophysics, finance, and the life sciences. The necessary matlab codes for running NMD are freely available for download.

  4. Nonlinear mode decomposition: a noise-robust, adaptive decomposition method.

    PubMed

    Iatsenko, Dmytro; McClintock, Peter V E; Stefanovska, Aneta

    2015-09-01

    The signals emanating from complex systems are usually composed of a mixture of different oscillations which, for a reliable analysis, should be separated from each other and from the inevitable background of noise. Here we introduce an adaptive decomposition tool-nonlinear mode decomposition (NMD)-which decomposes a given signal into a set of physically meaningful oscillations for any wave form, simultaneously removing the noise. NMD is based on the powerful combination of time-frequency analysis techniques-which, together with the adaptive choice of their parameters, make it extremely noise robust-and surrogate data tests used to identify interdependent oscillations and to distinguish deterministic from random activity. We illustrate the application of NMD to both simulated and real signals and demonstrate its qualitative and quantitative superiority over other approaches, such as (ensemble) empirical mode decomposition, Karhunen-Loève expansion, and independent component analysis. We point out that NMD is likely to be applicable and useful in many different areas of research, such as geophysics, finance, and the life sciences. The necessary matlab codes for running NMD are freely available for download. PMID:26465549

  5. Spectral decomposition of nonlinear systems with memory

    NASA Astrophysics Data System (ADS)

    Svenkeson, Adam; Glaz, Bryan; Stanton, Samuel; West, Bruce J.

    2016-02-01

    We present an alternative approach to the analysis of nonlinear systems with long-term memory that is based on the Koopman operator and a Lévy transformation in time. Memory effects are considered to be the result of interactions between a system and its surrounding environment. The analysis leads to the decomposition of a nonlinear system with memory into modes whose temporal behavior is anomalous and lacks a characteristic scale. On average, the time evolution of a mode follows a Mittag-Leffler function, and the system can be described using the fractional calculus. The general theory is demonstrated on the fractional linear harmonic oscillator and the fractional nonlinear logistic equation. When analyzing data from an ill-defined (black-box) system, the spectral decomposition in terms of Mittag-Leffler functions that we propose may uncover inherent memory effects through identification of a small set of dynamically relevant structures that would otherwise be obscured by conventional spectral methods. Consequently, the theoretical concepts we present may be useful for developing more general methods for numerical modeling that are able to determine whether observables of a dynamical system are better represented by memoryless operators, or operators with long-term memory in time, when model details are unknown.

  6. Regular Decompositions for H(div) Spaces

    SciTech Connect

    Kolev, Tzanio; Vassilevski, Panayot

    2012-01-01

    We study regular decompositions for H(div) spaces. In particular, we show that such regular decompositions are closely related to a previously studied “inf-sup” condition for parameter-dependent Stokes problems, for which we provide an alternative, more direct, proof.

  7. The Nonlinear Field Space Theory

    NASA Astrophysics Data System (ADS)

    Mielczarek, Jakub; Trześniewski, Tomasz

    2016-08-01

    In recent years the idea that not only the configuration space of particles, i.e. spacetime, but also the corresponding momentum space may have nontrivial geometry has attracted significant attention, especially in the context of quantum gravity. The aim of this letter is to extend this concept to the domain of field theories, by introducing field spaces (i.e. phase spaces of field values) that are not affine spaces. After discussing the motivation and general aspects of our approach we present a detailed analysis of the prototype (quantum) Nonlinear Field Space Theory of a scalar field on the Minkowski background. We show that the nonlinear structure of a field space leads to numerous interesting predictions, including: non-locality, generalization of the uncertainty relations, algebra deformations, constraining of the maximal occupation number, shifting of the vacuum energy and renormalization of the charge and speed of propagation of field excitations. Furthermore, a compact field space is a natural way to implement the "Principle of finiteness" of physical theories, which once motivated the Born-Infeld theory. Thus the presented framework has a variety of potential applications in the theories of fundamental interactions (e.g. quantum gravity), as well as in condensed matter physics (e.g. continuous spin chains), and can shed new light on the issue of divergences in quantum field theories.

  8. The Smooth Decomposition as a nonlinear modal analysis tool

    NASA Astrophysics Data System (ADS)

    Bellizzi, Sergio; Sampaio, Rubens

    2015-12-01

    The Smooth Decomposition (SD) is a statistical analysis technique for finding structures in an ensemble of spatially distributed data such that the vector directions not only keep the maximum possible variance but also the motions, along the vector directions, are as smooth in time as possible. In this paper, the notion of the dual smooth modes is introduced and used in the framework of oblique projection to expand a random response of a system. The dual modes define a tool that transforms the SD in an efficient modal analysis tool. The main properties of the SD are discussed and some new optimality properties of the expansion are deduced. The parameters of the SD give access to modal parameters of a linear system (mode shapes, resonance frequencies and modal energy participations). In case of nonlinear systems, a richer picture of the evolution of the modes versus energy can be obtained analyzing the responses under several excitation levels. This novel analysis of a nonlinear system is illustrated by an example.

  9. Nonlinear Decomposition of Climate Data: a New Method for Reconstruction of Dynamical Modes

    NASA Astrophysics Data System (ADS)

    Mukhin, D.; Gavrilov, A.; Loskutov, E. M.; Feigin, A. M.

    2014-12-01

    Modeling of multivariate time-series produced by complex systems requires efficient tools for compact data representation. In this report we consider this problem, relating to empirical modeling of climate, which implies an analysis of spatial-distributed time-series. The main goal is to establish the number of principal modes which have key contribution to data and actually governs the observed variability. Currently, the number of widely used methods based on PCA and factor analysis exists, which yield different data decompositions taking into consideration spatial/spatio-temporal correlations in observed dynamics: spatial empirical orthogonal functions, M-SSA, varimax rotation, empirical orthogonal teleconnections, and so on. However, the question about possibility of improving the decomposition by taking into account nonlinear couplings between variables often remains untouched. In the report the method for construction of principal dynamic modes on the basis of low-dimensional nonlinear parametric representation of observed multivariate time-series is suggested. It is aimed to extracting the set of latent modes that both explains an essential part of variability, and obeys the simplest evolution law. Thus, this approach can be used for optimal reconstruction of the phase space for empirical prognostic modeling of observed dynamics. Criterion of evidence of the nonlinearity, which allows estimating optimal parameters of data representation (including the number of parameters in mode definition and the number of principal modes) is proposed. The effectiveness of suggested method is firstly demonstrated on toy model example: two-dimensional strongly nonlinear problem of unknown dynamic mode retrieval from noisy data is considered. Next, the evidence of nonlinear couplings in SST space-distributed data covering the Globe is investigated by the proposed approach. It is demonstrated that the obtained principal modes capture more part of SST variability than

  10. Nonlinear Oscillators in Space Physics

    NASA Technical Reports Server (NTRS)

    Lester,Daniel; Thronson, Harley

    2011-01-01

    We discuss dynamical systems that produce an oscillation without an external time dependent source. Numerical results are presented for nonlinear oscillators in the Em1h's atmosphere, foremost the quasi-biennial oscillation (QBOl. These fluid dynamical oscillators, like the solar dynamo, have in common that one of the variables in a governing equation is strongly nonlinear and that the nonlinearity, to first order, has particular form. of 3rd or odd power. It is shown that this form of nonlinearity can produce the fundamental li'equency of the internal oscillation. which has a period that is favored by the dynamical condition of the fluid. The fundamental frequency maintains the oscillation, with no energy input to the system at that particular frequency. Nonlinearities of 2nd or even power could not maintain the oscillation.

  11. Nonlinear Elasticity in a Deforming Ambient Space

    NASA Astrophysics Data System (ADS)

    Yavari, Arash; Ozakin, Arkadas; Sadik, Souhayl

    2016-07-01

    In this paper, we formulate a nonlinear elasticity theory in which the ambient space is evolving. For a continuum moving in an evolving ambient space, we model time dependency of the metric by a time-dependent embedding of the ambient space in a larger manifold with a fixed background metric. We derive both the tangential and the normal governing equations. We then reduce the standard energy balance written in the larger ambient space to that in the evolving ambient space. We consider quasi-static deformations of the ambient space and show that a quasi-static deformation of the ambient space results in stresses, in general. We linearize the nonlinear theory about a reference motion and show that variation of the spatial metric corresponds to an effective field of body forces.

  12. Decomposition of Fuzzy Soft Sets with Finite Value Spaces

    PubMed Central

    Jun, Young Bae

    2014-01-01

    The notion of fuzzy soft sets is a hybrid soft computing model that integrates both gradualness and parameterization methods in harmony to deal with uncertainty. The decomposition of fuzzy soft sets is of great importance in both theory and practical applications with regard to decision making under uncertainty. This study aims to explore decomposition of fuzzy soft sets with finite value spaces. Scalar uni-product and int-product operations of fuzzy soft sets are introduced and some related properties are investigated. Using t-level soft sets, we define level equivalent relations and show that the quotient structure of the unit interval induced by level equivalent relations is isomorphic to the lattice consisting of all t-level soft sets of a given fuzzy soft set. We also introduce the concepts of crucial threshold values and complete threshold sets. Finally, some decomposition theorems for fuzzy soft sets with finite value spaces are established, illustrated by an example concerning the classification and rating of multimedia cell phones. The obtained results extend some classical decomposition theorems of fuzzy sets, since every fuzzy set can be viewed as a fuzzy soft set with a single parameter. PMID:24558342

  13. Decomposition of fuzzy soft sets with finite value spaces.

    PubMed

    Feng, Feng; Fujita, Hamido; Jun, Young Bae; Khan, Madad

    2014-01-01

    The notion of fuzzy soft sets is a hybrid soft computing model that integrates both gradualness and parameterization methods in harmony to deal with uncertainty. The decomposition of fuzzy soft sets is of great importance in both theory and practical applications with regard to decision making under uncertainty. This study aims to explore decomposition of fuzzy soft sets with finite value spaces. Scalar uni-product and int-product operations of fuzzy soft sets are introduced and some related properties are investigated. Using t-level soft sets, we define level equivalent relations and show that the quotient structure of the unit interval induced by level equivalent relations is isomorphic to the lattice consisting of all t-level soft sets of a given fuzzy soft set. We also introduce the concepts of crucial threshold values and complete threshold sets. Finally, some decomposition theorems for fuzzy soft sets with finite value spaces are established, illustrated by an example concerning the classification and rating of multimedia cell phones. The obtained results extend some classical decomposition theorems of fuzzy sets, since every fuzzy set can be viewed as a fuzzy soft set with a single parameter. PMID:24558342

  14. Decomposition of the complex system into nonlinear spatio-temporal modes: algorithm and application to climate data mining

    NASA Astrophysics Data System (ADS)

    Feigin, Alexander; Gavrilov, Andrey; Loskutov, Evgeny; Mukhin, Dmitry

    2015-04-01

    Proper decomposition of the complex system into well separated "modes" is a way to reveal and understand the mechanisms governing the system behaviour as well as discover essential feedbacks and nonlinearities. The decomposition is also natural procedure that provides to construct adequate and concurrently simplest models of both corresponding sub-systems, and of the system in whole. In recent works two new methods of decomposition of the Earth's climate system into well separated modes were discussed. The first method [1-3] is based on the MSSA (Multichannel Singular Spectral Analysis) [4] for linear expanding vector (space-distributed) time series and makes allowance delayed correlations of the processes recorded in spatially separated points. The second one [5-7] allows to construct nonlinear dynamic modes, but neglects delay of correlations. It was demonstrated [1-3] that first method provides effective separation of different time scales, but prevent from correct reduction of data dimension: slope of variance spectrum of spatio-temporal empirical orthogonal functions that are "structural material" for linear spatio-temporal modes, is too flat. The second method overcomes this problem: variance spectrum of nonlinear modes falls essentially sharply [5-7]. However neglecting time-lag correlations brings error of mode selection that is uncontrolled and increases with growth of mode time scale. In the report we combine these two methods in such a way that the developed algorithm allows constructing nonlinear spatio-temporal modes. The algorithm is applied for decomposition of (i) multi hundreds years globally distributed data generated by the INM RAS Coupled Climate Model [8], and (ii) 156 years time series of SST anomalies distributed over the globe [9]. We compare efficiency of different methods of decomposition and discuss the abilities of nonlinear spatio-temporal modes for construction of adequate and concurrently simplest ("optimal") models of climate systems

  15. Phase Space Structures Explain Hydrogen Atom Roaming in Formaldehyde Decomposition.

    PubMed

    Mauguière, Frédéric A L; Collins, Peter; Kramer, Zeb C; Carpenter, Barry K; Ezra, Gregory S; Farantos, Stavros C; Wiggins, Stephen

    2015-10-15

    We re-examine the prototypical roaming reaction--hydrogen atom roaming in formaldehyde decomposition--from a phase space perspective. Specifically, we address the question "why do trajectories roam, rather than dissociate through the radical channel?" We describe and compute the phase space structures that define and control all possible reactive events for this reaction, as well as provide a dynamically exact description of the roaming region in phase space. Using these phase space constructs, we show that in the roaming region, there is an unstable periodic orbit whose stable and unstable manifolds define a conduit that both encompasses all roaming trajectories exiting the formaldehyde well and shepherds them toward the H2···CO well. PMID:26499774

  16. Efficient variants of the vertex space domain decomposition algorithm

    SciTech Connect

    Chan, T.F.; Shao, J.P. . Dept. of Mathematics); Mathew, T.P. . Dept. of Mathematics)

    1994-11-01

    Several variants of the vertex space algorithm of Smith for two-dimensional elliptic problems are described. The vertex space algorithm is a domain decomposition method based on nonoverlapping subregions, in which the reduced Schur complement system on the interface is solved using a generalized block Jacobi-type preconditioner, with the blocks corresponding to the vertex space, edges, and a coarse grid. Two kinds of approximations are considered for the edge and vertex space subblocks, one based on Fourier approximation, and another based on an algebraic probing technique in which sparse approximations to these subblocks are computed. The motivation is to improve the efficiency of the algorithm without sacrificing the optimal convergence rate. Numerical and theoretical results on the performance of these algorithms, including variants of an algorithm of Bramble, Pasciak, and Schatz are presented.

  17. Parallel domain decomposition and the solution of nonlinear systems of equations

    SciTech Connect

    Gropp, W.D. ); Keyes, D.E. . Dept. of Mechanical Engineering)

    1990-01-01

    Many linear systems arise as subproblems in the solution of nonlinear equations, either as part of a simple fixed-point of a Newton's method iteration. This paper considers the use of domain decomposition techniques for the solution of these linear problems in the context of solving a multicomponent system of nonlinear equations on two types of parallel processors. One of the computations is drawn from fluid dynamics and includes locally refined grids. Such problems require great computational resources, and domain, decomposition seems to offer a way to efficiently solve these problems on computers with significant parallelism. The domain decomposition approach used is as in Gropp and Keyes, modified to achieve better parallelism and to reduce the computational work. 13 refs., 2 figs., 5 tabs.

  18. Control of nonlinear flexible space structures

    NASA Astrophysics Data System (ADS)

    Shi, Jianjun

    With the advances made in computer technology and efficiency of numerical algorithms over last decade, the MPC strategies have become quite popular among control community. However, application of MPC or GPC to flexible space structure control has not been explored adequately in the literature. The work presented in this thesis primarily focuses on application of GPC to control of nonlinear flexible space structures. This thesis is particularly devoted to the development of various approximate dynamic models, design and assessment of candidate controllers, and extensive numerical simulations for a realistic multibody flexible spacecraft, namely, Jupiter Icy Moons Orbiter (JIMO)---a Prometheus class of spacecraft proposed by NASA for deep space exploratory missions. A stable GPC algorithm is developed for Multi-Input-Multi-Output (MIMO) systems. An end-point weighting (penalty) is used in the GPC cost function to guarantee the nominal stability of the closed-loop system. A method is given to compute the desired end-point state from the desired output trajectory. The methodologies based on Fake Algebraic Riccati Equation (FARE) and constrained nonlinear optimization, are developed for synthesis of state weighting matrix. This makes this formulation more practical. A stable reconfigurable GPC architecture is presented and its effectiveness is demonstrated on both aircraft as well as spacecraft model. A representative in-orbit maneuver is used for assessing the performance of various control strategies using various design models. Different approximate dynamic models used for analysis include linear single body flexible structure, nonlinear single body flexible structure, and nonlinear multibody flexible structure. The control laws evaluated include traditional GPC, feedback linearization-based GPC (FLGPC), reconfigurable GPC, and nonlinear dissipative control. These various control schemes are evaluated for robust stability and robust performance in the presence of

  19. Nonlinear symmetries on spaces admitting Killing tensors

    NASA Astrophysics Data System (ADS)

    Visinescu, Mihai

    2010-04-01

    Nonlinear symmetries corresponding to Killing tensors are investigated. The intimate relation between Killing-Yano tensors and non-standard supersymmetries is pointed out. The gravitational anomalies are absent if the hidden symmetry is associated with a Killing-Yano tensor. In the case of the nonlinear symmetries the dynamical algebras of the Dirac-type operators is more involved and could be organized as infinite dimensional algebras or superalgebras. The general results are applied to some concrete spaces involved in theories of modern physics. As a first example it is considered the 4-dimensional Euclidean Taub-NUT space and its generalizations introduced by Iwai and Katayama. One presents the infinite dimensional superalgebra of Dirac type operators on Taub-NUT space that could be seen as a graded loop superalgebra of the Kac-Moody type. The axial anomaly, interpreted as the index of the Dirac operator, is computed for the generalized Taub-NUT metrics. Finally the existence of the conformal Killing-Yano tensors is investigated for some spaces with mixed Sasakian structures.

  20. A New Method for Nonlinear and Nonstationary Time Series Analysis: The Empirical Mode Decomposition Method

    NASA Technical Reports Server (NTRS)

    Huang, Norden E.; Zukor, Dorothy J. (Technical Monitor)

    2001-01-01

    A new method for analyzing nonlinear and nonstationary data has been developed. The key part of the method is the Empirical Mode Decomposition method with which any complicated data set can be decomposed into a finite and often small number of Intrinsic Mode Functions (IMF). An IMF is defined as any function having the same numbers of zero-crossing and extrema, and also having symmetric envelopes defined by the local maxima and minima respectively. The IMF also admits well-behaved Hilbert transform. This decomposition method is adaptive, and, therefore, highly efficient. Since the decomposition is based on the local characteristic time scale of the data, it is applicable to nonlinear and nonstationary processes. With the Hilbert transform, the Intrinsic Mode Functions yield instantaneous frequencies as functions of time that give sharp identifications of imbedded structures. The final presentation of the results is an energy-frequency-time distribution, designated as the Hilbert Spectrum. Classical nonlinear system models are used to illustrate the roles played by the nonlinear and nonstationary effects in the energy-frequency-time distribution.

  1. A Reconfigurable Sound Wave Decomposition Filterbank for Hearing Aids Based on Nonlinear Transformation.

    PubMed

    Huang, Shaoguang; Tian, Lan; Ma, Xiaojie; Wei, Ying

    2016-04-01

    Hearing impaired people have their own hearing loss characteristics and listening preferences. Therefore hearing aid system should become more natural, humanized and personalized, which requires the filterbank in hearing aids provides flexible sound wave decomposition schemes, so that patients are likely to use the most suitable scheme for their own hearing compensation strategy. In this paper, a reconfigurable sound wave decomposition filterbank is proposed. The prototype filter is first cosine modulated to generate uniform subbands. Then by non-linear transformation the uniform subbands are mapped to nonuniform subbands. By changing the control parameters, the nonlinear transformation changes which leads to different subbands allocations. It provides four different sound wave decomposition schemes without changing the structure of the filterbank. The performance of the proposed reconfigurable filterbank was compared with that of fixed filerbanks, fully customizable filterbanks and other existing reconfigurable filterbanks. It is shown that the proposed filterbank provides satisfactory matching performance as well as low complexity and delay, which make it suitable for real hearing aid applications. PMID:26168447

  2. Peculiar velocity decomposition, redshift space distortion, and velocity reconstruction in redshift surveys: The methodology

    NASA Astrophysics Data System (ADS)

    Zhang, Pengjie; Pan, Jun; Zheng, Yi

    2013-03-01

    Massive spectroscopic surveys will measure the redshift space distortion (RSD) induced by galaxy peculiar velocity to unprecedented accuracy and open a new era of precision RSD cosmology. We develop a new method to improve the RSD modeling and to carry out robust reconstruction of the 3D large scale peculiar velocity through galaxy redshift surveys, in light of RSD. (1) We propose a mathematically unique and physically motivated decomposition of peculiar velocity into three eigencomponents: an irrotational component completely correlated with the underlying density field (vδ), an irrotational component uncorrelated with the density field (vS), and a rotational (curl) component (vB). The three components have different origins, different scale dependences, and different impacts on RSD. (2) This decomposition has the potential to simplify and improve the RSD modeling. (i) vB damps the redshift space clustering. (ii) vS causes both damping and enhancement to the redshift space power spectrum Ps(k,u). Nevertheless, the leading order contribution to the enhancement has a u4 directional dependence, distinctively different from the Kaiser formula. Here, u≡kz/k, k is the amplitude of the wave vector, and kz is the component along the line of sight. (iii) vδ is of the greatest importance for the RSD cosmology. We find that the induced redshift clustering shows a number of important deviations from the usual Kaiser formula. Even in the limit of vS→0 and vB→0, the leading order contribution ∝(1+fW˜(k)u2)2. It differs from the Kaiser formula by a window function W˜(k). Nonlinear evolution generically drives W˜(k)≤1. We hence identify a significant systematical error causing underestimation of the structure growth parameter f by as much as O(10%) even at a relatively large scale k=0.1h/Mpc. (iv) The velocity decomposition reveals the three origins of the “finger-of-God” (FOG) effect and suggests how to simplify and improve the modeling of FOG by treating the

  3. Ship classification using nonlinear features of radiated sound: an approach based on empirical mode decomposition.

    PubMed

    Bao, Fei; Li, Chen; Wang, Xinlong; Wang, Qingfu; Du, Shuanping

    2010-07-01

    Classification for ship-radiated underwater sound is one of the most important and challenging subjects in underwater acoustical signal processing. An approach to ship classification is proposed in this work based on analysis of ship-radiated acoustical noise in subspaces of intrinsic mode functions attained via the ensemble empirical mode decomposition. It is shown that detection and acquisition of stable and reliable nonlinear features become practically feasible by nonlinear analysis of the time series of individual decomposed components, each of which is simple enough and well represents an oscillatory mode of ship dynamics. Surrogate and nonlinear predictability analysis are conducted to probe and measure the nonlinearity and regularity. The results of both methods, which verify each other, substantiate that ship-radiated noises contain components with deterministic nonlinear features well serving for efficient classification of ships. The approach perhaps opens an alternative avenue in the direction toward object classification and identification. It may also import a new view of signals as complex as ship-radiated sound. PMID:20649216

  4. Harmonic decomposition to describe the nonlinear evolution of stimulated Brillouin scattering

    SciTech Connect

    Hueller, S.; Masson-Laborde, P.E.; Pesme, D.; Casanova, M.; Detering, F.; Maximov, A.

    2006-02-15

    An efficient method to describe the nonlinear evolution of stimulated Brillouin scattering (SBS) in long scale-length plasmas is presented in the limit of a fluid description. The method is based on the decomposition of the various functions characterizing the plasma into their long- and short-wavelength components. It makes it possible to describe self-consistently the interplay between the plasma hydrodynamics, stimulated Brillouin scattering, and the generation of harmonics of the excited ion acoustic wave (IAW). This description is benchmarked numerically in one and two spatial dimensions [one dimensional (1D), two dimensional (2D)], by comparing the numerical results obtained along this method with those provided by a numerical code in which the decomposition into separate spatial scales is not made. The decomposition method proves to be very efficient in terms of computing time, especially in 2D, and very reliable, even in the extreme case of undamped ion acoustic waves. A novel picture of the SBS nonlinear behavior arises, in which the IAW harmonics generation gives rise to local defects appearing in the density and velocity hydrodynamics profiles. Consequently, SBS develops in various spatial domains which seem to be decorrelated one from each other, so that the backscattered Brillouin light is the sum of various backscattered waves generated in several independent spatial domains. It follows that the SBS reflectivity is chaotic in time and the resulting time-averaged value is significantly reduced as compared to the case when the IAW harmonics generation and flow modification are ignored. From the results of extensive numerical simulations carried out in 1D and 2D, we are able to infer the SBS reflectivity scaling law as a function of the plasma parameters and laser intensity, in the limit where the kinetic effects are negligible. It appears that this scaling law can be derived in the limit where the IAW harmonics generation is modeled simply by a

  5. Decomposition

    USGS Publications Warehouse

    Middleton, Beth A.

    2014-01-01

    A cornerstone of ecosystem ecology, decomposition was recognized as a fundamental process driving the exchange of energy in ecosystems by early ecologists such as Lindeman 1942 and Odum 1960). In the history of ecology, studies of decomposition were incorporated into the International Biological Program in the 1960s to compare the nature of organic matter breakdown in various ecosystem types. Such studies still have an important role in ecological studies of today. More recent refinements have brought debates on the relative role microbes, invertebrates and environment in the breakdown and release of carbon into the atmosphere, as well as how nutrient cycling, production and other ecosystem processes regulated by decomposition may shift with climate change. Therefore, this bibliography examines the primary literature related to organic matter breakdown, but it also explores topics in which decomposition plays a key supporting role including vegetation composition, latitudinal gradients, altered ecosystems, anthropogenic impacts, carbon storage, and climate change models. Knowledge of these topics is relevant to both the study of ecosystem ecology as well projections of future conditions for human societies.

  6. Nonlinear whistler wave scattering in space plasmas

    SciTech Connect

    Yukhimuk, V.; Roussel-Dupre, R.

    1997-04-01

    In this paper the evolution of nonlinear scattering of whistler mode waves by kinetic Alfven waves (KAW) in time and two spatial dimensions is studied analytically. The authors suggest this nonlinear process as a mechanism of kinetic Alfven wave generation in space plasmas. This mechanism can explain the dependence of Alfven wave generation on whistler waves observed in magnetospheric and ionospheric plasmas. The observational data show a dependence for the generation of long periodic pulsations Pc5 on whistler wave excitation in the auroral and subauroral zone of the magnetosphere. This dependence was first observed by Ondoh T.I. For 79 cases of VLF wave excitation registered by Ondoh at College Observatory (L=64.6 N), 52 of them were followed by Pc5 geomagnetic pulsation generation. Similar results were obtained at the Loparskaia Observatory (L=64 N) for auroral and subauroral zone of the magnetosphere. Thus, in 95% of the cases when VLF wave excitation occurred the generation of long periodic geomagnetic pulsations Pc5 were observed. The observations also show that geomagnetic pulsations Pc5 are excited simultaneously or insignificantly later than VLF waves. In fact these two phenomena are associated genetically: the excitation of VLF waves leads to the generation of geomagnetic pulsations Pc5. The observations show intensive generation of geomagnetic pulsations during thunderstorms. Using an electromagnetic noise monitoring system covering the ULF range (0.01-10 Hz) A.S. Fraser-Smith observed intensive ULF electromagnetic wave during a large thunderstorm near the San-Francisco Bay area on September 23, 1990. According to this data the most significant amplification in ULF wave activity was observed for waves with a frequency of 0.01 Hz and it is entirely possible that stronger enhancements would have been measured at lower frequencies.

  7. Decompositions of the polyhedral product functor with applications to moment-angle complexes and related spaces.

    PubMed

    Bahri, A; Bendersky, M; Cohen, F R; Gitler, S

    2009-07-28

    This article gives a natural decomposition of the suspension of a generalized moment-angle complex or partial product space which arises as the polyhedral product functor described below. The introduction and application of the smash product moment-angle complex provides a precise identification of the stable homotopy type of the values of the polyhedral product functor. One direct consequence is an analysis of the associated cohomology. For the special case of the complements of certain subspace arrangements, the geometrical decomposition implies the homological decomposition in earlier work of others as described below. Because the splitting is geometric, an analogous homological decomposition for a generalized moment-angle complex applies for any homology theory. Implied, therefore, is a decomposition for the Stanley-Reisner ring of a finite simplicial complex, and natural generalizations. PMID:19620727

  8. Reduced-order models for nonlinear vibrations of cylindrical shells via the proper orthogonal decomposition method

    NASA Astrophysics Data System (ADS)

    Amabili, M.; Sarkar, A.; Païdoussis, M. P.

    2003-09-01

    The nonlinear (large-amplitude) response of perfect and imperfect, simply supported circular cylindrical shells to harmonic excitation in the spectral neighbourhood of some of their lowest natural frequencies is investigated. The shell is assumed to be completely filled with an incompressible and inviscid fluid at rest. Donnell's nonlinear shallow-shell theory is used, and the solution is obtained by the Galerkin method. The proper orthogonal decomposition (POD) method is used to extract proper orthogonal modes that describe the system behaviour from time-series response data. These time series have been obtained via the conventional Galerkin approach (using normal modes as a projection basis) with an accurate model involving 16 degrees of freedom, validated in previous studies. The POD method, in conjunction with the Galerkin approach, permits a lower-dimensional model as compared to those obtainable via the conventional Galerkin approach. Different proper orthogonal modes computed from time series at different excitation frequencies are used and solutions are compared. Some of these sets of modes are capable of describing the system behaviour over the whole frequency range around the fundamental resonance with good accuracy and with only 3 degrees of freedom. They allow a drastic reduction in the computational effort, as compared to using the 16 degree-of-freedom model necessary when the conventional Galerkin approach is used.

  9. Nonlinear reduced order models for fluids systems using extended dynamic mode decomposition

    NASA Astrophysics Data System (ADS)

    Dawson, Scott; Rowley, Clarence

    2015-11-01

    The development of techniques that can extract simple, accurate, and computationally tractable models from fluids data is of importance for enhanced prediction, control, and fundamental understanding of such systems. Modeling approaches can take the form of identifying modes upon which to project the governing equations (e.g., Galerkin projection onto a set of POD modes), or in determining (or calibrating) the temporal dynamics from data, such as in dynamic mode decomposition (DMD), or various modifications to Galerkin projection. Here, we demonstrate that choosing appropriate observables (such as linear and quadratic monomials of POD coefficients) can allow for nonlinear behavior to be accurately captured using the recently proposed extended DMD algorithm. For cylinder wake data spanning the transient and vortex shedding limit cycle regimes, the identified nonlinear models show significant improvement in accuracy and robustness over standard DMD and Galerkin projection. Compared to traditional DMD, this approach should also allow for a better global approximation of the Koopman operator for the dynamical system. We make connections with other related model identification algorithms, and additionally investigate the performance of the method upon spatially sparse and noisy data. This work was supported by the Air Force Office of Scientific Research, under award No. FA9550-12-1-0075.

  10. Spinodal Decomposition for theCahn-Hilliard Equation in Higher Dimensions:Nonlinear Dynamics

    NASA Astrophysics Data System (ADS)

    Maier-Paape, Stanislaus; Wanner, Thomas

    This paper addresses the phenomenon of spinodal decomposition for the Cahn-Hilliard equation where Ω⊂n, n∈{1,2,3 }, is a bounded domain with sufficiently smooth boundary, and f is cubic-like, for example f(u) =u-u3. Based on the results of [26] the nonlinear Cahn-Hilliard equation will be discussed. This equation generates a nonlinear semiflow in certain affine subspaces of H2(Ω). In a neighborhood Uɛ with size proportional to ɛn around the constant solution , where μ lies in the spinodal region, we observe the following behavior. Within a local inertial manifold containing there exists a finite-dimensional invariant manifold which dominates the behavior of all solutions starting with initial conditions from a small ball around with probability almost 1. The dimension of is proportional to ɛ-n and the elements of exhibit a common geometric quantity which is strongly related to a characteristic wavelength proportional to ɛ.

  11. Detection of the Third Heart Sound Based on Nonlinear Signal Decomposition and Time-Frequency Localization.

    PubMed

    Barma, Shovan; Chen, Bo-Wei; Ji, Wen; Rho, Seungmin; Chou, Chih-Hung; Wang, Jhing-Fa

    2016-08-01

    This study presents a precise way to detect the third ( S3 ) heart sound, which is recognized as an important indication of heart failure, based on nonlinear single decomposition and time-frequency localization. The detection of the S3 is obscured due to its significantly low energy and frequency. Even more, the detected S3 may be misunderstood as an abnormal second heart sound with a fixed split, which was not addressed in the literature. To detect such S3, the Hilbert vibration decomposition method is applied to decompose the heart sound into a certain number of subcomponents while intactly preserving the phase information. Thus, the time information of all of the decomposed components are unchanged, which further expedites the identification and localization of any module/section of a signal properly. Next, the proposed localization step is applied to the decomposed subcomponents by using smoothed pseudo Wigner-Ville distribution followed by the reassignment method. Finally, based on the positional information, the S3 is distinguished and confirmed by measuring time delays between the S2 and S3. In total, 82 sets of cardiac cycles collected from different databases including Texas Heart Institute database are examined for evaluation of the proposed method. The result analysis shows that the proposed method can detect the S3 correctly, even when the normalized temporal energy of S3 is larger than 0.16, and the frequency of those is larger than 34 Hz. In a performance analysis, the proposed method demonstrates that the accuracy rate of S3 detection is as high as 93.9%, which is significantly higher compared with the other methods. Such findings prove the robustness of the proposed idea for detecting substantially low-energized S3 . PMID:26584485

  12. The geometry of the light-cone cell decomposition of moduli space

    SciTech Connect

    Garner, David Ramgoolam, Sanjaye

    2015-11-15

    The moduli space of Riemann surfaces with at least two punctures can be decomposed into a cell complex by using a particular family of ribbon graphs called Nakamura graphs. We distinguish the moduli space with all punctures labelled from that with a single labelled puncture. In both cases, we describe a cell decomposition where the cells are parametrised by graphs or equivalence classes of finite sequences (tuples) of permutations. Each cell is a convex polytope defined by a system of linear equations and inequalities relating light-cone string parameters, quotiented by the automorphism group of the graph. We give explicit examples of the cell decomposition at low genus with few punctures.

  13. Proper orthogonal decomposition method for analysis of nonlinear panel flutter with thermal effects in supersonic flow

    NASA Astrophysics Data System (ADS)

    Xie, Dan; Xu, Min; Dai, Honghua; Dowell, Earl H.

    2015-02-01

    The proper orthogonal decomposition (POD) method for analysis of nonlinear panel flutter subjected to supersonic flow is presented. Optimal POD modes are extracted from a chaotic Galerkin mode responses. The aeroelastic equations of motion are constructed using von Karman plate theory, first-order piston theory and quasi-steady thermal stress theory. A simply-supported plate with thermal loads from a uniformly distributed temperature is considered. Many types of panel behaviors, including stable flat, dynamically stable buckled, limit cycle oscillation, nonharmonic periodic motion, quasi-periodic motion and chaotic motion are observed. Our primary focus is on chaos and the route to chaos. It is found that a sudden transition from the buckled state to chaos occurs. Time history, phase portrait, Poincaré map, bifurcation diagram and Lyapunov exponent are employed to study chaos. The POD chaotic results obtained are compared with the traditional Galerkin solutions. It is shown that the POD method can obtain accurate chaotic solutions, using fewer modes and less computational effort than the Galerkin mode approach; additionally, the POD method converges faster in the analysis of chaotic transients. Effects of length-to-width ratios and thermal loads are presented. It is found that a smaller width for fixed length will produce more stable flutter response, while the thermal loads degrade the flutter boundary and result in a more complex evolution of dynamic motions. The numerical simulations show that the robustness of the POD modes depends on the dynamic pressure but not on temperature.

  14. Cardiopulmonary Resuscitation Pattern Evaluation Based on Ensemble Empirical Mode Decomposition Filter via Nonlinear Approaches.

    PubMed

    Sadrawi, Muammar; Sun, Wei-Zen; Ma, Matthew Huei-Ming; Dai, Chun-Yi; Abbod, Maysam F; Shieh, Jiann-Shing

    2016-01-01

    Good quality cardiopulmonary resuscitation (CPR) is the mainstay of treatment for managing patients with out-of-hospital cardiac arrest (OHCA). Assessment of the quality of the CPR delivered is now possible through the electrocardiography (ECG) signal that can be collected by an automated external defibrillator (AED). This study evaluates a nonlinear approximation of the CPR given to the asystole patients. The raw ECG signal is filtered using ensemble empirical mode decomposition (EEMD), and the CPR-related intrinsic mode functions (IMF) are chosen to be evaluated. In addition, sample entropy (SE), complexity index (CI), and detrended fluctuation algorithm (DFA) are collated and statistical analysis is performed using ANOVA. The primary outcome measure assessed is the patient survival rate after two hours. CPR pattern of 951 asystole patients was analyzed for quality of CPR delivered. There was no significant difference observed in the CPR-related IMFs peak-to-peak interval analysis for patients who are younger or older than 60 years of age, similarly to the amplitude difference evaluation for SE and DFA. However, there is a difference noted for the CI (p < 0.05). The results show that patients group younger than 60 years have higher survival rate with high complexity of the CPR-IMFs amplitude differences. PMID:27529068

  15. Cardiopulmonary Resuscitation Pattern Evaluation Based on Ensemble Empirical Mode Decomposition Filter via Nonlinear Approaches

    PubMed Central

    Ma, Matthew Huei-Ming

    2016-01-01

    Good quality cardiopulmonary resuscitation (CPR) is the mainstay of treatment for managing patients with out-of-hospital cardiac arrest (OHCA). Assessment of the quality of the CPR delivered is now possible through the electrocardiography (ECG) signal that can be collected by an automated external defibrillator (AED). This study evaluates a nonlinear approximation of the CPR given to the asystole patients. The raw ECG signal is filtered using ensemble empirical mode decomposition (EEMD), and the CPR-related intrinsic mode functions (IMF) are chosen to be evaluated. In addition, sample entropy (SE), complexity index (CI), and detrended fluctuation algorithm (DFA) are collated and statistical analysis is performed using ANOVA. The primary outcome measure assessed is the patient survival rate after two hours. CPR pattern of 951 asystole patients was analyzed for quality of CPR delivered. There was no significant difference observed in the CPR-related IMFs peak-to-peak interval analysis for patients who are younger or older than 60 years of age, similarly to the amplitude difference evaluation for SE and DFA. However, there is a difference noted for the CI (p < 0.05). The results show that patients group younger than 60 years have higher survival rate with high complexity of the CPR-IMFs amplitude differences. PMID:27529068

  16. Efficiently enclosing the compact binary parameter space by singular-value decomposition

    SciTech Connect

    Cannon, Kipp; Hanna, Chad; Keppel, Drew

    2011-10-15

    Gravitational-wave searches for the merger of compact binaries use matched filtering as the method of detecting signals and estimating parameters. Such searches construct a fine mesh of filters covering a signal parameter space at high density. Previously it has been shown that singular-value decomposition can reduce the effective number of filters required to search the data. Here we study how the basis provided by the singular-value decomposition changes dimension as a function of template-bank density. We will demonstrate that it is sufficient to use the basis provided by the singular-value decomposition of a low-density bank to accurately reconstruct arbitrary points within the boundaries of the template bank. Since this technique is purely numerical, it may have applications to interpolating the space of numerical relativity waveforms.

  17. A k-Space Method for Moderately Nonlinear Wave Propagation

    PubMed Central

    Jing, Yun; Wang, Tianren; Clement, Greg T.

    2013-01-01

    A k-space method for moderately nonlinear wave propagation in absorptive media is presented. The Westervelt equation is first transferred into k-space via Fourier transformation, and is solved by a modified wave-vector time-domain scheme. The present approach is not limited to forward propagation or parabolic approximation. One- and two-dimensional problems are investigated to verify the method by comparing results to analytic solutions and finite-difference time-domain (FDTD) method. It is found that to obtain accurate results in homogeneous media, the grid size can be as little as two points per wavelength, and for a moderately nonlinear problem, the Courant–Friedrichs–Lewy number can be as large as 0.4. Through comparisons with the conventional FDTD method, the k-space method for nonlinear wave propagation is shown here to be computationally more efficient and accurate. The k-space method is then employed to study three-dimensional nonlinear wave propagation through the skull, which shows that a relatively accurate focusing can be achieved in the brain at a high frequency by sending a low frequency from the transducer. Finally, implementations of the k-space method using a single graphics processing unit shows that it required about one-seventh the computation time of a single-core CPU calculation. PMID:22899114

  18. Performance of Scattering Matrix Decomposition and Color Spaces for Synthetic Aperture Radar Imagery

    NASA Astrophysics Data System (ADS)

    Terzuoli, Andrew; Arriagada, Manuel; Saville, Michael

    Polarimetrc Synthetic Aperture Radar (SAR) has been shown to be a powerful tool in re-mote sensing because uses up to four simultaneous measurements giving additional degrees of freedom for processing. Typically, polarization decomposition techniques are applied to the polarization-dependent data to form colorful imagery that is easy for operators systems to interpret. Yet, the presumption is that the SAR system operates with maximum bandwidth which requires extensive processing for near-or real-time application. In this research, color space selection is investigated when processing sparse polarimetric SAR data as in the case of the publicly available Volumetric SAR Data Set, Version 1:0". To improve information quality in resultant color imagery, three scattering matrix decompositions were investigated (linear, Pauli and Krogager) using two common color spaces (RGB, CMY) to deter-mine the best combination for accurate feature extraction. A mathematical model is presented for each de-composition technique and color space to the Cramer-Rao lower bound (CRLB) and quantify the performance bounds from an estimation perspective for given SAR system and processing parameters. After a deep literature review in color science, the mathematical model for color spaces was not able to be computed together with the mathematical model for decomposition techniques. The color spaces used for this research were functions of variables that are out of the scope of electrical engineering research and include factors such as the way humans sense color, envi-ronment inuences in the color stimulus and device technical characteristics used to display the SAR image. Hence, SAR imagery was computed for speci c combinations of decomposition technique and color space and allow the reader to gain an abstract view of the performance differences. The views expressed in this article are those of the authors and do not reflect the official policy of the U.S. Air Force, U.S. Department of Defense

  19. Nonlinear dynamic phenomena in the space shuttle thermal protection system

    NASA Technical Reports Server (NTRS)

    Housner, J. M.; Edighoffer, H. H.; Park, K. C.

    1981-01-01

    The development of an analysis for examining the nonlinear dynamic phenomena arising in the space shuttle orbiter tile/pad thermal protection system is presented. The tile/pad system consists of ceramic tiles bonded to the aluminum skin of the orbiter through a thin nylon felt pad. The pads are a soft nonlinear material which permits large strains and displays both hysteretic and nonlinear viscous damping. Application of the analysis to a square tile subjected to transverse sinusoidal motion of the orbiter skin is presented and the following nonlinear dynamic phenomena are considered: highly distorted wave forms, amplitude-dependent resonant frequencies which initially decrease and then increase with increasing amplitude of motion, magnification of substrate motion which is higher than would be expected in a similarly highly damped linear system, and classical parametric resonance instability.

  20. Nonlinear periodic space-charge waves in plasma

    SciTech Connect

    Kovalev, V. A.

    2009-05-15

    A solution is obtained in the form of coupled nonlinear periodic space-charge waves propagating in a magnetoactive plasma. The wave spectrum in the vicinity of the critical point, where the number of harmonics increases substantially, is found to fall with harmonic number as {proportional_to} s{sup -1/3}. Periodic space-charge waves are invoked to explain the zebra pattern in the radio emission from solar flares.

  1. Nonlinear Analysis of the Space Shuttle Superlightweight External Fuel Tank

    NASA Technical Reports Server (NTRS)

    Nemeth, Michael P.; Britt, Vicki O.; Collins, Timothy J.; Starnes, James H., Jr.

    1996-01-01

    Results of buckling and nonlinear analyses of the Space Shuttle external tank superlightweight liquid-oxygen (LO2) tank are presented. Modeling details and results are presented for two prelaunch loading conditions and for two full-scale structural tests that were conducted on the original external tank. The results illustrate three distinctly different types of nonlinear response for thin-walled shells subjected to combined mechanical and thermal loads. The nonlinear response phenomena consist of bifurcation-type buckling, short-wavelength nonlinear bending, and nonlinear collapse associated with a limit point. For each case, the results show that accurate predictions of non- linear behavior generally require a large-scale, high-fidelity finite-element model. Results are also presented that show that a fluid-filled launch-vehicle shell can be highly sensitive to initial geometric imperfections. In addition, results presented for two full-scale structural tests of the original standard-weight external tank suggest that the finite-element modeling approach used in the present study is sufficient for representing the nonlinear behavior of the superlightweight LO2 tank.

  2. A novel technique to solve nonlinear higher-index Hessenberg differential-algebraic equations by Adomian decomposition method.

    PubMed

    Benhammouda, Brahim

    2016-01-01

    Since 1980, the Adomian decomposition method (ADM) has been extensively used as a simple powerful tool that applies directly to solve different kinds of nonlinear equations including functional, differential, integro-differential and algebraic equations. However, for differential-algebraic equations (DAEs) the ADM is applied only in four earlier works. There, the DAEs are first pre-processed by some transformations like index reductions before applying the ADM. The drawback of such transformations is that they can involve complex algorithms, can be computationally expensive and may lead to non-physical solutions. The purpose of this paper is to propose a novel technique that applies the ADM directly to solve a class of nonlinear higher-index Hessenberg DAEs systems efficiently. The main advantage of this technique is that; firstly it avoids complex transformations like index reductions and leads to a simple general algorithm. Secondly, it reduces the computational work by solving only linear algebraic systems with a constant coefficient matrix at each iteration, except for the first iteration where the algebraic system is nonlinear (if the DAE is nonlinear with respect to the algebraic variable). To demonstrate the effectiveness of the proposed technique, we apply it to a nonlinear index-three Hessenberg DAEs system with nonlinear algebraic constraints. This technique is straightforward and can be programmed in Maple or Mathematica to simulate real application problems. PMID:27330880

  3. Domain decomposition methods for nonconforming finite element spaces of Lagrange-type

    NASA Technical Reports Server (NTRS)

    Cowsar, Lawrence C.

    1993-01-01

    In this article, we consider the application of three popular domain decomposition methods to Lagrange-type nonconforming finite element discretizations of scalar, self-adjoint, second order elliptic equations. The additive Schwarz method of Dryja and Widlund, the vertex space method of Smith, and the balancing method of Mandel applied to nonconforming elements are shown to converge at a rate no worse than their applications to the standard conforming piecewise linear Galerkin discretization. Essentially, the theory for the nonconforming elements is inherited from the existing theory for the conforming elements with only modest modification by constructing an isomorphism between the nonconforming finite element space and a space of continuous piecewise linear functions.

  4. Nonlinear shell analyses of the space shuttle solid rocket boosters

    NASA Technical Reports Server (NTRS)

    Knight, Norman F., Jr.; Gillian, Ronnie E.; Nemeth, Michael P.

    1989-01-01

    A variety of structural analyses have been performed on the Solid Rocket Boosters (SRB's) to provide information that would contribute to the understanding of the failure which destroyed the Space Shuttle Challenger. This paper describes nonlinear shell analyses that were performed to characterize the behavior of an overall SRB structure and a segment of the SRB in the vicinity of the External Tank Attachment (ETA) ring. Shell finite element models were used that would accurately reflect the global load transfer in an SRB in a manner such that nonlinear shell collapse and ovalization could be assessed. The purpose of these analyses was to calculate the overall deflection and stress distributions for these SRB models when subjected to mechanical loads corresponding to critical times during the launch sequence. Static analyses of these SRB models were performed using a snapshot picture of the loads. Analytical results obtained using these models show no evidence of nonlinear shell collapse for the pre-liftoff loading cases considered.

  5. Nonlinear shell analysis of the Space Shuttle Solid Rocket Boosters

    NASA Technical Reports Server (NTRS)

    Knight, N. F., Jr.; Gillian, R. E.; Nemeth, M. P.

    1990-01-01

    A variety of structural analyses have been performed on the Solid Rocket Boosters (SRB's) to provide information that would contribute to the understanding of the failure which destroyed the Space Shuttle Challenger. This paper describes nonlinear shell analyses that were performed to characterize the behavior of an overall SRB structure and a segment of the SRB in the vicinity of the External Tank Attachment (ETA) ring. Shell finite element models were used that would accurately reflect the global load transfer in an SRB in a manner such that nonlinear shell collapse and ovalization could be assessed. The purpose of these analyses was to calculate the overall deflection and stress distributions for these SRB models when subjected to mechanical loads corresponding to critical times during the launch sequence. Static analyses of these SRB models were performed using a 'snapshot picture' of the loads. Analytical results obtained using these models show no evidence of nonlinear shell collapse for the pre-liftoff loading cases considered.

  6. Geometrically Nonlinear Finite Element Analysis of a Composite Space Reflector

    NASA Technical Reports Server (NTRS)

    Lee, Kee-Joo; Leet, Sung W.; Clark, Greg; Broduer, Steve (Technical Monitor)

    2001-01-01

    Lightweight aerospace structures, such as low areal density composite space reflectors, are highly flexible and may undergo large deflection under applied loading, especially during the launch phase. Accordingly, geometrically nonlinear analysis that takes into account the effect of finite rotation may be needed to determine the deformed shape for a clearance check and the stress and strain state to ensure structural integrity. In this study, deformation of the space reflector is determined under static conditions using a geometrically nonlinear solid shell finite element model. For the solid shell element formulation, the kinematics of deformation is described by six variables that are purely vector components. Because rotational angles are not used, this approach is free of the limitations of small angle increments. This also allows easy connections between substructures and large load increments with respect to the conventional shell formulation using rotational parameters. Geometrically nonlinear analyses were carried out for three cases of static point loads applied at selected points. A chart shows results for a case when the load is applied at the center point of the reflector dish. The computed results capture the nonlinear behavior of the composite reflector as the applied load increases. Also, they are in good agreement with the data obtained by experiments.

  7. Phase space lattices and integrable nonlinear wave equations

    NASA Astrophysics Data System (ADS)

    Tracy, Eugene; Zobin, Nahum

    2003-10-01

    Nonlinear wave equations in fluids and plasmas that are integrable by Inverse Scattering Theory (IST), such as the Korteweg-deVries and nonlinear Schrodinger equations, are known to be infinite-dimensional Hamiltonian systems [1]. These are of interest physically because they predict new phenomena not present in linear wave theories, such as solitons and rogue waves. The IST method provides solutions of these equations in terms of a special class of functions called Riemann theta functions. The usual approach to the theory of theta functions tends to obscure the underlying phase space structure. A theory due to Mumford and Igusa [2], however shows that the theta functions arise naturally in the study of phase space lattices. We will describe this theory, as well as potential applications to nonlinear signal processing and the statistical theory of nonlinear waves. 1] , S. Novikov, S. V. Manakov, L. P. Pitaevskii and V. E. Zakharov, Theory of solitons: the inverse scattering method (Consultants Bureau, New York, 1984). 2] D. Mumford, Tata lectures on theta, Vols. I-III (Birkhauser); J. Igusa, Theta functions (Springer-Verlag, New York, 1972).

  8. An investigation of the use of temporal decomposition in space mission scheduling

    NASA Technical Reports Server (NTRS)

    Bullington, Stanley E.; Narayanan, Venkat

    1994-01-01

    This research involves an examination of techniques for solving scheduling problems in long-duration space missions. The mission timeline is broken up into several time segments, which are then scheduled incrementally. Three methods are presented for identifying the activities that are to be attempted within these segments. The first method is a mathematical model, which is presented primarily to illustrate the structure of the temporal decomposition problem. Since the mathematical model is bound to be computationally prohibitive for realistic problems, two heuristic assignment procedures are also presented. The first heuristic method is based on dispatching rules for activity selection, and the second heuristic assigns performances of a model evenly over timeline segments. These heuristics are tested using a sample Space Station mission and a Spacelab mission. The results are compared with those obtained by scheduling the missions without any problem decomposition. The applicability of this approach to large-scale mission scheduling problems is also discussed.

  9. Space Vehicle Pose Estimation via Optical Correlation and Nonlinear Estimation

    NASA Technical Reports Server (NTRS)

    Rakoczy, John; Herren, Kenneth

    2007-01-01

    A technique for 6-degree-of-freedom (6DOF) pose estimation of space vehicles is being developed. This technique draws upon recent developments in implementing optical correlation measurements in a nonlinear estimator, which relates the optical correlation measurements to the pose states (orientation and position). For the optical correlator, the use of both conjugate filters and binary, phase-only filters in the design of synthetic discriminant function (SDF) filters is explored. A static neural network is trained a priori and used as the nonlinear estimator. New commercial animation and image rendering software is exploited to design the SDF filters and to generate a large filter set with which to train the neural network. The technique is applied to pose estimation for rendezvous and docking of free-flying spacecraft and to terrestrial surface mobility systems for NASA's Vision for Space Exploration. Quantitative pose estimation performance will be reported. Advantages and disadvantages of the implementation of this technique are discussed.

  10. Space Vehicle Pose Estimation via Optical Correlation and Nonlinear Estimation

    NASA Technical Reports Server (NTRS)

    Rakoczy, John M.; Herren, Kenneth A.

    2008-01-01

    A technique for 6-degree-of-freedom (6DOF) pose estimation of space vehicles is being developed. This technique draws upon recent developments in implementing optical correlation measurements in a nonlinear estimator, which relates the optical correlation measurements to the pose states (orientation and position). For the optical correlator, the use of both conjugate filters and binary, phase-only filters in the design of synthetic discriminant function (SDF) filters is explored. A static neural network is trained a priori and used as the nonlinear estimator. New commercial animation and image rendering software is exploited to design the SDF filters and to generate a large filter set with which to train the neural network. The technique is applied to pose estimation for rendezvous and docking of free-flying spacecraft and to terrestrial surface mobility systems for NASA's Vision for Space Exploration. Quantitative pose estimation performance will be reported. Advantages and disadvantages of the implementation of this technique are discussed.

  11. Hubble space telescope: Pointing error effects on nonlinear ball joints

    NASA Technical Reports Server (NTRS)

    Farmer, J. E.; Grissett, F. R.

    1985-01-01

    The Hubble Space Telescope pointing error produced by optical benches mounted on free ball joints is examined. Spacecraft cable connections are assumed to produce translational and rotational damping and restoring forces which act through the optical bench center of mass. The nonlinear dynamics are modeled and then implemented using an existing computer program for simulating the vehicle dynamics and pointing control system algorithm. Results are presented for the test case which indicate acceptable performance.

  12. Multiple color-image authentication system using HSI color space and QR decomposition in gyrator domains

    NASA Astrophysics Data System (ADS)

    Rafiq Abuturab, Muhammad

    2016-06-01

    A new multiple color-image authentication system based on HSI (Hue-Saturation-Intensity) color space and QR decomposition in gyrator domains is proposed. In this scheme, original color images are converted from RGB (Red-Green-Blue) color spaces to HSI color spaces, divided into their H, S, and I components, and then obtained corresponding phase-encoded components. All the phase-encoded H, S, and I components are individually multiplied, and then modulated by random phase functions. The modulated H, S, and I components are convoluted into a single gray image with asymmetric cryptosystem. The resulting image is segregated into Q and R parts by QR decomposition. Finally, they are independently gyrator transformed to get their encoded parts. The encoded Q and R parts should be gathered without missing anyone for decryption. The angles of gyrator transform afford sensitive keys. The protocol based on QR decomposition of encoded matrix and getting back decoded matrix after multiplying matrices Q and R, enhances the security level. The random phase keys, individual phase keys, and asymmetric phase keys provide high robustness to the cryptosystem. Numerical simulation results demonstrate that this scheme is the superior than the existing techniques.

  13. Some Remarks on Space-Time Decompositions, and Degenerate Metrics, in General Relativity

    NASA Astrophysics Data System (ADS)

    Bengtsson, Ingemar

    Space-time decomposition of the Hilbert-Palatini action, written in a form which admits degenerate metrics, is considered. Simple numerology shows why D = 3 and 4 are singled out as admitting a simple phase space. The canonical structure of the degenerate sector turns out to be awkward. However, the real degenerate metrics obtained as solutions are the same as those that occur in Ashtekar's formulation of complex general relativity. An exact solution of Ashtekar's equations, with degenerate metric, shows that the manifestly four-dimensional form of the action, and its 3 + 1 form, are not quite equivalent.

  14. On the polar decomposition of right linear operators in quaternionic Hilbert spaces

    NASA Astrophysics Data System (ADS)

    G, Ramesh; P, Santhosh Kumar

    2016-04-01

    In this article, we prove the existence of the polar decomposition of densely defined closed right linear operators in quaternionic Hilbert spaces: If T is a densely defined closed right linear operator in a quaternionic Hilbert space H, then there exists a partial isometry U0 such that T = U 0 |" separators=" T | . In fact U0 is unique if N(U0) = N(T). In particular, if H is separable and U is a partial isometry with T = U |" separators=" T | , then we prove that U = U0 if and only if either N(T) = {0} or R(T)⊥ = {0}.

  15. Communication: Active space decomposition with multiple sites: Density matrix renormalization group algorithm

    SciTech Connect

    Parker, Shane M.; Shiozaki, Toru

    2014-12-07

    We extend the active space decomposition method, recently developed by us, to more than two active sites using the density matrix renormalization group algorithm. The fragment wave functions are described by complete or restricted active-space wave functions. Numerical results are shown on a benzene pentamer and a perylene diimide trimer. It is found that the truncation errors in our method decrease almost exponentially with respect to the number of renormalization states M, allowing for numerically exact calculations (to a few μE{sub h} or less) with M = 128 in both cases. This rapid convergence is because the renormalization steps are used only for the interfragment electron correlation.

  16. Parameter spaces for linear and nonlinear whistler-mode waves

    SciTech Connect

    Summers, Danny; Tang, Rongxin; Omura, Yoshiharu; Lee, Dong-Hun

    2013-07-15

    We examine the growth of magnetospheric whistler-mode waves which comprises a linear growth phase followed by a nonlinear growth phase. We construct time-profiles for the wave amplitude that smoothly match at the transition between linear and nonlinear wave growth. This matching procedure can only take place over a limited “matching region” in (N{sub h}/N{sub 0},A{sub T})-space, where A{sub T} is the electron thermal anisotropy, N{sub h} is the hot (energetic) electron number density, and N{sub 0} is the cold (background) electron number density. We construct this matching region and determine how the matching wave amplitude varies throughout the region. Further, we specify a boundary in (N{sub h}/N{sub 0},A{sub T})-space that separates a region where only linear chorus wave growth can occur from the region in which fully nonlinear chorus growth is possible. We expect that this boundary should prove of practical use in performing computationally expensive full-scale particle simulations, and in interpreting experimental wave data.

  17. Modeling diffusion control on organic matter decomposition in unsaturated soil pore space

    NASA Astrophysics Data System (ADS)

    Vogel, Laure; Pot, Valérie; Garnier, Patricia; Vieublé-Gonod, Laure; Nunan, Naoise; Raynaud, Xavier; Chenu, Claire

    2014-05-01

    Soil Organic Matter decomposition is affected by soil structure and water content, but field and laboratory studies about this issue conclude to highly variable outcomes. Variability could be explained by the discrepancy between the scale at which key processes occur and the measurements scale. We think that physical and biological interactions driving carbon transformation dynamics can be best understood at the pore scale. Because of the spatial disconnection between carbon sources and decomposers, the latter rely on nutrient transport unless they can actively move. In hydrostatic case, diffusion in soil pore space is thus thought to regulate biological activity. In unsaturated conditions, the heterogeneous distribution of water modifies diffusion pathways and rates, thus affects diffusion control on decomposition. Innovative imaging and modeling tools offer new means to address these effects. We have developed a new model based on the association between a 3D Lattice-Boltzmann Model and an adimensional decomposition module. We designed scenarios to study the impact of physical (geometry, saturation, decomposers position) and biological properties on decomposition. The model was applied on porous media with various morphologies. We selected three cubic images of 100 voxels side from µCT-scanned images of an undisturbed soil sample at 68µm resolution. We used LBM to perform phase separation and obtained water phase distributions at equilibrium for different saturation indices. We then simulated the diffusion of a simple soluble substrate (glucose) and its consumption by bacteria. The same mass of glucose was added as a pulse at the beginning of all simulations. Bacteria were placed in few voxels either regularly spaced or concentrated close to or far from the glucose source. We modulated physiological features of decomposers in order to weight them against abiotic conditions. We could evidence several effects creating unequal substrate access conditions for

  18. On the nonlinear dynamics and control of large space structures

    SciTech Connect

    Modi, V.J.

    1994-12-31

    The paper reviews, using the Lagrangian approach, dynamics of flexible multibody systems, of contemporary interest, and their control. To begin with, a relatively general formulation for studying the dynamics and control of an arbitrary spacecraft with interconnected flexible bodies is developed accounting for thermal deflection, transient system properties, shift in the center of mass, shear deformations, rotary inertias and geometric nonlinearities. A rather self-contained, comprehensive, numerical algorithm using system as well as component modes follows which is applicable to a large class of spacecraft configurations. In particular, versatility of the approach is demonstrated through the dynamics and control studies aimed at the proposed Space Station and the associated mobile manipulator system.

  19. Longitudinal emittance growth due to nonlinear space charge effect

    NASA Astrophysics Data System (ADS)

    Lau, Y. Y.; Yu, Simon S.; Barnard, John J.; Seidl, Peter A.

    2012-03-01

    Emittance posts limits on the key requirements of final pulse length and spot size on target in heavy ion fusion drivers. In this paper, we show studies on the effect of nonlinear space charge on longitudinal emittance growth in the drift compression section. We perform simulations, using the 3D PIC code WARP, for a high current beam under conditions of bends and longitudinal compression. The linear growth rate for longitudinal emittance turns out to depend only on the peak line charge density, and is independent of pulse length, velocity tilt, and/or the pipe and beam size. This surprisingly simple result is confirmed by simulations and analytic calculations.

  20. Nonlinear bulk viscosity in FRW cosmology: a phase space analysis.

    NASA Astrophysics Data System (ADS)

    Acquaviva, G.; Beesham, A.

    2015-11-01

    We consider a Friedmann-Robertson-Walker spacetime filled with both viscous radiation and nonviscous dust. The former has a bulk viscosity that is proportional to an arbitrary power of the energy density, i.e. \\zeta \\propto {ρ }{{v}}ν , and viscous pressure satisfying a nonlinear evolution equation. The analysis is carried out in the context of dynamical systems and the properties of solutions corresponding to the fixed points are discussed. For some ranges of the relevant parameter ν we find that the trajectories in the phase space evolve from a FRW singularity towards an asymptotic de Sitter attractor, confirming and extending previous analysis in the literature.

  1. Nonlinear rotordynamics analysis. [Space Shuttle Main Engine turbopumps

    NASA Technical Reports Server (NTRS)

    Noah, Sherif T.

    1991-01-01

    Effective analysis tools were developed for predicting the nonlinear rotordynamic behavior of the Space Shuttle Main Engine (SSME) turbopumps under steady and transient operating conditions. Using these methods, preliminary parametric studies were conducted on both generic and actual HPOTP (high pressure oxygen turbopump) models. In particular, a novel modified harmonic balance/alternating Fourier transform (HB/AFT) method was developed and used to conduct a preliminary study of the effects of fluid, bearing and seal forces on the unbalanced response of a multi-disk rotor in the presence of bearing clearances. The method makes it possible to determine periodic, sub-, super-synchronous and chaotic responses of a rotor system. The method also yields information about the stability of the obtained response, thus allowing bifurcation analyses. This provides a more effective capability for predicting the response under transient conditions by searching in proximity of resonance peaks. Preliminary results were also obtained for the nonlinear transient response of an actual HPOTP model using an efficient, newly developed numerical method based on convolution integration. Currently, the HB/AFT is being extended for determining the aperiodic response of nonlinear systems. Initial results show the method to be promising.

  2. Space-by-time decomposition for single-trial decoding of M/EEG activity.

    PubMed

    Delis, Ioannis; Onken, Arno; Schyns, Philippe G; Panzeri, Stefano; Philiastides, Marios G

    2016-06-01

    We develop a novel methodology for the single-trial analysis of multichannel time-varying neuroimaging signals. We introduce the space-by-time M/EEG decomposition, based on Non-negative Matrix Factorization (NMF), which describes single-trial M/EEG signals using a set of non-negative spatial and temporal components that are linearly combined with signed scalar activation coefficients. We illustrate the effectiveness of the proposed approach on an EEG dataset recorded during the performance of a visual categorization task. Our method extracts three temporal and two spatial functional components achieving a compact yet full representation of the underlying structure, which validates and summarizes succinctly results from previous studies. Furthermore, we introduce a decoding analysis that allows determining the distinct functional role of each component and relating them to experimental conditions and task parameters. In particular, we demonstrate that the presented stimulus and the task difficulty of each trial can be reliably decoded using specific combinations of components from the identified space-by-time representation. When comparing with a sliding-window linear discriminant algorithm, we show that our approach yields more robust decoding performance across participants. Overall, our findings suggest that the proposed space-by-time decomposition is a meaningful low-dimensional representation that carries the relevant information of single-trial M/EEG signals. PMID:27033682

  3. Spatiotemporal Domain Decomposition for Massive Parallel Computation of Space-Time Kernel Density

    NASA Astrophysics Data System (ADS)

    Hohl, A.; Delmelle, E. M.; Tang, W.

    2015-07-01

    Accelerated processing capabilities are deemed critical when conducting analysis on spatiotemporal datasets of increasing size, diversity and availability. High-performance parallel computing offers the capacity to solve computationally demanding problems in a limited timeframe, but likewise poses the challenge of preventing processing inefficiency due to workload imbalance between computing resources. Therefore, when designing new algorithms capable of implementing parallel strategies, careful spatiotemporal domain decomposition is necessary to account for heterogeneity in the data. In this study, we perform octtree-based adaptive decomposition of the spatiotemporal domain for parallel computation of space-time kernel density. In order to avoid edge effects near subdomain boundaries, we establish spatiotemporal buffers to include adjacent data-points that are within the spatial and temporal kernel bandwidths. Then, we quantify computational intensity of each subdomain to balance workloads among processors. We illustrate the benefits of our methodology using a space-time epidemiological dataset of Dengue fever, an infectious vector-borne disease that poses a severe threat to communities in tropical climates. Our parallel implementation of kernel density reaches substantial speedup compared to sequential processing, and achieves high levels of workload balance among processors due to great accuracy in quantifying computational intensity. Our approach is portable of other space-time analytical tests.

  4. The effect of space holder content and decomposition methods in fabrication of aluminum foams by powder metallurgy method using carbamide space holder

    NASA Astrophysics Data System (ADS)

    Amirah, A. H.; Nurulakmal, M. S.; Anasyida, A. S.

    2016-07-01

    The effect of space holder amount and decomposition methods on the morphology, density and porosity and compressive properties of aluminum foams were investigated. Aluminum foam was fabricated by powder metallurgy method using spherical carbamide as space holder using three different decomposition method of carbamide includes; dissolution process, normal sintering process and two step sintering process. Aluminum foam with 60 wt.% carbamide has the lowest density and exhibited the highest porosity for all the decomposition. The results indicated that Al foams produced by dissolution method have the highest compressive properties with acceptable density and porosity value.

  5. Non-linear stochastic growth rates and redshift space distortions

    SciTech Connect

    Jennings, Elise; Jennings, David

    2015-04-09

    The linear growth rate is commonly defined through a simple deterministic relation between the velocity divergence and the matter overdensity in the linear regime. We introduce a formalism that extends this to a non-linear, stochastic relation between θ = ∇ ∙ v(x,t)/aH and δ. This provides a new phenomenological approach that examines the conditional mean <θ|δ>, together with the fluctuations of θ around this mean. We also measure these stochastic components using N-body simulations and find they are non-negative and increase with decreasing scale from ~10 per cent at k < 0.2 h Mpc-1 to 25 per cent at k ~ 0.45 h Mpc-1 at z = 0. Both the stochastic relation and non-linearity are more pronounced for haloes, M ≤ 5 × 1012 M h-1, compared to the dark matter at z = 0 and 1. Non-linear growth effects manifest themselves as a rotation of the mean <θ|δ> away from the linear theory prediction -fLTδ, where fLT is the linear growth rate. This rotation increases with wavenumber, k, and we show that it can be well-described by second-order Lagrangian perturbation theory (2LPT) fork < 0.1 h Mpc-1. Furthermore, the stochasticity in the θ – δ relation is not so simply described by 2LPT, and we discuss its impact on measurements of fLT from two-point statistics in redshift space. Furthermore, given that the relationship between δ and θ is stochastic and non-linear, this will have implications for the interpretation and precision of fLT extracted using models which assume a linear, deterministic expression.

  6. Non-linear stochastic growth rates and redshift space distortions

    NASA Astrophysics Data System (ADS)

    Jennings, Elise; Jennings, David

    2015-06-01

    The linear growth rate is commonly defined through a simple deterministic relation between the velocity divergence and the matter overdensity in the linear regime. We introduce a formalism that extends this to a non-linear, stochastic relation between θ = nabla \\cdot v({x},t)/aH and δ. This provides a new phenomenological approach that examines the conditional mean <θ|δ>, together with the fluctuations of θ around this mean. We measure these stochastic components using N-body simulations and find they are non-negative and increase with decreasing scale from ˜10 per cent at k < 0.2 h Mpc-1 to 25 per cent at k ˜ 0.45 h Mpc-1 at z = 0. Both the stochastic relation and non-linearity are more pronounced for haloes, M ≤ 5 × 1012 M⊙ h-1, compared to the dark matter at z = 0 and 1. Non-linear growth effects manifest themselves as a rotation of the mean <θ|δ> away from the linear theory prediction -fLTδ, where fLT is the linear growth rate. This rotation increases with wavenumber, k, and we show that it can be well-described by second-order Lagrangian perturbation theory (2LPT) for k < 0.1 h Mpc-1. The stochasticity in the θ-δ relation is not so simply described by 2LPT, and we discuss its impact on measurements of fLT from two-point statistics in redshift space. Given that the relationship between δ and θ is stochastic and non-linear, this will have implications for the interpretation and precision of fLT extracted using models which assume a linear, deterministic expression.

  7. Non-linear stochastic growth rates and redshift space distortions

    DOE PAGESBeta

    Jennings, Elise; Jennings, David

    2015-04-09

    The linear growth rate is commonly defined through a simple deterministic relation between the velocity divergence and the matter overdensity in the linear regime. We introduce a formalism that extends this to a non-linear, stochastic relation between θ = ∇ ∙ v(x,t)/aH and δ. This provides a new phenomenological approach that examines the conditional mean <θ|δ>, together with the fluctuations of θ around this mean. We also measure these stochastic components using N-body simulations and find they are non-negative and increase with decreasing scale from ~10 per cent at k < 0.2 h Mpc-1 to 25 per cent at kmore » ~ 0.45 h Mpc-1 at z = 0. Both the stochastic relation and non-linearity are more pronounced for haloes, M ≤ 5 × 1012 M⊙ h-1, compared to the dark matter at z = 0 and 1. Non-linear growth effects manifest themselves as a rotation of the mean <θ|δ> away from the linear theory prediction -fLTδ, where fLT is the linear growth rate. This rotation increases with wavenumber, k, and we show that it can be well-described by second-order Lagrangian perturbation theory (2LPT) fork < 0.1 h Mpc-1. Furthermore, the stochasticity in the θ – δ relation is not so simply described by 2LPT, and we discuss its impact on measurements of fLT from two-point statistics in redshift space. Furthermore, given that the relationship between δ and θ is stochastic and non-linear, this will have implications for the interpretation and precision of fLT extracted using models which assume a linear, deterministic expression.« less

  8. Nonlinear instabilities driven by coherent phase-space structures

    NASA Astrophysics Data System (ADS)

    Lesur, Maxime

    2012-10-01

    Coherent phase-space (PS) structures are an important feature of plasma turbulence. They can drive nonlinear instabilities [1], intermittency in drift-wave turbulence [2], and transport [3]. We aim at a comprehensive understanding of turbulence, not just as an ensemble of waves, as quasilinear theory implies, but as a mixture of coupled waves and localized structures. This work, which focuses on isolated PS structures, is a fundamental advance in this direction. We analyze the effects of self-binding negative fluctuations (PS holes) on stability, intermittency and anomalous resistivity, both analytically and numerically. We present a new theory which describes the growth of a hole or clump [4]. We find that PS holes grow nonlinearly, independently of linear stability. Numerical simulations clarify the physics of nonlinear instabilities in both subcritical and supercritical conditions. When many resonances are unstable, several holes can coalesce into one main macro-scale structure, which survives much longer than a quasilinear diffusion time, suggesting that it may be crucial to resolve phase-space turbulence in analytical and numerical studies of transport. These findings are applied to two fundamental paradigms of plasma physics: bump-on-tail instabilities in 1D electronic plasma and current-driven ion-acoustic instabilities electron-ion plasma. Our results expose important limits of routinely-used linear and quasilinear theories.[4pt] [1] T.H. Dupree, Phys. Fluids 15, 334 (1972); R.H. Berman et al., Phys. Rev. Lett. 48, 1249 (1982).[0pt] [2] P.W. Terry, P.H. Diamond, and T.S. Hahm, Phys. Fluids B 2, 2048 (1990).[0pt] [3] H. Biglari et al., Phys. Fluids 31, 2644 (1988); Y. Kosuga et al., Phys. Plasmas 18, 122305 (2011).[0pt] [4] M. Lesur, P.H. Diamond, submitted to Phys. Rev. Lett.

  9. Nonlinear sigma models with compact hyperbolic target spaces

    NASA Astrophysics Data System (ADS)

    Gubser, Steven; Saleem, Zain H.; Schoenholz, Samuel S.; Stoica, Bogdan; Stokes, James

    2016-06-01

    We explore the phase structure of nonlinear sigma models with target spaces corresponding to compact quotients of hyperbolic space, focusing on the case of a hyperbolic genus-2 Riemann surface. The continuum theory of these models can be approximated by a lattice spin system which we simulate using Monte Carlo methods. The target space possesses interesting geometric and topological properties which are reflected in novel features of the sigma model. In particular, we observe a topological phase transition at a critical temperature, above which vortices proliferate, reminiscent of the Kosterlitz-Thouless phase transition in the O(2) model [1, 2]. Unlike in the O(2) case, there are many different types of vortices, suggesting a possible analogy to the Hagedorn treatment of statistical mechanics of a proliferating number of hadron species. Below the critical temperature the spins cluster around six special points in the target space known as Weierstrass points. The diversity of compact hyperbolic manifolds suggests that our model is only the simplest example of a broad class of statistical mechanical models whose main features can be understood essentially in geometric terms.

  10. Domain decomposition based iterative methods for nonlinear elliptic finite element problems

    SciTech Connect

    Cai, X.C.

    1994-12-31

    The class of overlapping Schwarz algorithms has been extensively studied for linear elliptic finite element problems. In this presentation, the author considers the solution of systems of nonlinear algebraic equations arising from the finite element discretization of some nonlinear elliptic equations. Several overlapping Schwarz algorithms, including the additive and multiplicative versions, with inexact Newton acceleration will be discussed. The author shows that the convergence rate of the Newton`s method is independent of the mesh size used in the finite element discretization, and also independent of the number of subdomains into which the original domain in decomposed. Numerical examples will be presented.

  11. Nonlinear Potential Model of Space-Charge Electron Beams

    NASA Astrophysics Data System (ADS)

    Litz, Marc Stuart

    1995-01-01

    This body of work is new and comprises theoretical analysis, numerical simulation and experimental investigations of the vircator, a tunable, compact, simply constructed, high power reflexing electron device that is used as a microwave source. A 1D theoretical model is formulated that is based on a time-varying, nonlinear potential to better understand the sensitivities of individual electron trajectories to macroscopic parameters of the system. Self -consistent 2D (space) and 3D(velocity) numerical PIC simulations are used to find common dynamical behavior that links the individual test-particle trajectories examined in the 1D nonlinear model with the experimental measurements of the final state of the electron trajectories. The experiment designed and built for these studies is unique and contrasts with other virtual cathode electron beam sources because it operates at lower voltage, and longer pulse width than previous forms of the device. The anode-cathode (AK) gap spacing is externally tunable expediting parametric studies. In addition, it is repetitively pulsed, which permits discrimination between single-shot anomalies and repeatable variation that often occurs in pulsed power experimental devices. These attributes distinguish this experimental apparatus from previous experimental platforms. This novel, repetitively pulsed device has particular relevance to the industrial and medical applications of rf. The approach followed in the present studies is (1) to investigate the dynamic behavior of individual test -particles in a simplified 1D, nonlinear, time-varying, potential; (2) utilize a 2D, self-consistent, electromagnetic PIC code to generate single particle trajectories subject to the collective effects of the electron beam; and (3) measure macroscopic quantities of voltage, current, electromagnetic fields, and electron flux in an experimental platform. Results of the 3D experiments are compared to predictions of the 1D and 2D models. Electron flux

  12. Image Restoration and Decomposition via Bounded Total Variation and Negative Hilbert-Sobolev Spaces

    SciTech Connect

    Lieu, Linh H. Vese, Luminita A.

    2008-10-15

    We propose a new class of models for image restoration and decomposition by functional minimization. Following ideas of Y. Meyer in a total variation minimization framework of L. Rudin, S. Osher, and E. Fatemi, our model decomposes a given (degraded or textured) image u{sub 0} into a sum u+v. Here u element of BV is a function of bounded variation (a cartoon component), while the noisy (or textured) component v is modeled by tempered distributions belonging to the negative Hilbert-Sobolev space H{sup -s}. The proposed models can be seen as generalizations of a model proposed by S. Osher, A. Sole, L. Vese and have been also motivated by D. Mumford and B. Gidas. We present existence, uniqueness and two characterizations of minimizers using duality and the notion of convex functions of measures with linear growth, following I. Ekeland and R. Temam, F. Demengel and R. Temam. We also give a numerical algorithm for solving the minimization problem, and we present numerical results of denoising, deblurring, and decompositions of both synthetic and real images.

  13. Domain Decomposition PN Solutions to the 3D Transport Benchmark over a Range in Parameter Space

    NASA Astrophysics Data System (ADS)

    Van Criekingen, S.

    2014-06-01

    The objectives of this contribution are twofold. First, the Domain Decomposition (DD) method used in the parafish parallel transport solver is re-interpreted as a Generalized Schwarz Splitting as defined by Tang [SIAM J Sci Stat Comput, vol.13 (2), pp. 573-595, 1992]. Second, parafish provides spherical harmonic (i.e., PN) solutions to the NEA benchmark suite for 3D transport methods and codes over a range in parameter space. To the best of the author's knowledge, these are the first spherical harmonic solutions provided for this demanding benchmark suite. They have been obtained using 512 CPU cores of the JuRoPa machine installed at the Jülich Computing Center (Germany).

  14. Localized spatially nonlinear matter waves in atomic-molecular Bose-Einstein condensates with space-modulated nonlinearity.

    PubMed

    Yao, Yu-Qin; Li, Ji; Han, Wei; Wang, Deng-Shan; Liu, Wu-Ming

    2016-01-01

    The intrinsic nonlinearity is the most remarkable characteristic of the Bose-Einstein condensates (BECs) systems. Many studies have been done on atomic BECs with time- and space- modulated nonlinearities, while there is few work considering the atomic-molecular BECs with space-modulated nonlinearities. Here, we obtain two kinds of Jacobi elliptic solutions and a family of rational solutions of the atomic-molecular BECs with trapping potential and space-modulated nonlinearity and consider the effect of three-body interaction on the localized matter wave solutions. The topological properties of the localized nonlinear matter wave for no coupling are analysed: the parity of nonlinear matter wave functions depends only on the principal quantum number n, and the numbers of the density packets for each quantum state depend on both the principal quantum number n and the secondary quantum number l. When the coupling is not zero, the localized nonlinear matter waves given by the rational function, their topological properties are independent of the principal quantum number n, only depend on the secondary quantum number l. The Raman detuning and the chemical potential can change the number and the shape of the density packets. The stability of the Jacobi elliptic solutions depends on the principal quantum number n, while the stability of the rational solutions depends on the chemical potential and Raman detuning. PMID:27403634

  15. Localized spatially nonlinear matter waves in atomic-molecular Bose-Einstein condensates with space-modulated nonlinearity

    PubMed Central

    Yao, Yu-Qin; Li, Ji; Han, Wei; Wang, Deng-Shan; Liu, Wu-Ming

    2016-01-01

    The intrinsic nonlinearity is the most remarkable characteristic of the Bose-Einstein condensates (BECs) systems. Many studies have been done on atomic BECs with time- and space- modulated nonlinearities, while there is few work considering the atomic-molecular BECs with space-modulated nonlinearities. Here, we obtain two kinds of Jacobi elliptic solutions and a family of rational solutions of the atomic-molecular BECs with trapping potential and space-modulated nonlinearity and consider the effect of three-body interaction on the localized matter wave solutions. The topological properties of the localized nonlinear matter wave for no coupling are analysed: the parity of nonlinear matter wave functions depends only on the principal quantum number n, and the numbers of the density packets for each quantum state depend on both the principal quantum number n and the secondary quantum number l. When the coupling is not zero, the localized nonlinear matter waves given by the rational function, their topological properties are independent of the principal quantum number n, only depend on the secondary quantum number l. The Raman detuning and the chemical potential can change the number and the shape of the density packets. The stability of the Jacobi elliptic solutions depends on the principal quantum number n, while the stability of the rational solutions depends on the chemical potential and Raman detuning. PMID:27403634

  16. FAST TRACK COMMUNICATION: Solvable nonlinear evolution PDEs in multidimensional space involving trigonometric functions

    NASA Astrophysics Data System (ADS)

    Calogero, F.; Françoise, J.-P.; Sommacal, M.

    2007-05-01

    A solvable nonlinear (system of) evolution PDEs in multidimensional space, involving trigonometric (or hyperbolic) functions, is identified. An isochronous version of this (system of) evolution PDEs in multidimensional space is also reported.

  17. Effect of Nonlinear Joints on Space Deployable Truss Structures

    NASA Astrophysics Data System (ADS)

    Guo, Hongwei; Deng, Zhongquan; Wu, Xiang; Liu, Rongqiang

    2012-07-01

    Joints nonlinearities with characteristics of freeplay and hysteresis are analyzed by describing joint nonlinear force-displacement based on describing function method. The nonlinear dynamic responses of the one- DOF system with joints under different exciting force levels are presented in the charts. The influence of the characterizing parameters, e.g., gaps, slipping forces of the joints on nonlinearities is analyzed. The nonlinear effects of freeplay and hysteresis present that the dynamic responses switch from one resonance frequency to another frequency when amplitude exceed the demarcation values. The hysteresis nonlinearity contributes nonlinear damping to the system. Dynamic responses of the modular beam-like deployable joint- dominated truss structure are tested under different sinusoidal exciting force levels which show obvious nonlinear behaviors. The nonlinear dynamic behaviors of the truss structure contributed by the joints shows a shift to lower resonance frequency and higher amplitude with the exciting force increases. The nonlinearity of the joints in the tested structure is identified to meet with the hysteresis nonlinearity. The experiment validates that describing method is an effective tool to model the joint nonlinearities.

  18. Color decomposition method for multiprimary display using 3D-LUT in linearized LAB space

    NASA Astrophysics Data System (ADS)

    Kang, Dong-Woo; Kim, Yun-Tae; Cho, Yang-Ho; Park, Kee-Hyon; Choe, Wonhee; Ha, Yeong-Ho

    2005-01-01

    This paper proposes a color decomposition method for a multi-primary display (MPD) using a 3-dimensional look-up-table (3D-LUT) in linearized LAB space. The proposed method decomposes the conventional three primary colors into multi-primary control values for a display device under the constraints of tristimulus matching. To reproduce images on an MPD, the color signals are estimated from a device-independent color space, such as CIEXYZ and CIELAB. In this paper, linearized LAB space is used due to its linearity and additivity in color conversion. First, the proposed method constructs a 3-D LUT containing gamut boundary information to calculate the color signals for the MPD in linearized LAB space. For the image reproduction, standard RGB or CIEXYZ is transformed to linearized LAB, then the hue and chroma are computed with reference to the 3D-LUT. In linearized LAB space, the color signals for a gamut boundary point are calculated to have the same lightness and hue as the input point. Also, the color signals for a point on the gray axis are calculated to have the same lightness as the input point. Based on the gamut boundary points and input point, the color signals for the input point are then obtained using the chroma ratio divided by the chroma of the gamut boundary point. In particular, for a change of hue, the neighboring boundary points are also employed. As a result, the proposed method guarantees color signal continuity and computational efficiency, and requires less memory.

  19. Color decomposition method for multiprimary display using 3D-LUT in linearized LAB space

    NASA Astrophysics Data System (ADS)

    Kang, Dong-Woo; Kim, Yun-Tae; Cho, Yang-Ho; Park, Kee-Hyon; Choe, Wonhee; Ha, Yeong-Ho

    2004-12-01

    This paper proposes a color decomposition method for a multi-primary display (MPD) using a 3-dimensional look-up-table (3D-LUT) in linearized LAB space. The proposed method decomposes the conventional three primary colors into multi-primary control values for a display device under the constraints of tristimulus matching. To reproduce images on an MPD, the color signals are estimated from a device-independent color space, such as CIEXYZ and CIELAB. In this paper, linearized LAB space is used due to its linearity and additivity in color conversion. First, the proposed method constructs a 3-D LUT containing gamut boundary information to calculate the color signals for the MPD in linearized LAB space. For the image reproduction, standard RGB or CIEXYZ is transformed to linearized LAB, then the hue and chroma are computed with reference to the 3D-LUT. In linearized LAB space, the color signals for a gamut boundary point are calculated to have the same lightness and hue as the input point. Also, the color signals for a point on the gray axis are calculated to have the same lightness as the input point. Based on the gamut boundary points and input point, the color signals for the input point are then obtained using the chroma ratio divided by the chroma of the gamut boundary point. In particular, for a change of hue, the neighboring boundary points are also employed. As a result, the proposed method guarantees color signal continuity and computational efficiency, and requires less memory.

  20. Null Space Imaging: Nonlinear Magnetic Encoding Fields Designed Complementary to Receiver Coil Sensitivities for Improved Acceleration in Parallel Imaging

    PubMed Central

    Tam, Leo K.; Galiana, Gigi; Stockmann, Jason P.; Constable, R. Todd

    2012-01-01

    To increase image acquisition efficiency, we develop alternative gradient encoding strategies designed to provide spatial encoding complementary to the spatial encoding provided by the multiple receiver coil elements in parallel image acquisitions. Intuitively, complementary encoding is achieved when the magnetic field encoding gradients are designed to encode spatial information where receiver spatial encoding is ambiguous, for example, along sensitivity isocontours. Specifically, the method generates a basis set for the null space of the coil sensitivities with the singular value decomposition (SVD) and calculates encoding fields from the null space vectors. A set of nonlinear gradients is used as projection imaging readout magnetic fields, replacing the conventional linear readout field and phase encoding. Multiple encoding fields are used as projections to capture the null space information, hence the term Null Space Imaging (NSI). The method is compared to conventional Cartesian SENSitivity Encoding (SENSE) as evaluated by mean squared error and robustness to noise. Strategies for developments in the area of nonlinear encoding schemes are discussed. The NSI approach yields a parallel imaging method that provides high acceleration factors with a limited number of receiver coil array elements through increased time efficiency in spatial encoding. PMID:22190380

  1. Coherent Control of Molecular Torsion and the Active-space Decomposition Method

    NASA Astrophysics Data System (ADS)

    Parker, Shane Matthew

    This dissertation discusses schemes and applications for the strong-field control of molecular torsions as well as introduces the active-space decomposition method. In the first part, a route to realize general control over the torsional motions of a class of biaryl compounds is proposed. Torsion in biaryl compounds--molecules with two aromatic moieties connected by a bond about which the barrier to rotation is small--mediates the electronic coupling between the two rings in the molecule. Thus, by controlling the torsion angle, one also controls the electron transfer and transport rates, the absorption and emission spectra, and the molecule's chirality. In our scheme, a non-resonant half-cycle pulse interacts with the permanent dipole of only one moiety of the pre-oriented biaryl compound. In the non-adiabatic regime, coherent motion is initiated by the half-cycle pulse. In the adiabatic regime, the torsion angle is tuned by the pulse. By properly choosing the parameters and polarization of the half-cycle pulse, we show that free internal rotation can be started or that the molecular chirality can be inverted. Then, with the aid of optimal control theory, we design "deracemizing" control pulses, i.e., control pulses that convert a racemic mixture into an enantiopure mixture. Finally, we explore the potential for this type of control in a single-molecule pulling experiment. In the second part, we describe the active space decomposition method for computing excited states of molecular dimers. In this method, the dimer's wavefunction is expressed as a linear combination of direct products of orthogonal localized monomer states. The adiabatic dimer states are found by diagonalizing the Hamiltonian in this direct product space. Matrix elements between direct product states are computed directly, without ever explicitly forming the dimer wavefunction, thus enabling calculations of dimers with active space sizes that would be otherwise impossible. The decomposed

  2. Rapid Transient Pressure Field Computations in the Nearfield of Circular Transducers using Frequency Domain Time-Space Decomposition

    PubMed Central

    Alles, E. J.; Zhu, Y.; van Dongen, K. W. A.; McGough, R. J.

    2013-01-01

    The fast nearfield method, when combined with time-space decomposition, is a rapid and accurate approach for calculating transient nearfield pressures generated by ultrasound transducers. However, the standard time-space decomposition approach is only applicable to certain analytical representations of the temporal transducer surface velocity that, when applied to the fast nearfield method, are expressed as a finite sum of products of separate temporal and spatial terms. To extend time-space decomposition such that accelerated transient field simulations are enabled in the nearfield for an arbitrary transducer surface velocity, a new transient simulation method, frequency domain time-space decomposition (FDTSD), is derived. With this method, the temporal transducer surface velocity is transformed into the frequency domain, and then each complex-valued term is processed separately. Further improvements are achieved by spectral clipping, which reduces the number of terms and the computation time. Trade-offs between speed and accuracy are established for FDTSD calculations, and pressure fields obtained with the FDTSD method for a circular transducer are compared to those obtained with Field II and the impulse response method. The FDTSD approach, when combined with the fast nearfield method and spectral clipping, consistently achieves smaller errors in less time and requires less memory than Field II or the impulse response method. PMID:23160476

  3. De Rham-Hodge decomposition and vanishing of harmonic forms by derivation operators on the Poisson space

    NASA Astrophysics Data System (ADS)

    Privault, Nicolas

    2016-05-01

    We construct differential forms of all orders and a covariant derivative together with its adjoint on the probability space of a standard Poisson process, using derivation operators. In this framewok we derive a de Rham-Hodge-Kodaira decomposition as well as Weitzenböck and Clark-Ocone formulas for random differential forms. As in the Wiener space setting, this construction provides two distinct approaches to the vanishing of harmonic differential forms.

  4. Application of 2D-Nonlinear Shallow Water Model of Tsunami by using Adomian Decomposition Method

    SciTech Connect

    Waewcharoen, Sribudh; Boonyapibanwong, Supachai; Koonprasert, Sanoe

    2008-09-01

    One of the most important questions in tsunami modeling is the estimation of tsunami run-up heights at different points along a coastline. Methods for numerical simulation of tsunami wave propagation in deep and shallow seas are well developed and have been widely used by many scientists (2001-2008). In this paper, we consider a two-dimensional nonlinear shallow water model of tsunami given by Tivon Jacobson is work [1]. u{sub t}+uu{sub x}+{nu}u{sub y} -c{sup 2}(h{sub x}+(h{sub b}){sub x}) {nu}{sub t}+u{nu}{sub x}+{nu}{nu}{sub y} = -c{sup 2}(h{sub y}+(h{sub b}){sub y}) h{sub t}+(hu){sub x}+(h{nu}){sub y} = 0 g-shore, h is surface elevation and s, t is time, u is velocity of cross-shore, {nu} is velocity of along-shore, h is surface elevation and h{sub b} is function of shore. This is a nondimensionalized model with the gravity g and constant reference depth H factored into c = {radical}(gH). We apply the Adomian Decompostion Method (ADM) to solve the tsunami model. This powerful method has been used to obtain explicit and numerical solutions of three types of diffusion-convection-reaction (DECR) equations. The ADM results for the tsunami model yield analytical solutions in terms of a rapidly convergent infinite power series. Symbolic computation, numerical results and graphs of solutions are obtained by Maple program.

  5. Verification of Multiphysics software: Space and time convergence studies for nonlinearly coupled applications

    SciTech Connect

    Jean C. Ragusa; Vijay Mahadevan; Vincent A. Mousseau

    2009-05-01

    High-fidelity modeling of nuclear reactors requires the solution of a nonlinear coupled multi-physics stiff problem with widely varying time and length scales that need to be resolved correctly. A numerical method that converges the implicit nonlinear terms to a small tolerance is often referred to as nonlinearly consistent (or tightly coupled). This nonlinear consistency is still lacking in the vast majority of coupling techniques today. We present a tightly coupled multiphysics framework that tackles this issue and present code-verification and convergence analyses in space and time for several models of nonlinear coupled physics.

  6. Dynamic mode decomposition of supersonic and transonic wakes of generic space launcher configurations

    NASA Astrophysics Data System (ADS)

    Statnikov, Vladimir; Sayadi, Taraneh; Meinke, Matthias; Schroeder, Wolfgang; Schmid, Peter

    2013-11-01

    Dynamic mode decomposition (DMD) is applied to supersonic and transonic wakes of generic space launcher configurations of Mach numbers 0.7 and 6 computed using a zonal RANS/LES approach. The axisymmetric geometry includes a backward facing step that causes the flow to separate. In addition to the separation bubble, acoustic waves are also radiated from the downstream region of the flow. Experimental and numerical observations clearly demonstrate the existence of peaks in the pressure spectra which can be attributed to both the flow inside the wake and the acoustic waves in the freestream. The objective of this work is to apply DMD to the set of numerical data in order to firstly, extract the spatial shape of the modes and secondly, identify their respective frequencies. This allows the dynamics associated to the separation bubble and those of the acoustic waves to be differentiated properly. In addition, with the help of DMD the modes responsible for pressure-loading on the backward face of the step are extracted and analyzed.

  7. Localized Nonlinear Waves in Systems with Time- and Space-Modulated Nonlinearities

    SciTech Connect

    Belmonte-Beitia, Juan; Perez-Garcia, Victor M.; Vekslerchik, Vadym; Konotop, Vladimir V.

    2008-04-25

    Using similarity transformations we construct explicit nontrivial solutions of nonlinear Schroedinger equations with potentials and nonlinearities depending both on time and on the spatial coordinates. We present the general theory and use it to calculate explicitly nontrivial solutions such as periodic (breathers), resonant, or quasiperiodically oscillating solitons. Some implications to the field of matter waves are also discussed.

  8. New Solutions of Three Nonlinear Space- and Time-Fractional Partial Differential Equations in Mathematical Physics

    NASA Astrophysics Data System (ADS)

    Yao, Ruo-Xia; Wang, Wei; Chen, Ting-Hua

    2014-11-01

    Motivated by the widely used ansätz method and starting from the modified Riemann—Liouville derivative together with a fractional complex transformation that can be utilized to transform nonlinear fractional partial differential equations to nonlinear ordinary differential equations, new types of exact traveling wave solutions to three important nonlinear space- and time-fractional partial differential equations are obtained simultaneously in terms of solutions of a Riccati equation. The results are new and first reported in this paper.

  9. Special class of nonlinear damping models in flexible space structures

    NASA Technical Reports Server (NTRS)

    Hu, Anren; Singh, Ramendra P.; Taylor, Lawrence W.

    1991-01-01

    A special class of nonlinear damping models is investigated in which the damping force is proportional to the product of positive integer or the fractional power of the absolute values of displacement and velocity. For a one-degree-of-freedom system, the classical Krylov-Bogoliubov 'averaging' method is used, whereas for a distributed system, both an ad hoc perturbation technique and the finite difference method are employed to study the effects of nonlinear damping. The results are compared with linear viscous damping models. The amplitude decrement of free vibration for a single mode system with nonlinear models depends not only on the damping ratio but also on the initial amplitude, the time to measure the response, the frequency of the system, and the powers of displacement and velocity. For the distributed system, the action of nonlinear damping is found to reduce the energy of the system and to pass energy to lower modes.

  10. Effects of nonlinear damping in flexible space structures

    NASA Technical Reports Server (NTRS)

    Hu, Anren; Taylor, Lawrence W.

    1988-01-01

    The classical Krylov-Bogoliubov "averaging" technique is used to study a class of nonlinear damping models, for which the damping force is proportional to the product of positive integer or fractional power of absolute values of displacement and that of velocity. The results are compared with linear viscous damping models. The amplitude decrement of free vibration for a single mode system with nonlinear models depends not only on damping ratio, but also on the initial amplitude, the time to measure the response, frequency of the system, and the powers of displacement and velocity. For the distributed system, the action of nonlinear damping is found to reduce energy of the system as well as to pass energy to higher modes. Experimental evidence such as in Spacecraft Control Laboratory Experiment seems to support the need for nonlinear models.

  11. Analytical Solution of the Space-Time Fractional Nonlinear Schrödinger Equation

    NASA Astrophysics Data System (ADS)

    Abdel-Salam, Emad A.-B.; Yousif, Eltayeb A.; El-Aasser, Mostafa A.

    2016-02-01

    The space-time fractional nonlinear Schrödinger equation is solved by mean of on the fractional Riccati expansion method. These solutions include generalized trigonometric and hyperbolic functions which could be useful for further understanding of mechanisms of the complicated nonlinear physical phenomena and fractional differential equations. Among these solutions, some are found for the first time.

  12. Nonlinear Dirac equation in two-spinor form: Separation in static RW space-time. Nonlinear Dirac equation in RW space-time

    NASA Astrophysics Data System (ADS)

    Zecca, Antonio

    2016-02-01

    The Dirac equation with nonlinear terms induced by torsion is studied in Robertson-Walker (RW) space-time. An extension of a separation method of the equation, based on the Newman-Penrose formalism and previously applied to the nonlinear case, is considered. Accordingly the angular dependence of the Dirac spinor solution is separated, under a special assumption, in the general time-dependent RW metric. In the case of static RW metric the time dependence of the Dirac spinor factors out and one is left with a pair of two coupled nonlinear radial equations. The radial equations are disentangled by a suitable substitution of the spinor solution. The problem amounts then to the solution of a single second-order highly nonlinear differential equation. Some elementary considerations are done on the asymptotic behavior of the solution of the equation.

  13. FAST TRACK COMMUNICATION: \\ {P}\\ {T}-symmetry, Cartan decompositions, Lie triple systems and Krein space-related Clifford algebras

    NASA Astrophysics Data System (ADS)

    Günther, Uwe; Kuzhel, Sergii

    2010-10-01

    Gauged \\ {P}\\ {T} quantum mechanics (PTQM) and corresponding Krein space setups are studied. For models with constant non-Abelian gauge potentials and extended parity inversions compact and noncompact Lie group components are analyzed via Cartan decompositions. A Lie-triple structure is found and an interpretation as \\ {P}\\ {T}-symmetrically generalized Jaynes-Cummings model is possible with close relation to recently studied cavity QED setups with transmon states in multilevel artificial atoms. For models with Abelian gauge potentials a hidden Clifford algebra structure is found and used to obtain the fundamental symmetry of Krein space-related J-self-adjoint extensions for PTQM setups with ultra-localized potentials.

  14. The nonlinear electromigration of analytes into confined spaces

    PubMed Central

    Chen, Zhen; Ghosal, Sandip

    2012-01-01

    We consider the problem of electromigration of a sample ion (analyte) within a uniform background electrolyte when the confining channel undergoes a sudden contraction. One example of such a situation arises in microfluidics in the electrokinetic injection of the analyte into a micro-capillary from a reservoir of much larger size. Here, the sample concentration propagates as a wave driven by the electric field. The dynamics is governed by the Nerst–Planck–Poisson system of equations for ionic transport. A reduced one-dimensional nonlinear equation, describing the evolution of the sample concentration, is derived. We integrate this equation numerically to obtain the evolution of the wave shape and determine how the injected mass depends on the sample concentration in the reservoir. It is shown that due to the nonlinear coupling of the ionic concentrations and the electric field, the concentration of the injected sample could be substantially less than the concentration of the sample in the reservoir. PMID:22977354

  15. Space formations and nonlinear properties of noncentrosymmetric germanates

    SciTech Connect

    Korotkov, Anton Sergeevich

    2014-10-15

    Space formations of Ge-O polyhedra have been analyzed for 114 noncentrosymmetric germanates. The type of Ge-O polyhedra space formations is dependent on the stoichiometric ratio SR=n(O)/n(Ge), where n(O) is the number of oxygen and n(Ge) is the number of germanium atoms in formal composition of the compound. Individual (Ge-O) polyhedra are found for SR≥3 chains of the polyhedra that exist in the range of SR=3–5. Only framework formations are possible at SR≤2.7. - graphical abstract: Examples of space formations of polyhedra. - Highlights: • Space formations of Ge-O polyhedra have been analyzed. • The dependence of the type of Ge-O polyhedra space formations on the stoichiometric ratio SR=n(O)/n(Ge) was established. • Set of experimental SHG data of noncentrosymmetric germanates was presented.

  16. Image enhancement by nonlinear extrapolation in frequency space

    NASA Astrophysics Data System (ADS)

    Greenspan, Hayit; Anderson, Charles H.

    1994-03-01

    A procedure for creating images with higher resolution than the sampling rate would allow is described. The enhancement algorithm augments the frequency content of the image using shape-invariant properties of edges across scale by using a non-linearity that generates phase- coherent higher harmonics. The procedure utilizes the Laplacian pyramid image representation. Results are presented depicting the power-spectra augmentation and the visual enhancement of several images. Simplicity of computations and ease of implementation allow for real-time applications such as high-definition television.

  17. NONLINEAR BEHAVIOR OF BARYON ACOUSTIC OSCILLATIONS IN REDSHIFT SPACE FROM THE ZEL'DOVICH APPROXIMATION

    SciTech Connect

    McCullagh, Nuala; Szalay, Alexander S.

    2015-01-10

    Baryon acoustic oscillations (BAO) are a powerful probe of the expansion history of the universe, which can tell us about the nature of dark energy. In order to accurately characterize the dark energy equation of state using BAO, we must understand the effects of both nonlinearities and redshift space distortions on the location and shape of the acoustic peak. In a previous paper, we introduced a novel approach to second order perturbation theory in configuration space using the Zel'dovich approximation, and presented a simple result for the first nonlinear term of the correlation function. In this paper, we extend this approach to redshift space. We show how to perform the computation and present the analytic result for the first nonlinear term in the correlation function. Finally, we validate our result through comparison with numerical simulations.

  18. Asynchronous space-time algorithm based on a domain decomposition method for structural dynamics problems on non-matching meshes

    NASA Astrophysics Data System (ADS)

    Subber, Waad; Matouš, Karel

    2016-02-01

    Large-scale practical engineering problems featuring localized phenomena often benefit from local control of mesh and time resolutions to efficiently capture the spatial and temporal scales of interest. To this end, we propose an asynchronous space-time algorithm based on a domain decomposition method for structural dynamics problems on non-matching meshes. The three-field algorithm is based on the dual-primal like domain decomposition approach utilizing the localized Lagrange multipliers along the space and time common-refinement-based interface. The proposed algorithm is parallel in nature and well suited for a heterogeneous computing environment. Moreover, two-levels of parallelism are embedded in this novel scheme. For linear dynamical problems, the algorithm is unconditionally stable, shows an optimal order of convergence with respect to space and time discretizations as well as ensures conservation of mass, momentum and energy across the non-matching grid interfaces. The method of manufactured solutions is used to verify the implementation, and an engineering application is considered, where a sandwich plate is impacted by a projectile.

  19. Image enhancement by non-linear extrapolation in frequency space

    NASA Technical Reports Server (NTRS)

    Anderson, Charles H. (Inventor); Greenspan, Hayit K. (Inventor)

    1998-01-01

    An input image is enhanced to include spatial frequency components having frequencies higher than those in an input image. To this end, an edge map is generated from the input image using a high band pass filtering technique. An enhancing map is subsequently generated from the edge map, with the enhanced map having spatial frequencies exceeding an initial maximum spatial frequency of the input image. The enhanced map is generated by applying a non-linear operator to the edge map in a manner which preserves the phase transitions of the edges of the input image. The enhanced map is added to the input image to achieve a resulting image having spatial frequencies greater than those in the input image. Simplicity of computations and ease of implementation allow for image sharpening after enlargement and for real-time applications such as videophones, advanced definition television, zooming, and restoration of old motion pictures.

  20. Effects of Transverse Physics on Nonlinear Evolution of Longitudinal Space-Charge Waves in Beams

    SciTech Connect

    K. Tian, I. Haber, R.A. Kishek, P.G. O'Shea, M. Reiser, D. Stratakis

    2009-05-01

    Longitudinal space-charge waves can introduce energy perturbations into charge particle beams and degrade the beam quality, which is critical to many modern applications of particle accelerators. Although many longitudinal phenomena arising from small perturbations can be explained by a one-dimensional cold fluid theory, nonlinear behavior of space-charge waves observed in experiments has not been well understood. In this paper, we summarize our recent investigation by means of more detailed measurements and self-consistent simulations. Combining the numerical capability of a PIC code, WARP, with the detailed initial conditions measured by our newly developed time resolved 6-D phase space mapping technique, we are able to construct a self consistent model for studying the complex physics of longitudinal dynamics of space-charge dominated beams. Results from simulation studies suggest that the unexplained nonlinear behavior of space-charge waves may be due to transverse mismatch or misalignment of beams.

  1. Task decomposition for a multilimbed robot to work in reachable but unorientable space

    NASA Technical Reports Server (NTRS)

    Su, Chau; Zheng, Yuan F.

    1991-01-01

    Robot manipulators installed on legged mobile platforms are suggested for enlarging robot workspace. To plan the motion of such a system, the arm-platform motion coordination problem is raised, and a task decomposition is proposed to solve the problem. A given task described by the destination position and orientation of the end effector is decomposed into subtasks for arm manipulation and for platform configuration, respectively. The former is defined as the end-effector position and orientation with respect to the platform, and the latter as the platform position and orientation in the base coordinates. Three approaches are proposed for the task decomposition. The approaches are also evaluated in terms of the displacements, from which an optimal approach can be selected.

  2. Study of nonlinear oscillations in a glow discharge plasma using empirical mode decomposition and Hilbert Huang transform

    SciTech Connect

    Wharton, A. M.; Sekar Iyengar, A. N.; Janaki, M. S.

    2013-02-15

    Hilbert Huang transform (HHT) based time series analysis was carried out on nonlinear floating potential fluctuations obtained from hollow cathode glow discharge plasma in the presence of anode glow. HHT was used to obtain contour plots and the presence of nonlinearity was studied. Frequency shift with time, which is a typical nonlinear behaviour, was detected from the contour plots. Various plasma parameters were measured and the concepts of correlation coefficients and the physical contribution of each intrinsic mode function have been discussed. Physically important quantities such as instantaneous energy and their uses in studying physical phenomena such as intermittency and non-stationary data have also been discussed.

  3. When Policy Structures Technology: Balancing upfront decomposition and in-process coordination in Europe's decentralized space technology ecosystem

    NASA Astrophysics Data System (ADS)

    Vrolijk, Ademir; Szajnfarber, Zoe

    2015-01-01

    This paper examines the decentralization of European space technology research and development through the joint lenses of policy, systems architecture, and innovation contexts. It uses a detailed longitudinal case history of the development of a novel astrophysics instrument to explore the link between policy-imposed institutional decomposition and the architecture of the technical system. The analysis focuses on five instances of collaborative design decision-making and finds that matching between the technical and institutional architectures is a predictor of project success, consistent with the mirroring hypothesis in extant literature. Examined over time, the instances reveal stability in the loosely coupled nature of institutional arrangements and a trend towards more integral, or tightly coupled, technical systems. The stability of the institutional arrangements is explained as an artifact of the European Hultqvist policy and the trend towards integral technical systems is related to the increasing complexity of modern space systems. If these trends persist, the scale of the mismatch will continue to grow. As a first step towards mitigating this challenge, the paper develops a framework for balancing upfront decomposition and in-process coordination in collaborative development projects. The astrophysics instrument case history is used to illustrate how collaborations should be defined for a given inherent system complexity.

  4. Suppression of space charge induced beam halo in nonlinear focusing channel

    NASA Astrophysics Data System (ADS)

    Batygin, Yuri K.; Scheinker, Alexander; Kurennoy, Sergey; Li, Chao

    2016-04-01

    An intense non-uniform particle beam exhibits strong emittance growth and halo formation in focusing channels due to nonlinear space charge forces of the beam. This phenomenon limits beam brightness and results in particle losses. The problem is connected with irreversible distortion of phase space volume of the beam in conventional focusing structures due to filamentation in phase space. Emittance growth is accompanied by halo formation in real space, which results in inevitable particle losses. A new approach for solving a self-consistent problem for a matched non-uniform beam in two-dimensional geometry is discussed. The resulting solution is applied to the problem of beam transport, while avoiding emittance growth and halo formation by the use of nonlinear focusing field. Conservation of a beam distribution function is demonstrated analytically and by particle-in-cell simulation for a beam with a realistic beam distribution.

  5. Near equilibrium distributions for beams with space charge in linear and nonlinear periodic focusing systems

    SciTech Connect

    Sonnad, Kiran G.; Cary, John R.

    2015-04-15

    A procedure to obtain a near equilibrium phase space distribution function has been derived for beams with space charge effects in a generalized periodic focusing transport channel. The method utilizes the Lie transform perturbation theory to canonically transform to slowly oscillating phase space coordinates. The procedure results in transforming the periodic focusing system to a constant focusing one, where equilibrium distributions can be found. Transforming back to the original phase space coordinates yields an equilibrium distribution function corresponding to a constant focusing system along with perturbations resulting from the periodicity in the focusing. Examples used here include linear and nonlinear alternating gradient focusing systems. It is shown that the nonlinear focusing components can be chosen such that the system is close to integrability. The equilibrium distribution functions are numerically calculated, and their properties associated with the corresponding focusing system are discussed.

  6. On the analytical modeling of the nonlinear vibrations of pretensioned space structures

    NASA Technical Reports Server (NTRS)

    Housner, J. M.; Belvin, W. K.

    1983-01-01

    Pretensioned structures are receiving considerable attention as candidate large space structures. A typical example is a hoop-column antenna. The large number of preloaded members requires efficient analytical methods for concept validation and design. Validation through analyses is especially important since ground testing may be limited due to gravity effects and structural size. The present investigation has the objective to present an examination of the analytical modeling of pretensioned members undergoing nonlinear vibrations. Two approximate nonlinear analysis are developed to model general structural arrangements which include beam-columns and pretensioned cables attached to a common nucleus, such as may occur at a joint of a pretensioned structure. Attention is given to structures undergoing nonlinear steady-state oscillations due to sinusoidal excitation forces. Three analyses, linear, quasi-linear, and nonlinear are conducted and applied to study the response of a relatively simple cable stiffened structure.

  7. Nonlinear Analysis of the Space Shuttle Super-Lightweight External Fuel Tank

    NASA Technical Reports Server (NTRS)

    Nemeth, Michael P.; Britt, Vicki O.; Collins, Timothy J.; Starnes, James H., Jr.

    1996-01-01

    The results of buckling and nonlinear analyses of the Space Shuttle External Tank super-lightweight liquid oxygen (LOX) tank are presented. Modeling details and results are presented for two prelaunch loading conditions and for two full-scale structural tests conducted on the original external tank. These results illustrate three distinctly different types of nonlinear responses for thin-walled shells subjected to combined mechanical and thermal loads. These nonlinear response phenomena consist of bifurcation-type buckling, short-wavelength nonlinear bending, and nonlinear collapse associated with a limit point. For each case, the results show that accurate predictions of nonlinear behavior generally require a large scale high-fidelity finite element model. Results are also presented that show that a fluid filled launch vehicle shell can be highly sensitive to initial geometric imperfections. In addition, results presented for two full scale structural tests of the original standard weight external tank suggest that the finite element modeling approach used in the present study is sufficient for representing the nonlinear behavior of the super lightweight LOX tank.

  8. Gombosi Receives 2013 Space Weather and Nonlinear Waves and Processes Prize: Citation

    NASA Astrophysics Data System (ADS)

    McCoy, Robert P.

    2014-08-01

    The Space Physics and Aeronomy section of AGU awards the Space Weather and Nonlinear Waves and Processes Prize to Tamas Gombosi of the University of Michigan. Gombosi, the founding director of the Center for Space Environment Modeling, has been a leader in space weather research, a visionary in space weather numerical modeling for several decades, and a pioneer of international space physics collaboration. Gombosi's major contributions to space weather include the development of the first time-dependent numerical model of the terrestrial polar wind and the creation of the BATS-R-US magnetohydrodynamic model, a powerful and versatile numerical tool widely used today for modeling the global geospace, the heliosphere, and the solar interior.

  9. Proper Criteria of Nonlinear Optical Crystals for Space Laser Systems and the Possible Causes for Space Laser Failures

    NASA Technical Reports Server (NTRS)

    Abdeldayem, Hossin A.; Dowdye, Edward; Jamison, Tracee; Canham, John

    2005-01-01

    NASA is striving to develop a scientific understanding of the universe and the Earth-Sun System and its response to natural or human-induced changes. Space lasers are vital tools for NASA's missions to advance our understanding of space research and improving our prediction capability for climate, weather, and natural hazards. Unfortunately, several past space missions that utilized lasers proved to be short-lived and unreliable. In this paper, we are reporting the results of our investigations on several nonlinear optical crystals, which are vital components in space lasers. Examples of these investigations are: The correlation of the phase diagrams of nonlinear crystals and its durability, the effect of radiating these crystals by high-energy beams of protons and gamma on their second harmonic efficiency, and measurements of the high-energy and low-energy thresholds for each crystal before and after irradiation. A set of proper criteria for these crystals will be presented. We will also discuss the possible causes of failures in a space laser and propose a solution to a contamination problem in all future space lasers.

  10. Photorefractive nonlinearities caused by the Dember space-charge field in undoped CdTe.

    PubMed

    Schroeder, W A; Stark, T S; Boggess, T F; Smirl, A L; Valley, G C

    1991-06-01

    The photorefractive nonlinearity associated with the Dember space-charge field between electrons and holes produced by two-photon absorption is unambiguously isolated and studied in undoped CdTe by using a nondegenerate, forward-probing, polarization-sensitive, transient-grating technique with a temporal resolution of <5 ps. PMID:19776789

  11. Periodic boundary value problems for nonlinear impulsive evolution equations on Banach spaces

    NASA Astrophysics Data System (ADS)

    Yu, Xiulan; Wang, JinRong

    2015-05-01

    In this paper, we consider periodic boundary value problems for two new type nonlinear impulsive evolution equations on Banach spaces. By using the theory of semigroup, a mixed type Gronwall inequality and fixed point methods, we establish several sufficient conditions on the existence of mild solutions for such problems. Finally, examples are given to illustrate our main results.

  12. Baseline-free damage visualization using noncontact laser nonlinear ultrasonics and state space geometrical changes

    NASA Astrophysics Data System (ADS)

    Liu, Peipei; Sohn, Hoon; Park, Byeongjin

    2015-06-01

    Damage often causes a structural system to exhibit severe nonlinear behaviors, and the resulting nonlinear features are often much more sensitive to the damage than their linear counterparts. This study develops a laser nonlinear wave modulation spectroscopy (LNWMS) so that certain types of damage can be detected without any sensor placement. The proposed LNWMS utilizes a pulse laser to generate ultrasonic waves and a laser vibrometer for ultrasonic measurement. Under the broadband excitation of the pulse laser, a nonlinear source generates modulations at various frequency values due to interactions among various input frequency components. State space attractors are reconstructed from the ultrasonic responses measured by LNWMS, and a damage feature called Bhattacharyya distance (BD) is computed from the state space attractors to quantify the degree of damage-induced nonlinearity. By computing the BD values over the entire target surface using laser scanning, damage can be localized and visualized without relying on the baseline data obtained from the pristine condition of a target structure. The proposed technique has been successfully used for visualizing fatigue crack in an aluminum plate and delamination and debonding in a glass fiber reinforced polymer wind turbine blade.

  13. Disentangling Redshift-Space Distortions and Nonlinear Bias using the 2D Power Spectrum

    SciTech Connect

    Jennings, Elise; Wechsler, Risa H.

    2015-08-07

    We present the nonlinear 2D galaxy power spectrum, P(k, µ), in redshift space, measured from the Dark Sky simulations, using galaxy catalogs constructed with both halo occupation distribution and subhalo abundance matching methods, chosen to represent an intermediate redshift sample of luminous red galaxies. We find that the information content in individual µ (cosine of the angle to the line of sight) bins is substantially richer then multipole moments, and show that this can be used to isolate the impact of nonlinear growth and redshift space distortion (RSD) effects. Using the µ < 0.2 simulation data, which we show is not impacted by RSD effects, we can successfully measure the nonlinear bias to an accuracy of ~ 5% at k < 0.6hMpc-1 . This use of individual µ bins to extract the nonlinear bias successfully removes a large parameter degeneracy when constraining the linear growth rate of structure. We carry out a joint parameter estimation, using the low µ simulation data to constrain the nonlinear bias, and µ > 0.2 to constrain the growth rate and show that f can be constrained to ~ 26(22)% to a kmax < 0.4(0.6)hMpc-1 from clustering alone using a simple dispersion model, for a range of galaxy models. Our analysis of individual µ bins also reveals interesting physical effects which arise simply from different methods of populating halos with galaxies. We also find a prominent turnaround scale, at which RSD damping effects are greater then the nonlinear growth, which differs not only for each µ bin but also for each galaxy model. These features may provide unique signatures which could be used to shed light on the galaxy–dark matter connection. Furthermore, the idea of separating nonlinear growth and RSD effects making use of the full information in the 2D galaxy power spectrum yields significant improvements in constraining cosmological parameters and may be a promising probe of galaxy formation models.

  14. A new 3D parallel high resolution electromagnetic nonlinear inversion based on new global magnetic integral and local differential decomposition (GILD)

    SciTech Connect

    Xie, G.; Li, J.

    1997-05-01

    A new 3D electromagnetic modeling and nonlinear inversion algorithm is presented based on global integral and local differential equations decomposition (GILD). The GILD parallel nonlinear inversion algorithm consists of five parts: (1) the domain is decomposed into subdomain SI and subdomain SII; (2) a new global magnetic integral equation in SI and the local magnetic differential equations IN SII will be used together to obtain the magnetic field in the modeling step; (3) the new global magnetic integral Jacobian equation in SI and the local magnetic differential Jacobian equations in SII will be used together to update the electric conductivity and permittivity from the magnetic field data in the inversion step; (4) the subdomain SII can naturally and uniformly be decomposed into 2{sup n} smaller sub-cubic-domains; the sparse matrix in each sub-cubic-domain can be eliminated separately, in parallel; (5) a new parallel multiple hierarchy substructure algorithm will be used to solve the smaller full matrices in SI, in parallel. The applications of the new 3D parallel GILD EM modeling and nonlinear inversion algorithm and software are: (1) to create high resolution controlled-source electric conductivity and permittivity imaging for interpreting electromagnetic field data acquired from cross hole, surface to borehole, surface to surface, single hole, and multiple holes; (2) to create the magnetotelluric high resolution imaging from the surface impedance and field data. The new GILD parallel nonlinear inversion will be a 3D/2.5D powerful imaging tool for the oil geophysical exploration and environmental remediation and monitoring.

  15. Analysis of the Nonlinear Trends and Non-Stationary Oscillations of Regional Precipitation in Xinjiang, Northwestern China, Using Ensemble Empirical Mode Decomposition.

    PubMed

    Guo, Bin; Chen, Zhongsheng; Guo, Jinyun; Liu, Feng; Chen, Chuanfa; Liu, Kangli

    2016-03-01

    Changes in precipitation could have crucial influences on the regional water resources in arid regions such as Xinjiang. It is necessary to understand the intrinsic multi-scale variations of precipitation in different parts of Xinjiang in the context of climate change. In this study, based on precipitation data from 53 meteorological stations in Xinjiang during 1960-2012, we investigated the intrinsic multi-scale characteristics of precipitation variability using an adaptive method named ensemble empirical mode decomposition (EEMD). Obvious non-linear upward trends in precipitation were found in the north, south, east and the entire Xinjiang. Changes in precipitation in Xinjiang exhibited significant inter-annual scale (quasi-2 and quasi-6 years) and inter-decadal scale (quasi-12 and quasi-23 years). Moreover, the 2-3-year quasi-periodic fluctuation was dominant in regional precipitation and the inter-annual variation had a considerable effect on the regional-scale precipitation variation in Xinjiang. We also found that there were distinctive spatial differences in variation trends and turning points of precipitation in Xinjiang. The results of this study indicated that compared to traditional decomposition methods, the EEMD method, without using any a priori determined basis functions, could effectively extract the reliable multi-scale fluctuations and reveal the intrinsic oscillation properties of climate elements. PMID:27007388

  16. Analysis of the Nonlinear Trends and Non-Stationary Oscillations of Regional Precipitation in Xinjiang, Northwestern China, Using Ensemble Empirical Mode Decomposition

    PubMed Central

    Guo, Bin; Chen, Zhongsheng; Guo, Jinyun; Liu, Feng; Chen, Chuanfa; Liu, Kangli

    2016-01-01

    Changes in precipitation could have crucial influences on the regional water resources in arid regions such as Xinjiang. It is necessary to understand the intrinsic multi-scale variations of precipitation in different parts of Xinjiang in the context of climate change. In this study, based on precipitation data from 53 meteorological stations in Xinjiang during 1960–2012, we investigated the intrinsic multi-scale characteristics of precipitation variability using an adaptive method named ensemble empirical mode decomposition (EEMD). Obvious non-linear upward trends in precipitation were found in the north, south, east and the entire Xinjiang. Changes in precipitation in Xinjiang exhibited significant inter-annual scale (quasi-2 and quasi-6 years) and inter-decadal scale (quasi-12 and quasi-23 years). Moreover, the 2–3-year quasi-periodic fluctuation was dominant in regional precipitation and the inter-annual variation had a considerable effect on the regional-scale precipitation variation in Xinjiang. We also found that there were distinctive spatial differences in variation trends and turning points of precipitation in Xinjiang. The results of this study indicated that compared to traditional decomposition methods, the EEMD method, without using any a priori determined basis functions, could effectively extract the reliable multi-scale fluctuations and reveal the intrinsic oscillation properties of climate elements. PMID:27007388

  17. Nonlinear spacing and frequency effects of an oscillating cylinder in the wake of a stationary cylinder

    NASA Astrophysics Data System (ADS)

    Yang, Xiaofan; Zheng, Zhongquan Charlie

    2010-04-01

    Nonlinear responses to a transversely oscillating cylinder in the wake of a stationary upstream cylinder are studied theoretically by using an immersed-boundary method at Re=100. Response states are investigated in the three flow regimes for a tandem-cylinder system: the "vortex suppression" regime, the critical spacing regime, and the "vortex formation" regime. When the downstream cylinder is forced to oscillate at a fixed frequency and amplitude, the response state of flow around the two cylinders varies with different spacing between the two cylinders, while in the same flow regime, the response state can change with the oscillating frequency and amplitude of the downstream cylinder. Based on velocity phase portraits, each of the nonlinear response states can be categorized into one of the three states in the order of increasing chaotic levels: lock-in, transitional, or quasiperiodic. These states can also be correlated with velocity spectral behaviors. The discussions are conducted using near-wake velocity phase portraits, spectral analyses, and related vorticity fields. A general trend in the bifurcation diagrams of frequency spacing shows the smaller the spacing, frequency, or amplitude, the less chaotic the response state of the system and more likely the downstream and upstream wakes are in the same response state. The system is not locked-in in any case when the spacing between the cylinders is larger than the critical spacing. The near-wake velocity spectral behaviors correspond to the nonlinear response states, with narrow-banded peaks shown at the oscillation frequency and its harmonics in the lock-in cases. High frequency harmonic peaks, caused by interactions between the upstream wake and the downstream oscillating cylinder, are reduced in the near-wake velocity spectra of the upstream cylinder when the spacing increases.

  18. A stable elemental decomposition for dynamic process optimization

    NASA Astrophysics Data System (ADS)

    Cervantes, Arturo M.; Biegler, Lorenz T.

    2000-08-01

    In Cervantes and Biegler (A.I.Ch.E.J. 44 (1998) 1038), we presented a simultaneous nonlinear programming problem (NLP) formulation for the solution of DAE optimization problems. Here, by applying collocation on finite elements, the DAE system is transformed into a nonlinear system. The resulting optimization problem, in which the element placement is fixed, is solved using a reduced space successive quadratic programming (rSQP) algorithm. The space is partitioned into range and null spaces. This partitioning is performed by choosing a pivot sequence for an LU factorization with partial pivoting which allows us to detect unstable modes in the DAE system. The system is stabilized without imposing new boundary conditions. The decomposition of the range space can be performed in a single step by exploiting the overall sparsity of the collocation matrix but not its almost block diagonal structure. In order to solve larger problems a new decomposition approach and a new method for constructing the quadratic programming (QP) subproblem are presented in this work. The decomposition of the collocation matrix is now performed element by element, thus reducing the storage requirements and the computational effort. Under this scheme, the unstable modes are considered in each element and a range-space move is constructed sequentially based on decomposition in each element. This new decomposition improves the efficiency of our previous approach and at the same time preserves its stability. The performance of the algorithm is tested on several examples. Finally, some future directions for research are discussed.

  19. Optical authentication based on moiré effect of nonlinear gratings in phase space

    NASA Astrophysics Data System (ADS)

    Liao, Meihua; He, Wenqi; Wu, Jiachen; Lu, Dajiang; Liu, Xiaoli; Peng, Xiang

    2015-12-01

    An optical authentication scheme based on the moiré effect of nonlinear gratings in phase space is proposed. According to the phase function relationship of the moiré effect in phase space, an arbitrary authentication image can be encoded into two nonlinear gratings which serve as the authentication lock (AL) and the authentication key (AK). The AL is stored in the authentication system while the AK is assigned to the authorized user. The authentication procedure can be performed using an optoelectronic approach, while the design process is accomplished by a digital approach. Furthermore, this optical authentication scheme can be extended for multiple users with different security levels. The proposed scheme can not only verify the legality of a user identity, but can also discriminate and control the security levels of legal users. Theoretical analysis and simulation experiments are provided to verify the feasibility and effectiveness of the proposed scheme.

  20. Some exact solutions of a system of nonlinear Schroedinger equations in three-dimensional space

    SciTech Connect

    Moskalyuk, S.S.

    1988-02-01

    Interactions that break the symmetry of systems of nonrelativistic Schroedinger equations but preserve their symmetry with respect to one-parameter subgroups of the Schroedinger group are described. Ansatzes for invariant solutions and the corresponding systems of reduced equations in invariant variables for Galileo-invariant Schroedinger equations are found. Exact solutions for the system of nonlinear Schroedinger equations in three-dimensional space for the generalized Hubbard model are obtained.

  1. Linear and nonlinear modulus surfaces in stress space, from stress-strain measurements on Berea sandstone

    NASA Astrophysics Data System (ADS)

    Boudjema, M.; Santos, I. B.; McCall, K. R.; Guyer, R. A.; Boitnott, G. N.

    The elastic response of many rocks to quasistatic stress changes is highly nonlinear and hysteretic, displaying discrete memory. Rocks also display unusual nonlinear response to dynamic stress changes. A model to describe the elastic behavior of rocks and other consolidated materials is called the Preisach-Mayergoyz (PM) space model. In contrast to the traditional analytic approach to stress-strain, the PM space picture establishes a relationship between the quasistatic data and a number density of hysteretic mesoscopic elastic elements in the rock. The number density allows us to make quantitative predictions of dynamic elastic properties. Using the PM space model, we analyze a complex suite of quasistatic stress-strain data taken on Berea sandstone. We predict a dynamic bulk modulus and a dynamic shear modulus surface as a function of mean stress and shear stress. Our predictions for the dynamic moduli compare favorably to moduli derived from time of flight measurements. We derive a set of nonlinear elastic constants and a set of constants that describe the hysteretic behavior of the sandstone.

  2. Evidence for a non-linear relationship between fracture spacing and layer thickness

    NASA Astrophysics Data System (ADS)

    Mandal, Nibir; Deb, Sanjitendra Krishna; Khan, Debdarpan

    1994-09-01

    In experiments, extension fractures were generated in rigid layers of Plaster of Paris resting on a viscous substratum (pitch). The experimental results predict a non-linear relationship between the spacing of fractures in uniform brittle layers and layer thickness for fractures generated by the tractional force of embeddins weak rocks. We derive an equation which relates the critical fracture spacing ( λc) with layer thickness ( b), tensile strength of the layer material/viscosity of the embedding medium ratio ( τo/ η) and bulk strain rate ( ɛ˙b). The equation shows that the spacing increases as a function of the square root of the layer thickness. The theory also predicts that the fracture spacing depends on the strain rate when the embedding weak medium is viscous.

  3. Contamination and Radiation Effects on Nonlinear Crystals for Space Laser Systems

    NASA Technical Reports Server (NTRS)

    Abdeldayem, Hossain A.; Dowdye, Edward; Jamison, Tracee; Canham, John; Jaeger, Todd

    2005-01-01

    Space Lasers are vital tools for NASA s space missions and military applications. Although, lasers are highly reliable on the ground, several past space laser missions proved to be short-lived and unreliable. In this communication, we are shedding more light on the contamination and radiation issues, which are the most common causes for optical damages and laser failures in space. At first, we will present results based on the study of liquids and subsequently correlate these results to the particulates of the laser system environment. We present a model explaining how the laser beam traps contaminants against the optical surfaces and cause optical damages and the role of gravity in the process. We also report the results of the second harmonic generation efficiency for nonlinear optical crystals irradiated with high-energy beams of protons. In addition, we are proposing to employ the technique of adsorption to minimize the presence of adsorbing molecules present in the laser compartment.

  4. Disentangling redshift-space distortions and non-linear bias using the 2D power spectrum

    NASA Astrophysics Data System (ADS)

    Jennings, Elise; Wechsler, Risa H.; Skillman, Samuel W.; Warren, Michael S.

    2016-03-01

    We present the 2D redshift-space galaxy power spectrum, P(k, μ), measured from the Dark Sky simulations, using catalogues constructed with halo occupation distribution and subhalo abundance matching methods, chosen to represent an intermediate redshift sample of luminous red galaxies. We find that the information content in individual μ (cosine of the angle to the line of sight) bins is substantially richer then multipole moments, and show that this can be used to isolate the impact of non-linear growth and redshift-space distortion (RSD) effects. Using the μ < 0.2 simulation data, which is not impacted by RSD, we can successfully measure the non-linear bias to ˜5 per cent at k < 0.6 h Mpc-1. Using the low μ simulation data to constrain the non-linear bias, and μ ≥ 0.2 to constrain the growth rate, we show that f can be constrained to ˜26(22) per cent to a kmax < 0.4(0.6) h Mpc-1 from clustering alone using a dispersion model, for a range of galaxy models. Our analysis of individual μ bins reveals interesting physical effects which arise from different methods of populating haloes with galaxies. We find a prominent turnaround scale, at which RSD damping effects are greater than the non-linear growth, which differs for each galaxy model. The idea of separating non-linear growth and RSD effects making use of the full information in the 2D galaxy power spectrum yields significant improvements in constraining cosmological parameters and may be a promising probe of galaxy formation models.

  5. Nonlinear alternatives of Schauder and Krasnosel'skij types with applications to Hammerstein integral equations in L spaces

    NASA Astrophysics Data System (ADS)

    Djebali, Smaïl; Sahnoun, Zahira

    This paper is devoted to establishing new variants of some nonlinear alternatives of Leray-Schauder and Krasnosel'skij type involving the weak topology of Banach spaces. The De Blasi measure of weak noncompactness is used. An application to solving a nonlinear Hammerstein integral equation in L spaces is given. Our results complement recent ones in [K. Latrach, M.A. Taoudi, A. Zeghal, Some fixed point theorems of the Schauder and the Krasnosel'skij type and application to nonlinear transport equations, J. Differential Equations 221 (2006) 256-2710] and [K. Latrach, M.A. Taoudi, Existence results for a generalized nonlinear Hammerstein equation on L spaces, Nonlinear Anal. 66 (2007) 2325-2333].

  6. Disentangling Redshift-Space Distortions and Nonlinear Bias using the 2D Power Spectrum

    DOE PAGESBeta

    Jennings, Elise; Wechsler, Risa H.

    2015-08-07

    We present the nonlinear 2D galaxy power spectrum, P(k, µ), in redshift space, measured from the Dark Sky simulations, using galaxy catalogs constructed with both halo occupation distribution and subhalo abundance matching methods, chosen to represent an intermediate redshift sample of luminous red galaxies. We find that the information content in individual µ (cosine of the angle to the line of sight) bins is substantially richer then multipole moments, and show that this can be used to isolate the impact of nonlinear growth and redshift space distortion (RSD) effects. Using the µ < 0.2 simulation data, which we show ismore » not impacted by RSD effects, we can successfully measure the nonlinear bias to an accuracy of ~ 5% at k < 0.6hMpc-1 . This use of individual µ bins to extract the nonlinear bias successfully removes a large parameter degeneracy when constraining the linear growth rate of structure. We carry out a joint parameter estimation, using the low µ simulation data to constrain the nonlinear bias, and µ > 0.2 to constrain the growth rate and show that f can be constrained to ~ 26(22)% to a kmax < 0.4(0.6)hMpc-1 from clustering alone using a simple dispersion model, for a range of galaxy models. Our analysis of individual µ bins also reveals interesting physical effects which arise simply from different methods of populating halos with galaxies. We also find a prominent turnaround scale, at which RSD damping effects are greater then the nonlinear growth, which differs not only for each µ bin but also for each galaxy model. These features may provide unique signatures which could be used to shed light on the galaxy–dark matter connection. Furthermore, the idea of separating nonlinear growth and RSD effects making use of the full information in the 2D galaxy power spectrum yields significant improvements in constraining cosmological parameters and may be a promising probe of galaxy formation models.« less

  7. Phase-space description of plasma waves: Linear and nonlinear theory

    NASA Astrophysics Data System (ADS)

    Biro, Thomas

    1992-11-01

    A (r,k) phase description of waves in plasmas is developed by introducing Gaussian window functions to separate short scale oscillations from long scale modulations of the wave fields and variations in the plasma parameters. To obtain a wave equation that unambiguously separates conservative dynamics from dissipation also in an inhomogeneous and time varying background plasma, the proper form of the current response function, is discussed. On the analogy of the particle distribution function f(v,r,t), a wave density N(k,r,t) is introduced on phase space. This function is proven to satisfy a simple continuity equation. Dissipation is also included, and this allows the damping or growth of wave density along rays to be described. Problems involving geometric optics of continuous media often appear simpler when viewed in phase space, since the flow of N in phase space is incompressible. Within the phase space representation, a very general formula for the second order nonlinear current is obtained in terms of the vector potential. This formula is a convenient starting point for studies of coherent as well as turbulent nonlinear processes. Kinetic equations for weakly inhomogeneous and turbulent plasmas are derived, including the effects of inhomogeneous turbulence, wave convection and refraction.

  8. Study of outgassing and decomposition of space shuttle heat protection tiles, fillers and adhesive

    NASA Technical Reports Server (NTRS)

    Proctor, B. L.; Hoffman, J. H.

    1982-01-01

    The purpose of this project was to determine the chemicals desorbing from the space shuttle heat protection tiles. The original protocol for this project involved direct insertion probe mass spectrometry (DIPMS) analysis of the outgassing products from the tiles. However, this method proved unsatisfactory due to the large number of compounds desorbing from the tiles. A purge and trap technique was then employed to collect and separate the chemicals desorbing from the tiles. The maximum temperature in this analysis was 180 C which is the gas chromatograph fused silica capillary column's temperature limit. The desorption was also carried out at atmospheric pressure with helium as the purge gas. A description of the modified protocol is given. All compounds are tentatively identified.

  9. The non-linear redshift-space power spectrum: Omega from redshift surveys

    NASA Astrophysics Data System (ADS)

    Fisher, Karl B.; Nusser, Adi

    1996-03-01

    We examine the anisotropies in the power spectrum by the mapping of real space to redshift space. Using the Zel'dovich approximation, we obtain an analytic expression for the non-linear redshift-space power spectrum in the distant observer limit. For a given unbiased galaxy distribution in redshift space, the anisotropies in the power spectrum depend on the parameter f(Omega)~=Omega^0.6, where Omega is the density parameter. We quantify these anisotropies by the ratio, R, of the quadrupole and monopole angular moments of the power spectrum. In contrast to linear theory, the Zel'dovich approximation predicts a decline in R with decreasing scale. This departure from linear theory is due to non-linear dynamics and is not a result of incoherent random velocities. The rate of decline depends strongly on Omega and the initial power spectrum. However, we find a scaling relation between the quantity R/R_lin (where R_lin is the linear theory value of R) and the dimensionless variable k/k_nl, where k_nl is a wavenumber determined by the scale of non-linear structures. The scaling is weakly dependent on the initial power spectrum and is in good agreement with a large N-body simulation. This universal scaling relation greatly extends the scales over which redshift distortions can be used as a probe of Omega. The scaling relation is in agreement with the observed quadrupole-to-monopole ratio from the 1.2-Jy IRAS survey, with a best estimate of Omega^0.6/b_lin=0.6+/-0.2 where b_lin is the linear bias parameter.

  10. Stable double spacing multiwavelength Brillouin-Erbium doped fiber laser based on highly nonlinear fiber

    NASA Astrophysics Data System (ADS)

    Ahmad, B. A.; Al-Alimi, A. W.; Abas, A. F.; Harun, S. W.; Mahdi, M. A.

    2012-05-01

    A double frequency spaced multiwavelength Brillouin-Erbium doped fiber laser (BEDFL) with figure-of-eight cavity have been successfully developed and tested. Double frequency spacing is achieved by using a piece of 2 km of highly nonlinear fiber (HNLF) as a gain medium. Figure-of-eight configuration removes the odd order Stokes signals via a four-port circulator. Fifteen Stokes channels are simultaneously generated with a spacing of 0.154 nm that is around 20 GHz, when the Brillouin pump and 980 nm pump powers are fixed at the optimized values of 6 dBm and 40 mW, respectively. Fourteen anti stoke channels are also obtained, which are generated through four wave mixing (FWM) process in the laser cavity. The output is smooth triangular comb. The BEFL can also be tuned from 1526.5 to 1567.5 nm.

  11. Picosecond dynamics of free carrier populations, space- charge fields, and photorefractive nonlinearities in zincblende semiconductors

    NASA Astrophysics Data System (ADS)

    Stark, Thomas Spencer

    Generally, nonlinear optics studies investigate optically-induced changes in refraction or absorption, and their application to spectroscopy or device fabrication. The photorefractive effect is a nonlinear optical effect that occurs in solids, where transport of an optically-induced free-carrier population results in an internal space-charge field, which produces an index change via the linear electrooptic effect. The photorefractive effect has been widely studied for a variety of materials and device applications, mainly because it allows large index changes to be generated with laser beams having only a few milliwatts of average power. Compound semiconductors are important photorefractive materials because they offer a near-infrared optical response, and because their carrier transport properties allow the index change to be generated quickly and efficiently. While many researchers have attempted to measure the fundamental temporal dynamics of the photorefractive effect in semiconductors using continuous-wave, nanosecond- and picosecond-pulsed laser beams, these investigations have been unsuccessful. However, studies with this goal are of clear relevance because they provide information about the fundamental physical processes that produce this effect, as well as the material's speed and efficiency limitations for device applications. In this dissertation, for the first time, we time-resolve the temporal dynamics of the photorefractive nonlinearities in two zincblende semiconductors, semi- insulating GaAs and undoped CdTe. While CdTe offers a lattice-match to the infrared material HgxCd1-xTe, semi-insulating GaAs has been widely used in optoelectronic and high- speed electronic applications. We use a novel transient- grating experimental method that allows picosecond temporal resolution and high sensitivity. Our results provide a clear and detailed picture of the picosecond photorefractive response of both materials, showing nonlinearities due to hot

  12. Efficient time-domain simulation of nonlinear, state-space, transmission-line models of the cochlea (L).

    PubMed

    Pan, Shuokai; Elliott, Stephen J; Teal, Paul D; Lineton, Ben

    2015-06-01

    Nonlinear models of the cochlea are best implemented in the time domain, but their computational demands usually limit the duration of the simulations that can reasonably be performed. This letter presents a modified state space method and its application to an example nonlinear one-dimensional transmission-line cochlear model. The sparsity pattern of the individual matrices for this alternative formulation allows the use of significantly faster numerical algorithms. Combined with a more efficient implementation of the saturating nonlinearity, the computational speed of this modified state space method is more than 40 times faster than that of the original formulation. PMID:26093443

  13. Analyzing the Variability and Teleconnection of Seasonal Precipitation over Western U.S. Using Empirical Orthogonal Function and Nonlinear Mode Decomposition

    NASA Astrophysics Data System (ADS)

    Liu, H.; Gao, X.; Sorooshian, S.

    2014-12-01

    Seasonal variations, decadal trends and spatial patterns of precipitation are crucial to water managements in western United States. We use the Empirical Orthogonal Function-Principal Component (EOF-PC) and Nonlinear Mode Decomposition (NMD) methods to analyze gridded monthly precipitation observation data in western US from 1971 to 2013.To differentiate trends, noises and physically-meaningful periodic patterns from the time series, we choose the wavelet-based NMD method. The results of our analysis show the leading EOFs and PC time series are capable of reconstructing the dynamic field of regional precipitation in high accuracy. The NMD method detects several amplitude-modulated periodic curves with complex waveforms. We found the first PC's varying amplitude is correlated to Pacific Decadal Oscillation index (R=0.41), and the fourth PC's varying amplitude is correlated to Nino3.4 index (R=0.57). These results suggest that the complex spatiotemporal variability of precipitation can be quantified with the dominant spatial patterns, periodicities, trends, and noises. And the EOF-NMD method can be used as a data-mining tool to investigate physically-meaningful information and teleconnections from a dataset of a dynamic field.

  14. Nonlinear spacing and frequency effects of an oscillating cylinder in the wake of a stationary cylinder

    NASA Astrophysics Data System (ADS)

    Zheng, Z.; Yang, Xiaofan

    2008-11-01

    Nonlinear responses to a transversely oscillating cylinder in the wake of a stationary upstream cylinder are studied theoretically by using an immersed-boundary method. It is found that flow around the two cylinders varies with different spacing between the two cylinders and the oscillation frequency of the downstream cylinder. As known in a stationary tandem-cylinder system, there exist the ``vortex suppression regime'' (VS) and the ``vortex formation regime'' (VF). These two regimes are divided by a critical spacing. When the downstream cylinder is forced to oscillate at a fixed amplitude but different frequency, different flow patterns appear in each of the regime. On the other hand, at the same oscillating frequency but different spacing, the response state (lock-in, transient or non-lock-in) changes. While each state has periodic or quasi-periodic behaviors, nonlinear responses appear. All of the analyses are based on vorticity contours, time histories of the velocities in the near wake regions, spectral analyses, and related phase portraits.

  15. Steady-state analysis of a nonlinear rotor-housing system. [Space Shuttle Main Engine

    NASA Technical Reports Server (NTRS)

    Noah, S. T.; Kim, Y. B.

    1990-01-01

    The periodic steady state response of a high pressure oxygen turbopump (HBOTP) of a Space Shuttle main engine (SSME), involving a clearance between the bearing and housing carrier, is sought. A harmonic balance method utilizig Fast Fourier Transform (FFT) algorithm is developed for the analysis. An impedance method is used to reduce the number of degrees of freedom to the displacements at the bearing clearance. Harmonic and subharmonic responses to imbalance for various system parameters are studied. The results show that the computational technique developed in this study is an effective and flexible method for determining the stable and unstable periodic response of complex rotor-housing systems with clearance type nonlinearity.

  16. 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.

  17. Inferring gene regulatory networks via nonlinear state-space models and exploiting sparsity.

    PubMed

    Noor, Amina; Serpedin, Erchin; Nounou, Mohamed; Nounou, Hazem N

    2012-01-01

    This paper considers the problem of learning the structure of gene regulatory networks from gene expression time series data. A more realistic scenario when the state space model representing a gene network evolves nonlinearly is considered while a linear model is assumed for the microarray data. To capture the nonlinearity, a particle filter-based state estimation algorithm is considered instead of the contemporary linear approximation-based approaches. The parameters characterizing the regulatory relations among various genes are estimated online using a Kalman filter. Since a particular gene interacts with a few other genes only, the parameter vector is expected to be sparse. The state estimates delivered by the particle filter and the observed microarray data are then subjected to a LASSO-based least squares regression operation which yields a parsimonious and efficient description of the regulatory network by setting the irrelevant coefficients to zero. The performance of the aforementioned algorithm is compared with the extended Kalman filter (EKF) and Unscented Kalman Filter (UKF) employing the Mean Square Error (MSE) as the fidelity criterion in recovering the parameters of gene regulatory networks from synthetic data and real biological data. Extensive computer simulations illustrate that the proposed particle filter-based network inference algorithm outperforms EKF and UKF, and therefore, it can serve as a natural framework for modeling gene regulatory networks with nonlinear and sparse structure. PMID:22350207

  18. Wave operators to a quadratic nonlinear Klein-Gordon equation in two space dimensions revisited

    NASA Astrophysics Data System (ADS)

    Hayashi, Nakao; Naumkin, Pavel I.; Tonegawa, Satoshi

    2012-08-01

    We continue to study the existence of the wave operators for the nonlinear Klein-Gordon equation with quadratic nonlinearity in two space dimensions {(partialt2-Δ+m2) u=λ u2,( t,x) in{R}×{R}2}. We prove that if u1+in{H}^{3/2+3γ,1}( {R}2),{ }u2+in{H}^{1/2+3γ,1}( {R} 2), where {γin( 0,1/4)} and the norm {Vert u1+Vert_{{H}1^{3/2+γ}}+Vert u2+Vert_{{H}1^{1/2+γ}}≤ρ,} then there exist ρ > 0 and T > 1 such that the nonlinear Klein-Gordon equation has a unique global solution {uin{C}( [ T,infty) ;{H}^{1/2}( {R}2) ) } satisfying the asymptotics Vert u( t) -u0 ( t) Vert _{{H}^{1/2}} ≤ Ct^{-1/2-γ} for all t > T, where u 0 denotes the solution of the free Klein-Gordon equation.

  19. Interferometric and nonlinear-optical spectral-imaging techniques for outer space and live cells

    NASA Astrophysics Data System (ADS)

    Itoh, Kazuyoshi

    2015-12-01

    Multidimensional signals such as the spectral images allow us to have deeper insights into the natures of objects. In this paper the spectral imaging techniques that are based on optical interferometry and nonlinear optics are presented. The interferometric imaging technique is based on the unified theory of Van Cittert-Zernike and Wiener-Khintchine theorems and allows us to retrieve a spectral image of an object in the far zone from the 3D spatial coherence function. The retrieval principle is explained using a very simple object. The promising applications to space interferometers for astronomy that are currently in progress will also be briefly touched on. An interesting extension of interferometric spectral imaging is a 3D and spectral imaging technique that records 4D information of objects where the 3D and spectral information is retrieved from the cross-spectral density function of optical field. The 3D imaging is realized via the numerical inverse propagation of the cross-spectral density. A few techniques suggested recently are introduced. The nonlinear optical technique that utilizes stimulated Raman scattering (SRS) for spectral imaging of biomedical targets is presented lastly. The strong signals of SRS permit us to get vibrational information of molecules in the live cell or tissue in real time. The vibrational information of unstained or unlabeled molecules is crucial especially for medical applications. The 3D information due to the optical nonlinearity is also the attractive feature of SRS spectral microscopy.

  20. New Space Weather and Nonlinear Waves and Processes Prize announced for 2013

    NASA Astrophysics Data System (ADS)

    Thompson, Victoria

    2012-01-01

    At the 2011 Fall Meeting in San Francisco, Calif., AGU announced the creation of a new award: the Space Weather and Nonlinear Waves and Processes Prize. The prize, which is being made possible by a generous contribution from longtime AGU members and NASA Jet Propulsion Laboratory (JPL), California Institute of Technology, scientists Bruce Tsurutani and Olga Verkhoglyadova, will recognize an AGU member scientist and will come with a $10,000 award. Tsurutani has served as a researcher with JPL since 1972 and is currently a senior research scientist. He was also the president of AGU's Space Physics and Aeronomy section from 1990 to 1992 and is a recipient of AGU's John Adam Fleming Medal, given “for original research and technical leadership in geomagnetism, atmospheric electricity, aeronomy, space physics, and related sciences.” Verkhoglyadova served as a professor of space physics in the Department of Astrophysics and Space Physics at Taras Shevchenko National University of Kyiv, in the Ukraine, prior to coming to the United States. Their leadership and dedication to AGU and to their field are apparent in their passion for this prize.

  1. Nonlinear Phase Mixing and Phase-Space Cascade of Entropy in Gyrokinetic Plasma Turbulence

    SciTech Connect

    Tatsuno, T.; Dorland, W.; Plunk, G. G.; Schekochihin, A. A.; Barnes, M.

    2009-07-03

    Electrostatic turbulence in weakly collisional, magnetized plasma can be interpreted as a cascade of entropy in phase space, which is proposed as a universal mechanism for dissipation of energy in magnetized plasma turbulence. When the nonlinear decorrelation time at the scale of the thermal Larmor radius is shorter than the collision time, a broad spectrum of fluctuations at sub-Larmor scales is numerically found in velocity and position space, with theoretically predicted scalings. The results are important because they identify what is probably a universal Kolmogorov-like regime for kinetic turbulence; and because any physical process that produces fluctuations of the gyrophase-independent part of the distribution function may, via the entropy cascade, result in turbulent heating at a rate that increases with the fluctuation amplitude, but is independent of the collision frequency.

  2. Efficient solution of two-sided nonlinear space-fractional diffusion equations using fast Poisson preconditioners

    NASA Astrophysics Data System (ADS)

    Moroney, Timothy; Yang, Qianqian

    2013-08-01

    We develop a fast Poisson preconditioner for the efficient numerical solution of a class of two-sided nonlinear space-fractional diffusion equations in one and two dimensions using the method of lines. Using the shifted Grünwald finite difference formulas to approximate the two-sided (i.e. the left and right Riemann-Liouville) fractional derivatives, the resulting semi-discrete nonlinear systems have dense Jacobian matrices owing to the non-local property of fractional derivatives. We employ a modern initial value problem solver utilising backward differentiation formulas and Jacobian-free Newton-Krylov methods to solve these systems. For efficient performance of the Jacobian-free Newton-Krylov method it is essential to apply an effective preconditioner to accelerate the convergence of the linear iterative solver. The key contribution of our work is to generalise the fast Poisson preconditioner, widely used for integer-order diffusion equations, so that it applies to the two-sided space-fractional diffusion equation. A number of numerical experiments are presented to demonstrate the effectiveness of the preconditioner and the overall solution strategy.

  3. Nonlinear model of space manipulator joint considering time-variant stiffness and backlash

    NASA Astrophysics Data System (ADS)

    Yang, Tianfu; Yan, Shaoze; Han, Zengyao

    2015-04-01

    Modeling of space manipulator joints has been studied for years but accurate positioning control is still unsatisfactory. One of the primary reasons is that, in the past researches, effects of the high-ratio reducers in the joints have usually been neglected. In this paper, a nonlinear dynamic model of the manipulator joint with planetary gear train transmission is developed by considering time-variant joint stiffness, backlash and reduction ratio. Based on the gear parameters and meshing phase relationship, the stiffness of the joint model is presented, in which the time-variant stiffness of 2K-H planetary gear train and the backlash are taken into consideration. The backlash effect is modeled as an alternate engagement mechanism, and the transmitted torque is defined as a dead zone function. This model is simulated on a two-link space manipulator system. The results show that the time-variant stiffness effect can be simplified as a constant value in most cases when other shafting parts are flexible, while if the total stiffness is approximate to the nonlinear stiffness, the positioning accuracy is reduced if neglecting the time-variant part. On the other hand, the backlash is the main source of positioning error and impact. Minimizing backlash is the most effective way to improve positioning accuracy and avoid the impact in the gearing system.

  4. Non-Linear Cosmological Power Spectra in Real and Redshift Space

    NASA Technical Reports Server (NTRS)

    Taylor, A. N.; Hamilton, A. J. S.

    1996-01-01

    We present an expression for the non-linear evolution of the cosmological power spectrum based on Lagrangian trajectories. This is simplified using the Zel'dovich approximation to trace particle displacements, assuming Gaussian initial conditions. The model is found to exhibit the transfer of power from large to small scales expected in self-gravitating fields. Some exact solutions are found for power-law initial spectra. We have extended this analysis into red-shift space and found a solution for the non-linear, anisotropic redshift-space power spectrum in the limit of plane-parallel redshift distortions. The quadrupole-to-monopole ratio is calculated for the case of power-law initial spectra. We find that the shape of this ratio depends on the shape of the initial spectrum, but when scaled to linear theory depends only weakly on the redshift-space distortion parameter, beta. The point of zero-crossing of the quadrupole, kappa(sub o), is found to obey a simple scaling relation and we calculate this scale in the Zel'dovich approximation. This model is found to be in good agreement with a series of N-body simulations on scales down to the zero-crossing of the quadrupole, although the wavenumber at zero-crossing is underestimated. These results are applied to the quadrupole-to-monopole ratio found in the merged QDOT plus 1.2-Jy-IRAS redshift survey. Using a likelihood technique we have estimated that the distortion parameter is constrained to be beta greater than 0.5 at the 95 percent level. Our results are fairly insensitive to the local primordial spectral slope, but the likelihood analysis suggests n = -2 un the translinear regime. The zero-crossing scale of the quadrupole is k(sub 0) = 0.5 +/- 0.1 h Mpc(exp -1) and from this we infer that the amplitude of clustering is sigma(sub 8) = 0.7 +/- 0.05. We suggest that the success of this model is due to non-linear redshift-space effects arising from infall on to caustic and is not dominated by virialized cluster cores

  5. Dynamic analysis of space-related linear and non-linear structures

    SciTech Connect

    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.

  6. A non-linear model predictive controller with obstacle avoidance for a space robot

    NASA Astrophysics Data System (ADS)

    Wang, Mingming; Luo, Jianjun; Walter, Ulrich

    2016-04-01

    This study investigates the use of the non-linear model predictive control (NMPC) strategy for a kinematically redundant space robot to approach an un-cooperative target in complex space environment. Collision avoidance, traditionally treated as a high level planning problem, can be effectively translated into control constraints as part of the NMPC. The objective of this paper is to evaluate the performance of the predictive controller in a constrained workspace and to investigate the feasibility of imposing additional constraints into the NMPC. In this paper, we reformulated the issue of the space robot motion control by using NMPC with predefined objectives under input, output and obstacle constraints over a receding horizon. An on-line quadratic programming (QP) procedure is employed to obtain the constrained optimal control decisions in real-time. This study has been implemented for a 7 degree-of-freedom (DOF) kinematically redundant manipulator mounted on a 6 DOF free-floating spacecraft via simulation studies. Real-time trajectory tracking and collision avoidance particularly demonstrate the effectiveness and potential of the proposed NMPC strategy for the space robot.

  7. 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.

  8. A decomposition theorem for the space of C{sup 1}-smooth skew products with complicated dynamics of the quotient map

    SciTech Connect

    Efremova, L S

    2013-11-30

    We use the notions of the Ω-function and functions suitable to it, to give a detailed proof of a decomposition theorem for the space of C{sup 1}-smooth skew products of interval maps whose quotient maps have complicated dynamics and satisfy the additional condition of Ω-stability with respect to the C{sup 1}-norm. In our theorem, the space of C{sup 1}-smooth skew products is decomposed into a union of four nonempty, pairwise disjoint subspaces. We give examples of maps contained in each of the four subspaces. Bibliography: 46 titles.

  9. Experiments in Nonlinear Adaptive Control of Multi-Manipulator, Free-Flying Space Robots

    NASA Technical Reports Server (NTRS)

    Chen, Vincent Wei-Kang

    1992-01-01

    Sophisticated robots can greatly enhance the role of humans in space by relieving astronauts of low level, tedious assembly and maintenance chores and allowing them to concentrate on higher level tasks. Robots and astronauts can work together efficiently, as a team; but the robot must be capable of accomplishing complex operations and yet be easy to use. Multiple cooperating manipulators are essential to dexterity and can broaden greatly the types of activities the robot can achieve; adding adaptive control can ease greatly robot usage by allowing the robot to change its own controller actions, without human intervention, in response to changes in its environment. Previous work in the Aerospace Robotics Laboratory (ARL) have shown the usefulness of a space robot with cooperating manipulators. The research presented in this dissertation extends that work by adding adaptive control. To help achieve this high level of robot sophistication, this research made several advances to the field of nonlinear adaptive control of robotic systems. A nonlinear adaptive control algorithm developed originally for control of robots, but requiring joint positions as inputs, was extended here to handle the much more general case of manipulator endpoint-position commands. A new system modelling technique, called system concatenation was developed to simplify the generation of a system model for complicated systems, such as a free-flying multiple-manipulator robot system. Finally, the task-space concept was introduced wherein the operator's inputs specify only the robot's task. The robot's subsequent autonomous performance of each task still involves, of course, endpoint positions and joint configurations as subsets. The combination of these developments resulted in a new adaptive control framework that is capable of continuously providing full adaptation capability to the complex space-robot system in all modes of operation. The new adaptive control algorithm easily handles free

  10. Variational space-time (dis)continuous Galerkin method for nonlinear free surface water waves

    NASA Astrophysics Data System (ADS)

    Gagarina, E.; Ambati, V. R.; van der Vegt, J. J. W.; Bokhove, O.

    2014-10-01

    A new variational finite element method is developed for nonlinear free surface gravity water waves using the potential flow approximation. This method also handles waves generated by a wave maker. Its formulation stems from Miles' variational principle for water waves together with a finite element discretization that is continuous in space and discontinuous in time. One novel feature of this variational finite element approach is that the free surface evolution is variationally dependent on the mesh deformation vis-à-vis the mesh deformation being geometrically dependent on free surface evolution. Another key feature is the use of a variational (dis)continuous Galerkin finite element discretization in time. Moreover, in the absence of a wave maker, it is shown to be equivalent to the second order symplectic Störmer-Verlet time stepping scheme for the free-surface degrees of freedom. These key features add to the stability of the numerical method. Finally, the resulting numerical scheme is verified against nonlinear analytical solutions with long time simulations and validated against experimental measurements of driven wave solutions in a wave basin of the Maritime Research Institute Netherlands.

  11. [Study on the nonlinear characteristics of the mixed pixel's reflectance in hyperspectral space].

    PubMed

    Zhu, Feng; Gong, Hui-Li; Sun, Tian-Lin; Zhao, Yun-Sheng

    2013-03-01

    Under the experimental condition of the 50 degree incidence zenith angle and 45 degree detection azimuth, 24 groups of reflectance spectral of the mixed pixel of lotus and water body acquired using the reflex platform and FieldSpec 3 Hi-Res portable spectrum instrument. The hyperspectral space was built based on the reflectance character. The relationship between similarity and the index of lotus area ratio was analyzed using the linear, logarithm and quadratic curve fitting, and the goodness of fitting is 63.6%, 76.2% and 82.9% respectively. According to the real relationship of the mixed pixel spectral vector and the reference spectral, the best fitting model has nonlinear characteristics. The idea that the mixed pixel may have the critical value was proposed on the base of the analysis. The research result will help understand the mixed pixel further, and provide a new direction for unmixing the mixed pixel. PMID:23705444

  12. Modal test/analysis correlation of Space Station structures using nonlinear sensitivity

    NASA Astrophysics Data System (ADS)

    Gupta, Viney K.; Newell, James F.; Berke, Laszlo; Armand, Sasan

    1992-09-01

    The modal correlation problem is formulated as a constrained optimization problem for validation of finite element models (FEM's). For large-scale structural applications, a pragmatic procedure for substructuring, model verification, and system integration is described to achieve effective modal correlations. The space station substructure FEM's are reduced using Lanczos vectors and integrated into a system FEM using Craig-Bampton component modal synthesis. The optimization code is interfaced with MSC/NASTRAN to solve the problem of modal test/analysis correlation; that is, the problem of validating FEM's for launch and on-orbit coupled loads analysis against experimentally observed frequencies and mode shapes. An iterative perturbation algorithm is derived and implemented to update nonlinear sensitivity (derivatives of eigenvalues and eigenvectors) during optimizer iterations, which reduced the number of finite element analyses.

  13. Modal test/analysis correlation of Space Station structures using nonlinear sensitivity

    NASA Astrophysics Data System (ADS)

    Gupta, Viney K.; Newell, James F.; Berke, Laszlo; Armand, Sasan

    1992-09-01

    The modal correlation problem is formulated as a constrained optimization problem for validation of finite element models (FEM's). For large-scale structural applications, a pragmatic procedure for substructuring, model verification, and system integration is described to achieve effective modal correlation. The space station substructure FEM's are reduced using Lanczos vectors and integrated into a system FEM using Craig-Bampton component modal synthesis. The optimization code is interfaced with MSC/NASTRAN to solve the problem of modal test/analysis correlation; that is, the problem of validating FEM's for launch and on-orbit coupled loads analysis against experimentally observed frequencies and mode shapes. An iterative perturbation algorithm is derived and implemented to update nonlinear sensitivity (derivatives of eigenvalues and eigenvectors) during optimizer iterations, which reduced the number of finite element analyses.

  14. Modal test/analysis correlation of Space Station structures using nonlinear sensitivity

    NASA Technical Reports Server (NTRS)

    Gupta, Viney K.; Newell, James F.; Berke, Laszlo; Armand, Sasan

    1992-01-01

    The modal correlation problem is formulated as a constrained optimization problem for validation of finite element models (FEM's). For large-scale structural applications, a pragmatic procedure for substructuring, model verification, and system integration is described to achieve effective modal correlations. The space station substructure FEM's are reduced using Lanczos vectors and integrated into a system FEM using Craig-Bampton component modal synthesis. The optimization code is interfaced with MSC/NASTRAN to solve the problem of modal test/analysis correlation; that is, the problem of validating FEM's for launch and on-orbit coupled loads analysis against experimentally observed frequencies and mode shapes. An iterative perturbation algorithm is derived and implemented to update nonlinear sensitivity (derivatives of eigenvalues and eigenvectors) during optimizer iterations, which reduced the number of finite element analyses.

  15. Modal Test/Analysis Correlation of Space Station Structures Using Nonlinear Sensitivity

    NASA Technical Reports Server (NTRS)

    Gupta, Viney K.; Newell, James F.; Berke, Laszlo; Armand, Sasan

    1992-01-01

    The modal correlation problem is formulated as a constrained optimization problem for validation of finite element models (FEM's). For large-scale structural applications, a pragmatic procedure for substructuring, model verification, and system integration is described to achieve effective modal correlation. The space station substructure FEM's are reduced using Lanczos vectors and integrated into a system FEM using Craig-Bampton component modal synthesis. The optimization code is interfaced with MSC/NASTRAN to solve the problem of modal test/analysis correlation; that is, the problem of validating FEM's for launch and on-orbit coupled loads analysis against experimentally observed frequencies and mode shapes. An iterative perturbation algorithm is derived and implemented to update nonlinear sensitivity (derivatives of eigenvalues and eigenvectors) during optimizer iterations, which reduced the number of finite element analyses.

  16. Phase space theory of quantum–classical systems with nonlinear and stochastic dynamics

    SciTech Connect

    Burić, Nikola Popović, Duška B.; Radonjić, Milan; Prvanović, Slobodan

    2014-04-15

    A novel theory of hybrid quantum–classical systems is developed, utilizing the mathematical framework of constrained dynamical systems on the quantum–classical phase space. Both, the quantum and classical descriptions of the respective parts of the hybrid system are treated as fundamental. Therefore, the description of the quantum–classical interaction has to be postulated, and includes the effects of neglected degrees of freedom. Dynamical law of the theory is given in terms of nonlinear stochastic differential equations with Hamiltonian and gradient terms. The theory provides a successful dynamical description of the collapse during quantum measurement. -- Highlights: •A novel theory of quantum–classical systems is developed. •Framework of quantum constrained dynamical systems is used. •A dynamical description of the measurement induced collapse is obtained.

  17. Entire solutions of nonlocal dispersal equations with monostable nonlinearity in space periodic habitats

    NASA Astrophysics Data System (ADS)

    Li, Wan-Tong; Wang, Jia-Bing; Zhang, Li

    2016-08-01

    This paper is concerned with the new types of entire solutions other than traveling wave solutions of nonlocal dispersal equations with monostable nonlinearity in space periodic habitats. We first establish the existence and properties of spatially periodic solutions connecting two steady states. Then new types of entire solutions are constructed by combining the rightward and leftward pulsating traveling fronts with different speeds and a spatially periodic solution. Finally, for a class of special heterogeneous reaction, we further establish the uniqueness of entire solutions and the continuous dependence of such an entire solution on parameters, such as wave speeds and the shifted variables. In other words, we build a five-dimensional manifold of solutions and the traveling wave solutions are on the boundary of the manifold.

  18. Fault detection and isolation of PEM fuel cell system based on nonlinear analytical redundancy. An application via parity space approach

    NASA Astrophysics Data System (ADS)

    Aitouche, A.; Yang, Q.; Ould Bouamama, B.

    2011-05-01

    This paper presents a procedure dealing with the issue of fault detection and isolation (FDI) using nonlinear analytical redundancy (NLAR) technique applied in a proton exchange membrane (PEM) fuel cell system based on its mathematic model. The model is proposed and simplified into a five orders state space representation. The transient phenomena captured in the model include the compressor dynamics, the flow characteristics, mass and energy conservation and manifold fluidic mechanics. Nonlinear analytical residuals are generated based on the elimination of the unknown variables of the system by an extended parity space approach to detect and isolate actuator and sensor faults. Finally, numerical simulation results are given corresponding to a faults signature matrix.

  19. Non-linear generalised minimum variance control state space design for a second-order Volterra series model

    NASA Astrophysics Data System (ADS)

    Maboodi, Mohsen; Camacho, Eduardo F.; Khaki-Sedigh, Ali

    2015-10-01

    This paper presents a non-linear generalised minimum variance (NGMV) controller for a second-order Volterra series model with a general linear additive disturbance. The Volterra series models provide a natural extension of a linear convolution model with the nonlinearity considered in an additive term. The design procedure is entirely carried out in the state space framework, which facilitates the application of other analysis and design methods in this framework. First, the non-linear minimum variance (NMV) controller is introduced and then by changing the cost function, NGMV controller is defined as an extended version of the linear cases. The cost function is used in the simplest form and can be easily extended to the general case. Simulation results show the effectiveness of the proposed non-linear method.

  20. A space-time collocation scheme for modified anomalous subdiffusion and nonlinear superdiffusion equations

    NASA Astrophysics Data System (ADS)

    Bhrawy, A. H.

    2016-01-01

    This paper reports a new spectral collocation technique for solving time-space modified anomalous subdiffusion equation with a nonlinear source term subject to Dirichlet and Neumann boundary conditions. This model equation governs the evolution for the probability density function that describes anomalously diffusing particles. Anomalous diffusion is ubiquitous in physical and biological systems where trapping and binding of particles can occur. A space-time Jacobi collocation scheme is investigated for solving such problem. The main advantage of the proposed scheme is that, the shifted Jacobi Gauss-Lobatto collocation and shifted Jacobi Gauss-Radau collocation approximations are employed for spatial and temporal discretizations, respectively. Thereby, the problem is successfully reduced to a system of algebraic equations. The numerical results obtained by this algorithm have been compared with various numerical methods in order to demonstrate the high accuracy and efficiency of the proposed method. Indeed, for relatively limited number of Gauss-Lobatto and Gauss-Radau collocation nodes imposed, the absolute error in our numerical solutions is sufficiently small. The results have been compared with other techniques in order to demonstrate the high accuracy and efficiency of the proposed method.

  1. New Facts on the Nature of Gravitational Force And Nonlinear Oscillations of Space

    NASA Astrophysics Data System (ADS)

    Kursunoglu, Behram N.

    2002-07-01

    This paper discusses the letters received by this author from Albert Einstein, Erwin Schrodinger, and Paul Adrian Maurice Dirac, about fifty years ago which comment on my nonsymmetrical generalization of Einstein's general relativistic theory of gravitation. The writing of this paper, because of the dates of the letters, seems to have been delayed by half a century. Of the three versions of the nonsymmetrical theory (Einstein, Schrodinger and Kursunoglu Theories) my own paper contains results obtained as solutions of Generalized Theory of Gravitation field equations. In this paper it is shown that the field equations yield space nonlinear oscillations; a quartet of gravitational forces, quintessence, and replace Einstein's Cosmological Constant by two invariant parameters r0 and q related according to r02 q2 = c4/2G, where r0 is a length varying between zero and infinity and where q2 has the dimensions of energy density. These parameters govern the expansion of the universe with increasing acceleration and their existence yield four different solutions at each space-time point.

  2. Nonlinear Aeroacoustics Computations by the Space-Time CE/SE Method

    NASA Technical Reports Server (NTRS)

    Loh, Ching Y.

    2003-01-01

    The Space-Time Conservation Element and Solution Element Method, or CE/SE Method for short, is a recently developed numerical method for conservation laws. Despite its second order accuracy in space and time, it possesses low dispersion errors and low dissipation. The method is robust enough to cover a wide range of compressible flows: from weak linear acoustic waves to strong discontinuous waves (shocks). An outstanding feature of the CE/SE scheme is its truly multi-dimensional, simple but effective non-reflecting boundary condition (NRBC), which is particularly valuable for computational aeroacoustics (CAA). In nature, the method may be categorized as a finite volume method, where the conservation element (CE) is equivalent to a finite control volume (or cell) and the solution element (SE) can be understood as the cell interface. However, due to its careful treatment of the surface fluxes and geometry, it is different from the existing schemes. Currently, the CE/SE scheme has been developed to a matured stage that a 3-D unstructured CE/SE Navier-Stokes solver is already available. However, in the present review paper, as a general introduction to the CE/SE method, only the 2-D unstructured Euler CE/SE solver is chosen and sketched in section 2. Then applications of the 2-D and 3-D CE/SE schemes to linear, and in particular, nonlinear aeroacoustics are depicted in sections 3, 4, and 5 to demonstrate its robustness and capability.

  3. Molecular Dynamics of Flexible Polar Cations in a Variable Confined Space: Toward Exceptional Two-Step Nonlinear Optical Switches.

    PubMed

    Xu, Wei-Jian; He, Chun-Ting; Ji, Cheng-Min; Chen, Shao-Li; Huang, Rui-Kang; Lin, Rui-Biao; Xue, Wei; Luo, Jun-Hua; Zhang, Wei-Xiong; Chen, Xiao-Ming

    2016-07-01

    The changeable molecular dynamics of flexible polar cations in the variable confined space between inorganic chains brings about a new type of two-step nonlinear optical (NLO) switch with genuine "off-on-off" second harmonic generation (SHG) conversion between one NLO-active state and two NLO-inactive states. PMID:27159779

  4. Coxeter decompositions of hyperbolic simplexes

    SciTech Connect

    Felikson, A A

    2002-12-31

    A Coxeter decomposition of a polyhedron in a hyperbolic space H{sup n} is a decomposition of it into finitely many Coxeter polyhedra such that any two tiles having a common facet are symmetric with respect to it. The classification of Coxeter decompositions is closely related to the problem of the classification of finite-index subgroups generated by reflections in discrete hyperbolic groups generated by reflections. All Coxeter decompositions of simplexes in the hyperbolic spaces H{sup n} with n>3 are described in this paper.

  5. A reduced basis localized orthogonal decomposition

    NASA Astrophysics Data System (ADS)

    Abdulle, Assyr; Henning, Patrick

    2015-08-01

    In this work we combine the framework of the Reduced Basis method (RB) with the framework of the Localized Orthogonal Decomposition (LOD) in order to solve parametrized elliptic multiscale problems. The idea of the LOD is to split a high dimensional Finite Element space into a low dimensional space with comparably good approximation properties and a remainder space with negligible information. The low dimensional space is spanned by locally supported basis functions associated with the node of a coarse mesh obtained by solving decoupled local problems. However, for parameter dependent multiscale problems, the local basis has to be computed repeatedly for each choice of the parameter. To overcome this issue, we propose an RB approach to compute in an "offline" stage LOD for suitable representative parameters. The online solution of the multiscale problems can then be obtained in a coarse space (thanks to the LOD decomposition) and for an arbitrary value of the parameters (thanks to a suitable "interpolation" of the selected RB). The online RB-LOD has a basis with local support and leads to sparse systems. Applications of the strategy to both linear and nonlinear problems are given.

  6. Very large space structures: Non-linear control and robustness to structural uncertainties

    NASA Astrophysics Data System (ADS)

    Gasbarri, Paolo; Monti, Riccardo; Sabatini, Marco

    2014-01-01

    Very Large Space Structures (VLSS) are challenging systems to be controlled, due to their high flexibility. In particular, rapid attitude maneuvers can determine great oscillations on the flexible elements of a spacecraft (solar wings, antennas, booms). On account of this, in the last decades many researchers have developed different strategies to effectively damp the elastic vibrations by means of active vibration devices (such as piezo-electric patches) or by means of robust control algorithms. The approach suggested in this paper is different, since neither additional devices nor complex control laws are introduced. In fact, the complete model of the system (including rigid, elastic and orbital dynamics, coupled with control actions) is controlled by the non-linear attitude controller named state dependent Riccati equation, which will be based on a simplified version of the spacecraft model. The task to reduce the mutual interaction between rigid attitude and flexible dynamics is entirely transferred to a modification of the desired trajectory that must be tracked. This command shaping technique is based on the knowledge of the parameters (inertial and elastics) of the VLSS. Unfortunately these parameters are not always exactly known and, however, they may change over the time. On account of this a Monte Carlo analysis has been also performed, showing the robustness of the proposed control strategy to the structural uncertainties. The numerical simulations prove that this strategy, based on the joint application of two well-known yet simple techniques, produces accurate and robust results.

  7. Properties of soil pore space regulate pathways of plant residue decomposition and community structure of associated bacteria.

    PubMed

    Negassa, Wakene C; Guber, Andrey K; Kravchenko, Alexandra N; Marsh, Terence L; Hildebrandt, Britton; Rivers, Mark L

    2015-01-01

    Physical protection of soil carbon (C) is one of the important components of C storage. However, its exact mechanisms are still not sufficiently lucid. The goal of this study was to explore the influence of soil structure, that is, soil pore spatial arrangements, with and without presence of plant residue on (i) decomposition of added plant residue, (ii) CO2 emission from soil, and (iii) structure of soil bacterial communities. The study consisted of several soil incubation experiments with samples of contrasting pore characteristics with/without plant residue, accompanied by X-ray micro-tomographic analyses of soil pores and by microbial community analysis of amplified 16S-18S rRNA genes via pyrosequencing. We observed that in the samples with substantial presence of air-filled well-connected large (>30 µm) pores, 75-80% of the added plant residue was decomposed, cumulative CO2 emission constituted 1,200 µm C g(-1) soil, and movement of C from decomposing plant residue into adjacent soil was insignificant. In the samples with greater abundance of water-filled small pores, 60% of the added plant residue was decomposed, cumulative CO2 emission constituted 2,000 µm C g(-1) soil, and the movement of residue C into adjacent soil was substantial. In the absence of plant residue the influence of pore characteristics on CO2 emission, that is on decomposition of the native soil organic C, was negligible. The microbial communities on the plant residue in the samples with large pores had more microbial groups known to be cellulose decomposers, that is, Bacteroidetes, Proteobacteria, Actinobacteria, and Firmicutes, while a number of oligotrophic Acidobacteria groups were more abundant on the plant residue from the samples with small pores. This study provides the first experimental evidence that characteristics of soil pores and their air/water flow status determine the phylogenetic composition of the local microbial community and directions and magnitudes of soil C

  8. Properties of Soil Pore Space Regulate Pathways of Plant Residue Decomposition and Community Structure of Associated Bacteria

    PubMed Central

    Negassa, Wakene C.; Guber, Andrey K.; Kravchenko, Alexandra N.; Marsh, Terence L.; Hildebrandt, Britton; Rivers, Mark L.

    2015-01-01

    Physical protection of soil carbon (C) is one of the important components of C storage. However, its exact mechanisms are still not sufficiently lucid. The goal of this study was to explore the influence of soil structure, that is, soil pore spatial arrangements, with and without presence of plant residue on (i) decomposition of added plant residue, (ii) CO2 emission from soil, and (iii) structure of soil bacterial communities. The study consisted of several soil incubation experiments with samples of contrasting pore characteristics with/without plant residue, accompanied by X-ray micro-tomographic analyses of soil pores and by microbial community analysis of amplified 16S–18S rRNA genes via pyrosequencing. We observed that in the samples with substantial presence of air-filled well-connected large (>30 µm) pores, 75–80% of the added plant residue was decomposed, cumulative CO2 emission constituted 1,200 µm C g-1 soil, and movement of C from decomposing plant residue into adjacent soil was insignificant. In the samples with greater abundance of water-filled small pores, 60% of the added plant residue was decomposed, cumulative CO2 emission constituted 2,000 µm C g-1 soil, and the movement of residue C into adjacent soil was substantial. In the absence of plant residue the influence of pore characteristics on CO2 emission, that is on decomposition of the native soil organic C, was negligible. The microbial communities on the plant residue in the samples with large pores had more microbial groups known to be cellulose decomposers, that is, Bacteroidetes, Proteobacteria, Actinobacteria, and Firmicutes, while a number of oligotrophic Acidobacteria groups were more abundant on the plant residue from the samples with small pores. This study provides the first experimental evidence that characteristics of soil pores and their air/water flow status determine the phylogenetic composition of the local microbial community and directions and magnitudes of soil C

  9. Properties of soil pore space regulate pathways of plant residue decomposition and community structure of associated bacteria

    SciTech Connect

    Negassa, Wakene C.; Guber, Andrey K.; Kravchenko, Alexandra N.; Marsh, Terence L.; Hildebrandt, Britton; Rivers, Mark L.

    2015-07-01

    Physical protection of soil carbon (C) is one of the important components of C storage. However, its exact mechanisms are still not sufficiently lucid. The goal of this study was to explore the influence of soil structure, that is, soil pore spatial arrangements, with and without presence of plant residue on (i) decomposition of added plant residue, (ii) CO₂ emission from soil, and (iii) structure of soil bacterial communities. The study consisted of several soil incubation experiments with samples of contrasting pore characteristics with/without plant residue, accompanied by X-ray micro-tomographic analyses of soil pores and by microbial community analysis of amplified 16S–18S rRNA genes via pyrosequencing. We observed that in the samples with substantial presence of air-filled well-connected large (>30 µm) pores, 75–80% of the added plant residue was decomposed, cumulative CO₂ emission constituted 1,200 µm C g⁻¹ soil, and movement of C from decomposing plant residue into adjacent soil was insignificant. In the samples with greater abundance of water-filled small pores, 60% of the added plant residue was decomposed, cumulative CO₂ emission constituted 2,000 µm C g⁻¹ soil, and the movement of residue C into adjacent soil was substantial. In the absence of plant residue the influence of pore characteristics on CO₂ emission, that is on decomposition of the native soil organic C, was negligible. The microbial communities on the plant residue in the samples with large pores had more microbial groups known to be cellulose decomposers, that is, Bacteroidetes, Proteobacteria, Actinobacteria, and Firmicutes, while a number of oligotrophic Acidobacteria groups were more abundant on the plant residue from the samples with small pores. This study provides the first experimental evidence that characteristics of soil pores and their air/water flow status determine the phylogenetic composition of the local microbial community and directions and magnitudes of

  10. Optimization by nonhierarchical asynchronous decomposition

    NASA Technical Reports Server (NTRS)

    Shankar, Jayashree; Ribbens, Calvin J.; Haftka, Raphael T.; Watson, Layne T.

    1992-01-01

    Large scale optimization problems are tractable only if they are somehow decomposed. Hierarchical decompositions are inappropriate for some types of problems and do not parallelize well. Sobieszczanski-Sobieski has proposed a nonhierarchical decomposition strategy for nonlinear constrained optimization that is naturally parallel. Despite some successes on engineering problems, the algorithm as originally proposed fails on simple two dimensional quadratic programs. The algorithm is carefully analyzed for quadratic programs, and a number of modifications are suggested to improve its robustness.

  11. Rice-Ramsperger-Kassel-Marcus theoretical prediction of high-pressure arrhenius parameters by nonlinear regression: application to silane and disilane decomposition

    SciTech Connect

    Roenigk, K.F.; Jensen, K.F.; Carr, R.W.

    1987-10-22

    Arrhenius parameters are estimated for silane and disilane thermal decomposition reactions by direct regression of RRKM predictions on published static and shock-tube data. For silane decomposition, they find E/sub infinity/ = 57.4-61.1 kcal/mol and log A/sub infinity/ = 14.9-16.3, while for disilane they find E/sub infinity/ = 51.1-52.5 kcal/mol and log A/sub infinity/ = 15.2-16.2. The lower limiting values correspond to inclusion of negative temperature dependence in the collision efficiency, while the higher values correspond to inclusion of weak or negligible temperature dependence. The Arrhenium parameters for both silane and disilane decomposition differ substantially from previously published values. For silane, they predict preexponentials approximately an order of magnitude greater than the previous values for the same activation energy. For disilane, they find A/sub infinity/ is roughly an order of magnitude higher than the literature values and E/sub infinity/ is greater by more than 2 kcal/mol. Falloff curves for both silane and disilane decomposition are given. Implications of these results for the activation energy of SiH/sub 2/ insertion into H/sub 2/ and SiH/sub 4/ and for ..delta..H/sub f//sup 0/(SiH/sub 2/) are discussed.

  12. Identification of Non-Linear Space Weather Models of the Van Allen Radiation Belts Using Volterra Networks

    NASA Astrophysics Data System (ADS)

    Taylor, M.; Daglis, I. A.; Anastasiadis, A.; Vassiliadis, D.

    2010-07-01

    Many efforts have been made to develop general dynamical models of the Van Allen radiation belts based on data alone. Early linear prediction filter studies focused on the response of daily-averaged relativistic electrons at geostationary altitudes Nagai 1988, Baker et al. 1990). Vassiliadis et al (2005) extended this technique spatially by incorporating SAMPEX electron flux data into linear prediction filters for a broad range of L-shells from 1.1 to 10.0 RE. Nonlinear state space models (Rigler & Baker 2008) have provided useful initial results on the timescales involved in modeling the impulse-response of the radiation belts. Here, we show how NARMAX models, in conjunction with nonlinear time-delay FIR neural networks (Volterra networks) hold great promise for the development of accurate and fully data-derived space weather specification and forecast tools.

  13. Properties of soil pore space regulate pathways of plant residue decomposition and community structure of associated bacteria

    DOE PAGESBeta

    Negassa, Wakene C.; Guber, Andrey K.; Kravchenko, Alexandra N.; Marsh, Terence L.; Hildebrandt, Britton; Rivers, Mark L.

    2015-07-01

    Physical protection of soil carbon (C) is one of the important components of C storage. However, its exact mechanisms are still not sufficiently lucid. The goal of this study was to explore the influence of soil structure, that is, soil pore spatial arrangements, with and without presence of plant residue on (i) decomposition of added plant residue, (ii) CO₂ emission from soil, and (iii) structure of soil bacterial communities. The study consisted of several soil incubation experiments with samples of contrasting pore characteristics with/without plant residue, accompanied by X-ray micro-tomographic analyses of soil pores and by microbial community analysis ofmore » amplified 16S–18S rRNA genes via pyrosequencing. We observed that in the samples with substantial presence of air-filled well-connected large (>30 µm) pores, 75–80% of the added plant residue was decomposed, cumulative CO₂ emission constituted 1,200 µm C g⁻¹ soil, and movement of C from decomposing plant residue into adjacent soil was insignificant. In the samples with greater abundance of water-filled small pores, 60% of the added plant residue was decomposed, cumulative CO₂ emission constituted 2,000 µm C g⁻¹ soil, and the movement of residue C into adjacent soil was substantial. In the absence of plant residue the influence of pore characteristics on CO₂ emission, that is on decomposition of the native soil organic C, was negligible. The microbial communities on the plant residue in the samples with large pores had more microbial groups known to be cellulose decomposers, that is, Bacteroidetes, Proteobacteria, Actinobacteria, and Firmicutes, while a number of oligotrophic Acidobacteria groups were more abundant on the plant residue from the samples with small pores. This study provides the first experimental evidence that characteristics of soil pores and their air/water flow status determine the phylogenetic composition of the local microbial community and directions and

  14. The soliton transform and a possible application to nonlinear Alfven waves in space

    NASA Technical Reports Server (NTRS)

    Hada, T.; Hamilton, R. L.; Kennel, C. F.

    1993-01-01

    The inverse scattering transform (IST) based on the derivative nonlinear Schroedinger (DNLS) equation is applied to a complex time series of nonlinear Alfven wave data generated by numerical simulation. The IST describes the long-time evolution of quasi-parallel Alfven waves more efficiently than the Fourier transform, which is adapted to linear rather than nonlinear problems. When dissipation is added, so the conditions for the validity of the DNLS are not strictly satisfied, the IST continues to provide a compact description of the wavefield in terms of a small number of decaying envelope solitons.

  15. Nonlinear dust phase-space vortices (holes) in charge-varying dusty plasmas

    SciTech Connect

    Tribeche, Mouloud; Aoutou, Kamel; Zerguini, Taha Houssine

    2005-03-01

    The recent analysis of dust voids [A. A. Mamun and P. K. Shukla, Phys. Plasmas 11, 1757 (2004)] is extended to include self-consistently the dust charge variation. Numerical solutions of highly nonlinear equations are carried out including dust charging and dust trapping. It is found that under certain conditions the effect of dust charge variation can be quite important. In particular, it may be noted that the dust charge variation leads to an additional enlargement of the nonlinear dust voids. The effects of ion/electron temperature, trapping parameter, and dust size on the properties of these nonlinear dust voids are briefly discussed.

  16. Unconditionally energy stable time stepping scheme for Cahn-Morral equation: Application to multi-component spinodal decomposition and optimal space tiling

    NASA Astrophysics Data System (ADS)

    Tavakoli, Rouhollah

    2016-01-01

    An unconditionally energy stable time stepping scheme is introduced to solve Cahn-Morral-like equations in the present study. It is constructed based on the combination of David Eyre's time stepping scheme and Schur complement approach. Although the presented method is general and independent of the choice of homogeneous free energy density function term, logarithmic and polynomial energy functions are specifically considered in this paper. The method is applied to study the spinodal decomposition in multi-component systems and optimal space tiling problems. A penalization strategy is developed, in the case of later problem, to avoid trivial solutions. Extensive numerical experiments demonstrate the success and performance of the presented method. According to the numerical results, the method is convergent and energy stable, independent of the choice of time stepsize. Its MATLAB implementation is included in the appendix for the numerical evaluation of algorithm and reproduction of the presented results.

  17. Infectious diseases in space and time: noise and nonlinearity in epidemiological dynamics

    NASA Astrophysics Data System (ADS)

    Grenfell, Bryan

    2005-03-01

    I illustrate the impact of noise and nonlinearity on the spatio-temporal dynamics and evolution of epidemics using mathematical models and analyses of detailed epidemiological data from childhood infections, such as measles.

  18. Low-dimensional models for the nonlinear vibration analysis of cylindrical shells based on a perturbation procedure and proper orthogonal decomposition

    NASA Astrophysics Data System (ADS)

    Gonçalves, P. B.; Silva, F. M. A.; Del Prado, Z. J. G. N.

    2008-08-01

    In formulating mathematical models for dynamical systems, obtaining a high degree of qualitative correctness (i.e. predictive capability) may not be the only objective. The model must be useful for its intended application, and models of reduced complexity are attractive in many cases where time-consuming numerical procedures are required. This paper discusses the derivation of discrete low-dimensional models for the nonlinear vibration analysis of thin cylindrical shells. In order to understand the peculiarities inherent to this class of structural problems, the nonlinear vibrations and dynamic stability of a circular cylindrical shell subjected to static and dynamic loads are analyzed. This choice is based on the fact that cylindrical shells exhibit a highly nonlinear behavior under both static and dynamic loads. Geometric nonlinearities due to finite-amplitude shell motions are considered by using Donnell's nonlinear shallow-shell theory. A perturbation procedure, validated in previous studies, is used to derive a general expression for the nonlinear vibration modes and the discretized equations of motion are obtained by the Galerkin method using modal expansions for the displacements that satisfy all the relevant boundary and symmetry conditions. Next, the model is analyzed via the Karhunen-Loève expansion to investigate the relative importance of each mode obtained by the perturbation solution on the nonlinear response and total energy of the system. The responses of several low-dimensional models are compared. It is shown that rather low-dimensional but properly selected models can describe with good accuracy the response of the shell up to very large vibration amplitudes.

  19. Ozone decomposition

    PubMed Central

    Batakliev, Todor; Georgiev, Vladimir; Anachkov, Metody; Rakovsky, Slavcho

    2014-01-01

    Catalytic ozone decomposition is of great significance because ozone is a toxic substance commonly found or generated in human environments (aircraft cabins, offices with photocopiers, laser printers, sterilizers). Considerable work has been done on ozone decomposition reported in the literature. This review provides a comprehensive summary of the literature, concentrating on analysis of the physico-chemical properties, synthesis and catalytic decomposition of ozone. This is supplemented by a review on kinetics and catalyst characterization which ties together the previously reported results. Noble metals and oxides of transition metals have been found to be the most active substances for ozone decomposition. The high price of precious metals stimulated the use of metal oxide catalysts and particularly the catalysts based on manganese oxide. It has been determined that the kinetics of ozone decomposition is of first order importance. A mechanism of the reaction of catalytic ozone decomposition is discussed, based on detailed spectroscopic investigations of the catalytic surface, showing the existence of peroxide and superoxide surface intermediates. PMID:26109880

  20. Ozone decomposition.

    PubMed

    Batakliev, Todor; Georgiev, Vladimir; Anachkov, Metody; Rakovsky, Slavcho; Zaikov, Gennadi E

    2014-06-01

    Catalytic ozone decomposition is of great significance because ozone is a toxic substance commonly found or generated in human environments (aircraft cabins, offices with photocopiers, laser printers, sterilizers). Considerable work has been done on ozone decomposition reported in the literature. This review provides a comprehensive summary of the literature, concentrating on analysis of the physico-chemical properties, synthesis and catalytic decomposition of ozone. This is supplemented by a review on kinetics and catalyst characterization which ties together the previously reported results. Noble metals and oxides of transition metals have been found to be the most active substances for ozone decomposition. The high price of precious metals stimulated the use of metal oxide catalysts and particularly the catalysts based on manganese oxide. It has been determined that the kinetics of ozone decomposition is of first order importance. A mechanism of the reaction of catalytic ozone decomposition is discussed, based on detailed spectroscopic investigations of the catalytic surface, showing the existence of peroxide and superoxide surface intermediates. PMID:26109880

  1. Optical ranked-order filtering using threshold decomposition

    DOEpatents

    Allebach, J.P.; Ochoa, E.; Sweeney, D.W.

    1987-10-09

    A hybrid optical/electronic system performs median filtering and related ranked-order operations using threshold decomposition to encode the image. Threshold decomposition transforms the nonlinear neighborhood ranking operation into a linear space-invariant filtering step followed by a point-to-point threshold comparison step. Spatial multiplexing allows parallel processing of all the threshold components as well as recombination by a second linear, space-invariant filtering step. An incoherent optical correlation system performs the linear filtering, using a magneto-optic spatial light modulator as the input device and a computer-generated hologram in the filter plane. Thresholding is done electronically. By adjusting the value of the threshold, the same architecture is used to perform median, minimum, and maximum filtering of images. A totally optical system is also disclosed. 3 figs.

  2. Optical ranked-order filtering using threshold decomposition

    DOEpatents

    Allebach, Jan P.; Ochoa, Ellen; Sweeney, Donald W.

    1990-01-01

    A hybrid optical/electronic system performs median filtering and related ranked-order operations using threshold decomposition to encode the image. Threshold decomposition transforms the nonlinear neighborhood ranking operation into a linear space-invariant filtering step followed by a point-to-point threshold comparison step. Spatial multiplexing allows parallel processing of all the threshold components as well as recombination by a second linear, space-invariant filtering step. An incoherent optical correlation system performs the linear filtering, using a magneto-optic spatial light modulator as the input device and a computer-generated hologram in the filter plane. Thresholding is done electronically. By adjusting the value of the threshold, the same architecture is used to perform median, minimum, and maximum filtering of images. A totally optical system is also disclosed.

  3. Extended phase space of AdS black holes in Einstein-Gauss-Bonnet gravity with a quadratic nonlinear electrodynamics

    NASA Astrophysics Data System (ADS)

    Hendi, S. H.; Panahiyan, S.; Momennia, M.

    2016-04-01

    In this paper, we consider quadratic Maxwell invariant as a correction to the Maxwell theory and study thermodynamic behavior of the black holes in Einstein and Gauss-Bonnet gravities. We consider cosmological constant as a thermodynamic pressure to extend phase space. Next, we obtain critical values in case of variation of nonlinearity and Gauss-Bonnet parameters. Although the general thermodynamical behavior of the black hole solutions is the same as usual Van der Waals system, we show that in special case of the nonlinear electromagnetic field, there will be a turning point for the phase diagrams and usual Van der Waals is not observed. This theory of nonlinear electromagnetic field provides two critical horizon radii. We show that this unusual behavior of phase diagrams is due to existence of second critical horizon radius. It will be pointed out that the power of the gravity and nonlinearity of the matter field modify the critical values. We generalize the study by considering the effects of dimensionality on critical values and make comparisons between our models with their special sub-classes. In addition, we examine the possibility of the existence of the reentrant phase transitions through two different methods.

  4. Nonlinear dynamic response of a uni-directional model for the tile/pad space shuttle thermal protection system

    NASA Technical Reports Server (NTRS)

    Housner, J. M.; Edighoffer, H. H.; Park, K. C.

    1980-01-01

    A unidirectional analysis of the nonlinear dynamic behavior of the space shuttle tile/pad thermal protection system is developed and examined for imposed sinusoidal and random motions of the shuttle skin and/or applied tile pressure. The analysis accounts for the highly nonlinear stiffening hysteresis and viscous behavior of the pad which joins the tile to the shuttle skin. Where available, experimental data are used to confirm the validity of the analysis. Both analytical and experimental studies reveal that the system resonant frequency is very high for low amplitude oscillations but decreases rapidly to a minimum value with increasing amplitude. Analytical studies indicate that with still higher amplitude the resonant frequency increases slowly. The nonlinear pad is also responsible for the analytically and experimentally observed distorted response wave shapes having high sharp peaks when the system is subject to sinusoidal loads. Furthermore, energy dissipation in the pad is studied analytically and it is found that the energy dissipated is sufficiently high to cause rapid decay of dynamic transients. Nevertheless, the sharp peaked nonlinear responses of the system lead to higher magnification factors than would be expected in such a highly damped linear system.

  5. Nonlinear space charge dynamics in mixed ionic-electronic conductors: Resistive switching and ferroelectric-like hysteresis of electromechanical response

    SciTech Connect

    Morozovska, Anna N.; Morozovsky, Nicholas V.; Eliseev, Eugene A.; Varenyk, Olexandr V.; Kim, Yunseok; Strelcov, Evgheni; Tselev, Alexander; Kalinin, Sergei V.

    2014-08-14

    We performed self-consistent modelling of nonlinear electrotransport and electromechanical response of thin films of mixed ionic-electronic conductors (MIEC) allowing for steric effects of mobile charged defects (ions, protons, or vacancies), electron degeneration, and Vegard stresses. We establish correlations between the features of the nonlinear space-charge dynamics, current-voltage, and bending-voltage curves for different types of the film electrodes. A pronounced ferroelectric-like hysteresis of the bending-voltage loops and current maxima on the double hysteresis current-voltage loops appear for the electron-transport electrodes. The double hysteresis loop with pronounced humps indicates a memristor-type resistive switching. The switching occurs due to the strong nonlinear coupling between the electronic and ionic subsystems. A sharp meta-stable maximum of the electron density appears near one open electrode and moves to another one during the periodic change of applied voltage. Our results can explain the nonlinear nature and correlation of electrical and mechanical memory effects in thin MIEC films. The analytical expression proving that the electrically induced bending of MIEC films can be detected by interferometric methods is derived.

  6. Development, analysis, and testing of robust nonlinear guidance algorithms for space applications

    NASA Astrophysics Data System (ADS)

    Wibben, Daniel R.

    This work focuses on the analysis and application of various nonlinear, autonomous guidance algorithms that utilize sliding mode control to guarantee system stability and robustness. While the basis for the algorithms has previously been proposed, past efforts barely scratched the surface of the theoretical details and implications of these algorithms. Of the three algorithms that are the subject of this research, two are directly derived from optimal control theory and augmented using sliding mode control. Analysis of the derivation of these algorithms has shown that they are two different representations of the same result, one of which uses a simple error state model (Delta r/Deltav) and the other uses definitions of the zero-effort miss and zero-effort velocity (ZEM/ZEV) values. By investigating the dynamics of the defined sliding surfaces and their impact on the overall system, many implications have been deduced regarding the behavior of these systems which are noted to feature time-varying sliding modes. A formal finite time stability analysis has also been performed to theoretically demonstrate that the algorithms globally stabilize the system in finite time in the presence of perturbations and unmodeled dynamics. The third algorithm that has been subject to analysis is derived from a direct application of higher-order sliding mode control and Lyapunov stability analysis without consideration of optimal control theory and has been named the Multiple Sliding Surface Guidance (MSSG). Via use of reinforcement learning methods an optimal set of gains has been found that make the guidance perform similarly to an open-loop optimal solution. Careful side-by-side inspection of the MSSG and Optimal Sliding Guidance (OSG) algorithms has shown some striking similarities. A detailed comparison of the algorithms has demonstrated that though they are nearly indistinguishable at first glance, there are some key differences between the two algorithms and they are indeed

  7. Decomposition techniques

    USGS Publications Warehouse

    Chao, T.T.; Sanzolone, R.F.

    1992-01-01

    Sample decomposition is a fundamental and integral step in the procedure of geochemical analysis. It is often the limiting factor to sample throughput, especially with the recent application of the fast and modern multi-element measurement instrumentation. The complexity of geological materials makes it necessary to choose the sample decomposition technique that is compatible with the specific objective of the analysis. When selecting a decomposition technique, consideration should be given to the chemical and mineralogical characteristics of the sample, elements to be determined, precision and accuracy requirements, sample throughput, technical capability of personnel, and time constraints. This paper addresses these concerns and discusses the attributes and limitations of many techniques of sample decomposition along with examples of their application to geochemical analysis. The chemical properties of reagents as to their function as decomposition agents are also reviewed. The section on acid dissolution techniques addresses the various inorganic acids that are used individually or in combination in both open and closed systems. Fluxes used in sample fusion are discussed. The promising microwave-oven technology and the emerging field of automation are also examined. A section on applications highlights the use of decomposition techniques for the determination of Au, platinum group elements (PGEs), Hg, U, hydride-forming elements, rare earth elements (REEs), and multi-elements in geological materials. Partial dissolution techniques used for geochemical exploration which have been treated in detail elsewhere are not discussed here; nor are fire-assaying for noble metals and decomposition techniques for X-ray fluorescence or nuclear methods be discussed. ?? 1992.

  8. On the Hodge-type decomposition and cohomology groups of k-Cauchy-Fueter complexes over domains in the quaternionic space

    NASA Astrophysics Data System (ADS)

    Chang, Der-Chen; Markina, Irina; Wang, Wei

    2016-09-01

    The k-Cauchy-Fueter operator D0(k) on one dimensional quaternionic space H is the Euclidean version of spin k / 2 massless field operator on the Minkowski space in physics. The k-Cauchy-Fueter equation for k ≥ 2 is overdetermined and its compatibility condition is given by the k-Cauchy-Fueter complex. In quaternionic analysis, these complexes play the role of Dolbeault complex in several complex variables. We prove that a natural boundary value problem associated to this complex is regular. Then by using the theory of regular boundary value problems, we show the Hodge-type orthogonal decomposition, and the fact that the non-homogeneous k-Cauchy-Fueter equation D0(k) u = f on a smooth domain Ω in H is solvable if and only if f satisfies the compatibility condition and is orthogonal to the set ℋ(k)1 (Ω) of Hodge-type elements. This set is isomorphic to the first cohomology group of the k-Cauchy-Fueter complex over Ω, which is finite dimensional, while the second cohomology group is always trivial.

  9. Nonlinear Electromagnetic Waves and Spherical Arc-Polarized Waves in Space Plasmas

    NASA Technical Reports Server (NTRS)

    Tsurutani, B.; Ho, Christian M.; Arballo, John K.; Lakhina, Gurbax S.; Glassmeier, Karl-Heinz; Neubauer, Fritz M.

    1997-01-01

    We review observations of nonlinear plasma waves detected by interplanetary spacecraft. For this paper we will focus primarily on the phase-steepened properties of such waves. Plasma waves at comet Giacobini-Zinner measured by the International Cometary Explorer (ICE), at comets Halley and Grigg-Skjellerup measured by Giotto, and interplanetary Alfven waves measured by Ulysses, will be discussed and intercompared.

  10. Stability study of the Large Space Telescope /LST/ system with nonlinear CMG gimbal friction. [Control Moment Gyroscopes

    NASA Technical Reports Server (NTRS)

    Kuo, B. C.; Singh, G.; Seltzer, S. M.

    1974-01-01

    The purpose of the investigation reported upon is to study the existence and characteristics of self-sustained oscillations in the dynamic behavior of the Large Space Telescope (LST) system due to the presence of nonlinear gimbal friction in the control moment gyroscopes (CMGs). A continuous-data single-axis model of the LST is considered. A solid friction model is used to represent CMG gimbal friction. A rigorous mathematical model is derived for use in a continuous describing function analysis. Conditions for self-sustained oscillations are then determined.

  11. Nonlinear Stability of Strong Planar Rarefaction Waves for the Relaxation Approximation of Conservation Laws in Several Space Dimensions

    NASA Astrophysics Data System (ADS)

    Zhao, Huijiang

    2000-04-01

    In this paper, we show that a strong planar rarefaction wave is nonlinear stable, namely it is an attractor for the relaxation approximation of the scalar conservation laws in several space dimensions. Compared with former results obtained by T. P. Liu (1987, Comm. Math. Phys.108, 153-175) and T. Luo (1997, J. Differential Equations133, 255-279), our main novelty lies in the fact that the planar rarefaction waves do not need to be small, and in the one-dimensional case, the initial disturbance can also be chosen arbitrarily large.

  12. Linear and non-linear studies of Alfven waves in space. Stationary and dynamic processes in magnetospheric plasmas

    NASA Technical Reports Server (NTRS)

    Bhattacharjee, A.; Hasegawa, A.

    1990-01-01

    The Final Technical Report on linear and non-linear studies of Alfven waves in space is presented. Areas of research included relaxation of magnetotail plasmas with field-aligned currents; the equilibrium dayside magnetosphere; macroscale particle simulation of kinetic Alfven wave physics; ballooning stability of plasmas with sheared equilibrium flows; theory of the drift-mirror instability; collisionless tearing instability in magnetotail plasmas; and nonadiabatic behavior of the magnetic moment of a charged particle in a dipole magnetic field and the development of stochastic webs.

  13. The temperature effect on the glycine decomposition induced by 2 keV electron bombardment in space analog conditions

    NASA Astrophysics Data System (ADS)

    Pilling, Sergio; Nair, Binu G.; Escobar, Antonio; Fraser, Helen; Mason, Nigel

    2014-03-01

    Glycine is the simplest proteinaceous amino acid that has been extensively detected in carbonaceous meteorites and was recently observed in the cometary samples returned to Earth by NASA's Stardust spacecraft. In space, such species is exposed to several radiation fields at different temperatures. In aqueous solutions, this species appears mainly as zwitterionic glycine (+NH3CH2COO-) however, in solid phase, it may be found in amorphous or crystalline forms. Here, we present an experimental study on the destruction of two zwitterionic glycine crystals ( α- and β-form) at two different temperatures (300 K and 14 K) by 2 keV electrons in an attempt to test the behavior and stability of this molecular species in different space environments. The samples were analyzed in situ by Fourier transform infrared spectrometry at electron fluences. The experiments were carried out under ultra-high vacuum conditions at the Molecular Physics Laboratory at the Open University at Milton Keynes, UK. The dissociation cross section of glycine is approximately 5 times higher for the 14 K samples when compared to the 300 K samples. In contrast, no significant differences emerged between the dissociation cross sections of α- and β-forms of glycine for fixed temperature experiments. We therefore conclude that the destruction cross section is more heavily dependent on temperature than the phase of the condensed glycine material. This may be associated with the opening of additional reaction routes in the frozen samples involving the trapped daughter species (e.g. CO2 and CO). The half-life of studied samples extrapolated to space conditions shows that glycine molecules on the surface of interstellar grains has less survivability and they are highly sensitive to ambient radiations, however, they can survive extended period of time in the solar system like environments. Survivability increases by a factor of 5 if the samples are at 300 K when compared to low temperature experiments at 14

  14. Kinematic dust viscosity effect on linear and nonlinear dust-acoustic waves in space dusty plasmas with nonthermal ions

    NASA Astrophysics Data System (ADS)

    El-Hanbaly, A. M.; Sallah, M.; El-Shewy, E. K.; Darweesh, H. F.

    2015-10-01

    Linear and nonlinear dust-acoustic (DA) waves are studied in a collisionless, unmagnetized and dissipative dusty plasma consisting of negatively charged dust grains, Boltzmann-distributed electrons, and nonthermal ions. The normal mode analysis is used to obtain a linear dispersion relation illustrating the dependence of the wave damping rate on the carrier wave number, the dust viscosity coefficient, the ratio of the ion temperature to the electron temperatures, and the nonthermal parameter. The plasma system is analyzed nonlinearly via the reductive perturbation method that gives the KdV-Burgers equation. Some interesting physical solutions are obtained to study the nonlinear waves. These solutions are related to soliton, a combination between a shock and a soliton, and monotonic and oscillatory shock waves. Their behaviors are illustrated and shown graphically. The characteristics of the DA solitary and shock waves are significantly modified by the presence of nonthermal (fast) ions, the ratio of the ion temperature to the electron temperature, and the dust kinematic viscosity. The topology of the phase portrait and the potential diagram of the KdV-Burgers equation is illustrated, whose advantage is the ability to predict different classes of traveling wave solutions according to different phase orbits. The energy of the soliton wave and the electric field are calculated. The results in this paper can be generalized to analyze the nature of plasma waves in both space and laboratory plasma systems.

  15. Nonlinear Analysis of the Space Shuttle Superlightweight LO2 Tank. Part 1; Bahavior Under Booster Ascent Loads

    NASA Technical Reports Server (NTRS)

    Young, Richard D.; Nemeth, Michael P.; Collins, Timothy J.; Starnes, James H., Jr.

    1998-01-01

    Results of linear bifurcation and nonlinear analyses of the Space Shuttle superlightweight (SLWT) external liquid-oxygen (LO2) tank for an important early booster ascent loading condition are presented. These results for thin-walled linear elastic shells that are subjected to combined mechanical and thermal loads illustrate an important type of response mode that may be encountered in the design of other liquid-fuel launch vehicles. Linear bifurcation analyses are presented that predict several nearly equal eigenvalues that correspond to local buckling modes in the forward ogive section of the LO2 tank. In contrast, the nonlinear response phenomenon is shown to consist of short-wavelength bending deformations in the forward ogive and barrel sections of the LO2 tank that growing amplitude in a stable manner increasing load. Imperfection sensitivity analyses are presented that show that the presence of several nearly equal eigenvalues does not lead to a premature general instability mode for the forward ogive section. For the linear bifurcation and nonlinear analyses, the results show that accurate predictions of the response of the shield generally require a large-scale, high-fidelity finite-element model. Results are also presented that show that the SLWT LO2 tank can support loads in excess of approximately 2.6 times the values of the operational loads considered.

  16. Higher-order corrections to nonlinear dust-ion-acoustic shock waves in a degenerate dense space plasma

    NASA Astrophysics Data System (ADS)

    El-Labany, S. K.; El-Taibany, W. F.; El-Samahy, A. E.; Hafez, A. M.; Atteya, A.

    2014-12-01

    A reductive perturbation technique is employed to investigate the contribution of higher-order nonlinearity and dissipation to nonlinear dust-ion-acoustic (DIA) shock waves in a three-component degenerate dense space plasma. The model consists of degenerate electron (being either ultrarelativistic or nonrelativistic), nonrelativistic ion fluid and stationary heavy dust grains. A nonlinear Burger equation and a linear inhomogeneous Burger-type equation are derived. The present model admits only compressive DIA shocks. Including these higher-order corrections results in creating new solitary wave structures " humped DIA shock" waves. For the case of ultrarelativistic (nonrelativistic) electrons, one (two) humped DIA shock is (are) created. The DIA shock wave amplitude and velocity is larger in case of ultrarelativistic electrons than of nonrelativistic electrons. It is shown that the effects of kinematic viscosity, heavy dust grains number density, and equilibrium ion number density have important roles in the basic features of the produced DIA shocks and the associated electric fields. The implications of our results to dense plasmas in astrophysical objects (e.g., non-rotating white dwarf stars) are briefly discussed.

  17. Kinematic dust viscosity effect on linear and nonlinear dust-acoustic waves in space dusty plasmas with nonthermal ions

    SciTech Connect

    El-Hanbaly, A. M.; Sallah, M.; El-Shewy, E. K.; Darweesh, H. F.

    2015-10-15

    Linear and nonlinear dust-acoustic (DA) waves are studied in a collisionless, unmagnetized and dissipative dusty plasma consisting of negatively charged dust grains, Boltzmann-distributed electrons, and nonthermal ions. The normal mode analysis is used to obtain a linear dispersion relation illustrating the dependence of the wave damping rate on the carrier wave number, the dust viscosity coefficient, the ratio of the ion temperature to the electron temperatures, and the nonthermal parameter. The plasma system is analyzed nonlinearly via the reductive perturbation method that gives the KdV-Burgers equation. Some interesting physical solutions are obtained to study the nonlinear waves. These solutions are related to soliton, a combination between a shock and a soliton, and monotonic and oscillatory shock waves. Their behaviors are illustrated and shown graphically. The characteristics of the DA solitary and shock waves are significantly modified by the presence of nonthermal (fast) ions, the ratio of the ion temperature to the electron temperature, and the dust kinematic viscosity. The topology of the phase portrait and the potential diagram of the KdV-Burgers equation is illustrated, whose advantage is the ability to predict different classes of traveling wave solutions according to different phase orbits. The energy of the soliton wave and the electric field are calculated. The results in this paper can be generalized to analyze the nature of plasma waves in both space and laboratory plasma systems.

  18. Woodland Decomposition.

    ERIC Educational Resources Information Center

    Napier, J.

    1988-01-01

    Outlines the role of the main organisms involved in woodland decomposition and discusses some of the variables affecting the rate of nutrient cycling. Suggests practical work that may be of value to high school students either as standard practice or long-term projects. (CW)

  19. Effect of joint damping and joint nonlinearity on the dynamics of space structures

    NASA Technical Reports Server (NTRS)

    Bowden, Mary; Dugundji, John

    1988-01-01

    Analyses of the effect of linear joint characteristics on the vibrations of a free-free, three-joint beam model show that increasing joint damping increases resonant frequencies and increases modal damping but only to the point where the joint gets 'locked up' by damping. This behavior is different from that predicted by modeling joint damping as proportional damping. Nonlinear analyses of the three-joint model with cubic springs at the joints show all the classical single DOF nonlinear response behavior at each resonance of the multiple DOF system: nondoubling of response for a doubling of forcing amplitude, multiple solutions, jump behavior, and resonant frequency shifts. These properties can be concisely quantified by characteristic backbone curves, which show the locus of resonant peaks for increasing forcing amplitude.

  20. Mixing of two collinear Rayleigh waves in an isotropic nonlinear elastic half-space

    SciTech Connect

    Morlock, Merlin B.; Kim, Jin-Yeon; Jacobs, Laurence J.; Qu, Jianmin

    2014-02-18

    Nonlinear mixing of two collinear, initially monochromatic, Rayleigh waves propagating in the same direction in an isotropic, nonlinear elastic solid is investigated analytically. A system of coupled ordinary differential equations is derived based on the Lagrange equations of the second kind to predict the evolution of the higher harmonic and combination frequency components of the fundamentals waves. Numerical results show that the energy transfer is larger for higher frequencies, and that the oscillation of the energy between the different frequency components depends on the amplitudes and frequencies of the fundamental waves. Furthermore, it is illustrated that the horizontal velocity component can form a shock wave while the vertical velocity component can form a pulse. The occurrence of these effects is independent of the specific fundamental frequencies and amplitudes that are mixed. However, the nonlinear interaction is more efficient when the mixing frequencies are close to each other which increases both effects. The analytical model is then extended by implementing diffraction effects in the parabolic approximation.

  1. Bridging Proper Orthogonal Decomposition methods and augmented Newton-Krylov algorithms: an adaptive model order reduction for highly nonlinear mechanical problems

    PubMed Central

    Kerfriden, P.; Gosselet, P.; Adhikari, S.; Bordas, S.

    2013-01-01

    This article describes a bridge between POD-based model order reduction techniques and the classical Newton/Krylov solvers. This bridge is used to derive an efficient algorithm to correct, “on-the-fly”, the reduced order modelling of highly nonlinear problems undergoing strong topological changes. Damage initiation problems are addressed and tackle via a corrected hyperreduction method. It is shown that the relevancy of reduced order model can be significantly improved with reasonable additional costs when using this algorithm, even when strong topological changes are involved. PMID:27076688

  2. Entropy solutions for a nonlinear parabolic problems with lower order term in Orlicz spaces

    NASA Astrophysics Data System (ADS)

    Mabdaoui, M.; Moussa, H.; Rhoudaf, M.

    2016-03-01

    We shall give the proof of existence results for the entropy solutions of the following nonlinear parabolic problem [Equation not available: see fulltext.]where A is a Leray-Lions operator having a growth not necessarily of polynomial type. The lower order term Φ :Ω × (0,T)× {R}→ {R}^N is a Carathéodory function, for a.e. (x,t)in Q_T and for all sin R , satisfying only a growth condition and the right hand side f belongs to L^1(Q_T).

  3. Numerical approximations to nonlinear conservation laws with locally varying time and space grids

    NASA Technical Reports Server (NTRS)

    Osher, S.; Sanders, R.

    1983-01-01

    Numerical approximations to the initial value problem for nonlinear systems of conservation laws are considered. The considered system is said to be hyperbolic when all eigenvalues of every real linear combination of the Jacobian matrices are real. Solutions may develop discontinuities in finite time, even when the initial data are smooth. In the investigation, explicit finite difference methods which use locally varying time grids are considered. The global CFL restriction is replaced by a local restriction. The numerical flux function is studied from a finite volume viewpoint, and a differencing technique is developed at interface points between regions of distinct time increments.

  4. Possible emittance growth induced by nonlinear space charge fields for arbitrary particle distributions

    NASA Astrophysics Data System (ADS)

    Kikuchi, Takashi; Horioka, Kazuhiko

    2016-06-01

    A procedure to obtain a ratio of beam radii at final and initial states in arbitrary particle distributions is proposed, and is applied to the estimation of possible emittance growth for Gaussian and thermal equilibrium distributions. The ratios are estimated for Gaussian and thermal equilibrium distributions as a function of tune depression. The possible emittance growth as a function of tune depression and nonlinear field energy factor is also estimated with and without a constant radius ratio approximation. It is confirmed that the possible emittance growths are almost the same in comparison to the cases with and without the constant radius ratio approximation at each distribution.

  5. Linear and nonlinear analysis of dust acoustic waves in dissipative space dusty plasmas with trapped ions

    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.

  6. A new nonlinear diffusion formalism in a magnetized plasma - Application to space physics and astrophysics

    NASA Technical Reports Server (NTRS)

    Karimbadi, H.; Krauss-Varban, D.

    1992-01-01

    A novel diffusion formalism that takes into account the finite width of resonances is presented. The resonance diagram technique is shown to reproduce the details of the particle orbits very accurately, and can be used to determine the acceleration/scattering in the presence of a given wave spectrum. Ways in which the nonlinear orbits can be incorporated into the diffusion equation are shown. The resulting diffusion equation is an extension of the Q-L theory to cases where the waves have large amplitudes and/or are coherent. This new equation does not have a gap at 90 deg in cases where the individual orbits can cross the gap. The conditions under which the resonance gap at 90-deg pitch angle exits are also examined.

  7. Nonlinear color space and spatiotemporal MRF for hierarchical segmentation of face features in video.

    PubMed

    Liévin, Marc; Luthon, Franck

    2004-01-01

    This paper deals with the low-level joint processing of color and motion for robust face analysis within a feature-based approach. To gain robustness and contrast under unsupervised viewing conditions, a nonlinear color transform relevant for hue segmentation is derived from a logarithmic model. A hierarchical segmentation scheme is based on Markov random field modeling, that combines hue and motion detection within a spatiotemporal neighborhood. Relevant face regions are segmented without parameter tuning. The accuracy of the label fields enables not only face detection and tracking but also geometrical measurements on facial feature edges, such as lips or eyes. Results are shown both on typical test sequences and on various sequences acquired from micro- or mobile cameras. The efficiency of the method makes it suitable for real-time applications aiming at audiovisual communication in unsupervised environments. PMID:15376958

  8. Nonlinear Dust Acoustic Waves in Dissipative Space Dusty Plasmas with Superthermal Electrons and Nonextensive Ions

    NASA Astrophysics Data System (ADS)

    El-Hanbaly, A. M.; El-Shewy, E. K.; Sallah, M.; Darweesh, H. F.

    2016-05-01

    The nonlinear characteristics of the dust acoustic (DA) waves are studied in a homogeneous, collisionless, unmagnetized, and dissipative dusty plasma composed of negatively charged dusty grains, superthermal electrons, and nonextensive ions. Sagdeev pseudopotential technique has been employed to study the large amplitude DA waves. It (Sagdeev pseudopotential) has an evidence for the existence of compressive and rarefractive solitons. The global features of the phase portrait are investigated to understand the possible types of solutions of the Sagdeev form. On the other hand, the reductive perturbation technique has been used to study small amplitude DA waves and yields the Korteweg-de Vries-Burgers (KdV-Burgers) equation that exhibits both soliton and shock waves. The behavior of the obtained results of both large and small amplitude is investigated graphically in terms of the plasma parameters like dust kinematic viscosity, superthermal and nonextensive parameters.

  9. Linear and nonlinear analysis of orbital telescope/space shuttle dynamics and control

    NASA Technical Reports Server (NTRS)

    Roberts, A. S., Jr.; Joshi, S. M.

    1976-01-01

    Work completed on the design and study of an annular suspension and pointing (ASP) system for the space shuttle was presented. This system makes use of a magnetically suspended vernier pointing assembly. The following objectives were pursued in this study: (1) development of a detailed mathematical model of the Space Shuttle/ASP system, (2) design of control laws in order to obtain the desired pointing performance, and (3) prediction of the statistical pointing accuracies in the presence of stochastic disturbances such as crew-motion, and sensor and actuator noise. The first two of these objectives are documented in this report.

  10. Nonlinear satellite wakes in planetary rings. I - Phase-space kinematics

    NASA Technical Reports Server (NTRS)

    Stewart, Glen R.

    1991-01-01

    The explicit expression presently derived for the phase-space density of a planetary ring subjected to perturbations by a proximate satellite recovers the usual perturbed-streamline equations by drawing first-order moments of the phase-space density. The surface density obtained is positive-definite in virtue of taking the ring particles' finite-velocity dispersion into account. The satellite-wake local mean velocity components deviate from the streamline equations' sinusoidal form; this deviation grows as the wake moves downstream from the shepherding satellite.

  11. Elastic-Plastic Nonlinear Response of a Space Shuttle External Tank Stringer. Part 2; Thermal and Mechanical Loadings

    NASA Technical Reports Server (NTRS)

    Knight, Norman F., Jr.; Warren, Jerry E.; Elliott, Kenny B.; Song, Kyongchan; Raju, Ivatury S.

    2012-01-01

    Elastic-plastic, large-deflection nonlinear thermo-mechanical stress analyses are performed for the Space Shuttle external tank s intertank stringers. Detailed threedimensional finite element models are developed and used to investigate the stringer s elastic-plastic response for different thermal and mechanical loading events from assembly through flight. Assembly strains caused by initial installation on an intertank panel are accounted for in the analyses. Thermal loading due to tanking was determined to be the bounding loading event. The cryogenic shrinkage caused by tanking resulted in a rotation of the intertank chord flange towards the center of the intertank, which in turn loaded the intertank stringer feet. The analyses suggest that the strain levels near the first three fasteners remain sufficiently high that a failure may occur. The analyses also confirmed that the installation of radius blocks on the stringer feet ends results in an increase in the stringer capability.

  12. Elastic-Plastic Nonlinear Response of a Space Shuttle External Tank Stringer. Part 1; Stringer-Feet Imperfections and Assembly

    NASA Technical Reports Server (NTRS)

    Knight, Norman F., Jr.; Song, Kyongchan; Elliott, Kenny B.; Raju, Ivatury S.; Warren, Jerry E.

    2012-01-01

    Elastic-plastic, large-deflection nonlinear stress analyses are performed for the external hat-shaped stringers (or stiffeners) on the intertank portion of the Space Shuttle s external tank. These stringers are subjected to assembly strains when the stringers are initially installed on an intertank panel. Four different stringer-feet configurations including the baseline flat-feet, the heels-up, the diving-board, and the toes-up configurations are considered. The assembly procedure is analytically simulated for each of these stringer configurations. The location, size, and amplitude of the strain field associated with the stringer assembly are sensitive to the assumed geometry and assembly procedure. The von Mises stress distributions from these simulations indicate that localized plasticity will develop around the first eight fasteners for each stringer-feet configuration examined. However, only the toes-up configuration resulted in high assembly hoop strains.

  13. Photorefractive nonlinearities caused by the Dember space-charge field in undoped CdTe

    NASA Astrophysics Data System (ADS)

    Schroeder, W. Andreas; Stark, Thomas S.; Boggess, Thomas F.; Smirl, Arthur L.; Valley, George C.

    1991-06-01

    The decay of a photorefractive grating in undoped CdTe produced by subband-gap picosecond excitation has been isolated and temporally resolved. The grating is observed to decay completely at the same rate as the free-carrier grating, a rate consistent with the decay of a 1.7 micron-period grating by ambipolar diffusion. This provides convincing evidence that the source of the photorefractive grating is the Dember space-charge field that is established between electron-hole pairs created by band-to-band two-photon absorption. This establishes the Dember field as a source for the photorefractive space-charge field that can be turned on in the dielectric relaxation time and turned off in the ambipolar diffusion time. It also allows the material to recover fully in the recombination lifetime which can be adjusted by doping or damaging the material.

  14. Nonlinear Attitude Control of Planar Structures in Space Using Only Internal Controls

    NASA Technical Reports Server (NTRS)

    Reyhanoglu, Mahmut; Mcclamroch, N. Harris

    1993-01-01

    An attitude control strategy for maneuvers of an interconnection of planar bodies in space is developed. It is assumed that there are no exogeneous torques and that torques generated by joint motors are used as means of control so that the total angular momentum of the multibody system is a constant, assumed to be zero. The control strategy utilizes the nonintegrability of the expression for the angular momentum. Large angle maneuvers can be designed to achieve an arbitrary reorientation of the multibody system with respect to an inertial frame. The theoretical background for carrying out the required maneuvers is summarized.

  15. Real-time human motion estimation using biomechanical models and non-linear state-space filters.

    PubMed

    Cerveri, P; Rabuffetti, M; Pedotti, A; Ferrigno, G

    2003-03-01

    In the field of sports biomechanics and rehabilitation engineering, the possibility of computing, in real time, the angular displacements and derivatives of human joints, from a video of motion sequences, represents an appealing goal. In particular, applications of biofeedback protocols in rehabilitation can benefit from this capability. The focus of the investigation was concerned with the application of biomechanical models, comprising of a kinematic chain and surface envelopes, and state-space filters, to the computation, in real time and with high accuracy, of the angular data and derivatives. By minimising the distances, measured with TV cameras, between the 2D marker projections and the corresponding back-projected markers located on the mannequin, the configuration of the biomechanical model was automatically updated. The use of state-space estimation allowed the computation of smooth derivatives of the orientation data. Owing to the non-linearity of the functions involved, the derivatives of the observation model were obtained through a multidimensional extension of Stirling's interpolation formula. Proper algorithms were developed to cope with the model calibration, initialisation and data labelling. Extensive experiments on real and simulated motions proved the reliability (maximum angular error less than 1 degree, maximum point reconstruction less than 1 mm) of the developed system, which is robust to false matching caused by marker occlusions. Moreover, orientation artifacts due to skin motion can be reduced by a factor of 50%. PMID:12691430

  16. An Introduction to Twisted Particle Filters and Parameter Estimation in Non-Linear State-Space Models

    NASA Astrophysics Data System (ADS)

    Ala-Luhtala, Juha; Whiteley, Nick; Heine, Kari; Piche, Robert

    2016-09-01

    Twisted particle filters are a class of sequential Monte Carlo methods recently introduced by Whiteley and Lee to improve the efficiency of marginal likelihood estimation in state-space models. The purpose of this article is to extend the twisted particle filtering methodology, establish accessible theoretical results which convey its rationale, and provide a demonstration of its practical performance within particle Markov chain Monte Carlo for estimating static model parameters. We derive twisted particle filters that incorporate systematic or multinomial resampling and information from historical particle states, and a transparent proof which identifies the optimal algorithm for marginal likelihood estimation. We demonstrate how to approximate the optimal algorithm for nonlinear state-space models with Gaussian noise and we apply such approximations to two examples: a range and bearing tracking problem and an indoor positioning problem with Bluetooth signal strength measurements. We demonstrate improvements over standard algorithms in terms of variance of marginal likelihood estimates and Markov chain autocorrelation for given CPU time, and improved tracking performance using estimated parameters.

  17. Tunable High-Intensity Electron Bunch Train Production Based on Nonlinear Longitudinal Space Charge Oscillation.

    PubMed

    Zhang, Zhen; Yan, Lixin; Du, Yingchao; Zhou, Zheng; Su, Xiaolu; Zheng, Lianmin; Wang, Dong; Tian, Qili; Wang, Wei; Shi, Jiaru; Chen, Huaibi; Huang, Wenhui; Gai, Wei; Tang, Chuanxiang

    2016-05-01

    High-intensity trains of electron bunches with tunable picosecond spacing are produced and measured experimentally with the goal of generating terahertz (THz) radiation. By imposing an initial density modulation on a relativistic electron beam and controlling the charge density over the beam propagation, density spikes of several-hundred-ampere peak current in the temporal profile, which are several times higher than the initial amplitudes, have been observed for the first time. We also demonstrate that the periodic spacing of the bunch train can be varied continuously either by tuning launching phase of a radio-frequency gun or by tuning the compression of a downstream magnetic chicane. Narrow-band coherent THz radiation from the bunch train was also measured with μJ-level energies and tunable central frequency of the spectrum in the range of ∼0.5 to 1.6 THz. Our results pave the way towards generating mJ-level narrow-band coherent THz radiation and driving high-gradient wakefield-based acceleration. PMID:27203327

  18. Tunable High-Intensity Electron Bunch Train Production Based on Nonlinear Longitudinal Space Charge Oscillation

    NASA Astrophysics Data System (ADS)

    Zhang, Zhen; Yan, Lixin; Du, Yingchao; Zhou, Zheng; Su, Xiaolu; Zheng, Lianmin; Wang, Dong; Tian, Qili; Wang, Wei; Shi, Jiaru; Chen, Huaibi; Huang, Wenhui; Gai, Wei; Tang, Chuanxiang

    2016-05-01

    High-intensity trains of electron bunches with tunable picosecond spacing are produced and measured experimentally with the goal of generating terahertz (THz) radiation. By imposing an initial density modulation on a relativistic electron beam and controlling the charge density over the beam propagation, density spikes of several-hundred-ampere peak current in the temporal profile, which are several times higher than the initial amplitudes, have been observed for the first time. We also demonstrate that the periodic spacing of the bunch train can be varied continuously either by tuning launching phase of a radio-frequency gun or by tuning the compression of a downstream magnetic chicane. Narrow-band coherent THz radiation from the bunch train was also measured with μ J -level energies and tunable central frequency of the spectrum in the range of ˜0.5 to 1.6 THz. Our results pave the way towards generating mJ-level narrow-band coherent THz radiation and driving high-gradient wakefield-based acceleration.

  19. The Global Nonlinear Stability of Minkowski Space for Self-gravitating Massive Fields - The Wave-Klein-Gordon Model

    NASA Astrophysics Data System (ADS)

    LeFloch, Philippe G.; Ma, Yue

    2016-01-01

    The Hyperboloidal Foliation Method (introduced by the authors in 2014) is extended here and applied to the Einstein equations of general relativity. Specifically, we establish the nonlinear stability of Minkowski spacetime for self-gravitating massive scalar fields, while existing methods only apply to massless scalar fields. First of all, by analyzing the structure of the Einstein equations in wave coordinates, we exhibit a nonlinear wave-Klein-Gordon model defined on a curved background, which is the focus of the present paper. For this model, we prove here the existence of global-in-time solutions to the Cauchy problem, when the initial data have sufficiently small Sobolev norms. A major difficulty comes from the fact that the class of conformal Killing fields of Minkowski space is significantly reduced in the presence of a massive scalar field, since the scaling vector field is not conformal Killing for the Klein-Gordon operator. Our method relies on the foliation (of the interior of the light cone) of Minkowski spacetime by hyperboloidal hypersurfaces and uses Lorentz-invariant energy norms. We introduce a frame of vector fields adapted to the hyperboloidal foliation and we establish several key properties: Sobolev and Hardy-type inequalities on hyperboloids, as well as sup-norm estimates, which correspond to the sharp time decay for the wave and the Klein-Gordon equations. These estimates allow us to control interaction terms associated with the curved geometry and the massive field by distinguishing between two levels of regularity and energy growth and by a successive use of our key estimates in order to close a bootstrap argument.

  20. Nonlinear state-space modeling of human motion using 2-D marker observations.

    PubMed

    Vartiainen, Paavo; Bragge, Timo; Arokoski, Jari P; Karjalainen, Pasi A

    2014-07-01

    A novel method for the estimation of human kinematics, based on state-space modeling, is proposed. The state consists of the positions, orientations, velocities, and accelerations of an articulated model. Estimation is performed using the unscented Kalman filter (UKF) algorithm with a fixed-interval smoother. Impulsive acceleration at floor contact of the foot is estimated by implementing a contact constraint in the UKF evolution model. The constraint inserts an acceleration impulse into the model state. The estimation method was applied to marker-based motion analysis in a motion laboratory. Validation measurements were performed with a rigid test device and with human gait. A triaxial accelerometer was used to evaluate acceleration estimates. Comparison between the proposed method and the extended Kalman smoother showed a clear difference in the quality of estimates during impulsive accelerations. The proposed approach enables estimation of human kinematics during both continuous and transient accelerations. The approach provides a novel way of estimating acceleration at foot initial contact, and thus enables more accurate evaluation of loading from the beginning of the floor contact. PMID:24760898

  1. Mode decomposition as a methodology for developing convective-scale representations in global models

    NASA Astrophysics Data System (ADS)

    Yano, Jun-Ichi; Redelsperger, Jean-Luc; Bechtold, Peter; Guichard, Françoise

    2005-07-01

    Mode decomposition is proposed as a methodology for developing subgrid-scale physical representations in global models by a systematic reduction of an originally full system such as a cloud-resolving model (CRM). A general formulation is presented, and also discussed are mathematical requirements that make this procedure possible. Features of this general methodology are further elucidated by the two specific examples: mass fluxes and wavelets.The traditional mass-flux formulation for convective parametrizations is derived as a special case from this general formulation. It is based on the decomposition of a horizontal domain into an approximate sum of piecewise-constant segments. Thus, a decomposition of CRM outputs on this basis is crucial for their direct verification. However, this decomposition is mathematically not well-posed nor unique due to the lack of admissibility. A classification into cloud types, primarily based on precipitation characteristics of the atmospheric columns, that has been used as its substitute, does not necessarily provide a good approximation for a piecewiseconstant segment decomposition. This difficulty with mass-flux decomposition makes a verification of the formulational details of parametrizations based on mass fluxes by a CRM inherently difficult.The wavelet decomposition is an alternative possibility that can more systematically decompose the convective system. Its completeness and orthogonality also allow a prognostic description of a CRM system in wavelet space in the same manner as is done in Fourier space. The wavelets can, furthermore, efficiently represent the various convective coherencies by a limited number of modes due to their spatial localizations. Thus, the degree of complexity of the wavelet-based prognostic representation of a CRM can be extensively reduced. Such an extensive reduction may allow its use in place of current cumulus parametrizations. This wavelet-based scheme can easily be verified from the full

  2. Uncertainty propagation in orbital mechanics via tensor decomposition

    NASA Astrophysics Data System (ADS)

    Sun, Yifei; Kumar, Mrinal

    2016-03-01

    Uncertainty forecasting in orbital mechanics is an essential but difficult task, primarily because the underlying Fokker-Planck equation (FPE) is defined on a relatively high dimensional (6-D) state-space and is driven by the nonlinear perturbed Keplerian dynamics. In addition, an enormously large solution domain is required for numerical solution of this FPE (e.g. encompassing the entire orbit in the x-y-z subspace), of which the state probability density function (pdf) occupies a tiny fraction at any given time. This coupling of large size, high dimensionality and nonlinearity makes for a formidable computational task, and has caused the FPE for orbital uncertainty propagation to remain an unsolved problem. To the best of the authors' knowledge, this paper presents the first successful direct solution of the FPE for perturbed Keplerian mechanics. To tackle the dimensionality issue, the time-varying state pdf is approximated in the CANDECOMP/PARAFAC decomposition tensor form where all the six spatial dimensions as well as the time dimension are separated from one other. The pdf approximation for all times is obtained simultaneously via the alternating least squares algorithm. Chebyshev spectral differentiation is employed for discretization on account of its spectral ("super-fast") convergence rate. To facilitate the tensor decomposition and control the solution domain size, system dynamics is expressed using spherical coordinates in a noninertial reference frame. Numerical results obtained on a regular personal computer are compared with Monte Carlo simulations.

  3. Rogue waves in electronegative space plasmas: The link between the family of the KdV equations and the nonlinear Schrödinger equation

    NASA Astrophysics Data System (ADS)

    El-Tantawy, S. A.

    2016-05-01

    We examine the likelihood of the ion-acoustic rogue waves propagation in a non-Maxwellian electronegative plasma in the framework of the family of the Korteweg-de Vries (KdV) equations (KdV/modified KdV/Extended KdV equation). For this purpose, we use the reductive perturbation technique to carry out this study. It is known that the family of the KdV equations have solutions of distinct structures such as solitons, shocks, kinks, cnoidal waves, etc. However, the dynamics of the nonlinear rogue waves is governed by the nonlinear Schrödinger equation (NLSE). Thus, the family of the KdV equations is transformed to their corresponding NLSE developing a weakly nonlinear wave packets. We show the possible region for the existence of the rogue waves and define it precisely for typical parameters of space plasmas. We investigate numerically the effects of relevant physical parameters, namely, the negative ion relative concentration, the nonthermal parameter, and the mass ratio on the propagation of the rogue waves profile. The present study should be helpful in understanding the salient features of the nonlinear structures such as, ion-acoustic solitary waves, shock waves, and rogue waves in space and in laboratory plasma where two distinct groups of ions, i.e. positive and negative ions, and non-Maxwellian (nonthermal) electrons are present.

  4. Study of nonlinear electron-acoustic solitary and shock waves in a dissipative, nonplanar space plasma with superthermal hot electrons

    SciTech Connect

    Han, Jiu-Ning He, Yong-Lin; Luo, Jun-Hua; Nan, Ya-Gong; Han, Zhen-Hai; Dong, Guang-Xing; Duan, Wen-Shan; Li, Jun-Xiu

    2014-01-15

    With the consideration of the superthermal electron distribution, we present a theoretical investigation about the nonlinear propagation of electron-acoustic solitary and shock waves in a dissipative, nonplanar non-Maxwellian plasma comprised of cold electrons, superthermal hot electrons, and stationary ions. The reductive perturbation technique is used to obtain a modified Korteweg-de Vries Burgers equation for nonlinear waves in this plasma. We discuss the effects of various plasma parameters on the time evolution of nonplanar solitary waves, the profile of shock waves, and the nonlinear structure induced by the collision between planar solitary waves. It is found that these parameters have significant effects on the properties of nonlinear waves and collision-induced nonlinear structure.

  5. Morphological Decomposition in Reading Hebrew Homographs

    ERIC Educational Resources Information Center

    Miller, Paul; Liran-Hazan, Batel; Vaknin, Vered

    2016-01-01

    The present work investigates whether and how morphological decomposition processes bias the reading of Hebrew heterophonic homographs, i.e., unique orthographic patterns that are associated with two separate phonological, semantic entities depicted by means of two morphological structures (linear and nonlinear). In order to reveal the nature of…

  6. Phase-Space Reconstruction: a Path Towards the Next Generation of Nonlinear Differential Equation Based Models and Its Implications Towards Non-Uniform Sampling Theory

    SciTech Connect

    Charles R. Tolle; Mark Pengitore

    2009-08-01

    This paper explores the overlaps between the Control community’s work on System Identification (SysID) and the Physics, Mathematics, Chaos, and Complexity communities’ work on phase-space reconstruction via time-delay embedding. There are numerous overlaps between the goals of each community. Nevertheless, the Controls community can gain new insight as well as some new very powerful tools for SysID from the latest developments within the Physics, Mathematics, Chaos, and Complexity communities. These insights are gained via the work on phase-space reconstruction of non-linear dynamics. New methods for discovering non-linear differential based equations that evolved from embedding operations can shed new light on hybrid-systems theory, Nyquest-Shannon’s Theories, and network based control theory. This paper strives to guide the Controls community towards a closer inspection of the tools and additional insights being developed within the Physics, Mathematics, Chaos, and Complexity communities for discovery of system dynamics, the first step in control system development. The paper introduces the concepts of phase-space reconstruction via time-delay embedding (made famous byWhitney, Takens, and Sauer’s Thoreoms), intergrate-and-fire embedding, and non-linear differential equation discovery based on Perona’s method.

  7. Fourier-Space Nonlinear Rayleigh-Taylor Growth Measurements of 3D Laser-Imprinted Modulations in Planar Targets

    SciTech Connect

    Smalyuk, V.A.; Sadot, O.; Delettrez, J.A.; Meyerhofer, D.D.; Regan, S.P.; Sangster, T.C.

    2005-12-05

    Nonlinear growth of 3-D broadband nonuniformities was measured near saturation levels using x-ray radiography in planar foils accelerated by laser light. The initial target modulations were seeded by laser nonuniformities and later amplified during acceleration by Rayleigh-Taylor instability. The nonlinear saturation velocities are measured for the first time and are found to be in excellent agreement with Haan predictions. The measured growth of long-wavelength modes is consistent with enhanced, nonlinear, long-wavelength generation in ablatively driven targets.

  8. Non-linearity's in proxy space: Three methods to deal with the non-linear behavior of proxies in calcareous marine skeletons

    NASA Astrophysics Data System (ADS)

    Bauwens, M.; Ohlsson, H.; Beelaerts, V.; Barbé, K.; Dehairs, F.; Schoukens, J.

    2009-04-01

    The reconstruction of seasonal variations in the paleo-environement is possible thanks to the bivalve shells that are commonly found back in archeological sites and that are sensitive environmental recorders. To do these reconstructions we look to the chemical composition of a shell along his growth axis. When a certain element or isotope can be related to an environmental parameter it is called a proxy for that parameter. Many elemental and isotopic have been proposed as potential temperature proxy. But the same problem is showing up over and over again: a certain elemental shows a good linear correlation with temperature, but at the moment that the intrinsic variation has to be explained, the incorporation of the proxy seems to be much more complex than assumed in the first instance. Two observations gave us the idea to use a new type of models to reconstruct temperature. In first instance we observed that all proxy incorporations were always influenced by more than one environmental parameter, which automatically implicate that the reconstruction of an environmental parameter should be done with more than one proxy input. A second observation is that some proxy-environment relations seemed to be non-linear. We propose to do climate reconstructions based on Non-linear multi-proxy models. We will present our own intuitive approach to do temperature reconstructions and we will compare our results to two recognized engineering techniques: Manifold learning and Supported vector machines. These three methods are validated on shell data of Scheldt estuarine environment using high resolution measurements of Mg, Mn, Ba, Sr and Pb on a LA-ICP-MS. The site specificity is tested for all methods using shell-data from 4 different sites along the Scheldt. Time specificity is checked by reconstructing the temperature corresponding to a shell that grew in another year. The best reconstructions are computed with the manifold learning algorithms, but the simplicity and high

  9. Cosmic flows and the expansion of the local Universe from non-linear phase-space reconstructions

    NASA Astrophysics Data System (ADS)

    Heß, Steffen; Kitaura, Francisco-Shu

    2016-03-01

    In this work, we investigate the impact of cosmic flows and density perturbations on Hubble constant H0 measurements using non-linear phase-space reconstructions of the Local Universe (LU). In particular, we rely on a set of 25 precise constrained N-body simulations based on Bayesian initial conditions reconstructions of the LU using the Two-Micron Redshift Survey galaxy sample within distances of about 90 h-1 Mpc. These have been randomly extended up to volumes enclosing distances of 360 h-1 Mpc with augmented Lagrangian perturbation theory (750 simulations in total), accounting in this way for gravitational mode coupling from larger scales, correcting for periodic boundary effects, and estimating systematics of missing attractors (σlarge = 134 s-1 km). We report on Local Group (LG) speed reconstructions, which for the first time are compatible with those derived from cosmic microwave background-dipole measurements: |vLG| = 685 ± 137 s-1 km. The direction (l, b) = (260.5° ± 13.3°, 39.1 ± 10.4°) is found to be compatible with the observations after considering the variance of large scales. Considering this effect of large scales, our local bulk flow estimations assuming a Λ cold dark matter model are compatible with the most recent estimates based on velocity data derived from the Tully-Fisher relation. We focus on low-redshift supernova measurements out to 0.01 < z < 0.025, which have been found to disagree with probes at larger distances. Our analysis indicates that there are two effects related to cosmic variance contributing to this tension. The first one is caused by the anisotropic distribution of supernovae, which aligns with the velocity dipole and hence induces a systematic boost in H0. The second one is due to the inhomogeneous matter fluctuations in the LU. In particular, a divergent region surrounding the Virgo Supercluster is responsible for an additional positive bias in H0. Taking these effects into account yields a correction of ΔH0 = -1

  10. Identification of the nonlinear state-space dynamics of the action-perception cycle for visually induced postural sway.

    PubMed

    Giese, M A; Dijkstra, T M; Schöner, G; Gielen, C C

    1996-05-01

    Human subjects standing in a sinusoidally moving visual environment display postural sway with characteristic dynamical properties. We analyzed the spatiotemporal properties of this sway in an experiment in which the frequency of the visual motion was varied. We found a constant gain near 1, which implies that the sway motion matches the spatial parameters of the visual motion for a large range of frequencies. A linear dynamical model with constant parameters was compared quantitatively with the data. Its failure to describe correctly the spatiotemporal properties of the system led us to consider adaptive and nonlinear models. To differentiate between possible alternative structures we directly fitted nonlinear differential equations to the sway and visual motion trajectories on a trial-by-trial basis. We found that the eigenfrequency of the fitted model adapts strongly to the visual motion frequency. The damping coefficient decreases with increasing frequency. This indicates that the system destabilizes its postural state in the inertial frame. This leads to a faster internal dynamics which is capable of synchronizing posture with fast-moving visual environments. Using an algorithm which allows the identification of essentially nonlinear terms of the dynamics we found small nonlinear contributions. These nonlinearities are not consistent with a limit-cycle dynamics, accounting for the robustness of the amplitude of postural sway against frequency variations. We interpret out results in terms of active generation of postural sway specified by sensory information. We derive also a number of conclusions for a behavior-oriented analysis of the postural system. PMID:8991458

  11. Nested Taylor decomposition in multivariate function decomposition

    NASA Astrophysics Data System (ADS)

    Baykara, N. A.; Gürvit, Ercan

    2014-12-01

    Fluctuationlessness approximation applied to the remainder term of a Taylor decomposition expressed in integral form is already used in many articles. Some forms of multi-point Taylor expansion also are considered in some articles. This work is somehow a combination these where the Taylor decomposition of a function is taken where the remainder is expressed in integral form. Then the integrand is decomposed to Taylor again, not necessarily around the same point as the first decomposition and a second remainder is obtained. After taking into consideration the necessary change of variables and converting the integration limits to the universal [0;1] interval a multiple integration system formed by a multivariate function is formed. Then it is intended to apply the Fluctuationlessness approximation to each of these integrals one by one and get better results as compared with the single node Taylor decomposition on which the Fluctuationlessness is applied.

  12. Effect of motor dynamics on nonlinear feedback robot arm control

    NASA Technical Reports Server (NTRS)

    Tarn, Tzyh-Jong; Li, Zuofeng; Bejczy, Antal K.; Yun, Xiaoping

    1991-01-01

    A nonlinear feedback robot controller that incorporates the robot manipulator dynamics and the robot joint motor dynamics is proposed. The manipulator dynamics and the motor dynamics are coupled to obtain a third-order-dynamic model, and differential geometric control theory is applied to produce a linearized and decoupled robot controller. The derived robot controller operates in the robot task space, thus eliminating the need for decomposition of motion commands into robot joint space commands. Computer simulations are performed to verify the feasibility of the proposed robot controller. The controller is further experimentally evaluated on the PUMA 560 robot arm. The experiments show that the proposed controller produces good trajectory tracking performances and is robust in the presence of model inaccuracies. Compared with a nonlinear feedback robot controller based on the manipulator dynamics only, the proposed robot controller yields conspicuously improved performance.

  13. Adaptive wavelet collocation methods for initial value boundary problems of nonlinear PDE's

    NASA Technical Reports Server (NTRS)

    Cai, Wei; Wang, Jian-Zhong

    1993-01-01

    We have designed a cubic spline wavelet decomposition for the Sobolev space H(sup 2)(sub 0)(I) where I is a bounded interval. Based on a special 'point-wise orthogonality' of the wavelet basis functions, a fast Discrete Wavelet Transform (DWT) is constructed. This DWT transform will map discrete samples of a function to its wavelet expansion coefficients in O(N log N) operations. Using this transform, we propose a collocation method for the initial value boundary problem of nonlinear PDE's. Then, we test the efficiency of the DWT transform and apply the collocation method to solve linear and nonlinear PDE's.

  14. Nonlinear Analysis of the Space Shuttle Superlightweight LO2 Tank. Part 2; Behavior Under 3g End-of-Flight Loads

    NASA Technical Reports Server (NTRS)

    Nemeth, Michael P.; Young, Richard D.; Collins, Timothy J.; Starnes, James H.,Jr.

    1998-01-01

    Results of linear bifurcation and nonlinear analyses of the Space Shuttle super lightweight (SLWT) external liquid-oxygen (LO2) tank are presented for an important end-of-flight loading condition. These results illustrate an important type of response mode for thin-walled shells, that are subjected to combined mechanical and thermal loads, that may be encountered in the design of other liquid-fuel launch vehicles. Linear bifurcation analyses are presented that predict several nearly equal eigenvalues that correspond to local buckling modes in the aft dome of the LO2 tank. In contrast, the nonlinear response phenomenon is shown to consist of a short-wavelength bending deformation in the aft elliptical dome of the LO2 tank that grows in amplitude in a stable manner with increasing load. Imperfection sensitivity analyses are presented that show that the presence of several nearly equal eigenvalues does not lead to a premature general instability mode for the aft dome. For the linear bifurcation and nonlinear analyses, the results show that accurate predictions of the response of the shell generally require a large-scale, high fidelity finite-element model. Results are also presented that show that the SLWT LO2 tank can support loads in excess of approximately 1.9 times the values of the operational loads considered.

  15. Diagnosis of nonlinear systems using time series analysis

    SciTech Connect

    Hunter, N.F. Jr.

    1991-01-01

    Diagnosis and analysis techniques for linear systems have been developed and refined to a high degree of precision. In contrast, techniques for the analysis of data from nonlinear systems are in the early stages of development. This paper describes a time series technique for the analysis of data from nonlinear systems. The input and response time series resulting from excitation of the nonlinear system are embedded in a state space. The form of the embedding is optimized using local canonical variate analysis and singular value decomposition techniques. From the state space model, future system responses are estimated. The expected degree of predictability of the system is investigated using the state transition matrix. The degree of nonlinearity present is quantified using the geometry of the transfer function poles in the z plane. Examples of application to a linear single-degree-of-freedom system, a single-degree-of-freedom Duffing Oscillator, and linear and nonlinear three degree of freedom oscillators are presented. 11 refs., 9 figs.

  16. Nonlinear mapping methods with adjustable computational complexity for hyperspectral image analysis

    NASA Astrophysics Data System (ADS)

    Myasnikov, E. V.

    2015-12-01

    Nonlinear mapping (Sammon mapping) is a well-known dimensionality reduction technique. Recently several nonlinear mapping methods with reduced computational complexity have been proposed but they do not provide a flexible control over a computational complexity. In this paper a nonlinear mapping method with adjustable computational complexity is proposed. The proposed method is based on the hierarchical decomposition of the multidimensional space, priority queues, and simple optimization procedure to provide fast and flexible dimensionality reduction process. The proposed method is compared to an alternative one based on stochastic optimization. The experiments are carried out on well-known hyperspectral images. Studied methods are evaluated in terms of the data mapping error and runtime. Experimental results for both two- and three-dimensional output spaces are presented.

  17. Analysis of pressure perturbation sources on a generic space launcher after-body in supersonic flow using zonal turbulence modeling and dynamic mode decomposition

    NASA Astrophysics Data System (ADS)

    Statnikov, Vladimir; Sayadi, Taraneh; Meinke, Matthias; Schmid, Peter; Schröder, Wolfgang

    2015-01-01

    A sparsity promoting dynamic mode decomposition (DMD) combined with a classical data-based statistical analysis is applied to the turbulent wake of a generic axisymmetric configuration of an Ariane 5-like launcher at Ma∞ = 6.0 computed via a zonal Reynolds-averaged Navier-Stokes/large-eddy simulation (RANS/LES) method. The objective of this work is to gain a better understanding of the wake flow dynamics of the generic launcher by clarification and visualization of initially unknown pressure perturbation sources on its after-body in coherent flow patterns. The investigated wake topology is characterized by a subsonic cavity region around the cylindrical nozzle extension which is formed due to the displacement effect of the afterexpanding jet plume emanating from the rocket nozzle (Mae = 2.52, pe/p∞ = 100) and the shear layer shedding from the main body. The cavity region contains two toroidal counter-rotating large-scale vortices which extensively interact with the turbulent shear layer, jet plume, and rocket walls, leading to the shear layer instability process to be amplified. The induced velocity fluctuations in the wake and the ultimately resulting pressure perturbations on the after-body feature three global characteristic frequency ranges, depending on the streamwise position inside the cavity. The most dominant peaks are detected at SrD r3 = 0.85 ± 0.075 near the nozzle exit, while the lower frequency peaks, in the range of SrD r2 = 0.55 ± 0.05 and SrD r1 = 0.25 ± 0.05, are found to be dominant closer to the rocket's base. A sparse promoting DMD algorithm is applied to the time-resolved velocity field to clarify the origin of the detected peaks. This analysis extracts three low-frequency spatial modes at SrD = 0.27, 0.56, and 0.85. From the three-dimensional shape of the DMD modes and the reconstructed modulation of the mean flow in time, it is deduced that the detected most dominant peaks of SrD r3 ≈ 0.85 are caused by the radial flapping motion of

  18. Conserved charges for black holes in Einstein-Gauss-Bonnet gravity coupled to nonlinear electrodynamics in AdS space

    NASA Astrophysics Data System (ADS)

    Mišković, Olivera; Olea, Rodrigo

    2011-01-01

    Motivated by possible applications within the framework of anti-de Sitter gravity/conformal field theory correspondence, charged black holes with AdS asymptotics, which are solutions to Einstein-Gauss-Bonnet gravity in D dimensions, and whose electric field is described by nonlinear electrodynamics are studied. For a topological static black hole ansatz, the field equations are exactly solved in terms of the electromagnetic stress tensor for an arbitrary nonlinear electrodynamic Lagrangian in any dimension D and for arbitrary positive values of Gauss-Bonnet coupling. In particular, this procedure reproduces the black hole metric in Born-Infeld and conformally invariant electrodynamics previously found in the literature. Altogether, it extends to D>4 the four-dimensional solution obtained by Soleng in logarithmic electrodynamics, which comes from vacuum polarization effects. Falloff conditions for the electromagnetic field that ensure the finiteness of the electric charge are also discussed. The black hole mass and vacuum energy as conserved quantities associated to an asymptotic timelike Killing vector are computed using a background-independent regularization of the gravitational action based on the addition of counterterms which are a given polynomial in the intrinsic and extrinsic curvatures.

  19. Conserved charges for black holes in Einstein-Gauss-Bonnet gravity coupled to nonlinear electrodynamics in AdS space

    SciTech Connect

    Miskovic, Olivera; Olea, Rodrigo

    2011-01-15

    Motivated by possible applications within the framework of anti-de Sitter gravity/conformal field theory correspondence, charged black holes with AdS asymptotics, which are solutions to Einstein-Gauss-Bonnet gravity in D dimensions, and whose electric field is described by nonlinear electrodynamics are studied. For a topological static black hole ansatz, the field equations are exactly solved in terms of the electromagnetic stress tensor for an arbitrary nonlinear electrodynamic Lagrangian in any dimension D and for arbitrary positive values of Gauss-Bonnet coupling. In particular, this procedure reproduces the black hole metric in Born-Infeld and conformally invariant electrodynamics previously found in the literature. Altogether, it extends to D>4 the four-dimensional solution obtained by Soleng in logarithmic electrodynamics, which comes from vacuum polarization effects. Falloff conditions for the electromagnetic field that ensure the finiteness of the electric charge are also discussed. The black hole mass and vacuum energy as conserved quantities associated to an asymptotic timelike Killing vector are computed using a background-independent regularization of the gravitational action based on the addition of counterterms which are a given polynomial in the intrinsic and extrinsic curvatures.

  20. Reconstruction of principal dynamical modes from climatic variability: nonlinear approach

    NASA Astrophysics Data System (ADS)

    Mukhin, Dmitry; Gavrilov, Andrey; Loskutov, Evgeny; Feigin, Alexander; Kurths, Juergen

    2015-04-01

    Analysis of multivariate time-series produced by complex systems requires efficient tools for reduction of data dimension. We consider this problem in relation to empirical modeling of climate, which implies an analysis of spatial-distributed time-series. The main goal is to establish the number of principal modes which have key contribution to data and actually governs the observed variability. Currently, the number of widely used linear methods based on PCA and factor analysis exists, which yield different data decompositions taking into consideration simultanious/time-lag correlations between spatial grid points. However, the question about possibility of improving the decomposition by taking into account nonlinear couplings between variables often remains untouched. In the report the method for constructing principal dynamic modes on the basis of low-dimensional nonlinear parametric representation of observed multivariate time-series is suggested. It is aimed to extracting the set of latent modes that both explains an essential part of variability, and obeys the simplest evolution law. Thus, this approach can be used for optimal reconstruction of the phase space for empirical prognostic modeling of observed dynamics. The evidence of nonlinear couplings in SST space-distributed data covering the Globe is investigated by the proposed approach. It is demonstrated that the obtained principal modes capture more part of SST variability than principal components (PCs) constructed by either EOF decomposition or its spatio-temporal extension. Relation of these modes to various climate phenomena is shown and discussed in the report. The application of the approach to data-driven forecast of climate bahavior is also discussed.

  1. The influence of crystallinity degree on the glycine decomposition induced by 1 MeV proton bombardment in space analog conditions.

    PubMed

    Pilling, Sergio; Mendes, Luiz A V; Bordalo, Vinicius; Guaman, Christian F M; Ponciano, Cássia R; da Silveira, Enio F

    2013-01-01

    Glycine is the simplest proteinaceous amino acid and is present in all life-forms on Earth. In aqueous solutions, it appears mainly as zwitterion glycine (+NH3CH2COO-); however, in solid phase, it may be found in amorphous or crystalline (α, β, and γ) forms. The crystalline forms differ from each other by the packing of zwitterions in the unitary cells and by the number of intermolecular hydrogen bonds. This molecular species has been extensively detected in carbonaceous meteorites and was recently observed in the cometary samples returned to Earth by NASA's Stardust spacecraft. In space, glycine is exposed to several radiation fields at different temperatures. We present an experimental study on the destruction of zwitterionic glycine crystals at room temperature by 1 MeV protons, in which the dependence of the destruction rates of the α-glycine and β-glycine crystals on bombardment fluence is investigated. The samples were analyzed in situ by Fourier transform infrared spectrometry at different proton fluences. The experiments occurred under ultrahigh vacuum conditions at the Van de Graaff accelerator lab at the Pontifical Catholic University at Rio de Janeiro (PUC-Rio), Brazil. For low fluences, the dissociation cross section of α-glycine was observed to be 2.5×10(-14) cm2, a value roughly 5 times higher than the dissociation cross section found for β-glycine. The estimated half-lives of α-glycine and β-glycine zwitterionic forms extrapolated to the Earth orbit environment are 9×10(5) and 4×10(6) years, respectively. In the diffuse interstellar medium the estimated values are 1 order of magnitude lower. These results suggest that pristine interstellar β-glycine is the one most likely to survive the hostile environments of space radiation. A small feature around 1650-1700 cm(-1), tentatively attributed to an amide functional group, was observed in the IR spectra of irradiated samples, suggesting that cosmic rays may induce peptide bond synthesis in

  2. Hydrological response to earthquakes in the Haibara well, central Japan - I. Groundwater level changes revealed using state space decomposition of atmospheric pressure, rainfall and tidal responses

    USGS Publications Warehouse

    Matsumoto, N.; Kitagawa, G.; Roeloffs, E.A.

    2003-01-01

    For the groundwater level observed at the Haibara well, Shizuoka Prefecture, central Japan, time series analysis using state-space modelling is applied to extract hydrological anomalies related to earthquakes. This method can decompose observed groundwater level time series into five components: atmospheric pressure, tidal, and precipitation responses, observation noise, and residual water level. The decomposed responses to atmospheric pressure and precipitation are independently determined and are consistent with the expected response to surface loading. In the groundwater level at the Haibara well, 28 coseismic changes can be discerned during the period from 1981 April to 1997 December. There is a threshold in the relationship between earthquake magnitude and the well-hypocentre distance, above which earthquakes cause coseismic changes in the residual water level. All of the coseismic water level changes at the Haibara well are decreases, although 33 per cent of the estimated coseismic volumetric strain steps are contraction, which would be expected to cause water level increases. The coseismic changes in groundwater level are more closely proportional to the estimated ground motion than to coseismic volumetric strain steps, suggesting that ground motion due to earthquakes is the major cause of the coseismic water level drops and that the contribution from static strain is rather small. Possible pre- or inter-earthquake water level changes have occurred at the Haibara well and may have been caused by local aseismic crustal deformation.

  3. A decomposition of irreversible diffusion processes without detailed balance

    NASA Astrophysics Data System (ADS)

    Qian, Hong

    2013-05-01

    As a generalization of deterministic, nonlinear conservative dynamical systems, a notion of canonical conservative dynamics with respect to a positive, differentiable stationary density ρ(x) is introduced: dot{x}=j(x) in which ∇.(ρ(x)j(x)) = 0. Such systems have a conserved "generalized free energy function" F[u] = ∫u(x, t)ln (u(x, t)/ρ(x))dx in phase space with a density flow u(x, t) satisfying ∂ut = -∇.(ju). Any general stochastic diffusion process without detailed balance, in terms of its Fokker-Planck equation, can be decomposed into a reversible diffusion process with detailed balance and a canonical conservative dynamics. This decomposition can be rigorously established in a function space with inner product defined as ⟨ϕ, ψ⟩ = ∫ρ-1(x)ϕ(x)ψ(x)dx. Furthermore, a law for balancing F[u] can be obtained: The non-positive dF[u(x, t)]/dt = Ein(t) - ep(t) where the "source" Ein(t) ⩾ 0 and the "sink" ep(t) ⩾ 0 are known as house-keeping heat and entropy production, respectively. A reversible diffusion has Ein(t) = 0. For a linear (Ornstein-Uhlenbeck) diffusion process, our decomposition is equivalent to the previous approaches developed by Graham and Ao, as well as the theory of large deviations. In terms of two different formulations of time reversal for a same stochastic process, the meanings of dissipative and conservative stationary dynamics are discussed.

  4. Analytical gradients of complete active space self-consistent field energies using Cholesky decomposition: Geometry optimization and spin-state energetics of a ruthenium nitrosyl complex

    SciTech Connect

    Delcey, Mickaël G.; Freitag, Leon; González, Leticia; Pedersen, Thomas Bondo; Aquilante, Francesco; Lindh, Roland

    2014-05-07

    We present a formulation of analytical energy gradients at the complete active space self-consistent field (CASSCF) level of theory employing density fitting (DF) techniques to enable efficient geometry optimizations of large systems. As an example, the ground and lowest triplet state geometries of a ruthenium nitrosyl complex are computed at the DF-CASSCF level of theory and compared with structures obtained from density functional theory (DFT) using the B3LYP, BP86, and M06L functionals. The average deviation of all bond lengths compared to the crystal structure is 0.042 Å at the DF-CASSCF level of theory, which is slightly larger but still comparable with the deviations obtained by the tested DFT functionals, e.g., 0.032 Å with M06L. Specifically, the root-mean-square deviation between the DF-CASSCF and best DFT coordinates, delivered by BP86, is only 0.08 Å for S{sub 0} and 0.11 Å for T{sub 1}, indicating that the geometries are very similar. While keeping the mean energy gradient errors below 0.25%, the DF technique results in a 13-fold speedup compared to the conventional CASSCF geometry optimization algorithm. Additionally, we assess the singlet-triplet energy vertical and adiabatic differences with multiconfigurational second-order perturbation theory (CASPT2) using the DF-CASSCF and DFT optimized geometries. It is found that the vertical CASPT2 energies are relatively similar regardless of the geometry employed whereas the adiabatic singlet-triplet gaps are more sensitive to the chosen triplet geometry.

  5. Analytical gradients of complete active space self-consistent field energies using Cholesky decomposition: geometry optimization and spin-state energetics of a ruthenium nitrosyl complex.

    PubMed

    Delcey, Mickaël G; Freitag, Leon; Pedersen, Thomas Bondo; Aquilante, Francesco; Lindh, Roland; González, Leticia

    2014-05-01

    We present a formulation of analytical energy gradients at the complete active space self-consistent field (CASSCF) level of theory employing density fitting (DF) techniques to enable efficient geometry optimizations of large systems. As an example, the ground and lowest triplet state geometries of a ruthenium nitrosyl complex are computed at the DF-CASSCF level of theory and compared with structures obtained from density functional theory (DFT) using the B3LYP, BP86, and M06L functionals. The average deviation of all bond lengths compared to the crystal structure is 0.042 Å at the DF-CASSCF level of theory, which is slightly larger but still comparable with the deviations obtained by the tested DFT functionals, e.g., 0.032 Å with M06L. Specifically, the root-mean-square deviation between the DF-CASSCF and best DFT coordinates, delivered by BP86, is only 0.08 Å for S0 and 0.11 Å for T1, indicating that the geometries are very similar. While keeping the mean energy gradient errors below 0.25%, the DF technique results in a 13-fold speedup compared to the conventional CASSCF geometry optimization algorithm. Additionally, we assess the singlet-triplet energy vertical and adiabatic differences with multiconfigurational second-order perturbation theory (CASPT2) using the DF-CASSCF and DFT optimized geometries. It is found that the vertical CASPT2 energies are relatively similar regardless of the geometry employed whereas the adiabatic singlet-triplet gaps are more sensitive to the chosen triplet geometry. PMID:24811621

  6. Analytical gradients of complete active space self-consistent field energies using Cholesky decomposition: Geometry optimization and spin-state energetics of a ruthenium nitrosyl complex

    NASA Astrophysics Data System (ADS)

    Delcey, Mickaël G.; Freitag, Leon; Pedersen, Thomas Bondo; Aquilante, Francesco; Lindh, Roland; González, Leticia

    2014-05-01

    We present a formulation of analytical energy gradients at the complete active space self-consistent field (CASSCF) level of theory employing density fitting (DF) techniques to enable efficient geometry optimizations of large systems. As an example, the ground and lowest triplet state geometries of a ruthenium nitrosyl complex are computed at the DF-CASSCF level of theory and compared with structures obtained from density functional theory (DFT) using the B3LYP, BP86, and M06L functionals. The average deviation of all bond lengths compared to the crystal structure is 0.042 Å at the DF-CASSCF level of theory, which is slightly larger but still comparable with the deviations obtained by the tested DFT functionals, e.g., 0.032 Å with M06L. Specifically, the root-mean-square deviation between the DF-CASSCF and best DFT coordinates, delivered by BP86, is only 0.08 Å for S0 and 0.11 Å for T1, indicating that the geometries are very similar. While keeping the mean energy gradient errors below 0.25%, the DF technique results in a 13-fold speedup compared to the conventional CASSCF geometry optimization algorithm. Additionally, we assess the singlet-triplet energy vertical and adiabatic differences with multiconfigurational second-order perturbation theory (CASPT2) using the DF-CASSCF and DFT optimized geometries. It is found that the vertical CASPT2 energies are relatively similar regardless of the geometry employed whereas the adiabatic singlet-triplet gaps are more sensitive to the chosen triplet geometry.

  7. Orthogonal tensor decompositions

    SciTech Connect

    Tamara G. Kolda

    2000-03-01

    The authors explore the orthogonal decomposition of tensors (also known as multi-dimensional arrays or n-way arrays) using two different definitions of orthogonality. They present numerous examples to illustrate the difficulties in understanding such decompositions. They conclude with a counterexample to a tensor extension of the Eckart-Young SVD approximation theorem by Leibovici and Sabatier [Linear Algebra Appl. 269(1998):307--329].

  8. Nonlinear ordinary difference equations

    NASA Technical Reports Server (NTRS)

    Caughey, T. K.

    1979-01-01

    Future space vehicles will be relatively large and flexible, and active control will be necessary to maintain geometrical configuration. While the stresses and strains in these space vehicles are not expected to be excessively large, their cumulative effects will cause significant geometrical nonlinearities to appear in the equations of motion, in addition to the nonlinearities caused by material properties. Since the only effective tool for the analysis of such large complex structures is the digital computer, it will be necessary to gain a better understanding of the nonlinear ordinary difference equations which result from the time discretization of the semidiscrete equations of motion for such structures.

  9. Consequences of using nonlinear particle trajectories to compute spatial diffusion coefficients. [for cosmic ray propagation in interstellar and interplanetary space

    NASA Technical Reports Server (NTRS)

    Goldstein, M. L.

    1977-01-01

    In a study of cosmic ray propagation in interstellar and interplanetary space, a perturbed orbit resonant scattering theory for pitch angle diffusion in a slab model of magnetostatic turbulence is slightly generalized and used to compute the diffusion coefficient for spatial propagation parallel to the mean magnetic field. This diffusion coefficient has been useful for describing the solar modulation of the galactic cosmic rays, and for explaining the diffusive phase in solar flares in which the initial anisotropy of the particle distribution decays to isotropy.

  10. Solution of a few nonlinear problems in aerodynamics by the finite elements and functional least squares methods. Ph.D. Thesis - Paris Univ.; [mathematical models of transonic flow using nonlinear equations

    NASA Technical Reports Server (NTRS)

    Periaux, J.

    1979-01-01

    The numerical simulation of the transonic flows of idealized fluids and of incompressible viscous fluids, by the nonlinear least squares methods is presented. The nonlinear equations, the boundary conditions, and the various constraints controlling the two types of flow are described. The standard iterative methods for solving a quasi elliptical nonlinear equation with partial derivatives are reviewed with emphasis placed on two examples: the fixed point method applied to the Gelder functional in the case of compressible subsonic flows and the Newton method used in the technique of decomposition of the lifting potential. The new abstract least squares method is discussed. It consists of substituting the nonlinear equation by a problem of minimization in a H to the minus 1 type Sobolev functional space.