Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “astrophysical fluid dynamics”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

X-Ray Spectra from General Relativistic Radiation Magnetohydrodynamic Simulations of Thin Disks

We compare X-ray emission from several general relativistic, multifrequency, radiation magnetohydrodynamic simulations of thin black hole accretion disks with different accretion rates and spins. The simulations were performed using the M1 closure scheme, resolved with 12 frequency (energy) bins logarithmically spaced from 5 × 10 −3 to 5 × 10 3 ​​​ keV. We apply a general relativistic Monte Carlo transport code to postprocess the simulation data with greater fidelity in frequency resolution and Compton scattering treatment. Despite the relatively few energy bins and Kompaneets approximation to Compton scattering utilized in the M1 method, we find generally good agreement between the methods. Both produce prominent thermal profiles with peaks around 2–2.5 keV, where agreement is particularly strong and representative of the soft state. Both also find weaker (lower-luminosity) thermally sourced emission extending out to 100 keV due to the hotter innermost regions of the disks. Inverse Compton scattering becomes increasingly effective at hardening spectral outputs with increasing black hole spin, and becomes the dominant mechanism for photons that escape with energies between 10 to several hundred keV. At very high rates of spin, the radiation flux in this upscattered component becomes comparable to the thermal flux, a phenomenon typically associated with intermediate states. Beyond 10 4 keV, we observe faint, free–free emission from hot, optically thin coronal regions developing near the horizon, common to both spinning and nonspinning black holes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The Multiple Paths to Merger of Unequal-mass Black Hole Binaries in the Disks of Active Galactic Nuclei

The accretion disks that power active galactic nuclei (AGN) are thought to house populations of stars and compact objects; after forming binaries these compact objects may merge, begetting gravitational waves (GWs) such as those detected by LIGO and VIRGO. We present a comprehensive study of the early evolution of binaries within AGN disks as their orbits are influenced by the surrounding gas, focusing on eccentric and unequal-mass binaries. Nearly equal-mass binaries behave similarly to their equal-mass counterparts: Prograde binaries inspiral, albeit somewhat slowly, and have their eccentricities damped; retrograde binaries inspiral ∼2–3 times faster than their prograde counterparts, and those with near-equal masses are driven quickly toward near-unity eccentricities. However, the primaries in retrograde binaries with mass ratios of m 2 /m 1 ≲ 0.4 experience significantly weaker headwinds and retain substantial accretion disks that help damp binary eccentricities, slowing binary inspirals. Additionally, we find that while accretion drives prograde binaries toward equal masses thanks to the exchange of material between the primary and secondary accretion disks, retrograde binaries are driven slowly toward more extreme mass ratios. Prograde binaries, and generally those with low mass ratios, likely accrete for multiple e-folding timescales before merger. On the other hand, high-mass-ratio retrograde binaries may merge before accreting substantially, potentially approaching merger with detectable eccentricity. Future ground-based GW observatories, with their broader frequency coverage, should be particularly useful for studying these populations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1997

Pentium Pro inside: A treecode at 430 Gigaflops on ASCI Red - 1

As an entry for the 1997 Gordon Bell performance prize, we present results from two methods of solving the gravitational N-body problem on the Intel Teraflops system at Sandia National Laboratory (ASCI Red). The first method, an O(N2) algorithm, obtained 635 Gigaflops for a 1 million particle problem on 6800 Pentium Pro processors. The second solution method, a tree-code which scales as O(N log N), sustained 170 Gigaflops over a continuous 9.4 hour period on 4096 processors, integrating the motion of 322 million mutually interacting particles in a cosmology simulation, while saving over 100 Gigabytes of raw data. Additionally, the tree-code sustained 430 Gigaflops on 6800 processors for the first 5 time-steps of that simulation. This tree-code solution is approximately 105 times more efficient than the O(N2) algorithm for this problem. As an entry for the 1997 Gordon Bell price/performance prize, we present two calculations from the disciplines of astrophysics and fluid dynamics. The simulations were performed on two 16 Pentium Pro processor Beowulf-class computers (Loki and Hyglac) constructed entirely from commodity personal computer technology, at a cost of roughly $50k each in September, 1996. The price of an equivalent system in August 1997 is less than $30. At Los Alamos, Loki performed a gravitational tree-code N-body simulation of galaxy formation using 9.75 million particles, which sustained an average of 879 Mflops over a ten day period, and produced roughly 10 Gbytes of raw data.

Warren, M. S.↗
DOE OSTI2023

Modeling Solids in Nuclear Astrophysics with Smoothed Particle Hydrodynamics

