Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLUID FLOW”

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 1,459 records · Page 81

Modeling Nanoconfinement Effects Using Active Learning

Predicting the spatial configuration of gas in nanopores of is relevant in applications such as fluid flow forecasting and hydrocarbon reserves estimation. For example, shale reservoirs have suffered from computationally intractable multiscale problems, since fluid properties such as viscosity, density, and adsorption must be calculated by using expensive molecular dynamics (MD) simulations within each nanopore, whereas flow through these connected nanopores must be simulated at the micrometer scale. We utilize machine learning techniques to quickly and accurately model nanoscale confinement effects as an important step toward bridging the nano and micro scales. Our workflow is based on building and training physics-based deep-neural-networks models by learning from a database of MD calculations. The model accounts for the adsorption phenomenon by predicting the statistical distribution of gas inside nanopores. Because large databases of MD calculations are expensive to create, we investigate active learning (AL) as a data set construction strategy. In this workflow, new data are selected based on the model uncertainty via the query-by-committee approach. We show that our workflow obtains accurate models that generalize to real scanning electron microscopy geometries with 1/10th of the number of MD calculations required vs random data set generation. Our method enables the possibility of modeling nanoconfinement effects at the mesoscale, where complex connected sets of nanopores affect flow.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On-Shore Central Hydraulic Power Generation for Wind and Tidal Energy

Tidal energy, offshore wind energy, and onshore wind energy can be converted to electricity at a central ground location by means of converting their respective energies into high-pressure hydraulic flows that are transmitted to a system of generators by high-pressure pipelines. The high-pressure flows are then efficiently converted to electricity by a central power plant, and the low-pressure outlet flow is returned. The Department of Energy (DOE) is presently supporting a project led by Sunlight Photonics to demonstrate a 15 kW tidal hydraulic power generation system in the laboratory and possibly later submerged in the ocean. All gears and submerged electronics are completely eliminated. A second portion of this DOE project involves sizing and costing a 15 MW tidal energy system for a commercial tidal energy plant. For this task, Atlantis Resources Corporation s 18-m diameter demonstrated tidal blades are rated to operate in a nominal 2.6 m/sec tidal flow to produce approximately one MW per set of tidal blades. Fifteen units would be submerged in a deep tidal area, such as in Maine s Western Passage. All would be connected to a high-pressure (20 MPa, 2900 psi) line that is 35 cm ID. The high-pressure HEPG fluid flow is transported 500-m to on-shore hydraulic generators. HEPG is an environmentally-friendly, biodegradable, watermiscible fluid. Hydraulic adaptations to ORPC s cross-flow turbines are also discussed. For 15 MW of wind energy that is onshore or offshore, a gearless, high efficiency, radial piston pump can replace each set of top-mounted gear-generators. The fluid is then pumped to a central, easily serviceable generator location. Total hydraulic/electrical efficiency is 0.81 at full rated wind or tidal velocities and increases to 0.86 at 1/3 rated velocities.

wind energy↗

Preliminary Design Study for an Eddy Current Flow Meter for Versatile Test Reactor

The eddy current flow meter (ECFM) serves the purpose of monitoring velocities and temperatures of conductive liquid coolant surrounding nuclear reactor cells. The ECFM can also detect voids in two-phase liquid metal flows. Measurement of this flow in reactor cells can greatly improve the accuracy of the empirical formulas used for fluid dynamics near the reactor core, as well as, aid in mapping any non-uniformities in fluid flow, its temperature estimation and phase. With modeling and experimental verification, an extensible ECFM for multiple facilities and environments can lead to an understanding of the liquid coolant flow; while saving cost of calibration and development. This can be done by characterizing the ECFM performance based on its parameters and then optimizing a sensor for a given use case. Simulation results using an ECFM and an aluminum rod medium are presented and discussed herein. Additionally, liquid metal simulations for Mercury and Galinstan will compare the ECFM behavior for different liquid velocities and excitation frequencies. After sensor calibration curves are generated and experimentally verified, ECFMs can be developed and deployed using the model introduced.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The Acoustic Air-Jet Generator

This paper describes the development history, operation, and working principle of an acoustic air-jet generator, which can be used to produce high intensity, high frequency acoustic waves. Flow imaging experiments of the system are shown, and applications of the system for various fluid flows are presented. While previous articles on this system have been written in English by the author, this article provides the only documented explanation and implementation of the focusing schlieren system used for flow imaging of the system. This paper has been cited in early works on focusing schlieren as one of the first implementations of such a system, but has not been widely cited in the bulk of the literature on the focusing schlieren method (aside from a more recent M.S. thesis), perhaps due to the difficulty of accessing the paper or that no translation was available. The purpose of this translation is to provide an English version of the article for researchers interested in this work of Dr. Hartmann, and in particular the description of his focusing schlieren system and its application.

