SEARCH · Engineering Papers

Results for “rotating fluid”

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 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Taylor-Couette Flow Modeling (Rev.1)

The Taylor-Couette Flow is a form of viscous flow in which fluid flows between two cylindrical walls with a prescribed rotation. A diagram is shown on the right (Figure 1) where Ω is the rotational velocity, H is the distance between the two cylinders, R is the radius the outer cylinder, r is the radius of the inner cylinder and L is the length of the cylinders. When the rotating inner cylinder reaches a critical velocity, the flow of the liquid becomes unstable and creates vortices in between the layers. In this experiment, the Taylor-Couette cylinders were modeled so that the vortices would be created after reaching a critical velocity. After the initial experiment, the velocity of the inner cylinder was decreased and increased to look at how the results changed.

42 ENGINEERING↗

DOE OSTI2021

The Polarized Image of a Synchrotron-emitting Ring of Gas Orbiting a Black Hole

Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov and conservation of the Walker–Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Computational Fluid Dynamics Simulation of Compressible Non-Newtonian Biomass in a Compression-Screw Feeder

Compression-screw feeders play a critical role in biorefineries to transport lignocellulosic biomass feedstocks from the feeding hoppers to biomass-conversion reactors in order to pretreat and convert biomass to hydrocarbon liquid biofuels and other power and energy resources. One of the main challenges in the operation of screw feeder is plugging and jamming of compressed biomass with high concentration of insoluble solids. The focus of this paper is to numerically investigate the screw feeder at these challenging operating conditions and help with the optimization of the screw feeder design to avoid operation failure. In this work a customized CFD model based on open-source OpenFOAM package [1] was developed to simulate the concentrated biomass as a highly viscous non-Newtonian fluid in the screw feeder. The biomass is modeled as a single-phase compressible Bingham fluid with a plastic viscosity as well as a density-dependent yield stress. The compressibility formulation (pressure-dependent density) and the density-dependent yield stress formulation in the governing equations follow the suggestions from a recent study by Duncan et al. [2]. A pilot-scale hopper/screw feeding system at NREL [3] is used to compare the experimental observations with our simulation results. The auger is 280 mm long and tapered with outer diameter changing from 80 mm to 35 mm. The auger rotates from 10 to 60 rpm in a conical throat which contains anti-rotational bars. The simulations predicted the required torque for the screw feeder and the pressure increase at the exit for biomass feedstocks with various fluid viscosity properties and auger rotating speeds. The analysis of the stress forces helped to identify the critical conditions were the screw feeder excessive wear or jamming could occur.

biofuels↗

DOE OSTI2020

Experimental Investigation of Hydraulic Fracturing and Stress Sensitivity of Fracture Permeability Under Changing Polyaxial Stress Conditions

Understanding and predicting fracture propagation and subsequent fluid flow characteristics is critical to geoenergy technologies that engineer and/or utilize favorable geological conditions to store or extract fluids from the subsurface. Fracture permeability decreases nonlinearly with increasing normal stress, but the relationship between shear displacement and fracture permeability is less well understood. We utilize the new Geo–Reservoir Experimental Analogue Technology (GREAT cell), which can apply polyaxial stress states and realistic reservoir temperatures and pressures to cylindrical samples and has the unique capability to alter both the magnitude and orientation of the radial stress field by increments of 11.25° during an experiment. We load synthetic analogue materials and real rock samples to stress conditions representative of 500–1,000 m depth, investigate the hydraulic stimulation process, and then conduct flow experiments while changing the fluid pressure and the orientation of the intermediate and minimum principal stresses. High–resolution circumferential strain measurements combined with fluid pressure data indicate fracture propagation can be both stable (no fluid pressure drop) and unstable (fluid pressure drop). The induced fractures exhibit both opening and shear displacements during their creation and/or during fluid flow with changing radial stress states. Flow tests during radial stress field rotation reveal that fracture normal effective stress has first–order control on fracture permeability but increasing fracture offset can lead to elevated permeabilities at maximum shear stress. Here, the results have implications for our conceptual understanding of fracture propagation as well as fluid flow and deformation around fractures.

58 GEOSCIENCES↗

DOE OSTI2020

Solvation Energy of Ions in a Stockmayer Fluid

