Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamics methods”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics

In this work, we present a conservative numerical method for radiation magnetohydrodynamics with frequency-dependent full transport in stationary spacetimes. This method is stable and accurate for both large and small optical depths and radiation pressures. The radiation stress–energy tensor is evolved in flux-conservative form, and closed with a swarm of samples that each transport a multigroup representation of the invariant specific intensity along a null geodesic. In each zone, the enclosed samples are used to efficiently construct a Delaunay triangulation of the unit sphere in the comoving frame, which in turn is used to calculate the Eddington tensor, average source terms, and adaptively refine the sample swarm. Furthermore, radiation four-forces are evaluated in the moment sector in a semi-implicit fashion. The radiative transfer equation is solved in invariant form deterministically for each sample. Since each sample carries a discrete representation of the full spectrum, the cost of evaluating the transport operator is independent of the number of frequency groups, representing a significant reduction of algorithmic complexity for transport in frequency-dependent problems. The major approximation we make in this work is performing scattering in an angle-averaged way. Local adaptivity in samples also makes this scheme more amenable to nonuniform meshes than a traditional Monte Carlo method. We describe the method and present results on a suite of test problems. We find that Method of Characteristics Moment Closure converges at least as ~N -1 , rather than the canonical Monte Carlo N -1/2 , where N is the number of samples per zone.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1973

An economical method of analyzing transient motion of gas-lubricated rotor-bearing systems.

A method of economically evaluating the hydrodynamic forces generated in a gas-lubricated tilting-pad bearing is presented. The numerical method consists of solving the case of the infinite width bearing and then converting this solution to the case of the finite bearing by accounting for end leakage. The approximate method is compared to the finite-difference solution of Reynolds equation and yields acceptable accuracy while running about one-hundred times faster. A mathematical model of a gas-lubricated tilting-pad vertical rotor systems is developed. The model is capable of analyzing a two-bearing-rotor system in which the rotor center of mass is not at midspan by accounting for gyroscopic moments. The numerical results from the model are compared to actual test data as well as analytical results of other investigators.

Falkenhagen, G. L.↗
DOE OSTI2020

Efficient smoothed particle radiation hydrodynamics I: Thermal radiative transfer

This work presents efficient solution techniques for radiative transfer in the smoothed particle hydrodynamics discretization. Two choices that impact efficiency are how the material and radiation energy are coupled, which determines the number of iterations needed to converge the emission source, and how the radiation diffusion equation is solved, which must be done in each iteration. The coupled material and radiation energy equations are solved using an inexact Newton iteration scheme based on nonlinear elimination, which reduces the number of Newton iterations needed to converge within each time step. During each Newton iteration, the radiation diffusion equation is solved using Krylov iterative methods with a multigrid preconditioner, which abstracts and optimizes much of the communication when running in parallel. The code is verified for an infinite medium problem, a one-dimensional Marshak wave, and a two and three-dimensional manufactured problem, and exhibits first-order convergence in time and second-order convergence in space. For these problems, the number of iterations needed to converge the inexact Newton scheme and the diffusion equation is independent of the number of spatial points and the number of processors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Comparison of two pressure–temperature equilibration methods

We compare and contrast the traditionally used method of solving the pressure–temperature equilibration problem in hydrodynamics, where specific internal energy and density are considered independent variables, with a different method where pressure and temperature are independent variables. With the goal of examining the robustness of the two methods as the number of components increases, we examine 2-, 4-, 6-, and 8-component systems. After equilibrating more than 10 4 initial conditions for each system using both methods, we demonstrate that the latter method constrains the search space by lowering its dimensionality and forces a better initial guess, resulting in a higher probability of convergence to solution with fewer, cheaper iterations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Advancing our understanding of instabilities in high-energy-density systems through the study of benchmark datasets with advanced modeling tools

