Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLUID MECHANICS”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Learning the boundary-to-domain mapping using Lifting Product Fourier Neural Operators for partial differential equations

Neural operators such as the Fourier Neural Operator (FNO) have been shown to provide resolution-independent deep learning models that can learn mappings between function spaces. For example, an initial condition can be mapped to the solution of a partial differential equation (PDE) at a future time-step using a neural operator. Despite the popularity of neural operators, their use to predict solution functions over a domain given only data over the boundary (such as a spatially varying Dirichlet boundary condition) remains unexplored. In this paper, we refer to such problems as boundary-to-domain problems; they have a wide range of applications in areas such as fluid mechanics, solid mechanics, heat transfer etc. We present a novel FNO-based architecture, named Lifting Product FNO (or LP-FNO) which can map arbitrary boundary functions defined on the lower-dimensional boundary to a solution in the entire domain. Specifically, two FNOs defined on the lower-dimensional boundary are lifted into the higher dimensional domain using our proposed lifting product layer. We demonstrate the efficacy and resolution independence of the proposed LP-FNO for the 2D Poisson equation.

Kashi, Aditya↗
DOE OSTI2016

High order ADER schemes for a unified first order hyperbolic formulation of continuum mechanics: Viscous heat-conducting fluids and elastic solids

Highlights: • High order schemes for a unified first order hyperbolic formulation of continuum mechanics. • The mathematical model applies simultaneously to fluid mechanics and solid mechanics. • Viscous fluids are treated in the frame of hyper-elasticity as generalized visco-plastic solids. • Formal asymptotic analysis reveals the connection with the Navier–Stokes equations. • The distortion tensor A in the model appears to be well-suited for flow visualization. - Abstract: This paper is concerned with the numerical solution of the unified first order hyperbolic formulation of continuum mechanics recently proposed by Peshkov and Romenski [110], further denoted as HPR model. In that framework, the viscous stresses are computed from the so-called distortion tensor A, which is one of the primary state variables in the proposed first order system. A very important key feature of the HPR model is its ability to describe at the same time the behavior of inviscid and viscous compressible Newtonian and non-Newtonian fluids with heat conduction, as well as the behavior of elastic and visco-plastic solids. Actually, the model treats viscous and inviscid fluids as generalized visco-plastic solids. This is achieved via a stiff source term that accounts for strain relaxation in the evolution equations of A. Also heat conduction is included via a first order hyperbolic system for the thermal impulse, from which the heat flux is computed. The governing PDE system is hyperbolic and fully consistent with the first and the second principle of thermodynamics. It is also fundamentally different from first order Maxwell–Cattaneo-type relaxation models based on extended irreversible thermodynamics. The HPR model represents therefore a novel and unified description of continuum mechanics, which applies at the same time to fluid mechanics and solid mechanics. In this paper, the direct connection between the HPR model and the classical hyperbolic–parabolic Navier–Stokes–Fourier theory is established for the first time via a formal asymptotic analysis in the stiff relaxation limit. From a numerical point of view, the governing partial differential equations are very challenging, since they form a large nonlinear hyperbolic PDE system that includes stiff source terms and non-conservative products. We apply the successful family of one-step ADER–WENO finite volume (FV) and ADER discontinuous Galerkin (DG) finite element schemes to the HPR model in the stiff relaxation limit, and compare the numerical results with exact or numerical reference solutions obtained for the Euler and Navier–Stokes equations. Numerical convergence results are also provided. To show the universality of the HPR model, the paper is rounded-off with an application to wave propagation in elastic solids, for which one only needs to switch off the strain relaxation source term in the governing PDE system. We provide various examples showing that for the purpose of flow visualization, the distortion tensor A seems to be particularly useful.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Computational thermal, chemical, fluid, and solid mechanics for geosystems management.

This document summarizes research performed under the SNL LDRD entitled - Computational Mechanics for Geosystems Management to Support the Energy and Natural Resources Mission. The main accomplishment was development of a foundational SNL capability for computational thermal, chemical, fluid, and solid mechanics analysis of geosystems. The code was developed within the SNL Sierra software system. This report summarizes the capabilities of the simulation code and the supporting research and development conducted under this LDRD. The main goal of this project was the development of a foundational capability for coupled thermal, hydrological, mechanical, chemical (THMC) simulation of heterogeneous geosystems utilizing massively parallel processing. To solve these complex issues, this project integrated research in numerical mathematics and algorithms for chemically reactive multiphase systems with computer science research in adaptive coupled solution control and framework architecture. This report summarizes and demonstrates the capabilities that were developed together with the supporting research underlying the models. Key accomplishments are: (1) General capability for modeling nonisothermal, multiphase, multicomponent flow in heterogeneous porous geologic materials; (2) General capability to model multiphase reactive transport of species in heterogeneous porous media; (3) Constitutive models for describing real, general geomaterials under multiphase conditions utilizing laboratory data; (4) General capability to couple nonisothermal reactive flow with geomechanics (THMC); (5) Phase behavior thermodynamics for the CO2-H2O-NaCl system. General implementation enables modeling of other fluid mixtures. Adaptive look-up tables enable thermodynamic capability to other simulators; (6) Capability for statistical modeling of heterogeneity in geologic materials; and (7) Simulator utilizes unstructured grids on parallel processing computers.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2019

