Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLUID BOUNDARY”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

An immersed boundary method for fluid–structure interaction with compressible multiphase flows

This paper presents a two-dimensional immersed boundary method for fluid–structure interaction with compressible multiphase flows involving large structure deformations. This method involves three important parts: flow solver, structure solver and fluid–structure interaction coupling. In the flow solver, the compressible multiphase Navier–Stokes equations for ideal gases are solved by a finite difference method based on a staggered Cartesian mesh, where a fifth-order accuracy Weighted Essentially Non-Oscillation (WENO) scheme is used to handle spatial discretization of the convective term, a fourth-order central difference scheme is employed to discretize the viscous term, the third-order TVD Runge–Kutta scheme is used to discretize the temporal term, and the level-set method is adopted to capture the multi-material interface. In this work, the structure considered is a geometrically non-linear beam which is solved by using a finite element method based on the absolute nodal coordinate formulation (ANCF). The fluid dynamics and the structure motion are coupled in a partitioned iterative manner with a feedback penalty immersed boundary method where the flow dynamics is defined on a fixed Lagrangian grid and the structure dynamics is described on a global coordinate. We perform several validation cases (including fluid over a cylinder, structure dynamics, flow induced vibration of a flexible plate, deformation of a flexible panel induced by shock waves in a shock tube, an inclined flexible plate in a hypersonic flow, and shock-induced collapse of a cylindrical helium cavity in the air), and compare the results with experimental and other numerical data. The present results agree well with the published data and the current experiment. Finally, we further demonstrate the versatility of the present method by applying it to a flexible plate interacting with multiphase flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Exocortex Network for AI-Augmented Human-Led Scientific Expedition

AI advances in science can be viewed along two main directions with a fluid boundary: enhancing efficiency through automation and smart tools to accelerate tasks that humans can already perform; and enabling exploration into uncharted territories and potentially toward AGI. These advances manifest in the AI cognitive core through the development and explainability of foundation models; in the physical embodiment of instruments and facilities; and in the integrated agency of AI workflows exemplified by the science exocortex. To address the role of humans in this evolving landscape, in this Perspective, we suggest a third direction: the development of personalized agents that form human-centered networks, supporting both efficiency and exploration while ensuring that AI remains aligned with human vision.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Towards self-consistent integrated modeling of the tokamak pedestal, scrape-off layer, and divertor using SOLPS-ITER and EPED

A predictive, self-consistent integrated model of the tokamak edge plasma and neutral system has been developed and validated against experimental data from DIII-D, spanning from the divertor target and scrape-off layer to the top of the pedestal. The model uses the SOLPS-ITER fluid boundary plasma and neutral code to constrain the EPED prediction of the pedestal height and width by taking the separatrix temperature and density calculated from SOLPS as inputs to EPED. An empirical proxy function based on data from a subset of these discharges is then used to infer the ratio of density at the pedestal top to that at the separatrix based on the electron temperature at the divertor target. This allows the density at the top of the pedestal to be provided to EPED as well, so that knowledge of this quantity is no longer required a priori, while SOLPS radial transport coefficients are needed instead. The gas puff rate provided to SOLPS is then a primary model input for this forward model coupling (along with radial transport coefficients, which are kept fixed through each scan in this study). Comparisons with measurements from DIII-D density scans through a range of target conditions find significant pedestal pressure degradation as the density is increased towards detachment in three different divertor configurations, consistent with the experimental data in these ballooning-limited pedestal regimes.

core-edge integration↗
DOE OSTI2022

MAGNETICALLY SUSPENDED CANNED ROTOR PUMPS FOR THE INTEGRAL MOLTEN SALT REACTOR

