Search NASA⌕ Search

SEARCH · Search NASA

Results for “two fluid model”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Coupling Thermal-Hydraulics with Reactor Pressure Vessel Fracture Models

Because of the importance of maintaining the integrity of the reactor pressure vessel (RPV) in light-water reactors, significant effort has been devoted over the last several decades to understand the potential failure mechanisms, characterize how material properties evolve due to environmental conditions, and develop tools to assess the probability of fracture in aged RPVs during potential transient events. One important aspect of this problem that has not yet been fully addressed is the effects of nonuniform cooling on the RPV inner wall due to the way low-temperature water injected into the primary coolant system is distributed within the RPV during potential loss-of-coolant accident transients. Colder temperatures are expected in the regions near the inlets, causing a “cold-plume” effect. Rapid cooling of the RPV wall increases the likelihood of fracture initiation, and localized cooling due to cold-plume effects could potentially exacerbate this problem. The Grizzly code can perform a three-dimensional simulation of the RPV thermomechanical response as well as the probabilistic assessment of fracture for a population of pre-existing flaws, which is necessary to account for spatial variations in coolant temperature. This report documents a first effort to couple Grizzly RPV fracture models with computational fluid dynamics models performed using the Cardinal code, which allows for a realistic representation of the temperature distribution. This report demonstrates this coupling approach on two test cases: a simplified proof-of-concept scenario and a realistic small-break loss-of-coolant accident scenario. In both cases, the effects of the spatially varying temperature are evident; although, because the small-break loss-of-coolant accident scenario was far from challenging the RPV’s margins, it is still inconclusive how important cold-plume effects are in realistic accident scenarios. However, this effort represents an important step toward performing realistic analyses to determine to what extent spatially nonuniform cooling might affect RPV integrity under more aggressive accident scenarios.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Layer pseudospin magnetism in a transition metal dichalcogenide double-moiré system

Spontaneous order of layer pseudospins in two-dimensional bilayers is common in quantum Hall systems, where it is responsible for hysteretic responses to gate fields in states with Ising order and giant drag voltages in states with XY (spontaneous interlayer phase coherence) order. In this article we predict that layer pseudospin order will also occur in double-moiré strongly correlated two-dimensional electron systems. Here, we comment on similarities and differences in the competition between the two types of order in quantum Hall and double-moiré systems, and relate our findings to previous work on Falicov-Kimball models of electronic ferroelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Additive Manufacturing Interface Profiles for Jet Studies on the Vertical Shock Tube

Understanding the mixing between two fluids of disparate densities (i.e., variable density mixing) is important for the lab's mission. Experiments are required to validate our codes, such as RAGE, FLAG, and the BHR mix model. The role initial conditions (ICs) have on instability growth and the transition to turbulence is not completely understood and comparisons with experiments are needed. The Extreme Fluids team focuses on experiments at Los Alamos that utilize gases and low shock velocities (Mach <2 in air). These experiments have the advantage of high fidelity diagnostics capable of temporally and spatially resolved velocity and density measurements. One main difficulty is setting up a specific interface profile between the two gases. While attempts have been made to do so without the aid of a solid structure (i.e., a membrane), the results are limited in the range of possible shapes and reproducibility. Membranes provide the best way to provide a deterministic initial condition. However, fragments from the membrane can interfere with the flow and instability growth, as well as the diagnostics. In this project, we focused on exploring the materials and techniques available via modern additive manufacturing (i.e., 3D printing) to test if they could address some of the shortcomings of prior membrane designs. Our goal was to find a material that would have minimal impact on our experiment while also providing us with a means to produce complex initial conditions. This report details our work to explore the range of additive manufacturing capabilities and the jet-like structures we were able to generate using three-dimensional initial conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Interpreting Mass and Radius Measurements of Neutron Stars with Dark Matter Halos

Abstract The high densities of neutron stars (NSs) could provide astrophysical locations for dark matter (DM) to accumulate. Depending on the DM model, these DM admixed NSs (DANSs) could have significantly different properties than pure baryonic NSs, accessible through X-ray observations of rotation-powered pulsars. We adopt the two-fluid formalism in general relativity to numerically simulate stable configurations of DANSs, assuming a fermionic equation of state (EOS) for the DM with repulsive self-interaction. The distribution of DM in the DANS as a halo affects the path of X-rays emitted from hot spots on the visible baryonic surface, causing notable changes in the pulse profile observed by telescopes such as NICER, compared to pure baryonic NSs. We explore how various DM models affect the DM mass distribution, leading to different types of dark halos. We quantify the deviation in observed X-ray flux from stars with each of these halos. We identify the pitfalls in interpreting mass and radius measurements of NSs inferred from electromagnetic radiation and constraining the baryonic matter EOS if these dark halos exist.

