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,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Multibody for Everybody (M4E) - A Linearization Approach to Enable Frequency Domain Analysis, Time Integration and Control Co-Design

1.1 Background/Objectives: Marine energy represents a promising yet underexploited source of power. To increase the harvested power, significant efforts have been made to improve wave energy converter (WEC) modeling capabilities and optimize power take-off (PTO) performance; however, these efforts have often treated WEC dynamics, PTO design, and controller development sequentially. In contrast, control co-design (CCD) is emerging as a promising strategy to address these issues directly, creating a growing need for fast analysis tools suitable for repeated simulation and parametric studies [1]. To support this need, this work presents the Multibody for Everybody (M4E) [2] linearization module, which employs a symbolic toolbox to provide deeper insight of WEC design parameters. The objective is to demonstrate that a minimal-coordinate linearization of articulated WEC dynamics can provide accurate wave response predictions and substantial computational savings relative to nonlinear time-domain simulation, while preserving compatibility with broader wave-energy analysis workflows, enabling CCD. 1.2 Approach/Activities: The proposed approach linearizes the equations of motion, generated by M4E, in minimal coordinates about a selected operating point and combines the resulting system with frequencydomain hydrodynamic terms to incorporate the reduced mass, damping, stiffness, and forcing operators. The linearized model is used for both impedance-based response amplitude operator (RAO) prediction and rapid regular-wave time integration. The methodology is demonstrated on a single-flap device and a FOSWEC configuration, with linearized M4E responses compared against the corresponding nonlinear M4E simulations and WEC-Sim results. Regular-wave time histories, RAO trends, and runtime differences are assessed. The framework is also compatible with broader wave-energy workflows, including coupling to WecOptTool, although that capability is not the focus of this work [3]. 1.3 Results/Lessons: The linearized M4E model reproduces key regularwave response characteristics such as integration and Response Amplitude over multiple frequencies. This module matches nonlinear M4E and WEC-Sim results while substantially reducing integration cost. Thus, the proposed framework can serve as a rapid analysis layer for articulated WEC design, parameter studies, and controls-oriented workflows. The analysis is most appropriate in the near-equilibrium regime, about the linearization point.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

CFD Data on a Vertical Axis Wave Turbine

In this study January to July, 2023, different variations of the original geometry of a vertical-axis wave turbine (VAWT) were generated and evaluated for hydrodynamic power efficiency using computational fluid dynamics (CFD). The key geometrical parameters considered in this parametric study included the chord length of the rotor blades and the horizontal semi-axis length. The immersion depth of the rotor was also examined as a key deployment parameter for the wave turbine. The CFD simulation results revealed that a medium chord length of the blade (i.e., the same as that of the baseline design) and a shorter horizontal semi-axis for the guide curve of the blade than that of the baseline design resulted in higher hydrodynamic power to extract. With the most efficient turbine rotor geometry identified in this study, a deployment depth that could assure full submergence of the rotor in waves but as close to the free surface as possible led to a higher hydrodynamic power. These findings revealed a pathway for the improvement of the wave turbine energy efficiency.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Efficiently resolving the terrestrial-aquatic interface in E3SM with sub-grid methods to improve coastal simulations (Final Technical Report)

The broad goal of this project was to extend the capabilities of MPAS‐Ocean, the ocean component of DOE’s Energy Exascale Earth System Model (E3SM), to better resolve the hydrodynamics of the terrestrial‐aquatic interface without significantly increasing computing resource requirements. Specifically, our goal was to better resolve the small-scale features across the terrestrial‐aquatic interface which includes the low lying coastal floodplain and a dendritic network of hydraulically efficient channels connected to estuarine systems and back bays, barrier islands and high energy inlets. The processes driving flow across the terrestrial‐aquatic interface include upland and floodplain hydrology, tides, wind and atmospheric pressure (including coastally generated storm surge across continental shelves and within estuaries and bays), wind waves and water temperature and salinity gradients.

58 GEOSCIENCES↗
NASA NTRS1980

Dynamic simulation of coronal mass ejections

