Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics simulations”

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

Modeling the impact of extreme weather events and future climate on the radiologically contaminated sites of Enewetak Atoll

Enewetak Atoll underwent 43 historical nuclear tests from 1948 to 1958, including the first hydrogen bomb test, resulting in a substantial nuclear material fallout contaminating the Atoll and the lagoon waters. The radionuclide fallout material deposited in lagoon sediments and land soil will remain for decades to come. With intensifying climate and extreme weather events, the possibility of redistribution of deposited radionuclide material has become a great concern. This study uses a numerical modeling approach to estimate the potential elevated radionuclide concentrations that can be distributed during storm events under current and future climates. We simulated three historical storm scenarios that are most likely to impact Atoll’s environment and remobilize the radionuclide-bound sediments. WRF-ARW was used to reconstruct these storm scenarios under current year (2015) and future year (2090) climates. Storm-induced ocean hydrodynamics conditions were generated using FVCOM. FVCOM-ICM was externally coupled to simulate the fate and transport of radionuclides. Given that the 239 Pu is the largest inventory of the lagoon and Atoll islands, the model results show the highest average 239 Pu concentration that an island may be exposed to is 3.25E-4 Bq/m 3 (becquerel per cubic meters), which is an increase of 84 times the average baseline/existing 239 Pu concentration without the storm conditions. The overall increase in 239 Pu average over all the islands of Atoll is about 20 folds relative to the baseline concentration. Despite the high relative increase ratios, the significantly low activity concentrations may not pose an immediate exposure risk. However, due to the limitations of the study and uncertainties/biases in the historical data used, further research supported by field surveys to better characterize the current contamination level may be needed to make more accurate predictions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Shock wave formation in radiative plasmas

The temporal evolution of weak shocks in radiative media is theoretically investigated in this work. The structure of radiative shocks has traditionally been studied in a stationary framework. Their systematic classification is complex because layers of optically thick and thin regions alternate to form a radiatively-driven precursor and a temperature-relaxation layer, between which the hydrodynamic shock is embedded. In this work, we analyze the formation of weak shocks when two radiative plasmas with different pressures are put in contact. Applying a reductive perturbative method yields a Burgers-type equation that governs the temporal evolution of the perturbed variables including the radiation field. The conditions upon which optically thick and thin solutions exist have been derived and expressed as a function of the shock strength and Boltzmann number. Below a certain Boltzmann number threshold, weak shocks always become optically thick asymptotically in time, while thin solutions appear as transitory structures. The existence of an optically thin regime is related to the presence of an overdense layer in the compressed material. Scaling laws for the characteristic formation time and shock width are provided for each regime. The theoretical analysis is supported by FLASH simulations, and a comprehensive testcase has been designed to benchmark radiative hydrodynamic codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Validated Finite Element Modeling Tool for Hydrodynamic Loading and Structural Analysis of Ocean Deployed Macroalgae Farms

The objective of this project was to develop and validate a fine-tuned, 3D computational modeling tool for high-fidelity simulations of macroalgae cultivation and harvest systems. To develop this tool, a full understanding of the geometric, material and hydrodynamic properties associated with farming kelp in exposed environments was necessary. Experiments and numerical modeling were done at increasing scales from 1) tank tests with a 1 m model, 2) to field tests with a 61 m culture line at a sheltered site, 3) to field tests with a 122 m culture line at an exposed site and to 4) a 5- line array at the exposed site with 186 m of culture length. The results of the experimental tank tests yielded a set of hydrodynamic drag-area coefficients for a 1 m aggregate based on densely grown kelp at 16 kg/m with a length up to 3 m as described in Fredriksson et al. This work provided the basis to model kelp as a 1 m aggregate based on in-situ geometric and material properties to include measurements of kelp blade length and width, number of blades, yield per m, material mass density and cantilever tests to obtain modulus of elasticity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Energy conserving quadrature based dimensionality reduction for nonlinear hydrodynamics problems

