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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Simulating Lindbladian evolution with non-Abelian symmetries: Ballistic front propagation in the SU(2) Hubbard model with a localized loss

In this work, we develop a non-Abelian time evolving block decimation (NA-TEBD) approach to study open systems governed by Lindbladian time evolution, while exploiting an arbitrary number of Abelian or non-Abelian symmetries. We illustrate this method in a one-dimensional fermionic SU(2) Hubbard model on a semi-infinite lattice with localized particle loss at one end. We observe a ballistic front propagation with strongly renormalized front velocity, and a hydrodynamic current density profile. For large loss rates, a suppression of the particle current is observed, as a result of the quantum Zeno effect. Operator entanglement is found to propagate faster than the depletion profile, preceding the latter.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Fluctuating hydrodynamics of multi-species reactive mixtures

We formulate and study computationally the fluctuating compressible Navier-Stokes equations for reactive multi-species fluid mixtures. We contrast two different expressions for the covariance of the stochastic chemical production rate in the Langevin formulation of stochastic chemistry, and compare both of them to predictions of the chemical master equation for homogeneous well-mixed systems close to and far from thermodynamic equilibrium. We develop a numerical scheme for inhomogeneous reactive flows, based on our previous methods for non-reactive mixtures [Balakrishnan , Phys. Rev. E 89, 013017 (2014)]. We study the suppression of non-equilibrium long-ranged correlations of concentration fluctuations by chemical reactions, as well as the enhancement of pattern formation by spontaneous fluctuations. Good agreement with available theory demonstrates that the formulation is robust and a useful tool in the study of fluctuations in reactive multi-species fluids. At the same time, several problems with Langevin formulations of stochastic chemistry are identified, suggesting that future work should examine combining Langevin and master equation descriptions of hydrodynamic and chemical fluctuations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Tempus v. 1.0

Tempus is a software framework for the time integration of ordinary differential equations (ODEs), partial differential equations (PDEs), and differential-algebraic equations (DAEs). Examples of time-integration methods available are Forward Euler, Backward Euler, Explicit Runge-Kutta and Implicit Runge-Kutta methods. Tempus can be used from small systems of equations (e.g., single ODEs for the time evolution of plasticity models) to large-scale transient simulations requiring exascale computing (e.g., flow fields around reentry vehicles and magneto-hydrodynamics). Tempus will be released as part of Trilinos under a BSD open source license.

Obert, CurtisC.↗
DOE OSTI2023

REYNO: A Reactive Hydrodynamics Modeling Suite

REYNO is a Python library that is primarily designed as a platform for the agile development and testing of novel reactive flow models. Equations of state and multi-step rate laws can be implemented using the provided class hierarchy with different thermodynamic closure conditions. The library includes a one-dimensional Lagrangian hydrodynamic solver to simulate multi-material reactive or inert problems in planar, cylindrical, or spherical coordinates. Custom time-dependent boundary conditions are also available for problems such as ramp compression. This document describes the different numerical methods and algorithms that are included in REYNO.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1990

Modeling the dust size distribution in comets with dust fragmentation

A hydrodynamic model was developed of a spherically symmetric dusty gas flow in a cometary atmosphere assuming a single fluid, inviscid, perfect gas. The hydrodynamics for gas and dust, which involves the gas drag force (momentum transfer), heat exchange between gas and dust, photodissociation energy for H2O gas, and radiative heating and cooling terms for dust particles, are solved using the Gear method for stiff, coupled differential equations. Calculations were done with a dust size distribution for radii alpha = 0.01 micron to 10 cm with densities variable with the size. A nucleus size of 4.0 km radius with a density of 0.5 g/cu cm and a total dust-to-gas mass ratio chi = 1 was adopted. There are indications from in situ observations that dust particle fragment into smaller ones. Fragmentation of dust particles was incorporated into the model. This is done by adding source and sink terms in the continuity equations for the dust. Lifetimes for the decay of dust particles were assumed as a function of particle size. It is also assumed that dust particles always fragment only into the next smaller size.