This project developed a molten salt pump design for small modular reactors, concentrated solar, and Gen IV nuclear reactors with a magnetically levitated rotor that can operate at temperatures of up to 700 °C. This eliminates the need for rotating seals and roller element bearings which require maintenance and are prone to failure. The pump design is also more compact than existing molten salt pump designs which require long shafts to thermally isolate the electric motor. The aims of the project are to: 1) design and test the molten salt pump, 2) develop and test a high temperature position sensor needed for magnetic bearing control, and 3) develop and test high temperature coil fabrication methods for electric motor and magnetic bearing fabrication. To develop the high temperature coils, conductor materials were analyzed for oxidation, diffusion, cladding, and resistivity. Of the candidate conductor materials that met the requirements, silver was chosen. The long-term impact of radiation on silver resistivity was analyzed by developing a model of a generic molten salt fast reactor to model the conductor irradiation in worst case scenarios. The analysis showed that increased resistivity due to transmutation of silver decreases exponentially with the reactor reflector thickness. Different electrical insulation materials were tested and the final high temperature coils developed were able to withstand conductor to ground voltages up to 1800 V before insulation breakdown. Magnetic bearings require high-precision high-speed measurements of the rotor position to operate. Current commercial position sensors can only operate up to 550 °C so a novel position sensor was developed that can operate up to 800 °C. These sensors were also designed to measure the position of the metallic rotor through a thin metallic corrosion barrier that protects that coils, insulation, and magnetic cores so that the salt containment barrier does not have any penetrations. The high-temperature sensor was tested and had comparable performance to commercially available inductive position sensors. The sensor design and signal processing has led to several invention disclosures. The magnetically levitated molten salt pump design was completed and analyzed. The chloride salt test loop at Oak Ridge National Laboratory was chosen to test the pump and the pump hydraulics and power were designed for the test loop flow and pressure requirements. Classes of motor designs that do not utilize permanent magnetic were considered for the high-temperature motor. After extensive analysis of their performance, a synchronous reluctance motor was chosen for the pump. The design was refined and analyzed to optimize performance. The magnetic bearings were designed to meet the force and frequency requirements necessary to levitate the rotor and reject disturbances from rotor imbalance, motor forces, and hydraulic forces. Compliant mounts for the position sensors were designed that would not damage the ceramic sensor components and maintain sensor alignment over a wide temperature range. A thermal hydraulic analysis of the molten salt coolant flow used to cool the rotor and stator was performed to ensure that the maximum temperature in the rotor and stator will not exceed 725 °C. The corrosion barrier and casing were designed to maintain the salt fluid boundary by using compliant features to absorb large expansion mismatches due to different coefficients of thermal expansion (CTE). Finally, the salt drain and fill systems were designed along with the pump support structure, integration with the test loop, and thermal insulation and heating. Corrosion barrier fabrication experiments were performed to assess using laser cladding to apply a sub-millimeter corrosion barrier to the surface of a Fe-Co-V magnetic core material. The cladding method showed minimal mixing between materials at their interface and withstood temperature cycling without damage to the cladding. The technology commercialization is focused on the high-temperature position sensors and there has been significant interest from the nuclear industry in utilizing these sensors along with some interest from the aerospace industry. The project also led to a research collaboration with a fusion energy industrial partner to study barriers to scaling the pump design from the kW size to the MW size. This research collaboration will also study methods to create high-temperature coils that have conductor to ground breakdown voltages above 5 kV. The successful development of high-temperature electromagnetic coils and high-temperature position sensors greatly increases the temperature limits for electromagnetic devices that don’t rely on permanent magnets. This includes electric motors, magnetic bearings, linear actuators, rotary position and velocity sensing, actuated valves, and generators to name a few. This will fill a current need in molten salt reactors, concentrated solar, and fusion energy for these critical peripheral devices needed for practical reactor designs that do not currently have commercial solutions.

Hines, J Wesley↗
DOE OSTI2019

Universal molecular-kinetic scaling relation for slip of a simple fluid at a solid boundary

Using the observation that slip in simple fluids at low and moderate shear rates is a thermally activated process driven by the shear stress in the fluid close to the solid boundary, we develop a molecular-kinetic model for simple fluid slip at solid boundaries. Here, the proposed model, which is in the form of a universal scaling relation that connects slip and shear rate, reduces to the well-known Navier slip condition under low shear conditions, providing a direct connection between molecular parameters and the slip length. Molecular-dynamics simulations are in very good agreement with the predicted dependence of slip on system parameters, including the temperature and fluid-solid interaction strength. Connections between our model and previous work, as well as simulation and experimental results, are explored and discussed.

42 ENGINEERING↗
DOE OSTI2018

A novel algorithm of immersed moving boundary scheme for fluid-particle interactions in DEMLBM

This paper presents an efficient and accurate Immersed Moving Boundary (IMB) algorithm for solving fluid-particle interactions in the framework of the Lattice Boltzmann Method (LBM). Although the IMB scheme has been widely employed in many fluid-particle coupling problems in a wide range of applications, the algorithm of its implementation, especially in identifying both fluid and solid boundary nodes for particles, is seldom reported. Besides, the computational cost of handling fluid-particle coupling is very expensive. To provide a bridge between theory and application and improve the computing efficiency of IMB, a novel boundary tracing procedure and an efficient method for computing the solid nodal ratio using Gaussian quadrature are proposed in this paper. Both accuracy and efficiency of the proposed algorithm are examined by two benchmark tests. It is also found that the IMB scheme are more efficient and stable compared to another widely used the Immersed Boundary Method (IBM) in LBM.