Julius Hartmann↗

Fully quantum algorithm for mesoscale fluid simulations with application to partial differential equations

Fluid flow simulations marshal our most powerful computational resources. In many cases, even this is not enough. Quantum computers provide an opportunity to speed up traditional algorithms for flow simulations. We show that lattice-based mesoscale numerical methods can be executed as efficient quantum algorithms due to their statistical features. This approach revises a quantum algorithm for lattice gas automata to reduce classical computations and state preparation at every time step. For this, the algorithm approximates the qubit relative phases and subtracts them at the end of each time step. Phases are evaluated using the iterative phase estimation algorithm and subtracted using single-qubit rotation phase gates. Further, this method optimizes the quantum resource required and makes it more appropriate for near-term quantum hardware. We also demonstrate how the checkerboard deficiency that the D1Q2 scheme presents can be resolved using the D1Q3 scheme. The algorithm is validated by simulating two canonical partial differential equations: the diffusion and Burgers' equations on different quantum simulators. We find good agreement between quantum simulations and classical solutions for the presented algorithm.

97 MATHEMATICS AND COMPUTING↗

Computational Analysis of the Impact of Boundary Conditions on a Particle-Laden Flow: A Case Study in a Pressurized Oxy-Coal Combustor

Designing an effective burner is vital for the development of coal combustion technologies. Because of high pressure, the volumetric fraction of the coal particles in the injected fuel in a pressurized oxy-combustion (POC) burner approaches or even exceeds the limitations allowed by the commercial computational fluid dynamics codes (e.g., Ansys Fluent). Consequently, for such high particle volumetric fractions, the interplay between the particles, the fluid flow, and the burner wall needs to be re-evaluated. The present computational work is a first step in a systematic analysis of the roles of various characteristics involved in the POC process, such as the method of particle release, its location, and the particle size. Specifically, pulverized coal is burned under an elevated pressure of 15 bar in an O 2 /CO 2 environment. A 100 kW, a POC combustor, is modeled with Ansys Fluent using the Reynolds-averaged Navier−Stokes approach. It is revealed that for this pilot-scale, pressurized burner, the gas phase flow velocity in the near-wall region exhibits anomalies. With the major focus on POC, this work aims to eliminate/reduce the impact of high particle loading on the gas-phase flow. To scrutinize the role of particle loading in the near-wall region and eliminate the impact of this velocity on POC downstream, the particle−gas interplay in the boundary layer is investigated by means of the computational simulations incorporating the coupling between the turbulent flow and the particles. It is found that the tuning of the particle release location makes the gas-phase flow velocity in the presence of particles consistent with the pure gas flow velocity profile. The particles size is also found to have a significant impact on the particle trajectory.

CFD↗

User's guide of TOUGH2-EGS-MP: A Massively Parallel Simulator with Coupled Geomechanics for Fluid and Heat Flow in Enhanced Geothermal Systems VERSION 1.0

TOUGH2-EGS-MP is a parallel numerical simulation program coupling geomechanics with fluid and heat flow in fractured and porous media, and is applicable for simulation of enhanced geothermal systems (EGS). TOUGH2-EGS-MP is based on the TOUGH2-MP code, the massively parallel version of TOUGH2. In TOUGH2-EGS-MP, the fully-coupled flow-geomechanics model is developed from linear elastic theory for thermo-poro-elastic systems and is formulated in terms of mean normal stress as well as pore pressure and temperature. Reservoir rock properties such as porosity and permeability depend on rock deformation, and the relationships between these two, obtained from poro-elasticity theories and empirical correlations, are incorporated into the simulation. This report provides the user with detailed information on the TOUGH2-EGS-MP mathematical model and instructions for using it for Thermal-Hydrological-Mechanical (THM) simulations. The mathematical model includes the fluid and heat flow equations, geomechanical equation, and discretization of those equations. In addition, the parallel aspects of the code, such as domain partitioning and communication between processors, are also included. Although TOUGH2-EGS-MP has the capability for simulating fluid and heat flows coupled with geomechanical effects, it is up to the user to select the specific coupling process, such as THM or only TH, in a simulation. There are several example problems illustrating applications of this program. These example problems are described in detail and their input data are presented. Their results demonstrate that this program can be used for field-scale geothermal reservoir simulation in porous and fractured media with fluid and heat flow coupled with geomechanical effects.

