A conservative level set method for N-phase flows with a free-energy-based surface tension model
Not Available
Howard, Amanda A. (ORCID:0000000264116198)↗
SEARCH · Search NASA
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.
A method to generate small-scale, high-resolution sedimentary bedform architecture models representing realistic geologic facies
Small-scale (mm to m) sedimentary structures (e.g. ripple lamination, cross-bedding) have received a great deal of attention in sedimentary geology. The influence of depositional heterogeneity on subsurface fluid flow is now widely recognized, but incorporating these features in physically-rational bedform models at various scales remains problematic. The current investigation expands the capability of an existing set of open-source codes, allowing generation of high-resolution 3D bedform architecture models. The implemented modifications enable the generation of 3D digital models consisting of laminae and matrix (binary field) with characteristic depositional architecture. The binary model is then populated with petrophysical properties using a textural approach for additional analysis such as statistical characterization, property upscaling, and single and multiphase fluid flow simulation. One example binary model with corresponding threshold capillary pressure field and the scripts used to generate them are provided, but the approach can be used to generate dozens of previously documented common facies models and a variety of property assignments. An application using the example model is presented simulating buoyant fluid (CO 2 ) migration and resulting saturation distribution.
58 GEOSCIENCES↗
Simplified Pore Network Model for Analysis of Flow Capacity and Residual Oil Saturation Distribution Based on Computerized Tomography
The pore network (connected tunnels) can be used not only to describe pore geometry inside a rock but also to calculate various fluid flow properties on a micro scale. Based on a “digital core” extracted from micro computerized tomography (μ-CT), we have developed a simplified pore network model. This study simplifies the pore structure as a lattice of cubes. Pore size, pore distribution, tunnel connectivity, fluid viscosity, interfacial tension, and the external driving pressure gradient are considered for each adjacent cube to evaluate the flow capacity of the pore network. We calculated the permeability, effective permeability, and effective sweep efficiency. The results show that fluid flow initially starts in the largest tunnel, which corresponds to the threshold pressure of the rock sample (fluid flow onset pressure). Permeability changes nonlinearly before all the tunnels are involved in the flow. Changes in the microscopic residual oil saturation and the sweep efficiency were visualized based on the three-dimensional pore distribution. Visualization shows that more oil remains in the formation, the heterogeneity is greater, and the oil viscosity is higher.
58 GEOSCIENCES↗
2017 GTO Project review Laboratory Evaluation of EGS Shear Stimulation.
The objectives and purpose of this research has been to produce laboratory-based experimental and numerical analyses to provide a physics-based understanding of shear stimulation phenomena (hydroshearing) and its evolution during stimulation. Water was flowed along fractures in hot and stressed fractured rock, to promote slip. The controlled laboratory experiments provide a high resolution/high quality data resource for evaluation of analysis methods developed by DOE to assess EGS “behavior” during this stimulation process. Segments of the experimental program will provide data sets for model input parameters, i.e., material properties, and other segments of the experimental program will represent small scale physical models of an EGS system, which may be modeled. The coupled lab/analysis project has been a study of the response of a fracture in hot, water-saturated fractured rock to shear stress experiencing fluid flow. Under this condition, the fracture experiences a combination of potential pore pressure changes and fracture surface cooling, resulting in slip along the fracture. The laboratory work provides a means to assess the role of “hydroshearing” on permeability enhancement in reservoir stimulation. Using the laboratory experiments and results to define boundary and input/output conditions of pore pressure, thermal stress, fracture shear deformation and fluid flow, and models were developed and simulations completed by the University of Oklahoma team. The analysis methods are ones used on field scale problems. The sophisticated numerical models developed contain parameters present in the field. The analysis results provide insight into the role of fracture slip on permeability enhancement-“hydroshear” is to be obtained. The work will provide valuable input data to evaluate stimulation models, thus helping design effective EGS.
58 GEOSCIENCES↗
Convection in grain refining
The relationship between fluid flow phenomena, nucleation, and grain refinement in solidifying metals both in the presence and in the absence of a gravitational field was investigated. The reduction of grain size in hard-to-process melts; the effects of undercooling on structure in solidification processes, including rapid solidification processing; and control of this undercooling to improve structures of solidified melts are considered. Grain refining and supercooling thermal modeling of the solidification process, and heat and fluid flow phenomena in the levitated metal droplets are described.
Flemings, M. C.↗
Basin-wide gradients in dolomite d26Mg and 87Sr/86Sr identify the structural centre of the Williston Basin as a source of magnesium for dolomitization of the Middle Devonian Winnipegosis Formation
Dolomitized carbonates in the Williston Basin are common below the Middle Devonian Prairie Evaporite Formation, indicating a plausible role for brine-reflux as the source of Mg for dolomitization. However, the Lower Member Winnipegosis Formation, which directly underlies the evaporite sequence, exhibits dolomite 87Sr/86Sr ratios that are higher than expected for Middle Devonian seawater, instead suggesting that the dolomitizing fluids ascended from the deep basin where present-day formation waters are known to be more radiogenic. This interpretation is supported in this study by Mg isotope data showing a basin-scale gradient of increasing d26Mg values from the deep center (–2.03‰) to the edge (–1.34‰) of the Williston Basin along the lower Member Winnipegosis bed. Dolomite favors light isotopes of Mg when it forms, thus driving dolomitizing solutions, and the produced dolomite, towards heavier isotopic compositions in the direction of fluid flow. However, dolomite is also Mg-deficient when it forms, and its initial Mg isotope composition is susceptible to alteration by fluid-rock interactions occurring later in the depositional history of the deposit. We argue that the Mg isotope gradient in the Lower Member Winnipegosis was created by an ascending flow of isotopically light Mg-bearing fluid that differentially altered an earlier formed brine-reflux dolomite, with more alteration occurring deeper in the basin and less towards to edges. Alteration of pre-existing dolomite by an Mg isotope exchange mechanism requires large supplies of Mg. This fact, and the direction of fluid flow point to the underlying crystalline basement as an important Mg source for dolomitization in the Williston Basin. We speculate that ultramafic rocks supplied the Mg when they were hydrothermally altered during an anomalous heating event that affected the Williston Basin in the Late Devonian / Early Carboniferous. Seismic pumping drove these fluids upwards through fracture networks in the crust, and eventually into the deep center of the Williston Basin along vertically-oriented down-to-the-basement faults, pressuring deep aquifers like the Winnipegosis, which resulted in fluids migrating up the stratigraphic dip of the beds towards the edges of the basin. Plotting and contouring the dolomite d26Mg data on a map of the Williston Basin reveals the pathways taken by the migrating fluids. One prominent pathway extending outwards from the center of the basin towards the northeast bears a strong resemblance to petroleum migration fairways deduced for the underlying Late Ordovician Red River Formation. This study demonstrates that Mg isotopic mapping of large dolomite bodies is a useful tool for deciphering fluid flow histories in sedimentary basins and may assist in locating migrated petroleum and Mississippi Valley Type ore deposits.
Magnesium Isotopes, Strontium Isotopes, Dolomitiza↗
Integrated Task And Data Parallel Programming: Language Design
his research investigates the combination of task and data parallel language constructs within a single programming language. There are an number of applications that exhibit properties which would be well served by such an integrated language. Examples include global climate models, aircraft design problems, and multidisciplinary design optimization problems. Our approach incorporates data parallel language constructs into an existing, object oriented, task parallel language. The language will support creation and manipulation of parallel classes and objects of both types (task parallel and data parallel). Ultimately, the language will allow data parallel and task parallel classes to be used either as building blocks or managers of parallel objects of either type, thus allowing the development of single and multi-paradigm parallel applications. 1995 Research Accomplishments In February I presented a paper at Frontiers '95 describing the design of the data parallel language subset. During the spring I wrote and defended my dissertation proposal. Since that time I have developed a runtime model for the language subset. I have begun implementing the model and hand-coding simple examples which demonstrate the language subset. I have identified an astrophysical fluid flow application which will validate the data parallel language subset. 1996 Research Agenda Milestones for the coming year include implementing a significant portion of the data parallel language subset over the Legion system. Using simple hand-coded methods, I plan to demonstrate (1) concurrent task and data parallel objects and (2) task parallel objects managing both task and data parallel objects. My next steps will focus on constructing a compiler and implementing the fluid flow application with the language. Concurrently, I will conduct a search for a real-world application exhibiting both task and data parallelism within the same program m. Additional 1995 Activities During the fall I collaborated with Andrew Grimshaw and Adam Ferrari to write a book chapter which will be included in Parallel Processing in C++ edited by Gregory Wilson. I also finished two courses, Compilers and Advanced Compilers, in 1995. These courses complete my class requirements at the University of Virginia. I have only my dissertation research and defense to complete.
Grimshaw, Andrew S.↗
Integrated Task and Data Parallel Programming
This research investigates the combination of task and data parallel language constructs within a single programming language. There are an number of applications that exhibit properties which would be well served by such an integrated language. Examples include global climate models, aircraft design problems, and multidisciplinary design optimization problems. Our approach incorporates data parallel language constructs into an existing, object oriented, task parallel language. The language will support creation and manipulation of parallel classes and objects of both types (task parallel and data parallel). Ultimately, the language will allow data parallel and task parallel classes to be used either as building blocks or managers of parallel objects of either type, thus allowing the development of single and multi-paradigm parallel applications. 1995 Research Accomplishments In February I presented a paper at Frontiers 1995 describing the design of the data parallel language subset. During the spring I wrote and defended my dissertation proposal. Since that time I have developed a runtime model for the language subset. I have begun implementing the model and hand-coding simple examples which demonstrate the language subset. I have identified an astrophysical fluid flow application which will validate the data parallel language subset. 1996 Research Agenda Milestones for the coming year include implementing a significant portion of the data parallel language subset over the Legion system. Using simple hand-coded methods, I plan to demonstrate (1) concurrent task and data parallel objects and (2) task parallel objects managing both task and data parallel objects. My next steps will focus on constructing a compiler and implementing the fluid flow application with the language. Concurrently, I will conduct a search for a real-world application exhibiting both task and data parallelism within the same program. Additional 1995 Activities During the fall I collaborated with Andrew Grimshaw and Adam Ferrari to write a book chapter which will be included in Parallel Processing in C++ edited by Gregory Wilson. I also finished two courses, Compilers and Advanced Compilers, in 1995. These courses complete my class requirements at the University of Virginia. I have only my dissertation research and defense to complete.
Grimshaw, A. S.↗
Developing a system for in vivo imaging of maize roots containing iodinated contrast media in soil using synchrotron XCT and XRF
We sought to develop a novel experimental system which enabled application of iodinated contrast media to in vivo plant roots intact in soil and was compatible with time-resolved synchrotron X-ray computed tomography imaging. The system was developed to overcome issues of low contrast to noise within X-ray computed tomography images of plant roots and soil environments, the latter of which can complicate image processing and result in the loss of anatomical information. To demonstrate the efficacy of the system we employ the novel use of both synchrotron X-ray computed tomography and synchrotron X-ray fluorescence mapping to capture the translocation of the contrast media through root vasculature into the leaves. With the application of contrast media we identify fluid flow in root vasculature and visualise anatomical features, which are otherwise often only observable in ex vivo microscopy, including: the xylem, metaxylem, pith, fibres in aerenchyma and leaf venation. We are also able to observe interactions between aerenchyma cross sectional area and solute transport in the root vasculature with depth. Our novel system was capable of successfully delivering sufficient contrast media into root and leaf tissues such that anatomical features could be visualised and internal fluid transport observed. We propose that our system could be used in future to study internal plant transport mechanisms and parameterise models for fluid flow in plants.
59 BASIC BIOLOGICAL SCIENCES↗
Multidisciplinary System Reliability Analysis
The objective of this study is to develop a new methodology for estimating the reliability of engineering systems that encompass multiple disciplines. The methodology is formulated in the context of the NESSUS probabilistic structural analysis code, developed under the leadership of NASA Glenn Research Center. The NESSUS code has been successfully applied to the reliability estimation of a variety of structural engineering systems. This study examines whether the features of NESSUS could be used to investigate the reliability of systems in other disciplines such as heat transfer, fluid mechanics, electrical circuits etc., without considerable programming effort specific to each discipline. In this study, the mechanical equivalence between system behavior models in different disciplines are investigated to achieve this objective. A new methodology is presented for the analysis of heat transfer, fluid flow, and electrical circuit problems using the structural analysis routines within NESSUS, by utilizing the equivalence between the computational quantities in different disciplines. This technique is integrated with the fast probability integration and system reliability techniques within the NESSUS code, to successfully compute the system reliability of multidisciplinary systems. Traditional as well as progressive failure analysis methods for system reliability estimation are demonstrated, through a numerical example of a heat exchanger system involving failure modes in structural, heat transfer and fluid flow disciplines.
Mahadevan, Sankaran↗
Multi-Disciplinary System Reliability Analysis
The objective of this study is to develop a new methodology for estimating the reliability of engineering systems that encompass multiple disciplines. The methodology is formulated in the context of the NESSUS probabilistic structural analysis code developed under the leadership of NASA Lewis Research Center. The NESSUS code has been successfully applied to the reliability estimation of a variety of structural engineering systems. This study examines whether the features of NESSUS could be used to investigate the reliability of systems in other disciplines such as heat transfer, fluid mechanics, electrical circuits etc., without considerable programming effort specific to each discipline. In this study, the mechanical equivalence between system behavior models in different disciplines are investigated to achieve this objective. A new methodology is presented for the analysis of heat transfer, fluid flow, and electrical circuit problems using the structural analysis routines within NESSUS, by utilizing the equivalence between the computational quantities in different disciplines. This technique is integrated with the fast probability integration and system reliability techniques within the NESSUS code, to successfully compute the system reliability of multi-disciplinary systems. Traditional as well as progressive failure analysis methods for system reliability estimation are demonstrated, through a numerical example of a heat exchanger system involving failure modes in structural, heat transfer and fluid flow disciplines.
Mahadevan, Sankaran↗
Connect the Dots: In Situ 4-D Seismic Monitoring of CO 2 Storage With Spatio-Temporal CNNs
4-D seismic imaging has been widely used in CO 2 sequestration projects to monitor the fluid flow in the volumetric subsurface region that is not sampled by wells. Ideally, real-time monitoring and near-future forecasting would provide site operators with great insights to understand the dynamics of the subsurface reservoir and assess any potential risks. However, due to obstacles such as high deployment cost, availability of acquisition equipment, exclusion zones around surface structures, only very sparse seismic imaging data can be obtained during monitoring. That leads to an unavoidable and growing knowledge gap over time. The operator needs to understand the fluid flow throughout the project lifetime and the seismic data are only available at a limited number of times. This is insufficient for understanding reservoir behavior. To overcome those challenges, we have developed spatio-temporal neural-network-based models that can produce high-fidelity interpolated or extrapolated images effectively and efficiently. Specifically, our models are built on an autoencoder, and incorporate the long short-term memory (LSTM) structure with a new loss function regularized by optical flow. We validate the performance of our models using real 4-D post-stack seismic imaging data acquired at the Sleipner CO 2 sequestration field. We employ two different strategies in evaluating our models. Numerically, we compare our models with different baseline approaches using classic pixel-based metrics. We also conduct a blind survey and collect a total of 20 responses from domain experts to evaluate the quality of data generated by our models. Finally, via both numerical and expert evaluation, we conclude that our models can produce high-quality 2-D/3-D seismic imaging data at a reasonable cost, offering the possibility of real-time monitoring or even near-future forecasting of the CO 2 storage reservoir.
4-D seismic imaging↗
Measuring higher-order moments of neutron-time-of-flight data for cryogenic inertial confinement fusion implosions on OMEGA
Ion temperatures serve as an important diagnostic for inertial confinement fusion (ICF) implosions. In direct-drive ICF experiments on OMEGA, neutron-time-of-flight (nTOF) data are used to infer the ion temperature of the fusing plasma produced in the implosion experiment. The analysis of the nTOF data requires an assumption about the shape of the underlying source signal. Since the source nTOF signal is a near-replica of the neutron energy spectrum, an ideal Gaussian shape, corresponding to the neutron energy spectrum of a uniform temperature plasma, is routinely employed. However, spatial and temporal variations of the ion temperature in the plasma give rise to higher-order moments, which were first described by Munro [Nucl. Fusion 56, 036001 (2016)]. In this work, we show a simpler alternative analysis to derive moments of the neutron energy spectrum for a plasma with variations in ion temperature. We also present a revised analysis of measured nTOF signals that uses a model with an additional degree of freedom to take into account the effect of ion temperature variations on the shape of the spectrum. Compared to presently used nTOF analysis, the revised analysis yields on average ≈2× more accurate fits to the data and up to 15% higher ion temperatures for cryogenic experiments. As a result, we quantify the ion temperature inflation caused by radially symmetric fluid flows, which are present even in a symmetric implosion, and which serve as a lower bound on the ion temperature inflation in real implosions.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Evolution of anisotropy with saturation and its implications for the elastoplastic responses of clay rocks
Many clay rocks have distinct bedding planes. Experimental studies have shown that their mechanical properties evolve with the degree of saturation (DOS), often with higher stiffness and strength after drying. For transversely isotropic rocks, the effects of saturation can differ between the bed-normal (BN) and bed-parallel (BP) directions, which gives rise to saturation-dependent stiffness and strength anisotropy. Accurate prediction of the mechanical behavior of clay rocks under partially saturated conditions requires numerical models that can capture the evolving elastic and plastic anisotropy with DOS. In this study, we present an anisotropy framework for coupled solid deformation-fluid flow in unsaturated elastoplastic media. Here we incorporate saturation-dependent strength anisotropy into an anisotropic modified Cam-Clay (MCC) model and consider the evolving anisotropy in both the elastic and plastic responses. The model was calibrated using experimental data from triaxial tests to demonstrate its capability in capturing strength anisotropy at various levels of saturation. Through numerical simulations, we demonstrate the role of evolving stiffness and strength anisotropy in the mechanical behavior of clay rocks. Plane strain simulations of triaxial compression tests were also conducted to demonstrate the impacts of material anisotropy and DOS on the mechanical and fluid flow responses.
36 MATERIALS SCIENCE↗
Depressurization-Induced Fines Migration in Sediments Containing Methane Hydrate: X-Ray Computed Tomography Imaging Experiments
Depressurization of hydrate-bearing sediments (HBS) can cause the movement of fine particles, and in turn, such fines migration affects fluid flow and mechanical behavior of sediments, ultimately affecting long-term hydrocarbon production and wellbore stability. In this work, we investigate how and to what extent depressurization of HBS causes fines migration using X-ray computed tomography (CT) imaging. Methane hydrate was synthesized in sediments with 10% fines content (FC), composed of sands with silt and/or clay, and the hydrate-bearing samples were stepwisely depressurized while acquiring CT images. The CT images were analyzed to quantify the spatial changes in FC in the host sediment and thus to capture the fines migration during depressurization. It was found that the FC changes began occurring from the hydrate dissociation regions. This confirms that the multiphase flow caused by depressurization accompanies fines migration. Depressurization of HBS with a hydrate saturation of ~20–40% caused FC reduction from ~10% to ~6–9%, and the extent of fines migration differed with the particle sizes of the host sands and the types of fines. It was found that fines migration was more pronounced with coarse sands and with silty fines. Such observed level of FC reduction is estimated to increase sediment permeability by several factors based on the Kozeny-type permeability model. In conclusion, our results support the notion that the extent of fines migration and its effect on fluid flow behavior need to be assessed in consideration of physical properties of host sediment and fine particles to identify optimum depressurization strategies.
42 ENGINEERING↗
Inlets, ducts and nozzles
The internal fluid mechanics research program in inlets, ducts, and nozzles is described. The program consists of a balanced effort between the development of computational tools and the conduct of experimental research. The experiments are designed to better understand the fluid flow physics, to develop new or improved flow models, and to provide benchmark quality data sets for validation of the computational methods. The inlet, duct, and nozzle research program is described according to three major classifications of flow phenomena: highly three-dimensional flow fields; shock-boundary layer interactions; and shear layer control. Specific examples of current and future elements of the research program are described for each of these phenomena. In particular, the highly three-dimensional flow field phenomena is highlighted by describing the computational and experiemental research program in transition ducts having a round-to-rectangular area variation. In the case of shock-boundary layer interactions, the specific details of research for normal shock-boundary layer interactions are described. For shear layer control research in vortex generators and the use of aerodynamic excitation for enhancement of the jet mixing process are described. Future research in inlets, ducts, and nozzles will include more emphasis on three-dimensional full Navier-Stokes methods and corresponding experiments designed to concentrate on the appropriate three-dimensional fluid flow physics.
Abbott, John M.↗
An integrated development of the equations of motion for elastic hypersonic flight vehicles
An integrated, consistent analytical framework is developed for modeling the dynamics of elastic hypersonic flight vehicles. A Lagrangian approach is used in order to capture the dynamics of rigid-body motion, elastic deformation, fluid flow, rotating machinery, wind, and a spherical rotating earth model, and to account for their interactions with each other. A vector form of the force, moment and elastic-deformation equations is developed from Lagrange's equation; a useable scalar form of these equations is also presented. The appropriate kinematic equations are developed, and are presented in a useable form. A preliminary study of the significance of selected terms in the equations of motion is conducted. Using generic data for a single-stage-to-orbit vehicle, it was found that the Coriolis force can reach values of up to 6 percent of the vehicle weight, and that the forces and moments attributable to fluid-flow terms can be significant.
Bilimoria, Karl D.↗
Effective permeability of an immiscible fluid in porous media determined from its geometric state
Based on the phenomenological extension of Darcy's law, two-fluid flow is dependent on a relative permeability function of saturation only that is process or path dependent with an underlying dependency on pore structure and wettability. For applications (fuel cells to underground CO 2 storage), it is imperative to determine the effective phase permeability relationships where the traditional approach is based on the inverse modeling of time-consuming experiments. The underlying reason is that the fundamental upscaling step from pore to Darcy scale, which links the pore structure of the porous medium to the continuum hydraulic conductivities, is not solved. Herein, we develop an artificial neural network (ANN) that relies on fundamental geometrical relationships to determine the mechanical energy dissipation during creeping immiscible two-fluid flow. Finally, the developed ANN is based on a prescribed set of state variables based on physical insights that predict the effective permeability of 4500 unseen pore-scale geometrical states with R 2 = 0.98.
42 ENGINEERING↗