42 ENGINEERING↗
DOE OSTI2025

Platinum liquid-vapor phase boundary mapped by fluid flyer experiments

We report a direct measurement of the temperature and density of a metal along its liquid-vapor coexistence (L-V) curve. By shocking platinum to a high-pressure liquid, we imparted sufficient heat for subsequent isentropic release to place it in a state on the boundary between the liquid and vapor phases. Released material in the liquid phase acted as a high velocity flyer pinned to the L-V curve. We measured velocity and radiant emission of the flyer as well as the interface motion and transiting shock states induced in a downstream window material by its impact. We used these measurements to calculate temperature and density of the L-V curve state which we compare to density functional theory predictions.

Equations of state↗
DOE OSTI2024

Pushing the Boundaries of Relativistic Fluid Dynamics: From Particle Colliders to Neutron Star Mergers

This project aimed at establishing connections between fluid dynamics, kinetic theory, and string theory to develop new tools to characterize the far-from-equilibrium properties of relativistic systems, which are instrumental in understanding the emergence of liquid-like behavior of the quark-gluon plasma formed in heavy-ion collisions. This project investigated the fundamental properties of viscous fluids in the full nonlinear regime to pave the way for the systematic inclusion of viscous effects in future numerical simulations of neutron star mergers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Inertial dynamics of an interface with interfacial mass flux: Stability and flow fields’ structure, inertial stabilization mechanism, degeneracy of Landau’s solution, effect of energy fluctuations, and chemistry-induced instabilities

This work focuses on the long-standing problem of inertial dynamics of an interface with interfacial mass flux and reports new mechanisms for the interface stabilization and destabilization. The interface is a phase boundary separating fluids of different densities and having interfacial mass flux. To analyze the interface dynamics from a far field, we develop and apply the general matrix method to rigorously solve the boundary value problem involving the governing equations in the fluid bulk and the boundary conditions at the interface and at the outside boundaries of the domain. We find the fundamental solutions for the linearized system of equations and analyze the interplay of interface stability with flow fields’ structure by directly linking rigorous mathematical attributes to physical observables. We find that the interface is stable when the dynamics conserves the fluxes of mass, momentum, and energy; the stabilization is due to an inertial mechanism causing small oscillations of the interface velocity. In the classic Landau’s dynamics, the postulate of perfect constancy of the interface velocity leads to the development of Landau–Darrieus instability. This destabilization is also linked to the imbalance of the perturbed energy at the interface. The classic Landau’s solution is found to have degeneracy; lifting of the degeneracy may lead to singularity and self-similar dynamics. Our results compare well with traditional theories of combustion and propose new experiments to study the dynamics of the interface and the flow fields in combustible systems. We further conduct reactive molecular dynamics simulations to elucidate the complexity of chemical processes, to study the destabilizing effect of energy fluctuations on the interface stability, and to illustrate the chemistry-induced instabilities. In summary, we identify the extreme sensitivity of the interface dynamics to the interfacial boundary conditions, including the formal properties of fundamental solutions and the qualitative and quantitative properties of the flow fields. This provides new opportunities for studies, diagnostics, and control of multiphase flows in a broad range of processes in nature and technology.

42 ENGINEERING↗
DOE OSTI2025

Underwater thermomagnetic generator for remote marine thermal energy harvesting and sensing

Thermomagnetic generators offer a promising approach for sustainable power generation in remote marine environments. Here, this study presents the design, thermal modeling, and experimental validation of a passively driven underwater thermomagnetic generator developed for powering ocean observation and monitoring sensors. The generator was evaluated under varying working fluids, thermal boundary conditions, and extended operation to assess real-world applicability. Two fluids, deionized water and silicone-based Thermal C5, were tested under simulated shallow- and deep-ocean conditions. Deionized water outperformed Thermal C5, especially in colder environments (~5°C), achieving a peak output of 2.7 mW due to larger temperature gradients and enhanced convective-evaporative heat transfer. Long-duration tests revealed a transient evaporation-condensation cycle that temporarily reduced rotor immersion and performance before stabilizing. The generator powered commercial marine sensors for over 6 h without external batteries, demonstrating the viability of compact, passively cooled thermomagnetic systems for autonomous, off-grid marine sensing.

