Search NASA⌕ Search

SEARCH · Search NASA

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

Exploring Whether Subsurface Fluid Production Can Minimize Triggered Seismicity in Geothermal Fields

Fluid injection and production related to energy recovery and other industrial operations alter the pressure and stress state of subsurface reservoirs which can lead to induced seismicity. Here, the primary goal is to investigate the hypothesis that the modulation of stress state in subsurface reservoirs through fluid injection/production operations can reduce the likelihood of inducing seismicity. Validation of this mitigation strategy will provide an active operational method to control induced seismicity in subsurface reservoir exploitations such as geothermal energy recovery and carbon storage. In this project, we analyze the relationship between extensive fluid production and the paucity of aftershock activity at the Coso Geothermal plant (CGP) following the 2019 Ridgecrest earthquake sequence, where high rates of aftershock triggering were expected. We developed a high-fidelity multiphase coupled thermo-hydro-mechanical (THM) model to simulate fluid injection/production activities at the CGP between 1986 and 2020. THM results of surface subsidence due to high rates of fluid production agree well with field observations from global positioning system (GPS) and interferometric synthetic aperture radar (InSAR). Subsequently, the variations of pore pressure, temperature and stress state were used to drive numerical earthquake simulations of the aftershock response to the 2019 Ridgecrest earthquakes. The earthquake simulations show that seismic quiescence may occur following the Ridgecrest event, depending on the initial stress state at the time of the mainshock and at the initiation of geothermal energy production. Seismic quiescence occurs in 20% of the cases we explored, where rates are decreased by 50% of the background rate in the two years following the mainshock. In circumstances where aftershock rates increase near the CGP following Ridgecrest, the average rate change is a factor of two larger than background rates from simulations with no operations. These findings indicate that unlike many other operations that bring faults closer to failure, operations at the CGP are acting to stabilize faults such that triggering in the current stress state is minimal.

58 GEOSCIENCES↗

Infrared emission spectra of elastohydrodynamic contacts

A small diamond disk mounted as a window in a steel plate was covered with test fluid, and a weighted steel ball was rotated over the window so as to form a sliding elastohydrodynamic (EHD) contact region. Some of the radiant energy generated in this region, both in the fluid and at the boundaries, passed through the window into an infrared interferometer, giving rise to an emission spectrum. This spectrum could be separated into contributions from the fluid and from the ball surface, making it possible, by appropriate calibrations, to estimate their temperatures separately under operating conditions. Moreover, the shape of the discrete spectral bands of the fluid permitted some inferences on its state. Two fluids were studied under identical mechanical conditions, a polyester and a naphthenic oil, each containing an equal amount of polymethylstyrene as a spectral indicator. Differences of band intensity, band width, and frequency could, therefore, be attributed to differences in the behavior of the base fluid.

Lauer, J. L.↗

High-performance parallel analysis of coupled problems for aircraft propulsion

This research program deals with the application of high-performance computing methods to the numerical simulation of complete jet engines. The program was initiated in 1993 by applying two-dimensional parallel aeroelastic codes to the interior gas flow problem of a by-pass jet engine. The fluid mesh generation, domain decomposition and solution capabilities were successfully tested. Attention was then focused on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by these structural displacements. The latter is treated by an ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field fluid elements. New partitioned analysis procedures to treat this coupled 3-component problem were developed in 1994. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers. For the global steady-state axisymmetric analysis of a complete engine we have decided to use the NASA-sponsored ENG10 program, which uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor for parallel versions of ENG10 has been developed. It is planned to use the steady-state global solution provided by ENG10 as input to a localized three-dimensional FSI analysis for engine regions where aeroelastic effects may be important.

Felippa, C. A.↗

High-Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

This research program dealt with the application of high-performance computing methods to the numerical simulation of complete jet engines. The program was initiated in January 1993 by applying two-dimensional parallel aeroelastic codes to the interior gas flow problem of a bypass jet engine. The fluid mesh generation, domain decomposition and solution capabilities were successfully tested. Attention was then focused on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by these structural displacements. The latter is treated by a ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field fluid elements. New partitioned analysis procedures to treat this coupled three-component problem were developed during 1994 and 1995. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers, including the iPSC-860, Paragon XP/S and the IBM SP2. For the global steady-state axisymmetric analysis of a complete engine we have decided to use the NASA-sponsored ENG10 program, which uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor tor parallel versions of ENG10 was developed. During 1995 and 1996 we developed the capability tor the first full 3D aeroelastic simulation of a multirow engine stage. This capability was tested on the IBM SP2 parallel supercomputer at NASA Ames. Benchmark results were presented at the 1196 Computational Aeroscience meeting.

