Theoretical models and numerical studies of waves in a three-fluid medium
Mathematical models of wave propagation in three-fluid plasma and numerical solutions to dispersion and propagation properties
Murphree, D. L.↗
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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Computational procedures for mixed equations with shock waves
This paper discusses the procedures we have developed to treat a canonical problem involving a mixed nonlinear equation with boundary data that imply a discontinuous solution. This equation arises in various physical contexts and is basic to the description of the nonlinear acoustic behavior of a shock wave near a caustic. The numerical scheme developed is of second order, treats discontinuities as such by applying the appropriate jump conditions across them, and eliminates the numerical dissipation and dispersion associated with large gradients. Our results are compared with the results of a first-order scheme and with those of a second-order scheme we have developed. The algorithm used here can easily be generalized to more complicated problems, including transonic flows with imbedded shocks.
Yu, N. J.↗
The influence of the directional energy distribution on the nonlinear dispersion relation in a random gravity wave field
The influence of the directional distribution of wave energy on the dispersion relation is calculated numerically using various directional wave spectrum models. The results indicate that the dispersion relation varies both as a function of the directional energy distribution and the direction of propagation of the wave component under consideration. Furthermore, both the mean deviation and the random scatter from the linear approximation increase as the energy spreading decreases. Limited observational data are compared with the theoretical results. The agreement is favorable.
Huang, N. E.↗
On the stability of nongyrotropic ion populations - A first (analytic and simulation) assessment
The wave and dispersion equations for perturbations propagating parallel to an ambient magnetic field in magnetoplasmas with nongyrotropic ion populations show, in general, the occurrence of coupling between the parallel (left- and right-hand circularly polarized electromagnetic and longitudinal electrostatic) eigenmodes of the associated gyrotropic medium. These interactions provide a means to driving linearly one mode with free-energy sources of other modes in homogeneous media. Different types of nongyrotropy bring about distinct classes of coupling. The stability of a hydrogen magnetoplasma with anisotropic, nongyrotropic protons that only couple the electromagnetic modes to each other is investigated analytically (via solution of the derived dispersion equation) and numerically (via simulation with a hybrid code). Nongyrotropy enhances growth and enlarges the unstable spectral range relative to the corresponding gyrotropic situation. The relevance of the properties of nongyrotropic populations to space plasma environments is also discussed.
Brinca, A. L.↗
Application of the Space-Time Conservation Element and Solution Element Method to One-Dimensional Advection-Diffusion Problems
Test problems are used to examine the performance of several one-dimensional numerical schemes based on the space-time conservation and solution element (CE/SE) method. Investigated in this paper are the CE/SE schemes constructed previously for solving the linear unsteady advection-diffusion equation and the schemes derived here for solving the nonlinear viscous and inviscid Burgers equations. In comparison with the numerical solutions obtained using several traditional finite-difference schemes with similar accuracy, the CE/SE solutions display much lower numerical dissipation and dispersion errors.
Wang, Xiao-Yen↗
Generalized Langmuir Waves in Magnetized Kinetic Plasmas
The properties of unmagnetized Langmuir waves and cold plasma magnetoionic waves (x, o, z and whistler) are well known. However, the connections between these modes in a magnetized kinetic plasma have not been explored in detail. Here, wave properties are investigated by numerically solving the dispersion equation derived from the Vlasov equations both with and without a beam instability present. For omega(sub p)>Omega(sub e), it is shown that the generalized Langmuir mode at oblique propagation angles has magnetic z-mode characteristics at low wave numbers and thermal Langmuir mode characteristics at high wave numbers. For omega(sub p)<Omega(sub e), it is shown that the (oblique) Langmuir mode instead connects to the whistler mode at low wave numbers. The transition from the Langmuir/z mode to the Langmuir/whistler mode near omega(sub p) = Omega(sub e) is rapid. In addition, the effects on wave dispersion and polarization after adding a beam are investigated. Applications of this theory to magnetized Langmuir waves in Earth's foreshock and the solar wind, to waves observed near the plasma frequency in the auroral regions, and to solar type III bursts are discussed.
Willes, A. J.↗
On Efficient Multigrid Methods for Materials Processing Flows with Small Particles
Multiscale modeling of materials requires simulations of multiple levels of structural hierarchy. The computational efficiency of numerical methods becomes a critical factor for simulating large physical systems with highly desperate length scales. Multigrid methods are known for their superior efficiency in representing/resolving different levels of physical details. The efficiency is achieved by employing interactively different discretizations on different scales (grids). To assist optimization of manufacturing conditions for materials processing with numerous particles (e.g., dispersion of particles, controlling flow viscosity and clusters), a new multigrid algorithm has been developed for a case of multiscale modeling of flows with small particles that have various length scales. The optimal efficiency of the algorithm is crucial for accurate predictions of the effect of processing conditions (e.g., pressure and velocity gradients) on the local flow fields that control the formation of various microstructures or clusters.
Thomas, James↗
Hyperbolic Method for Dispersive PDEs: Same High-Order of Accuracy for Solution, Gradient, and Hessian
In this paper, we introduce a new hyperbolic first-order system for general dispersive partial differential equations (PDEs). We then extend the proposed system to general advection-diffusion-dispersion PDEs. We apply the fourth-order RD scheme of Ref. 1 to the proposed hyperbolic system, and solve time-dependent dispersive equations, including the classical two-soliton KdV and a dispersive shock case. We demonstrate that the predicted results, including the gradient and Hessian (second derivative), are in a very good agreement with the exact solutions. We then show that the RD scheme applied to the proposed system accurately captures dispersive shocks without numerical oscillations. We also verify that the solution, gradient and Hessian are predicted with equal order of accuracy.
Mazaheri, Alireza↗
Exascale Multiphysics Nuclear Reactor Simulations for Advanced Designs
ENRICO is a coupled application developed under the U.S. Department of Energy's Exascale Computing Project (ECP) targeting the modeling of advanced nuclear reactors. It couples radiation transport with heat and fluid simulation, including the high-fidelity, highresolution Monte-Carlo code Shift and the Computational fluid dynamics code NekRS. NekRS is a highly-performant open-source code for simulation of incompressible and low-Mach fluid flow, heat transfer, and combustion with a particular focus on turbulent flows in complex domains. It is based on rapidly convergent high-order spectral element discretizations that feature minimal numerical dissipation and dispersion. State-of-the-art multilevel preconditioners, efficient high-order time-splitting methods, and runtime-adaptive communication strategies are built on a fast OCCA-based kernel library, libParanumal, to provide scalability and portability across the spectrum of current and future high-performance computing platforms. On Frontier, Nek5000/RS has recently achieved an unprecedented milestone in breaching over 1 billion spectral elements and 350 billion degrees of freedom. Shift has demonstrated the capability to transport upwards of 1 billion particles per second in full core nuclear reactor simulations featuring complete temperature-dependent, continuous-energy physics on Frontier. Shift achieved a weak-scaling efficiency of 97.8% on 8192 nodes of Frontier and calculated 6 reactions in 214,896 fuel pin regions below 1% statistical error yielding first-of-a-kind resolution for a Monte Carlo transport application.
Hamilton, Steven P.↗
Investigating Sources of Variability and Error in Simulations of Carbon Dioxide in an Urban Region
Greenhouse gas (GHG) emissions estimation methods that use atmospheric trace gas observations, including inverse modeling techniques, perform better when carbon dioxide (CO2) fluxes are more accurately transported and dispersed in the atmosphere by a numerical model. In urban areas, transport and dispersion is particularly difficult to simulate using current mesoscale meteorological models due, in part, to added complexity from surface heterogeneity and fine spatial/temporal scales. It is generally assumed that the errors in GHG estimation methods in urban areas are dominated by errors in transport and dispersion. Other significant errors include, but are not limited to, those from assumed emissions magnitude and spatial distribution. To assess the predictability of simulated trace gas mole fractions in urban observing systems using a numerical weather prediction model, we employ an Eulerian model that combines traditional meteorological variables with multiple passive tracers of atmospheric CO2 from anthropogenic inventories and a biospheric model. The predictability of the Eulerian model is assessed by comparing simulated atmospheric CO2 mole fractions to observations from four in situ tower sites (three urban and one rural) in the Washington DC/Baltimore, MD area for February 2016. Four different gridded fossil fuel emissions inventories along with a biospheric flux model are used to create an ensemble of simulated atmospheric CO2 observations within the model. These ensembles help to evaluate whether the modeled observations are impacted more by the underlying emissions or transport. The spread of modeled observations using the four emission fields indicates the model's ability to distinguish between the different inventories under various meteorological conditions. Overall, the Eulerian model performs well; simulated and observed average CO2 mole fractions agree within 1% when averaged at the three urban sites across the month. However, there can be differences greater than 10% at any given hour, which are attributed to complex meteorological conditions rather than differences in the inventories themselves. On average, the mean absolute error of the simulated compared to actual observations is generally twice as large as the standard deviation of the modeled mole fractions across the four emission inventories. This result supports the assumption, in urban domains, that the predicted mole fraction error relative to observations is dominated by errors in model meteorology rather than errors in the underlying fluxes in winter months. As such, minimizing errors associated with atmospheric transport and dispersion may help improve the performance of GHG estimation models more so than improving flux priors in the winter months. We also find that the errors associated with atmospheric transport in urban domains are not restricted to certain times of day. This suggests that atmospheric inversions should use CO2 observations that have been filtered using meteorological observations rather than assuming that meteorological modeling is most accurate at certain times of day (such as using only mid-afternoon observations).
CO2↗
Explicit dispersion relations for warm fluid waves in a uniform plasma (invited)
Classical dispersion relations for waves in a fluid plasma are expressed as implicit functions of wave frequency and wave- number. The implicit dispersion relation must in general be solved numerically. This work introduces an astute method for obtaining an explicit dispersion relation for waves in a fluid plasma. The explicit expression has an advantage of providing eigenmodes without numerical computations and giving a dispersion relation more detailed than the implicit expression. As in the case of cold waves, the wavenumber for an arbitrary frequency can be directly obtained from the explicit formula. The explicit formula also enables an accurate evaluation of finite-temperature effects on a dispersion relation even near the cyclotron resonance, where a numerical analysis of the implicit relation is impractical. As a result, the analytic formula can be used to investigate temperature effects on electromagnetic waves.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
BASS. XXVI. DR2 Host Galaxy Stellar Velocity Dispersions
We present new central stellar velocity dispersions for 484 Sy 1.9 and Sy 2 from the second data release of the Swift/BAT AGN Spectroscopic Survey (BASS DR2). This constitutes the largest study of velocity dispersion measurements in X-ray-selected obscured active galactic nuclei (AGN) with 956 independent measurements of the Ca II H and K λ3969, 3934 and Mg I λ5175 region (3880–5550 Å) and the calcium triplet region (8350–8730 Å) from 642 spectra mainly from VLT/X-Shooter or Palomar/DoubleSpec. Our sample spans velocity dispersions of 40–360 km s 1 , corresponding to 4–5 orders of magnitude in black hole mass (M BH = 10 5.5-9.6 M ⊙ ), bolometric luminosity (L bol ~ 10 42–46 erg s -1 ), and Eddington ratio (L/L Edd ~ 10 -5 to 2). For 281 AGN, our data and analysis provide the first published central velocity dispersions, including six AGN with low-mass black holes (M BH = 10 5.5-6.5 M ⊙ ), discovered thanks to high spectral resolution observations (σ inst ~ 25 km s -1 ). The survey represents a significant advance with a nearly complete census of velocity dispersions of hard X-ray–selected obscured AGN with measurements for 99% of nearby AGN (z < 0.1) outside the Galactic plane (|b| > 10°). The BASS AGN have much higher velocity dispersions than the more numerous optically selected narrow-line AGN (i.e., ~150 versus ~100 km s -1 ) but are not biased toward the highest velocity dispersions of massive ellipticals (i.e., >250 km s -1 ). Despite sufficient spectral resolution to resolve the velocity dispersions associated with the bulges of small black holes (~10 4–5 M ⊙ ), we do not find a significant population of super-Eddington AGN. Using estimates of the black hole sphere of influence from velocity dispersion, direct stellar and gas black hole mass measurements could be obtained with existing facilities for more than ~100 BASS AGN.
79 ASTRONOMY AND ASTROPHYSICS↗
[Kelvin-Helmholtz instability in protostellar jets]
NASA grant NAG 5 2866, funded by the Astrophysics Theory Program, enabled the study the Kelvin-Helmholtz instability in protostellar jets. In collaboration with co-investigator Philip Hardee, the PI derived the analytic dispersion relation for the instability in including a cooling term in the energy equation which was modeled as one of two different power laws. Numerical solutions to this dispersion relation over a wide range of perturbation frequencies, and for a variety of parameter values characterizing the jet (such as Mach number, and density ratio) were found It was found that the growth rates and wavelengths associated with unstable roots of the dispersion relation in cooling jets are significantly different than those associated with adiabatic jets, which have been studied previously. In collaboration with graduate student Jianjun Xu (funded as a research associate under this grant), hydrodynamical simulations were used to follow the growth of the instability into the nonlinear regime. It was found that asymmetric surface waves lead to large amplitude, sinusoidal distortions of the jet, and ultimately to disruption Asymmetric body waves, on the other hand, result in the formation of shocks in the jet beam in the nonlinear regime. In cooling jets, these shocks lead to the formation of dense knots and filaments of gas within the jet. For sufficiently high perturbation frequencies, however, the jet cannot respond and it remains symmetric. Applying these results to observed systems, such as the Herbig-Haro jets HH34, HH111 and HH47 which have been observed with the Hubble Space Telescope, we predicted that some of the asymmetric structures observed in these systems could be attributed to the K-H modes, but that perturbations on timescales associated with the inner disk (about 1 year) would be too rapid to cause disruption. Moreover, it was discovered that weak shock 'spurs' in the ambient gas produced by ripples in the jet surface due to nonlinear, modes of surface and/or body waves could accelerate the ambient gas to low velocity. This latter effect represents a new mechanism by which supersonic jets can accelerate low velocity outflows.
Stone, James↗
Numerical simulation of advection fog formation on multi-disperse aerosols due to combustion-related pollutants
The effects of multi-disperse distribution of the aerosol population are presented. Single component and multi-component aerosol species on the condensation/nucleation processes which affect the reduction in visibility are described. The aerosol population with a high particle concentration provided more favorable conditions for the formation of a denser fog than the aerosol population with a greater particle size distribution when the value of the mass concentration of the aerosols was kept constant. The results were used as numerical predictions of fog formation. Two dimensional observations in horizontal and vertical coordinates, together with time-dependent measurements were needed as initial values for the following physical parameters: (1)wind profiles; (2) temperature profiles; (3) humidity profiles; (4) mass concentration of aerosol particles; (5) particle size distribution of aerosols; and (6) chemical composition of aerosols. Formation and dissipation of advection fog, thus, can be forecasted numerically by introducing initial values obtained from the observations.
Hung, R. J.↗
Electrostatic wave generation above and below the plasma frequency by electron beams
The unmagnetized electrostatic dispersion equation is solved numerically to study the growth of electrostatic waves near the plasma frequency resulting from a unstable electron beam. An attempt is made to find the dispersion relations, frequencies of maximum growth, and the resonant or nonresonant character of the waves for particular beam parameters. It is found that the unstable waves do not have the Langmuir dispersion relation except in the limit of a very dilute beam with growth on the connected mode of O'Neil and Malmberg. The properties of the unstable mode depend strongly on beam parameters such as beam density, speed, and temperature.
Cairns, Iver H.↗
Evolution of planetesimals. I - Dynamics: Relaxation in a thin disk. II - Numerical simulations
The study examines the effects of density inhomogeneity and differential rotation as well as inelastic collisions on the dynamical evolution of planetesimals. Consideration is given to a three-step analysis: the dynamical evolution of the planetesimals, collisions and mass accumulation, and interaction with gas. It is shown that the velocity dispersion of a cold system of planetesimals increases rapidly due to elastic gravitational scattering. When the dispersion in the epicycle amplitude becomes comparable to the planetesimals' Roche radius, energy is transferred from the systematic Keplerian shear to the dispersive motion. With a numerical N-body scheme, gravitational scattering and physical collisions among a system of planetesimals is simulated. It is shown that dynamical equilibrium is attained with a velocity dispersion comparable to the surface escape velocity of those planetesimals which contribute most of the system mass.
Palmer, P. L.↗
An analysis of scattered light in low dispersion IUE spectra
A detailed numerical simulation of light scattering from the low-resolution grating in the short wavelength spectrograph of the IUE Observatory was developed, in order to quantitatively analyze the effects of scattering on both continuum and line emission spectra. It is found that: (1) the redistribution of light by grating scattering did not appreciably alter either the shape or the absolute flux level of continuum spectra for A-F stars; (2) late-type stellar continua showed a tendency to flatten when observed in scattered light toward the shorter wavelengths; and (3) the effect of grating scattering on emission lines is to decrease measured line intensities by an increasing percentage toward the shorter wavelengths. The spectra obtained from scattering experiments for solar-type and late type stars are reproduced in graphic form.
Basri, G.↗
Solitary-wave interaction
The interaction of solitary-wave solutions of a model equation for long waves in dispersive media is examined numerically. It is found that the waves do not emerge from the interaction unscathed. Instead, two new solitary waves, having slightly different amplitudes from the original waves, together with a small dispersive tail are generated as a result of the interaction.
Bona, J. L.↗