13 HYDRO ENERGY↗
DOE OSTI2026

Nanoconfined Grain Boundaries Increase the Conductivity of Polycrystalline Molecular Crystals

Soft-solid molecular crystals consist of crystalline grains and fluid grain boundaries (GBs) that enhance the grain binding and transport of Li + ions between the grains. The total ionic conductivity consists of ion migration in both the grains and GBs. To unravel these contributions in adiponitrile (Adpn):LiPF 6 molecular crystals, the GB volume fraction was varied by changing the size of the crystals and the Adpn:LiPF 6 molar ratio. Molecular dynamics (MD) simulations indicate that ion motion was subdiffusive in the grains and “well-diffusive” in the GBs, with GBs characterized as disordered nanoconfined regions of higher charge carrier concentration (∼1 M) than in saturated Adpn:LiPF 6 solutions (0.04 M), and Li + ions predominantly solvated by cyano groups with few contact ion pairs. The diffusivity in the GBs is at least an order of magnitude higher than that in the crystalline grains. The emergent picture is the grains as a reservoir of ions that migrate to faster-conducting GBs.

Crystals↗
DOE OSTI2024

Recent developments in mathematical aspects of relativistic fluids

Abstract We review some recent developments in mathematical aspects of relativistic fluids. The goal is to provide a quick entry point to some research topics of current interest that is accessible to graduate students and researchers from adjacent fields, as well as to researches working on broader aspects of relativistic fluid dynamics interested in its mathematical formalism. Instead of complete proofs, which can be found in the published literature, here we focus on the proofs’ main ideas and key concepts. After an introduction to the relativistic Euler equations, we cover the following topics: a new wave-transport formulation of the relativistic Euler equations tailored to applications; the problem of shock formation for relativistic Euler; rough (i.e., low-regularity) solutions to the relativistic Euler equations; the relativistic Euler equations with a physical vacuum boundary; relativistic fluids with viscosity. We finish with a discussion of open problems and future directions of research.

Disconzi, Marcelo (ORCID:0000000234497778)↗
DOE OSTI2023

Gaussian process hydrodynamics

Abstract We present a Gaussian process (GP) approach, called Gaussian process hydrodynamics (GPH) for approximating the solution to the Euler and Navier-Stokes (NS) equations. Similar to smoothed particle hydrodynamics (SPH), GPH is a Lagrangian particle-based approach that involves the tracking of a finite number of particles transported by a flow. However, these particles do not represent mollified particles of matter but carry discrete/partial information about the continuous flow. Closure is achieved by placing a divergence-free GP prior ξ on the velocity field and conditioning it on the vorticity at the particle locations. Known physics (e.g., the Richardson cascade and velocity increment power laws) is incorporated into the GP prior by using physics-informed additive kernels. This is equivalent to expressing ξ as a sum of independent GPs ξ l , which we call modes, acting at different scales (each mode ξ l self-activates to represent the formation of eddies at the corresponding scales). This approach enables a quantitative analysis of the Richardson cascade through the analysis of the activation of these modes, and enables us to analyze coarse-grain turbulence statistically rather than deterministically. Because GPH is formulated by using the vorticity equations, it does not require solving a pressure equation. By enforcing incompressibility and fluid-structure boundary conditions through the selection of a kernel, GPH requires significantly fewer particles than SPH. Because GPH has a natural probabilistic interpretation, the numerical results come with uncertainty estimates, enabling their incorporation into an uncertainty quantification (UQ) pipeline and adding/removing particles (quanta of information) in an adapted manner. The proposed approach is suitable for analysis because it inherits the complexity of state-of-the-art solvers for dense kernel matrices and results in a natural definition of turbulence as information loss. Numerical experiments support the importance of selecting physics-informed kernels and illustrate the major impact of such kernels on the accuracy and stability. Because the proposed approach uses a Bayesian interpretation, it naturally enables data assimilation and predictions and estimations by mixing simulation data and experimental data.

Mathematics↗
DOE OSTI2017

Thermal-hydraulic performance of metal foam heat exchangers under dry operating conditions

