Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow physics”

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 919 records · Page 51

The mechanisms, dynamics, and implications of self-sealing and CO 2 resistance in wellbore cements

Here, this study analyzes the dynamics and mechanisms of the interactions of carbonated brine with hydrated-Portland-cement; in particular, the study focuses on self-sealing, a process whereby hydrated-Portland cement reacts with carbonated brine to for silica and calcium carbonate in sufficient quantities to seal the flow pathway. The analysis is based on a comprehensive set of reactive-transport simulations that explore the complex coupled dynamics between the fluid flow and mineral reactions that underlie self-sealing, and it relies heavily on the synthesis of the extensive body of work on wellbore integrity that has been conducted over the past decade. The analysis explores a large chemical and mineralogical diversity and a wide range in physical conditions and flow regimes, attempting to assess the robustness of the analysis. Self-sealing conditions arise over a wide range in cement properties and reservoir conditions. Although some properties and conditions promote a stronger self-sealing response, self-sealing occurs for a wide range of Ca:Si ratios in cement and for various reservoir fluid compositions. Self-sealing conditions move along a wellbore proportional to the flux of the leaking carbonated brine, and the reaction zone spreads out proportional to the fluid velocity, where volumetric flux and velocity are related by porosity (flux = velocity * porosity). However, self-sealing conditions can be maintained in a specific section of a wellbore by controlling the pressure drive and/or effective wellbore permeability, which in turn can limit the flux and velocity of any leaking fluid. Finally, the phases produced by hydrating Portland cement represent a carbonic cement that will react with a carbonated brine to produce end products (calcium carbonate and silica) that can maintain integrity in the presence of carbonic acid. Lastly, the attributes that make hydrated Portland cement phases a carbonic cement are required for self-sealing.

58 GEOSCIENCES↗

Validation of the space shuttle main engine steady state performance model

The primary objective was to present methods for validating predictions of Rocketdyne's most current version of the Space Shuttle Main Engine (SSME) Power Balance Model (PBM) with respect to physical relations governing flow systems. This required the development and implementation of postprocessors to check results of PBM computations for satisfaction of conservation relations. A cursory uncertainty analysis of PBM predictions with respect to mass and energy balances was performed. In addition, an effort to identify the empirical relations and physical assumptions within PBM which impact the ability of the model to attain rigorous balance was begun.

Santi, L. Michael↗

Joint state-parameter estimation for the reduced fracture model via the united filter

Here, in this paper, we introduce an effective United Filter method for jointly estimating the solution state and physical parameters in flow and transport problems within fractured porous media. Fluid flow and transport in fractured porous media are critical in subsurface hydrology, geophysics, and reservoir geomechanics. Reduced fracture models, which represent fractures as lower-dimensional interfaces, enable efficient multi-scale simulations. However, reduced fracture models also face accuracy challenges due to modeling errors and uncertainties in physical parameters such as permeability and fracture geometry. To address these challenges, we propose a United Filter method, which integrates the Ensemble Score Filter (EnSF) for state estimation with the Direct Filter for parameter estimation. EnSF, based on a score-based diffusion model framework, produces ensemble representations of the state distribution without deep learning. Meanwhile, the Direct Filter, a recursive Bayesian inference method, estimates parameters directly from state observations. The United Filter combines these methods iteratively: EnSF estimates are used to refine parameter values, which are then fed back to improve state estimation. Numerical experiments demonstrate that the United Filter method surpasses the state-of-the-art Augmented Ensemble Kalman Filter, delivering more accurate state and parameter estimation for reduced fracture models. This framework also provides a robust and efficient solution for PDE-constrained inverse problems with uncertainties and sparse observations.

Bayesian inference↗

Electron heat flux and propagating fronts in plasma thermal quench via ambipolar transport

