Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical stabilization”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Thermodynamic stability of a spin microemulsion in Rashba spin-orbit-coupled bosons

Recent finite-temperature numerical simulations have unveiled a quantum “spin” microemulsion analog, found by raising the temperature of a stripe supersolid phase in a Rashba spin-orbit coupled Bose gas. This microemulsion state is a highly correlated, isotropic normal fluid where atoms self-arrange based on their internal pseudospin into patterns that resemble bicontinuous microemulsions. This finding leaves several open questions regarding the broader accessibility of this phase in experiments. Here, we use equilibrium finite-temperature numerical simulations based on a coherent-state path integral representation to perform a computational investigation into the thermodynamic stability of the spin microemulsion state. Numerical simulations emphasize the requirement of a nearly, but not perfectly, isotropic spin-orbit coupling in order to achieve the microemulsion phase in cold-atom experiments. Moreover, the microemulsion state exists independent of miscibility of the pseudospin components and for a wide range of pseudospin population imbalance, suggesting a high degree of flexibility in choosing the atom and hyperfine states in an experimental realization. Lastly, we demonstrate this feasibility by mimicking a Rashba spin-orbit-coupled 87 Rb experiment in an isotropic harmonic trap, where we confirm the microemulsion's existence via its density profile and equilibrium quasimomentum distribution.

Complex Langevin dynamics↗
DOE OSTI2022

Construction of wave dark matter halos: Numerical algorithm and analytical constraints

Here we present a wave generalization of the classic Schwarzschild method for constructing self-consistent halos—such a halo consists of a suitable superposition of waves instead of particle orbits, chosen to yield a desired mean density profile. As an illustration, the method is applied to spherically symmetric halos. We derive an analytic relation between the particle distribution function and the wave superposition amplitudes and show how it simplifies in the high-energy (WKB) limit. We verify the stability of such constructed halos by numerically evolving the Schrödinger-Poisson system. The algorithm provides an efficient and accurate way to simulate the time-dependent halo substructures from wave interference. We use this method to construct halos with a variety of density profiles, all of which have a core from the ground-state wave function, though the core-halo relation need not be the standard one.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Nonlinear localized modes in PT-symmetric optical media with competing gain and loss

The existence and stability of the nonlinear spatial localized modes are investigated in parity-time symmetric optical media characterized by a generic complex hyperbolic refractive index distribution with competing gain and loss profile. The exact analytical expression of the localized modes are found for all values of the competing parameter and in the presence of both the self-focusing and self-defocusing Kerr nonlinearity. The effects of competing gain/loss profile on the stability structure of these localized modes are discussed with the help of linear stability analysis followed by the direct numerical simulation of the governing equation. The spatial localized modes in two-dimensional geometry as well as the transverse power-flow density associated with these localized modes are also examined. -- Highlights: • Existence of localized modes is investigated in PT-symmetric complex potentials. • Exact analytical expression of the localized modes is obtained. • Effect of gain/loss profile on the stability of these localized modes is discussed. • Localized modes in 2D and associated transverse power-flow density are also examined.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1981

Cylindrical annulus models

The stability diagram for the atmospheric general circulation experiment (AGCE) was numerically determined. The cylindrical code of Warn-Varnas modified for linear stability and adapted to spherical geometry. The code computes axisymmetric basic states and then test their stability to zonal disturbances. The Warn-Varnas code and some of the basic state results are reported.

Leslie, F.↗
DOE OSTI2021

On Distributed Model-Free Reinforcement Learning Control with Stability Guarantee

Distributed learning can enable scalable and effective decision making in numerous complex cyber-physical systems such as smart transportation, robotics swarm, power systems, etc. However, the stability of the system is usually not guaranteed in most existing learning paradigms; and this limitation can hinder the wide deployment of machine learning in decision making of safety-critical systems. This paper presents a stability guaranteed distributed reinforcement learning (SGDRL) framework for interconnected linear subsystems, without knowing the subsystem models. While the learning process requires data from a peer-to-peer (p2p) communication architecture, the control implementation of each subsystem is only based on its local state. The stability of the interconnected subsystems will be ensured by a diagonally dominant eigenvalue condition, which will then be used in a model-free RL algorithm to learn the feedback control gains. The RL algorithm structure follows an off-policy iterative framework, with interleaved policy evaluation and policy update steps. We numerically validate our theoretical results by performing simulations on four interconnected sub-systems.

