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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Unifying disparate rate-dependent rheological regimes in non-Brownian suspensions

A typical dense non-Brownian particulate suspension exhibits shear thinning (decreasing viscosity) at low shear rate/stress followed by a Newtonian plateau (constant viscosity) at intermediate shear rate/stress values which transitions to shear thickening (increasing viscosity) beyond a critical shear rate/stress value and finally, undergoes a second shear thinning transition at extremely high shear rate/stress values. In this study, we unify and quantitatively reproduce all the disparate rate dependent regimes and the corresponding transitions for a dense non-Brownian suspension with increasing shear rate/stress. Here, we employ discrete particle dynamics simulations based on the proposed mechanism to elucidate its accuracy. We find that a competition between inter-particle interactions of hydrodynamic and non-hydrodynamic origins and the switching in the dominant stress scale with increasing the shear rate/stress lead to each of the above transitions. Inclusion of traditional hydrodynamic interactions, attractive/repulsive DLVO (Derjaguin and Landau, Verwey and Overbeek) interactions, the inter-particle contact interactions and a constant friction (or other constraint mechanism) reproduce the initial thinning as well as the shear thickening transition. However, to quantitatively capture the intermediate Newtonian plateau and the second shear thinning, an additional non-hydrodynamic interaction of non-DLVO origin and a decreasing coefficient of friction, respectively, are essential; thus, providing the first explanation for the presence of the intermediate Newtonian plateau along with reproducing the second shear thinning in a single framework. Expressions utilized for various interactions and friction are determined from experimental measurements and hence, result in an excellent quantitative agreement between the simulations and previous experiments.

42 ENGINEERING↗
DOE OSTI2022

Design of a high energy density experiment to measure the suppression of hydrodynamic instability in an applied magnetic field

A transverse magnetic field can suppress the hydrodynamic instability growth of an unstable plasma interface. This effect, of interest to inertial confinement fusion and astrophysics, has mostly been studied via simulation. Here, we present the design of an experiment at the National Ignition Facility (NIF) to demonstrate this effect in a laboratory. Simulations indicate that the timescale for the diffusion of the magnetic field across the mixing region should be at least comparable to the timescale of the instability growth in order to have a measurable suppression effect. This motivates the use of lower density target materials than usual high energy density (HED) hydrodynamics experiments to permit faster hydrodynamics and higher plasma conductivities (through higher temperature), for a given laser drive and magnetic field. We discuss a target design for creating a Rayleigh–Taylor unstable HED plasma interface that uses 320 mg/cc iodine-doped carbon foam as the heavy material, 20 mg/cc carbon foam as the light material, and a 6 lm amplitude, a 120 lm wavelength ripple machined at the interface, which shows a measurable suppression effect with a nominal NIF drive and a 30 T magnetic field (the present facility limit). Here, models indicate lower density foams that may display even larger suppression effects, as the Hohlraum drive also radiatively preheats the foam to permit even higher temperatures and, hence, higher conductivities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modelling a Heaving Point-Absorber with a Closed-Loop Control System Using the DualSPHysics Code

The present work addresses the need for an efficient, versatile, accurate and open-source numerical tool to be used during the design stage of wave energy converters (WECs). The device considered here is the heaving point-absorber developed and tested by Sandia National Laboratories. The smoothed particle hydrodynamics (SPH) method, as implemented in DualSPHysics, is proposed since its meshless approach presents some important advantages when simulating floating devices. The dynamics of the power take-off system are also modelled by coupling DualSPHysics with the multi-physics library Project Chrono. A satisfactory matching between experimental and numerical results is obtained for: (i) the heave response of the device when forced via its actuator; (ii) the vertical forces acting on the fixed device under regular waves and; (iii) the heave response of the WEC under the action of both regular waves and the actuator force. This proves the ability of the numerical approach proposed to simulate accurately the fluid–structure interaction along with the WEC’s closed-loop control system. In addition, radiation models built from the experimental and WAMIT results are compared with DualSPHysics by plotting the intrinsic impedance in the frequency domain, showing that the SPH method can be also employed for system identification.