We calculate the solvation energy of monovalent and divalent ions in various liquids with coarse-grained molecular dynamics simulations. Our theory treats the solvent as a Stockmayer fluid, which accounts for the intrinsic dipole moment of molecules and the rotational dynamics of the dipoles. Despite the simplicity of the model, we obtain qualitative agreement between the simulations and experimental data for the free energy and enthalpy of ion solvation, which indicates that the primary contribution to the solvation energy arises mainly from the first and possibly second solvation shells near the ions. Our results suggest that a Stockmayer fluid can serve as a reference model that enables direct comparison between theory and experiment and may be invoked to scale up electrostatic interactions from the atomic to the molecular length scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Global Λ-hyperon polarization in Au + Au collisions at $\sqrt{S_{NN}}$ = 3 GeV

Global hyperon polarization, $\overline{P}$ H , in Au + Au collisions over a large range of collision energy, $\sqrt{S_{NN}}$, was recently measured and successfully reproduced by hydrodynamic and transport models with intense fluid vorticity of the quark-gluon plasma. While na¨ıve extrapolation of data trends suggests a large $\overline{P}$ H as the collision energy is reduced, the behavior of $\overline{P}$ H at small $\sqrt{S_{NN}}$ < 7.7 GeV is unknown. Operating the STAR experiment in fixed-target mode, we measured the polarization of Λ hyperons along the direction of global angular momentum in Au + Au collisions at $\sqrt{S_{NN}}$ = 3 GeV. As a result, the observation of substantial polarization of 4.91 ± 0.81 (stat.) ± 0.15 (syst.) % in these collisions may require a reexamination of the viscosity of any fluid created in the collision, of the thermalization timescale of rotational modes, and of hadronic mechanisms to produce global polarization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Dispersion Validation for Flow Involving a Large Structure Revisited: 45 Degree Rotation

The atmospheric dispersion of contaminants in the wake of a large urban structure is a challenging fluid mechanics problem of interest to the scientific and engineering communities. Magnetic Resonance Velocimetry (MRV) and Magnetic Resonance Concentration (MRC) are relatively new techniques that leverage diagnostic equipment used primarily by the medical field to make 3D engineering measurements of flow and contaminant dispersal. SIERRA/Fuego, a computational fluid dynamics (CFD) code at Sandia National Labs is employed to make detailed comparisons to the dataset to evaluate the quantitative and qualitative accuracy of the model. This work is the second in a series of scenarios. In the prior work, a single large building in an array of similar buildings was considered with the wind perpendicular to a building face. In this work, the geometry is rotated by 45 degrees and improved studies are performed for simulation credibility. The comparison exercise shows conditionally good comparisons between the model and experiment. Model uncertainties are assessed through parametric variations. Various methods of quantifying the accuracy between experiments and data are examined Three-dimensional analysis of accuracy is performed. The effort helped identify deficiencies in the techniques used to make these comparisons, and further methods development therefore becomes one of the main recommendations for follow-on work.

42 ENGINEERING↗

DOE OSTI2024

Axisymmetric hybrid Vlasov equilibria with applications to tokamak plasmas

We derive axisymmetric equilibrium equations in the context of the hybrid Vlasov model with kinetic ions and massless fluid electrons, assuming isothermal electrons and deformed Maxwellian distribution functions for the kinetic ions. The equilibrium system comprises a Grad–Shafranov partial differential equation and an integral equation. These equations can be utilized to calculate the equilibrium magnetic field and ion distribution function, respectively, for given particle density or given ion and electron toroidal current density profiles. The resulting solutions describe states characterized by toroidal plasma rotation and toroidal electric current density. Additionally, due to the presence of fluid electrons, these equilibria also exhibit a poloidal current density component. This is in contrast to the fully kinetic Vlasov model, where axisymmetric Jeans equilibria can only accommodate toroidal currents and flows, given the absence of a third integral of the microscopic motion.

Physics↗

DOE OSTI2026

Anisotropic gravitational waves from inhomogeneous axion rotation

Gravitational waves (GWs) provide a powerful probe of the early universe due to their ability to free-stream across cosmic history. We study GW production in a compelling scenario where a rotating axion(-like) field becomes relevant for a brief period in the early universe before transitioning into a kination fluid and rapidly dissipating its energy through cosmic expansion. During this short epoch, the curvature perturbation can be predominantly sourced by the rotating axion and may significantly exceed the adiabatic component. Moreover, axion field perturbations grow on superhorizon scales during this phase. These effects can generate a strong stochastic background of induced GWs. This GW background also exhibits a pronounced large-scale anisotropy inherited from the axion fluctuations, serving as a distinctive signature of the scenario. Importantly, the transient nature of axion relevance enables this scenario to evade stringent bounds on large-scale perturbations. We analyze various observational constraints and find that both the amplitude and anisotropy of the resulting GW signal could be accessible to future detectors.

