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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The aluminium-26 distribution in a cosmological simulation of a Milky Way-type Galaxy

Context . The 1.8 MeV γ -rays corresponding to the decay of the radioactive isotope 26 Al (with a half-life of 0.72 Myr ) have been observed by the SPI detector on the INTEGRAL spacecraft and extensively used as a tracer of star formation and current nucleosynthetic activity in the Milky Way Galaxy. Further information is encoded in the observation related to the higher 26 Al content found in regions of the Galaxy with the highest line-of-sight (LoS) velocity relative to an observer located in the Solar System. However, this feature remains unexplained. Aims . We ran a cosmological “zoom-in” chemodynamical simulation of a Milky Way-type galaxy, including the production and decays of radioactive nuclei in a fully self-consistent way. We then analyzed the results to follow the evolution of 26 Al throughout the lifetime of the simulated galaxy to provide a new method for interpreting the 26 Al observations. Methods . We included the massive star sources of 26 Al in the Galaxy and its radioactive decay into a state-of-the-art galactic chemical evolution model, coupled with cosmological growth and hydrodynamics. This approach allowed us to follow the spatial and temporal evolution of the 26 Al content in the simulated galaxy. Results . Our results are in agreement with the observations with respect to the fact that gas particles in the simulation with relatively higher 26 Al content also have the highest LoS velocities. On the other hand, gas particles with relatively lower 26 Al content (i.e., not bright enough to be observed) generally display the lowest LoS velocities. However, this result is not conclusive because the overall rotational velocity of our simulated galaxy is higher than that observed for cold CO gas in the Milky Way Galaxy. Furthermore, we found no significant correlation between gas temperature, rotational velocity, and 26 Al content at any given radius. We also found the presence of transient 26 Al-rich spots at low LoS velocities and we show that one such spot had been captured by the INTEGRAL/SPI data. Based on our model, we present a prediction for the detection of 1.8 MeV γ -rays by the future COSI mission. We find that according to our model, the new instrument will be able to observe similar 26 Al-emission patterns to those seen by INTEGRAL/SPI.

Wehmeyer, B. (ORCID:0009000254149292)↗
DOE OSTI2017

Analysis of the Conformally Flat Approximation for Binary Neutron Star Initial Conditions

The spatially conformally flat approximation (CFA) is a viable method to deduce initial conditions for the subsequent evolution of binary neutron stars employing the full Einstein equations. Here in this paper, we analyze the viability of the CFA for the general relativistic hydrodynamic initial conditions of binary neutron stars. We illustrate the stability of the conformally flat condition on the hydrodynamics by numerically evolving ~100 quasicircular orbits. We illustrate the use of this approximation for orbiting neutron stars in the quasicircular orbit approximation to demonstrate the equation of state dependence of these initial conditions and how they might affect the emergent gravitational wave frequency as the stars approach the innermost stable circular orbit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Understanding the compound flood risk along the coast of the contiguous United States

Abstract. Compound flooding is a type of flood event caused by multiple flood drivers. The associated risk has usually been assessed using statistics-based analyses or hydrodynamics-based numerical models. This study proposes a compound flood (CF) risk assessment (CFRA) framework for coastal regions in the contiguous United States (CONUS). In this framework, a large-scale river model is coupled with a global ocean reanalysis dataset to (a) evaluate the CF exposure related to the coastal backwater effects on river basins, and (b) generate spatially distributed data for analyzing the CF hazard using a bivariate statistical model of river discharge and storm surge. The two kinds of risk are also combined to achieve a holistic understanding of the continental-scale CF risk. The estimated CF risk shows remarkable inter- and intra-basin variabilities along the CONUS coast with more variabilities in the CF hazard over the US west and Gulf coastal basins. Different risk assessment methods present significantly different patterns in a few key regions such as the San Francisco Bay area, the lower Mississippi River, and Puget Sound. Our results highlight the need to weigh different CF risk measures and avoid using single statistics-based or hydrodynamics-based CFRAs. Uncertainty sources in these CFRAs include the use of gauge observations, which cannot account for the flow physics or resolve the spatial variability of risks, and underestimations of the flood extremes and the dependence of CF drivers in large-scale models, highlighting the importance of understanding the CF risks for developing a more robust CFRA.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1990

Thermal instabilities in protogalactic clouds

