SEARCH · Engineering Papers

Results for “Fluids”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 883 records · Page 49

All sources NASA NTRS DOE OSTI

DOE OSTI2015

Fluid vs. kinetic magnetic reconnection with strong guide fields

The fast rates of magnetic reconnection found in both nature and experiments are important to understand theoretically. Recently, it was demonstrated that two-fluid magnetic reconnection remains fast in the strong guide field regime, regardless of the presence of fast-dispersive waves. This conclusion is in agreement with recent results from kinetic simulations, and is in contradiction to the findings in an earlier two-fluid study, where it was suggested that fast-dispersive waves are necessary for fast reconnection. In this paper, we give a more detailed derivation of the analytic model presented in a recent letter and present additional simulation results to support the conclusions that the magnetic reconnection rate in this regime is independent of both collisional dissipation and system-size. In particular, we present a detailed comparison between fluid and kinetic simulations, finding good agreement in both the reconnection rate and overall length of the current layer. Finally, we revisit the earlier two-fluid study, which arrived at different conclusions, and suggest an alternative interpretation for the numerical results presented therein.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Unsteady hydromagnetic flow of dusty fluid and heat transfer over a vertical stretching sheet with thermal radiation

In this paper, the hydromagnetic flow of dusty fluid over a vertical stretching sheet with thermal radiation is investigated. The governing partial differential equations are reduced to nonlinear ordinary differential equations using similarity transformation. These nonlinear ordinary differential equations are solved numerically using Runge-Kutta Fehlberg fourth-fifth order method (RKF45 Method). The behavior of velocity and temperature profiles of hydromagnetic fluid flow of dusty fluid is analyzed and discussed for different parameters of interest such as unsteady parameter, fluid-particle interaction parameter, the magnetic parameter, radiation parameter and Prandtl number on the flow.

42 ENGINEERING↗

DOE OSTI2026

Filament extension atomizers: A study of integration potential into spray dryers for high-viscosity and non-Newtonian fluids

Filament extension atomizers (FEAs) are an emerging class of spray nozzles designed to atomize high-viscosity and non-Newtonian fluids that are challenging for conventional pressure or two-fluid nozzles. In this study, we investigate the influence of roller geometry, surface velocity, roller material, and fluid rheology on the atomization performance of FEA systems using concentrated whey protein suspensions (50–70 wt. %). Extensional and shear rheology experiments, along with high-speed imaging and particle image velocimetry, reveal that filament breakup dynamics are governed by competition between inertial, capillary, and viscoelastic stresses. High roller rotational velocity leads to narrower spray cones, contradicting rheology experiments and suggests a significant inertial contribution to filament breakup. Smaller rollers operating at the same rotational velocity led to broader spray cones consistent with expectations. An FEA nozzle was integrated into a conventional dryer producing particles in the 100 μm range. Results suggest that FEA technology enables atomization of highly viscous fluids at industrially relevant spray cone angles similar to those generated by pressure nozzles, offering a pathway to improve energy efficiency in spray drying by enabling higher solids loading feedstocks. Furthermore, these insights provide critical guidance for optimizing FEA nozzle designs and process parameters across a range of applications.

Energy efficiency↗

DOE OSTI2025

Thermal-Fluid and Thermal-Structural Response of the T-Tube Modular Divertor to Spatiotemporally Varying Heat Loads

