Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Hydrodynamically Controlled Self‐Organization in Mixtures of Active and Passive Colloids

Abstract Active particles are known to exhibit collective behavior and induce structure in a variety of soft‐matter systems. However, many naturally occurring complex fluids are mixtures of active and passive components. The authors examine how activity induces organization in such multi‐component systems. Mixtures of passive colloids and colloidal micromotors are investigated and it is observed that even a small fraction of active particles induces reorganization of the passive components in an intriguing series of phenomena. Experimental observations are combined with large‐scale simulations that explicitly resolve the near‐ and far‐field effects of the hydrodynamic flow and simultaneously accurately treat the fluid–colloid interfaces. It is demonstrated that neither conventional molecular dynamics simulations nor the reduction of hydrodynamic effects to phoretic attractions can explain the observed phenomena, which originate from the flow field that is generated by the active colloids and subsequently modified by the aggregating passive units. These findings not only offer insight into the organization of biological or synthetic active–passive mixtures, but also open avenues to controlling the behavior of passive building blocks by means of small amounts of active particles.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Simulation of dynamic crystal plasticity with a Lagrangian discontinuous Galerkin hydrodynamic method

Here we present a new Lagrangian modal discontinuous Galerkin (DG) hydrodynamic method that supports a dynamic dislocation based crystal plasticity model for simulating the mechanical behavior of crystallographic materials, both single crystal and polycrystalline, under dynamic conditions. A modal DG approach is used to evolve fields relevant to conservation laws. These fields are approximated by Taylor series polynomials of varying degree. These polynomials describe macro-scale hydrodynamic behavior while their evolution is determined by evaluating the dynamic crystal plasticity model at material points within the element. The dynamic crystal plasticity model is sensitive to the time increment size, with too large of time increments leading to instability in the model. To mitigate this, the temporal evolution of the dynamic crystal plasticity model is achieved with the combination of a sub-incrementing scheme with Heun’s third-order time integration scheme, which is also used to temporally evolve the governing equations. The implementation of the dynamic crystal plasticity model within the DG framework is tested using a 2D approximation of the Taylor impact test with a single crystal material, using quadratic elements that have faces that can bend. In addition to the standard continuous material modeling, we propose a new simulation method that would represent the heterogeneous behavior of polycrystalline microstructures within an element by varying the position and material properties of the material points within the element. This method is demonstrated using random orientation distributions on materials points that are arranged in both structured and random configurations.

42 ENGINEERING↗
DOE OSTI2020

Locally adaptive artificial viscosity strategies for Lagrangian hydrodynamics

To accurately model inviscid flow with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity is introduced to convert kinetic energy into internal energy, thereby providing a mechanism to generate the required entropy increase across shocks. In this paper, we propose a new method for constructing an adaptive, artificial viscosity in the context of one-dimensional, staggered-grid Lagrangian hydrodynamics. Our adaptive, artificial viscosity is defined in terms of two parameters that depend on density locally in a neighborhood around a shock, and hence, vary cell-by-cell. Our methodology is based on building a reference set of pre-computed optimal, globally constant artificial viscosity coefficients for a family of isolated shock test problems. For arbitrary flows, the evaluation of the unknown coefficients is automated by first estimating shock intensity locally, and second computing the corresponding adaptive parameter value via interpolation over the data from the pre-computed reference set of optimal, globally constant coefficient values. To illustrate the performance of our new approach, we compare our results against two existing methods. The first method is a limiter-based approach, which relies on estimating velocity gradients of the flow, and the second method is utilized in several commercial codes. Finally, we demonstrate that our new adaptive methodology produces more accurate results for a variety of tests with propagating shock waves, as well as for the aforementioned family of isolated shock problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

3D Cell-centered hydrodynamics with subscale closure model and multi-material remap

We extend a higher-order finite volume cell-centered hydrodynamic (CCH) formulation to include an interface-aware subscale closure model and a multi-material remap for simulating 3D compressible hydrodynamic problems within an arbitrary Lagrangian-Eulerian (ALE) framework. This CCH formulation involves a multidirectional approximate Riemann solution using quadratic polynomial reconstructions of the stress tensor and the velocity. At the subscale level, we determine pair-wise material interactions by solving a distinct approximate Riemann problem at the common interface, using the volume of fluids (VOF) method to find the interface. Material interactions are constrained to ensure smooth pressure equilibration among materials. The accuracy and robustness of the ALE method is demonstrated by simulating a suite of 3D Cartesian multi-material problems covering both gas and solid dynamics, where each test case has two or more materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