Fluid-based gating mechanism with tunable multiphase selectivity and antifouling behavior

A gating mechanism that uses a capillary stabilized liquid as a reversible, reconfigurable gate that fills and seals pores in the closed state, and creates a non-fouling, liquid-lined pore in the open state is disclosed. Theoretical modeling and experiments demonstrate that for each transport substance, the gating threshold—the pressure needed to open the pores—can be rationally tuned over a wide pressure range. This enables realizing in one system differential response profiles for a variety of liquids and gases, even letting liquids flow through the pore while preventing gas from escaping. These capabilities allow dynamic modulation of gas-liquid sorting and to separate multi-phase mixtures, with the liquid lining ensuring sustained antifouling behavior. Because the liquid gating strategy enables efficient short-term and long-term operation and can be applied to a variety of pore structures and membrane materials, and to nano, micro as well as macroscale fluid systems, the gating systems is useful in a wide range of applications.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Multiscale Simulation Framework for Coupled Fluid Flow and Mechanical Deformation

A multiscale linear-solver framework for the pressure equation associated with flow in highly heterogeneous porous formations was developed. The multiscale based approach is cast in a general algebraic form, which facilitates integration of the new scalable linear solver in existing flow simulators. The Algebraic Multiscale Solver (AMS) is employed as a preconditioner within a multi-stage strategy. The formulations investigated include the standard MultiScale Finite-Element (MSFE) andMultiScale Finite-Volume (MSFV) methods. The local-stage solvers include incomplete factorization and the so-called Correction Functions (CF) associated with the MSFV approach. Extensive testing of AMS, as an iterative linear solver, indicate excellent convergence rates and computational scalability. AMS compares favorably with advanced Algebraic MultiGrid (AMG) solvers for highly detailed three-dimensional heterogeneous models. Moreover, AMS is expected to be especially beneficial in solving time-dependent problems of coupled multiphase flow and transport in large-scale subsurface formations.

02 PETROLEUM↗
DOE OSTI2016

Multiscale Simulation Framework for Coupled Fluid Flow and Mechanical Deformation

Our work in this project is aimed at making fundamental advances in multiscale methods for flow and transport in highly heterogeneous porous media. The main thrust of this research is to develop a systematic multiscale analysis and efficient coarse-scale models that can capture global effects and extend existing multiscale approaches to problems with additional physics and uncertainties. A key emphasis is on problems without an apparent scale separation. Multiscale solution methods are currently under active investigation for the simulation of subsurface flow in heterogeneous formations. These procedures capture the effects of fine-scale permeability variations through the calculation of specialized coarse-scale basis functions. Most of the multiscale techniques presented to date employ localization approximations in the calculation of these basis functions. For some highly correlated (e.g., channelized) formations, however, global effects are important and these may need to be incorporated into the multiscale basis functions. Other challenging issues facing multiscale simulations are the extension of existing multiscale techniques to problems with additional physics, such as compressibility, capillary effects, etc. In our project, we explore the improvement of multiscale methods through the incorporation of additional (single-phase flow) information and the development of a general multiscale framework for flows in the presence of uncertainties, compressible flow and heterogeneous transport, and geomechanics. We have considered (1) adaptive local-global multiscale methods, (2) multiscale methods for the transport equation, (3) operator-based multiscale methods and solvers, (4) multiscale methods in the presence of uncertainties and applications, (5) multiscale finite element methods for high contrast porous media and their generalizations, and (6) multiscale methods for geomechanics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A High-Fidelity Approach for the Simulation of Flow-Induced Vibration

Flow-induced vibration (FIV) is a widespread problem in energy systems becausethey rely on fluid movement for energy conversion. Vibrating structures may be damaged as fatigue or wear occur. Given the importance of reliable components in the nuclear industry, flow-induced vibrations have long been a major concern in the safety and operation of nuclear reactors. In particular,nuclear fuel and steam generators have been known to suffer from flow-induced vibrations and related failures. Over the past five years, the Nuclear EnergyAdvanced Modeling and Simulation program has developed the integrated multiphysics code suite SHARP. The goal of developing such a tool is to perform multiphysics modeling of the components inside a reactor core, the full reactor core or portions of it, and be able to achieve that with various levelsof fidelity. This flexibility allows users to select the appropriate level of fidelity for their computational resources and design constraints.In particular SHARP contains high-fidelity single-physics codes for structural mechanics and fluid mechanics calculations: the structural mechanics implicit code Diablo and the computational fluid dynamics spectral element code Nek5000. Both codes are state-of-the-art.highly scalable (up to millions of processors in the case of Nek5000) tools that have been extensively validated. These tools form a strong basis on which to build an FIV modeling capability.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Shock tube experiments on the three-layer Richtmyer–Meshkov instability

