Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid Structure Interaction”

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 865 records · Page 48

Water Imbibition and Oil Recovery in Shale: Dynamics and Mechanisms Using Integrated Centimeter-to-Nanometer-Scale Imaging

Water imbibition, and the associated oil displacement, is an important process in shale oil reservoirs after hydraulic fracturing and in water-based enhanced oil recovery (EOR). Current techniques for water imbibition measurement are mostly “black-box”-type methods. A more explicit understanding of the water imbibition/oil recovery dynamics and geological controls is in demand. In this paper, a multiscale imaging technique that covers centimeter to nanometer scale (i.e., core to pore scale), integrating neutron radiography, microcomputed tomography (micro-CT), and scanning electron microscope (SEM) is applied to investigate the water imbibition depth and rate and the cause of heterogeneity of imbibition in shale samples. The dynamic processes of water imbibition in the 1-in. (25.4-mm) core sample were explicitly demonstrated, and the imbibition along the matrix and imbibition through microfractures are distinguished through neutron radiography image analysis. The causes of observed imbibition heterogeneity were further investigated through micro-CT and SEM image analysis for 1.5-mm diameter miniplug samples from different laminas of the 1-in. core samples. Imbibition depth and rate were calculated on the basis of image analysis as well. Estimation of oil recovery through water imbibition in shale matrix was performed for an example shale field. Here this innovative and integrated multiscale imaging technique provides a “white/gray-box” method to understand water imbibition and water-oil displacement in shale. The wide span of the length scale (from centimeter to nanometer) of this technique enables a more comprehensive, accurate, and specific understanding of both the core-scale dynamics and pore-scale mechanisms of water imbibition, oil recovery, and matrix-fracture interaction.

04 OIL SHALES AND TAR SANDS↗

Subsurface Characterization of Hydraulic Fracture Test Site-2 (HFTS-2), Delaware Basin

Hydraulic Fracturing Test Site-2 (HFTS-2) is a field-based research experiment performed in the Wolfcamp Formation of the Permian (Delaware) Basin. This paper focuses on integration, advanced geological characterization, and 3D subsurface modeling of the comprehensive HFTS-2 dataset. The study showcases a multidisciplinary reservoir characterization approach that incorporates geology, petrophysics, geochemistry, geomechanics, microseismic, and subsurface engineering analysis. Subsurface characterization of organic-rich mudstone formations requires understanding complex hydraulic fracture network growth in relation to inherent lithology, geomechanical properties, and interaction with pre-existing natural fractures. This paper presents a characterization workflow incorporating pre- and post-stimulation subsurface data, unique to the HFTS-2 dataset. The study integrated: (1) rock properties from logs, cores, and thin sections; (2) natural and hydraulic fracture descriptions from cores and image logs; (3) local and regional stresses; (4) geomechanics; (5) microseismic; (6) fiber optic (FO) and bottomhole pressure gauge (BHPG) response; and (7) produced fluids analysis. During a stimulation treatment, creation of the stimulated rock volume (SRV) is influenced by several subsurface factors. Key contributing factors include structural context, stress conditions, lithology, facies architecture, pre-existing natural fractures, and geomechanical properties. The HFTS-2 subsurface data integration indicates that the SRV is comprised of a complex juxtaposition of hydraulic fracture swarms, as evidenced by image logs analysis, core description, and microseismic monitoring. The HFTS-2 microseismic event density was used to generate 3D heat maps that serve as a representative SRV footprint, corroborated by secondary datasets. These maps were further integrated with petrophysical and geomechanical characteristics, as well as responses from FO and BHPG, to estimate the lateral and vertical dimensions of the effective fractures. The geological characterization for the HFTS-2 dataset combined with 3D modeling for petrophysical and geomechanical properties provides a strong foundation for subsurface simulation and optimization studies. Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/1-21URTC/D011S005R001/2477501/urtec-2021-5243-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5243 The workflow improved our understanding of HFTS-2 hydraulic fracture propagation and characteristics in relation to offset pressure depletion and interaction with pre-existing natural fractures. Analysis showed that fracture geometry varies by stage and by well, and a complex fracture network is generated with varying fracture density. The multidisciplinary workflow presented herein for integration and characterization serves as a foundation to evaluate completion efficiency and estimate areal and vertical stimulation and depletion extent for the project. Furthermore, the workflow and learnings can also be transferred to other unconventional plays.

58 GEOSCIENCES↗

Carbonate formation in Wyoming montmorillonite under high pressure carbon dioxide