An adaptive model-free robotic force control strategy for hydrodynamic real-time hybrid simulation of floating offshore wind turbines

Real-time hybrid simulation (RTHS) - a cyber-physical testing approach - promises to enhance the simulation fidelity of the model-scale experiments used to prototype floating offshore wind turbines (FOWTs). In hydrodynamic RTHS (hydro-RTHS), actuators emulate aerodynamic forces on model-scale FOWT specimens subjected to physical waves in a hydrodynamic laboratory. Robotic arms are promising candidates for actuation in hydro-RTHS due to their compact multi-degree-of-freedom (DOF) capabilities. Unlike classical RTHS for seismic applications, which typically relies on displacement control, hydro-RTHS requires 6-DOF force control on newly designed floating prototypes in a model-scale setting, which presents significant challenges, including modeling uncertainties, directional asymmetry, configuration drift, bandwidth limitations, and time-varying delays. To mitigate these constraints without extensive pre-test calibration, this study proposes an adaptive model-free robotic force control strategy that combines task-space explicit force control with a secondary joint-space pose-keeping task. The Adaptive Feedforward Compensator (AFC) is integrated into the force control loop to compensate for time-varying delay. Experimental testing was conducted using a Franka Emika Panda robotic arm with a 1:50 scale FOWT specimen under operational wind and wave conditions. Results demonstrate stable and consistent 6-DOF force tracking. Effective delay compensation was observed, with low-frequency delay reductions ranging from 71.4% to 91.8% and improvements in low-frequency surge force tracking of 25.0% to 52.1%. This study enhances robotic actuation performance in hydro-RTHS and introduces a force control strategy that supports reliable robotic operation in uncertain floating environments. Future work will explore disturbance-observer mechanisms to further enhance wave rejection capabilities under extreme wind and wave conditions.

17 WIND ENERGY↗
DOE OSTI2025

Towards smarter green infrastructure: Fusing bark ecology and stemflow hydrodynamics on tree stems

A wide array of bark surfaces sheath wooded plants in rural and urban areas alike. Much work has examined the function and role of bark in different contexts and different environs, including urban areas, finding that it is rich in life and can play a role in the transfer of water and matter to the ground surface. Accordingly, this paper presents a first step to weld and fuse bark ecology and stemflow hydrodynamics. It is an effort to develop a physically-based understanding of the transport of water and matter (e.g., solutes, particulates, microorganisms) along tree stems using relevant equations to allow a more informed consideration of bark in green infrastructure initiatives. In particular, the hydrodynamical equations are based on the conservation of water mass, conservation of momentum, and conservation of scalar mass. These equations, coupled with contemplation of corticular life, underpin and substantiate bark’s unifying role as a modulator and cultivator. By elucidating the ‘black box’ of the tree stem and utilizing the formulations set forth in this paper, urban foresters and planners can develop green infrastructure to help advance ecosystem services and sustainability development goals (SDG), especially SDG 11 and SDG 15.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Development of spatially and temporally resolved electron density measurements for the assessment of radiation hydrodynamics simulations of laboratory X-ray photoionized plasmas

The photoionized plasma gas cell experiment is an established platform we use to make at-parameter (ξ > > 1 ergs cm s –1 ) measurements of plasma properties with application to high-energy astrophysical systems. We model the experiments with 1D radiation hydrodynamics simulations using the HELIOS-CR code to inform our understanding and assist in the interpretation of results. The simulations predict that the bulk of the plasma is in a quasi-uniform and hydrodynamically unperturbed state throughout the duration of the experiment. To evaluate this prediction, we introduced a photonic Doppler velocimetry (PDV) diagnostic to measure spatially and temporally resolved plasma electron density. The initial measurements were successful but had limitations that made model-data comparisons challenging. To address this, we re-designed the gas cell PDV diagnostic and doubled the number of measurement locations to sample across two thirds of the depth of the cell. We also present a comparison of the results from the upgraded PDV diagnostic to the HELIOS-CR simulations for the first time. As a result, the experimental data confirms the prediction of an unperturbed region in the bulk of the plasma but reveals discrepancies in the time evolution and spatial distribution of the simulated electron density.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2021

Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star