The means by which a protogalaxy can fragment to form the first generation of stars and globular clusters remains an important problem in astrophysics. Gravitational instabilities grow on timescales too long to drive fragmentation before the background density grows by many orders of magnitude (see Murray and Lin 1989a, and references therein). Thermal instability provides a much more likely mechanism. After its initial collapse, a protogalactic cloud is expected to be shock heated to its virial temperature approx. 10(exp 6) K. Cooling by H and He+ below 10(exp 6) K has a negative slope, so that the cloud is subject to strong thermal instabilities. Density enhancements may then grow rapidly, fragmenting the protogalaxy as it cools to lower temperatures. The role of dynamical effects upon the growth of perturbations is considered here. The method used is similar to that used in Murray and Lin (1989a; see also the Erratum to appear September 15), which examined the growth of thermal instabilities with a one-dimensional Lagrangian hydrodynamics code, written for spherical symmetry. Perturbed regions therefore take the form of shells. The dynamical variables are integrated explicitly, while the temperature, ionization fraction, and molecular fraction are integrated implicitly, and account is taken for non-equilibrium values of these quantities.

Murray, S. D.↗
NASA NTRS2001

Computational Relativistic Astrophysics Using the Flowfield-Dependent Variation Theory

Theoretical models, observations and measurements have preoccupied astrophysicists for many centuries. Only in recent years, has the theory of relativity as applied to astrophysical flows met the challenges of how the governing equations can be solved numerically with accuracy and efficiency. Even without the effects of relativity, the physics of magnetohydrodynamic flow instability, turbulence, radiation, and enhanced transport in accretion disks has not been completely resolved. Relativistic effects become pronounced in such cases as jet formation from black hole magnetized accretion disks and also in the study of Gamma-Ray bursts (GRB). Thus, our concern in this paper is to reexamine existing numerical simulation tools as to the accuracy and efficiency of computations and introduce a new approach known as the flowfield-dependent variation (FDV) method. The main feature of the FDV method consists of accommodating discontinuities of shock waves and high gradients of flow variables such as occur in turbulence and unstable motions. In this paper, the physics involved in the solution of relativistic hydrodynamics and solution strategies of the FDV theory are elaborated. The general relativistic astrophysical flow and shock solver (GRAFSS) is introduced, and some simple example problems for Computational Relativistic Astrophysics (CRA) are demonstrated.

Richardson, G. A.↗
DOE OSTI2026

A cell-centered AMR-ALE framework for 3D multi-material hydrodynamics. Part I: Lagrangian and indirect Euler AMR algorithms

Many applications of physics and engineering involve wide ranges of time and spatial scales. The numerical simulation of localized small scales such as shock waves and material interfaces requires a large number of computational cells in these regions. For these applications, Lagrangian and Arbitrary-Lagrangian-Eulerian (ALE) related methods are engaging since the moving mesh feature naturally brings mesh cells on shock discontinuities and material interfaces are carefully captured. In addition, Adaptive-Mesh-Refinement (AMR) strategies aim to optimize computational resources by concentrating finer mesh cells only in areas of interest while using coarser cells elsewhere. A key but challenging AMR requirement consists in efficiently distributing the computational effort to achieve high accuracy without the prohibitive computational costs associated with uniformly fine grids. Here, in this document, the coupling of the p4est AMR library with a cell-centered Lagrangian scheme is presented with the goal to perform reliable 3D Lagrangian-AMR and indirect Euler-AMR multi-material simulations. In particular, it is shown that starting from a 3D indirect ALE code, the memory management and load balancing requirements can be delegated to an external library (here the p4est library) to unlock ALE-AMR capabilities. First, we present a strategy to transcribe the octant-based connectivity of the 3D AMR framework with that of an unstructured mesh of polygonal cells used in Lagrangian hydrodynamics. Then, we show how refinement and coarsening operations must be adapted to the particular Lagrangian framework to ensure the conservation of volume during those steps. Finally, several numerical test cases are presented that demonstrate the capabilities of the Lagrangian-AMR and indirect Euler-AMR algorithms.

3D cell-centered Lagrangian numerical scheme↗
NASA NTRS1994

Accretion disk coronae in high-luminosity systems