The thermal collapse of a nearly collisionless plasma interacting with a cooling spot, in which the electron parallel heat flux plays an essential role, is both theoretically and numerically investigated. Here, we show that such thermal collapse, which is known as thermal quench in tokamaks, comes about in the form of propagating fronts, originating from the cooling spot, along magnetic field lines. The slow fronts, propagating with local ion sound speed, limit the aggressive cooling of plasma, which is accompanied by a plasma cooling flow toward the cooling spot. The extraordinary physics underlying such a cooling flow is that the fundamental constraint of ambipolar transport along the field line limits the spatial gradient of electron thermal conduction flux to the much weaker convective scaling, as opposed to the free-streaming scaling, so that a large electron temperature and, hence, pressure gradient can be sustained. The last ion front for a radiative cooling spot is a shock front where cold but flowing ions meet hot ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Turbine stator flow field simulations

The increased capability and accessibility of modern computers, coupled with increasingly sophisticated and accurate numerical and physical modeling, has led to a marked impact of numerical simulations upon current turbine design and research problems. The turbine section represents a considerable challenge as it contains significant regions of complex three-dimensional flow, including both aerodynamic and heat transfer phenomena. The focus of the present effort is the development of an efficient and accurate three-dimensional Navier-Stokes calculation procedure for application to the turbine stator and rotor problems. In particular, an effective procedure is sought which: (1) adequately represents the flow physics, (2) allows for sufficient resolution in regions of small length scale, and (3) has sufficiently good convergence properties so as to allow use on a regular basis.

Buggeln, R. C.↗

Laser velocimetry applied to transonic and supersonic aerodynamics

Measurements obtained with laser velocimetry in a Mach 2.9 separated turbulent boundary layer and in the transonic flow past a two-dimensional airfoil section are presented and compared to data realized by conventional techniques. Agreement in mean velocities was realized where the pressure measurements could be considered reliable; however, in regions of instantaneous reverse velocities, the laser results were found to be consistent with the physics of the flow whereas the pressure data were not. Streamwise turbulence intensities are also presented. In the transonic airfoil study, velocity measurements obtained immediately outside the upper surface boundary layer of a 6-inch chord NACA 64A010 airfoil are compared to edge velocities inferred from surface pressure measurements. For free-stream Mach numbers of 0.6 and 0.8, the agreement in results was very good. "Dual scatter" optical arrangements in conjunction with a single particle, counter-type signal processor were employed in these investigations.

Johnson, D. A.↗

Ionization and Deionization Processes in Low-Density Plasma Flows

Various plasma relaxation processes that are significant in laboratory plasma flows for aerospace physics studies are discussed on the basis of a survey of time decay studies of immobile plasmas. Some of the problems analyzed and discussed are: the relative importance of the various recombination mechanisms, catalytic deionization by electro- negative particles, delayed ionization produced by metastable action, and heating of the flow by deionization reactions.

Barger, Raymond L.↗

Development of a Navier-Stokes algorithm for parallel-processing supercomputers

An explicit flow solver, applicable to the hierarchy of model equations ranging from Euler to full Navier-Stokes, is combined with several techniques designed to reduce computational expense. The computational domain consists of local grid refinements embedded in a global coarse mesh, where the locations of these refinements are defined by the physics of the flow. Flow characteristics are also used to determine which set of model equations is appropriate for solution in each region, thereby reducing not only the number of grid points at which the solution must be obtained, but also the computational effort required to get that solution. Acceleration to steady-state is achieved by applying multigrid on each of the subgrids, regardless of the particular model equations being solved. Since each of these components is explicit, advantage can readily be taken of the vector- and parallel-processing capabilities of machines such as the Cray X-MP and Cray-2.

Swisshelm, Julie M.↗

Development of an aerodynamics algorithm for parallel-processing supercomputers

An explicit flow solver, applicable to the hierarchy of model equations ranging from Euler to full Navier-Stokes, is combined with several techniques designed to reduce computational expense. The computational domain consists of local grid refinements embedded in a global coarse mesh, where the locations of these refinements are defined by the physics of the flow. Flow characteristics are also used to determine which set of model equations is appropriate for solution in each region, thereby reducing not only the number of grid points at which the solution must be obtained, but also the computational effort required to get that solution. Acceleration to steady-state is achieved by applying multigrid on each of the subgrids, regardless of the particular model equations being solved. Since each of these components is explicit, advantage can readily be taken of the vector- and parallel-processing capabilities of machines such as the Cray X-MP and Cray-2.

Swisshelm, Julie M.↗

Gas Flow Detection System