Abstract Smoothed particle hydrodynamics (SPH) is a frequently applied tool in computational astrophysics to solve the fluid dynamics equations governing the systems under study. For some problems, for example when involving asteroids and asteroid impacts, the additional inclusion of material strength is necessary in order to accurately describe the dynamics. In compact stars, that is white dwarfs and neutron stars, solid components are also present. Neutron stars have a solid crust, which is the strongest material known in nature. However, their dynamical evolution, when modeled via SPH or other computational fluid dynamics codes, is usually described as a purely fluid dynamics problem. Here, we present the first 3D simulations of neutron star crustal toroidal oscillations including material strength with the Los Alamos National Laboratory SPH code FleCSPH. In the first half of the paper, we present the numerical implementation of solid material modeling together with standard tests. The second half is on the simulation of crustal oscillations in the fundamental toroidal mode. Here, we dedicate a large fraction of the paper to approaches that can suppress numerical noise in the solid. If not minimized, the latter can dominate the crustal motion in the simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1987

Spacelab 3 Mission Science Review

Papers and abstracts of the presentations made at the symposium are given as the scientific report for the Spacelab 3 mission. Spacelab 3, the second flight of the National Aeronautics and Space Administration's (NASA) orbital laboratory, signified a new era of research in space. The primary objective of the mission was to conduct applications, science, and technology experiments requiring the low-gravity environment of Earth orbit and stable vehicle attitude over an extended period (e.g., 6 days) with emphasis on materials processing. The mission was launched on April 29, 1985, aboard the Space Shuttle Challenger which landed a week later on May 6. The multidisciplinary payload included 15 investigations in five scientific fields: material science, fluid dynamics, life sciences, astrophysics, and atmospheric science.

Fichtl, George H.↗
NASA NTRS1987

Supercomputing in Aerospace

Topics addressed include: numerical aerodynamic simulation; computational mechanics; supercomputers; aerospace propulsion systems; computational modeling in ballistics; turbulence modeling; computational chemistry; computational fluid dynamics; and computational astrophysics.

Kutler, Paul↗
NASA NTRS2003

Finite Element Method for Capturing Ultra-relativistic Shocks

While finite element methods are used extensively by researchers solving computational fluid dynamics in fields other than astrophysics, their use in astrophysical fluid simulations has been predominantly overlooked. Current simulations using other methods such as finite difference and finite volume (based on finite difference) have shown remarkable results, but these methods are limited by their fundamental properties in aspects that are important for simulations with complex geometries and widely varying spatial and temporal scale differences. We have explored the use of finite element methods for astrophysical fluids in order to establish the validity of using such methods in astrophysical environments. We present our numerical technique applied to solving ultra-relativistic (Lorentz Factor Gamma >> 1) shocks which are prevalent in astrophysical studies including relativistic jets and gamma-ray burst studies. We show our finite element formulation applied to simulations where the Lorentz factor ranges up to 2236 and demonstrate its stability in solving ultra-relativistic flows. Our numerical method is based on the Flowfield Dependent Variation (FDV) Method, unique in that numerical diffusion is derived from physical parameters rather than traditional artificial viscosity methods. Numerical instabilities account for most of the difficulties when capturing shocks in this regime. Our method results in stable solutions and accurate results as compared with other methods.

Richardson, G. A.↗
DOE OSTI2023

Fluid dynamic mathematical aspects of supernova remnants

Supernovae—explosions of stars—are a central problem in astrophysics since they contain information on the entire process of stellar evolution and nucleosynthesis. Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities, developing during the supernova blast, lead to intense interfacial RT/RM mixing of the star's materials and couple astrophysical to atomic scales. This work analyzes some fluid dynamic mathematical aspects of the titanic task of supernova's blast. We handle mathematical challenges of RT/RM dynamics in supernova relevant conditions by directly linking the conservation laws governing RT/RM dynamics to symmetry-based momentum model, by exactly deriving the model parameters in the scale-dependent and scale-invariant regimes, and by exploring the special self-similar class for RT/RM interfacial mixing with variable accelerations. Here we reveal that RT/RM dynamics is strongly influenced by deterministic (the initial and the flow) conditions in the scale-dependent linear and nonlinear regimes and in the self-similar mixing regime. The theory outcomes are consistent with the observations of supernova remnants, explain the results of the scaled laboratory experiments in high energy density plasmas, and yield the design of future experiments for the accurate quantification of RT/RM dynamics in supernova relevant conditions. We find that from fluid dynamic mathematical perspectives, supernovae can be regarded as an astrophysical initial value problem. Along with the guidance of what explodes at microscopic scales, supernova remnants encapsulate information on the explosion hydrodynamics and the associated deterministic conditions at macroscopic scales. We urge such effects be considered in interpretations of the observational data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Spectra of Magnetic Turbulence in a Relativistic Plasma