Numerous types of pulsed power driven inertial confinement fusion (ICF) and high energy density (HED) systems rely on implosion stability to achieve desired temperatures, pressures, and densities. Sandia National Laboratories Pulsed Power Sciences Center’s main ICF platform, Magnetized Liner Inertial Fusion (MagLIF), suffers from implosion instabilities which limit attainable fuel conditions and can compromise fuel confinement. This Truman Fellowship research primarily focused on computationally exploring (a) methods for improving our understanding of hydrodynamic and magnetohydrodynamic instabilities that form during cylindrical liner implosions, (b) methods for mitigating implosion instabilities, particularly those that degrade performance of MagLIF targets, and (c) novel MagLIF target designs intended to improve target performance primarily via enhanced implosion stability. Several multi-dimensional computational tools were used, including the magnetohydrodynamics code ALEGRA, the radiation-magnetohydrodynamics code HYDRA, and the magnetohydrodynamics code KRAKEN. This research succeeded in executing and analyzing simulations of automagnetizing liner implosions, shockless MagLIF implosions, dynamic screw pinch driven cylindrical liner implosions, and cylindrically convergent HED instability studies. The methods and tools explored and developed in this Truman Fellowship research have been published in several peer-reviewed journal articles and will serve as useful contributions to the fields of pulsed power science and engineering, particularly pertaining to pulsed power ICF and HED science.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Coupling laser physics to radiation-hydrodynamics

In order to accurately model implosion hydrodynamics in a radiation-hydrodynamics code, it is essential to include accurate accounting for energy deposition physics. In inertial confinement fusion (ICF), where capsules are driven by lasers or laser-driven x-rays, energy deposition profiles and energy transport have a strong impact on the development and evolution of capsule dynamics and hydrodynamic instabilities. Nevertheless, accurately modeling laser beam propagation in radiation-hydrodynamics codes presents unique challenges associated with disparate resolution requirements, the potential to seed spurious noise in highly unstable systems, and computational expense. We discuss a new method for coupling laser ray-tracing physics to a radiation hydrodynamics code, developed in the process of implementing the Mazinisin laser ray-trace into the xRAGE radiation hydrodynamics code. In contrast to previous approaches, in which laser ray-tracing is performed on the radiation-hydrodynamics mesh, our method involves a mesh generation and evolution strategy that addresses the unique requirements of the laser ray-trace in a separate mesh, enabling performance enhancements and strategies to reduce noise seeded by the discretization of beams into computational rays. In addition, we have employed several methods to ensure that spurious mesh imprinting is minimized. These involved optimizing the laser and radiation-hydrodynamics meshes as well as interpolation between them and requires the use of an exact initialization method for the radiation-hydrodynamics mesh. Furthermore, these techniques have enabled efficient computation of laser-driven implosions and other experiments with minimal introduction of spurious noise.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Observational signatures of circumbinary discs – I. Kinematics

We present five morphological and kinematic criteria to aid in asserting the binary nature of a protoplanetary disc, based on 3D hydrodynamical simulations of circumbinary discs post-processed with Monte Carlo radiative transfer. We find that circumbinary discs may be identified by (i) a central cavity, (ii) spiral arms both in and outside of their central cavities, (iii) non-localized perturbations in their iso-velocity curves, (iv) asymmetry between the lines of maximum speed of the blue shifted and redshifted wings, and (v) asymmetry between the area of the blue shifted and redshifted wings. In this study, we provide quantitative metrics for the last two criteria that can be used, in conjunction with the morphological criteria, to signal whether a protoplanetary disc is likely to be a circumbinary disc.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Deformations, relaxation, and broken symmetries in liquids, solids, and glasses: A unified topological field theory

