SEARCH · Search NASA

Results for “numerical dispersion”

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.

SourceResearch connectionFrom yearThrough year

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

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS1987

On the thermal instability of galactic and cluster halos

The paper presents a detailed study of thermal instabilities in cooling flows associated with galaxies and clusters of galaxies. In the case of purely radiation-driven accretion onto a central object such as the cD galaxy M87, it is found that the gas is largely subject to overstability, rather than to monotonic instability. If thermal conductivity is taken into account, the flow is stabilized on scales of several kiloparsecs, even if the conductivity is appreciably reduced (e.g., about 1 percent) with respect to the Spitzer value. In no case are the globular perturbations (i.e., perturbations with comparable radial and azimuthal dimensions) found to be monotonically unstable. The paper presents numerical solutions of the local dispersion relation for the cooling flow in M87 and discusses the possible consequences of the results for a correct understanding of cooling flows.

Malagoli, A.↗

NASA NTRS1988

Comet-solar wind interaction through ion-proton beam instability

The acceleration of cometary ions in the distant cometary tail by an electromagnetic beam instability, which is caused by the relative drift velocity between the cometary ions and solar protons is studied. The linear phase of the instability is analyzed by numerical evaluation of the Vlasov dispersion relation. A particle simulation is used to study the nonlinear phase of the instability. Right-hand polarized cyclotron waves are made unstable by the presence of ion beams. The cometary tail ions are accelerated mainly during the reduced-growth phase of the instability. The cyclotron waves cause momentum transfer from solar wind protons to cometary tail ions. The acceleration of cometary ions through the electromagnetic beam instability is found to be 500-5000 cm/s-squared.

Price, Channon P.↗

NASA NTRS1988

Nature and the nonlinear evolution of electrostatic waves associated with the AMPTE solar wind releases

The nonlinear evolution of the electrostatic waves observed during the AMPTE (Active Magnetosphere Particle Tracer Explorers) solar wind releases is investigated. Previous linear studies indicate two distinct sets of instabilities may be responsible for the generation. One set consists of ion-acoustic type instabilities which are insensitive to the presence of a background magnetic field, while the other group corresponds to the modified two stream instabilities and requires the solar wind flow to be across the ambient magnetic field. To establish which set of instabilities is more viable a detailed linear Vlasov theory has been conducted by numerically solving the full electromagnetic dispersion relation. In addition both the plasma wave as well as the magnetic field measurements by the IRM (Ion Release Module) spacecraft were used to correlate the frequency and the power of the observed waves with the magnitude and the direction of the solar wind magnetic field. The results of these analyses indicate that the ion-acoustic type instabilities have growth rates that are an order of magnitude or more larger than those of the modified two stream instabilities. Results show that both the solar wind protons and the released ions may be heated and accelerated in the directions oblique to the solar wind flow velocity.

Omidi, N.↗

NASA NTRS1989

The oblique behavior of low-frequency electromagnetic waves excited by newborn cometary ions

The free energy in oxygen or hydrogen ions freshly created in the solar wind stimulates low-frequency electromagnetic waves whose growth does not always maximize at parallel propagation. Exploration of the wave vector plane discloses the frequent occurrence of islets of oblique growth unconnected to the unstable parallel modes. Contour plots of the growth rate, real frequency, polarization, and magnetic compression characterize the oblique wave behavior for large values of the initial pitch angle of the cometary particles. Although wave-particle (Landau and cyclotron) resonances feed most of the surveyed oblique instabilities, some are seemingly fluidlike. The results, obtained from the numerical solution of the kinetic dispersion and wave equations, imply that newborn ions can easily excite significant oblique hydromagnetic wave activity. Cometary environments provide the adopted plasma model, but the study is helpful in the interpretation of other low-frequency wave observations in space.

Brinca, Armando L.↗

NASA NTRS1991

Flux-corrected transport techniques for multidimensional compressible magnetohydrodynamics

A prescription is given for conservatively integrating generalized hydromagnetic equations using flux-corrected transport (FCT) techniques. By placing the magnetic-field components at the interface locations of the finite-difference grid, the field is kept divergence-free to within machine roundoff error. The use of FCT techniques allows an integration scheme of high accuracy to be employed, while the numerical ripples associated with large dispersion errors are avoided. The method is particularly well suited for problems involving magnetohydrodynamic shocks and other discontinuities.

Devore, C. Richard↗

NASA NTRS1991

Feedback instability of the ionospheric resonant cavity