This system provides a portable means to detect gas flow through a thin-walled tube without breaking into the tubing system. The flow detection system was specifically designed to detect flow through two parallel branches of a manifold with only one inlet and outlet, and is a means for verifying a space shuttle program requirement that saves time and reduces the risk of flight hardware damage compared to the current means of requirement verification. The prototype Purge Vent and Drain Window Cavity Conditioning System (PVD WCCS) Flow Detection System consists of a heater and a temperature-sensing thermistor attached to a piece of Velcro to be attached to each branch of a WCCS manifold for the duration of the requirement verification test. The heaters and thermistors are connected to a shielded cable and then to an electronics enclosure, which contains the power supplies, relays, and circuit board to provide power, signal conditioning, and control. The electronics enclosure is then connected to a commercial data acquisition box to provide analog to digital conversion as well as digital control. This data acquisition box is then connected to a commercial laptop running a custom application created using National Instruments LabVIEW. The operation of the PVD WCCS Flow Detection System consists of first attaching a heater/thermistor assembly to each of the two branches of one manifold while there is no flow through the manifold. Next, the software application running on the laptop is used to turn on the heaters and to monitor the manifold branch temperatures. When the system has reached thermal equilibrium, the software application s graphical user interface (GUI) will indicate that the branch temperatures are stable. The operator can then physically open the flow control valve to initiate the test flow of gaseous nitrogen (GN2) through the manifold. Next, the software user interface will be monitored for stable temperature indications when the system is again at thermal equilibrium with the test flow of GN2. The temperature drop of each branch from its "no flow" stable temperature peak to its stable "with flow" temperature will allow the operator to determine whether a minimum level of flow exists. An alternative operation has the operator turning on the software only long enough to record the ambient temperature of the tubing before turning on the heaters and initiating GN2 flow. The stable temperature of the heated tubing with GN2 flow is then compared with the ambient tubing temperature to determine if flow is present in each branch. To help quantify the level of flow in the manifolds, each branch will be bench calibrated to establish its thermal properties using the flow detection system and different flow rates. These calibration values can then be incorporated into the software application to provide more detailed flow rate information.

Moss, Thomas↗

Implementation of Radiation, Ablation, and Free Energy Minimization Modules for Coupled Simulations of Hypersonic Flow

A description of models and boundary conditions required for coupling radiation and ablation physics to a hypersonic flow simulation is provided. Chemical equilibrium routines for varying elemental mass fraction are required in the flow solver to integrate with the equilibrium chemistry assumption employed in the ablation models. The capability also enables an equilibrium catalytic wall boundary condition in the non-ablating case. The paper focuses on numerical implementation issues using FIRE II, Mars return, and Apollo 4 applications to provide context for discussion. Variable relaxation factors applied to the Jacobian elements of partial equilibrium relations required for convergence are defined. Challenges of strong radiation coupling in a shock capturing algorithm are addressed. Results are presented to show how the current suite of models responds to a wide variety of conditions involving coupled radiation and ablation.

Gnoffo, Peter A.↗

A Review of Lattice-Boltzmann Models Coupled with Geochemical Modeling Applied for Simulation of Advanced Waterflooding and Enhanced Oil Recovery Processes

To maintain economic profit and improve the oil production efficiency after the primary and secondary production phase, advanced waterflooding techniques such as low salinity waterflooding in carbonate reservoirs have been investigated in numerical simulations, laboratory experiments, and field pilot tests. Multiple underlying mechanisms have been proposed based on these studies, and they are still under debate. Various numerical modeling approaches are introduced, but there exists a lack of a pore-scale comprehensive modeling scheme to fully understand the processes. Lattice-Boltzmann method (LBM) is a type of numerical fluid flow modeling technique that shows capabilities and flexibilities in modeling pore-scale fluid flow to integrate physical–chemical processes within complex structures. The intrinsic feature of LBM makes it a promising framework for simulating advanced waterflooding due to its flexibility, accuracy, and parallel efficiency. LBM works either by itself for solving reactive transport problems or by coupling with a third-party reaction solver. This review mainly introduces the LBM fluid flow and reactive transport capabilities and the concept and modeling approaches to simulate advanced waterflooding techniques. Meanwhile, an evaluation of the coupled LBM models for enhanced oil recovery (EOR) simulations is discussed with future research challenges and directions concluded.