Here, we combine hydrodynamic and field theoretic methods to develop a general theory of phonons as Goldstone bosons in crystals, glasses, and liquids based on nonaffine displacements and the consequent Goldstone phase relaxation. We relate the conservation, or lack thereof, of specific higher-form currents with properties of the underlying deformation field—nonaffinity—which dictates how molecules move under an applied stress or deformation. In particular, the single-valuedness of the deformation field is associated with conservation of higher-form charges that count the number of topological defects. Our formalism predicts, from first principles, the presence of propagating shear waves above a critical wave vector in liquids, thus giving a formal derivation of the phenomenon in terms of fundamental symmetries. The same picture provides also a theoretical explanation of the corresponding “positive sound dispersion” phenomenon for longitudinal sound. Importantly, accordingly to our theory, the main collective relaxation timescale of a liquid or a glass (known as the α relaxation for the latter) is given by the phase relaxation time, which is not necessarily related to the Maxwell time. Finally, we build a nonequilibrium effective action using the in-in formalism defined on the Schwinger-Keldysh contour, that further supports the emerging picture. In summary, our work suggests that the fundamental difference between solids, fluids, and glasses has to be identified with the associated generalized higher-form global symmetries and their topological structure, and that the Burgers vector for the displacement fields serves as a suitable topological order parameter distinguishing the solid (ordered) phase and the amorphous ones (fluids, glasses).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The coexistence of the streaming instability and the vertical shear instability in protoplanetary disks

The streaming instability is a promising mechanism to induce the formation of planetesimals. Nonetheless, this process has been found in previous studies to require either a dust-to-gas surface density ratio or a dust size that is enhanced compared to observed values. Employing two-dimensional global simulations of protoplanetary disks, we show that the vertical shear instability and the streaming instability in concert can cause dust concentration that is sufficient for planetesimal formation for lower surface density ratios and smaller dust sizes than the streaming instability in isolation, and in particular under conditions that are consistent with observational constraints. This is because dust overdensities forming in pressure bumps induced by the vertical shear instability act as seeds for the streaming instability and are enhanced by it. While our two-dimensional model does not include self-gravity, we find that strong dust clumping and the formation (and dissolution) of gravitationally unstable overdensities can be robustly inferred from the evolution of the maximum or the mean dust-to-gas volume density ratio. The vertical shear instability puffs up the dust layer to an average mid-plane dust-to-gas density ratio that is significantly below unity. We therefore find that reaching a mid-plane density ratio of one is not necessary to trigger planetesimal formation via the streaming instability when it acts in unison with the vertical shear instability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Structure of a chaotic tangle of quantum vortices in turbulent superfluid liquids and in a Bose – Einstein condensate

Based on the theory of the thermodynamic equilibrium in a system of quantum vortices in superfluid liquids and in a Bose – Einstein condensate in the presence of a counterflow of normal and superfluid components, we investigate the structure of a chaotic tangle of quantum vortices in turbulent superfluid liquids. Using the characteristic functional method, the properties of hydrodynamic vortex filaments are examined. It is shown that the average curvature of the vortex lines is on the order of the intervortex distance, with the proportionality coefficient being independent of the counterflow velocity. It is found that the degree of anisotropy of the vortex loops does not depend on the applied counterflow velocity. The obtained results explain the origin of the anisotropy as well as the relationship between the curvature of the lines and the intervortex space and their dependence on the parameters of the problem. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

The impact of the $\text{WHIM}$ on the $\text{IGM}$ thermal state determined from the low- z Lyman $\alpha$ forest