A model is developed for the formation and propagation through the lower corona of the loop-like coronal transients in which mass is ejected from near the solar surface to the outer corona. It is assumed that the initial state for the transient is a coronal streamer. The initial state for the streamer is a polytropic, hydrodynamic solution to the steady-state radial equation of motion coupled with a force-free dipole magnetic field. The numerical solution of the complete time-dependent equations then gradually approaches a stationary coronal streamer configuration. The streamer configuration becomes the initial state for the coronal transient. The streamer and transient simulations are performed completely independent of each other. The transient is created by a sudden increase in the pressure at the base of the closed-field region in the streamer configuration. Both coronal streamers and coronal transients are calculated for values of the plasma beta (the ratio of thermal to magnetic pressure) varying from 0.1 to 100.

Steinolfson, R. S.↗
DOE OSTI2026

Influence of the drive design on ablation front hydrodynamic instability growth in a capsule implosion at the National Ignition Facility

A critical aspect of inertial confinement fusion experiments lies in the control of ablation front instabilities during the implosion process. The growth of Rayleigh–Taylor and Richtmyer–Meshkov instabilities, seeded by target defects, can significantly degrade the performance of the implosion. Here, this study explores the influence of drive design on ablation front stability using the hydro-growth radiography platform at the National Ignition Facility. This platform allows the measurement of the ablation front hydrodynamic instability growth. Because these measurements are restricted to early convergence, their primary role is to constrain the initial instability growth and benchmark the simulations used to infer shell integrity and performance at peak velocity and ignition. Three ignition designs, Hybrid-E (HyE), SQ-n (“S” for scaling and “Q” for quality), and High temperature, High thickness (HiT) were analyzed. The results demonstrated a high dependence of the growth factor on the choice of design, capsule scale as well as the hohlraum conditions, mainly the radiative temperature and the gold M-band emission. HyE and SQ-n display a similar growth factor range, but their evolution is different due to the design differences in hohlraum conditions. HiT has the lowest growth factor of the studied designs, which is a result of a higher radiative temperature. These results highlight the importance of design choices in controlling instability dynamics. These insights inform future drive design strategies to enhance stability and efficiency in fusion ignition experiments.

Physics↗
DOE OSTI2023

Rheology of dense fiber suspensions: Origin of yield stress, shear thinning, and normal stress differences

We explain the origins of yield stress, shear thinning, and normal stress differences in rigid fiber suspensions. We investigate the interplay between the hydrodynamic, noncontact attractive and repulsive, and interfiber contact interactions. The shear-thinning viscosity and finite yield stress obtained from the Immerse Boundary Method simulations are in quantitative agreement with experiential results from the literature. In this study, we show that attractive interactions result in yield stress and shear thinning rheology in the suspensions of rigid fibers. This is an important finding, given the ongoing discussion regarding the origin of the yield stress for suspensions of fibers. Finally, the ability of the proposed model to quantitatively predict the rheology is not limited to only shear thinning and yield stress but also extends to normal stresses.

42 ENGINEERING↗
DOE OSTI2023

Rheology of dense fiber suspensions: Origin of yield stress, shear thinning, and normal stress differences

We use high-fidelity computational model of fiber suspension to investigate suspension rheology. We investigate the interplay between the hydrodynamic, noncontact attractive and repulsive, and interfiber contact interactions. The shear-thinning viscosity and finite yield stress obtained from the Immerse Boundary Method simulations align quantitatively with experimental findings from the literature. The study demonstrates that attractive interactions lead to both yield stress and shear thinning behavior in rigid fiber suspensions. This discovery holds significance as it contributes to the ongoing debate on the source of yield stress in fiber suspensions. The proposed model is used to quantify normal stresses in addition to shear thinning and yield stress.

Ardekani, Arezoo↗
NASA NTRS1972

Hydrodynamic water impact

The hydrodynamic impact of a falling body upon a viscous incompressible fluid was investigated by numerically solving the equations of motion. Initially the mathematical model simulated the axisymmetric impact of a rigid right circular cylinder upon the initially quiescent free surface of a fluid. A compressible air layer exists between the falling cylinder and the liquid free surface. The mathematical model was developed by applying the Navier-Stokes equations to the incompressible air layer and the incompressible fluid. Assuming the flow to be one dimensional within the air layer, the average velocity, pressure and density distributions were calculated. The liquid free surface was allowed to deform as the air pressure acting on it increases. For the liquid the normalized equations were expressed in two-dimensional cylindrical coordinates. The governing equations for the air layer and the liquid were expressed in finite difference form and solved numerically. For the liquid a modified version of the Marker-and-Cell method was used. The mathematical model has been reexamined and a new approach has recently been initiated. Essentially this consists of examining the impact of an inclined plate onto a quiesent water surface with the equations now formulated in cartesian coordinates.