Shawqi, Shafayat (ORCID:0000000210956183)↗
DOE OSTI2020

Coupled Electro-Thermal Analysis of Permanent Magnet Synchronous Motor for Electric Vehicles

Permanent magnet synchronous motors (PMSM) are extensively used in electric vehicles. However, high internal heat generation and inefficient heat dissipation often limit the operational reliability, and longevity of the PMSM. Therefore, proper quantification of heat generation in electric motor and advanced embedded motor cooling techniques remain topics of immense interest. In order to accurately predict electro-magnetic performance, i.e., efficiency, component-wise heat losses and the corresponding temperature distribution of a jacket cooled machine, this paper presents a two-way iteratively coupled electro-thermal modeling framework. Finite element based software Motor-CAD has been utilized for electro-magnetic performance calculation of BMW i3 PMSM. Finite volume based computational fluid dynamics/heat transfer (CFD/HT) software ANSYS® FLUENT® has been employed to simulate the temperature distribution of the PMSM, using the electro-magnetic losses as heat input. Computed heat losses, stator, winding, rotor, and magnet temperatures are utilized as coupling parameters between the electro-magnetic and thermal models. A conventional one-way coupling algorithm has also been developed and compared to the newly proposed two-way coupling algorithm. Numerical results confirm that at high current density, one-way coupling algorithm significantly over-predicts the motor temperature compared to the two-way algorithm. A comprehensive analysis has been carried out to characterize the influences of current density, speed, and forced convection heat transfer coefficient on the heat losses, overall efficiency, and maximum temperature of the PMSM. Lastly, an efficiency map has been interpreted from the coupled electro-magnetic simulation.

42 ENGINEERING↗
DOE OSTI2024

Light-wave-controlled Haldane model in monolayer hexagonal boron nitride

In recent years, the stacking and twisting of atom-thin structures with matching crystal symmetry has provided a unique way to create new superlattice structures in which new properties emerge. In parallel, control over the temporal characteristics of strong light fields has allowed researchers to manipulate coherent electron transport in such atom-thin structures on sublaser-cycle timescales. Here we demonstrate a tailored light-wave-driven analogue to twisted layer stacking. Tailoring the spatial symmetry of the light waveform to that of the lattice of a hexagonal boron nitride monolayer and then twisting this waveform result in optical control of time-reversal symmetry breaking and the realization of the topological Haldane model in a laser-dressed two-dimensional insulating crystal. Further, the parameters of the effective Haldane-type Hamiltonian can be controlled by rotating the light waveform, thus enabling ultrafast switching between band structure configurations and allowing unprecedented control over the magnitude, location and curvature of the bandgap. This results in an asymmetric population between complementary quantum valleys that leads to a measurable valley Hall current, which can be detected by optical harmonic polarimetry. The universality and robustness of our scheme paves the way to valley-selective bandgap engineering on the fly and unlocks the possibility of creating few-femtosecond switches with quantum degrees of freedom.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Atwood effects on nonlocality of the scalar transport closure in Rayleigh-Taylor mixing

The importance of nonlocality is assessed in modeling mean scalar transport for turbulent Rayleigh-Taylor (RT) mixing at different Atwood numbers. Building on the two-dimensional incompressible work of Lavacot et al. [J. Fluid Mech. 985, A47 (2024)], the present work extends the macroscopic forcing method to variable density problems in three-dimensional space to measure moments of the generalized eddy diffusivity kernel in RT mixing for increasing Atwood numbers (𝐴 = 0.05, 0.3, 0.5, 0.8). It is found that as 𝐴 increases, (1) the eddy diffusivity moments become asymmetric and (2) the higher-order eddy diffusivity moments become larger relative to the leading-order diffusivity, indicating that nonlocality becomes more important at higher 𝐴. There is a particularly strong temporal nonlocality at higher 𝐴, suggesting stronger history effects. In conclusion, the implications of these findings for closure modeling for finite-Atwood RT are discussed.

general physics↗
DOE OSTI2023

Sucrose transport inside the phloem: Bridging hydrodynamics and geometric characteristics