Abstract We present a phenomenological and numerical study of strong Alfvénic turbulence in a magnetically dominated collisionless relativistic plasma with a strong background magnetic field. In contrast with the nonrelativistic case, the energy in such turbulence is contained in magnetic and electric fluctuations. We argue that such turbulence is analogous to turbulence in a strongly magnetized nonrelativistic plasma in the regime of broken quasi-neutrality. Our 2D particle-in-cell numerical simulations of turbulence in a relativistic pair plasma find that the spectrum of the total energy has the scaling k −3/2 , while the difference between the magnetic and electric energies, the so-called residual energy, has the scaling k −2.4 . The electric and magnetic fluctuations at scale ℓ exhibit dynamic alignment with the alignment angle scaling close to cos ϕ ℓ ∝ ℓ 1 / 4 . At scales smaller than the (relativistic) plasma inertial scale, the energy spectrum of relativistic inertial Alfvén turbulence steepens to k −3.5 .

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2015

Marshall Space Flight Center Faculty Fellowship Program

The Faculty Fellowship program was revived in the summer of 2015 at NASA Marshall Space Flight Center, following a period of diminished faculty research activity here since 2006 when budget cuts in the Headquarters' Education Office required realignment. Several senior Marshall managers recognized the need to involve the Nation's academic research talent in NASA's missions and projects to the benefit of both entities. These managers invested their funds required to establish the renewed Faculty Fellowship program in 2015, a 10-week residential research involvement of 16 faculty in the laboratories and offices at Marshall. These faculty engineers and scientists worked with NASA collaborators on NASA projects, bringing new perspectives and solutions to bear. This Technical Memorandum is a compilation of the research reports of the 2015 Marshall Faculty Fellowship program, along with the Program Announcement (appendix A) and the Program Description (appendix B). The research touched on seven areas-propulsion, materials, instrumentation, fluid dynamics, human factors, control systems, and astrophysics. The propulsion studies included green propellants, gas bubble dynamics, and simulations of fluid and thermal transients. The materials investigations involved sandwich structures in composites, plug and friction stir welding, and additive manufacturing, including both strength characterization and thermosets curing in space. The instrumentation projects involved spectral interfero- metry, emissivity, and strain sensing in structures. The fluid dynamics project studied the water hammer effect. The human factors project investigated the requirements for close proximity operations in confined spaces. Another team proposed a controls system for small launch vehicles, while in astrophysics, one faculty researcher estimated the practicality of weather modification by blocking the Sun's insolation, and another found evidence in satellite data of the detection of a warm-hot intergalactic medium filament. Our goal is to continue the Faculty Fellowship effort with Center funds in succeeding summers.

Six, N. F.↗
DOE OSTI2023

Spatial Intermittency of Particle Distribution in Relativistic Plasma Turbulence

Abstract Relativistic magnetically dominated turbulence is an efficient engine for particle acceleration in a collisionless plasma. Ultrarelativistic particles accelerated by interactions with turbulent fluctuations form nonthermal power-law distribution functions in the momentum (or energy) space, f ( γ ) d γ ∝ γ − α d γ , where γ is the Lorenz factor. We argue that in addition to exhibiting non-Gaussian distributions over energies, particles energized by relativistic turbulence also become highly intermittent in space. Based on particle-in-cell numerical simulations and phenomenological modeling, we propose that the bulk plasma density has lognormal statistics, while the density of the accelerated particles, n , has a power-law distribution function, P ( n ) dn ∝ n − β dn . We argue that the scaling exponents are related as β ≈ α + 1, which is broadly consistent with numerical simulations. Non-space-filling, intermittent distributions of plasma density and energy fluctuations may have implications for plasma heating and for radiation produced by relativistic turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Turbulence and Particle Acceleration in a Relativistic Plasma

Abstract In a collisionless plasma, the energy distribution function of plasma particles can be strongly affected by turbulence. In particular, it can develop a nonthermal power-law tail at high energies. We argue that turbulence with initially relativistically strong magnetic perturbations (magnetization parameter σ ≫ 1) quickly evolves into a state with ultrarelativistic plasma temperature but mildly relativistic turbulent fluctuations. We present a phenomenological and numerical study suggesting that in this case, the exponent α in the power-law particle-energy distribution function, f ( γ ) d γ ∝ γ − α d γ , depends on magnetic compressibility of turbulence. Our analytic prediction for the scaling exponent α is in good agreement with the numerical results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Theories of Relativistic Dissipative Fluid Dynamics