Kettleborough, C. F.↗
NASA NTRS1991

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative, manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. These results are highlighted. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable. Another one of our goals is to examine the stability of our predicted flows. Because time-dependent three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulations). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, Robert W.↗
NASA NTRS1990

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program was to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, with the funding from this NASA program, we concentrated on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we developed unique global models that allowed us to study the coupling between the different regions. These results are highlighted in the next section. Another important aspect of our NASA Theory Program concerned the effect that localized 'structure' had on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkland current patterns) or time variations in these input due to storms and substorms. Also, some of the plasma flows that we predicted with our macroscopic models could be unstable, and another one of our goals was to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulations). Therefore, another goal of our NASA Theory Program was to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This could involve a detailed comparison of kinetic, semi-kinetic, and hydrodynamic predictions for a given polar wind scenario or it could involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations provides insight into when the various models can be used with confidence.

Schunk, Robert W.↗
NASA NTRS1992

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable, and another one of our goals is to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulation). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, R. W.↗
DOE OSTI2022

A Parameterized Neutrino Emission Model to Study Mass Ejection in Failed Core-collapse Supernovae

Some massive stars end their lives as failed core-collapse supernovae (CCSNe) and become black holes (BHs). Although in this class of phenomena the stalled supernova (SN) shock is not revived, the outer stellar envelope can still be partially ejected. This occurs because the hydrodynamic equilibrium of the star is disrupted by the gravitational mass loss of the protoneutron star (PNS) due to neutrino emission. We develop a simple parameterized model that emulates PNS evolution and its neutrino emission and use it to simulate failed CCSNe in spherical symmetry for a wide range of progenitor stars. Our model allows us to study mass ejection of failed CCSNe where the PNS collapses into a BH within ~100 ms and up to ~10 6 s. We perform failed CCSNe simulations for 262 different pre-SN progenitors and determine how the energy and mass of the ejecta depend on progenitor properties and the equation of state (EOS) of dense matter. In the case of a future failed CCSN observation, the trends obtained in our simulations can be used to place constraints on the pre-SN progenitor characteristics, the EOS, and on PNS properties at BH formation time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Baseline predictions of elliptic flow and fluctuations for the RHIC Beam Energy Scan using response coefficients

Currently the RHIC Beam Energy Scan is exploring a new region of the Quantum Chromodynamic phase diagram at large baryon densities that approaches nuclear astrophysics regimes. This provides an opportunity to study relativistic hydrodynamics in a regime where the net conserved charges of baryon number, strangeness, and electric charge play a role, which will significantly change the theoretical approach to simulating the baryon-dense Quark-Gluon Plasma. In this work, we detail many of the important changes needed to adapt both initial conditions and the medium to baryon-rich matter. Then, we make baseline predictions for the elliptical flow and fluctuations based on extrapolating the physics at LHC and top RHIC energies to support future analyses of where and how the new baryon-dense physics causes these extrapolations to break down. First we compare eccentricities across beam energies, exploring their underlying assumptions; we find the the extrapolated initial state is predicted to be nearly identical to that at AuAu $\sqrt{s_{NN}}$ = 200 GeV. Then the final flow harmonic predictions are based on linear+cubic response. We discuss preliminary STAR results in order to determine the implications that they have for linear+cubic response coefficients at the lowest beam energy of AuAu $\sqrt{s_{NN}}$ = 7 GeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Effects of dissipative baryon current in heavy-ion collisions at RHIC-BES energies

The CLVisc (3+1)D viscous hydrodynamic model is extended to include the equation of net baryon conservation and the Israel-Stewart-like equations for dissipative baryon current. Using the NEOSB equation of state, we simulate the dynamical evolution and collectivity of the quark-gluon plasma with finite chemical potential, assuming smooth energy density and net baryon density distributions at the initial proper time. Numerical results are shown for the impact of net-baryon dissipation on particle yields and p T spectra in heavy-ion collisions at beam energy scan energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Virtual Flow Solver - Geophysics: A 3D Incompressible Navier-Stokes Solver