Felippa, C. A.↗

A posteriori error analysis of fluid–structure interactions: Time dependent error

A posteriori error analysis is a technique to quantify the error in particular simulations of a numerical approximation method. In this article, we use such an approach to analyze how various error components propagate in certain moving boundary problems. We study quasi-steady state simulations where slowly moving boundaries remain in mechanical equilibrium with a surrounding fluid. Such problems can be numerically approximated with the Method of Regularized Stokelets (MRS), a popular method used for studying viscous fluid–structure interactions, especially in biological applications. Our approach to monitoring the regularization error of the MRS is novel, along with the derivation of linearized adjoint equations to the governing equations of the MRS with a elastic elements. Our main numerical results provide a clear illustration of how the error evolves over time in several MRS simulations.

97 MATHEMATICS AND COMPUTING↗

Progress in the Physisorption Characterization of Nanoporous Gas Storage Materials

Assessing the adsorption properties of nanoporous materials and determining their structural characterization is critical for progressing the use of such materials for many applications, including gas storage. Gas adsorption can be used for this characterization because it assesses a broad range of pore sizes, from micropore to mesopore. In the past 20 years, key developments have been achieved both in the knowledge of the adsorption and phase behavior of fluids in ordered nanoporous materials and in the creation and advancement of state-of-the-art approaches based on statistical mechanics, such as molecular simulation and density functional theory. Together with high-resolution experimental procedures for the adsorption of subcritical and supercritical fluids, this has led to significant advances in physical adsorption textural characterization. In this short, selective review paper, we discuss a few important and central features of the underlying adsorption mechanisms of fluids in a variety of nanoporous materials with well-defined pore structure. In conclusion, the significance of these features for advancing physical adsorption characterization and gas storage applications is also discussed.

42 ENGINEERING↗

Structural and preliminary thermal performance testing of a pressure activated contact heat exchanger

A contact heat exchanger concept is being developed for use onboard Space Station as an interface device between external thermal bus and pressurized modules. The concept relies on mechanical contact activated by the fluid pressure inside thin-walled tubes. Structural testings were carried out to confirm the technology feasibility of using such thin-walled tubes. The test results also verified the linear elastic stress analysis which was used to predict the tube mechanical behaviors. A preliminary thermal testing was also performed with liquid Freon-11 flowing inside tubes and heat being supplied by electrical heating from the bottom of the contact heat exchanger baseplate. The test results showed excellent agreement of test data with analytical prediction for all thermal resistances except for the two-phase flow characteristics. Testing with two-phase flow inside tubes will, however, be performed on the NASA-JSC test bed.

Lee, C. Y.↗

Technical evaluation report of AGARD Technical Evaluation Meeting on Unsteady Aerodynamics: Fundamentals and Applications to Aircraft Dynamics

From May 6 to 9, 1985, the Fluid Dynamics Panel and Flight Mechanics Panel of AGARD jointly arranged a Symposium on Unsteady Aerodynamics-Fundamentals and Applications to Aircraft Dynamics at the Stadthall, Goettingen, West Germany. This Symposium was organized by an international program committee chaired by Dr. K. J. Orlik-Ruckemann of the Fluid Dynamics Panel. The program consisted of five sessions grouped in two parts: (1) Fundamentals of Unsteady Aerodynamics; and (2) Applications to Aircraft Dynamics. The 35 papers presented at the 4 day meeting are published in AGARD CP 386 and listed in the Appendix. As the papers are already available and cover a very wide field, the evaluators have offered brief comments on every paper, followed by an overall evaluation of the meeting, together with some general conclusions and recommendations.

Mabey, D. G.↗

Modeling of Transient Heat Pipe Operation

The major goal of this project is to develop mathematical models of heat pipes which can be used to predict transient behavior under normal and adverse conditions. The models and solution techniques are to be formulated so that they can be incorporated into existing NASA structural design codes. The major parameters of interest are heat flux distribution, temperature distribution, working fluid pressure distribution, fluid and containment thermal and mechanical properties and geometry. Normal transient operation is taken to be operating conditions where the capillary structure remains fully wetted. Adverse transient operation occurs when drying, re-wetting, choking, non-continuum flow, thawing, freezing, etc., occur in the internal heat pipe working fluid.