Projection based dimensionality reduction is used to approximate the full discretiza tion of a PDE problem in order to reduce the cost of its numerical simulation while staying faithful to the full dynamics. For nonlinear problems, hyperreduction methods are additionally needed to remove the dependence of the nonlinear terms on the full problem’s size. In this project, we develop an energy conserving hyperreduction method and apply it to Eulerian-Lagrangian hydrodynamics problems, demonstrating that the reduced model successfully accelerates the numerical simulations while conserving the problem’s total energy with high accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Demonstration of a new unstructured mesh IMC x-ray transport capability in LAP codes

The Advanced Simulation and Computing (ASC) Transport project’s Jayenne Implicit Monte Carlo (IMC) transport library now includes an unstructured mesh capability and is available in a Lagrangian Applications Project (LAP) code. In this presentation, we discuss recent work by the LAP and Transport projects that provides an IMC transport capability for radiation hydrodynamics in the Lagrangian frame. Verification problems testing the new capabilities have been simulated and analyzed, i.e. Marshak wave, Mach 45, Su-Olsen, picket fence, and crooked pipe, both in one and two dimensions. We also present results on two stretch goal problems: a simplified COAX high energy density physics experiment and a supernova shock simulation. Finally, we identify current limitations and future work needed to bring a full capability to the user community.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Demonstrating the viability of Lagrangian in situ reduction on supercomputers

Performing exploratory analysis and visualization of large-scale time-varying computational science applications is challenging due to inaccuracies that arise from under-resolved data. In recent years, Lagrangian representations of the vector field computed using in situ processing are being increasingly researched and have emerged as a potential solution to enable exploration. However, prior works have offered limited estimates of the encumbrance on the simulation code as they consider “theoretical” in situ environments. Further, the effectiveness of this approach varies based on the nature of the vector field, benefitting from an in-depth investigation for each application area. With this study, an extended version of Sane et al. (2021), we contribute an evaluation of Lagrangian analysis viability and efficacy for simulation codes executing at scale on a supercomputer. We investigated previously unexplored cosmology and seismology applications as well as conducted a performance benchmarking study by using a hydrodynamics mini-application targeting exascale computing. Here, to inform encumbrance, we integrated in situ infrastructure with simulation codes, and evaluated Lagrangian in situ reduction in representative homogeneous and heterogeneous HPC environments. To inform post hoc accuracy, we conducted a statistical analysis across a range of spatiotemporal configurations as well as a qualitative evaluation. Additionally, our study contributes cost estimates for distributed-memory post hoc reconstruction. In all, we demonstrate viability for each application — data reduction to less than 1% of the total data via Lagrangian representations, while maintaining accurate reconstruction and requiring under 10% of total execution time in over 90% of our experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Fracton Hydrodynamics without Time-Reversal Symmetry

Here, we present an effective field theory for the nonlinear fluctuating hydrodynamics of a single conserved charge with or without time-reversal symmetry, based on the Martin-Siggia-Rose formalism. Applying this formalism to fluids with only charge and multipole conservation, and with broken time-reversal symmetry, we predict infinitely many new dynamical universality classes, including some with arbitrarily large upper critical dimensions. Using large scale simulations of classical Markov chains, we find numerical evidence for a breakdown of hydrodynamics in quadrupole-conserving models with broken time-reversal symmetry in one spatial dimension. Our framework can be applied to the hydrodynamics around stationary states of open systems, broadening the applicability of previously developed ideas and methods to a wide range of systems in driven and active matter.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

OC7 phase I: Toward practical sea-state-dependent modeling of hydrodynamic viscous drag and damping

Here, this article presents a collaborative research campaign under the OC7 project on refining the engineering modeling approach for hydrodynamic viscous drag and damping on floating wind platforms, focusing on the adjustment of hydrodynamic drag and damping coefficients for different sea states. The participant simulation results show significant improvements over the previous OC6 project in predicting the low-frequency resonance motion under nonoperational conditions. The improvements are mainly due to enhanced modeling, including the adoption of wave stretching, and directly tuning the coefficients to measured platform motion in waves instead of free decay. For accurate predictions of mean- and slow-drift motion, the better performing models use a decreasing column splash zone drag coefficient and increasing surge damping/drag with increasing wave height. The model tuning for heave and pitch resonance shows less consistency. Generally, both quadratic drag and additional heave or pitch damping are needed for accurate predictions. Alternatively, a quadratic drag formulation with velocity filtering for the rectangular pontoons leads to improved predictions without additional damping. This model is also potentially more predictive, requiring minimal adjustment to its parameters for different conditions.

