Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics modeling”

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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Simulations of the Mg II K and Ca II 8542 Lines From an Alfvén Wave-Heated Flare Chromosphere

We use radiation hydrodynamic simulations to examine two models of solar flare chromospheric heating: Alfven wave dissipation and electron beam collisional losses. Both mechanisms are capable of strong chromospheric heating, and we show that the distinctive atmospheric evolution in the mid-to-upper chromosphere results in Mg II k-line emission that should be observably different between wave-heated and beam-heated simulations. We also present Ca II 8542 A profiles that are formed slightly deeper in the chromosphere. The Mg II k-line profiles from our wave-heated simulation are quite different from those from a beam-heated model and are more consistent with Interface Region Imaging Spectrograph observations. The predicted differences between the Ca II 8542 A in the two models are small. We conclude that careful observational and theoretical study of lines formed in the mid-to-upper chromosphere holds genuine promise for distinguishing between competing models for chromospheric heating inflares.

Kerr, Graham S.↗
NASA NTRS2020

Tidal Disruptions of Main Sequence Stars - I. Observable Quantities and their Dependence on Stellar and Black Hole Mass

This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main-sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses(MM0.1510)and six black hole masses( MM105 105BH7). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of Må for Må/Me3. Independent of Må, relativistic effects increase the critical pericenter distance for full disruption events by up to a factor of∼3 relative to the Newtonian prediction. The probability of a direct capture is also independent of Må; at MBH/Me;5×106 this probability is equal to the probability of a complete disruption. The breadth of the debris energy distribution ΔE can differ from the standard estimate by factors of0.35−2, depending on Må and MBH, implying a corresponding change(∝(ΔE)−3/2)in the characteristic mass-return timescale. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly universal relation between the star’s angular momentum and the fraction of Må remaining. Within the “empty loss-cone” regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions.

Taeho Ryu↗
NASA NTRS2020

Tidal Disruptions of Main Sequence Stars - II. Simulation Methodology and Stellar Mass Dependence of the Character of Full Tidal Disruptions

This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main-sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses (0.15 < M_/M<10)and six black hole masses( ́10<MM105 105BH7). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of Må for Må/Me3. Independent of Må, relativistic effects increase the critical pericenter distance for full disruption events by up to a factor of ∼3 relative to the Newtonian prediction. The probability of a direct capture is also independent of Må; at MBH/Me5×106this probability is equal to the probability of a complete disruption. The breadth of the debris energy distribution ΔE can differ from the standard estimate by factors of 0.35−2, depending on Må and MBH, implying a corresponding change(∝(ΔE)−3/2)in the characteristic mass-return timescale. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly universal relation between the star’s angular momentum and the fraction of Må remaining. Within the “empty loss-cone” regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions

Taeho Ryu↗
DOE OSTI2022

1.2.4.404 - Data-Driven Approach for Hydropower Plant Controller Prototyping Using Remote Hardware in the Loop (DR-HIL)

Real-time prototyping of hydropower plant controls is important for reducing the cost and the risk of field deployment. This project will 1) collect design and operational data from actual hydro plants and 2) use a physics-informed machine learning approach for real-time emulation of hydropower plants, including hydro turbine and hydrodynamics. The data-driven models will be interfaced with digital real-time simulation at NREL’s Flatirons campus for hardware-in-the-loop (HIL) testing of the governor hardware device or controller-HIL (CHIL). The proposed approach will also establish the connectivity based remote CHIL testing capability using real-time data streams from an actual hydro plant. This integrated hydro-plant emulation with CHIL will be used to prototype hydro-governor controls and eventually provide an opportunity to test hydropower integrated with various technologies (e.g. conventional and renewable generation, energy conversion, etc.) as HIL.

controls prototyping↗
DOE OSTI2023

Prediction of Breakdown in Air and Solid Dielectrics: A Complete Plasma Model from Discharge Initiation to Flashover (Final Technical Report)

