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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

A subzone reconstruction algorithm for efficient staggered compatible remapping

Staggered-grid Lagrangian hydrodynamics algorithms frequently make use of subzonal discretization of state variables for the purposes of improved numerical accuracy, generality to unstructured meshes, and exact conservation of mass, momentum, and energy. For Arbitrary Lagrangian–Eulerian (ALE) methods using a geometric overlay, it is difficult to remap subzonal variables in an accurate and efficient manner due to the number of subzone–subzone intersections that must be computed. This becomes prohibitive in the case of 3D, unstructured, polyhedral meshes. A new procedure is outlined in this paper to avoid direct subzonal remapping. The new algorithm reconstructs the spatial profile of a subzonal variable using remapped zonal and nodal representations of the data. The reconstruction procedure is cast as an under-constrained optimization problem. Enforcing conservation at each zone and node on the remapped mesh provides the set of equality constraints; the objective function corresponds to a quadratic variation per subzone between the values to be reconstructed and a set of target reference values. Numerical results for various pure-remapping and hydrodynamics tests are provided. Ideas for extending the algorithm to staggered-grid radiation-hydrodynamics are discussed as well as ideas for generalizing the algorithm to include inequality constraints.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Halted-pendulum Relaxation: Application to White Dwarf Binary Initial Data

Studying compact-star binaries and their mergers is integral to determining progenitors for observable transients. Today, compact-star mergers are typically studied via state-of-the-art computational fluid dynamics codes. One such numerical technique, smoothed particle hydrodynamics (SPH), is frequently chosen for its excellent mass, energy, and momentum conservation. The natural treatment of vacuum and the ability to represent highly irregular morphologies make SPH an excellent tool for the study of compact-star binaries and mergers. For many scenarios, including binary systems, the outcome of simulations is only as accurate as the initial conditions. For SPH, it is essential to ensure that the particles are distributed regularly, representing the initial density profile but without long-range correlations. Particle noise in the form of high-frequency local motion and low-frequency global dynamics must be damped out. Damping the latter can be as computationally intensive as the actual simulation. We discuss a new and straightforward relaxation method, halted-pendulum relaxation (HPR), to remove global oscillation modes of SPH particle configurations. In combination with effective external potentials representing gravitational and orbital forces, we show that HPR has an excellent performance in efficiently relaxing SPH particles to the desired density distribution and removing global oscillation modes. We compare the method to frequently used relaxation approaches and test it on a white dwarf binary model at its Roche-lobe overflow limit. We highlight the importance of our method in achieving accurate initial conditions and its effect on achieving circular orbits and realistic accretion rates when compared with other general relaxation methods.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

thornado+FLASH-X: A Hybrid Discontinuous Galerkin–Implicit-explicit and Finite-volume Framework for Neutrino-radiation Hydrodynamics in Core-collapse Supernovae

We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement–based multiphysics simulation framework FLASH-X. thornado, developed primarily for simulations of core-collapse supernovae (CCSNe), employs a spectral, six-species two-moment formulation with algebraic closure and special-relativistic observer corrections accurate to $\mathcal{O}(v/c)$, and uses discontinuous Galerkin (DG) methods for phase-space discretization combined with implicit-explicit time stepping. A key development is a nonlinear neutrino–matter coupling algorithm based on nested fixed-point iteration with Anderson acceleration, enabling fully implicit treatment of collisional processes, including energy-coupling interactions such as neutrino–electron scattering and pair production. Coupling to finite-volume (FV) hydrodynamics is achieved through a hybrid DG-FV representation of the fluid variables and operator-split evolution within FLASH-X. The implementation is verified using basic transport tests with idealized opacities and relaxation and deleptonization problems with tabulated microphysics. Spherically symmetric CCSN simulations demonstrate accuracy and robustness of the coupled scheme, including close agreement with the CCSN simulation code Chimera. An axisymmetric CCSN simulation further demonstrates the viability of DG-based neutrino transport for multidimensional supernova modeling within FLASH-X. thornado’s neutrino-transport solver is GPU-enabled using OpenMP offloading or OpenACC, and all CCSN applications included in this work use the GPU implementation. Together, these results establish a foundation for future enhancements in physics fidelity, numerical algorithms, and computational performance, for increasingly realistic large-scale CCSN simulations.