15 GEOTHERMAL ENERGY↗

User's Guide of TOUGH2-EGS. A Coupled Geomechanical and Reactive Geochemical Simulator for Fluid and Heat Flow in Enhanced Geothermal Systems Version 1.0

TOUGH2-EGS is a numerical simulation program coupling geomechanics and chemical reactions for fluid and heat flows in porous media and fractured reservoirs of enhanced geothermal systems. The simulator includes the fully-coupled geomechanical (THM) module, the fully-coupled geochemical (THC) module, and the sequentially coupled reactive geochemistry (THMC) module. The fully-coupled flow-geomechanics model is developed from the linear elastic theory for the thermo-poro-elastic system and is formulated with the mean normal stress as well as pore pressure and temperature. The chemical reaction is sequentially coupled after solution of flow equations, which provides the flow velocity and phase saturation for the solute transport calculation at each time step. In addition, reservoir rock properties, such as porosity and permeability, are subjected to change due to rock deformation and chemical reactions. The relationships between rock properties and geomechanical and chemical effects from poro-elasticity theories and empirical correlations are incorporated into the simulator. This report provides the user with detailed information on both mathematical models and instructions for using TOUGH2-EGS for THM, THC or THMC simulations. The mathematical models include the fluid and heat flow equations, geomechanical equation, reactive geochemistry equations, and discretization methods. Although TOUGH2-EGS has the capability for simulating fluid and heat flows coupled with both geomechanical and chemical effects, it is up to the users to select the specific coupling process, such as THM, THC, or THMC in a simulation. There are several example problems illustrating the applications of this program. These example problems are described in details and their input data are presented. The results demonstrate that this program can be used for field-scale geothermal reservoir simulation with fluid and heat flow, geomechanical effect, and chemical reaction in porous and fractured media.

15 GEOTHERMAL ENERGY↗

Using Magnetic Fields to Control Convection during Protein Crystallization: Analysis and Validation Studies

The effect of convection during the crystallization of proteins is not very well understood. In a gravitational field, convection is caused by crystal sedimentation and by solutal buoyancy induced flow and these can lead to crystal imperfections. While crystallization in microgravity can approach diffusion limited growth conditions (no convection), terrestrially strong magnetic fields can be used to control fluid flow and sedimentation effects. In this work, we develop the analysis for magnetic flow control and test the predictions using analog experiments. Specifically, experiments on solutal convection in a paramagnetic fluid were conducted in a strong magnetic field gradient using a dilute solution of Manganese Chloride. The observed flows indicate that the magnetic field can completely counter the settling effects of gravity locally and are consistent with the theoretical predictions presented. This phenomenon suggests that magnetic fields may be useful in mimicking the microgravity environment of space for some crystal growth ana biological applications where fluid convection is undesirable.

Ramachandran, N.↗

Acoustic transmission matrix of a variable area duct or nozzle carrying a compressible subsonic flow

The differential equations governing the propagation of sound in a variable area duct or nozzle carrying a one dimensional subsonic compressible fluid flow are derived and put in state variable form using acoustic pressure and particle velocity as the state variables. The duct or nozzle is divided into a number of regions. The region size is selected so that in each region the Mach number can be assumed constant and the area variation can be approximated by an exponential area variation. Consequently, the state variable equation in each region has constant coefficients. The transmission matrix for each region is obtained by solving the constant coefficient acoustic state variable differential equation. The transmission matrix for the duct or nozzle is the product of the individual transmission matrices of each region. Solutions are presented for several geometries with and without mean flow.

Miles, J. H.↗

Thermocapillary flow and melt/solid interfaces in floating-zone crystal growth under microgravity

Computer simulation of steady-state axisymmetrical heat transfer and fluid flow was conducted to study thermocapillary flow and melt/solid interfaces in floating-zone crystal growth under microgravity. The effects of key variables on the extent of thermocapillary flow in the melt zone, the shapes of melt/solid interfaces and the length of the melt zone were discussed. These variables are: (1) the temperature coefficient of surface tension (or the Marangoni number), (2) the pulling speed (or the Peclet number), (3) the feed rod radius, (4) the ambient temperature distribution, (5) the heat transfer coefficient (or the Biot number), and (6) the thermal diffusivity of the material (or the Prandtl number).

Lan, C. W.↗

Algorithms for parallel flow solvers on message passing architectures