Here, a vertical shock tube is used for experiments on the three-layer Richtmyer–Meshkov instability. Two closely spaced membrane-less interfaces are formed by the flow of two different sects of three gases: one with air above CO 2 above SF 6 and the other with helium above air above SF 6 . The lightest of the three gases enters the shock tube at the top of the driven section and flows downward. Conversely, the heaviest gas enters at the bottom of the shock tube and flows upward while the intermediate density gas enters at the middle through porous plates. All three gases are allowed to escape through holes at the layer location, leaving an approximately 30-mm layer of intermediate-density gas suspended between the lightest gas from above and the heaviest gas from below. A single-mode, two-dimensional initial perturbation is then imposed on the lower interface by oscillating the shock tube in the horizontal direction. The flow is visualized by seeding the intermediate gas with particles and illuminating it with a pulsed laser. Image sequences are then captured using high-speed video cameras. Perturbation amplitude measurements are made from the three-layer system and compared with measurements from 2, two-layer systems. It is observed that the presence of the upper, initially flat interface produces a decrease in growth of instability amplitude in the nonlinear phase over an equivalent single-interface configuration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Effect of pore fluid chemistry on the mechanical behavior of a divalent compacted bentonite, an experimental and constitutive study

Ongoing research in isolating high-level nuclear waste and spent fuel has highlighted compacted bentonite as a suitable material for engineered barrier systems in deep geological repositories due to its extraordinary swelling and retention properties. This research focuses on the chemo-mechanical behavior of compacted bentonite exposed to different pore fluids with different concentrations and loading conditions. The study involves swelling pressure and compressibility experiments along with mineralogy analysis employing X-ray diffraction (XRD) and Cation exchange. The tests were conducted on BCV (a Mg/Ca- bentonite) compacted at a dry density of 1.48 ± .02 Mg/m 3 . An advanced chemical-mechanical constitutive model for unsaturated highly expansive clays was adopted to simulate the material response and better understand its behavior. The model is able to account for the main phenomena at both macro and microstructural levels and the interactions between them. The model successfully replicated experimental observations. The XRD analyses support the macroscopic observation, indicating that salinity impacts crystalline swelling as demonstrated by the reduction of basal spacing from 19.27 Å to 15.68 Å when the osmotic suction increases from 0 MPa to 33 MPa. The results suggested that the osmotic pressure generated by the concentration in the pore fluids promotes a reduction in swelling pressures, swelling strains, and crystalline swelling of clay minerals. Also, it affects the pre-consolidation stress and the compressibility of the compacted samples. In conclusion, it was also observed that both solution type and solution concentration impact the clay swelling pressure.

Chemo-mechanical constitutive model↗
DOE OSTI2012

THERMO-HYDRO-MECHANICAL MODELING OF WORKING FLUID INJECTION AND THERMAL ENERGY EXTRACTION IN EGS FRACTURES AND ROCK MATRIX

Development of enhanced geothermal systems (EGS) will require creation of a reservoir of sufficient volume to enable commercial-scale heat transfer from the reservoir rocks to the working fluid. A key assumption associated with reservoir creation/stimulation is that sufficient rock volumes can be hydraulically fractured via both tensile and shear failure, and more importantly by reactivation of naturally existing fractures (by shearing), to create the reservoir. The advancement of EGS greatly depends on our understanding of the dynamics of the intimately coupled rock-fracture-fluid-heat system and our ability to reliably predict how reservoirs behave under stimulation and production. Reliable performance predictions of EGS reservoirs require accurate and robust modeling for strongly coupled thermal-hydrological-mechanical (THM) processes. Conventionally, these types of problems have been solved using operator-splitting methods, usually by coupling a subsurface flow and heat transport simulators with a solid mechanics simulator via input files. An alternative approach is to solve the system of nonlinear partial differential equations that govern multiphase fluid flow, heat transport, and rock mechanics simultaneously, using a fully coupled, fully implicit solution procedure, in which all solution variables (pressure, enthalpy, and rock displacement fields) are solved simultaneously. This paper describes numerical simulations used to investigate the poro- and thermal- elastic effects of working fluid injection and thermal energy extraction on the properties of the fractures and rock matrix of a hypothetical EGS reservoir, using a novel simulation software FALCON (Podgorney et al., 2011), a finite element based simulator solving fully coupled multiphase fluid flow, heat transport, rock deformation, and fracturing using a global implicit approach. Investigations are also conducted on how these poro- and thermal-elastic effects are related to fracture permeability evolution.