A model is developed that provides a theoretical basis for previous numerical results showing a feedback instability with frequencies characteristic of Alfven travel times within the region of the large increase of Alfven speed above the ionosphere. These results have been extended to arbitrary ionospheric conductivity by developing a numerical solution of the cavity dispersion relation that involves Bessel functions of complex order and argument. It is concluded that the large contrast between the magnetospheric and ionospheric Alfven speed leads to the formation of resonant cavity modes with frequencies ranging from 0.1 to 1 Hz. The presence of the cavity leads to a modification of the reflection characteristics of Alfven waves with frequencies that compare to the cavity's normal modes.

Lysak, Robert L.↗

NASA NTRS1991

ISEE 1 observations of electrostatic ion cyclotron waves in association with ion beams on auroral field lines from about 2.5 to 4.5 R(E)

Quasi-monochromatic waves at about the hydrogen cyclotron frequency were observed as the ISEE 1 satellite traversed auroral field lines at radial distances of about 2.5-4.5 R(E) near midnight on June 19, 1981. Waves and both lower and higher frequencies were observed at higher altitudes, and possible electrostatic helium cyclotron and oxygen cyclotron waves occurred at lower altitudes. Upflowing hydrogen and oxygen beams and field-aligned currents occurred simultaneously. The features of the waves are most consistent with the current-driven mode. In addition, numerical studies of the linear dispersion relation, using parameters based on the observations, show that both the parallel and oblique two-stream modes and the ion-beam-driven modes were stable while oblique current-driven modes were unstable. The O(+) and H(+) distributions provide evidence for interactions with local electrostatic ion cyclotron waves and for the H(+)-O(+) two-stream instability at altitudes below the satellite.

Catell, C. A.↗

NASA NTRS1992

Computing high-speed flows past an oscillating cylinder near a vertical wall

A computational method to simulate unsteady flows involving moving rigid boundaries and interference has been developed. The method is used to solve inviscid equations governing the fluid flow and the dynamic equations governing the motion of rigid bodies. A second-order accurate, upwind-biased, and alterating-direction-implicit method is used to solve the governing equations of the flow. A kinematic domain decomposition (KDD) procedure is extended to treat 3D problems with a high degree of accuracy and generality. The method under consideration is applied to both transonic and supertransonic flows. Both cases involve flow past a cylinder which is forced to pitch sinusoidally near a vertical wall. Benefits of the proposed approach include accurate calculation of the flow around 3D moving multiple bodies with interference; reduction of a numerical error; in particular, the dispersion error which strongly affects wave propagation; and minimization of the phase error which is accumulated according to the time advance procedure.

Yen, Guan-Wei↗

NASA NTRS1993

Strongly nonlinear evolution of low-frequency wave packets in a dispersive plasma

The evolution of strongly nonlinear, strongly modulated wave packets is investigated in a dispersive plasma using a hybrid numerical code. These wave packets have amplitudes exceeding the strength of the external magnetic field, along which they propagate. Alfven (left helicity) wave packets show strong steepening for p < 1, while fast (fight heIicity) wave packets hardly steepen for any beta. Substantial regions of opposite helicity form on the leading side of steepened Alfven wave packets. This behavior differs qualitatively from that exhibited by the solutions to the derivative nonlinear Schrodinger (DNLS) equation.

Vasquez, Bernard J.↗

NASA NTRS2023

Oblique Propagation and Temperature Effects on the Resonant Right-Hand Ion Beam Instability

The resonant right-hand instability (RHI) is often the dominant mode driven by reflected ions upstream of Earth’s quasi-parallel bow shock. In the tradition of Peter Gary, this paper further explores the right-hand instability using numerical solutions of the plasma dispersion relation and non-linear kinetic simulations, with parameters inspired by observations from NASA’s Magnetospheric Multiscale (MMS) mission. Agreement is found between the ion distributions in the particle-in-cell simulations and Magnetospheric Multiscale spacecraft data, which show the gyrophase bunching characteristic of the instability. The non-linear structures created by right-hand instability tend to be stronger when the plasma beta is lower. These structures have sizes of around 100 to 200 ion inertial lengths perpendicular to the magnetic field, presenting planet-sized disturbances to the magnetosphere. 2d and 3D hybrid particle-in-cell simulations show that modes with a range of propagation angles oblique to the magnetic field are excited, providing a ground to understand previous statistical studies of observed foreshock waves.

Ari Le↗

DOE OSTI2023

An evaluation of meteorological data prediction over Washington, D.C.: Comparison of DCNet observations and NAM model outputs