Colwell, G. T.↗

Physical property characteristics of gas hydrate-bearing reservoir and associated seal sediments collected during NGHP-02 in the Krishna-Godavari Basin, in the offshore of India

India’s National Gas Hydrate Program Expedition 02 (NGHP-02) was conducted to better understand geologic controls on gas hydrate occurrence and morphology, targeting coarse-grained sediments along the lower continental slope offshore eastern India. This study combines seismic, logging-while-drilling data, and a petroleum system approach to provide a regional geologic and lithologic context for: 1) gas hydrate morphology and distribution, and 2) effects of fine-grained sediments (clays and other grains smaller than 63µm) on gas production in NGHP-02 Area B in the Krishna-Godavari Basin. Area B seismic data show a buried anticline/syncline structure with strong reflectors, R1 and R2, that delineate two of the five lithologic units: Unit I (shallowest), II (R1), III, IV (R2) and V (below the bas of gas hydrate stability). Throughout Area B, gas hydrate morphology depends on its placement within these units. Specifically, core- and grain-scale measurements indicate fines content exerts a primary control on the gas hydrate distribution and morphology. Units I, II and III are generally fine-grained. On the anticline crest, these units host primarily grain-displacing gas hydrate veins linked to pore-occupying gas hydrate in thin, localized, coarser-grained deposits. Diatoms in Unit III increase porosity with depth, reaching ~70% where it contacts Unit IV, the gas hydrate reservoir associated with reflector R2. The Unit III lithology and porosity allow fluid and dissolved-phase methane to escape Unit IV and form gas hydrate in the fine-grained overburden. Within Unit IV, fine-grained layers are interbedded with coarser-grained gas hydrate reservoir sands, and the fines content even in the sands is high enough to impact hydraulic and mechanical properties during gas production. Fluid motion during gas production can mobilize fines, which can then clog pore throats, limiting production rates. In conclusion, pore-water freshening during gas hydrate dissociation can increase fines mobilization, particularly given the smectite identified in the fine-grained interbeds.

03 NATURAL GAS↗

Visualization of flow interaction in a channel flow produced by a pair of helically coiled rods

Helically coiled tubes have been considered in the design of steam generators or heat exchangers of next generation nuclear power plants. Its compactness, efficiency, and the capability to absorb thermal expansion are desirable for nuclear reactor designs. However, the complex flow generated by the helically coiled geometry can result in flow induced vibration (FIV). In order to predict the FIV, there have been efforts using Computational Fluid dynamics (CFD)4, 5 to be coupled with multi-mechanics simulations. However, it is required to validate CFD to have confidence in capturing the flow structure characteristics of the transient complex flow behavior and fluctuation. In the present study, a simplified helically coiled steam generator (HCSG) design experiment was built acquiring data for CFD validation purpose. It is not perfect representation of the multiple rods, however, the current test facility produced transient data of the flow interaction between the rods and flow fluctuations near the rods. Particle Image Velocimetry (PIV) was applied to visualize the turbulent flow in a channel with slant rods. The two dimensional turbulent flow fields were analyzed. The velocity components were defined as below. u = U + u' (1) v = V + v' (2) where u is the lateral velocity on x-axis, v is the vertical velocity on y-axis, U and V are the mean velocities, and u' and v' are the fluctuation. Reynolds number, Re, was calculated using the rod diameter, d = 15.9 mm, and the upstream approach velocity, V{sub 0} = 0.5 m/s. Re = V{sub 0} d/η (3) where η is the kinematic viscosity of water. The experiment was conducted at Re ≅ 9,000 considering the currently reasonable computing power. Flow visualization in a flow channel produced by helically coiled rods was performed using PIV. The mean velocity field and streamlines were presented in this paper. The vortical flow structures in the cavity between rods were observed. Because of the asymmetric geometry, the flow profile was skewed. Additional measurements and transient analysis will be conducted to produce benchmark data for FIV study using CFD. Also, a test facility with a 5-rod heat exchanger design is operating to produce more data as a further work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Ionization effects in three-dimensional solar granulation simulations