Axions and ALPs↗

DOE OSTI2023

Review of rotating wing dynamic stall: Experiments and flow control

Dynamic stall has been a technical challenge and a fluid dynamical subject of interest for more than fifty years; but in the last decade significant advances have been made in the understanding, prediction, modeling, and control of dynamic stall on rotors. This paper provides a summary of the state of the art of dynamic stall experiments and future directions in the understanding of dynamic stall on rotors. Experimental data sets are discussed, as well the direction of future research for control of dynamic stall. Coordinated testing between airfoils and rotating blades, as well as close integration between computational and experimental studies were found to be productive approaches. Advanced analysis methods, including statistical methods, modal representations, and artificial intelligence methods have led to significant advances in the understanding of dynamic stall. Investigations of dynamic stall control devices have allowed many useful targeted investigations of the transition to separated flow, but have not yet resulted in a commercially implemented device.

42 ENGINEERING↗

DOE OSTI2022

A Non-perturbative Approach to Computing Seismic Normal Modes in Rotating Planets

In this work, a continuous Galerkin method based approach is presented to compute the seismic normal modes of rotating planets. Special care is taken to separate out the essential spectrum in the presence of a fluid outer core using a polynomial filtering eigensolver. The relevant elastic-gravitational system of equations, including the Coriolis force, is subjected to a mixed finite-element method, while self-gravitation is accounted for with the fast multipole method. Our discretization utilizes fully unstructured tetrahedral meshes for both solid and fluid regions. The relevant eigenvalue problem is solved by a combination of several highly parallel and computationally efficient methods. We validate our three-dimensional results in the non-rotating case using analytical results for constant elastic balls, as well as numerical results for an isotropic Earth model from standard “radial” algorithms. We also validate the computations in the rotating case, but only in the slowly-rotating regime where perturbation theory applies, because no other independent algorithms are available in the general case. The algorithm and code are used to compute the point spectra of eigenfrequencies in several Earth and Mars models studying the effects of heterogeneity on a large range of scales.

58 GEOSCIENCES↗

DOE OSTI2021

Physics-informed machine learning of the Lagrangian dynamics of velocity gradient tensor

Reduced models describing the Lagrangian dynamics of the velocity gradient tensor (VGT) in homogeneous isotropic turbulence (HIT) are developed under the physics-informed machine learning (PIML) framework. We consider the VGT at both Kolmogorov scale and coarse-grained scale within the inertial range of HIT. Building reduced models requires resolving the pressure Hessian and subfilter contributions, which is accomplished by constructing them using the integrity bases and invariants of the VGT. The developed models can be expressed using the extended tensor basis neural network (TBNN) introduced by Ling et al. [J. Fluid Mech. 807, 155 (2016)]. Physical constraints, such as Galilean invariance, rotational invariance, and incompressibility condition, are thus embedded in the models explicitly. Our PIML models are trained on the Lagrangian data from a high-Reynolds number direct numerical simulation (DNS). To validate the results, we perform a comprehensive out-of-sample test. We observe that the PIML model provides an improved representation for the magnitude and orientation of the small-scale pressure Hessian contributions. Statistics of the flow, as indicated by the joint PDF of second and third invariants of the VGT, show good agreement with the “ground-truth” DNS data. A number of other important features describing the structure of HIT are reproduced by the model successfully. We have also identified challenges in modeling inertial range dynamics, which indicates that a richer modeling strategy is required. This helps us identify important directions for future research, in particular towards including inertial range geometry into the TBNN.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Rocking response of liquid-filled cylindrical tanks

The hydrodynamic response of liquid-filled storage tanks subjected to translational earthquake shaking has been the focus of studies for more than 75 years. In service, tanks experience six components of seismic input, with rocking due to support flexibility and/or rotational components of ground motion. Published theory for predicting hydrodynamic responses in rigid cylindrical tanks due to rotational seismic input is examined. Analytical solutions for pressure in the fluid and on the tank wall, and base moment are modified to account for the effects of angular displacements at the base. A finite element model of a water-filled cylindrical tank using the Arbitrary Lagrangian and Eulerian (ALE) solver in LS-DYNA is used to demonstrate the efficacy of the proposed expressions for three base rocking motions. Results from the modified analytical solutions and the finite element model are in excellent agreement.

58 GEOSCIENCES↗

DOE OSTI2021

Rotational and reflectional equivariant convolutional neural network for data-limited applications: Multiphase flow demonstration