This study presents an example of how outputs of operational and readily-available mesoscale numerical models can be adapted to initialize dispersion calculations within the urban surface roughness layer. Three years of urban meteorological observations from central Washington, DC, are compared against forecast outputs of the North American Mesoscale (NAM) model. NAM wind speed predictions underestimate the observations in light winds and overestimate the measurements in high winds. Average wind directions are consistent. However, an adjustment of the predicted direction of the plume by –20° is needed. The uncertainty associated with this adjustment is large in light NAM wind speed with no evident variation by season. The values of the standard deviation of the wind direction, σθ derived from NAM model outputs underestimate the observations by a small amount (about –1.5 to –2.5°). The results presented here indicate that mesoscale numerical model outputs can provide information adequate for dispersion calculations. However, levels of uncertainty associated with implementation of the suggested procedures increase with decreasing wind speed, causing considerable uncertainty in the implementation of adjustments as low wind speed conditions are approached. Results and recommendations reported here should not be extended to other numerical models or other cities without further testing.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1993

Dispersed ion structures at the poleward edge of the auroral oval: Low-altitude observations and numerical modeling

We have compared the AUREOL 3 (A3) observations of auroral ion precipitation, particularly ion beams, with the results from the global kinetic model of magnetotail plasma of Ashour-Abdalla et al. (1993). We have identified 101 energetic keV H(+) velocity dispersed precipitating ion structures (VDIS) with fluxes above 10(exp -3) ergs./sq cm./s in the A3 record between the end of 1981 and mid-1984. These beams display a systematic increase in energy with increasing latitude and were observed in a narrow region within less than 1 deg in latitude of the polar cap boundary. The VDIS are the most distinctive feature in the auroral zone of the plasma sheet boundary layer. We report first on a statistical analysis of the possible ralationships between magnetic activity or substorm phase and the VDIS properties. Our particle simulations of the precipitating ions have been extended by using a series of modified versions of the Tsyganenko (1989) magnetic field model and by varying the cross-magnetosphere electric field. In the simulations, plasma from a mantle source is subject to strong nonlinear acceleration, forming beams which flow along the PSBL. Only 3 to 4% of these beams precipitate into the ionosphere to form the VDIS while the majority return to the equatorial plane after mirroring and form the thermalized central plasma sheet. The final energy and the dispersion of the beams in the model depend on the amplitude of the cross-tail electric field. Two unsual observations of low-energy (less than 5 keV) O(+) VDIS, shifted by 4 deg 5 deg in invariant latitude equatorward of H(+) VDIS are analyzed in detail. The sparsity of such O(+) events and the absence of the changes in the flux and frequency of occurrence indicate a solar wind origin for the plasma. Finally, large-scale kinetic modeling, even with its simplifications and assumptions (e.g., static magnetic field, solar wind source), reproduces low-altitude auroral ion features fairly well; it may therefore be presented as an appropriate framework into which data on energization and transport of the hot plasma, obtained in the equatorial plane, could be inserted in the near future.

Bosqued, J. M.↗

DOE OSTI2021

Tracing the Mott-Hubbard transition in one-dimensional Hubbard models without Umklapp scattering

Here we apply the density-matrix renormalization group (DMRG) method to a one-dimensional Hubbard model that lacks Umklapp scattering and thus provides an ideal case to study the Mott-Hubbard transition analytically and numerically. The model has a linear dispersion and displays a metal-to-insulator transition when the Hubbard interaction U equals the bandwidth U c = W, where the single-particle gap opens linearly, Δ(U ≥ W) = U - W. The simple nature of the elementary excitations permits us to determine numerically with high accuracy the critical interaction strength and the gap function in the thermodynamic limit. The jump discontinuity of the momentum distribution $n_k$ at the Fermi wave number $k_F$ = 0 cannot be used to locate accurately U c from finite-size systems. However, the slope of $n_k$ at the band edges $k_B$ = ±π, reveals the formation of a single-particle bound state which can be used to determine U c reliably from $n_k$ using accurate finite-size data.

1-dimensional systems↗

NASA NTRS1997

Computations of Complex Three-Dimensional Turbulent Free Jets

Three-dimensional, incompressible turbulent jets with rectangular and elliptical cross-sections are simulated with a finite-difference numerical method. The full Navier- Stokes equations are solved at low Reynolds numbers, whereas at high Reynolds numbers filtered forms of the equations are solved along with a sub-grid scale model to approximate the effects of the unresolved scales. A 2-N storage, third-order Runge-Kutta scheme is used for temporary discretization and a fourth-order compact scheme is used for spatial discretization. Although such methods are widely used in the simulation of compressible flows, the lack of an evolution equation for pressure or density presents particular difficulty in incompressible flows. The pressure-velocity coupling must be established indirectly. It is achieved, in this study, through a Poisson equation which is solved by a compact scheme of the same order of accuracy. The numerical formulation is validated and the dispersion and dissipation errors are documented by the solution of a wide range of benchmark problems. Three-dimensional computations are performed for different inlet conditions which model the naturally developing and forced jets. The experimentally observed phenomenon of axis-switching is captured in the numerical simulation, and it is confirmed through flow visualization that this is based on self-induction of the vorticity field. Statistical quantities such as mean velocity, mean pressure, two-point velocity spatial correlations and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stress equations are presented to aid in the turbulence modeling of complex jets. Simulations of circular jets are used to quantify the effect of the non-uniform curvature of the non-circular jets.