Konno, Ichishiro↗
DOE OSTI2014

Heating the intergalactic medium by X-rays from population III binaries in high-redshift galaxies

Due to their long mean free path, X-rays are expected to have an important impact on cosmic reionization by heating and ionizing the intergalactic medium (IGM) on large scales, especially after simulations have suggested that Population III (Pop III) stars may form in pairs at redshifts as high as 20-30. We use the Pop III distribution and evolution from a self-consistent cosmological radiation hydrodynamic simulation of the formation of the first galaxies and a simple Pop III X-ray binary model to estimate their X-ray output in a high-density region larger than 100 comoving (Mpc){sup 3}. We then combine three different methods—ray tracing, a one-zone model, and X-ray background modeling—to investigate the X-ray propagation, intensity distribution, and long-term effects on the IGM thermal and ionization state. The efficiency and morphology of photoheating and photoionization are dependent on the photon energies. The sub-kiloelectronvolt X-rays only impact the IGM near the sources, while the kiloelectronvolt photons contribute significantly to the X-ray background and heat and ionize the IGM smoothly. The X-rays just below 1 keV are most efficient in heating and ionizing the IGM. We find that the IGM might be heated to over 100 K by z = 10 and the high-density source region might reach 10{sup 4} K, limited by atomic hydrogen cooling. This may be important for predicting the 21 cm neutral hydrogen signals. On the other hand, the free electrons from X-ray ionizations are not enough to contribute significantly to the optical depth of the cosmic microwave background to the Thomson scattering.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Data-Driven Discovery of Active Nematic Hydrodynamics

Active nematics can be modeled using phenomenological continuum theories that account for the dynamics of the nematic director and fluid velocity through partial differential equations (PDEs). While these models provide a statistical description of the experiments, the relevant terms in the PDEs and their parameters are usually identified indirectly. We adapt a recently developed method to automatically identify optimal continuum models for active nematics directly from spatio-temporal data, via sparse regression of the coarse-grained fields onto generic low order PDEs. After extensive benchmarking, we apply the method to experiments with microtubule-based active nematics, finding a surprisingly minimal description of the system. Furthermore, our approach can be generalized to gain insights into active gels, microswimmers, and diverse other experimental active matter systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Feasibility of flow studies at NICA/MPD

In the light of recent developments in heavy ion physic, anisotropic flow measurements play a key role in a better understanding of the hot and dense barionic matter. In the presented article a short introduction to the proposed NICA/MPD project is given, as well as a brief description of the event-plane method used to estimate the elliptic flow of reconstructed and identified hadrons (p, π, Λ)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Multi-frame, ultrafast, x-ray microscope for imaging shockwave dynamics

Inertial confinement fusion (ICF) holds increasing promise as a potential source of abundant, clean energy, but has been impeded by defects such as micro-voids in the ablator layer of the fuel capsules. It is critical to understand how these micro-voids interact with the laser-driven shock waves that compress the fuel pellet. At the Matter in Extreme Conditions (MEC) instrument at the Linac Coherent Light Source (LCLS), we utilized an x-ray pulse train with ns separation, an x-ray microscope, and an ultrafast x-ray imaging (UXI) detector to image shock wave interactions with micro-voids. To minimize the high- and low-frequency variations of the captured images, we incorporated principal component analysis (PCA) and image alignment for flat-field correction. After applying these techniques we generated phase and attenuation maps from a 2D hydrodynamic radiation code (xRAGE), which were used to simulate XPCI images that we qualitatively compare with experimental images, providing a one-to-one comparison for benchmarking material performance. Moreover, we implement a transport-of-intensity (TIE) based method to obtain the average projected mass density (areal density) of our experimental images, yielding insight into how defect-bearing ablator materials alter microstructural feature evolution, material compression, and shock wave propagation on ICF-relevant time scales.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Tracer Particles for Core-collapse Supernova Nucleosynthesis: The Advantages of Moving Backward