The purpose of this project has been to identify and test suitable technologies for implementation of fluid flow solvers -- possibly coupled with structures and heat equation solvers -- on MIMD parallel computers. In the course of this investigation much attention has been paid to efficient domain decomposition strategies for ADI-type algorithms. Multi-partitioning derives its efficiency from the assignment of several blocks of grid points to each processor in the parallel computer. A coarse-grain parallelism is obtained, and a near-perfect load balance results. In uni-partitioning every processor receives responsibility for exactly one block of grid points instead of several. This necessitates fine-grain pipelined program execution in order to obtain a reasonable load balance. Although fine-grain parallelism is less desirable on many systems, especially high-latency networks of workstations, uni-partition methods are still in wide use in production codes for flow problems. Consequently, it remains important to achieve good efficiency with this technique that has essentially been superseded by multi-partitioning for parallel ADI-type algorithms. Another reason for the concentration on improving the performance of pipeline methods is their applicability in other types of flow solver kernels with stronger implied data dependence. Analytical expressions can be derived for the size of the dynamic load imbalance incurred in traditional pipelines. From these it can be determined what is the optimal first-processor retardation that leads to the shortest total completion time for the pipeline process. Theoretical predictions of pipeline performance with and without optimization match experimental observations on the iPSC/860 very well. Analysis of pipeline performance also highlights the effect of uncareful grid partitioning in flow solvers that employ pipeline algorithms. If grid blocks at boundaries are not at least as large in the wall-normal direction as those immediately adjacent to them, then the first processor in the pipeline will receive a computational load that is less than that of subsequent processors, magnifying the pipeline slowdown effect. Extra compensation is needed for grid boundary effects, even if all grid blocks are equally sized.

Vanderwijngaart, Rob F.↗

On-Shore Central Hydraulic Power Generation for Wind and Tidal Energy

Tidal energy, offshore wind energy, and onshore wind energy can be converted to electricity at a central ground location by means of converting their respective energies into high-pressure hydraulic flows that are transmitted to a system of generators by high-pressure pipelines. The high-pressure flows are then efficiently converted to electricity by a central power plant, and the low-pressure outlet flow is returned. The Department of Energy (DOE) is presently supporting a project led by Sunlight Photonics to demonstrate a 15 kilowatt tidal hydraulic power generation system in the laboratory and possibly later submerged in the ocean. All gears and submerged electronics are completely eliminated.A second portion of this DOE project involves sizing and costing a 15 megawatt tidal energy system for a commercial tidal energy plant. For this task, Atlantis Resources Corporation's 18-m diameter demonstrated tidal blades are rated to operate in a nominal 2.6 m/sec tidal flow to produce approximately one megawatt per set of tidal blades. Fifteen units would be submerged in a deep tidal area, such as in Maine's Western Passage. All would be connected to a high-pressure (20 megapascals, 2900 pounds per square inch) line that is 35 cm ID. The high-pressure HEPG fluid flow is transported 500-m to on-shore hydraulic generators. HEPG is an environmentally-friendly, biodegradable, water-miscible fluid. Hydraulic adaptations to ORPC's cross-flow turbines are also discussed.For 15 megawatt of wind energy that is onshore or offshore, a gearless, high efficiency, radial piston pump can replace each set of top-mounted gear-generators. The fluid is then pumped to a central, easily serviceable generator location. Total hydraulic/electrical efficiency is 0.81 at full rated wind or tidal velocities and increases to 0.86 at 1/3 rated velocities.

wind energy↗

Simulating many-engine spacecraft: Exceeding 1 quadrillion degrees of freedom via information geometric regularization

We present an optimized implementation of the recently proposed information geometric regularization (IGR) for unprecedented scale simulation of compressible fluid flows applied to multi-engine spacecraft boosters. We improve upon state-of-the-art computational fluid dynamics (CFD) techniques in terms of computational cost, memory footprint, and energy-to-solution metrics. Unified memory on coupled CPU–GPU or APU platforms increases problem size with negligible overhead. Mixed half/single-precision storage and computation are used on well-conditioned numerics. We simulate flow at 200 trillion grid points and 1 quadrillion degrees of freedom, exceeding the current record by a factor of 20. A factor of 4 wall-time speedup is achieved over optimized baselines. Ideal weak scaling is observed on OLCF Frontier, LLNL El Capitan, and CSCS Alps using the full systems. Strong scaling is near ideal at extreme conditions, including 80% efficiency on CSCS Alps with an 8 node baseline and stretching to the full system.