This article deals with approximating steady-state particle-resolved fluid flow around a fixed particle of interest under the influence of randomly distributed stationary particles in a dispersed multiphase setup using convolutional neural network (CNN). The considered problem involves rotational symmetry about the mean velocity (streamwise) direction. Thus, this work enforces this symmetry using SE(3)-equivariant, special Euclidean group of dimension 3, CNN architecture, which is translation and three-dimensional rotation equivariant. This study mainly explores the generalization capabilities and benefits of a SE(3)-equivariant network. Accurate synthetic flow fields for Reynolds number and particle volume fraction combinations spanning over a range of [86.22, 172.96] and [0.11, 0.45], respectively, are produced with careful application of symmetry-aware data-driven approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Numerical Analysis of Combustion Dynamics in a Full-Scale Rotating Detonation Rocket Engine using Large Eddy Simulations

Large eddy simulations (LES) using detailed chemistry and leveraging adaptive mesh refinement (AMR) are performed to gain insights into the combustion dynamics within a full-scale methane-oxygen non-premixed rotating detonation rocket engine (RDRE) employing impinging discrete injection schemes. In particular, a comparative analysis of two operating conditions corresponding to the same global equivalence ratio but different mass flow rates is carried out to investigate the resultant impact on detonation wave characteristics and RDRE global performance. Multiple co-rotating detonation waves with spatially-distributed wave structure and preferential alignment with the inner wall of the annulus (due to asymmetry in fuel distribution) are encountered under both conditions. Both cases exhibit pre-detonation deflagrative burning in the fill region, while one of the cases shows higher susceptibility to backflow into the feed plenums due to lower plenum pressures. Furthermore, heat release analysis shows that the thrust obtained from the RDRE is closely linked to the distribution of total heat release between detonative and deflagrative combustion. On the other hand, combustion efficiency is associated with the fraction of heat release occurring in fuel-rich versus fuel-lean regions within the RDRE.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2022

Fluid-powered linear motor with rotary pistons and motion rectifier and synthetic diamond bearing assemblies

A fluid-powered linear motor with rotary pistons is disclosed. An application is for a downhole motor but it could be used in other applications. Rotational pistons provide increased torque generation as the torque generated is proportional to motor length. Since downhole drills are long (generally up to a maximum length of 30 ft.), a high-torque motor can be produced using this method. A pressurized fluid is used to drive the piston assemblies to produce bit shaft power. This concept employs rotary pistons and a mechanical rectifier to convert the rotational reciprocation produced by a reciprocating rotational piston into continuous rotary motion. Integration of synthetic diamond bearing assemblies into the motor allows for reduced friction and abrasive drilling fluid compatibility.

Raymond, David W.↗

DOE OSTI2021

Prediction of Detonation-Induced Disturbances Propagating Upstream into Inlets of Rotating Detonation Combustors

Disturbances caused by the detonation wave in a rotating detonation combustor (RDC) propagate upstream through the inlet, and can potentially affect and couple to upstream components, such as turbomachinery or isolators. These disturbances can potentially also affect the operation of the RDC itself. By drawing from the analogy of a detonation wave bounded by an inert gas, the pressure disturbances observed upstream of the inlet are explained as the consequence of the passage of an upstream propagating oblique shock. In this study, the pressure rise in the plenum from the oblique shock is measured in an axial air inlet RDC. The speed of the upstream propagating wave is estimated to be moderately above the acoustic speed of the oxidizer in the plenum. The wave propagates into the plenum despite local regions of choking in the inlet. It is estimated that the time it takes a fluid particle to transit from the plenum to the detonation channel through the inlet is much larger than the rotational time of the detonation wave. This implies that a fluid particle experiences multiple shocks prior to entering the detonation channel. The oblique shock propagating upstream through the inlet area change is modeled by leveraging an analogy with a quasi-1D shock wave moving in a variable area duct with mean (incoming) flow. Due to flow expansion along the area change, fluid particles are found to experience stronger shocks in the inlet than in the plenum, thereby creating different thermodynamic states within the fill region as the oxidizer emerges from the inlet.

Feleo, Alexander↗

DOE OSTI2025

A High-Order Discontinuous Galerkin Spectral Element Method for Compressible Reacting Flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large eddy simulations due to their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reactive Euler equations encountered in high-speed combustion. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of DG approach. Thus, the framework is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. The numerical method is implemented within the spectral element solver Nek5000. Validation cases are conducted for both non-reactive and reactive discontinuous flows to demonstrate the solver capability. In particular, canonical one-dimensional and two-dimensional detonation simulations are performed and the high-order numerical results are validated against available literature data.

computational fluid dynamics (CFD)↗