17 WIND ENERGY↗
DOE OSTI2020

Modeling hydrodynamics, magnetic fields, and synthetic radiographs for high-energy-density plasma flows in shock-shear targets

Three-dimensional FLASH radiation-magnetohydrodynamics (radiation-MHD) modeling is carried out to study the hydrodynamics and magnetic fields in the shock-shear derived platform. Here, simulations indicate that fields of tens of Tesla can be generated via the Biermann battery effect due to vortices and mix in the counterpropagating shock-induced shear layer. Synthetic proton radiography simulations using MPRAD and synthetic X-ray image simulations using SPECT3D are carried out to predict the observable features in the diagnostics. Quantifying the effects of magnetic fields in inertial confinement fusion and high-energy-density plasmas represents frontier research that has far-reaching implications in basic and applied sciences.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Two-dimensional turbulence in a massless fluid with a relativistic lattice Boltzmann model

We investigate a relativistic adaptation of the lattice Boltzmann method that reproduces the equations of motion for a turbulent, two-dimensional, massless hydrodynamic system. The classical lattice Boltzmann method and its extension to relativistic fluid dynamics are described. The numeric formulation is evaluated using a zero-averaged stirring force introduced into the numerics to induce turbulence, and the flow characteristics produced are compared to properties of a classical turbulent hydrodynamic flow. The model can reasonably be expected to offer quantitative simulations of electron fluid flows in graphene or Kagome lattices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling electrokinetic flows with the discrete ion stochastic continuum overdamped solvent algorithm

In this article we develop an algorithm for the efficient simulation of electrolytes in the presence of physical boundaries. In previous work the discrete ion stochastic continuum overdamped solvent (DISCOS) algorithm was derived for triply periodic domains, and was validated through ion-ion pair correlation functions and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In extending this approach to include an accurate treatment of physical boundaries we must address several important issues. First, the modifications to the spreading and interpolation operators necessary to incorporate interactions of the ions with the boundary are described. Next we discuss the modifications to the electrostatic solver to handle the influence of charges near either a fixed potential or dielectric boundary. An additional short-ranged potential is also introduced to represent interaction of the ions with a solid wall. Finally, the dry diffusion term is modified to account for the reduced mobility of ions near a boundary, which introduces an additional stochastic drift correction. Herein, several validation tests are presented confirming the correct equilibrium distribution of ions in a channel. Additionally, the methodology is demonstrated using electro-osmosis and induced-charge electro-osmosis, with comparison made to theory and other numerical methods. Notably, the DISCOS approach achieves greater accuracy than a continuum electrostatic simulation method. We also examine the effect of under-resolving hydrodynamic effects using a “dry diffusion” approach, and find that considerable computational speedup can be achieved with a negligible impact on accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Compressible pairwise interaction extended point-particle model for force prediction of shock-particle bed interaction