42 ENGINEERING↗
DOE OSTI2022

Impacts of climate change on estuarine stratification and implications for hypoxia within a shallow subtropical system

Vertical density stratification often plays an important role in the formation and expansion of coastal hypoxic zones through its effect on near-bed circulation and vertical oxygen flux. However, the impact of future climate change on estuarine circulation is widely unknown. Here, we developed and calibrated a three-dimensional hydrodynamic model for Pensacola Bay, a shallow subtropical estuary in the northeastern Gulf of Mexico. Model simulations based on years 2013–2017 were applied to examine changes in salinity, temperature, and density under future climate scenarios, including increased radiative forcing (IR) and temperature (T), increased freshwater discharge (D), sea level rise (SLR), and wind intensification (W). Simulations showed that the impacts of climate change on modeled state variables varied over time with external forcing conditions. The model demonstrated the potential for sea level rise and increased freshwater discharge to episodically increase vertical density gradients in the Bay. However, increased wind forcing destabilized vertical gradients, at times reducing the spatial extent and duration of stable stratification. For time periods with low freshwater discharge, moderate increases in wind speed (10%) can destabilize density gradients strengthened by increased discharge (10%) and sea level rise (0.48 m). In contrast, destruction of strong density gradients that form near the mid-Bay channel following flood events requires stronger wind forcing. Furthermore, these results highlight the importance of considering natural variability in freshwater and wind forcing, as well as local phenomena that are generally unresolved by global climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Studying QGP with flow: A theory overview

I review recent developments in the phenomenological study of the quark-gluon plasma (QGP) transport properties based on a personal selection of results that were presented at Quark Matter 2019. Additionally, the constraints on the temperature dependence of QGP shear and bulk viscosity are summarized. I discuss new theory advancements towards more realistic 3D dynamical simulations of heavy-ion collisions at finite baryon density. The challenges and opportunities of applying hydrodynamics to small collision systems are highlighted.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Thermal fluctuations in the dissipation range of homogeneous isotropic turbulence

In this work, using fluctuating hydrodynamics we investigate the effect of thermal fluctuations in the dissipation range of homogeneous isotropic turbulence. Simulations confirm theoretical predictions that the energy spectrum is dominated by these fluctuations at length scales comparable to the Kolmogorov length. We also find that the extreme intermittency in the far-dissipation range predicted by Kraichnan is replaced by Gaussian thermal equipartition.

42 ENGINEERING↗
DOE OSTI2025

Shear-Induced Conformations of Salt-Free Polyelectrolytes in Semidilute Solutions

Here, we present coarse-grained molecular dynamics simulations of salt-free polyelectrolyte chains in semidilute solutions under simple shear flow, with full hydrodynamic interactions and explicit dipolar solvent. At equilibrium, chain orientation statistics follow a pseudo-Voigt distribution, and the structural correlation length and chain end-to-end vector autocorrelation function exhibit scaling behavior consistent with theoretical predictions for polyelectrolytes. Under shear, chains transition from coiled to stretched states and the end-to-end vector autocorrelation function reveals oscillatory dynamics at high Weissenberg numbers. Analysis of the gyration tensor and shear strain distributions identified three distinct chain populations with directional alignment along and against the shear gradient. Compared with their neutral polymer counterparts, polyelectrolytes exhibit stronger shear thinning and enhanced chain alignment under the same Weissenberg number, which is attributed to electrostatic interactions and shear-induced counterion release. These findings provide molecular insight into the distinct flow response of charged polymers and have implications for tailoring the rheological properties of polyelectrolyte-based materials.

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

Fluctuating hydrodynamics of chiral active fluids