Abstract After decades, the theoretical study of core-collapse supernova explosions is moving from parameterized, spherically symmetric models to increasingly realistic multidimensional simulations. However, obtaining nucleosynthesis yields based on such multidimensional core-collapse supernova simulations is not straightforward. Frequently, tracer particles are employed. Tracer particles may be tracked in situ during the simulation, but often they are reconstructed in a post-processing step based on the information saved during the hydrodynamic simulation. Reconstruction can be done in a number of ways, and here we compare the approaches of backward and forward integration of the equations of motion to the results based on inline particle trajectories. We find that both methods agree reasonably well with the inline results for isotopes for which a large number of particles contribute. However, for rarer isotopes that are produced only by a small number of particle trajectories, deviations can be large. For our setup, we find that backward integration leads to better agreement with the inline particles by more accurately reproducing the conditions following freeze-out from nuclear statistical equilibrium, because the establishment of nuclear statistical equilibrium erases the need for detailed trajectories at earlier times. Based on our results, if inline tracers are unavailable, we recommend backward reconstruction to the point when nuclear statistical equilibrium was last applied, with an interval between simulation snapshots of at most 1 ms for nucleosynthesis post-processing.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

ON THE MISALIGNMENT BETWEEN CHROMOSPHERIC FEATURES AND THE MAGNETIC FIELD ON THE SUN

Observations of the upper chromosphere show an enormous amount of intricate fine structure. Much of this comes in the form of linear features, which are most often assumed to be well aligned with the direction of the magnetic field in the low plasma β regime that is thought to dominate the upper chromosphere. We use advanced radiative magnetohydrodynamic simulations, including the effects of ion-neutral interactions (using the generalized Ohm’s law) in the partially ionized chromosphere, to show that the magnetic field is often not well aligned with chromospheric features. This occurs where the ambipolar diffusion is large, i.e., ions and neutral populations decouple as the ion-neutral collision frequency drops, allowing the field to slip through the neutral population; where currents perpendicular to the field are strong; and where thermodynamic timescales are longer than or similar to those of ambipolar diffusion. We find this often happens in dynamic spicule or fibril-like features at the top of the chromosphere. This has important consequences for field extrapolation methods, which increasingly use such upper chromospheric features to help constrain the chromospheric magnetic field: our results invalidate the underlying assumption that these features are aligned with the field. In addition, our results cast doubt on results from 1D hydrodynamic models, which assume that plasma remains on the same field lines. Finally, our simulations show that ambipolar diffusion significantly alters the amount of free energy available in the coronal part of our simulated volume, which is likely to have consequences for studies of flare initiation.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2022

The extended atmosphere and circumstellar environment of the cool evolved star VX Sagittarii as seen by MATISSE

Context. VX Sgr is a cool, evolved, and luminous red star whose stellar parameters are difficult to determine, which affects its classification. Aims. We aim to spatially resolve the photospheric extent as well as the circumstellar environment. Methods. We used interferometric observations obtained with the MATISSE instrument in the L (3–4 μm), M (4.5–5 μm), and N (8–13 μm) bands. We reconstructed monochromatic images using the MIRA software. We used 3D radiation-hydrodynamics simulations carried out with CO5BOLD and a uniform disc model to estimate the apparent diameter and interpret the stellar surface structures. Moreover, we employed the radiative transfer codes OPTIM3D and RADMC3D to compute the spectral energy distribution for the L, M, and N bands, respectively. Results. MATISSE observations unveil, for the first time, the morphology of VX Sgr across the L, M, and N bands. The reconstructed images show a complex morphology with brighter areas whose characteristics depend on the wavelength probed. We measured the angular diameter as a function of the wavelength and showed that the photospheric extent in the L and M bands depends on the opacity through the atmosphere. In addition to this, we also concluded that the observed photospheric inhomogeneities can be interpreted as convection-related surface structures. The comparison in the N band yielded a qualitative agreement between the N-band spectrum and simple dust radiative transfer simulations. However, it is not possible to firmly conclude on the interpretation of the current data because of the difficulty in constraing the model parameters using the limited accuracy of our absolute flux calibration. Conclusions. MATISSE observations and the derived reconstructed images unveil the appearance of VX Sgr’s stellar surface and circumstellar environment across a very large spectral domain for the first time.