Endeve, Eirik [Oak Ridge National Laboratory (ORNL↗
DOE OSTI2017

A solution algorithm for fluid–particle flows across all flow regimes

Many fluid–particle flows occurring in nature and in technological applications exhibit large variations in the local particle volume fraction. For example, in circulating fluidized beds there are regions where the particles are closepacked as well as very dilute regions where particle–particle collisions are rare. Thus, in order to simulate such fluid–particle systems, it is necessary to design a flow solver that can accurately treat all flow regimes occurring simultaneously in the same flow domain. In this work, a solution algorithm is proposed for this purpose. The algorithm is based on splitting the free-transport flux solver dynamically and locally in the flow. In close-packed to moderately dense regions, a hydrodynamic solver is employed, while in dilute to very dilute regions a kinetic-based finite-volume solver is used in conjunction with quadrature-based moment methods. To illustrate the accuracy and robustness of the proposed solution algorithm, it is implemented in OpenFOAM for particle velocity moments up to second order, and applied to simulate gravity-driven, gas–particle flows exhibiting cluster-induced turbulence. By varying the average particle volume fraction in the flow domain, it is demonstrated that the flow solver can handle seamlessly all flow regimes present in fluid–particle flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Tuning the reactivity of Al/Fe{sub 2}O{sub 3} nanoenergetic materials via an approach combining soft template self-assembly with sol–gel process process

A bottom-up approach combining soft template self-assembly with sol–gel process, was adopted to prepare the assembled Al/Fe{sub 2}O{sub 3} nanoenergetic materials, assembly-Al/Fe{sub 2}O{sub 3} sample. The other two unassembled Al/Fe{sub 2}O{sub 3}a nanoenergetic materials, sol–gel–Al/Fe{sub 2}O{sub 3} sample and mixing-Al/Fe{sub 2}O{sub 3} sample, were prepared by sol–gel method and physical mixing method respectively. The assembly process within the preparation of the assembly-Al/Fe{sub 2}O{sub 3} sample was analyzed through the changes in the average hydrodynamic diameters of the particles and the micelles in solution. SEM, EDS and TEM tests were performed to demonstrate a significant improvement regarding to dispersity and arrangements of the Al and Fe{sub 2}O{sub 3} particles in the assembled samples, compared to that of the unassembled Al/Fe{sub 2}O{sub 3} samples. DSC test was employed to characterize the reactivity of the samples. The heat release of the assembled Al/Fe{sub 2}O{sub 3} sample was 2088 J/g, about 400 and 990 J/g more than that of the sol–gel–Al/Fe{sub 2}O{sub 3} sample and mixing-Al/Fe{sub 2}O{sub 3} sample, respectively. - Graphical abstract: Modified aluminum (Al) nanoparticles with hydrophobic surface assembled into the Brij S10 micelle in Fe(III) sol, then the well dispersed system was transformed into Al/Fe{sub 2}O{sub 3} nanoenergetic materials with high reactivity. - Highlights: • An approach combining soft template self-assembly with sol–gel process was adopted. • The aggregation of Al nanoparticles in the final product was reduced significantly. • The reactivity of Al/Fe{sub 2}O{sub 3} nanoenergetic materials was improved to a large extent.

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

Perspective on ab initio phonon thermal transport

Coupling of the Peierls-Boltzmann equation with density functional theory paved the way for predictive thermal materials discovery and a variety of new physical insights into vibrational transport behaviors. Here, rapid theoretical and numerical developments have generated a wealth of thermal conductivity data and understanding of a wide variety of materials—1D, 2D, and bulk—for thermoelectric and thermal management applications. Nonetheless, modern ab initio descriptions of phonon thermal transport face challenges regarding the effects of defects, disorder, structural complexity, strong anharmonicity, quasiparticle couplings, and time and spatially varying perturbations. Highlighting recent research on these issues, this perspective explores opportunities to expand current ab initio phonon transport techniques beyond the paradigm of weakly perturbed crystals, to the wider variety of materials possible. Recent developments in phonon-defect interactions, complexity, disorder and anharmonicity, hydrodynamic transport, and the rising roles of molecular dynamics simulations, high throughput, and machine learning tools are included in this perspective. Lastly, as more sophisticated theoretical and computational methods continue to advance thermal transport predictions, novel vibrational physics and thermally functional materials will be discovered for improved energy technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Quantitative x ray phase contrast imaging of oblique shock wave–interface interactions

Oblique shock wave–interface interactions of gases and liquids have been extensively studied in shock tubes using optical methods to measure equation-of-state (EOS) parameters. However, this is difficult with solids due to their opaqueness to visible light. X ray phase contrast imaging (XPCI) has the penetrative strength to probe solids while still being sensitive to mass density and enhancing the visibility of material boundaries. We investigate the accuracy and repeatability of measuring the mean value of the average mass density (areal density divided by thickness) over region S (B S ) and flow deflection angle (θ) from XPCI images of a sample. To that end, a Hough transform-based method for measuring θ is developed. To measure B S , the XPCI image intensity probability density function (PDF) is modeled accounting for the spatial distribution of x ray energy, scintillator response, and pulse-to-pulse variation in the x ray intensity. In addition, a Monte Carlo-based algorithm for computing the B S PDF is developed. Both methods are validated on an impact-generated oblique shock wave interacting at a solid polymer-to-polymer interface. Furthermore this is accomplished through a comparison to hydrodynamic simulations using well-established EOS. Under the modeling framework for the XPCI image intensity, B S is computed with an accuracy of < 0.1% and precision of 3%–5%, while θ has an uncertainty of 0.2°, respectively. This shows that the XPCI-based model that is developed here could be an invaluable tool for high-fidelity testing of hydrodynamic models in shock polar configurations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Optimizing the Smoothness and Thickness Uniformity of Thin-Film Parylene-N Vapor-Deposited Coatings for Inertial Confinement Fusion Experiments

Polymer coatings with submicrometer smoothness and constant thickness are a required component in a variety of inertial confinement fusion experiments. Smoothness is important for minimizing Rayleigh-Taylor-driven hydrodynamic instabilities, and uniform thickness is important for uniform shock propagation and shell convergence, both of which are critical phenomena that affect the experiment. The preferred polymer coating method is to vapor deposit the parylene-N polymer because it provides nominally smooth conformal coatings. As the coating thickness exceeds 5 µm, however, dome-shaped nodular growth defects develop and the thickness will vary by up to 17% over a distance of 3 cm. This study presents a deterministic method for achieving uniform film thicknesses with ±2% variability over 3 cm and a predictive method to control the thickness to within 5% of the desired value. A coating smoothness of ∼50 nm rms, measured over 40 000 µm 2 , was achieved by adding additional surfaces near the substrates. This additional area improved the thickness uniformity, an effect that is attributed to the low sticking coefficient of the parylene monomer.

chemical vapor deposition (CVD)↗
DOE OSTI2014

Images of gravitational and magnetic phenomena derived from two-dimensional back-projection Doppler tomography of interacting binary stars

We have used two-dimensional back-projection Doppler tomography as a tool to examine the influence of gravitational and magnetic phenomena in interacting binaries that undergo mass transfer from a magnetically active star onto a non-magnetic main-sequence star. This multitiered study of over 1300 time-resolved spectra of 13 Algol binaries involved calculations of the predicted dynamical behavior of the gravitational flow and the dynamics at the impact site, analysis of the velocity images constructed from tomography, and the influence on the tomograms of orbital inclination, systemic velocity, orbital coverage, and shadowing. The Hα tomograms revealed eight sources: chromospheric emission, a gas stream along the gravitational trajectory, a star-stream impact region, a bulge of absorption or emission around the mass-gaining star, a Keplerian accretion disk, an absorption zone associated with hotter gas, a disk-stream impact region, and a hot spot where the stream strikes the edge of a disk. We described several methods used to extract the physical properties of the emission sources directly from the velocity images, including S-wave analysis, the creation of simulated velocity tomograms from hydrodynamic simulations, and the use of synthetic spectra with tomography to sequentially extract the separate sources of emission from the velocity image. In summary, the tomography images have revealed results that cannot be explained solely by gravitational effects: chromospheric emission moving with the mass-losing star, a gas stream deflected from the gravitational trajectory, and alternating behavior between stream state and disk state. Our results demonstrate that magnetic effects cannot be ignored in these interacting binaries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Features of the Upsurge in Drift Velocity of Electrons in DA-pHEMT

A simple phenomenological model for estimating the upsurge in drift velocity of electrons in transistor heterostructures is proposed. This model is based on a self-consistent solution of Schrödinger and Poisson equations and the hydrodynamic system of equations of energy and momentum conservation. It is demonstrated that the conditions in the layer channel of DA-pHEMT structures with additional potential barriers, which are produced by donor–acceptor doping and enhance the localization of hot electrons, are even more conducive to reducing the time of flight of electrons under the gate than those established in heterostructures with deeper quantum wells produced by increasing the conduction band offset at the heterojunction interface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

A Novel Approach to Fast Calculation of High-Order Q-Cumulants

The method of Q-cumulants has been shown as a powerful tool to study the fine details of the azimuthal anisotropies in high-energy nucleus-nucleus collisions. A new method for the fast calculation of arbitrary order Q-cumulant $v_{n}\{2k\}$ values, based on the partition of a non-negative integer l $\le$ m for calculation of the 2m-particle azimuthal correlations is presented in this paper. Unlike the standard Q-cumulants method in which the calculation of high-order multi-particle calculations is impractical, the newly proposed method enables easy calculation. The validity of the method is proven via a toy model that uses the elliptic power distribution to simulate anisotropic emission of particles. The method enables the study of fine details of the $v_{2}$ distribution, such as higher-order central moments of the $v_{2}$ distribution, as well as the hydrodynamic behavior of the Quark-Gluon Plasma.

FOS: Physical sciences↗
NASA NTRS2005

Computation of Tone Noise From Supersonic Jet Impinging on Flat Plates

A supersonic jet impinging normally on a flat plate has both practical importance and theoretical interests. The physical phenomenon is not fully understood yet. Research concentrates either on the hydrodynamics (e.g., lift loss for STOVL) or on the aeroacoustic loading. In this paper, a finite volume scheme - the space-time conservation element and solution element (CE/SE) method - is employed to numerically study the near-field noise of an underexpanded supersonic jet from a converging nozzle impinging normally on a flat plate. The numerical approach is of the MILES type (monotonically integrated large eddy simulation). The computed results compare favorably with the experimental findings.

Loh, Ching Y.↗
DOE OSTI2021

Modification of the magneto-hydrodynamic equilibrium by the lower-hybrid wave driven fast electrons on the TST-2 spherical tokamak

Presence of a large population of fast particles may qualitatively modify the tokamak equilibrium from the one that can be described by standard magneto-hydrodynamics (MHD). The kinetic modification of the MHD equilibrium was studied for a lower-hybrid (LH) wave driven plasma on the TST-2 spherical tokamak. The analysis was performed using an equilibrium reconstruction method based on an extended MHD model that considers a two-component plasma of bulk MHD fluid and kinetic fast electrons. Scrape-off-layer current appeared naturally in the extended MHD analysis because of the finite electron orbit excursion from the flux surfaces. This reduced the bulk pressure contribution to the toroidal current substantially. Furthermore, the resulting poloidal flux profile was more consistent with that inferred from the Thomson scattering measurement.

extended magneto-hydrodynamics↗
DOE OSTI2016

A unified gas-kinetic scheme for continuum and rarefied flows IV: Full Boltzmann and model equations

Fluid dynamic equations are valid in their respective modeling scales, such as the particle mean free path scale of the Boltzmann equation and the hydrodynamic scale of the Navier–Stokes (NS) equations. With a variation of the modeling scales, theoretically there should have a continuous spectrum of fluid dynamic equations. Even though the Boltzmann equation is claimed to be valid in all scales, many Boltzmann solvers, including direct simulation Monte Carlo method, require the cell resolution to the order of particle mean free path scale. Therefore, they are still single scale methods. In order to study multiscale flow evolution efficiently, the dynamics in the computational fluid has to be changed with the scales. A direct modeling of flow physics with a changeable scale may become an appropriate approach. The unified gas-kinetic scheme (UGKS) is a direct modeling method in the mesh size scale, and its underlying flow physics depends on the resolution of the cell size relative to the particle mean free path. The cell size of UGKS is not limited by the particle mean free path. With the variation of the ratio between the numerical cell size and local particle mean free path, the UGKS recovers the flow dynamics from the particle transport and collision in the kinetic scale to the wave propagation in the hydrodynamic scale. The previous UGKS is mostly constructed from the evolution solution of kinetic model equations. Even though the UGKS is very accurate and effective in the low transition and continuum flow regimes with the time step being much larger than the particle mean free time, it still has space to develop more accurate flow solver in the region, where the time step is comparable with the local particle mean free time. In such a scale, there is dynamic difference from the full Boltzmann collision term and the model equations. This work is about the further development of the UGKS with the implementation of the full Boltzmann collision term in the region where it is needed. The central ingredient of the UGKS is the coupled treatment of particle transport and collision in the flux evaluation across a cell interface, where a continuous flow dynamics from kinetic to hydrodynamic scales is modeled. The newly developed UGKS has the asymptotic preserving (AP) property of recovering the NS solutions in the continuum flow regime, and the full Boltzmann solution in the rarefied regime. In the mostly unexplored transition regime, the UGKS itself provides a valuable tool for the non-equilibrium flow study. The mathematical properties of the scheme, such as stability, accuracy, and the asymptotic preserving, will be analyzed in this paper as well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1967

Computer analysis of large structural systems.

Numerical procedure for structural systems analysis, discussing computer application to hydrodynamic, electric, magnetic, thermodynamic, elastostatic and elastodynamic problems

NUMERICAL ANALYSIS↗
DOE OSTI2015

Development of an Efficient Meso- scale Multi-phase Flow Solver in Nuclear Applications

The proposed research aims at formulating a predictive high-order Lattice Boltzmann Equation for multi-phase flows relevant to nuclear energy related application - namely, saturated and sub-cooled boiling in reactors, and liquid- liquid mixing and extraction for fuel cycle separation. An efficient flow solver will be developed based on the Finite Element based Lattice Boltzmann Method (FE- LBM), accounting for phase-change heat transfer and capable of treating multiple phases over length scales from the submicron to the meter. A thermal LBM will be developed in order to handle adjustable Prandtl number, arbitrary specific heat ratio, a wide range of temperature variations, better numerical stability during liquid-vapor phase change, and full thermo-hydrodynamic consistency. Two-phase FE-LBM will be extended to liquid–liquid–gas multi-phase flows for application to high-fidelity simulations building up from the meso-scale up to the equipment sub-component scale. While several relevant applications exist, the initial applications for demonstration of the efficient methods to be developed as part of this project include numerical investigations of Critical Heat Flux (CHF) phenomena in nuclear reactor fuel bundles, and liquid-liquid mixing and interfacial area generation for liquid-liquid separations. In addition, targeted experiments will be conducted for validation of this advanced multi-phase model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Weak boundary conditions for Lagrangian shock hydrodynamics: A high-order finite element implementation on curved boundaries

Here, we propose a new Nitsche-type approach for weak enforcement of normal velocity boundary conditions for a Lagrangian discretization of the compressible shock-hydrodynamics equations using high-order finite elements on curved boundaries. Specifically, the variational formulation is appropriately modified to enforce free-slip wall boundary conditions, without perturbing the structure of the function spaces used to represent the solution, with a considerable simplification with respect to traditional approaches. Total energy is conserved and the resulting mass matrices are constant in time. The robustness and accuracy of the proposed method are validated with an extensive set of tests involving nontrivial curved boundaries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Global anomalies, discrete symmetries and hydrodynamic effective actions

We derive effective actions for parity-violating fluids in both (3 + 1) and (2 + 1) dimensions, including those with anomalies. As a corollary we confirm the most general constitutive relations for such systems derived previously using other methods. We discuss in detail connections between parity-odd transport and underlying discrete symmetries. In (3+1) dimensions we elucidate connections between anomalous transport coefficients and global anomalies, and clarify a previous puzzle concerning transports and local gravitational anomalies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