A realistic parametrization of convection and convective boundary mixing in conventional stellar evolution codes is still the subject of ongoing research. Furthermore, to improve the current situation, multidimensional hydrodynamic simulations are used to study convection in stellar interiors. Such simulations are numerically challenging, especially for flows at low Mach numbers which are typical for convection during early evolutionary stages. We explore the benefits of using a low-Mach hydrodynamic flux solver and demonstrate its usability for simulations in the astrophysical context. Simulations of convection for a realistic stellar profile are analyzed regarding the properties of convective boundary mixing. The time-implicit Seven-League Hydro (SLH) code was used to perform multidimensional simulations of convective helium shell burning based on a 25 M ⊙ star model. The results obtained with the low-Mach AUSM + -up solver were compared to results when using its non low-Mach variant AUSM B + -up. We applied well-balancing of the gravitational source term to maintain the initial hydrostatic background stratification. The computational grids have resolutions ranging from 180 × 90 2 to 810 × 540 2 cells and the nuclear energy release was boosted by factors of 3 × 10 3 , 1 × 10 4 , and 3 × 10 4 to study the dependence of the results on these parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Hydrodynamic simulations of shock-driven chemistry in polyimide

Here, the development for equations of state both for polyimide and its reaction products is presented along with hydrodynamic simulations linking the two EOS through an Arrhenius rate law. The equations of state compare favorably to available data, and the hydrodynamic simulations are able to qualitatively reproduce many of the features seen in both legacy shock data and more modern embedded gauge data. However, quantitative agreement is not achieved, suggesting that new rate laws are needed to fully capture the effect of the reaction over a wide range of conditions.

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

Coupling solidification kinetics with phase-behavior computations in hydrodynamic simulations of high-pressure, dynamic-compression processes

In this study, we report a numerical scheme to integrate models for the kinetics of solidification processes together with phase-behavior computations in the context of continuum-scale hydrodynamic simulations. The objective of the phase-behavior computations is to determine the pressure and temperature, given the following three sets of inputs: (1) an appropriate equation of state to describe our system, (2) the phase fraction(s) produced by the kinetic models, (3) and the volume and internal energy obtained by solving the conservation equations that govern the hydrodynamic behavior. The kinetics are assumed to be governed by the Kolmogorov–Johnson–Mehl–Avrami equation, and the nucleation and growth rates that enter into that equation are functions of the pressure and temperature produced by the phase-behavior computations. Our formulation allows for the fluid and solid phases to be at different temperatures (thermal nonequilibrium) and pressures (arising from surface-tension-induced Laplace contributions). The formulation is presented in a fairly general setting that is independent of any particular material, although we demonstrate it in some examples related to high-energy-density science applications where materials are rapidly compressed to pressures exceeding several gigapascals in less than a microsecond. We conclude with a critical evaluation of our approach and provide suggestions for future work to improve the predictive capabilities and generality of the models.

42 ENGINEERING↗
DOE OSTI2024

Measurements of dense fuel hydrodynamics in the NIF burning plasma experiments using backscattered neutron spectroscopy

The hydrodynamics of the dense confining fuel shell is of great importance in defining the behavior of the burning plasma and burn propagation regimes of inertial confinement fusion experiments. However, it is difficult to probe due to its low emissivity in comparison with the central fusion core. In this work, we utilize the backscattered neutron spectroscopy technique to directly measure the hydrodynamic conditions of the dense fuel during fusion burn. Experimental data are fit to obtain dense fuel velocities and apparent ion temperatures. Trends of these inferred parameters with yield and velocity of the burning plasma are used to investigate their dependence on alpha heating and low mode drive asymmetry. It is shown that the dense fuel layer has an increased outward radial velocity as yield increases, showing that burn has continued into re-expansion, a key signature of hotspot ignition. A comparison with analytic and simulation models shows that the observed dense fuel parameters are displaying signatures of burn propagation into the dense fuel layer, including a rapid increase in dense fuel apparent ion temperature with neutron yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Reaction-in-flight neutrons as a diagnostic for hydrodynamical mixing in double shell inertial confinement fusion capsules

We examine reaction-in-flight (RIF) neutrons as diagnostics for hydrodynamical mixing of high-Z shell material into the hotspot of double shell capsules that are designed for the National Ignition Facility. In particular, we consider the effects of different levels of mixing of tungsten shell material into the DT gas on RIF spectra. Using a set of 1D simulations from the radiation hydrodynamic code xRAGE to determine the temperature and density profiles of the mixed W-DT gas, we find that increasing the mass of mix systematically reduces the ratio of RIF neutrons to primary 14 MeV neutrons. The shape of the RIF spectrum also changes with mix, with the predicted spectrum softening in energy with increasing mix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

ETHOS: An automated framework to generate multi-fidelity constitutive data tables and propagate uncertainties to hydrodynamic simulations