We propose a pairwise influence framework for the complex unsteady compressible particle-laden flow problem by accounting for the scattered hydrodynamic waves emitting from neighboring particles in a Euler-Lagrange simulation. It has been observed from particle-resolved (PR) simulations of randomly dispersed particle beds under a loading shock that the compressible pseudoturbulence dominates the flow system even after the primary shock has passed, which causes fluctuations observed in the forces experienced by the particle. Moreover, the fact that each particle exists in the vicinity of a random arrangement of other particles modifies the time history of the drag force experienced by each particle during and after the passage of the shock. First, the scattering flow field due to an incoming shock interacting with a single sphere is constructed using an analysis of the flow in the acoustic limit. Then we examine the validity of the compressible Maxey-Riley-Gatignol force model by comparing the force prediction against a PR simulation of two interacting particles for various particle arrangements and incoming shock strength. Subsequently, the neighboring influences are stored as a library of maps that can be used readily in the calculation of the perturbation force. Lastly, the pairwise interaction assumption is evaluated by comparing the force predicted with the model with PR simulations of a randomly packed particle bed of 10% volume fraction for both water and air as the fluid medium for an incoming shock Mach number 1.22. With a considerably lower cost for the implementation of the model compared to PR simulations, it is verified that the model is reasonably accurate in pinpointing particles whose peak force is significantly larger or smaller than the mean drag but also to capture the prolonged fluctuations after the initial shock.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A New Method for Ice-Ice Aggregation in the Adaptive Habit Model

In this paper, a novel methodology for modeling ice-ice aggregation is presented. This methodology combines a modified hydrodynamic collection algorithm with bulk aggregate characteristic information from an offline simulator that collects ice particles, namely, the Ice Particle and Aggregate Simulator, and has been implemented into the Adaptive Habit Microphysics scheme in the Weather Research and Forecasting Model. Aggregates, or snow, are formed via collection of cloud ice particles, where initial ice characteristics and the resulting geometry determine aggregate characteristics. Upon implementation, idealized squall-line simulations are performed to examine the new methodology in comparison with commonly used bulk microphysics schemes. It is found that the adaptive habit aggregation parameterization develops snow and reduces ice mass and number concentrations compared to other schemes. The development of aggregates through the new methodology cascades into other interesting effects, including enhancements in ice and snow growth, as well as homogeneous freezing. Further microphysical analyses reveal varying sensitivities, where snow processes are most sensitive to the new parameterization, followed by ice, then cloud, rain, and graupel processes. Further, the new scheme results in enhancements in surface precipitation due to the persistence of snow at lower altitudes. This persistence is a result of shape-dependent melting and sublimation, increasing the residence time. Moreover, these low-level enhancements are reflected in increases in radar reflectivity at the surface and its spatial distribution. Finally, the ability to predict snow shape and density allows for the simulation of polarimetric radar quantities, resulting in signature enhancements compared to schemes that do not consider spatial and temporal variations in snow shape and density.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Understanding Transients: Needs from the Laboratory Astrophysics Community

A growing number of astrophysical transients are pushing astronomers to develop increasingly complex computational tools to model both radiation hydrodynamics and electron transport. Here we review the physics needs and simulation uncertainties associated with modeling these transients. Although this review will focus on theory and simulation aspects of this problem, we also review some experiments designed to study this physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Larval connectivity for European green crab management in the Salish Sea and surrounding waters

The presence of invasive species is a growing concern in coastal marine ecosystems because of their adverse effects on biodiversity. The European green crab Carcinus maenas (EGC) is a small crab inhabiting inshore areas. Although it is native to the Northeast Atlantic Ocean and Baltic Sea, its distribution has expanded to North America, where it is an invasive species. Its main food sources are small invertebrate species that support valuable fisheries in the USA. The first presence of EGC in northern Washington was observed around 20 yr ago, along the Pacific Coast of the USA. Recently, EGC has been detected throughout the Salish Sea (Washington, USA, and British Columbia, Canada) wherein spread dynamics are unknown. The overall distribution of EGC is mainly driven by larval dispersal and, in the Salish Sea and surrounding waters, the assessment of EGC population dynamics is essential to understand its migration patterns and prevent its future expansion. To investigate the dispersal patterns of EGC larvae, a larval dispersal model was developed which couples a regional model of hydrodynamic circulation with an individual-based model of ichthyoplankton dynamics. Simulations were performed over 9 yr (2013-2022) to analyze average larval transport trends in the Salish Sea and surrounding waters, interannual variability of EGC larval connectivity, and the influence of larval behavior on connectivity patterns. Lastly, areas were identified to inform invasive species management moving forward. The prediction of likely sources and settlement locations of EGC larvae from the model will help improve the management of the population in the Salish Sea and surrounding waters.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Cosmic recombination in the presence of primordial magnetic fields