A major challenge for most of the technologies of today, and the technologies of the future, is the compact and reliable design of electrical insulation systems for their electrical components. Although controlled low-temperature plasmas (LTPs) or non-thermal plasmas in air at atmospheric pressure are of interest for a wide range of technologies, the appearance of such discharges is considered highly undesirable in insulation systems because they are precursors of a breakdown or flashover, where complete insulation failure takes place. For non-thermal plasmas, often only one mechanism, or the transition from one discharge mechanism to another, has typically been studied to date. On the other hand, in the case of thermal plasma applications, the conditions following the occurrence of an electric arc or flashover have always been considered. Here, there is a big technical gap: there is no paper describing a complete plasma model simulating all discharge mechanisms and their transitions from discharge initiation to electrical arc and flashover. To address this gap, a complete plasma model for both thermal and non-thermal plasmas should be developed using a hydrodynamic approach. The project aimed to build the model. The work has the potential to address industrial applications and is of great importance to higher voltage systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Measurement and Modeling of Electron Temperature in Laboratory Photoionized Plasmas Relevant to Astrophysics

The electron temperature of photoionized plasmas characterizes the thermalization of photoelectrons, impacts the charge-state distribution, emissivity, and opacity through atomic recombination processes, and is needed to perform detailed comparisons with theory predictions. We discuss temperature measurements in laboratory photoionized plasmas and a comparison with model calculations done with several theory approximations and codes. These include a radiation-hydrodynamics simulation and two nonequilibrium heating and ionization models that tracked the evolution of the internal energy of the electrons. Furthermore, for the same physics model and X-ray flux time history, calculations were performed assuming steady-state or time-dependent conditions. The time history of steady-state results correlates with that of the X-ray flux, while that of the time-dependent cases does not, and it is qualitatively and quantitatively different from the steady-state case. Steady-state results significantly overestimated temperature measurements, while time-dependent results produced better approximations, which suggests the importance of transient effects in the experiment and also the need for time-resolved measurements.

Laboratory astrophysics↗
DOE OSTI2020

Experiments and simulations of isochorically heated warm dense carbon foam at the Texas Petawatt Laser

An experimental and simulation study of warm dense carbon foams at ambient density (n e ~10 21 cm –3 ) is presented. This study of isochorically heated foams is motivated by their potential application in carbon-atmosphere white-dwarf envelopes, where there are modeling uncertainties due to the equation of state. The foams are heated on an approximately picosecond time scale with a laser-accelerated proton beam. The cooling and expansion of the heated foams can be modeled with appropriately initialized radiation-hydrodynamics codes; xRAGE code is used in this work. The primary experimental diagnostic is the streaked optical pyrometer, which images a narrow band of radiation from the rear surface of the heated material. Presented are xRAGE modeling results for both solid aluminum targets and carbonized resorcinol-formaldehyde foam targets, showing that the foam appears to cool slowly on the pyrometer because of partial transparency. So that simulations of cooling foam are processed properly, it is necessary to account for finite optical depth in the photosphere calculation, and the methods for performing that calculation are presented in depth.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Learning robust parameter inference and density reconstruction in flyer plate impact experiments

Estimating physical parameters or material properties from experimental observations is a common objective in many areas of physics and material science. In many experiments, especially in shock physics, radiography is the primary means of observing the system of interest. However, radiography does not provide direct access to key state variables, such as density, which prevents the application of traditional parameter estimation approaches. Here we focus on flyer plate impact experiments on porous materials, and resolving the underlying parameterized equation of state (EoS) and crush porosity model parameters given radiographic observation(s). We use machine learning as a tool to demonstrate with high confidence that using only high impact velocity data does not provide sufficient information to accurately infer both EoS and crush model parameters, even with fully resolved density fields or a dynamic sequence of images. We thus propose an observable data set consisting of low and high impact velocity experiments/simulations that capture different regimes of compaction and shock propagation, and proceed to introduce a generative machine learning approach which produces a posterior distribution of physical parameters directly from radiographs. We demonstrate the effectiveness of the approach in estimating parameters from simulated flyer plate impact experiments, and show that the obtained estimates of EoS and crush model parameters can then be used in hydrodynamic simulations to obtain accurate and physically admissible density reconstructions. Finally, we examine the robustness of the approach to model mismatches, and find that the learned approach can provide useful parameter estimates in the presence of out-of-distribution radiographic noise and previously unseen physics, thereby promoting a potential breakthrough in estimating material properties from experimental radiographic images.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

T2DM v1.0

T2DM adds two-dimensional hydrodynamic dispersion to TOUGH2/EOS7C. T2DM can be used when modeling flow and transport of brine in porous media to account for hydrodynamic dispersion.

