Search NASA⌕ Search

SEARCH · Search NASA

Results for “Analytic linearization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

A site-selective antiferromagnetic ground state in layered pnictide-oxide BaTi{sub 2}As{sub 2}O

The electronic and magnetic properties of BaTi{sub 2}As{sub 2}O have been investigated using both the first-principles and analytical methods. The full-potential linearized augmented plane-wave calculations show that the most stable state is a site-selective antiferromagnetic (AFM) metal with a 2×1×1 magnetic unit cell containing two nonmagnetic Ti atoms and two other Ti atoms with antiparallel moments. Further analysis to Fermi surface and spin susceptibility shows that the site-selective AFM ground state is driven by the Fermi surface nesting and the Coulomb correlation. Meanwhile, the charge density distribution remains uniform, suggesting that the phase transition at 200 K in experiment is a spin-density-wave transition.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Modulating resonance behaviors by noise recycling in bistable systems with time delay

In this paper, the impact of noise recycling on resonance behaviors is studied theoretically and numerically in a prototypical bistable system with delayed feedback. According to the interior cooperating and interacting activity of noise recycling, a theory has been proposed by reducing the non-Markovian problem into a two-state model, wherein both the master equation and the transition rates depend on not only the current state but also the earlier two states due to the recycling lag and the feedback delay. By virtue of this theory, the formulae of the power spectrum density and the linear response function have been found analytically. And the theoretical results are well verified by numerical simulations. It has been demonstrated that both the recycling lag and the feedback delay play a crucial role in the resonance behaviors. In addition, the results also suggest an alternative scheme to modulate or control the coherence or stochastic resonance in bistable systems with time delay.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Numerical simulation of scattering of acoustic waves by inelastic bodies using hypersingular boundary integral equation

A numerical method for solving a problem of diffraction of acoustic waves by system of solid and thin objects based on the reduction the problem to a boundary integral equation in which the integral is understood in the sense of finite Hadamard value is proposed. To solve this equation we applied piecewise constant approximations and collocation methods numerical scheme. The difference between the constructed scheme and earlier known is in obtaining approximate analytical expressions to appearing system of linear equations coefficients by separating the main part of the kernel integral operator. The proposed numerical scheme is tested on the solution of the model problem of diffraction of an acoustic wave by inelastic sphere.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Three species one-dimensional kinetic model for weakly ionized plasmas

A three species one-dimensional kinetic model is presented for a spatially homogeneous weakly ionized plasma subjected to the action of a time varying electric field. Planar geometry is assumed, which means that the plasma evolves in the privileged direction of the field. The energy transmitted to the electric charges is channelized to the neutrals thanks to collisions, a mechanism that influences the plasma dynamics. Charge-charge interactions have been designed as a one-dimensional collision term equivalent to the Landau operator used for fully ionized plasmas. Charge-neutral collisions are modelled by a conservative drift-diffusion operator in the Dougherty's form. The resulting set of coupled integro-differential equations is solved with the stable and robust propagator integral method. This semi–analytical method feasibility accounts for non–linear effects without appealing to linearisation or simplifications, providing conservative physically meaningful solutions even for initial or emerging sharp velocity distribution function profiles. It is found that charge-neutral collisions exert a significant effect since a quite different plasma evolution arises if compared to the collisionless limit. In addition, substantial differences in the system motion are found for constant and temperature dependent collision frequencies cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Algebraic motion of vertically displacing plasmas

In this paper, the vertical motion of a tokamak plasma is analytically modelled during its non-linear phase by a free-moving current-carrying rod inductively coupled to a set of fixed conducting wires or a cylindrical conducting shell. The solutions capture the leading term in a Taylor expansion of the Green's function for the interaction between the plasma column and the surrounding vacuum vessel. The plasma shape and profiles are assumed not to vary during the vertical drifting phase such that the plasma column behaves as a rigid body. In the limit of perfectly conducting structures, the plasma is prevented to come in contact with the wall due to steep effective potential barriers created by the induced Eddy currents. Resistivity in the wall allows the equilibrium point to drift towards the vessel on the slow timescale of flux penetration. The initial exponential motion of the plasma, understood as a resistive vertical instability, is succeeded by a non-linear “sinking” behaviour shown to be algebraic and decelerating. Finally, the acceleration of the plasma column often observed in experiments is thus concluded to originate from an early sharing of toroidal current between the core, the halo plasma, and the wall or from the thermal quench dynamics precipitating loss of plasma current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Coupling to rotational manifolds to improve gas-phase pump–probe spectroscopic models