We present the results of self-consistent models of Compton-heated accretion disk coronae. The models are calculated using a new method for computing monochromatic radiative transfer n two dimensions. The method splits the radiation into direct and scattered components. The direct radiation is computed by calculating the optical depth along rays, while transfer of the scattered radiation is approximated by flux-limited diffusion. The resulting code agrees with more accurate treatments to within 50%, and is highly efficient, making it practical for use in large hydrodynamic simulations. The coronal models are used to confirm the results of earlier work, and to extend it to higher luminosities. In contrast to earlier work, which found the outer disks to be shadowed by the inner corona at high luminosities, we find our results to form an almost continuous extension of the models at lower luminosities. This is due to the presence of multiply scattered radiation, which acts to partially offset the loss of direct radiation from the central source. Although the analytic methods derived at lower luminosities cannot be used to derive the coronal structure for L/L(sub Edd) approx. greater than 0.1, the results of the models are amenable to semiempirical fits. We also discuss possible observational consequences of the results for coronal veiling and line fluorescence from the disk.

Murray, Stephen D.↗
NASA NTRS1998

Controlling the Mobility of the Fluid Interface of Moving Gas Bubbles or Liquid Drops by Using Micellar Solutions of Surfactants

Microgravity processes must rely on mechanisms other than buoyancy to move bubbles or droplets from one region to another in a continuous liquid phase. One suggested method is thermocapillary migration in which a temperature gradient is applied to the continuous phase. A significant and as yet unresolved impediment to the use of thermocapillary migration to direct bubble or drop motion is that these migrations can be significantly retarded by the adsorption onto the fluid particle surface of surface active impurities unavoidably present in the continuous or (if the particle is a liquid) droplet phases. The focus of our research was to develop a theory for remobilizing fluid particle interfaces retarded by a surfactant impurity in an effort to make more viable the use of thermocapillary migrations for the management of bubbles and drops in microgravity. We postulated that a surfactant at high bulk concentration which kinetically exchanges rapidly with the surface can restore interface mobility. The scaling arguments along with a discussion of the previous literature is reviewed in the context of the scaling framework. The specific objectives of the research were twofold. The first was to prove the remobilization theory by studying a model problem. As the mechanism for remobilization is independent of the force which drives the particle, the fluid particle shape and the presence of fluid inertia, we chose the simplest model consisting of a spherical bubble rising steadily by buoyancy in creeping flow. We solved the hydrodynamic and surfactant transport equations for rapid kinetic exchange to demonstrate that as the concentration increases, the Marangoni retardation at first increases (the low k behavior) and then decreases (the high k behavior). The second objective was to develop a method to determine the kinetic rate constants of a surfactant molecule, since this information is necessary to select surfactants which will exchange rapidly enough relative to the convective rate in the thermocapillary process of interest. To measure the kinetic rate, we measure the dynamic tension change accompanying adsorption onto an initially clean interface, or the re-equilibration in tension when an equilibrium interface is compressed. The dynamic tension measurements are made by a pendant bubble method, in which surfactant adsorbs onto a pendant bubble, and the tension is measured by analyzing the shape change in the bubble. We conclude this report by detailing the publications, presentations and doctoral thesis completed under the auspices of this grant.

Maldarelli, Charles↗
DOE OSTI2024

A modern concept of Lagrangian hydrodynamics

Here, we offer a modern interpretation of Lagrangian hydrodynamics as employed in Lagrangian simulations of compressible fluid flow. Our main result is to show that artificial viscosity, traditionally viewed as a numerical artifice to control unphysical oscillations in flows with shocks, actually represents a physical process and is necessary to derive accurate simulations in any compressible flow. We begin by reviewing the origins of two numerical devices, artificial viscosity and finite-volume methods. We proceed to construct a mathematical (PDE) model that incorporates those numerics and in which a new length scale, the observer, arises representing the discretization. Associated with that length scale, there are new inviscid fluxes that are the artificial viscosity as first formulated by Richtmyer and an artificial heat flux postulated by Noh but typically not included in Lagrangian codes. We discuss the connection of our results to bivelocity hydrodynamics. We conclude with some speculation as to the direction of future developments in multidimensional Lagrangian codes as computers get faster and have larger memories.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

AN OPEN-SOURCE NEUTRINO RADIATION HYDRODYNAMICS CODE FOR CORE-COLLAPSE SUPERNOVAE