Carbonation reaction with silicate minerals that are common components of the host rock and cap rock within geological storage reservoirs and the associated structural deformation were investigated for better understanding of the geochemical reactions associated with geologic CO 2 storage. Exposure of a model expanding clay, Wyoming montmorillonite, SWy-2, to high-pressure CO 2 resulted in the formation of a mineral carbonate phase via dry CO 2 –clay mineral interactions at two different temperatures. The experimental evidence suggests that the properties of CO 2 fluid at 70 °C provide more favorable conditions for carbonate formation at the clay surface less accessible to CO 2 at 22 °C. The carbonation reaction occurred predominantly within the first couple of days of exposure to the fluid and then proceeded slower with continuing exposure. As compared to the as-received clay under the same ambient conditions, the (0 0 1) basal spacing of the clay bearing carbonates (after the CO 2 exposure) was slightly expanded at a relative humidity (RH) level of 12% but it was slightly collapsed at the RH level of 40%. Finally, experimental observations suggest that the carbonation reaction occurs at the external surface as well as internal surface (interlayer) of the clay particles.

Hur, Tae-Bong↗

Collaborative Research: Explosive reconnection in relativistic magnetically-dominated plasmas (Final Report)

Magnetic fields can be the dominant component of astrophysical plasmas, so that the magnetic energy density might exceed even the rest-mass energy density of matter. In this extreme (and largely unexplored) regime the magnetic field controls the overall plasma evolution, dissipation, and acceleration of non-thermal particles. This plasma regime, applicable to a variety of astrophysical sources - magnetars, pulsars and pulsar wind nebulae (PWNe), jets of Active Galactic Nuclei (AGNs) and Gamma-Ray Bursters (GRBs) - is dramatically different from laboratory plasmas, the magnetospheres of planets, and the interplanetary plasma. Relativistic astrophysical sources then provide an unique opportunity to study the fundamental plasma physics of magnetically-dominated plasmas; a novel and fast-evolving field of theoretical research which, by investigating energy conversion and particle energization processes in plasmas, is of vital importance to the Fusion Energy Sciences DoE program. Data coming from astrophysical high-energy missions, especially the Crab Nebula flares recently observed by the Fermi and AGILE satellites, suggest that the acceleration of non-thermal particles to the highest energies occurs in magnetic reconnection events - a major change of paradigm in high-energy astrophysics. Most importantly, observations demand that particle acceleration should proceed extremely fast (with accelerating electric field of the order of the magnetic field) and on macroscopic scales (much larger, e.g., than the microscopic plasma skin depth). We are conducting studies of the microphysics of magnetically-dominated plasmas focussing in particular on the highly dynamic regime of explosive reconnection and associated particle acceleration in relativistic plasmas. We are studying the stability and explosive plasma dynamics, particle acceleration and radiation production in a number of idealized plasma configurations that approximate relevant astrophysical sources (like the magnetic ABC structures and interacting flux tubes, as well as generalizations of analytical models of X-point collapse to relativistic plasmas). The well-studied case of the Crab Nebula is taken as a prototypical example for the application of the model. We are combining analytical studies of explosive magnetic dynamics and dissipation in relativistic plasmas with particle-in-cell (PIC) simulations and fluid simulations. The theoretical model, fluid and particle-in-cell simulations are cross-checked for agreement and convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SPAR improved structure/fluid dynamic analysis capability

The capability of analyzing a coupled dynamic system of flowing fluid and elastic structure was added to the SPAR computer code. A method, developed and adopted for use in SPAR utilizes the existing assumed stress hybrid plan element in SPAR. An operational mode was incorporated in SPAR which provides the capability for analyzing the flaw of a two dimensional, incompressible, viscous fluid within rigid boundaries. Equations were developed to provide for the eventual analysis of the interaction of such fluids with an elastic solid.

Oden, J. T.↗

Calculations of inlet distortion induced compressor flow field instability

Calculations of the onset of flow instability are carried out for low-speed multistage axial compressors operating with asymmetric inlet flow. The modeling of the fluid dynamic interaction between the spoiled and unspoiled sectors of the compressor is the most important feature of the calculation procedure. The calculations show that annulus average slope of the compressor pressure rise characteristic equal to zero is a useful approximate stability criterion for situations where the dynamics of the compressor flow field do not couple strongly to the compression system or the structure of the imposed distortion is not similar to that of the eigenmodes of the flow in the compressor annulus. This criterion is employed to investigate the relationship between the present model and the 'parallel compressor' model. Calculations are also presented for cases when compressor and compressor system are closely coupled, and situations in which the compressor is subjected to a rotating distortion. These first-of-a-kind computations, and the accompanying description of the physical mechanisms, show that the stability of the flow in the compressor can be adversely affected if the temporal or spatial structure of the distortion is such that resonant type responses can be evoked either from the compressor or from compressor/compression system interactions.