The physical picture of gas-phase optical transitions is normally presented as an isolated two-level system balanced by upward and downward processes. Isolated models assume a phenomenological treatment of collisional dephasing but do not strictly account for collisional population exchange with the rotational baths. While this assumption is valid under low-intensity conditions, where excitation is rate-limiting, isolated models can deviate from Beer’s Law at sufficient pressures and monochromatic intensities when both collisional broadening and power broadening are comparable to (or greater than) lifetime broadening, which are not uncommon conditions for cavity enhanced spectroscopies in the mid-IR spectral range. Although this problem has been addressed by rate-equation models for linear absorption measurements, a general treatment for multi-level quantum mechanical models suitable for non-linear absorption measurements (two-photon/two-color/pump–probe) is lacking. Isolated models require physical parameter inputs that disagree with expected values by at least an order of magnitude. These non-physical models undermine the ability to predict non-linear signal strengths under untested conditions and thereby limit the potential to optimize the sensitivity of non-linear spectroscopies and to expand their analytical applications (e.g., new analytes and/or buffer gases, changes in cavity free-spectral-range, changes in intracavity powers or wavelengths, and accurate investigation of physical phenomena). In this study, we derive bath-coupled models for gaseous pump–probe spectroscopy by application of the quantum Lindblad equation and detailed balance. Bath-coupled models are shown to fit data consistently across variations in intensity and agree with all physically expected values.

Cavity ring-down spectroscopy↗
DOE OSTI2019

Formation of spectra of saturated-absorption resonances on closed transitions in the spectroscopy of unidirectional waves

The formation of spectra of saturated-absorption resonances on the atomic transitions with level momenta J = 0 → J = 1 and J = 1 → J = 2 are investigated, both numerically and analytically, by the method of unidirectional linearly polarised laser waves in order to determine how the transition openness (radiation branching) affects these spectra. It is shown that, along with quantitative changes, determined by the relaxation constants of levels, value of their splitting, polarisation orientation, and wave intensities, the spectra undergo qualitative changes with a change in the degree of atomic transition openness. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

THE THERMAL INSTABILITY OF SOLAR PROMINENCE THREADS

The fine structure of solar prominences and filaments appears as thin and long threads in high-resolution images. In H{alpha} observations of filaments, some threads can be observed for only 5-20 minutes before they seem to fade and eventually disappear, suggesting that these threads may have very short lifetimes. The presence of an instability might be the cause of this quick disappearance. Here, we study the thermal instability of prominence threads as an explanation of their sudden disappearance from H{alpha} observations. We model a prominence thread as a magnetic tube with prominence conditions embedded in a coronal environment. We assume a variation of the physical properties in the transverse direction so that the temperature and density continuously change from internal to external values in an inhomogeneous transitional layer representing the particular prominence-corona transition region (PCTR) of the thread. We use the nonadiabatic and resistive magnetohydrodynamic equations, which include terms due to thermal conduction parallel and perpendicular to the magnetic field, radiative losses, heating, and magnetic diffusion. We combine both analytical and numerical methods to study linear perturbations from the equilibrium state, focusing on unstable thermal solutions. We find that thermal modes are unstable in the PCTR for temperatures higher than 80,000 K, approximately. These modes are related to temperature disturbances that can lead to changes in the equilibrium due to rapid plasma heating or cooling. For typical prominence parameters, the instability timescale is of the order of a few minutes and is independent of the form of the temperature profile within the PCTR of the thread. This result indicates that thermal instability may play an important role for the short lifetimes of threads in the observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

PERTURBATION THEORY OF THE COSMOLOGICAL LOG-DENSITY FIELD

The matter density field exhibits a nearly lognormal probability density distribution after entering into the nonlinear regime. Recently, it has been shown that the shape of the power spectrum of a logarithmically transformed density field is very close to the linear density power spectrum, motivating an analytic study of it. In this paper, we develop cosmological perturbation theory for the power spectrum of this field. Our formalism is developed in the context of renormalized perturbation theory, which helps to regulate the convergence behavior of the perturbation series, and of the Taylor series expansion we use for the logarithmic mapping. This approach allows us to handle the critical issue of density smoothing in a straightforward way. We also compare our perturbative results with simulation measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

UNCERTAINTIES IN ATOMIC DATA AND THEIR PROPAGATION THROUGH SPECTRAL MODELS. I