We present an open-source update to the spherically symmetric, general-relativistic hydrodynamics, core-collapse supernova (CCSN) code GR1D. The source code is available at http://www.GR1Dcode.org. We extend its capabilities to include a general-relativistic treatment of neutrino transport based on the moment formalisms of Shibata et al. and Cardall et al. We pay special attention to implementing and testing numerical methods and approximations that lessen the computational demand of the transport scheme by removing the need to invert large matrices. This is especially important for the implementation and development of moment-like transport methods in two and three dimensions. A critical component of neutrino transport calculations is the neutrino–matter interaction coefficients that describe the production, absorption, scattering, and annihilation of neutrinos. In this article we also describe our open-source neutrino interaction library NuLib (available at http://www.nulib.org). We believe that an open-source approach to describing these interactions is one of the major steps needed to progress toward robust models of CCSNe and robust predictions of the neutrino signal. We show, via comparisons to full Boltzmann neutrino-transport simulations of CCSNe, that our neutrino transport code performs remarkably well. Furthermore, we show that the methods and approximations we employ to increase efficiency do not decrease the fidelity of our results. We also test the ability of our general-relativistic transport code to model failed CCSNe by evolving a 40-solar-mass progenitor to the onset of collapse to a black hole.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Scale bridging damage model for quasi-brittle metals informed with crack evolution statistics

Computationally efficient methods for bridging length scales, from highly resolved micro/meso-scale models that can explicitly model crack growth, to macro-scale continuum models that are more suitable for modeling large parts, have been of interest to researchers for decades. In this work, an improved brittle damage model is presented for the simulation of dynamic fracture in continuum scale quasi-brittle metal components. Crack evolution statistics, including the number, length, and orientation of individual cracks, are extracted from high-fidelity, finite discrete element method (FDEM) simulations and used to generate effective material moduli that reflect the material’s damaged state over time. This strategy allows for the retention of small-scale physical behaviors such as crack growth and coalescence in continuum scale hydrodynamic simulations. However, the high-fidelity simulations required to generate the crack statistics are computationally expensive. Thus, steps were taken to produce a flexible constitutive model to reduce the number of costly high-fidelity simulations needed to produce accurate results. A new stress based degradation criterion is introduced for the degradation of individual material zones. This allows for the development of a heterogeneous damage distribution within the bulk material. Then a flow stress model is added to the hydrodynamic simulation to account for plasticity in quasi-brittle materials. As a result, the effective moduli model can be applied to a larger range of materials. The effective moduli constitutive model is used to simulate beryllium flyer plate experiments. The results from the continuum scale simulations using statistics from a single high-fidelity simulation are found to be in excellent agreement with numerical and experimental velocity interferometer data. The same set of crack statistics are used to extrapolate the results of a higher rate flyer plate case using the effective moduli model. In conclusion, the extension of this model to higher rate cases shows promise for further reducing the number of costly high-fidelity simulations needed to generate crack statistics.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Convergence of shock waves between conical and parabolic boundaries

Convergence of shock waves, generated by underwater electrical explosions of cylindrical wire arrays, between either parabolic or conical bounding walls is investigated. A high-current pulse with a peak of ∼550 kA and rise time of ∼300 ns was applied for the wire array explosion. Strong self-emission from an optical fiber placed at the origin of the implosion was used for estimating the time of flight of the shock wave. 2D hydrodynamic simulations coupled with the equations of state of water and copper showed that the pressure obtained in the vicinity of the implosion is ∼7 times higher in the case of parabolic walls. However, comparison with a spherical wire array explosion showed that the pressure in the implosion vicinity in that case is higher than the pressure in the current experiment with parabolic bounding walls because of strong shock wave reflections from the walls. It is shown that this drawback of the bounding walls can be significantly minimized by optimization of the wire array geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Design, Development, and Operation of an Integrated Fluidized Carbon Capture Unit Using Polyethylenimine Sorbents

In this paper, we present the design, development, and operation of a reactor system for CO 2 capture. Modifications were implemented to address differences in sorbent from 180 μm Geldart group B to 115 μm Geldart group A material; operational issues were discovered during experimental trials. The major obstacle in system operation was the ability to maintain a constant circulation of a solid sorbent stemming from this change in sorbent material. The system consisted of four fluid beds, through which a polyamine impregnated sorbent was circulated and adsorption, preheat, regeneration, and cooling processes occurred. Pressure transducers, thermocouples, gas flow meters, and gas composition instrumentation were used to characterize thermal, hydrodynamic, and gas adsorption performance in this integrated unit. A series of shakedown tests were performed and the configuration altered to meet the needs of the sorbent performance and achieve desired target capture efficiencies. Methods were identified, tested, and applied to continuously monitor critical operating parameters including solids circulation rate, adsorbed and desorbed CO 2 , solids inventories, and pressures. The working capacity and CO 2 capture efficiency were used to assess sorbent performance while CO 2 closure was used to define data quality and approach to steady-state. Testing demonstrated >90% capture efficiencies and identified the regenerator to be the process step limiting throughput. Sorbent performance was found to be related to the reactant stoichiometry. Finally, a stochastic model with an exponential dependence on the relative CO 2 /amine concentration was used to describe 90% of the variance in the data.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2026

