Search NASA⌕ Search

SEARCH · Search NASA

Results for “Flow Simulator”

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 181 records · Page 10

Modeling and simulations of hydrodynamic shocks in a plasma flowing across randomized ICF scale laser beams

High-energy laser beams interacting with flowing plasmas can produce a plasma response that leads to deflection of the beam, beam bending. Such beams have usually a speckle structure generated by optical smoothing techniques that reduce the spatial and temporal coherence in the laser field pattern. The cumulative plasma response from laser speckles slows down the velocity of the incoming flow by momentum conservation. For slightly super-sonic flow the cumulative plasma response to the ponderomotive force exerted by the beam speckle ensemble is the strongest, such that slowing down the flow to subsonic velocities leads eventually to the generation of a shock around the cross section of the beam. This scenario has been predicted theoretically and is confirmed here by our hydrodynamic simulations in two dimensions with speckled beams and in one dimension with a reduced model. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed laser beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging from those produced by beams with only a random phase plates in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser plasma coupling in inertial confinement fusion (ICF) experiments are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A code-to-code benchmark for magneto-convection in a horizontal duct

Liquid metals and magnetic fields are used in many technical applications such as metallurgy, crystal growth and nuclear fusion reactors. When an electrically conducting fluid moves in a magnetic environment, electric currents and electromagnetic forces are generated that affect velocity and pressure losses in the flow. These magnetohydrodynamic (MHD) interactions have to be investigated to optimize the engineering processes. The characteristics of MHD flows depend on the geometrical configuration, the strength of the applied magnetic field, the electrical properties of fluid and structural materials and the thermal conditions. In the so-called blankets for fusion reactors, where liquid metals are used to breed the plasma fuel component tritium and to extract the generated heat, magneto-convective flows play a crucial role in determining heat and mass transfer. Therefore, the availability of numerical codes to simulate this type of flow is mandatory and their validation is a necessary step to guarantee the reliability of the results. For that reason, a benchmark problem has been defined to simulate liquid metal flows in a horizontal rectangular duct heated from below and exposed to a non-uniform magnetic field. Results obtained by five research groups using different codes are compared.

benchmark↗

Hydrodynamic characterization of the coastal pioneer array ocean observing system

Ocean observation buoys require relatively small amounts of power, yet traditionally necessitate costly resupply trips for battery replacement. With the offshore location of the buoys and small power requirements, wave energy may be an effective solution for providing consistent and reliable power to support the buoy instrumentation. The US National Science Foundation Ocean Observatories Initiative (OOI) includes arrays of point absorber-like buoy systems used for ocean observation that have been deployed at multiple locations including the Southern Mid-Atlantic Bight. A study is currently underway to design a pitch resonator wave energy converter to supplement existing renewable energy generation for powering observation instrumentation. This paper details field measurements from surface moorings of the OOI Coastal Pioneer Array, which informs the subsequent development of a numerical model for the moored observation system. The model is developed in Wave Energy Converter Simulator (WEC-Sim), which leverages the Simscape multibody solver within the MATLAB/Simulink framework and linear potential flow theory to simulate the hydrodynamic interactions and multibody dynamics in 6 degrees of freedom. Multiple tuning variables are considered to produce a model for the system that matches well with empirical data (about 8% error). In conclusion, the WEC-Sim model will serve as a platform for integrating the pitch resonator wave energy converter concept and deployment preparation (detailed design including power take-off and control systems, response evaluation, etc.).

hydrodynamic modeling↗

The Effect of Flow on CO2 Corrosion of Self-Healing Metallic Coatings