We present a method for computing uncertainties in spectral models, i.e., level populations, line emissivities, and emission line ratios, based upon the propagation of uncertainties originating from atomic data. We provide analytic expressions, in the form of linear sets of algebraic equations, for the coupled uncertainties among all levels. These equations can be solved efficiently for any set of physical conditions and uncertainties in the atomic data. We illustrate our method applied to spectral models of O III and Fe II and discuss the impact of the uncertainties on atomic systems under different physical conditions. As to intrinsic uncertainties in theoretical atomic data, we propose that these uncertainties can be estimated from the dispersion in the results from various independent calculations. This technique provides excellent results for the uncertainties in A-values of forbidden transitions in [Fe II].

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2010

N-body simulations with generic non-Gaussian initial conditions I: power spectrum and halo mass function

We address the issue of setting up generic non-Gaussian initial conditions for N-body simulations. We consider inflationary-motivated primordial non-Gaussianity where the perturbations in the Bardeen potential are given by a dominant Gaussian part plus a non-Gaussian part specified by its bispectrum. The approach we explore here is suitable for any bispectrum, i.e. it does not have to be of the so-called separable or factorizable form. The procedure of generating a non-Gaussian field with a given bispectrum (and a given power spectrum for the Gaussian component) is not univocal, and care must be taken so that higher-order corrections do not leave a too large signature on the power spectrum. This is so far a limiting factor of our approach. We then run N-body simulations for the most popular inflationary-motivated non-Gaussian shapes. The halo mass function and the non-linear power spectrum agree with theoretical analytical approximations proposed in the literature, even if they were so far developed and tested only for a particular shape (the local one). We plan to make the simulations outputs available to the community via the non-Gaussian simulations comparison project web site http://icc.ub.edu/∼liciaverde/NGSCP.html.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Detection of the large-scale tidal field with galaxy multiplet alignment in the DESI Y1 spectroscopic survey

We explore correlations between the orientations of small galaxy groups, or ‘multiplets’, and the large-scale gravitational tidal field. Using data from the Dark Energy Spectroscopic Instrument (DESI) Y1 survey, we detect the intrinsic alignment (IA) of multiplets to the galaxy-traced matter field out to separations of $100\,h^{-1}$ Mpc. Unlike traditional IA measurements of individual galaxies, this estimator is not limited by imaging of galaxy shapes and allows for direct IA detection beyond redshift $z=1$. Multiplet alignment is a form of higher order clustering, for which the scale-dependence traces the underlying tidal field and amplitude is a result of small-scale ($\lt 1h^{-1}$ Mpc) dynamics. Within samples of bright galaxies, luminous red galaxies (LRG) and emission-line galaxies, we find similar scale-dependence regardless of intrinsic luminosity or colour. This is promising for measuring tidal alignment in galaxy samples that typically display no IA. DESI’s LRG mock galaxy catalogues created from the A BACUS S UMMIT N -body simulations produce a similar alignment signal, though with a 33 per cent lower amplitude at all scales. An analytic model using a non-linear power spectrum (NLA) only matches the signal down to 20 $h^{-1}$ Mpc. Our detection demonstrates that galaxy clustering in the non-linear regime of structure formation preserves an interpretable memory of the large-scale tidal field. Multiplet alignment complements traditional two-point measurements by retaining directional information imprinted by tidal forces, and contains additional line-of-sight information compared to weak lensing. This is a more effective estimator than the alignment of individual galaxies in dense, blue, or faint galaxy samples.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

On the dual advantage of placing observations through forward sensitivity analysis

The four-dimensional variational data assimilation methodology for assimilating noisy observations into a deterministic model has been the workhorse in forecasting centres for over three decades. While this method provides a computationally efficient framework for dynamic data assimilation, it is largely silent on the important question concerning the minimum number and placement of observations. To answer this question, we demonstrate the dual advantage of placing the observations where the square of the sensitivity of the model solution with respect to the unknown control variables, called forward sensitivities, attains its maximum. By following this approach, we can force the observability Gramian to be of full rank, which in turn guarantees efficient recovery of the optimal values of the control variables, which is the first of the two advantages of this strategy. We further show that the proposed strategy of placing observations has another inherent optimality: the square of the sensitivity of the optimal estimates of the control with respect to the observations (used to obtain these estimates) attains its minimum value, a second advantage that is a direct consequence of the above strategy for placing observations. Our analytical framework and numerical experiments on linear and nonlinear systems confirm the effectiveness of our proposed strategy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

Gravitational lens equation for embedded lenses; magnification and ellipticity