Tungsten (W) is the leading candidate for divertor target plates because of its high melting point (>3000°C), thermal conductivity, and ultimate tensile stress. While W and its alloys are the only solid materials that can survive the high heat fluxes incident on the divertor, W’s low-ductility high ductile-to-brittle transition temperature of ~600°C and relatively low recrystallization temperature (RT) of ~1300°C pose structural (among other) challenges. The objective of this work is to estimate the thermal-fluid and thermal-structural performance of the helium (He)-cooled T-tube divertor, which was originally developed by the Advanced Reactor Innovation and Evaluation Study (ARIES) using numerical simulations. Here, predictions of temperature distributions across the plasma-facing structural component and surface pressures from computational fluid dynamics simulations are used to determine stress distributions using commercial structural finite element modeling software over a range of fusion-relevant conditions. The maximum allowable incident heat fluxes are determined based on the temperature limits imposed by the ITER elastic Structural Design Criteria for In-vessel Components (SDC-IC) and the maximum RT over a range of He mass flow rates and presented in the form of performance design charts. Our recent work found that thermal- structural criteria accounting for the low ductility of W in a finger-type modular divertor constrain the maximum incident heat fluxes to values well below the ITER specifications, and those based on considering only the RT demonstrate that integrated thermal-fluid and elastic structural performance evaluation are required for accurate assessment of divertor performance. This novel analysis of the T-tube considers how nonuniform and transient incident heat fluxes affect its thermal-fluid and thermal-structural performance, as well as the effect of volumetric heating, which can be as great as 27% of the power incident on the divertor surface. The W tile of the T-tube, with its relatively large plasma-facing area of ~15 cm 2 , will likely experience significant spatial variations in incident heat flux. This work therefore assesses whether steady-state incident heat flux profiles with a peak of 10 MW/m 2 and maximum heat flux gradients of 200 MW/m 2 per m exceed the structural limits imposed by the ITER elastic SDC-IC and the maximum RT over a range of fusion-relevant conditions. The effect of transient heat fluxes typical of plasma detachment and reattachment from the target plate due, for example, to gas injection are also evaluated

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Local well-posedness and singularity formation in non-Newtonian compressible fluids

We investigate the initial value problem of a very general class of 3 + 1 non-Newtonian compressible fluids in which the viscous stress tensor with shear and bulk viscosity relaxes to its Navier–Stokes values. These fluids correspond to the non-relativistic limit of well-known Israel–Stewart-like theories used in the relativistic fluid dynamic simulations of high-energy nuclear and astrophysical systems. After establishing the local well-posedness of the Cauchy problem, we show for the first time in the literature that there exists a large class of initial data for which the corresponding evolution breaks down in finite time due to the formation of singularities. This implies that a large class of non-Newtonian fluids do not have finite solutions defined at all times.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Parametric electroconvection in a weakly conducting fluid in a horizontal parallel-plate capacitor

We study the flows of a nonuniformly heated weakly conducting fluid in an ac electric field of a horizontal parallel-plate capacitor. Analysis is carried out for fluids in which the charge formation is governed by electroconductive mechanism associated with the temperature dependence of the electrical conductivity of the medium. Periodic and chaotic regimes of fluid flow are investigated in the limiting case of instantaneous charge relaxation and for a finite relaxation time. Bifurcation diagrams and electroconvective regimes charts are constructed. The regions where fluid oscillations synchronize with the frequency of the external field are determined. Hysteretic transitions between electroconvection regimes are studied. The scenarios of transition to chaotic oscillations are analyzed. Depending on the natural frequency of electroconvective system and the external field frequency, the transition from periodic to chaotic oscillations can occur via quasiperiodicity, a subharmonic cascade, or intermittence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

Minkowski and Galilei/Newton Fluid Dynamics: A Geometric 3 + 1 Spacetime Perspective

A kinetic theory of classical particles serves as a unified basis for developing a geometric 3 + 1 spacetime perspective on fluid dynamics capable of embracing both Minkowski and Galilei/Newton spacetimes. Parallel treatment of these cases on as common a footing as possible reveals that the particle four-momentum is better regarded as comprising momentum and inertia rather than momentum and energy; and, consequently, that the object now known as the stress-energy or energy-momentum tensor is more properly understood as a stress-inertia or inertia-momentum tensor. In dealing with both fiducial and comoving frames as fluid dynamics requires, tensor decompositions in terms of the four-velocities of observers associated with these frames render use of coordinate-free geometric notation not only fully viable, but conceptually simplifying. A particle number four-vector, three-momentum ( 1 , 1 ) tensor, and kinetic energy four-vector characterize a simple fluid and satisfy balance equations involving spacetime divergences on both Minkowski and Galilei/Newton spacetimes. Reduced to a fully 3 + 1 form, these equations yield the familiar conservative formulations of special relativistic and non-relativistic fluid dynamics as partial differential equations in inertial coordinates, and in geometric form will provide a useful conceptual bridge to arbitrary-Lagrange–Euler and general relativistic formulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

STATISTICAL DECOUPLING OF A LAGRANGIAN FLUID PARCEL IN NEWTONIAN COSMOLOGY