At z ≲ 1, shock heating caused by large-scale velocity flows and possibly violent feedback from galaxy formation, converts a significant fraction of the cool gas (T ~ 10 4 K) in the intergalactic medium (IGM) into warm–hot phase (WHIM) with T > 10 5 K, resulting in a significant deviation from the previously tight power-law IGM temperature–density relationship, T = T 0 (P/$\bar{p}$) γ-1 ⁠. This study explores the impact of the WHIM on measurements of the low-z IGM thermal state, [T 0 , γ], based on the b–NH 1 distribution of the Ly α forest. Exploiting a machine learning-enabled simulation-based inference method trained on Nyx hydrodynamical simulations, we demonstrate that [T 0 , γ] can still be reliably measured from the b–NH 1 distribution at z = 0.1, notwithstanding the substantial WHIM in the IGM. To investigate the effects of different feedback, we apply this inference methodology to mock spectra derived from the IllustrisTNG and Illustris simulations at z = 0.1. The results suggest that the underlying [T 0 , γ] of both simulations can be recovered with biases as low as |Δlog(T 0 /K)| ≲ 0.05 dex, |Δγ| ≲ 0.1, smaller than the precision of a typical measurement. Given the large differences in the volume-weighted WHIM fractions between the three simulations (Illustris 38 percent, IllustrisTNG 10 percent, and Nyx 4 per cent), we conclude that the b–N H1 distribution is not sensitive to the WHIM under realistic conditions. Finally, we investigate the physical properties of the detectable Ly α absorbers, and discover that although their T and Δ distributions remain mostly unaffected by feedback, they are correlated with the photoionization rate used in the simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Error Study of FlecSim Using the Method of Manufactured Solutions [Slides]

Error was calculated using exact solution chosen with the method of manufactured solutions. Hydrodynamics match expected order of accuracy for spatial error within 1.5% for all variables tested. Radiation transport matches expected order of temporal accuracy within 0.5% for all variables tested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Ring Formation in Protoplanetary Disks Driven by an Eccentric Instability

We find that, under certain conditions, protoplanetary disks may spontaneously generate multiple, concentric gas rings without an embedded planet through an eccentric cooling instability. Using both linear theory and nonlinear hydrodynamics simulations, we show that a variety of background states may trap a slowly processing, one-armed spiral mode that becomes unstable when a gravitationally stable disk rapidly cools. The angular momentum required to excite this spiral comes at the expense of nonuniform mass transport that generically results in multiple rings. For example, one long-term hydrodynamics simulation exhibits four long-lived, axisymmetric gas rings. We verify the instability evolution and ring-formation mechanism from first principles with our linear theory, which shows remarkable agreement with the simulation results. Dust trapped in these rings may produce observable features consistent with observed disks. Additionally, direct detection of the eccentric gas motions may be possible when the instability saturates, and any residual eccentricity left over in the rings at later times may also provide direct observational evidence of this mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The AGORA High-resolution Galaxy Simulations Comparison Project. V. Satellite Galaxy Populations in a Cosmological Zoom-in Simulation of a Milky Way–Mass Halo

We analyze and compare the satellite halo populations at z ~ 2 in the high-resolution cosmological zoom-in simulations of a 10 12 M ⊙ target halo (z = 0 mass) carried out on eight widely used astrophysical simulation codes (Art-I, Enzo, Ramses, Changa, Gadget-3, Gear, Arepo-t, and Gizmo) for the AGORA High-resolution Galaxy Simulations Comparison Project. We use slightly different redshift epochs near z = 2 for each code (hereafter "z ~ 2") at which the eight simulations are in the same stage in the target halo's merger history. After identifying the matched pairs of halos between the CosmoRun simulations and the DMO simulations, we discover that each CosmoRun halo tends to be less massive than its DMO counterpart. When we consider only the halos containing stellar particles at z ~ 2, the number of satellite galaxies is significantly fewer than that of dark matter halos in all participating AGORA simulations and is comparable to the number of present-day satellites near the Milky Way or M31. The so-called "missing satellite problem" is fully resolved across all participating codes simply by implementing the common baryonic physics adopted in AGORA and the stellar feedback prescription commonly used in each code, with sufficient numerical resolution (≲100 proper pc at z = 2). We also compare other properties such as the stellar mass–halo mass relation and the mass–metallicity relation. Our work highlights the value of comparison studies such as AGORA, where outstanding problems in galaxy formation theory are studied simultaneously on multiple numerical platforms.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1982

X-rays from quasars and active galaxies