Towards Robust Calibration of the AWSD Reactive Burn Model

Calibration of a reactive burn model for detonation of high explosive is an important step towards predictive hy drodynamic simulations of detonation. A typical calibration consists of varying model parameters (e.g., rate constants, activation energies) until results of hydrodynamic simulations match the experimental data for a certain set of ex periments. Hydrodynamic simulations of the dependence of steady detonation velocity on the radius of a cylindrical high-explosive charge - often used in such calibrations - can be computationally expensive. In this work, we propose a method where such expensive simulations are performed infrequently, and only to parameterize and refine a surrogate model for the dependence of the detonation velocity on calibrated parameters. The method is developed, implemented and applied to an example problem - calibration of the AWSD reactive burn model for important high explosive PBX 9502. Two different flavors of the surrogate model are investigated, and the calibration is performed successfully.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Cholla-MHD: An Exascale-capable Magnetohydrodynamic Extension to the Cholla Astrophysical Simulation Code

Abstract We present an extension of the massively parallel, GPU native, astrophysical hydrodynamics code Cholla to magnetohydrodynamics (MHD). Cholla solves the ideal MHD equations in their Eulerian form on a static Cartesian mesh utilizing the Van Leer + constrained transport integrator, the HLLD Riemann solver, and reconstruction methods at second and third order. Cholla’s MHD module can perform ≈260 million cell updates per GPU-second on an NVIDIA A100 while using the HLLD Riemann solver and second order reconstruction. The inherently parallel nature of GPUs combined with increased memory in new hardware allows Cholla’s MHD module to perform simulations with resolutions ∼500 3 cells on a single high-end GPU (e.g., an NVIDIA A100 with 80 GB of memory). We employ GPU direct Message Passing Interface to attain excellent weak scaling on the exascale supercomputer Frontier, while using 74,088 GPUs and simulating a total grid size of over 7.2 trillion cells. A suite of test problems highlights the accuracy of Cholla’s MHD module and demonstrates that zero magnetic divergence in solutions is maintained to round off error. We also present new testing and CI tools using GoogleTest, GitHub Actions, and Jenkins that have made development more robust and accurate and ensure reliability in the future.

Astronomy & Astrophysics↗
NASA NTRS2017

Overview of the AVT-191 Project to Assess Sensitivity Analysis and Uncertainty Quantification Methods for Military Vehicle Design

A NATO symposium held in 2008 identified many promising sensitivity analysis and un-certainty quantification technologies, but the maturity and suitability of these methods for realistic applications was not known. The STO Task Group AVT-191 was established to evaluate the maturity and suitability of various sensitivity analysis and uncertainty quantification methods for application to realistic problems of interest to NATO. The program ran from 2011 to 2015, and the work was organized into four discipline-centric teams: external aerodynamics, internal aerodynamics, aeroelasticity, and hydrodynamics. This paper presents an overview of the AVT-191 program content.

Benek, John A.↗
DOE OSTI2016

SANTA BARBARA CLUSTER COMPARISON TEST WITH DISPH

The Santa Barbara cluster comparison project revealed that there is a systematic difference between entropy profiles of clusters of galaxies obtained by Eulerian mesh and Lagrangian smoothed particle hydrodynamics (SPH) codes: mesh codes gave a core with a constant entropy, whereas SPH codes did not. One possible reason for this difference is that mesh codes are not Galilean invariant. Another possible reason is the problem of the SPH method, which might give too much “protection” to cold clumps because of the unphysical surface tension induced at contact discontinuities. In this paper, we apply the density-independent formulation of SPH (DISPH), which can handle contact discontinuities accurately, to simulations of a cluster of galaxies and compare the results with those with the standard SPH. We obtained the entropy core when we adopt DISPH. The size of the core is, however, significantly smaller than those obtained with mesh simulations and is comparable to those obtained with quasi-Lagrangian schemes such as “moving mesh” and “mesh free” schemes. We conclude that both the standard SPH without artificial conductivity and Eulerian mesh codes have serious problems even with such an idealized simulation, while DISPH, SPH with artificial conductivity, and quasi-Lagrangian schemes have sufficient capability to deal with it.

79 ASTRONOMY AND ASTROPHYSICS↗