The Lagrangian dynamics of a single fluid element within a self-gravitational matter field is intrinsically non-local due to the presence of the tidal force. This complicates the theoretical investigation of the nonlinear evolution of various cosmic objects, e.g., dark matter halos, in the context of Lagrangian fluid dynamics, since fluid parcels with given initial density and shape may evolve differently depending on their environments. In this paper, we provide a statistical solution that could decouple this environmental dependence. After deriving the evolution equation for the probability distribution of the matter field, our method produces a set of closed ordinary differential equations whose solution is uniquely determined by the initial condition of the fluid element. Mathematically, it corresponds to the projected characteristic curve of the transport equation of the density-weighted probability density function (ρPDF). Consequently it is guaranteed that the one-point ρPDF would be preserved by evolving these local, yet nonlinear, curves with the same set of initial data as the real system. Physically, these trajectories describe the mean evolution averaged over all environments by substituting the tidal tensor with its conditional average. For Gaussian distributed dynamical variables, this mean tidal tensor is simply proportional to the velocity shear tensor, and the dynamical system would recover the prediction of the Zel’dovich approximation (ZA) with the further assumption of the linearized continuity equation. For a weakly non-Gaussian field, the averaged tidal tensor could be expanded perturbatively as a function of all relevant dynamical variables whose coefficients are determined by the statistics of the field.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Experimental and numerical studies of seismic fluid-structure interaction in a base-supported cylindrical vessel

Seismic design and qualification of advanced reactors will rely heavily on the use of verified and validated numerical models capable of capturing the interaction of the vessel, its contained fluid, and the internal equipment: fluid-structure interaction (FSI) analysis. Analytical solutions can be used for preliminary sizing and design of such vessels but their application is limited to simple geometries and boundary conditions, and small amplitude, translational (and rotational) inputs. To validate numerical models for seismic FSI analysis in finite element codes, a comprehensive set of experiments was performed on a liquid-filled cylindrical vessel, using a 6 degree-of-freedom earthquake simulator. Results in terms of sloshing frequency, damping ratio in sloshing modes, and hydrodynamic responses (wave height, hydrodynamic pressure, base shear, and base moment) for multi-directional earthquake simulator inputs are reported and compared with analytical solutions for liquid-filled vessels. The impact of seismic (base) isolation on hydrodynamic responses was studied using earthquake simulator inputs generated using a virtual isolation system. Data from the experiments are used to validate a numerical model of the fluid-filled vessel using the Arbitrary Lagrangian Eulerian (ALE) solver in the commercial finite element program LS-DYNA. Validation studies are presented for multi-directional seismic inputs, including rocking motions. Lagrangian modeling approaches using an elastic material formulation for the fluid are also investigated and their limitations and possible applications are identified. Here, the results are broadly applicable to the seismic response of base supported, liquid-filled vessels.

58 GEOSCIENCES↗

DOE OSTI2022

Development, verification, and validation of comprehensive acoustic fluid-structure interaction capabilities in an open-source computational platform

The acoustic fluid-structure interaction (FSI) formulation is a practical numerical approach for the seismic analysis of fluid-filled tanks. However, there are no verification and validation studies reported in the literature that demonstrate the ability of an acoustic FSI numerical model to predict responses important to structural and mechanical design for intense translational and rotational earthquake inputs. Herein, an acoustic FSI formulation is implemented in the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE), and is formally verified and validated using analytical solutions and code-to-code verification, and experimental data, respectively. The analytical solutions are for small amplitude, unidirectional seismic inputs. The code-to-code verification utilizes a previously verified and validated Arbitrary Lagrangian-Eulerian (ALE) numerical model in the commercial finite element code LS-DYNA. The validation studies utilize a comprehensive data set assembled from results of 3D earthquake-simulator tests of a fluid-filled vessel. The acoustic numerical model in MOOSE is verified and validated for hydrodynamic pressures and support reactions except for cases that involve significant convective response. For small amplitude inputs, numerically predicted wave heights match those of the analytical solutions. The numerical model is not verified and validated for wave height calculations under intense 3D seismic inputs. The run times for the acoustic FSI simulations in MOOSE are an order of magnitude, or more, shorter than for the corresponding ALE simulations in LS-DYNA. The utility of the MOOSE acoustic FSI implementation is demonstrated by seismic analysis of a building equipped with a fluid-filled, advanced nuclear reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗

DOE OSTI2022

A two-fluid single-column model of turbulent shallow convection. Part II: Single-column model formulation and numerics

The two-fluid single-column model of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) is extended to include moisture and horizontal wind shear. Turbulent kinetic energy is introduced as a prognostic variable, dependence on a diagnosed boundary-layer height is removed, and subfilter fluxes are approximated using a two-fluid version of a Mellor–Yamada scheme. Three mechanisms for entrainment and detrainment processes are introduced, which represent entrainment of unstable air at the surface, forced detrainment of air at the top of the boundary/cloud layers, and turbulent mixing that relaxes the convective fluid to a reference profile. A semi-implicit Eulerian discretization replaces the semi-implicit semi-Lagrangian implementation of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) to improve numerical stability and conservation. The equations for the implicit time step are solved using a quasi-Newton method, which is shown to perform well in numerical tests for conservation and convergence. The two-fluid single-column model presented in this article will be applied to simulations of shallow cumulus convection in Part III.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Ensemble Kalman filter for data assimilation coupled with low-resolution computations techniques applied in fluid dynamics

This paper presents an innovative Reduced-order model (ROM) for merging experimental and simulation data using data assimilation (DA) to estimate the "True" state of a fluid dynamics system, leading to more accurate predictions. Our methodology introduces a novel approach by implementing the ensemble Kalman filter (EnKF) within a reduced-dimensional framework, grounded in a robust theoretical foundation and applied to fluid dynamics. To address the substantial computational demands of DA, the proposed ROM employs low-resolution (LR) techniques to drastically reduce computational costs. This innovative approach involves downsampling datasets for DA computations, followed by an advanced reconstruction technique based on low-cost singular value decomposition (lcSVD). The lcSVD method, a key innovation in this paper, has never been applied to DA before and offers a highly efficient way to enhance resolution with minimal computational resources. Our results demonstrate significant reductions in both computation time and RAM usage through these LR techniques without compromising the accuracy of the estimations. For instance, in a turbulent test case, for a data compression rate of 15.9, the LR approach can achieve a speed-up of 13.7 and a RAM compression of 90.9% while maintaining a low relative root mean square error (RRMSE) of 2.6%, compared to 0.8% in the high-resolution (HR) reference. Furthermore, we highlight the effectiveness of the EnKF in estimating and predicting the state of fluid flow systems based on limited observations and given low-fidelity numerical data. This paper highlights the potential of the proposed DA method in fluid dynamics applications, particularly for improving computational efficiency in CFD and related fields. Its ability to balance accuracy with low computational and memory costs makes it especially suitable for large-scale and real-time applications, such as environmental monitoring or engineering design. This method will be incorporated into ModelFLOWs-app.

Data Assimilation↗

DOE OSTI2014

Axial and transverse acoustic radiation forces on a fluid sphere placed arbitrarily in Bessel beam standing wave tweezers

The axial and transverse radiation forces on a fluid sphere placed arbitrarily in the acoustical field of Bessel beams of standing waves are evaluated. The three-dimensional components of the time-averaged force are expressed in terms of the beam-shape coefficients of the incident field and the scattering coefficients of the fluid sphere using a partial-wave expansion (PWE) method. Examples are chosen for which the standing wave field is composed of either a zero-order (non-vortex) Bessel beam, or a first-order Bessel vortex beam. It is shown here, that both transverse and axial forces can push or pull the fluid sphere to an equilibrium position depending on the chosen size parameter ka (where k is the wave-number and a the sphere’s radius). The corresponding results are of particular importance in biophysical applications for the design of lab-on-chip devices operating with Bessel beams standing wave tweezers. Moreover, potential investigations in acoustic levitation and related applications in particle rotation in a vortex beam may benefit from the results of this study. -- Highlights: •The axial and transverse forces on a fluid sphere in acoustical Bessel beams tweezers are evaluated. •The attraction or repulsion to an equilibrium position in the standing wave field is examined. •Potential applications are in particle manipulation using standing waves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Gravitational collapse of dissipative fluid as a source of gravitational waves