02 PETROLEUM↗

Discrete velocity computations with stochastic variance reduction of the Boltzmann equation for gas mixtures

We extend a variance reduced discrete velocity method developed at UT Austin [1, 2] to gas mixtures with large mass ratios and flows with trace species. The mixture is stored as a collection of independent velocity distribution functions, each with a unique grid in velocity space. Different collision types (A-A, A-B, B-B, etc.) are treated independently, and the variance reduction scheme is formulated with different equilibrium functions for each separate collision type. The individual treatment of species enables increased focus on species important to the physics of the flow, even if the important species are present in trace amounts. The method is verified through comparisons to Direct Simulation Monte Carlo computations and the computational workload per time step is investigated for the variance reduced method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

In Search of the Physics: The Interplay of Experiment and Computation in Airframe Noise Research: Flap-Edge Noise

The LaRC investigative process for airframe noise has proven to be a useful guide for elucidation of the physics of flow-induced noise generation over the last five years. This process, relying on a close interplay between experiment and computation, is described and demonstrated here on the archetypal problem of flap-edge noise. Some detailed results from both experiment and computation are shown to illustrate the process, and a description of the multi-source physics seen in this problem is conjectured.

Streett, C. L.↗

From Fluid Flow to Coupled Processes in Fractured Rock: Recent Advances and New Frontiers

Abstract Quantitative predictions of natural and induced phenomena in fractured rock is one of the great challenges in the Earth and Energy Sciences with far‐reaching economic and environmental impacts. Fractures occupy a very small volume of a subsurface formation but often dominate fluid flow, solute transport and mechanical deformation behavior. They play a central role in CO 2 sequestration, nuclear waste disposal, hydrogen storage, geothermal energy production, nuclear nonproliferation, and hydrocarbon extraction. These applications require predictions of fracture‐dependent quantities of interest such as CO 2 leakage rate, hydrocarbon production, radionuclide plume migration, and seismicity; to be useful, these predictions must account for uncertainty inherent in subsurface systems. Here, we review recent advances in fractured rock research covering field‐ and laboratory‐scale experimentation, numerical simulations, and uncertainty quantification. We discuss how these have greatly improved the fundamental understanding of fractures and one's ability to predict flow and transport in fractured systems. Dedicated field sites provide quantitative measurements of fracture flow that can be used to identify dominant coupled processes and to validate models. Laboratory‐scale experiments fill critical knowledge gaps by providing direct observations and measurements of fracture geometry and flow under controlled conditions that cannot be obtained in the field. Physics‐based simulation of flow and transport provide a bridge in understanding between controlled simple laboratory experiments and the massively complex field‐scale fracture systems. Finally, we review the use of machine learning‐based emulators to rapidly investigate different fracture property scenarios and accelerate physics‐based models by orders of magnitude to enable uncertainty quantification and near real‐time analysis.

58 GEOSCIENCES↗

Wickability-optimized textured liquid-desiccant air dehumidifiers for independent moisture management in energy-efficient buildings

Liquid-desiccant-based air conditioning systems are envisioned to enable independent humidity management, thereby improving the energy efficiency of future buildings. Existing liquid-desiccant-based air conditioning concepts, however, suffer from a poor liquid flow distribution deteriorating moisture removal rate. They are consequently flooded with the liquid-desiccant solution, which significantly degrades the energy efficiency of the dehumidification process. Here, in this work, the capillary forces and wickability effect of textured air dehumidifier surfaces are altered to minimize the liquid-desiccant flow rate of the fully wetted state, thereby transforming the physics of interfacial desiccant flow distribution. Consequently, the wickability-optimized air dehumidifier surface maximizes both moisture removal rate and dehumidification energy efficiency. It was interestingly found that the length scale of a textured air dehumidifier surface concept is optimized at an intermediate pattern density. Dry solid-air menisci appear at length scales exceeding the optimum pattern distance while the effective liquid–air interfacial area is reduced at smaller length scales, both of which degrade the moisture removal rate. At the optimum pattern density, the effective liquid–air interfacial area increases with the solution flow rate, thereby increasing the dehumidification rate. At a water vapor pressure potential of 3 kPa and a solution flow rate of 2.8 g/s, experimental results indicated a moisture removal rate of 0.16 g/m2-s for a textured surface concept with a capillary length scale of 3 mm, a 28% improvement compared with that of smooth-plate dehumidifier surfaces. A high moisture removal rate of the textured surface at a low desiccant flow rate led to a high thermal efficiency of 0.75 at a water vapor pressure potential of 5.6 kPa and a LiBr flow rate of 2.8 g/s. The insights gained from the present study accelerate the development of advanced textured surface concepts for next-generation liquid-desiccant-based air dehumidification systems offering independent humidity management for future energy-efficient buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Flow field characteristics of a 127-pin rod bundle with hexagonal spacer grids