Chue, R.↗

Group-kinetic theory and modeling of atmospheric turbulence

A group kinetic method is developed for analyzing eddy transport properties and relaxation to equilibrium. The purpose is to derive the spectral structure of turbulence in incompressible and compressible media. Of particular interest are: direct and inverse cascade, boundary layer turbulence, Rossby wave turbulence, two phase turbulence; compressible turbulence, and soliton turbulence. Soliton turbulence can be found in large scale turbulence, turbulence connected with surface gravity waves and nonlinear propagation of acoustical and optical waves. By letting the pressure gradient represent the elementary interaction among fluid elements and by raising the Navier-Stokes equation to higher dimensionality, the master equation was obtained for the description of the microdynamical state of turbulence.

Tchen, C. M.↗

Turbine blade tip/outer airseal/platform flow analysis

The interaction between the main passage flow through a high pressure turbine and the secondary flows due to leakage through the wheelspace rim seals is reviewed. Various experimental and analytical studies of these interactions are outlined and a numerical investigation of the interaction between the main passage flow and a wheelspace cavity seal flow is described. The numerical investigation uses a structured grid method to study the overall interaction between the turbine stage components and the wheelspace seal flow, and an unstructured grid method to resolve the detailed flow features within the geometrically complex cavity seal. The numerical results agree with various observations from experimental studies under similar flow conditions. As the flow rate through the rim cavity seal is increased, the ingestion of fluid from the main passage flow into the rim seal area decreases rapidly. A small amount of main flow gas is ingested into the rim seal area, even at high seal flow rates.

Hah, C.↗

Hybrid Vlasov-Maxwell simulations of two-dimensional turbulence in plasmas

Turbulence in plasmas is a very challenging problem since it involves wave-particle interactions, which are responsible for phenomena such as plasma dissipation, acceleration mechanisms, heating, temperature anisotropy, and so on. In this work, a hybrid Vlasov-Maxwell numerical code is employed to study local kinetic processes in a two-dimensional turbulent regime. In the present model, ions are treated as a kinetic species, while electrons are considered as a fluid. As recently reported in [S. Servidio, Phys. Rev. Lett. 108, 045001 (2012)], nearby regions of strong magnetic activity, kinetic effects manifest through a deformation of the ion velocity distribution function that consequently departs from the equilibrium Maxwellian configuration. Here, the structure of turbulence is investigated in detail in phase space, by evaluating the high-order moments of the particle velocity distribution, i.e., temperature, skewness, and kurtosis. This analysis provides quantitative information about the non-Maxwellian character of the system dynamics. This departure from local thermodynamic equilibrium triggers several processes commonly observed in many astrophysical and laboratory plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stern Layer Structure and Energetics at Mica-Water Interfaces

The screening of surface charge by dissolved ions at solid liquid interfaces in the region of interfacial fluid known as the electrical double layer (EDL)-plays a recurrent role in surface science, from ion adsorption to colloidal mechanics to the transport properties of nanoporous media. A persistent unknown in theories of EDL-related phenomena is the structure of the Stern layer, the near-surface portion of the EDL where water molecules and adsorbed ions form specific short-range interactions with surface atoms. Here, we describe a set of synchrotron X-ray reflectivity (XRR) experiments and molecular dynamics (MD) simulations carried out under identical conditions for a range of 0.1 M alkali chloride (Li-, Na-, K-, Rb-, or CsCl) solutions on the basal surface of muscovite mica, a mineral isostructural to phyllosilicate clay minerals and one of the most widely studied reference surfaces in interfacial science. Our XRR and MD simulation results provide a remarkably consistent view of the structure and energetics of the Stern layer, with some discrepancy on the fraction of the minor outer-sphere component of Rb and on the adsorption energetics of Li. The results of both techniques, along with surface complexation model calculations, provide insight into the sensitivity of water structure and ion adsorption to surface topography and the type of adsorbed counterion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Design and Preparation of Nematic Colloidal Particles