High porosity metal foams with novel thermal, mechanical, electrical, and acoustic properties are being more widely adopted for application. Due to their large surface-area-to-volume ratio and complex structure which induces better fluid mixing, boundary layer restarting and wake destruction, they hold promise for heat transfer applications. In this study, the thermal-hydraulic performance of open-cell aluminum metal foam heat exchanger has been evaluated. The impact of flow conditions and metal foam geometry on the heat transfer coefficient and gradient have been investigated. Metal foam heat exchanger with same geometry (face area, flow depth and fin dimensions) consisting of four different type of metal foams have been built for the study. Experiments are conducted in a closed-loop wind tunnel at different flow rate under dry operating condition. Metal foams with a smaller pore size (40 PPI) have a larger heat transfer coefficient compared to foams with a larger pore size (5 PPI). However, foams with larger pores result in relatively smaller pressure gradients. Current thermal-hydraulic modeling practices have been reviewed and potential issues have been identified. Permeability and inertia coefficients are determined and compared to data reported in open literature. Finally, on the basis of the new experimental results, correlations are developed relating the foam characteristics and flow conditions through the friction factor f and the Colburn j factor.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Evolution of structure and dynamics of thermo-reversible nanoparticle gels—A combined XPCS and rheology study

A combined X-ray photon correlation spectroscopy and rheology report is carried out to capture the evolution of structure, fast particle-scale dynamics, and moduli (elastic and loss) at early times of gel formation near the fluid-gel boundary of a suspension of nanoparticles. The system is made up of moderately concentrated suspensions of octadecyl silica in decalin (φ = 0.2) undergoing thermoreversible gelation. Near the gel boundary, the rate of gel formation is very sensitive to changes in attraction strength. Yet, we find that at different attraction strengths, the system goes through identical intermediate states of microscopic and macroscopic behavior, even though the absolute time needed to form a gel varies by orders of magnitude. We identify a single dimensionless time parameter, t w /t g , where t w is the wait time following the quench and t g is the rheologically determined gel time, that captures the similarity in gel formation at a range of attraction strengths. Following a temperature quench below the gel boundary, the system is initially fluidlike and forms diffusive clusters (similar to 8.5 times the particle diameter). After a lag-time, t L , clusters aggregate to form a network like structure which is characterized by the onset of mechanical rigidity and a rapid growth in microscopic relaxation times. At t g , the Baxter parameter obtained from adhesive hard sphere fits of the structure factor attains a constant value corresponding to the theoretical percolation boundary, thus demonstrating that gelation is percolation driven.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

An assessment of thermodynamic merits for current and potential future engine operating strategies

The present work compares the fundamental thermodynamic underpinnings (i.e., working fluid properties and heat release profile) of various combustion strategies with engine measurements. The approach employs a model that separately tracks the impacts on efficiency due to differences in rate of heat addition, volume change, mass addition, and molecular weight change for a given combination of working fluid, heat release profile, and engine geometry. Comparative analysis between measured and modeled efficiencies illustrates fundamental sources of efficiency reductions or opportunities inherent to various combustion regimes. Engine operating regimes chosen for analysis include stoichiometric spark-ignited combustion and lean compression-ignited combustion including HCCI, SA-HCCI, RCCI, GCI, and CDC. Within each combustion regime, effects such as engine load, combustion duration, combustion phasing, combustion chamber geometry, fuel properties, and charge dilution are explored. Model findings illustrate that even in the absence of losses such as heat transfer or incomplete combustion, the maximum possible thermal efficiency inherent to each operating strategy varies to a significant degree. Additionally, the experimentally measured losses are observed to be unique within a given operating strategy. The findings highlight the fact that in order to create a roadmap for future directions in ICE technologies, it is important to not only compare the absolute real-world efficiency of a given combustion strategy, but to also examine the measured efficiency in context of what is thermodynamically possible with the working fluid and boundary conditions prescribed by a strategy.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Reservoir Dynamics in Proposed Operational Scenarios Where CO 2 -EOR Fields are Transitioned From CO 2 -Flood Enhanced Oil Recovery to Dedicated Carbon Storage: A Field Case Study

This study models the transition of CO 2 -enhanced oil recovery (CO 2 -EOR) fields to dedicated carbon storage using a generalized reservoir model informed by real field data, SACROC. Through scenario-based simulations, the work evaluates how reservoir depletion levels, boundary conditions, fluid properties, and domain size affect pressure buildup, CO 2 plume migration, and long-term containment performance. The results inform practices for repurposing oil fields into reliable CO 2 storage sites.

02 PETROLEUM↗