Wilfong, Benjamin [Georgia Institute of Technology↗

Thermal Analysis of Magnetically-Coupled Pump for Cryogenic Applications

Magnetically-coupled pump is under evaluation at Kennedy Space Center for possible cryogenic applications. A major concern is the impact of low temperature fluid flows on the pump performance. As a first step toward addressing this and related issues, a computational fluid dynamics and heat transfer tool has been adopted in a pump geometry. The computational tool includes (i) a commercial grid generator to handle multiple grid blocks and complicated geometric definitions, and (ii) an in-house computational fluid dynamics and heat transfer software developed in the Principal Investigator's group at the University of Florida. Both pure-conduction and combined convection-conduction computations have been conducted. A pure-conduction analysis gives insufficient information about the overall thermal distribution. Combined convection-conduction analysis indicates the significant influence of the coolant over the entire flow path. Since 2-D simulation is of limited help, future work on full 3-D modeling of the pump using multi-materials is needed. A comprehensive and accurate model can be developed to take into account the effect of multi-phase flow in the cooling flow loop, and the magnetic interactions.

Senocak, Inanc↗

A massively-parallel, unstructured overset method to simulate moving bodies in turbulent flows

An unstructured overset method capable of performing direct numerical simulation (DNS) and large eddy simulation (LES) of many (O(10 5 )) moving bodies, utilizing many computational cores (O(10 5 )), in turbulent, incompressible fluid flow is presented. Unstructured meshes are attached to bodies and placed within a fixed background domain. Body meshes are allowed to arbitrarily overlap and move throughout the domain. Within each mesh a high resolution, unstructured, non-dissipative finite volume method is used to solve for the flow field. Boundary conditions for each mesh are provided by interpolation from flow solutions on overlapping meshes. When many unstructured meshes of different resolution overlap, care is required in the connection between the different flow solutions. An interpolant is created which seeks to preserve volume conservation of flow quantities between meshes regardless of overlapping mesh differences. An implicit fractional step method is used for time advancement, requiring the calculation of a predicted fluid velocity and corrector pressure field. For the predictor step, the resulting interpolation is directly introduced into the implicit equations for the predicted flow field. For the corrected pressure field, the continuity between meshes is weakly enforced using a penalty formulation. The pressure formulation is symmetric, positive-definite and non-singular resulting in a formulation which is readily solvable using traditional iterative matrix inversion techniques. An Arbitrary Euler-Lagrangian (ALE) method coupled to a 6 degrees of freedom rigid body equation system (6-DOF) is used for body motion. For rotation, a quaternion representation is used to solve Euler's equations of rigid body motion. A linear spring damper model, which uses geometry information readily available from the overset assembly process, is used for collisions. Validation of the method for canonical flow fields is presented including assessment of order of accuracy and kinetic energy conservation properties. Particle-resolved direct numerical simulation (PR-DNS) of single particles in various flow fields are presented for validation. PR-DNS results of 50,000 spherical particles freely moving within turbulent channel flow are shown as a demonstration of the method at full scale. LES results of a marine propeller under crashback conditions are shown to demonstrate the ability to simulate highly unsteady turbulent flows over complex, moving geometries.

97 MATHEMATICS AND COMPUTING↗

Dynamic Earthquake Triggering in Southern California in High Resolution: Intensity, Time Decay, and Regional Variability

Abstract Earthquake triggering by seismic waves has been recognized as a phenomenon for nearly 30 years. However, our ability to study dynamic triggering has been limited by our ability to capture the triggering stresses accurately and record the resultant earthquakes. Here we use full waveforms from a dense seismic network and a modern, high‐resolution seismic catalog to measure triggering in Southern California from 2008 to 2017 based on interevent time ratios. We find that the fractional seismicity rate change, which we term triggering intensity or triggerability, as a function of peak strain change for the period of ∼20 s due to distant earthquakes is monotonically increasing and compatible with earlier measurements made with a disjoint data set from 1984 to 2008. A triggering strain of 1 microstrain is equivalent to the local productivity generated by an M 1.8 earthquakes. This result implies that a prediction of seismicity rate changes can be made based on recorded ground shaking using the same formalism as currently used for aftershock prediction. For a teleseismic event, this small level of triggering occurs throughout the region and thus aggregates to a regional effect. We find that the triggering rate decays after the triggerer follows an Omori‐Utsu law, but at a much slower rate than a typical aftershock sequence. The slow decay rate suggests that an ancillary process such as creep or fluid flow must be part of dynamic triggering. The prevalence of triggering in areas of creep or fluid involvement reinforces this inference. A triggering cascade of secondary earthquakes is insufficient to explain the data.

Miyazawa, Masatoshi↗