Colloidal systems are abundant in technology, in biomedical settings, and in our daily life. The so-called “colloidal atoms” paradigm exploits interparticle interactions to self-assemble colloidal analogs of atomic and molecular crystals, liquid crystal glasses, and other types of condensed matter from nanometer- or micrometer-sized colloidal building blocks. Nematic colloids, which comprise colloidal particles dispersed within an anisotropic nematic fluid host medium, provide a particularly rich variety of physical behaviors at the mesoscale, not only matching but even exceeding the diversity of structural and phase behavior in conventional atomic and molecular systems. This feature article, using primarily examples of works from our own group, highlights recent developments in the design, fabrication, and self-assembly of nematic colloidal particles, including the capabilities of preprogramming their behavior by controlling the particle’s surface boundary conditions for liquid crystal molecules at the colloidal surfaces as well as by defining the shape and topology of the colloidal particles. Here recent progress in defining particle-induced defects, elastic multipoles, self-assembly, and dynamics is discussed along with open issues and challenges within this research field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigation of the interaction between competing types of nondiffusive transport in drift wave turbulence

Radial transport in turbulence dominated tokamak plasmas has been observed to deviate from classical diffusion in certain regimes relevant for magnetic confinement fusion. These situations at least include near-marginal turbulence, where radial transport becomes superdiffusive and mediated by elongated radial structures (or avalanches) and transport across radially sheared poloidal flows, where radial subdiffusion often ensues. In this paper, the interaction between very different physical ingredients responsible for these two types of nondiffusive dynamics (namely, turbulent profile relaxation close to a local threshold and the interaction with radially sheared zonal flows) is studied in detail in the context of a simple two-dimensional electrostatic plasma fluid turbulence model based on the dissipative trapped electron mode. It is shown that, depending on the relative relevance of each of these ingredients, which can be tuned in various ways, a variety of non-diffusive radial transport behaviors can be found in the system. The results also illustrate the fact that the classical diffusion paradigm is often insufficient to describe turbulent transport in systems with self-generated flows and turbulent profile relaxations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mechanotransduction and the functional response of bone to mechanical strain

Mechanotransduction plays a crucial role in the physiology of many tissues including bone. Mechanical loading can inhibit bone resorption and increase bone formation in vivo. In bone, the process of mechanotransduction can be divided into four distinct steps: (1) mechanocoupling, (2) biochemical coupling, (3) transmission of signal, and (4) effector cell response. In mechanocoupling, mechanical loads in vivo cause deformations in bone that stretch bone cells within and lining the bone matrix and create fluid movement within the canaliculae of bone. Dynamic loading, which is associated with extracellular fluid flow and the creation of streaming potentials within bone, is most effective for stimulating new bone formation in vivo. Bone cells in vitro are stimulated to produce second messengers when exposed to fluid flow or mechanical stretch. In biochemical coupling, the possible mechanisms for the coupling of cell-level mechanical signals into intracellular biochemical signals include force transduction through the integrin-cytoskeleton-nuclear matrix structure, stretch-activated cation channels within the cell membrane, G protein-dependent pathways, and linkage between the cytoskeleton and the phospholipase C or phospholipase A pathways. The tight interaction of each of these pathways would suggest that the entire cell is a mechanosensor and there are many different pathways available for the transduction of a mechanical signal. In the transmission of signal, osteoblasts, osteocytes, and bone lining cells may act as sensors of mechanical signals and may communicate the signal through cell processes connected by gap junctions. These cells also produce paracrine factors that may signal osteoprogenitors to differentiate into osteoblasts and attach to the bone surface. Insulin-like growth factors and prostaglandins are possible candidates for intermediaries in signal transduction. In the effector cell response, the effects of mechanical loading are dependent upon the magnitude, duration, and rate of the applied load. Longer duration, lower amplitude loading has the same effect on bone formation as loads with short duration and high amplitude. Loading must be cyclic to stimulate new bone formation. Aging greatly reduces the osteogenic effects of mechanical loading in vivo. Also, some hormones may interact with local mechanical signals to change the sensitivity of the sensor or effector cells to mechanical load.

Non-NASA Center↗

Bayesian estimations of orientation distribution functions from small-angle scattering enable direct prediction of mechanical stress in anisotropic materials

Properties of soft materials are influenced by their anisotropic structuring under nonequilibrium fields. Although anisotropic structure-property relationships have been extensively explored theoretically, comparison to experiments requires determination of the microstructural orientation probability distribution function (OPDF) of microstructural elements. Small angle scattering (SAS) measurements encode information about the OPDF, but tools to navigate this connection are incomplete. Here, we develop and validate an explicit framework to link arbitrary OPDFs to SAS measurements. Specifically, we propose, validate, and apply a method, maximum a posteriori scattering inference (MAPSI), whereby the OPDF may be obtained from SAS measurements using a Bayesian estimation method. Using this method, we obtain estimates of the full 3D OPDF for two model semidilute fd-virus (rodlike) dispersions at concentrations that are approximately equal to and twice the overlap concentration. From the OPDF, we calculate its second and fourth moments and compare these to predictions for a dilute suspension of rigid rods and to a recent theory for semidilute suspensions. Finally, we use both the theoretical and measured moments to calculate the stress, both for dilute and semidilute suspensions. These predictions are not only compared to each other, but also to measured values of the shear stress, and point to new insights into the behavior of suspensions of highly elongated particles in the transition between dilute and semidilute behavior. We also use this new framework to provide perspective on the connection between scalar parameterizations of scattering and the OPDF that have frequently been used in the past. Here, the new tools developed in this work provide an unprecedented path toward experimental validation of dynamical theories of rodlike colloids and polymers, and for measurement of nonequilibrium structures and stresses of other complex fluids and soft materials with SAS.