Virtual Flow Solver - Geophysics (VFS-Geophysics) is a three-dimensional (3D) incompressible Navier-Stokes solver based on the Curvilinear Immersed Boundary (CURVIB) method. The CURVIB is a sharp interface type of immersed boundary (IB) method that enables the simulation of fluid flows in the presence of geometrically complex moving bodies. The CURVIB method can be applied to wind/MHK turbine simulations and energy applications. VFS-Geophysics is the result of many years of research work by several graduate students, post-docs, and research associates that have been involved in the Computational Hydrodynamics and Biofluids Laboratory directed by Professor Fotis Sotiropoulos. The preparation of the present manual has been supported by the U.S. Department of Energy (DE-EE 0005482).

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Validation of magnetized gas-jet experiments to investigate the effects of an external magnetic field on laser-plasma instabilities

Laser-plasma instabilities (LPI) play a detrimental role in energy coupling to the target in inertial confinement fusion (ICF). The recent development of applied strong magnetic fields for use in ICF and laboratory astrophysics experiments has opened opportunities to investigate the role of external magnetic fields on LPIs. Recent numerical studies have shown that stimulated Raman scattering (SRS) can be mitigated by external magnetic fields in the kinetic regime of the instability and warrant systematic experimental studies to validate modelling. To this end, we design experiments at the OMEGA-EP laser facility to investigate the effect of an external perpendicular B-field of 5–30 T on the backscattered light spectrum from a gas-jet target. We present measurements from a proof-of-principle experiment, where the backward-SRS (BSRS) is in the kinetic regime, for which the magnetization is expected to produce the greatest effects on instability growth. New simulations motivated by this experiment are used to inform the proposal of an upgraded experimental design. Our simulation predictions show that the new design is suited to experimentally demonstrating BSRS mitigation by an external magnetic field, despite the magnetization effects on the hydrodynamics, as well as the inherent temporal and spatial variations in plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Laser-heated capillary discharge plasma waveguides for electron acceleration to 8 GeV

A plasma channel created by the combination of a capillary discharge and inverse Bremsstrahlung laser heating enabled the generation of electron bunches with energy up to 7.8 GeV in a laser-driven plasma accelerator. The capillary discharge created an initial plasma channel and was used to tune the plasma temperature, which optimized laser heating. Although optimized colder initial plasma temperatures reduced the ionization degree, subsequent ionization from the heater pulse created a fully ionized plasma on-axis. The heater pulse duration was chosen to be longer than the hydrodynamic timescale of ≈ 1 ns, such that later temporal slices were more efficiently guided by the channel created by the front of the pulse. Here, simulations are presented which show that this thermal self-guiding of the heater pulse enabled channel formation over 20 cm. The post-heated channel had lower on-axis density and increased focusing strength compared to relying on the discharge alone, which allowed for guiding of relativistically intense laser pulses with a peak power of 0.85 PW and wakefield acceleration over 15 diffraction lengths. Electrons were injected into the wake in multiple buckets and times, leading to several electron bunches with different peak energies. To create single electron bunches with low energy spread, experiments using localized ionization injection inside a capillary discharge waveguide were performed. A single injected bunch with energy 1.6 GeV, charge 38 pC, divergence 1 mrad, and relative energy spread below 2% full-width half-maximum was produced in a 3.3 cm-long capillary discharge waveguide. This development shows promise for mitigation of energy spread and future high efficiency staged acceleration experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Coarsening dynamics of ternary polymer solutions with mobility and viscosity contrasts

Using phase-field simulations, we investigate the bulk coarsening dynamics of ternary polymer solutions undergoing a glass transition for two models of phase separation: diffusion only and with hydrodynamics. The glass transition is incorporated in both models by imposing mobility and viscosity contrasts between the polymer-rich and polymer-poor phases of the evolving microstructure. For microstructures composed of polymer-poor clusters in a polymer-rich matrix, the mobility and viscosity contrasts significantly hinder coarsening, effectively leading to structural arrest. For microstructures composed of polymer-rich clusters in a polymer-poor matrix, the mobility and viscosity contrasts do not impede domain growth; rather, they change the transient concentration of the polymer-rich phase, altering the shape of the discrete domains. This effect introduces several complexities to the coarsening process, including percolation inversion of the polymer-rich and polymer-poor phases—a phenomenon normally attributed to viscoelastic phase separation.

Chemistry↗