Gravitational collapse of cylindrical anisotropic fluid has been considered in analogy with the work of Misner and Sharp. Using Darmois matching conditions, the interior cylindrical dissipative fluid (in the form of shear viscosity and heat flux) is matched to an exterior vacuum Einstein–Rosen space–time. It is found that on the bounding 3-surface the radial pressure of the anisotropic perfect fluid is linearly related to the shear viscosity and the heat flux of the dissipative fluid on the boundary. This non-zero radial pressure on the bounding surface may be considered as the source of gravitational waves outside the collapsing matter distribution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

Modeling Concept and Numerical Simulation of Ultrasonic Wave Propagation in a Moving Fluid-Structure Domain based on a Monolithic Approach

Here in the present study, we propose a novel multiphysics model that merges two time-dependent problems – the Fluid-Structure Interaction (FSI) and the ultrasonic wave propagation in a fluid-structure domain with a one directional coupling from the FSI problem to the ultrasonic wave propagation problem. This model is referred to as the “eXtended fluid-structure interaction (eXFSI)” problem. This model comprises isothermal, incompressible Navier-Stokes equations with nonlinear elastodynamics using the Saint-Venant Kirchhoff solid model. The ultrasonic wave propagation problem comprises monolithically coupled acoustic and elastic wave equations. To ensure that the fluid and structure domains are conforming, we use the ALE technique. The solution principle for the coupled problem is to first solve the FSI problem and then to solve the wave propagation problem. Accordingly, the boundary conditions for the wave propagation problem are automatically adopted from the FSI problem at each time step. The overall problem is highly nonlinear, which is tackled via a Newton-like method. The model is verified using several alternative domain configurations. To ensure the credibility of the modeling approach, the numerical solution is contrasted against experimental data.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Estimation of the Thermal-Fluids and Thermal-Structural Performance of Helium-Cooled Modular Finger-Type Divertors

Over the past decade, our group has investigated the thermal-fluid performance of the helium-cooled modular divertor with multiple jets (HEMJ) and a simplified “flat” design of the HEMJ for long-pulse magnetic fusion energy (MFE) reactors. Experimental studies were performed in a helium (He) loop at the prototypical pressure of 10 MPa, nearly prototypical He temperatures and incident heat fluxes using test sections made from stainless steel and tungsten alloys. Correlations for average Nusselt numbers and pressure loss coefficients were developed from the data and are used to validate computational fluid dynamics (CFD) models. Here, this work presents updated thermal-fluids performance curves based on these correlations that estimate the maximum heat flux that can be accommodated by the plasma-facing surface and coolant pumping power requirements at prototypical operating conditions. Thermal-structural performance curves developed from ITER structural design criteria are introduced, which include protection against ductile and non-ductile failure, ratcheting fatigue, and creep fatigue. The performance design curves for these finger-type divertors demonstrate that the “flat” design, with a significantly less complicated geometry than the HEMJ, has thermal-fluid and thermal-structural performance comparable to the original HEMJ concept, and may be superior in terms of non-ductile failure criteria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Pool boiling heat transfer evaluation of next-generation dielectric fluid: Opteon™ 2P50

The growing use of artificial intelligence has led to heavy thermal loads and high heat dissipation rates in data centers. Conventional air-cooled technologies are not able to fulfill these requirements. To overcome these challenges, two-phase immersion cooling (2PIC) has emerged as one of the leading technologies for high power-density chips. 2PIC increases the heat dissipation rate and efficiency of the system while reducing the footprint of the cooling equipment. A fluid with adequate dielectric properties, a suitable normal boiling temperature to maintain chip temperatures, and good material compatibility, is desired for 2PIC system. In this study, the pool boiling heat transfer of a new developmental dielectric fluid, Opteon™ 2P50, was experimentally investigated. The heat transfer coefficients at various heat fluxes (20–150 kW/m 2 ) and the critical heat flux were measured using a smooth aluminum surface. Compared with HFE-7100, Opteon™ 2P50 shows higher heat transfer coefficient (up to 59% higher) and a slightly lower value of critical heat flux (around 5.9% lower). The modified Cooper correlation with the optimized leading constant resulted in reliable prediction accuracy with a 5.3% mean absolute error percentage. Overall, these results indicate that the new dielectric fluid provides similar thermal performance to some legacy fluids.

2P50↗