36 MATERIALS SCIENCE↗

Overview of the Helios Version 2.0 Computational Platform for Rotorcraft Simulations

This article summarizes the capabilities and development of the Helios version 2.0, or Shasta, software for rotary wing simulations. Specific capabilities enabled by Shasta include off-body adaptive mesh refinement and the ability to handle multiple interacting rotorcraft components such as the fuselage, rotors, flaps and stores. In addition, a new run-mode to handle maneuvering flight has been added. Fundamental changes of the Helios interfaces have been introduced to streamline the integration of these capabilities. Various modifications have also been carried out in the underlying modules for near-body solution, off-body solution, domain connectivity, rotor fluid structure interface and comprehensive analysis to accommodate these interfaces and to enhance operational robustness and efficiency. Results are presented to demonstrate the mesh adaptation features of the software for the NACA0015 wing, TRAM rotor in hover and the UH-60A in forward flight.

Sankaran, Venkateswaran↗

Cold Flow Propulsion Test Complex Pulse Testing

When the propellants in a liquid rocket engine burn, the rocket not only launches and moves in space, it causes forces that interact with the vehicle itself. When these interactions occur under specific conditions, the vehicle's structures and components can become unstable. One instability of primary concern is termed pogo (named after the movement of a pogo stick), in which the oscillations (cycling movements) cause large loads, or pressure, against the vehicle, tanks, feedlines, and engine. Marshall Space Flight Center (MSFC) has developed a unique test technology to understand and quantify the complex fluid movements and forces in a liquid rocket engine that contribute strongly to both engine and integrated vehicle performance and stability. This new test technology was established in the MSFC Cold Flow Propulsion Test Complex to allow injection and measurement of scaled propellant flows and measurement of the resulting forces at multiple locations throughout the engine.

McDougal, Kris↗

Formation and Topology of Foreshock Bubbles

We use global and local hybrid (kinetic ions and fluid electrons) simulations to investigate the conditions under which foreshock bubbles (FBs) form and how their topology changes with solar wind conditions. FBs form as a result of the interaction between solar wind discontinuities and backstreaming ion beams in the foreshock. They consist of an outer shock and its associated sheath plasma and a low density high temperature core with low magnetic field strength. The structure of FBs is determined by the angle between the interplanetary magnetic field and the normal to the solar wind discontinuity. We show that interaction of rotational discontinuities with the foreshock during small angles between the interplanetary magnetic field and discontinuity normal results in the formation of a nearly spherical bubble with a radius that scales with the width of the foreshock. As this angle increases, FBs become more elongated and eventually become nearly planar structures with dimensions that scale with the length of the foreshock. Despite this transformation, the signatures of FBs in spacecraft time series data remain the same in agreement with the observations. Global simulation results show that FBs form when the solar wind flow speed corresponds to high or intermediate Alfvén Mach numbers (approximately >7 MA). In general, this is tied to the relative speed between the solar wind and ion beams and drop in density of the back streaming ions.

Nick Omidi↗

Computational Evaluation of an OML-based Heat Exchanger Concept for HEATheR

Feasibility of outer mold line (OML) heat exchangers for electrified aircraft is computationally explored for three different aircraft concepts within The High-efficiency Electrified Aircraft Thermal Research (HEATheR) project. OML surface temperature limits were selected based on aluminum-alloy and carbon-fiber composite materials commonly used in transport aircraft. Heat flux distribution over the aircraft OML at the surface temperature limit was evaluated with computational fluid dynamics (CFD) analysis. The candidate OML locations and the size for the implementation of the OML-based heat exchanger is down-selected using this initial analysis and other considerations including the concern for structural integrity and proximity to heat-generating electrical components. The sensitivity of the heat transfer analysis to the computational grid, aircraft angle of attack, surface temperature, and the interaction of individual OML heat exchanger patches were investigated.

CAS↗