stars: atmospheres↗
DOE OSTI2015

Physical Processes in the MAGO/MFT Systems

The Monograph is devoted to theoretical discussion of the physical effects, which are most significant for the alternative approach to the problem of controlled thermonuclear fusion (CTF): the MAGO/MTF approach. The book includes the description of the approach, its difference from the major CTF systems—magnetic confinement and inertial confinement systems. General physical methods of the processes simulation in this approach are considered, including plasma transport phenomena and radiation, and the theory of transverse collisionless shock waves, the surface discharges theory, important for such kind of research. Different flows and magneto-hydrodynamic plasma instabilities occurring in the frames of this approach are also considered. In virtue of the general physical essence of the considered phenomena the presented results are applicable to a wide range of plasma physics and hydrodynamics processes. The book is intended for the plasma physics and hydrodynamics specialists, post-graduate students, and senior students-physicists.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Galaxy cluster matter profiles - I. Self-similarity, mass calibration, and observable-mass relation validation employing cluster mass posteriors

We present a study of the weak lensing inferred matter profiles ΔΣ(R) of 698 South Pole Telescope (SPT) thermal Sunyaev-Zel’dovich effect (tSZE) selected and MCMF optically confirmed galaxy clusters in the redshift range 0.25 < z < 0.94 that have associated weak gravitational lensing shear profiles from the Dark Energy Survey (DES). Rescaling these profiles to account for the mass dependent size and the redshift dependent density produces average rescaled matter profiles ΔΣ(R/R200c)/(ρcritR200c) with a lower dispersion than the unscaled ΔΣ(R) versions, indicating a significant degree of self-similarity. Galaxy clusters from hydrodynamical simulations also exhibit matter profiles that suggest a high degree of self-similarity, with RMS variation among the average rescaled matter profiles with redshift and mass falling by a factor of approximately six and 23, respectively, compared to the unscaled average matter profiles. We employed this regularity in a new Bayesian method for weak lensing mass calibration that employs the so-called cluster mass posterior P(M200|ζ̂, λ̂, z), which describes the individual cluster masses given their tSZE (ζ̂) and optical (λ̂, z) observables. This method enables simultaneous constraints on richness λ-mass and tSZE detection significance ζ-mass relations using average rescaled cluster matter profiles. We validated the method using realistic mock datasets and present observable-mass relation constraints for the SPT×DES sample, where we constrained the amplitude, mass trend, redshift trend, and intrinsic scatter. Our observable-mass relation results are in agreement with the mass calibration derived from the recent cosmological analysis of the SPT×DES data based on a cluster-by-cluster lensing calibration. Our new mass calibration technique offers a higher efficiency when compared to the single cluster calibration technique. We present new validation tests of the observable-mass relation that indicate the underlying power-law form and scatter are adequate to describe the real cluster sample but that also suggest a redshift variation in the intrinsic scatter of the λ-mass relation may offer a better description. In addition, the average rescaled matter profiles offer high signal-to-noise ratio (S/N) constraints on the shape of real cluster matter profiles, which are in good agreement with available hydrodynamical ΛCDM simulations. This high S/N profile contains information about baryon feedback, the collisional nature of dark matter, and potential deviations from general relativity.Key words: gravitational lensing: weak / galaxies: clusters: general / large-scale structure of Universe

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Direct photon elliptic flow in Pb–Pb collisions at s NN = 2.76 TeV

The elliptic flow of inclusive and direct photons was measured at mid-rapidity in two centrality classes 0–20% and 20–40% in Pb–Pb collisions at s NN =2.76 TeV by ALICE. Photons were detected with the highly segmented electromagnetic calorimeter PHOS and via conversions in the detector material with the e + e - pairs reconstructed in the central tracking system. The results of the two methods were combined and the direct-photon elliptic flow was extracted in the transverse momentum range 0.9 < p T < 6.2 GeV/ c . A comparison to RHIC data shows a similar magnitude of the measured direct-photon elliptic flow. Hydrodynamic and transport model calculations are systematically lower than the data, but are found to be compatible.