In plants, the delivery of the products of photosynthesis is achieved through a hydraulic system labeled as phloem. This semi-permeable plant tissue consists of living cells that contract and expand in response to fluid pressure and flow velocity fluctuations. The Münch pressure flow theory, which is based on osmosis providing the necessary pressure gradient to drive the mass flow of carbohydrates, is currently the most accepted model for such sucrose transport. When this hypothesis is combined with the conservation of fluid mass and momentum as well as sucrose mass, many simplifications must be invoked to mathematically close the problem and to resolve the flow. Here this study revisits such osmotically driven flows by developing a new two-dimensional numerical model in cylindrical coordinates for an elastic membrane and a concentration-dependent viscosity. It is demonstrated that the interaction between the hydrodynamic and externally supplied geometrical characteristic of the phloem has a significant effect on the front speed of sucrose transport. These results offer a novel perspective about the evolutionary adaptation of plant hydraulic traits to optimize phloem soluble compounds transport efficiency.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Radial drift of plasma blobs in a toroidal magnetic field with fully kinetic and reduced fluid models

In curved magnetic geometries, field-aligned regions of enhanced plasma pressure and density, termed ‘blobs,’ move as coherent filaments across the magnetic field lines. Coherent blobs account for a significant fraction of transport at the edges of magnetic fusion experiments and arise in naturally-occurring space plasmas. This work examines the dynamics of blobs with a fully kinetic electromagnetic particle-in-cell code and with a drift-reduced fluid code. In low-beta regimes with moderate blob speeds, good agreement is found in the maximum blob velocity between the two simulation schemes and simple analytical estimates. The fully kinetic code demonstrates that blob speeds saturate near the initial sound speed, which is a regime outside the validity of the reduced fluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Energy transport analysis of NSTX plasmas with the TGLF turbulent and NEO neoclassical transport models

This work presents a study of plasma transport at low aspect ratio on the National Spherical Torus Experiment tokamak, where the turbulent and neoclassical energy fluxes calculated by the quasilinear Trapped Gyro Landau Fluid (TGLF) model and the multi species drift-kinetic Neoclassical solver (NEO) are validated against experimental data. The turbulent energy transport of two plasma discharges, one in the L-mode confinement regime and another in the H-mode regime, is dominated by electrostatic drift-wave instabilities, while the ion heat transport has a significant neoclassical contribution. The data analysis workflow is described in detail to understand how the variations of mapping and fitting of experimental data affect the power balance solution and subsequent flux-matching plasma profile predictions with the TGYRO solver. On average, the predicted plasma profiles are consistent with experimental data. However, the solutions are sensitive to various input parameters, including boundary conditions, and the electron-ion coupling. Linear gyrokinetic stability analysis demonstrates close agreement of the real frequencies of unstable modes between TGLF and CGYRO gyrokinetic simulations, but higher growth rates are predicted by TGLF, especially for the H-mode case. Estimates of the low-k modes' contributions to the total flux are consistent with linear stability analysis and the E × B suppression of turbulence in TGLF simulations with the SAT1 saturation model, while the SAT2 saturation model over-predicts the low-k modes' contribution in the H-mode case. Moreover, the results with SAT1 model are consistent with power balance analysis, which indicates only neoclassical ion energy fluxes inside ρ < 0.4 in the L-mode case and $\rho \unicode{x2A7D} 0.7$ in the H-mode case. The presence of multi-scale turbulence and ion-scale driven zonal flow mixing effects are also observed in TGLF scans of the electron turbulent heat flux over a range of temperature gradients and the electron-ion temperature ratio, which could explain the strong model sensitivity to variations of input parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A quantitative comparison between velocity dependent SIDM cross-sections constrained by the gravothermal and isothermal models

One necessary step for probing the nature of self-interacting dark matter (SIDM) particles with astrophysical observations is to pin down any possible velocity dependence in the SIDM cross-section. Major challenges for achieving this goal include eliminating, or mitigating, the impact of the baryonic components and tidal effects within the dark matter halos of interest – the effects of these processes can be highly degenerate with those of dark matter self-interactions at small scales. In this work, we select 9 isolated galaxies and brightest cluster galaxies (BCGs) with baryonic components small enough such that the baryonic gravitational potentials do not significantly influence the halo gravothermal evolution processes. We then constrain the parameters of Rutherford and Møller scattering cross-section models with the measured rotation curves and stellar kinematics through the gravothermal fluid formalism and isothermal method. Cross-sections constrained by the two methods are consistent at $1\sigma$ confidence level, but the isothermal method prefers cross-sections greater than the gravothermal approach constraints by a factor of ~ 3.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Chimera F-Series Gravitational Wave Emission Sourced from Matter and Neutrino Anisotropy