Oldenburg, CurtisM↗
DOE OSTI2022

Validation of ray-based cross-beam energy transfer models

Ray-based cross-beam energy transfer (CBET) models have become a common feature of the radiation-hydrodynamic codes used to simulate inertial confinement fusion experiments. Here, these models are necessary for achieving better agreement with experimental measurements, but their detailed implementation can vary widely between the codes and often rely on artificial multipliers. To address this, a series of 2D and 3D test cases has been developed with validated solutions from wave-based calculations. Comparisons of various ray-based CBET models to the wave-based calculations highlight the essential physics that is required for accurate ray-based CBET modeling. Quantitative comparison metrics and/or field data from the wave-based calculations have been made available for use in the validation of other ray-based CBET codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1954

Hydrodynamic Investigation of the the Take-Off Characteristics of a 1/10-Scale Dynamic Model of the Convair XF2Y-1 Airplane

An investigation was made of the take-off characteristics of a 1/10-scale dynamic model of the Convair XF2Y-1 airplane. This airplane is a water-based, jet-propelled, delta-wing fighter incorporating a hydro-ski landing gear. Tests were made with the original configuration, with the beaching wheels removed, and with the wheels installed and fairings added in front of the wheels. Each configuration was tested at weight and balance conditions simulating 17,000 pounds gross weight with the moment due t o 7,600 pounds of thrust, 17,300 pounds gross weight with a 9,500-pound thrust condition, and 23,000 pounds gross weight with a 9,300-pound thrust condition. Constant-speed runs were made at various elevon settings and vertical ski-strut positions; and trim, rise, and resistance were measured. Accelerated runs were made with controlled elevons and scale shock struts which could be extended as desired, and the longitudinal stability and spray characteristics were observed and photographed.

McBridge, Ellis E.↗
DOE OSTI2025

Accelerating computational fluid dynamics simulation of post-combustion carbon capture modeling with MeshGraphNets

Packed columns are commonly used in post-combustion processes to capture CO 2 emissions by providing enhanced contact area between a CO 2 -laden gas and CO 2 -absorbing solvent. To study and optimize solvent-based post-combustion carbon capture systems (CCSs), computational fluid dynamics (CFD) can be used to model the liquid–gas countercurrent flow hydrodynamics in these columns and derive key determinants of CO 2 -capture efficiency. However, the large design space of these systems hinders the application of CFD for design optimization due to its high computational cost. In contrast, data-driven modeling approaches can produce fast surrogates to study large-scale physics problems. We build our surrogates using MeshGraphNets (MGN), a graph neural network framework that efficiently learns and produces mesh-based simulations. We apply MGN to a random packed column modeled with over 160K graph nodes and a design space consisting of three key input parameters: solvent surface tension, inlet velocity, and contact angle. Our models can adapt to a wide range of these parameters and accurately predict the complex interactions within the system at rates over 1700 times faster than CFD, affirming its practicality in downstream design optimization tasks. This underscores the robustness and versatility of MGN in modeling complex fluid dynamics for large-scale CCS analyses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Advancement of hybrid fluid-kinetic modeling for HEDP and ICF science

We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Recent Developments in the WEC-Sim Open-Source Design Tool: Preprint

WEC-Sim (Wave Energy Converter SIMulator) is an open-source code for simulating wave energy converters, which has been actively developed and applied to simulate a wide variety of device archetypes and has become a popular tool since its initial release in 2014. WEC-Sim is developed jointly by the National Renewable Energy Laboratory (NREL) and Sandia National Laboratories (SNL) within the MATLAB/SIMULINK environment. Figure 1 illustrates a general wave-to-wire model which begins with a deployment site resource characterization, which is used to complete the hydrodynamic simulation of a single WEC (or array), with the power generation profile imported to a grid simulator to understand the influence on the local electrical network. While modelling the entire wave-to-wire is difficult and encompass multiple time scales and physics, WEC-Sim is focused on the hydrodynamics simulation to predict, analyze and optimize WEC dynamics and power performance. WEC-Sim simulations are performed in the time domain based on the radiation and diffraction method using hydrodynamics coefficients derived from boundary element method (BEM) based-frequency-domain potential flow solvers (e.g., WAMIT, NEMOH, Capytaine, or ANSYS-AQWA). Within this level of modeling fidelity, WEC-Sim can handle floating body hydrodynamics, mechanical and electrical power generation methods, advanced control implementation, mooring systems, and other unique applications such as desalination. Table 1 lists additional WEC-Sim functionalities, which are created using prebuilt Simulink blocks and MATLAB scripts that can simulate a wide range of floating systems and the corresponding auxiliary subsystems.

