Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Modeling of injected current stream-induced 3D perturbations in local helicity injection plasmas

Solenoid-free tokamak startup techniques are essential for spherical tokamaks and offer a pathway to cost reduction and design simplification in fusion energy systems. Local helicity injection (LHI) is one such approach, employing compact edge current sources to drive open field line current that initiates and sustains tokamak plasmas. The recently commissioned Pegasus-III Experiment provides a platform for advancing this and other solenoid-free startup methods. This study investigates the effect of LHI on magnetic topology in Pegasus-III plasmas. A helical filament model represents the injected current, and the linear plasma response to its three-dimensional field is calculated with M3D-C1. Poincaré mapping reveals substantial flux surface degradation in all modeled cases. The onset of overlapping magnetic structures and large-scale surface deformation begins near Ψ N ≈ 0.37, indicating a broad region of perturbed topology extending toward the edge. In rotating plasmas, both single-fluid and two-fluid models exhibit partial screening of the n = 1 perturbation, with two-fluid calculations showing stronger suppression near the edge. In contrast, the absence of rotation leads to strong resonant field amplification in the single-fluid case, while the two-fluid case with zero electron rotation mitigates this amplification and preserves edge screening. Magnetic probe measurements indicate that modeling the stream with spatial spreading—representing distributed current and/or oscillatory motion—better reproduces measured magnetic power profiles than a rigid filament model. The results underscore the role of rotation and two-fluid physics in screening stream perturbations and point to plasma flow measurements and refined stream models as key steps toward improving predictive fidelity.