Internal corrosion is an issue that affects natural gas pipelines, a significant part of the United States’ energy infrastructure. Over time, this corrosion has worn away longstanding pipelines due to the original construction materials used for the lines and impurities in the gas and liquid streams flowing through them. The main impurities in natural gas are H2O, CO2, H2S, and O2. One of the leading causes of corrosion is the CO2 dissolved in water, giving rise to carbonic acid formation, which can further dissolve the steel pipe. The National Energy Technology Laboratory (NETL) has been studying different solutions to this problem. One potential answer is a self-healing sacrificial metallic coating applied by a novel cold spray technique. This study explored the impact of flow on corrosion with this coating and carbon steel when exposed to a saturated CO2 environment by simulating the pipeline flow profile in a small-scale lab setting. The samples were evaluated using multiple electrochemical techniques that found corrosion rates and were backed up with surface analysis to conclude the behavior of the corrosion mechanisms. The cold spray coating exhibited steel corrosion protection under flow conditions. This study has significantly advanced our understanding of CO2 corrosion, particularly in the context of natural gas pipelines and coating design.

carbon steel↗

The Effect of Flow on CO2 Corrosion of Self-Healing Metallic Coatings

Internal corrosion is an issue that affects natural gas pipelines, a significant part of the United States energy infrastructure. Over time, this corrosion has worn away longstanding pipelines due to the original construction materials used for the lines and impurities in the gas and liquid streams flowing through them. The main impurities in natural gas are H2O, CO2, H2S, and O2. One of the leading causes of corrosion is the CO2 dissolved in water, giving rise to carbonic acid formation, which can further dissolve the steel pipe. The National Energy Technology Laboratory (NETL) has been studying different solutions to this problem. One potential answer is a self-healing sacrificial metallic coating applied by a novel cold spray technique. Previous studies at NETL have shown that the zinc-rich self-healing coating can withstand corrosion through passivation and create a barrier to the diffusion of corrosive species. This study explored the impact of flow on corrosion with this coating and carbon steel when exposed to a saturated CO2 environment by simulating the pipeline flow profile in a small-scale lab setting. The samples were evaluated using multiple electrochemical techniques that found corrosion rates and were backed up with surface analysis to conclude the behavior of the corrosion mechanisms.

CO2 corrosion↗

Optimization and Experimental Validation of Annular Finned PCM-HX for a Domestic Hot Water Heater Application

The load profile for domestic water heating is time-dependent and can result in high energy demand during peak operating times. Shifting this peak load can have significant environmental and economic impacts. Phase change material (PCM)-based thermal energy storage (TES) is a potentially useful technology for peak load shifting in domestic hot water (DHW) applications thanks to its high latent heat and energy density. In this study, an annular finned-tube PCM-HX design concept was optimized for a load-shifting TES unit to meet the Department of Energy standard for a medium-usage DHW heater using a resistance-capacitance model (RCM) integrated with a Multi-Objective Genetic Algorithm. The optimized design comprised 70 identical annular finned-tube PCM-HX units connected in parallel and utilizing RT62HC as the PCM. A single PCM-HX unit was prototyped and tested in a vertically oriented setup with upward heat transfer fluid (HTF) flow. The hot water supply time was defined based on a cutoff temperature of 51.7°C. The as-designed mass flow rate (1.5 g/s) was tested to assess the performance of the prototyped PCM-HX unit for RCM validation. For the experimental investigation, RTD sensor bundles measured HTF temperature at the PCM-HX inlet and outlet, and a Coriolis flow meter accurately measured the HTF mass flow rate. The simulated discharging power underpredicted the experimental result by about 12%, and the simulated hot water supply time underpredicted the experimental result by approximately 13% for the as-designed mass flow rate (1.5 g/s). The average deviation of the hot water supply temperature between the experimental and RCM results during the complete PCM solidification process was 1.3 K for the as-designed mass flow rate. The overall good agreement between the experimental and RCM results provides confidence that computationally efficient models such as RCM can be utilized for design optimization of PCM-HXs.

42 ENGINEERING↗

Discrete-Element and Material-Point Method (DEM and MPM) Based Solvers for Sustainable Technologies