15 GEOTHERMAL ENERGY↗
DOE OSTI2015

Geomechanical and Fluid Transport Properties

Complex changes to fluid and mechanical properties of a pyrolyzing formation could lead to significant subsidence. Capturing subsidence requires accounting for changes in porosity due to generation and drainage of fluid. Pores, permeability and pore pressure are generated during pyrolysis of kerogen into oil and gas. In a series reaction scheme, oil is an intermediate product of pyrolysis. A large portion of oil-range fluids may be lost in an in-situ operation if the oil produced cannot be drained quickly enough. Petroleum fluid quality and production are highly dependent on the relative rates of creation, destruction, and transport of oil and gas in a given volume. All of these are dependent on the heating rate of oil shale. If oil shale is heated quickly, fluid flow may increase, but at a cost to the amount of oil that is finally extracted. To account for the effect of pore fluid on mechanical deformation, an elastic model was combined with a Mohr-Coulomb failure criterion and Biot’s linear consolidation model. The overall model is poroelastic, because it accounts for rock deforming by considering the response of the bulk system and the pore fluids it contains. Biot’s linear model includes all ways that pore fluid volume can change in a system, including stress-strain interactions between pores and the bulk material. Additionally, the model captures the distinction between fluid and solid matrix compression and describes how bulk deformation is affected by fluid drainage. The fluid response to changing stress is affected by changes in pore pressure, while the bulk material responds to effective stress changes. To capture all of these effects, a number of theoretical and empirical models were employed in FLAC3D™, a rock-mechanics software package produced by Itasca Consulting Group. The basic capabilities of FLAC3D™ have been supplemented by customized algorithms for the simulation package. The most important contribution has been to develop a methodology for describing poro-thermo-mechanical properties of oil shale across the complete chronology of in situ pyrolysis. This report highlights those methods and results for inferring realistic oil shale properties.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2016

Resolved-particle simulation by the Physalis method: Enhancements and new capabilities

We present enhancements and new capabilities of the Physalis method for simulating disperse multiphase flows using particle-resolved simulation. The current work enhances the previous method by incorporating a new type of pressure-Poisson solver that couples with a new Physalis particle pressure boundary condition scheme and a new particle interior treatment to significantly improve overall numerical efficiency. Further, we implement a more efficient method of calculating the Physalis scalar products and incorporate short-range particle interaction models. We provide validation and benchmarking for the Physalis method against experiments of a sedimenting particle and of normal wall collisions. We conclude with an illustrative simulation of 2048 particles sedimenting in a duct. In the appendix, we present a complete and self-consistent description of the analytical development and numerical methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

On plasma detachment in propulsive magnetic nozzles

Three detachment mechanisms proposed in the literature (via resistivity, via electron inertia, and via induced magnetic field) are analyzed with an axisymmetric model of the expansion of a small-beta, weakly collisional, near-sonic plasma in a diverging magnetic nozzle. The model assumes cold, partially magnetized ions and hot, isothermal, fully magnetized electrons. Different conditions of the plasma beam at the nozzle throat are considered. A central feature is that a positive thrust gain in the nozzle of a plasma thruster is intimately related to the azimuthal current in the plasma being diamagnetic. Then, and contrary to existing expectations, the three aforementioned detachment mechanisms are divergent, that is, the plasma beam diverges outwards of the guide nozzle, further hindering its axial expansion and the thrust efficiency. The rate of divergent detachment is quantified for the small-parameter range of the three mechanisms. Alternative mechanisms for a convergent detachment of the plasma beam are suggested.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Modeling of shock-induced force on an isolated particle in water and air

The prediction of force on an isolated particle, while a shock is passing over it, is an important problem in many natural and industrial applications. Although the flow monotonically changes from the pre-shock to the post-shock state, the particle's force has been observed to behave nonmonotonically with a sharp peak when the shock is located halfway across the particle. This nonmonotonic behavior is due to the unsteady nature of the compression and rarefaction waves that radiate as the shock diffracts around the particle and, therefore, cannot be predicted by a quasi-steady model. An accurate force model must account for the unsteady nature of the flow and the sharp discontinues in the flow properties across the shock. In this work, we test four different inviscid models and observe that the compressible Maxey–Riley–Gatignol (C-MRG) model is the most accurate based on comparison with results from particle-resolved inviscid simulations at two different Mach numbers for both water and air as the medium. The C-MRG model is first demonstrated to predict the force on a stationary particle accurately and then extended to capture the force on a moving particle. Numerical complexities regarding the implementation of the C-MRG model are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