We give the lens equation for light deflections caused by point mass condensations in an otherwise spatially homogeneous and flat universe. We assume the signal from a distant source is deflected by a single condensation before it reaches the observer. We call this deflector an embedded lens because the deflecting mass is part of the mean density. The embedded lens equation differs from the conventional lens equation because the deflector mass is not simply an addition to the cosmic mean. We prescribe an iteration scheme to solve this new lens equation and use it to compare our results with standard linear lensing theory. We also compute analytic expressions for the lowest order corrections to image amplifications and distortions caused by incorporating the lensing mass into the mean. We use these results to estimate the effect of embedding on strong lensing magnifications and ellipticities and find only small effects, <1%, contrary to what we have found for time delays and for weak lensing, {approx}5%.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2011

Temporal enhancement of super-horizon curvature perturbations from decays of two curvatons and its cosmological consequences

If more than one curvaton dominate the Universe at different epochs from each other, curvature perturbations can be temporarily enhanced to a value much larger than the observed one 10{sup -5}. The traces of the enhancement may be left as higher order correlation functions, that is, as non-Gaussianity, the stochastic gravitational waves that are sourced by scalar-scalar mode couplings, as well as the primordial black holes that are formed by the gravitational collapse of the enhanced curvature perturbations. We first confirm that such a temporal enhancement indeed occurs by solving the linearized perturbation equations both numerically and analytically. We then derive an analytic expression of the full-order curvature perturbation which does not rely on the frequently used sudden decay approximation and is exact on superhorizon scales. By using this analytic formula, we provide expressions of the nonlinearity parameters f{sub NL}, {tau}{sub NL} and g{sub NL}. If both two curvatons contribute to the final curvature perturbations, the strong non-Gaussianity appears in the trispectrum rather than in the bispectrum. We also find a unique consistency relation between {tau}{sub NL} and g{sub NL} without f{sub NL}. By using the second-order perturbation theory, we numerically show that the spectrum of the induced gravitational waves has a plateau corresponding to duration of the enhancement and such gravitational waves can be probed by ultimate-DECIGO and space-based atomic interferometers. We finally calculate the abundance of the primordial black holes and put a constraint on the amplitude of the enhanced curvature perturbations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Entanglement Entropy in Lifshitz Theories

We discuss and compute entanglement entropy (EE) in (1+1)-dimensional free Lifshitz scalar field theories with arbitrary dynamical exponents. We consider both the subinterval and periodic sublattices in the discretized theory as subsystems. In both cases, we are able to analytically demonstrate that the EE grows linearly as a function of the dynamical exponent. Furthermore, for the subinterval case, we determine that as the dynamical exponent increases, there is a crossover from an area law to a volume law. Lastly, we deform Lifshitz field theories with certain relevant operators and show that the EE decreases from the ultraviolet to the infrared fixed point, giving evidence for a possible c c -theorem for deformed Lifshitz theories.

He, Temple↗
NASA NTRS1972

Determination of decay coefficients for combustors with acoustic absorbers

An analytical technique for the calculation of linear decay coefficients in combustors with acoustic absorbers is presented. Tuned circumferential slot acoustic absorbers were designed for the first three transverse modes of oscillation, and decay coefficients for these absorbers were found as a function of backing distance for seven different chamber configurations. The effectiveness of the absorbers for off-design values of the combustion response and acoustic mode is also investigated. Results indicate that for tuned absorbers the decay coefficient increases approximately as the cube of the backing distance. For most off-design situations the absorber still provides a damping effect. However, if an absorber designed for some higher mode of oscillation is used to damp lower mode oscillations, a driving effect is frequently found.

Mitchell, C. E.↗
NASA NTRS1973

Mixed finite-difference scheme for analysis of simply supported thick plates.

A mixed finite-difference scheme is presented for the stress and free vibration analysis of simply supported nonhomogeneous and layered orthotropic thick plates. The analytical formulation is based on the linear, three-dimensional theory of orthotropic elasticity and a Fourier approach is used to reduce the governing equations to six first-order ordinary differential equations in the thickness coordinate. The governing equations possess a symmetric coefficient matrix and are free of derivatives of the elastic characteristics of the plate. In the finite difference discretization two interlacing grids are used for the different fundamental unknowns in such a way as to reduce both the local discretization error and the bandwidth of the resulting finite-difference field equations. Numerical studies are presented for the effects of reducing the interior and boundary discretization errors and of mesh refinement on the accuracy and convergence of solutions. It is shown that the proposed scheme, in addition to a number of other advantages, leads to highly accurate results, even when a small number of finite difference intervals is used.

Noor, A. K.↗