Schaefer, Carolyn E. [Univ. of Wisconsin, Madison,↗
DOE OSTI2021

Comprehensive framework for data-driven model form discovery of the closure laws in thermal-hydraulics codes

The two-phase two-fluid model is a basis of many thermal-hydraulics codes used in design, licensing, and safety considerations of nuclear power plants. Thermal-hydraulics codes rely on the closure laws to close the system of conservation equations and describe the interactions between phases. These laws, derived from years of experimental investigations, are semi-empirical correlations that lack generality and have a limited range of applicability. Increase of computational power, availability of new experiments, and development of high-fidelity simulations has increased the number of validation data. The discrepancies between the code predictions and the validation data are a great source of knowledge. Missing physics that are not included in the model but are important for the considered phenomena can be discovered by propagating the information from the experimental results through the model. Furthermore, physics-discovered data-driven model form (P3DM) methodology integrates available integral effect tests and separate effects tests to determine the necessary corrections to the model form of the closure laws. In contrast to existing calibration techniques, the methodology modifies the functional form of the closure laws. Based on the functional form of the correction, the missing physics that were not included in the original model can be discovered. The methodology provides the alternative to the machine learning approach, in which the model is discovered in the form of the intractable black-box relation. In this work, the methodology was applied to the CTF subchannel code to improve the prediction of the two-phase flow phenomena.

42 ENGINEERING↗
DOE OSTI2020

Non-local model for surface tension in fluid-fluid simulations

We propose a non-local model for surface tension obtained in the form of an integral of a molecular-force-like function with support 3.5e added to the Navier-Stokes momentum conservation equation. We demonstrate analytically and numerically that with the non-local model interfaces with a radius of curvature larger than the support length behave macroscopically and microscopically, otherwise. For static droplets, the pressure difference Pe,in - Pe,out satisfies the Young-Laplace law for droplet radius greater than 3.5e and other- wise deviates from the Young-Laplace law. The latter indicates that the surface tension in the proposed model decreases with decreasing radius of curvature, which agrees with molecular dynamics and experimental studies of nanodroplets. Using the non-local model we perform numerical simulations of droplets under dynamic conditions, including a rising droplet, a droplet in shear flow, and two colliding droplets in shear flow, and compare results with a standard Navier-Stokes model subject to the Young-Laplace boundary condition at the fluid-fluid interface implemented via the Conservative Level Set (CLS) method. We find good agreement with existing numerical methods and analytical results for a rising macroscopic droplet and a droplet in a shear flow. For colliding droplets in shear flow, the non- local model converges (with respect to the grid size) to the correct behavior, including sliding, coalescence, and merging and breaking of two droplets depending on the capillary number. In contrast, we find that the results of the CLS model are highly grid-size dependent.

Howard, Amanda A.↗
DOE OSTI2020

HYDRODYNAMIC ANALYSIS OF MACROALGAE LOCAL MODEL USING COMPUTATIONAL FLUID DYNAMICS

In this article, a local scale, fully nonlinear coupled fluid-structural interaction (FSI) sugar kelp model has been developed using a computational fluid dynamics (CFD) method. In this model, to be consistent with available experimental data, the sugar kelp is approximated as elongated rectangles with smoothed isosceles triangles at the ends and a single kelp model with one end fixed in a channel with constant current model is developed. Several different current speeds are simulated, and the resulting drag forces and calculated drag coefficients are validated by comparison with experimental data from the literature. In a previous study, a global scale model was developed using a computational structural dynamics (CSD) method to simulate macroalgae farming system and guide the system configuration design. In the global scale model, the hydrodynamic forces are calculated using Morison’s equation and the kinematics and dynamics of the sugar kelp are simplified and the group of kelps attached to the long line is modeled as a slender structure with the same length and an effective diameter such that the volumes are consistent with the real physical system. This simplified model matches the weight and buoyancy but adjusting the hydrodynamic properties when the general hydrodynamic coefficients are employed. Therefore, optimal hydrodynamic coefficients used in global scale model were determined to obtain the hydrodynamic force more accurately. The validated local scale model is then be applied to determine the hydrodynamic coefficients of the simplified sugar kelp model for global dynamic analysis.

hydrodynamic modeling, macroalgae, CFD model↗
DOE OSTI2025

Speed of Sound Measurements of Binary Mixtures of trans -1,2-Difluoroethylene (R-1132(E)) with Difluoromethane (R-32) or 2,3,3,3-Tetrafluoropropene (R-1234yf)

The speed of sound of two binary mixtures containing 0.336 mole fraction trans-1,2- difluoroethylene (R-1132(E)) with difluoromethane (R-32) and 0.435 mole fraction R-1132(E) with 2,3,3,3-tetrafluoropropene (R-1234yf) were measured with a dual-path pulse-echo instrument. The speed of sound was measured along several pseudo-isochores for each blend at temperatures ranging from 230 K to 325 K for blends with R-32 and to 342 K for blends with R- 1234yf. Measurements started at pressures just above each mixtures bubble point pressure and were limited to 8 MPa to avoid potential disproportionation reactions of the R-1132(E). The data were compared to multi-fluid models incorporating Helmholtz-energy-explicit equations of state (EOS) for each pure fluid.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1989

The generalized drift flux approach: Identification of the void-drift closure law

The main characteristics and the potential advantages of generalized drift flux models are presented. In particular it is stressed that the issue on the propagation properties and on the mathematical nature (hyperbolic or not) of the model and the problem of closure are easier to tackle than in two fluid models. The problem of identifying the differential void-drift closure law inherent to generalized drift flux models is then addressed. Such a void-drift closure, based on wave properties, is proposed for bubbly flows. It involves a drift relaxation time which is of the order of 0.25 s. It is observed that, although wave properties provide essential closure validity tests, they do not represent an easily usable source of quantitative information on the closure laws.

Boure, J. A.↗
NASA NTRS2018

Thermal Modeling and Correlation of the Space Environments Complex Vacuum Chamber and Cryoshroud

A thermal model of the Space Environments Complex (SEC) vacuum chamber and cryoshroud has been developed in support of upcoming thermal vacuum/thermal balance testing for Orion EM-1. The model was developed in Thermal Desktop and includes the vacuum chamber itself, a fluid model of the gaseous nitrogen flowing through the cryoshroud and the chamber piping, the Heat Flux System (HFS) within the cryoshroud, and the mechanical ground support equipment (MGSE) that interfaces with the vehicle. It has been correlated with steady state data from three tests. Two tests ran the cryoshroud in hot mode at 170 °F, and one test ran the cryoshroud in cold mode at -263 °F. Correlation was done using an optimization algorithm to find values of unknown contact conductances that minimized the RMS error between the model predictions and the test data. Overall, model quality was very good with a total RMS error of 2.6 °F.

Stalcup, Erik J.↗
NASA NTRS2017

Fusion-Enabled Pluto Orbiter and Lander

The Pluto orbiter mission proposed here is credible and exciting. The benefits to this and all outer-planet and interstellar-probe missions are difficult to overstate. The enabling technology, Direct Fusion Drive, is a unique fusion engine concept based on the Princeton Field-Reversed Configuration (PFRC) fusion reactor under development at the Princeton Plasma Physics Laboratory. The truly game-changing levels of thrust and power in a modestly sized package could integrate with our current launch infrastructure while radically expanding the science capability of these missions. During this Phase I effort, we made great strides in modeling the engine efficiency, thrust, and specific impulse and analyzing feasible trajectories. Based on 2D fluid modeling of the fusion reactors outer stratum, its scrape-off-layer (SOL), we estimate achieving 2.5 to 5 N of thrust for each megawatt of fusion power, reaching a specific impulse, Isp, of about 10,000 s. Supporting this model are particle-in-cell calculations of energy transfer from the fusion products to the SOL electrons. Subsequently, this energy is transferred to the ions as they expand through the magnetic nozzle and beyond. Our point solution for the Pluto mission now delivers 1000 kg of payload to Pluto orbit in 3.75 years using 7.5 N constant thrust. This could potentially be achieved with a single 1 MW engine. The departure spiral from Earth orbit and insertion spiral to Pluto orbit require only a small portion of the total delta-V. Departing from low Earth orbit reduces mission cost while increasing available mission mass. The payload includes a lander, which utilizes a standard green propellant engine for the landing sequence. The lander has about 4 square meters of solar panels mounted on a gimbal that allows it to track the orbiter, which beams 30 to 50 kW of power using a 1080 nm laser. Optical communication provides dramatically high data rates back to Earth. Our mass modeling investigations revealed that if current high-temperature superconductors are utilized at liquid nitrogen temperatures, they drive the mass of the engine, partly because of the shielding required to maintain their critical temperature. Second generation materials are thinner but the superconductor is a very thin layer deposited on a substrate with additional layers of metallic classing. Tremendous research is being performed on a variety of these superconducting materials, and new irradiation data is now available. This raises the possibility of operating nearfuture high-temperature superconductors at a moderately low temperature to dramatically reduce the amount of shielding required. At the same time, a first generation space engine may require low-temperature superconductors, which are higher TRL and have been designed for space coils before (AMS-02 experiment for the ISS). We performed detailed analysis of the startup system and thermal conversion system components. The ideal working fluid was determined to be a blend of Helium and Xenon. No significant problems were identified with these subsystems. For the RF system, we conceived of a new, more efficient design using state-of-the-art switch amplifiers, which have the potential for 100% efficiency. This report presents details of our engine and trajectory analyses, mass modeling efforts, and updated vehicle designs.

temperature↗
DOE OSTI2020

Probabilistic neural networks for fluid flow surrogate modeling and data recovery

We consider the use of probabilistic neural networks for fluid flow surrogate modeling and data recovery. This framework is constructed by assuming that the target variables are sampled from a Gaussian distribution conditioned on the inputs. Consequently, the overall formulation sets up a procedure to predict the hyperparameters of this distribution which are then used to compute an objective function given training data. We demonstrate that this framework has the ability to provide for prediction confidence intervals based on the assumption of a probabilistic posterior, given an appropriate model architecture and adequate training data. The applicability of the present framework to cases with noisy measurements and limited observations is also assessed. To demonstrate the capabilities of this framework, we consider canonical regression problems of fluid dynamics from the viewpoint of reduced-order modeling and spatial data recovery for four canonical data sets. The examples considered in this study arise from (i) the shallow-water equations, (ii) a two-dimensional cylinder flow, (iii) the wake of a NACA0012 airfoil with a Gurney flap, and (iv) the NOAA sea surface temperature data set. Furthermore, the present results indicate that the probabilistic neural network not only produces a machine-learning-based fluid flow surrogate model but also systematically quantifies the uncertainty therein to assist with model interpretability.

42 ENGINEERING↗
DOE OSTI2021

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1973

Magnetospheric current sheets

A review is presented of theories and observations of current sheets in the magnetosphere, emphasizing the magnetopause and tail current sheets. Theoretical models of the magnetopause current sheet can be divided into two types: specular reflection and fluid models containing no external magnetic field, and models with an interplanetary/magnetosheath magnetic field. Our understanding of the first type is much better than of the second, although magnetospheric observations indicate that the external magnetic field plays an important role. Most of the theoretical models of the tail current sheet attempt to arrive at an understanding of the reconnection process or start with the assumption that reconnection is important. Observations are reviewed and implications for the models are discussed.

Speiser, T. W.↗
NASA NTRS1994

The evolution of cosmic-ray-mediated magnetohydrodynamic shocks: A two-fluid approach

We study the shock structure and acceleration efficiency of cosmic-ray mediated Magnetohydrodynamic (MHD) shocks both analytically and numerically by using a two-fluid model. Our model includes the dynamical effect of magnetic fields and cosmic rays on a background thermal fluid. The steady state solution is derived by following the technique of Drury & Voelk (1981) and compared to numerical results. We explore the time evolution of plane-perpendicular, piston-driven shocks. From the results of analytical and numerical studies, we conclude that the mean magnetic field plays an important role in the structure and acceleration efficiency of cosmic-ray mediated MHD shocks. The acceleration of cosmic-ray particles becomes less efficient in the presence of strong magnetic pressure since the field makes the shock less compressive. This feature is more prominent at low Mach numbers than at high Mach numbers.

Jun, Byung-Il↗
NASA NTRS2017

Validation of a Solid Rocket Motor Internal Environment Model

In a prior effort, a thermal/fluid model of the interior of Penn State University's laboratory-scale Insulation Test Motor (ITM) was constructed to predict both the convective and radiative heat transfer to the interior walls of the ITM with a minimum of empiricism. These predictions were then compared to values of total and radiative heat flux measured in a previous series of ITM test firings to assess the capabilities and shortcomings of the chosen modeling approach. Though the calculated fluxes reasonably agreed with those measured during testing, this exercise revealed means of improving the fidelity of the model to, in the case of the thermal radiation, enable direct comparison of the measured and calculated fluxes and, for the total heat flux, compute a value indicative of the average measured condition. By replacing the P1-Approximation with the discrete ordinates (DO) model for the solution of the gray radiative transfer equation, the radiation intensity field in the optically thin region near the radiometer is accurately estimated, allowing the thermal radiation flux to be calculated on the heat-flux sensor itself, which was then compared directly to the measured values. Though the fully coupling the wall thermal response with the flow model was not attempted due to the excessive computational time required, a separate wall thermal response model was used to better estimate the average temperature of the graphite surfaces upstream of the heat flux gauges and improve the accuracy of both the total and radiative heat flux computations. The success of this modeling approach increases confidence in the ability of state-of-the-art thermal and fluid modeling to accurately predict SRM internal environments, offers corrections to older methods, and supplies a tool for further studies of the dynamics of SRM interiors.

Martin, Heath T.↗
NASA NTRS2017

Validation of a Solid Rocket Motor Internal Environment Model

In a prior effort, a thermal/fluid model of the interior of Penn State University's laboratory-scale Insulation Test Motor (ITM) was constructed to predict both the convective and radiative heat transfer to the interior walls of the ITM with a minimum of empiricism. These predictions were then compared to values of total and radiative heat flux measured in a previous series of ITM test firings to assess the capabilities and shortcomings of the chosen modeling approach. Though the calculated fluxes reasonably agreed with those measured during testing, this exercise revealed means of improving the fidelity of the model to, in the case of the thermal radiation, enable direct comparison of the measured and calculated fluxes and, for the total heat flux, compute a value indicative of the average measured condition. By replacing the P1-Approximation with the discrete ordinates (DO) model for the solution of the gray radiative transfer equation, the radiation intensity field in the optically thin region near the radiometer is accurately estimated, allowing the thermal radiation flux to be calculated on the heat-flux sensor itself, which was then compared directly to the measured values. Though the fully coupling the wall thermal response with the flow model was not attempted due to the excessive computational time required, a separate wall thermal response model was used to better estimate the average temperature of the graphite surfaces upstream of the heat flux gauges and improve the accuracy of both the total and radiative heat flux computations. The success of this modeling approach increases confidence in the ability of state-of-the-art thermal and fluid modeling to accurately predict SRM internal environments, offers corrections to older methods, and supplies a tool for further studies of the dynamics of SRM interiors.

Martin, Heath T.↗
DOE OSTI2021

A modified many-body dissipative particle dynamics model for mesoscopic fluid simulation: methodology, calibration, and application for hydrocarbon and water

The many-body dissipative particle dynamics (mDPD) is a prominent mesoscopic multiphase model for fluid transport in mesoconfinement. However, it has been a long-standing challenge for mDPD (and other multiphase-enabled DPD models) to accurately predict real-fluid static and dynamic properties simultaneously. We have developed a modified mDPD model to overcome the issue and a rigorous calibration approach that uses reference data, including experimental and/or molecular dynamics (MD) simulations to parameterise the modified mDPD for real fluids. We choose heptane as a representative example of hydrocarbon in source rocks to demonstrate the model's capability to accurately predict the equation of state (EOS), free surface tension, diffusivity, and viscosity. Our timing test shows that the modified mDPD is 400–500 times faster than its MD counterpart for simulating bulk heptane in equivalent volumes. To further demonstrate the robustness of the model, we revisited the benchmark problem of mesoscopic modelling of water, in which all the previous works on DPD reported only a limited portion of the water properties. Here, we show that the modified mDPD can provide accurate modelling of water static and dynamic properties and an EOS that matches the experimental data to a large range of confinement pressure.

42 ENGINEERING↗
DOE OSTI2021

Development of a Filtered CFD-DEM Drag Model with Multiscale Markers Using an Artificial Neural Network and Nonlinear Regression

Here, the accuracy of coarse-grained Euler-Lagrangian simulations of fluidized beds heavily depends on the mesoscale drag models to account for the influences of the unresolved sub-grid structures. Traditional filtered drag models are regressed with mesoscale markers such as voidage and slip velocities. In this research, a filtered drag was regressed with both mesoscale and macro-scale markers using fine grid Computational Fluid Dynamics - Discrete Element Method (CFD-DEM) simulations. The traditional non-linear regression method was compared with machine learning regression using an Artificial Neural Network (ANN) implemented in PyTorch and coupled with MFiX. The new drag showed higher accuracy than the Wen-Yu drag and another filtered drag derived from the two-fluid model. The nonlinear regression shows slightly better results than ANN regression in cases with similar R 2 values. The utilization of the gas inlet velocity as an additional macro-scale marker reduced the errors by up to 55.3% in the tested cases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2004

Development of a Computationally Efficient, High Fidelity, Finite Element Based Hall Thruster Model

This report documents the development of a two dimensional finite element based numerical model for efficient characterization of the Hall thruster plasma dynamics in the framework of multi-fluid model. Effect of the ionization and the recombination has been included in the present model. Based on the experimental data, a third order polynomial in electron temperature is used to calculate the ionization rate. The neutral dynamics is included only through the neutral continuity equation in the presence of a uniform neutral flow. The electrons are modeled as magnetized and hot, whereas ions are assumed magnetized and cold. The dynamics of Hall thruster is also investigated in the presence of plasma-wall interaction. The plasma-wall interaction is a function of wall potential, which in turn is determined by the secondary electron emission and sputtering yield. The effect of secondary electron emission and sputter yield has been considered simultaneously, Simulation results are interpreted in the light of experimental observations and available numerical solutions in the literature.

Jacobson, David↗
NASA NTRSDate not supplied

Thermal Design Challenges Posed by the Four Bed CO2 Scrubber COTS Air-Save Pump

The Four Bed Carbon Dioxide (4BCO2) scrubber Air-Save Pump (ASP) operates as part of the adsorbent bed regeneration cycle. The ASP removes residual air from the bed for return to the cabin prior to heat and vacuum exposure which drives out the CO2 regenerating the bed. 4BCO2 employs a Commercial Off-the-Shelf (COTS) scroll type air pump, repackaged in an acoustically insulated enclosure to reduce noise and mounted to a cold plate. The International Space Station (ISS) Low Temperature Loop (LTL) flow (operating between 38°F and 50°F) first cools the process air, and then flows through the cold plate, cooling the pump. This results in competing ASP thermal design goals: (1) to keep the pump and motor sufficiently cool and (2) to avoid forming condensation due to over-cooling. Surfaces below 60°F typically warrant careful consideration of condensation. A test-calibrated thermal model demonstrates such a balanced design is feasible with temperatures above 60°F. A separate, coupled fluid model predicts the potential for condensation formation, allowing risk assessment of flying with the unmodified design.

Daniel P Popok↗