Active materials are characterized by continuous injection of energy at the microscopic level and typically cannot be adequately described by equilibrium thermodynamics. Here we study a class of active fluids in which equilibrium-like properties emerge when fluctuating and activated degrees of freedom are statistically decoupled, such that their mutual information is negligible. We analyse three paradigmatic systems: chiral active fluids composed of spinning frictional particles that are free to translate, oscillating granular gases and active Brownian rollers. In all of these systems, a single effective temperature generated by activity parameterizes both the equation of state and the emergent Boltzmann statistics. The same effective temperature, renormalized by velocity correlations, relates viscosities to steady-state stress fluctuations via a Green–Kubo relation. To rationalize these observations, we develop a theory for the fluctuating hydrodynamics of these non-equilibrium fluids and validate it through large-scale molecular dynamics simulations. Our work sheds light on the microscopic origin of odd viscosities and stress fluctuations characteristic of parity-violating fluids, in which mirror symmetry and detailed balance are broken.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

DOE-ICoM/RIFT

Rapid Infrastructure Flood Tool (RIFT) is a two-dimensional hydrodynamic model based on the complete shallow water equations. RIFT has specifically been designed with rapid simulation in mind by utilizing commodity high performance computing technology and best-available nation-wide data. RIFT is used to predict the movement of water over land and resolve the spatial and temporal variability of flood depths, extent, and velocity. RIFT can be applied to many flood situations and has primarily been used to quantify flood extents from dam/levee failure or inland rainfall flooding.