We present the use of discrete element method (DEM) and material point method (MPM) in three relevant green technology applications that include biomass feedstock handling, lithium-ion battery manufacturing, and high-pressure reverse osmosis. Our open-source DEM and MPM solvers are developed using performance portable grid and particle management library, AMReX, thus enabling superior performance on NVIDIA and AMD GPUs with > 100 million particles. Our DEM solver resolves the motion of individual particles in a granular system and includes a bonded sphere method for modeling non-spherical particles along with Hertzian and liquid bridge-based contact models. We simulate highly variable biomass feedstock flows in large-scale hoppers for biofuel production and electrode calendering in battery manufacturing using DEM. Our simulations predict flow blockage in large scale biomass hoppers and electrode microstructure variations, thus providing valuable information for biofuel and battery manufacturers, respectively. The second half of the talk will be on MPM and its application towards pore resolved simulations of reverse osmosis membranes under compressive loads. We present a validation study of our MPM simulations with membrane microscopy imaging thus providing useful insights on membrane stability under high pressure conditions. We also present a spectral stability analysis of using linear hat, quadratic and cubic spline basis in MPM indicating regions of numerical stability.

BIOMASS FUELS,MATHEMATICS AND COMPUTING↗

Liquid crystal torons in Poiseuille-like flows

Three-dimensional (3D) simulations of the structure of liquid crystal (LC) torons, topologically protected distortions of the LC director field, under material flows are rare but essential in microfluidic applications. Here, we show that torons adopt a steady-state configuration at low flow velocity before disintegrating at higher velocities, in line with experimental results. Furthermore, we show that under partial slip conditions at the boundaries, the flow induces a reversible elongation of the torons, also consistent with the experimental observations. These results are in contrast with previous simulation results for 2D skyrmions under similar flow conditions, highlighting the need for a 3D description of this LC soliton in relation to its coupling to the material flow. These findings pave the way for future studies of other topological solitons, like hopfions and heliknotons, in flowing soft matter systems.

Science & Technology - Other Topics↗

Flow-Driven Stretch Fluctuations Govern the Nonlinear Viscoelasticity of Elongating Associative Polymer Networks

We use nonequilibrium molecular dynamics simulations to verify recent tube-model predictions that associative polymer networks exhibit broad stretch fluctuations during elongational flow. Simulations further show that these fluctuating dynamics give rise to the rate-dependent extensional viscosity 𝜂 𝐸 measured in filament stretching experiments on H-bonding networks. Simulations model bivalent associative networks with a reactive bead-spring model for varying association strength and extensional strain rate. We observe that stretch fluctuations are driven by a new form of chain tumbling, where chains continually collapse and elongate as their associations break and reform within the advecting network. This produces a broad, nearly uniform distribution of chain stretch over a wide range of strain rates, manifesting as a rate-independent plateau in the extensional stress. Our results show that the nonlinear viscoelasticity of associative networks is dominated by large fluctuations in molecular response, which cannot be captured by current mean-field models.

biomimetic & bio-inspired materials↗

A modern concept of Lagrangian hydrodynamics

Here, we offer a modern interpretation of Lagrangian hydrodynamics as employed in Lagrangian simulations of compressible fluid flow. Our main result is to show that artificial viscosity, traditionally viewed as a numerical artifice to control unphysical oscillations in flows with shocks, actually represents a physical process and is necessary to derive accurate simulations in any compressible flow. We begin by reviewing the origins of two numerical devices, artificial viscosity and finite-volume methods. We proceed to construct a mathematical (PDE) model that incorporates those numerics and in which a new length scale, the observer, arises representing the discretization. Associated with that length scale, there are new inviscid fluxes that are the artificial viscosity as first formulated by Richtmyer and an artificial heat flux postulated by Noh but typically not included in Lagrangian codes. We discuss the connection of our results to bivelocity hydrodynamics. We conclude with some speculation as to the direction of future developments in multidimensional Lagrangian codes as computers get faster and have larger memories.