Mukherjee, Sayak↗
DOE OSTI2016

Nonlinear Elastic Plate in a Flow of Gas: Recent Results and Conjectures

We give a survey of recent results on flow-structure interactions modeled by a modified wave equation coupled at an interface with equations of nonlinear elasticity. Both subsonic and supersonic flow velocities are considered. The focus of the discussion here is on the interesting mathematical aspects of physical phenomena occurring in aeroelasticity, such as flutter and divergence. This leads to a partial differential equation treatment of issues such as well-posedness of finite energy solutions, and long-time (asymptotic) behavior. The latter includes theory of asymptotic stability, convergence to equilibria, and to global attracting sets. We complete the discussion with several well known observations and conjectures based on experimental/numerical studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1990

The inviscid stability of supersonic flow past a sharp cone

The laminar boundary layer which forms on a sharp cone in a supersonic freestream, where lateral curvature plays a key role in the physics of the problem is considered. This flow is then analyzed from the point of view of linear, temporal, inviscid stability. The basic, non-axisymmetric disturbance equations are derived for general flows of this class, and a so called triply generalized inflexion condition is found for the existence of subsonic neutral modes of instability. This condition is analogous to the well-known generalized inflexion condition found in planar flows, although in the present case the condition depends on both axial and aximuthal wavenumbers. Extensive numerical results are presented for the stability problem at a freestream Mach number of 3.8, for a range of streamwise locations. These results reveal that a new mode of instability may occur, peculiar to flows of this type involving curvature. Additionally, asymptotic analyses valid close to the tip of the cone, far downstream of the cone are presented, and these give a partial (asymptotic) description of this additional mode of instability.

Duck, Peter W.↗
NASA NTRS2003

Technical Evaluation Report, Part A - Vortex Flow and High Angle of Attack

A symposium entitled Vortex Flow and High Angle of Attack was held in Loen, Norway, from May 7 through May 11, 2001. The Applied Vehicle Technology (AVT) panel, under the auspices of the Research and Technology Organization (RTO), sponsored this symposium. Forty-eight papers, organized into nine sessions, addressed computational and experimental studies of vortex flows pertinent to both aircraft and maritime applications. The studies also ranged from fundamental fluids investigations to flight test results, and significant results were contributed from a broad range of countries. The principal emphasis of this symposium was on "the understanding and prediction of separation-induced vortex flows and their effects on military vehicle performance, stability, control, and structural design loads." It was further observed by the program committee that "separation- induced vortex flows are an important part of the design and off-design performance of conventional fighter aircraft and new conventional or unconventional manned or unmanned advanced vehicle designs (UAVs, manned aircraft, missiles, space planes, ground-based vehicles, and ships)." The nine sessions addressed the following topics: vortical flows on wings and bodies, experimental techniques for vortical flows, numerical simulations of vortical flows, vortex stability and breakdown, vortex flows in maritime applications, vortex interactions and control, vortex dynamics, flight testing, and vehicle design. The purpose of this paper is to provide brief reviews of these papers along with some synthesizing perspectives toward future vortex flow research opportunities. The paper includes the symposium program. (15 refs.)

VORTEX BREAKDOWN↗
NASA NTRS1993

Mapping the stability region of the 3:2 Neptune-Pluto resonance