Primordial magnetic fields (PMFs) may explain observations of magnetic fields on extragalactic scales. They are most cleanly constrained by measurements of cosmic microwave background radiation (CMB) anisotropies. Their effects on cosmic recombination may even be at the heart of the resolution of the Hubble tension. We present the most detailed analysis of the effects of PMFs on cosmic recombination to date. To this end we extend the public magneto-hydrodynamic code ENZO with a new cosmic recombination routine, Monte-Carlo simulations of Lyman-α photon transport, and a Compton drag term in the baryon momentum equation. The resulting code allows us, for the first time, to realistically predict the impact of PMFs on the cosmic ionization history and the clumping of baryons during cosmic recombination. Our results identify the importance of mixing of Lyman-α photons between overdense- and underdense- regions for small PMF strength. This mixing speeds up recombination beyond the speed-up due to clumping. We also investigate the effects of pecuilar flows on the recombination rate and find it to be small for small PMF strengths. For non-helical PMFs with a Batchelor spectrum we find a surprising dependency of results on ultra-violet magnetic modes. We further show that the increase in the ionization fraction at low redshift by hydrodynamic baryon heating due to PMF dissipation is completely compensated by the faster recombination from baryon clumping. In conclusion, the present study shall serve as a theoretical foundation for a future precise comparison of recombination with PMFs to CMB data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

2D end-to-end modelling of kilonovae from binary neutron star merger remnants

We investigate the kilonova emission resulting from outflows produced in a 3D general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. We map the outflows into the flash hydrodynamics code to model their expansion in axisymmetry, and study the effects of employing different r-process heating rates. Except for the highest heating rate prescription, we find no significant differences with respect to overall ejecta dynamics and morphology compared to the simulation without heating. Once homologous expansion is attained, typically after ~2s for these ejecta, we map the outflows to the sedona radiative transfer code and compute the spectral evolution of the kilonova and broad-band light curves in various Legacy Survey of Space and Time (LSST) bands. The kilonova properties depend on the remnant lifetime, with peak luminosities and peak time-scales increasing for longer lived remnants that produce more massive ejecta. For all models, there is a strong dependence of both the bolometric and broad-band light curves on the viewing angle. While the short-lived (12 ms) remnant produces higher luminosities when viewed from angles closer to the pole, longer lived remnants (240 ms and 2.5 s) are more luminous when viewed from angles closer to the equator. Our results highlight the importance of self-consistent, long-term modelling of merger ejecta, and taking viewing-angle dependence into account when interpreting observed kilonova light curves. We find that magnetized outflows from an HMNS – if it survives long enough – could explain blue kilonovae, such as the blue emission seen in AT2017gfo.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Manipulation of self-organized multi-vortical states in active magnetic roller suspensions

Ensembles of active magnetic colloids demonstrate complex collective behavior and self-organization when driven out of equilibrium by external magnetic fields. The interplay between magnetic and hydrodynamic interactions often result in spontaneous self-assembly and the emergence of localized vortices in ensembles of ferromagnetic rollers. We report experimental and computational study of the self-organized multi-vortical state in ensembles of magnetic rollers under the imposed external confining potential realized by patterned substrates. We explore the behavior of the system on a substrate patterned with shallow wells ordered in a lattice with square symmetry and lattice spacing smaller than the wells’ diameter. Here we demonstrate, that ensembles of active magnetic rollers evolve in response to changes in the geometry of the confining lattice from a globally correlated state characterized by anti-ferromagnetic vortex ordering to a state with rapidly changing particle flows and self-organized vortical states with a lattice constant larger than one of the confining pattern. Our experimental observations are accompanied by numerical simulations based on phenomenological coarse grained particle dynamics coupled to Navier–Stokes hydrodynamics in shallow water approximation. The reported results provide insights into the collective behavior of active magnetic rollers under confinement and suggest strategies for the manipulation of collective dynamic states in active magnetic liquids.

36 MATERIALS SCIENCE↗