Relativistic dissipative fluid dynamics finds widespread applications in high-energy nuclear physics and astrophysics. However, formulating a causal and stable theory of relativistic dissipative fluid dynamics is far from trivial; efforts to accomplish this reach back more than 50 years. In this review, we give an overview of the field and attempt a comparative assessment of (at least most of) the theories for relativistic dissipative fluid dynamics proposed until today and used in applications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A simple model of a strong shock driven by a spherical or cylindrical piston

A simple model of piston-driven spherical and cylindrical shocks is suggested. The model is based on a consistent use of two factors: (a) an almost uniform pressure across the shocked layer and, (b) continuous geometrical stretching of the surface elements of the expanding piston. It turns out that for a uniform pre-shock medium the gas between the piston and the shock behaves essentially as an incompressible fluid. An algebraic equation for the shock vs piston position is obtained. Detailed evaluation of the accuracy of the proposed solution shows that its accuracy is a few percent for the adiabatic index γ = 5/3 (as in ideal plasma). A closed-form solution describing enhancement of a weak ambient magnetic field by the shock is presented. The proposed model of piston-driven shocks goes beyond the classical self-similar solutions in that it: (1) naturally covers an early, non-asymptotic dynamics and its transition to asymptotic regime; (2) allows for smooth radial density variation of an ambient gas of the form of bumps, dimples or ramps between two constant values; and (3) allows for smooth temporal variation of piston velocity of the form of bumps, dimples, or ramps. Furthermore, this simple and versatile model provides some new insights into a classical hydrodynamical problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

A comparative study of computational methods in cosmic gas dynamics

Many theoretical investigations of fluid flows in astrophysics require extensive numerical calculations. The selection of an appropriate computational method is, therefore, important for the astronomer who has to solve an astrophysical flow problem. The present investigation has the objective to provide an informational basis for such a selection by comparing a variety of numerical methods with the aid of a test problem. The test problem involves a simple, one-dimensional model of the gas flow in a spiral galaxy. The numerical methods considered include the beam scheme, Godunov's method (G), the second-order flux-splitting method (FS2), MacCormack's method, and the flux corrected transport methods of Boris and Book (1973). It is found that the best second-order method (FS2) outperforms the best first-order method (G) by a huge margin.

Van Albada, G. D.↗
DOE OSTI2025

A review of high order strong stability preserving two-derivative explicit, implicit, and IMEX methods

High order strong stability preserving time discretizations ensure the nonlinear non-inner-product strong stability properties of spatial discretizations suited for the stable simulation of hyperbolic PDEs in a wide variety of application areas including fluid dynamics, magnetohydrodynamics, semiconductor devices, electromagnetics, and astrophysics. Over the past decade multiderivative time-stepping have been increasingly used for the time-evolution hyperbolic PDEs, so that the strong stability properties of these methods have become important. In this work we review sufficient conditions for a two-derivative multistage method to preserve the strong stability properties of spatial discretizations in a forward Euler and different conditions on the second derivative. In particular we present the strong stability preserving theory for explicit and implicit two-derivative Runge–Kutta schemes, including a special condition on the second derivative under which these implicit methods may be unconditionally strong stability preserving. This special condition is natural for the stiff component of wide range of plasma physics problems, and can be useful in the context of strong stability preserving implicit-explicit multi-derivative Runge–Kutta schemes, where the time-step restriction is then independent of the stiff term. Lastly, we present the strong stability preserving theory for implicit-explicit multi-derivative general linear methods, and some novel second and third order methods where the time-step restriction is independent of the stiff term.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Local well-posedness and singularity formation in non-Newtonian compressible fluids

We investigate the initial value problem of a very general class of 3 + 1 non-Newtonian compressible fluids in which the viscous stress tensor with shear and bulk viscosity relaxes to its Navier–Stokes values. These fluids correspond to the non-relativistic limit of well-known Israel–Stewart-like theories used in the relativistic fluid dynamic simulations of high-energy nuclear and astrophysical systems. After establishing the local well-posedness of the Cauchy problem, we show for the first time in the literature that there exists a large class of initial data for which the corresponding evolution breaks down in finite time due to the formation of singularities. This implies that a large class of non-Newtonian fluids do not have finite solutions defined at all times.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