This dataset contains gravitational wave data sourced from the the time-dependent fluid quadrupole motion as well as neutrino anisotropy, in the Chimera F-Series two-dimensional core collapse supernova simulations. Data from two models initiated from two different progenitors are presented: F15.78 and F15.79. Please see the README for more information about the data structure and progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Radial drift of plasma blobs in a toroidal magnetic field with fully kinetic and reduced fluid models

Abstract In curved magnetic geometries, field-aligned regions of enhanced plasma pressure and density, termed ‘blobs,’ move as coherent filaments across the magnetic field lines. Coherent blobs account for a significant fraction of transport at the edges of magnetic fusion experiments and arise in naturally-occurring space plasmas. This work examines the dynamics of blobs with a fully kinetic electromagnetic particle-in-cell code and with a drift-reduced fluid code. In low-beta regimes with moderate blob speeds, good agreement is found in the maximum blob velocity between the two simulation schemes and simple analytical estimates. The fully kinetic code demonstrates that blob speeds saturate near the initial sound speed, which is a regime outside the validity of the reduced fluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A Comparison of Interface Growth Models Applied to Rayleigh–Taylor and Richtmyer–Meshkov Instabilities

Sophisticated numerical models that contain fluid interfaces rely upon interface evolution models to approximate the transition to turbulence near interfaces, in the presence of Rayleigh–Taylor (RTI) or Richtmyer–Meshkov (RMI) instability. Semi-analytical models have been developed in recent decades to predict the interface growth from an initial state into the nonlinear regime. Two of these models are considered in this study. They are the Goncharov and the z-models. Both of these interface models have strengths and weaknesses, which are examined here. Both of them have been implemented in the xRAGE compressible flow solver, which models fluid interfaces. The flow solver provides the fluids acceleration as a body force to the interface model. The purpose of such interface model is to evolve the early times of interface position as a subgrid model within a compressible flow simulation in order to then initialize a turbulence model. In this work, the interface models are assessed and compared for their evolution of RTI and RMI. The z-model performed better than the Goncharov model for all cases that were explored.

42 ENGINEERING↗
DOE OSTI2024

Thermal Fluid Modeling Approaches in SAM for High Temperature Gas-Cooled Reactor Applications

This paper provides an overview of the unique modeling and simulation challenges and needs in the system and safety analysis of two common types of High Temperature Gas-cooled Reactor (HTGR) designs. The challenges are associated with the complex core geometry configurations and the change of the dominant heat removal mechanisms between normal operating conditions and decay heat removal transients. This requires that simulation tools utilize models with reasonable length scales and capture both the localized heat transfer and the core-scale heat transfer at the same time. To meet such analysis needs, different methodologies were developed. In conclusion, this includes adopting existing methods and further developing/improving them, as well as proposing innovative methods that leverage advanced computational frameworks such as the Multiphysics Object Oriented Simulation Environment (MOOSE).

HTGR↗
DOE OSTI2022

A level set approach for the computational study of a yield stress fluid filling a thin mold

Many important engineering and scientific applications such as cement slurries, foams, crude oil, and granular avalanches involve the concept of yield stress. Therefore, modeling yield stress fluids in different flow configurations, including the accurate prediction of the yield surface, is important. In this paper, we present a computational model based on the finite element method to study the flow of yield stress fluids in a thin mold and compare the results with data from flow visualization experiments. Here, we use the level set method to describe the interface between the filling fluid and air. We use polypropylene glycol as a model Newtonian fluid and Carbopol for the model yield stress fluid, as the Carbopol solution demonstrates yielding without thixotropy. To describe the yielding and shear-thinning behavior, we use a generalized Newtonian constitutive equation with a Bingham–Carreau–Yasuda form. We compare the results obtained from the mold filling experiments with the results from the three-dimensional (3D) model and from a reduced-order Hele-Shaw (HS) model that is two-dimensional, including the effect of shear-thinning along the thin direction only approximately. We show that both the 3D and the HS model can capture the experimental meniscus shape reasonably well for all the fluids considered at three different flow rates. This indicates that the shape evolution is insensitive to the dimensionality of the model. However, the viscosity and yield surfaces predicted by the 3D and HS models are different. The HS model underestimates the high viscosity and unyielded regions compared to the estimation by the 3D model.

42 ENGINEERING↗
DOE OSTI2024

Relative permeabilities for two-phase flow through wellbore cement fractures