97 MATHEMATICS AND COMPUTING↗

Electron Heating in the Transrelativistic Perpendicular Shocks of Tilted Accretion Flows

Abstract General relativistic magnetohydrodynamic (GRMHD) simulations of black hole tilted disks—where the angular momentum of the accretion flow at large distances is misaligned with respect to the black hole spin—commonly display standing shocks within a few to tens of gravitational radii from the black hole. In GRMHD simulations of geometrically thick, optically thin accretion flows, applicable to low-luminosity sources like Sgr A* and M87*, the shocks have transrelativistic speed, moderate plasma beta (the ratio of ion thermal pressure to magnetic pressure is β pi1 ∼ 1–8), and low sonic Mach number (the ratio of shock speed to sound speed is M s ∼ 1–6). We study such shocks with 2D particle-in-cell simulations, and we quantify the efficiency and mechanisms of electron heating for the special case of preshock magnetic fields perpendicular to the shock direction of propagation. We find that the postshock electron temperature T e2 exceeds the adiabatic expectation T e2,ad by an amount T e 2 / T e 2 , ad − 1 ≃ 0.0016 M s 3.6 , nearly independent of the plasma beta and of the preshock electron-to-ion temperature ratio T e1 / T i1 , which we vary from 0.1 to unity. We investigate the heating physics for M s ∼ 5–6 and find that electron superadiabatic heating is governed by magnetic pumping at T e1 / T i1 = 1, whereas heating by B -parallel electric fields (i.e., parallel to the local magnetic field) dominates at T e1 / T i1 = 0.1. Our results provide physically motivated subgrid prescriptions for electron heating at the collisionless shocks seen in GRMHD simulations of black hole accretion flows.

Astronomy & Astrophysics↗

Modeling a divertor with mid-leg pumping for high-power H-mode scenarios in DIII-D considering E x B drift flows

Edge-plasma simulations of a baffled, long-legged divertor in DIII-D, performed using the multi-fluid code UEDGE, indicate that the position of the detachment front is constrained to the location of the pump duct along the low-field side (LFS) baffle. Simulations including magnetic and E x B drifts were performed for 12.5 MW deuterium plasmas including intrinsic carbon and seeded neon to assess the optimal location of the LFS divertor pump to create a stable detachment front between the target and the X-point. The radiation front position in the simulations, taken to be indicative of the detachment front, can be controlled between the pump and X-point in the favorable magnetic field direction for H-mode access by moving the pump duct location upstream of the target along the LFS baffle. In the unfavorable magnetic field direction, the radial E θ x B drift flows are directed towards the pumping surface, efficiently removing the injected deuterium gas and limiting the sensitivity of the radiation front location to the gas injection rate. The role of pumping rate and drift direction on the pumping efficiency are also found to affect the divertor plasma conditions and detachment front location in UEDGE simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High temporal frequency data from a four turbine, blade-resolved wind farm simulation with ExaWind