These numerical studies show that ionization influences both the transport and dynamical properties of compressible convection near the surface of the Sun. About two-thirds of the enthalpy transported by convective motions in the region of partial hydrogen ionization is carried as latent heat. The role of fast downflow plumes in total convective transport is substantially elevated by this contribution. Instability of the thermal boundary layer is strongly enhanced by temperature sensitive variations in the radiative properties of the fluid, and this provides a mechanism for plume initiation and cell fragmentation in the surface layers. As the plumes descend, temperature fluctuations and associated buoyancy forces are maintained because of the increased specific heat of the partially ionized material. This can result is supersonic vertical flows. At greater depths, ionization effects diminish, and the plumes are decelerated by significant entrainment of surrounding fluid.

Rast, Mark P.↗

Surveillance of Ocular Parameters and Visual Function in Bed Rest Subjects

Recent visual changes in astronauts have raised concern about ocular health during long duration spaceflight. Seven cases have been documented in astronauts who spent 6 months aboard the International Space Station. These astronauts were male ranging in age from 45 to 55 years old. All astronauts exhibited pre- to post flight refractive changes. Decreased intraocular pressure (IOP) post flight was observed in 3 cases. Fundoscopic exams revealed post flight findings of choroidal folds in 4 cases, optic disc edema in 5 cases and the presence of cotton wool spots in 3 cases. Optical coherence tomography (OCT) confirmed findings of choroidal folds and disc edema, and also documented retinal nerve fiber layer thickening (5 cases). Findings from MRI examinations showed posterior globe flattening (5 cases), optic nerve sheath distention (6 cases) and torturous optic nerves (2 cases). Of the 7 cases, intracranial pressure was measured on 4 astronauts. These 4 showed elevated ICP post-flight that remained elevated for as long as 19 months in one case. While the etiology remains unknown, hypotheses speculate that venous insufficiency or hypertension in the brain caused by cephalad fluid shifts during spaceflight are possible mechanisms for ocular changes seen in astronauts. Head-down tilt bed rest is a spaceflight analog that induces cephalad fluid shifts. This study is designed to provide ocular monitoring of bed rest subjects and determine whether clinically relevant changes are found. Ocular Changes

Cromwell, Ronita L.↗

Liquid-induced topological transformations of cellular microstructures

The fundamental topology of cellular structures—the location, number and connectivity of nodes and compartments—can profoundly affect their acoustic, electrical5, chemical, mechanical and optical properties, as well as heat, fluid and particle transport. Approaches that harness swelling, electromagnetic actuation and mechanical instabilities in cellular materials have enabled a variety of interesting wall deformations and compartment shape alterations, but the resulting structures generally preserve the defining connectivity features of the initial topology. Achieving topological transformation presents a distinct challenge for existing strategies: it requires complex reorganization, repacking, and coordinated bending, stretching and folding, particularly around each node, where elastic resistance is highest owing to connectivity. Here we introduce a two-tiered dynamic strategy that achieves systematic reversible transformations of the fundamental topology of cellular microstructures, which can be applied to a wide range of materials and geometries. Our approach requires only exposing the structure to a selected liquid that is able to first infiltrate and plasticize the material at the molecular scale, and then, upon evaporation, form a network of localized capillary forces at the architectural scale that ‘zip’ the edges of the softened lattice into a new topological structure, which subsequently restiffens and remains kinetically trapped. Reversibility is induced by applying a mixture of liquids that act separately at the molecular and architectural scales (thus offering modular temporal control over the softening–evaporation–stiffening sequence) to restore the original topology or provide access to intermediate modes. Guided by a generalized theoretical model that connects cellular geometries, material stiffness and capillary forces, we demonstrate programmed reversible topological transformations of various lattice geometries and responsive materials that undergo fast global or localized deformations. In conclusion, we then harness dynamic topologies to develop active surfaces with information encryption, selective particle trapping and bubble release, as well as tunable mechanical, chemical and acoustic properties.

42 ENGINEERING↗

CFD simulations of electric motor end ring cooling for improved thermal management