Acharya, S.↗
DOE OSTI2017

Fast heating of fuel assembled in a spherical deuterated polystyrene shell target by counter-irradiating tailored laser pulses delivered by a HAMA 1 Hz ICF driver

Fast heating is a method of heating an assembled high-density plasma into a hot state by irradiating it with short-duration (sub-picosecond), high-intensity ($>$ $10^{18}$ W ${\rm cm}^{-2}$ ) laser pulses before the plasma expands and dissolves hydrodynamically. In this paper, we present detailed experimental results of fast heating fuel assembled in a spherical deuterated polystyrene shell target of 500 ${\mu}$ m diameter and 7 μm thickness with counterbeam illumination by using a HAMA 1 Hz, 5.9 J inertial confinement fusion laser driver with pulse tailoring. These tailored pulses contain three pulses in sequence: a 'foot' pulse of 2.4 J/25 ns, a 'spike' pulse of 0.5 J/300 ps and a 'heater' pulse of 0.4 J/110 fs; these pulses are designed to assemble the fuel and heat it. By varying the energy of the foot pulse, we find that fast heating the fuel is achieved only if the fuel is weakly ablated by the foot pulse and then shock-assembled by the spike pulse into the target centre so that the heater pulse can access the fuel with a focal intensity greater than $10^{18}$ W ${\rm cm}^{-2}$ . Without a foot pulse, the heater pulse contributes to assembling the fuel. For higher foot-pulse energies, the heater pulse drives a hydrodynamic motion with speeds of the order 107 cm ${\rm s}^{-1}$ with intensities of the order $10^{17}$ W ${\rm cm}^{-2}$ , resulting in re-assembling and additional heating of the pre-assembled fuel. Once a shock-assembled core is achieved at the target centre, we succeed qualitatively in fast heating the core for shots in sequence with variations of laser energy within 18%. The coupling efficiency from the heating laser to the core is inferred to be $(10 \pm 2)$ % in total: $(8 \pm 1.6)$ % for the ionized bulk electrons and $(2 \pm 0.4)$ % for the bulk ions. The fusion neutron spectrum detected on the laser axis exhibits peaks at 1.0 MeV, 1.7 MeV and 3.8 MeV. These peaks are attributed to the $C(d, n){\hspace{0pt}}^{13}$ N and $d(d, n){\hspace{0pt}}^3$ He reactions induced by counterpropagating fast deuterons accelerated by the photon pressure of the heating pulses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Virtual Flow Solver - Geophysics: A 3D Incompressible Navier-Stokes Solver

Virtual Flow Solver - Geophysics (VFS-Geophysics) is a three-dimensional (3D) incompressible Navier-Stokes solver based on the Curvilinear Immersed Boundary (CURVIB) method. The CURVIB is a sharp interface type of immersed boundary (IB) method that enables the simulation of fluid flows in the presence of geometrically complex moving bodies. The CURVIB method can be applied to wind/MHK turbine simulations and energy applications. VFS-Geophysics is the result of many years of research work by several graduate students, post-docs, and research associates that have been involved in the Computational Hydrodynamics and Biofluids Laboratory directed by Professor Fotis Sotiropoulos. The preparation of the present manual has been supported by the U.S. Department of Energy (DE-EE 0005482).

16 TIDAL AND WAVE POWER↗
DOE OSTI2013

Examination of Beam Profile Effects on DARHT Axis II Forward Modeling

The reconstruction of the density field for hydrodynamic tests at DARHT is of paramount importance to the LANL hydrodynamic program. Arguably, the most accurate means by which the density field may be obtained is via the development of a forward model. The forward model incorporates all of the salient elements of the radiographic process and consequently is the method by which density reconstructions have been performed by LANL for over a decade utilizing the BIE.

42 ENGINEERING↗