The data was generated with ExaWind (https://github.com/Exawind) which couples AMR-Wind (https://github.com/Exawind/amr-wind/), Nalu-Wind (https://github.com/Exawind/nalu-wind), TIOGA (https://github.com/Exawind/tioga), and OpenFAST (https://github.com/OpenFAST/openfast). This is a large-scale simulation of a blade-resolved wind farm using the ExaWind software stack. ExaWind couples together a background flow solver, AMR-Wind, and a near-body solver, Nalu-Wind, through an overset technique from the TIOGA application. Another application, OpenFAST, handles the structural dynamics of the turbine blades and towers, which informs the fluid-structure interaction of the wind turbines with the flow solvers. This particular simulation includes four blade-resolved wind turbines operating in a turbulent atmospheric boundary layer. The AMR-Wind solver uses 500 million cells and is being solved on 256 AMD GPUs of the Oakridge Leadership Computing Facility Frontier supercomputer. Each turbine is assigned its own Nalu-Wind solver with over 13 million elements per turbine and solved using 448 CPU cores, for a total of 1792 CPU cores. For each node, 56 cores contain Nalu-Wind, while 8 cores correspond to AMR-Wind operations on the GPUs. Consequently, ExaWind is entirely utilizing the CPUs and the GPUs of the nodes concurrently. The data used in the visualization is full flow field data output from the simulation. It is lossy-compressed to a specific accuracy using ZFP and written to disk every 16 time-steps to enable real-time flow visualization. The flow fields are sampled at a high temporal frequency to enable real-time, 24fps visualization. The flow fields are sampled every 12 simulation time steps (every 0.04132s).

17 WIND ENERGY↗

Posterior comparison of model dynamics in several hybrid turbulence model forms

Hybrid turbulence models that can accurately reproduce unsteady three-dimensional flow physics across the entire range of grid scales and turbulence dynamics from Reynolds-averaged Navier–Stokes (RANS), through large-eddy simulation (LES), down to direct numerical simulations (DNS) are of increasing interest to the turbulence modeling community. However, despite decades of research and development, the basic tasks of eliminating poor-performing hybrid RANS-LES models and accelerating adoption of superior models through well-designed validation and verification have yet to occur. As a step in this direction, in this work we evaluate thirteen different hybrid RANS-LES models via systematic grid refinement of decaying homogeneous isotropic turbulence. We further derive a novel mathematical framework for assessing the energy partitioning dynamics of each Hybrid RANS-LES model, wherein model-to-model variations in energy partitioning can be interpreted as different feedback mechanisms operating on a low-dimensional nonlinear dynamical system. We found that model forms similar to the flow simulation methodology—also often termed very-large eddy simulation—are dynamically inconsistent with DNS at all resolutions. Additionally, we found a strong dynamical similarity in the feedback mechanisms of all models related to detached eddy simulation and partially averaged Navier–Stokes that is inherent to their general model forms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Online randomized interpolative decomposition with a posteriori error estimator for temporal PDE data reduction

Traditional low-rank approximation is a powerful tool for compressing large data matrices that arise in simulations of partial differential equations (PDEs), but suffers from high computational cost and requires several passes over the PDE data. The compressed data may also lack interpretability thus making it difficult to identify feature patterns from the original data. Here, to address these issues, we present an online randomized algorithm to compute the interpolative decomposition (ID) of large-scale data matrices in situ. Compared to previous randomized IDs that used the QR decomposition to determine the column basis, we adopt a streaming ridge leverage score-based column subset selection algorithm that dynamically selects proper basis columns from the data and thus avoids an extra pass over the data to compute the coefficient matrix of the ID. In particular, we adopt a single-pass error estimator based on the non-adaptive Hutch++ algorithm to provide real-time error approximation for determining the best coefficients. As a result, our approach only needs a single pass over the original data and thus is suitable for large and high-dimensional matrices stored outside of core memory or generated in PDE simulations. A strategy to improve the accuracy of the reconstructed data gradient, when desired, within the ID framework is also presented. We provide numerical experiments on turbulent channel flow and ignition simulations, and on the NSTX Gas Puff Image dataset, comparing our algorithm with the offline ID algorithm to demonstrate its utility in real-world applications.

Column subset selection↗

Thermal and Well Flow Performance of Closed-Loop Geothermal in the Wattenberg Area

Closed-loop geothermal systems provide an alternative to resource-constrained hydrothermal systems and stimulation-intensive enhanced geothermal systems. In this work, we apply the slender-body theory (SBT) model, to simulate the well flow and heat transfer performance of U-loop well designs in the Wattenberg area of the Denver-Julesburg Basin. Three U-loop well patterns are investigated, including a single-, double-, and multi-lateral design. The subsurface within the area of interest is characterized by deep, hot (> 200 degrees C) igneous/metamorphic basement rock underlying multiple sedimentary formations. The lateral section(s) of the U-loop lie(s) within a target depth of 6 km, where temperatures are estimated to approach 300 degrees C. As a base case, conduction-only heat transfer is investigated through simulations with the SBT model within U-loops with open-hole laterals that exchange heat directly with the hot, dry rock using water as a working fluid. The utilization of supercritical CO2 as a heat transfer fluid is also considered. For each scenario, the system performance in terms of annual heat production and temperature profile over a 20-year project lifetime are assessed. Also, the levelized costs of heat and electricity (LCOH and LCOE) are determined using a top-down techno-economic analysis model. The results show that the performance- and cost-optimized U-loop design is one having an injection-production well spacing of 1,000 meters with ten 50-meter-spaced laterals that traverse a subsurface system with a temperature gradient of 60 degrees C/km. By injecting 20 degrees C-water at a rate of 60 kg/s through this loop, an average heat production of 19 MWth (i.e., 2.2 MWe net plant output) can be achieved, resulting in an LCOE and LCOH of $136/MWhe and $1.53/GJ, respectively, over a 20-year project life.

closed-loop geothermal↗

Thermal and Well Flow Performance of Closed-Loop Geothermal in the Wattenberg Area: Preprint

Closed-loop geothermal systems provide an alternative to resource-constrained hydrothermal systems and stimulation-intensive enhanced geothermal systems. In this work, we apply the slender-body theory (SBT) model, to simulate the well flow and heat transfer performance of U-loop well designs drilled in the Wattenberg area of the Denver-Julesburg Basin. Three U-loop well patterns are investigated including a single, double, and multi-lateral design. The subsurface within area is characterized by deep, hot (> 200degreesC) igneous/metamorphic basement rock underlying multiple sedimentary formations. The lateral section(s) of the U-loop lies within a target depth of 6 km where temperatures are estimated to approach 300degreesC. As a base case, conduction-only heat transfer is investigated through simulations with the SBT model within U-loops with open-hole laterals that exchange heat directly with the hot dry rock using water as a working fluid. The utilization of supercritical CO2 as a heat transfer fluid is also considered. For each scenario, the system performance in terms of annual heat production and temperature profile over a 20-year project lifetime are assessed. Also, the levelized costs of heat and electricity (LCOH and LCOE) are determined using a top-down technoeconomic analysis model. The results show that the performance and cost optimized U-loop design is one having an injection-production well spacing of 1,000 meters with ten 50-meter-spaced laterals that traverse a subsurface system with a temperature gradient of 60degreesC/km. By injecting 20 degreesC-water at a rate of 60 kg/s through this loop, an average heat production of 19 MWth can be achieved, resulting in an LCOE and LCOH of $136/MWh and $1.53/GJ, respectively, over a 20-year project life.

closed-loop geothermal↗

Probabilistic Forward Modeling of Galaxy Catalogs with Normalizing Flows

Abstract Evaluating the accuracy and calibration of the redshift posteriors produced by photometric redshift (photo- z ) estimators is vital for enabling precision cosmology and extragalactic astrophysics with modern wide-field photometric surveys. Evaluating photo- z posteriors on a per-galaxy basis is difficult, however, as real galaxies have a true redshift but not a true redshift posterior. We introduce PZFlow, a Python package for the probabilistic forward modeling of galaxy catalogs with normalizing flows. For catalogs simulated with PZFlow, there is a natural notion of “true” redshift posteriors that can be used for photo- z validation. We use PZFlow to simulate a photometric galaxy catalog where each galaxy has a redshift, noisy photometry, shape information, and a true redshift posterior. We also demonstrate the use of an ensemble of normalizing flows for photo- z estimation. We discuss how PZFlow will be used to validate the photo- z estimation pipeline of the Dark Energy Science Collaboration, and the wider applicability of PZFlow for statistical modeling of any tabular data.

Astronomy & Astrophysics↗