Proper thermal management of an electric motor for vehicle applications extends its operating range. One cooling approach is to impinge Automatic Transmission Fluid (ATF) onto the rotor end ring. Increased ATF coverage correlates to enhanced heat transfer. Computational Fluid Dynamics (CFD) analytical tools provide a mechanism to assess motor thermal management prior to hardware fabrication. The complexity of the fluid flow (e.g., jet atomization, interface tracking, wall impingement) and heat transfer makes these simulations challenging. Computational costs are high when solving these flows on high-speed rotating meshes. Typically, a Volume-of Fluid (VOF) technique (i.e., two-fluid system) is used to resolve ATF dynamics within this rotating framework. Suitable numerical resolution of the relevant physics for thin films under strong inertial forces at high rotor speeds is computationally expensive, further increasing the run times. In this work, a numerical study of rotor-ring cooling by ATF is presented using a patent automated Cartesian cut-cell based method coupled with Automatic Mesh Refinement (AMR). This approach automatically creates the Cartesian mesh on-the-fly and can effectively handle complex rotating geometries by adaptively refining the mesh based on local gradients in the flow field which results in better resolution of the air-ATF interface. A Single non-inertial Reference Frame (SRF) approach is used to account for the rotating geometry and to further improve the overall computational efficiency. Quasi-steady state conditions are targeted in the analysis of the results. Important physics such as ATF jet structure, velocity detail near the air-jet interface, ATF coverage/accumulation on the ring surface, and cooling capacity are presented for a low-resolution Reynolds averaged Navier-Stokes (RANS), high-resolution RANS, and high-resolution Large-Eddy Simulation (LES) models. Computations are scaled over hundreds of cores on a supercomputer to maximize turnaround time. Each numerical approach is shown to capture the general trajectory of the oil jet prior to surface impingement. The high-resolution LES simulation, however, is superior in capturing small scale details and heat transfer between the free jet and surrounding air.

42 ENGINEERING↗

Micro- and nano-porosity of the active Alpine Fault zone, New Zealand

Porosity reduction in rocks from a fault core can cause elevated pore fluid pressures and consequently influence the recurrence time of earthquakes. We investigated the porosity distribution in the New Zealand's Alpine Fault core in samples recovered during the first phase of the Deep Fault Drilling Project (DFDP-1B) by using two-dimensional nanoscale and three-dimensional microscale imaging. Synchrotron X-ray microtomography-derived analyses of open pore spaces show total microscale porosities in the range of 0.1 %–0.24 %. These pores have mainly non-spherical, elongated, flat shapes and show subtle bipolar orientation. Scanning and transmission electron microscopy reveal the samples' microstructural organization, where nanoscale pores ornament grain boundaries of the gouge material, especially clay minerals. Our data imply that (i) the porosity of the fault core is very small and not connected; (ii) the distribution of clay minerals controls the shape and orientation of the associated pores; (iii) porosity was reduced due to pressure solution processes; and (iv) mineral precipitation in fluid-filled pores can affect the mechanical behavior of the Alpine Fault by decreasing the already critically low total porosity of the fault core, causing elevated pore fluid pressures and/or introducing weak mineral phases, and thus lowering the overall fault frictional strength. We conclude that the current state of very low porosity in the Alpine Fault core is likely to play a key role in the initiation of the next fault rupture.

58 GEOSCIENCES↗

Simulating Fracture Sealing by Granular LCM Particles in Geothermal Drilling

Lost circulation occurs when the returned fluid is less than what is pumped into the well due to loss of fluid to pores or fractures. A lost-circulation event is a common occurrence in a geothermal well. Typical geothermal reservoirs are often under-pressured and have larger fracture apertures. A severe lost-circulation event is costly and may lead to stuck pipe, well instability, and well abandonment. One typical treatment is adding lost-circulation materials (LCMs) to seal fractures. Conventional LCMs fail to properly seal fractures because their mechanical limit is exceeded at elevated temperatures. In this paper, parametric studies in numerical simulations are conducted to better understand different thermal effects on the sealing mechanisms of LCMs. The computational fluid dynamics (CFDs) and the discrete element method (DEM) are coupled to accurately capture the true physics of sealing by granular materials. Due to computational limits, the traditional Eulerian–Eulerian approach treats solid particles as a group of continuum matter. With the advance of modern computational power, particle bridging is achievable with DEM to track individual particles by modeling their interactive forces between each other. Particle–fluid interactions can be modeled by coupling CFD algorithms. Fracture sealing capability is investigated by studying the effect of four individual properties including fluid viscosity, particle size, friction coefficient, and Young’s modulus. It is found that thermally degraded properties lead to inefficient fracture sealing.

15 GEOTHERMAL ENERGY↗