Perkins, Bill [Pacific Northwest National Laborato↗
DOE OSTI2021

Machine Learning for Conservative-to-Primitive in Relativistic Hydrodynamics

The numerical solution of relativistic hydrodynamics equations in conservative form requires root-finding algorithms that invert the conservative-to-primitive variables map. These algorithms employ the equation of state of the fluid and can be computationally demanding for applications involving sophisticated microphysics models, such as those required to calculate accurate gravitational wave signals in numerical relativity simulations of binary neutron stars. This work explores the use of machine learning methods to speed up the recovery of primitives in relativistic hydrodynamics. Artificial neural networks are trained to replace either the interpolations of a tabulated equation of state or directly the conservative-to-primitive map. The application of these neural networks to simple benchmark problems shows that both approaches improve over traditional root finders with tabular equation-of-state and multi-dimensional interpolations. In particular, the neural networks for the conservative-to-primitive map accelerate the variable recovery by more than an order of magnitude over standard methods while maintaining accuracy. Neural networks are thus an interesting option to improve the speed and robustness of relativistic hydrodynamics algorithms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Error quantification and comparison of binary neutron star gravitational waveforms from numerical relativity codes

Future gravitational wave detections of merging binary neutron star systems have the possibility to tightly constrain the equation of state of dense nuclear matter. In order to extract such constraints, gravitational waveform models need to be calibrated to accurate numerical relativity simulations of the late inspiral and merger. In this work, we take an essential step toward classifying the error and potential systematics in current generation numerical relativity simulations of merging binary neutron stars. To this end, we perform a direct comparison of two codes (fil, spec), which differ in many aspects, including the numerical methods and discretizations used and equations solved. We find that despite these different approaches, the codes are—within current numerical resolution bounds—fully consistent, and broadly comparable in cost for a given accuracy level. Our results indicate that the error in the waveforms is primarily dominated by the hydrodynamic evolution, consistent with earlier findings in the literature. Finally, we also discuss current limitations and cost estimates for numerical relativity simulations to reach the accuracies required in the era of next-generation gravitational detectors.

gravitational waves↗
DOE OSTI2023

SPH modeling of biomass granular flow: Theoretical implementation and experimental validation

The commercialization of biomass-derived energy is impeded by flowability challenges arising from the feeding and handling of granular biomass materials in full-scale biorefineries. To overcome these obstacles, a robust and accurate model to simulate the flow of granular biomass is indispensable. However, conventional mesh-based numerical codes are limited by inherent mesh distortion in simulating large deformation that commonly occurs in granular biomass handling. Here, in this study, we propose a graphics processing unit (GPU)-accelerated meshless Smoothed Particle Hydrodynamics (SPH) code to model the flow of granular biomass materials. A modified void ratio-based mass conversation, a hybrid particle-to-particle/surface frictional boundary treatment, and a hypoplastic constitutive model are implemented. Four numerical examples, an elastic block sliding on inclined planes, sand column collapse, Angle of Repose, and axial compression tests for pine chips, were simulated using the developed SPH code. The results demonstrate good agreement between numerical predictions and analytical and experimental data for all four examples, validating the SPH code and increasing confidence that it can be applied to simulate more complex granular biomass handling processes, such as hopper feeding or auger conveyance.

09 BIOMASS FUELS↗
DOE OSTI2022

Can Cooling and Heating Functions Be Modeled with Homogeneous Radiation Fields?

Abstract Cooling and heating functions describe how radiative processes impact the thermal state of a gas as a function of its temperature and other physical properties. In a most general case the functions depend on the detailed distributions of ionic species and on the radiation spectrum. Hence, these functions may vary on a very wide range of spatial and temporal scales. In this paper, we explore cooling and heating functions between 5 ≤ z ≤ 10 in simulated galaxies from the Cosmic Reionization On Computers project. We compare three functions: (1) the actual cooling and heating rates of hydrodynamic cells as a function of cell temperature, (2) the median cooling and heating functions computed using median interstellar medium (ISM) properties (median ISM), and (3) the median of the cooling and heating functions of all gas cells (instantaneous). We find that the median ISM and instantaneous approaches to finding a median cooling and heating function give identical results within the spread due to cell-to-cell variation. However, the actual cooling (heating) rates experienced by the gas at different temperatures in the simulations do not correspond to either summarized cooling (heating) functions. In other words, the thermodynamics of the gas in the simulations cannot be described by a single set of a cooling plus a heating function with a spatially constant radiation field that could be computed with common tools, such as CLOUDY.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Control Design for a Marine Hydrokinetic Cycloturbine Vehicle

Abstract Marine Hydrokinetic (MHK) cycloturbines generate sustainable power by exploiting tidal currents. By powering the turbines and using pitching foils for control, a vehicle comprised of MHK cycloturbines also has the ability to station keep and maneuver. The vehicle consists of four counter-rotating cycloturbines, with hydrofoils oriented perpendicular to the flow in a paddlewheel configuration. Lift and drag generated from these foils sum together to produce thrust. An experimentally tuned simulation model that solves the six-degrees-of-freedom rigid body equations of motion for the MHK vehicle subject to hydrodynamic, hydrostatic, and propulsive forces is used to aid the design of vehicle controllers. Global feedback controllers are initially designed by applying classical control methods to an approximate linear model of the system dynamics. A higher performing nonlinear controller is designed using the nonlinear dynamic inversion (NDI) method. NDI accounts for the nonlinearities of the MHK system and therefore is suitable for a wide range of operating conditions. The response of the classical and NDI controllers to speed, depth, roll, pitch, and yaw commands are evaluated and compared in simulation. The classical controller outperforms the NDI controller for small amplitude maneuvers, although the degradation with NDI is minor. However, in the nonlinear operating regime the NDI controller outperforms the classical controller and the classical controller exhibits instability.

Automation & Control Systems↗
DOE OSTI2021

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Future Response of Coastal Wetlands to Environmental Stresses: Algorithm Comparison of Numerical Models

Coastal wetlands are a critical component of the earth system that strongly influence the global water and biogeochemical cycles. They are also likely important sentinel of climate change. Because eco-geomorphological processes have long been recognized to be important for coastal wetland survival under accelerated sea-level rise (SLR), many eco-geomorphology models have been developed to assess the impact of climate change on coastal wetlands. Although these models differ substantially in complexity and numerical methods, few studies have investigated the algorithm-level uncertainties in these models. In this study, we developed a multiple-algorithm model framework of coastal wetlands that represents coastal hydrodynamics (such as water level, significant wave height and bottom shear stress) and four eco-geomorphological processes: mineral accretion, organic matter accretion, storm surge erosion and landward migration. We validated the model at three representative coastal wetland sites (Venice Lagoon, Plum Island Estuary and Hunter Estuary) for hydrodynamics, mineral accretion and organic matter accretion. Through model-data comparison, we showed that the model can well capture the dynamics of hydrodynamical and eco-geomorphological conditions in the study sites. Importantly, analysis of the multiple-algorithm simulations suggests that differences in the process representation of mineral and organic matter accretion may contribute to the recent contradicting predictions of coastal wetland evolution under accelerated SLR.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Spatial Mapping of Riverbed Grain-Size Distribution Using Machine Learning

Recent alluvial sediments in riverbeds play a significant role in controlling hydrologic exchange flows (HEFs) in river systems. The alluvial layer is usually associated with strong heterogeneity in physical properties (e.g., permeability and hydraulic conductivity), which affects local HEFs and therefore biogeochemical processes. The spatial distribution of these physical properties needs to be determined to inform the numerical models used to reveal the realistic hydro-biogeochemical behaviors. Such information can be obtained based on the intrinsic link between sediment grain-size distribution and hydraulic properties where sediment texture information is available. However, grain-size measurements are usually spatially sparse and do not have adequate coverage and resolution, particularly for a relatively large domain such as the Hanford Reach of the Columbia River. In this paper, we adopted machine learning (ML) approaches for categorizing and mapping the spatial distributions of riverbed substrate grain size and filling in missing areas of substrate data using the ML models along the reach. Such ML models for substrate size mapping were trained at 13,372 locations using measured substrate sizes along with observed and simulated attributes, including bathymetric attributes (e.g., elevation, slope, and aspect ratio) from LIDAR and bathymetric surveys, and hydrodynamic properties (e.g., water depth, velocity, shear stress, and their statistical moments). An ensemble bagging-based ML technique, Random Forest, was adopted to identify the most influential factors as predictors to develop the predictive models with over-fitting issues addressed. The models were evaluated with respect to each individual substrate size class and the lumped group, and then used to generate the final substrate size maps covering all the grid cells in the numerical modeling domain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Investigating coastal backwater effects and flooding in the coastal zone using a global river transport model on an unstructured mesh

Abstract. Coastal backwater effects are caused by the downstream water level increase as a result of elevated sea level, high river discharge and their compounding influence. Such effects have crucial impacts on floods in densely populated regions but have not been well represented in large-scale river models used in Earth system models (ESMs), partly due to model mesh deficiency and oversimplifications of river hydrodynamics. Using two mid-Atlantic river basins as a testbed, we perform the first attempt to simulate the backwater effects comprehensively over a coastal region using the MOSART river transport model under an ESM framework, i.e., Energy Exascale Earth System Model (E3SM) configured on a regionally refined unstructured mesh, with a focus on understanding the backwater drivers and their long-term variations. By including sea level variations at the river downstream boundary, the model performance in capturing backwaters is greatly improved. We also propose a new flood event selection scheme to facilitate the decomposition of backwater drivers into different components. Our results show that while storm surge is a key driver, the influence of extreme discharge cannot be neglected, particularly when the river drains to a narrow river-like estuary. Compound flooding, while not necessarily increasing the flood peaks, exacerbates the flood risk by extending the duration of multiple coastal and fluvial processes. Furthermore, our simulations and analysis highlight the increasing strength of backwater effects due to sea level rise and more frequent storm surge during 1990–2019. Thus, backwaters need to be properly represented in ESMs to improve the predictive understanding of coastal flooding.

54 ENVIRONMENTAL SCIENCES↗