hydrodynamics modeling↗
NASA NTRS1994

Thermohydrodynamic analysis of cryogenic liquid turbulent flow fluid film bearings, phase 2

The Phase 2 (1994) Annual Progress Report presents two major report sections describing the thermal analysis of tilting- and flexure-pad hybrid bearings, and the unsteady flow and transient response of a point mass rotor supported on fluid film bearings. A literature review on the subject of two-phase flow in fluid film bearings and part of the proposed work for 1995 are also included. The programs delivered at the end of 1994 are named hydroflext and hydrotran. Both codes are fully compatible with the hydrosealt (1993) program. The new programs retain the same calculating options of hydrosealt plus the added bearing geometries, and unsteady flow and transient forced response. Refer to the hydroflext & hydrotran User's Manual and Tutorial for basic information on the analysis and instructions to run the programs. The Examples Handbook contains the test bearing cases along with comparisons with experimental data or published analytical values. The following major tasks were completed in 1994 (Phase 2): (1) extension of the thermohydrodynamic analysis and development of computer program hydroflext to model various bearing geometries, namely, tilting-pad hydrodynamic journal bearings, flexure-pad cylindrical bearings (hydrostatic and hydrodynamic), and cylindrical pad bearings with a simple elastic matrix (ideal foil bearings); (2) improved thermal model including radial heat transfer through the bearing stator; (3) calculation of the unsteady bulk-flow field in fluid film bearings and the transient response of a point mass rotor supported on bearings; and (4) a literature review on the subject of two-phase flows and homogeneous-mixture flows in thin-film geometries.

Sanandres, Luis↗
NASA NTRS1952

Hydrodynamic Impact of a System with a Single Elastic Mode II : Comparison of Experimental Force and Response with Theory

Hydrodynamic impact tests were made on an elastic model approximating a two-mass spring system to determine experimentally the effects of structural flexibility on the hydrodynamic loads encountered during seaplane landing impacts and to correlate the results with theory. A flexible seaplane was represented by a two-mass spring system consisting of a rigid prismatic float connected to a rigid upper mass by an elastic structure. The model had a ratio of sprung mass to hull mass of 0.6 and a natural frequency of 3.0 cycles per second. The tests were conducted in smooth water at fixed trims and included both high and low flight-path angles and a range of velocity. Theoretical and experimental comparisons indicated that the theoretical results agreed well with the experimental results.

Miller, Robert W↗
DOE OSTI2025

How Do Uncertainties in Galaxy Formation Physics Impact Field-level Galaxy Bias?

Our ability to extract cosmological information from galaxy surveys is limited by uncertainties in the galaxy–dark matter halo relationship for a given galaxy population, which are governed by the intricacies of galaxy formation. To quantify these uncertainties, we examine quenched and star-forming galaxies using two distinct approaches to modeling galaxy formation: U$\scriptsize{NIVERSE}$M$\scriptsize{ACHINE}$, an empirical semianalytic model, and the I$\scriptsize{llustris}$TNG hydrodynamical simulation. We apply a second-order hybrid N-body perturbative bias expansion to each galaxy sample, enabling direct comparison of modeling approaches and revealing how uncertainties in the galaxy–halo connection affect bias parameters and non-Poisson noise across number densities and redshifts. Notably, we find that quenched and star-forming galaxies occupy distinct parts of the bias parameter space, and that the scatter induced from these different galaxy formation models is small when conditioned on similar selections of galaxies. We also detect the signature of assembly bias in our samples; this leads to small but significant deviations from analytic bias predictions, while assembly bias-removed samples match these predictions well. This work indicates that galaxy samples from a spectrum of reasonable, physically motivated models for galaxy formation give a relatively small range of field-level galaxy bias parameters. We estimate a set of priors from these models that should be useful in extracting cosmological constraints from luminous red galaxy– and emission line galaxy–like samples. Looking forward, careful estimates of the range of impacts of galaxy formation, for a given sample and cosmological analysis, will be an essential ingredient for extracting the most precise cosmological information from current and future large galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