Multiple fluids are likely to exist in fractures and flow paths associated with leaky wellbores, including liquids (e.g., crude oil) and gases (e.g., gas exsolved from liquid). These fluids occupy and move through different portions of the pore spaces within the fractures depending on many factors, including fluid properties, fracture size, and the amount of the different fluids. Upward leakage of any phase, through the fracture, can contaminate water-bearing formations, create hazardous surface conditions, and compromise the functionality of the wellbore. Early signs of wellbore leaks may be expressed by anomalous pressure behavior at surface monitoring points on cavern storage wells. These pressure anomalies are difficult to interpret, necessitating knowledge of the factors that affect the multiphase flow in fractures and porous media. These parameters are critical to modeling multiphase flow in fractures. This insight can guide further diagnosis and maximize leak remediation. Here, our study focuses on the relationship of the liquid–gas relative permeabilities for representative variable-aperture wellbore cement fracture. To obtain the relative permeability of each phase, two-phase flow tests were conducted where both fluids were flowing simultaneously through a fractured wellbore cement specimen under a range of factors, namely (1) aperture size, (2) capillary numbers, and (3) viscosity ratio. The flow experiments were conducted under a range of confining stresses and flow velocities, using nitrogen gas and silicone oils (of different viscosities) in a specially designed pressure vessel. The sum of gas and oil relative permeabilities were found to be less than one under all conditions, which indicates that the presence of one phase affects the permeability of the other phase, and vice versa. Since the gas phase flow conditions include a significant inertial flow component in addition to viscous flow, the inertial flow coefficients at different saturation states are presented. The factors affecting the relationship between the relative permeabilities are discussed in detail. A new mathematical model for estimating the relative permeability of wellbore cement fracture is presented and experimentally validated.

58 GEOSCIENCES↗
DOE OSTI2020

3D CFD Model Validation Using Benchmark Data of 1/16th Scaled VHTR Upper Plenum and Development of Wall Heat-Transfer Correlation For Laminar Flow

With support from the U.S. Department of Energy-Office of Nuclear Energy’s (DOE-NE’s) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, an effort has been pursued to support high-temperature gas-cooled reactor (HTGR) technology development and its modeling and simulation needs. There is a particular need for advanced modeling and simulation tools to predict thermal-fluid behavior in the nuclear reactor primary system, especially in the core and the lower and upper plena, during safety-related transients. In this report, two main such activities are presented relevant to the HTGRs: (1) three-dimensional (3D) computational fluid dynamics (CFD) validation using benchmark data from the upper plenum of Texas A&M University’s 1/16th scaled very-high-temperature gas-cooled reactor (VHTR), and (2) development of wall heat-transfer correlation for laminar flow in a wall-heated pipe. The CFD tool validation exercises can be helpful to choose the models and CFD tools to simulate and design specific components of the HTRGs such as upper plenum where jet mixing is a complex phenomenon. In a loss of forced circulation event, the laminar flow can be observed during the development of natural circulation flow. This work includes the development and validation of heat transfer correlations for laminar flow using the Nek5000 CFD code due to limited available experimental data for laminar flow conditions to guide low-order models (1D). In this report, the flow characteristics of a single isothermal jet discharging into the upper plenum was investigated using the Nek5000 Large-Eddy Simulation (LES) CFD tool. Several numerical simulations were performed for various jet-discharged Reynolds numbers ranging from 3,413 to 12,819. A grid-independent study was performed. The numerical results of mean velocity, root-mean-square fluctuating velocity, and Reynolds stress were compared against the benchmark data. Good agreement was obtained between simulated and measured data for axial mean velocities, except near the upper plenum hemisphere. The maximum predicted errors for axial mean velocities at various normalized coolant channel diameter heights of 1, 5, and 10 are 1.56%, 1.88%, and 3.82%, respectively. In addition, the predicted root-mean-square fluctuating velocity and Reynolds stress are qualitatively in agreement with the experimental data. The Nek5000 code was used to develop wall-heat transfer correlation for laminar flow in a cylindrical tube. Several simulations were performed for various Reynolds flow and wall-heat fluxes. A new heat transfer correlation was developed using data from Nek5000 simulation results and regression functions in Matlab. The developed heat transfer correlation is valid for various Reynolds flows from 200 to 2000. The predicted R² value for model fit was 0.875, which ensures that 87.5% of the model data lies on the Nek5000 data. Moreover, a machine learning (ML) tool was used to train and test the Nek5000 data. A good fit of the ML-based model was observed with the test data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