A determination of nominal flow phenomena in liquid metal fast reactor (LMFR) fuel assemblies is critical toward generation-IV reactor development. Axially positioned spacer grids are used to maintain the geometry of hexagonal rod bundles and simultaneously introduce perturbations in the flow. Three-dimensional (3D) printed asymmetric honeycomb spacer grids were installed in a prototypical 127-pin LMFR fuel assembly model to study complex fluid dynamics interactions induced by the spacer grid and rods. To characterize flow dynamics in this intricate geometry, time-resolved particle image velocimetry (TR-PIV) using the matched-index-of-refraction method was employed to obtain non-intrusive velocity measurements for three axial planes (one near-wall and two interior planes) at a Reynolds number of 6000. The statistical TR-PIV results compared sub-channel-dependent normalized time-averaged velocity, velocity fluctuations, Reynolds stress, vorticity, and turbulence kinetic energy distributions. TR-PIV line profiles characterized downstream spacer grid flow dynamics. Two-point spatial and spatial–temporal cross-correlation fields revealed local coherent structures and quantified convection velocities of traveling vortices. Spatial–temporal decomposition using dynamic mode decomposition (DMD) applied to the near-wall vorticity fields extracted turbulent structures and flow instabilities in the wake region of the spacer grid, along with their decay and frequency rates. Reduced-order velocity fields from DMD reconstructions identified the most energy-containing coherent structures persistent in the near-wall region. This research provides experimental data sets and analyses of flow behavior in rod bundles with hexagonal spacer grids. In conclusion, the results are critical toward LMFR design and geometry optimization, crucial for the validation of computational fluid dynamics and reduced-order flow models.

42 ENGINEERING↗

Unified first-principles theory of thermal properties of insulators

The conventional first-principles theory for the thermal and thermodynamic properties of insulators is based on the perturbative treatment of the anharmonicity of crystal bonds. While this theory has been a successful predictive tool for strongly bonded solids such as diamond and silicon, here we show that it fails dramatically for strongly anharmonic (weakly bonded) materials, and that the conventional quasiparticle picture breaks down at relatively low temperatures. To address this failure, we present a unified first-principles theory of the thermodynamic and thermal properties of insulators that captures multiple thermal properties within the same framework across the full range of anharmonicity from strongly bonded to weakly bonded insulators. This theory features a new phonon renormalization approach derived from many-body physics that creates well-defined quasiparticles even at relatively high temperatures, and it accurately captures the effects of strongly anharmonic bonds on phonons and thermal transport. Using a prototypical strongly anharmonic material, sodium chloride (NaCl), as an example, we demonstrate that our new first-principles framework simultaneously captures the apparently contradictory experimental observations of large thermal expansion and low thermal conductivity of NaCl on the one hand, and anomalously weak temperature dependence of phonon modes on the other, while the conventional theory fails in all three cases. Here we demonstrate that four-phonon scattering due to higher-order anharmonicity significantly lowers the thermal conductivity of NaCl and is required for a proper comparison to experiment. Furthermore, we show that our renormalization framework, along with four-phonon scattering, also successfully predicts the measured phonon frequencies and thermal properties of a weakly anharmonic material, diamond, indicating universal applicability for thermal properties of insulators. Our work gives new insights into the physics of heat flow in solids, and presents a computationally efficient and rigorous framework that captures the thermal and thermodynamic properties of both weakly and strongly bonded insulators simultaneously.

36 MATERIALS SCIENCE↗