Pluto and Charon are most likely the remnants of a large number of objects that existed in the Uranus-Neptune region at early epochs of the solar system. Numerical integrations have shown that, in general, such objects were ejected from the planetary region on timescales of approximately 10(exp 7) years after Neptune and Uranus reached their current masses. It is thought that the Pluto-Charon system survived to current times without being dynamically removed in this way because it is trapped in a set of secular and mean motion resonances with Neptune. The best-known Pluto-Neptune orbit coupling is the 3:2 mean motion resonance discovered almost 30 years ago by C. Cohen and E. Hubbard. These workers showed that the resonance angle, delta is equivalent to 3(lambda(sub P)) - 2(lambda(sub N)) - omega-bar(sub P) where omega-bar(sub P) is the longitude of perihelion of the Pluto-Charon system, and lambda(sub N) and lambda(sub P) are the mean longitude of Neptune and Pluto-Charon respectively, librates about 180 deg with an amplitude, A(sub delta), of 76 deg. A numerical simulation project to map out the stability region of the 3:2 resonance is reported. The results of these simulations are important to understanding whether Pluto's long-term heliocentric stability requires only the 3:2 resonance, or whether it instead requires one or more of the other Pluto-Neptune resonances. Our study also has another important application. By investigating stability timescales as a function of orbital elements, we gain insight into the fraction of orbital phase space which the stable 3:2 resonance occupies. This fraction is directly related to the probability that the Pluto-Charon system (and possibly other small bodies) could have been captured into this resonance.

Levison, H. F.↗
NASA NTRS2010

Receptivity of Hypersonic Boundary Layers over Straight and Flared Cones

The effects of adverse pressure gradients on the receptivity and stability of hypersonic boundary layers were numerically investigated. Simulations were performed for boundary layer flows over a straight cone and two flared cones. The steady and the unsteady flow fields were obtained by solving the two-dimensional Navier-Stokes equations in axi-symmetric coordinates using the 5th order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The mean boundary layer profiles were analyzed using local stability and non-local parabolized stability equations (PSE) methods. After the most amplified disturbances were identified, two-dimensional plane acoustic waves were introduced at the outer boundary of the computational domain and time accurate simulations were performed. The adverse pressure gradient was found to affect the boundary layer stability in two important ways. Firstly, the frequency of the most amplified second-mode disturbance was increased relative to the zero pressure gradient case. Secondly, the amplification of first- and second-mode disturbances was increased. Although an adverse pressure gradient enhances instability wave growth rates, small nose-tip bluntness was found to delay transition due to the low receptivity coefficient and the resulting weak initial amplitude of the instability waves. The computed and measured amplitude-frequency spectrums in all three cases agree very well in terms of frequency and the shape except for the amplitude.

Balakumar, Ponnampalam↗
DOE OSTI2013

Numerical integration of the extended variable generalized Langevin equation with a positive Prony representable memory kernel

Generalized Langevin dynamics (GLD) arise in the modeling of a number of systems, ranging from structured fluids that exhibit a viscoelastic mechanical response, to biological systems, and other media that exhibit anomalous diffusive phenomena. Molecular dynamics (MD) simulations that include GLD in conjunction with external and/or pairwise forces require the development of numerical integrators that are efficient, stable, and have known convergence properties. In this article, we derive a family of extended variable integrators for the Generalized Langevin equation with a positive Prony series memory kernel. Using stability and error analysis, we identify a superlative choice of parameters and implement the corresponding numerical algorithm in the LAMMPS MD software package. Salient features of the algorithm include exact conservation of the first and second moments of the equilibrium velocity distribution in some important cases, stable behavior in the limit of conventional Langevin dynamics, and the use of a convolution-free formalism that obviates the need for explicit storage of the time history of particle velocities. Capability is demonstrated with respect to accuracy in numerous canonical examples, stability in certain limits, and an exemplary application in which the effect of a harmonic confining potential is mapped onto a memory kernel.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1982

On the stability of Schwarzschild's triaxial galaxy model

The numerical model for a triaxial galaxy constructed by Schwartzschild (1979) has been tested by means of fully self-consistent three-dimensional n-body numerical experiments that use a representation of this model to start the integration. The integration, study models, initial load, and correcting potential of the experiment are detailed. Preliminary checks included density values, potential values, and integration in potential of an ideal model. The results showed some particle loss and no significant shape change. The model is rugged and robust, and no growing disturbances with growth rate in excess of 0.5 per crossing time could be detected. Runs were done at 50 percent excess halo and with tube orbits reversed.

Smith, B. F.↗
NASA NTRS1999

Experimental and Numerical Optimization of a High-Lift System to Improve Low-Speed Performance, Stability, and Control of an Arrow-Wing Supersonic Transport