Some of the most important theoretical ideas and observations for quasars and the nuclei of active galaxies are reviewed, and areas of future research are suggested. Emphasis is on the nature of the power source, the radiation processes, and the mechanism for formation and collimation of jets. Phenomena that produce X-rays are of particular concern. Particular topics discussed are the observed and expected time variabilility, the gas supply mechanisms and luminosity evolution, thermal and nonthermal radiation processes, observed and theoretical spectra, criteria for thermalization of electrons and ions, effects of electron-positron pairs on relativistic plasmas, hydrodynamic, electrodynamic and inertial methods for producing and confining jets. The paper concludes with a list of needed observations.

Lightman, A. P.↗
NASA NTRS1993

Galaxy tracers in N-body simulations

Using the method of smoothed particle hydrodynamics, we have modeled the formation of a compact group of galaxies with sufficient resolution to trace galaxies. Radiative cooling allows the baryons to dissipate their thermal energy and collapse to overdensities characteristic of real galaxies. With their cross section greatly reduced, these galaxy tracers remain distinct during cluster formation while their dark matter halos merge. In addition, the number density, the mass distribution function, and even the morphology of these objects are similar to those of observed galaxies. A viable population of galaxy tracers can be unambiguously defined.

Summers, F. J.↗
NASA NTRS2015

Multiscale and Multiphysics Modeling of Additive Manufacturing of Advanced Materials

The objective of this proposed project is to research and develop a prediction tool for advanced additive manufacturing (AAM) processes for advanced materials and develop experimental methods to provide fundamental properties and establish validation data. Aircraft structures and engines demand materials that are stronger, useable at much higher temperatures, provide less acoustic transmission, and enable more aeroelastic tailoring than those currently used. Significant improvements in properties can only be achieved by processing the materials under nonequilibrium conditions, such as AAM processes. AAM processes encompass a class of processes that use a focused heat source to create a melt pool on a substrate. Examples include Electron Beam Freeform Fabrication and Direct Metal Deposition. These types of additive processes enable fabrication of parts directly from CAD drawings. To achieve the desired material properties and geometries of the final structure, assessing the impact of process parameters and predicting optimized conditions with numerical modeling as an effective prediction tool is necessary. The targets for the processing are multiple and at different spatial scales, and the physical phenomena associated occur in multiphysics and multiscale. In this project, the research work has been developed to model AAM processes in a multiscale and multiphysics approach. A macroscale model was developed to investigate the residual stresses and distortion in AAM processes. A sequentially coupled, thermomechanical, finite element model was developed and validated experimentally. The results showed the temperature distribution, residual stress, and deformation within the formed deposits and substrates. A mesoscale model was developed to include heat transfer, phase change with mushy zone, incompressible free surface flow, solute redistribution, and surface tension. Because of excessive computing time needed, a parallel computing approach was also tested. In addition, after investigating various methods, a Smoothed Particle Hydrodynamics Model (SPH Model) was developed to model wire feeding process. Its computational efficiency and simple architecture makes it more robust and flexible than other models. More research on material properties may be needed to realistically model the AAM processes. A microscale model was developed to investigate heterogeneous nucleation, dendritic grain growth, epitaxial growth of columnar grains, columnar-to-equiaxed transition, grain transport in melt, and other properties. The orientations of the columnar grains were almost perpendicular to the laser motion's direction. Compared to the similar studies in the literature, the multiple grain morphology modeling result is in the same order of magnitude as optical morphologies in the experiment. Experimental work was conducted to validate different models. An infrared camera was incorporated as a process monitoring and validating tool to identify the solidus and mushy zones during deposition. The images were successfully processed to identify these regions. This research project has investigated multiscale and multiphysics of the complex AAM processes thus leading to advanced understanding of these processes. The project has also developed several modeling tools and experimental validation tools that will be very critical in the future of AAM process qualification and certification.

Liou, Frank↗