Accurate constitutive data, such as equations of state and plasma transport coefficients, are necessary for reliable hydrodynamic simulations of plasma systems such as fusion targets, planets, and stars. Here, we develop a framework for automatically generating transport-coefficient tables using a parameterized model that incorporates data from both high-fidelity sources (e.g., density functional theory calculations and reference experiments) and lower-fidelity sources (e.g., average-atom and analytic models). The framework incorporates uncertainties from these multi-fidelity sources, generating ensembles of optimally diverse tables that are suitable for uncertainty quantification of hydrodynamic simulations. We illustrate the utility of the framework with magnetohydrodynamic simulations of magnetically launched flyer plates, which are used to measure material properties in pulsed-power experiments. We explore how changes in the uncertainties assigned to the multi-fidelity data sources propagate to changes in simulation outputs and find that our simulations are most sensitive to uncertainties near the melting transition. The presented framework enables computationally efficient uncertainty quantification that readily incorporates new high-fidelity measurements or calculations and identifies plasma regimes where additional data will have high impact.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A fluctuating hydrodynamics model for nanoscale surfactant-laden interfaces

A multispecies diffuse interface model is formulated in a fluctuating hydrodynamics framework for the purpose of simulating surfactant interfaces at the nanoscale. The model generalizes previous work to ternary mixtures, employing a Cahn-Hilliard free energy density combined with incompressible, isothermal fluctuating hydrodynamics where dissipative fluxes include both deterministic and stochastic terms. The intermolecular parameters in the free energy are chosen such that one species acts as a partially miscible surfactant. From Laplace pressure measurements, we show that in this model the surface tension decreases linearly with surfactant concentration, leading to Marangoni convection for interfaces with concentration gradients. In the capillary wave spectrum for interfaces with and without surfactant, we find that for the former, the spectrum deviates significantly from classical capillary wave theory, presumably due to Gibbs elasticity. In non-equilibrium simulations of the Rayleigh-Plateau instability, deterministic simulations showed that the surfactant delays pinching of a fluid cylinder into droplets. However, stochastic simulations indicate that thermal fluctuations disrupt the surfactant's stabilizing effect. Similarly, the spreading of a patch of surfactant, driven by Marangoni convection, was found to be partially suppressed by thermal fluctuations.

Capillary waves↗
DOE OSTI2026

Non-local transport in radiation-hydrodynamics codes for ICF by efficient coupling to an external Vlasov–Fokker–Planck code

Accurately incorporating non-local transport into radiation-hydrodynamics codes, and indeed any fluid system, has long been elusive. To date, a simplified and accurate theory that can be easily integrated has not been available. This limitation affects modeling in inertial confinement fusion (ICF) and magnetic confinement fusion systems, among others, where non-local transport is well-known to be present. Here, we present a coupling methodology between a full Vlasov–Fokker–Planck (VFP) electron kinetic code and radiation-hydrodynamics (rad-hydro) codes. The VFP code is used to adjust native electron transport in the rad-hydro code, thus enabling improved transport without the need to integrate a full electron VFP solver into the rad-hydro code. This approach necessitates only occasional invocation of the VFP code, reducing computational intensity compared to following the dynamic evolution entirely with the VFP code on fluid time scales. We illustrate that the methodology is more accurate than other simplified methods in thermal decay systems relevant to ICF and can replicate standard theoretical results with high accuracy.

Electronic transport↗
DOE OSTI2020

Hydrodynamic and frictional modulation of deformations in switchable colloidal crystallites

Displacive transformations in colloidal crystals may offer a pathway for increasing the diversity of accessible configurations without the need to engineer particle shape or interaction complexity. To date, binary crystals composed of spherically symmetric particles at specific size ratios have been formed that exhibit floppiness and facile routes for transformation into more rigid structures that are otherwise not accessible by direct nucleation and growth. There is evidence that such transformations, at least at the micrometer scale, are kinetically influenced by concomitant solvent motion that effectively induces hydrodynamic correlations between particles. Here, we study quantitatively the impact of such interactions on the transformation of binary bcc-CsCl analog crystals into close-packed configurations. We first employ principal-component analysis to stratify the explorations of a bcc-CsCl crystallite into orthogonal directions according to displacement. We then compute diffusion coefficients along the different directions using several dynamical models and find that hydrodynamic correlations, depending on their range, can either enhance or dampen collective particle motions. These two distinct effects work synergistically to bias crystallite deformations toward a subset of the available outcomes.

42 ENGINEERING↗