An investigation was performed to evaluate leading-and trailing-edge flap deflections for optimal aerodynamic performance of a High-Speed Civil Transport concept during takeoff and approach-to-landing conditions. The configuration used for this study was designed by the Douglas Aircraft Company during the 1970's. A 0.1-scale model of this configuration was tested in the Langley 30- by 60-Foot Tunnel with both the original leading-edge flap system and a new leading-edge flap system, which was designed with modem computational flow analysis and optimization tools. Leading-and trailing-edge flap deflections were generated for the original and modified leading-edge flap systems with the computational flow analysis and optimization tools. Although wind tunnel data indicated improvements in aerodynamic performance for the analytically derived flap deflections for both leading-edge flap systems, perturbations of the analytically derived leading-edge flap deflections yielded significant additional improvements in aerodynamic performance. In addition to the aerodynamic performance optimization testing, stability and control data were also obtained. An evaluation of the crosswind landing capability of the aircraft configuration revealed that insufficient lateral control existed as a result of high levels of lateral stability. Deflection of the leading-and trailing-edge flaps improved the crosswind landing capability of the vehicle considerably; however, additional improvements are required.

Hahne, David E.↗
NASA NTRS2009

Mean Element Propagations Using Numerical Averaging

The long-term evolution characteristics (and stability) of an orbit are best characterized using a mean element propagation of the perturbed two body variational equations of motion. The averaging process eliminates short period terms leaving only secular and long period effects. In this study, a non-traditional approach is taken that averages the variational equations using adaptive numerical techniques and then numerically integrating the resulting EOMs. Doing this avoids the Fourier series expansions and truncations required by the traditional analytic methods. The resultant numerical techniques can be easily adapted to propagations at most solar system bodies.

orbit propagations↗
DOE OSTI2018

Analysis of fractional-order models for hepatitis B

This paper presents two models for hepatitis B, both given by fractional differential equations. The first model is formulated without parameters that indicate drug therapy, while the second one considers the drug therapy. The basic reproduction number and the stability analysis are considered for both models. Moreover, some numerical simulations by nonstandard finite difference schemes are presented. The numerical results show that the solutions converges to an equilibrium point as predicted in the stability analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Reliable Devices Yield Stable Quantum Computations

Stable quantum computation requires noisy results to remain bounded even in the presence of noise fluctuations. Yet non-stationary noise processes lead to drift in the varying characteristics of a quantum device that can greatly influence the circuit outcomes. Here we address how temporal and spatial variations in noise relate device reliability to quantum computing stability. First, our approach quantifies the differences in statistical distributions of characterization metrics collected at different times and locations using Hellinger distance. We then validate an analytical bound that relates this distance directly to the stability of a computed expectation value. Our demonstration uses numerical simulations with models informed by the washington superconducting transmon device. We find that the stability metric is consistently bounded from above by the corresponding Hellinger distance, which can be cast as a specified tolerance level. These results underscore the significance of reliable quantum computing devices and the impact for stable quantum computation.

Dasgupta, Samudra↗
DOE OSTI2013

Ion temperature gradient instability at sub-Larmor radius scales with non-zero ballooning angle

Linear gyro-kinetic stability calculations predict unstable toroidal ion temperature gradient modes (ITGs) with normalised poloidal wave vectors well above one (k{sub {theta}}{rho}{sub i}>1) for standard tokamak parameters with adiabatic electron response. These modes have a maximum amplitude at a poloidal angle {theta} that is shifted away from the low field side ({theta}{ne}0). The physical mechanism is clarified through the use of a fluid model. It is shown that the shift of the mode away from the low field side reduces the effective drift frequency which allows for the instability to develop. Numerical tests using the gyro-kinetic model confirm this physical mechanism. Furthermore, it is shown that modes localized away from the low field side can be important also for k{sub {theta}}{rho}{sub i}<1 close to the threshold of the ITG. In fact, modes with maximum amplitude at {theta}{ne}0 can exist for normalised temperature gradient lengths below the threshold of the ITG obtained for the case with the maximum at {theta}=0.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