Wilson, Robert V.↗

NASA NTRS1989

Assessment of numerical methods for the solution of fluid dynamics equations for nonlinear resonance systems

The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

Przekwas, A. J.↗

DOE OSTI2021

Measurement and control of optical nonlinearities in dispersive dielectric multilayers

Dispersive dielectric multilayer mirrors, high-dispersion chirped mirrors in particular, are widely used in modern ultrafast optics to manipulate spectral chirps of ultrashort laser pulses. Dispersive mirrors are routinely designed for dispersion compensation in ultrafast lasers and are assumed to be linear optical components. In this work, we report the experimental characterization of an unexpectedly strong nonlinear response in these chirped mirrors. At modest peak intensities <2 TW/cm 2 —well below the known laser-induced damage threshold of these dielectric structures—we observed a strong reflectivity decrease, local heating, transient spectral modifications, and time-dependent absorption of the incident pulse. Through computational analysis, we found that the incident laser field can be enhanced by an order of magnitude in the dielectric layers of the structure. The field enhancement leads to a wavelength-dependent nonlinear absorption, that shows no signs of cumulative damage before catastrophic failure. The nonlinear absorption is not a simply two-photon process but instead is likely mediated by defects that facilitate two-photon absorption. To mitigate this issue, we designed and fabricated a dispersive multilayer design that strategically suppresses the field enhancement in the high-index layers, shifting the high-field regions to the larger-bandgap, low-index layers. This strategy significantly increases the maximum peak intensity that the mirror can sustain. However, our finding of an onset of nonlinear absorption even at ‘modest’ fluence and peak intensity has significant implications for numerous past published experimental works employing dispersive mirrors. Additionally, our results will guide future ultrafast experimental work and ultrafast laser design.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRSDate not supplied

The Spacecraft Mass Balance as a Diagnostic Tool for Cabin Air Quality

The spacecraft-level mass balance is a powerful tool for diagnosing and decoupling factors influencing cabin air quality, especially with respect to trace contaminant propagation. Successful implementation of this approach relies upon accurate and temporally relevant air quality measurements, which are often challenging to attain for emerging trace contaminants. In addition to the proper interpretation of air quality data trends, a thorough understanding of subsystem-level mass transfer is required to characterize process performance. For many Environmental Control and Life Support System processes, subsystem mass transfer may be predicted based on physicochemical properties and classic unit operation design approaches. Gaps in understanding of process performance should be supplemented with thorough ground testing. A high confidence in the fidelity of one of these two aspects, air quality data or process performance, may help offset uncertainties in the other. At minimum, thresholding assumptions and numerical methods can be imposed to constrain unknown parameters within the physical envelope of the situation at hand. Finally, the complexity of cabin-wide integrated mass balances may be reduced by implementing a phenomenological approach towards subsystem discretization.

Cabin Air-quality↗

DOE OSTI2021

On the Reliability of Parameter Inferences in a Multiscale Model for Transport in Stream Corridors

Nonreacting tracer tests capture information about physical processes in transient storage zones including the hyporheic zone (HZ). However, reliably extracting this information from breakthrough curves (BTCs) and distinguishing the effects of in-channel dispersion and transient storage are well-known challenges. Using BTCs from a nonreacting tracer test monitored at multiple locations, we explore ways for reliable parameter estimations. The identifiability of parameters is greatly influenced by the choice of forward and inverse modeling frameworks in addition to the quality of the data. Our forward model is a recently proposed multiscale model that uses subgrid transport models written in the Lagrangian form to represent transport along a diverse set of HZ pathways with a shape-free distribution of travel times. Joint distributions of HZ and channel parameters are estimated using the Markov Chain Monte Carlo technique. Numerical experiments show ambiguity between channel dispersion and HZ transport when the reach length is too short to allow significant solute-HZ interaction, the observation period is too brief to observe the tailing behavior, or the solute source is spread in time. In contrast, we obtained reliable parameter estimates by simultaneously fitting BTCs observed at different locations in the test reach using a single set of HZ parameters and section-specific channel areas and dispersion coefficients. Furthermore, this study demonstrates the estimation of travel time distributions, HZ exchange rates, and channel parameters in a new multiscale approach and offers guidance for extracting reliable parameter estimates from multiple BTCs.

54 ENVIRONMENTAL SCIENCES↗