Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid Network”

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 451 records · Page 25

Orion Active Thermal Control System Dynamic Modeling Using Simulink/MATLAB

This paper presents dynamic modeling of the crew exploration vehicle (Orion) active thermal control system (ATCS) using Simulink (Simulink, developed by The MathWorks). The model includes major components in ATCS, such as heat exchangers and radiator panels. The mathematical models of the heat exchanger and radiator are described first. Four different orbits were used to validate the radiator model. The current model results were compared with an independent Thermal Desktop (TD) (Thermal Desktop, PC/CAD-based thermal model builder, developed in Cullimore & Ring (C&R) Technologies) model results and showed good agreement for all orbits. In addition, the Orion ATCS performance was presented for three orbits and the current model results were compared with three sets of solutions- FloCAD (FloCAD, PC/CAD-based thermal/fluid model builder, developed in C&R Technologies) model results, SINDA/FLUINT (SINDA/FLUINT, a generalized thermal/fluid network-style solver ) model results, and independent Simulink model results. For each case, the fluid temperatures at every component on both the crew module and service module sides were plotted and compared. The overall agreement is reasonable for all orbits, with similar behavior and trends for the system. Some discrepancies exist because the control algorithm might vary from model to model. Finally, the ATCS performance for a 45-hr nominal mission timeline was simulated to demonstrate the capability of the model. The results show that the ATCS performs as expected and approximately 2.3 lb water was consumed in the sublimator within the 45 hr timeline before Orion docked at the International Space Station.

Wang, Xiao-Yen J.↗

Virtual Flight through a Catalyst Particle

This visualization of the experimental data shows how scientists mapped the distribution of chemical elements in a single fluid catalytic cracking (FCC) particle and merged it with structural information about the pore networks. Because of the high resolution at which they mapped the catalyst, they were able to look deep into the pores and learn more about the metal poisoning reaction. The changing colors of the "fog" inside the pores reflect the changing chemistry.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Heat Transfer Experiments of a 1st Stage Blade Cascade for Supercritical CO2 Oxy-Combustion Turbine Application

The results of internally cooled 1st stage blade (S1B) cascade testing in a supercritical CO2 environment is presented. The turbine blade design has been previously established for the end application of an oxy-combustion turbine operating in the Allam-Fetvedt cycle with turbine inlet conditions of 305 bar and 1150°C. The internally cooled blade features leading edge (LE) region impingement cooling, mid-section ribbed serpentine passages, and a pin-finned trailing edge (TE) region before cooling ejection holes. The geometry for the tested blade cascade has a cooled central blade with un-cooled blades on either side to match flowpath areas of the actual turbine. The flowpath reuses internal components previously employed for mid-section region ribbed serpentine passage experiments that established Nusselt number enhancement ratios over a range of Reynolds numbers from 100,000-400,000. New components include flow conditioning plates upstream and downstream of the blade cascade to adequately represent the flow field and blade external heat transfer coefficient profiles for the actual turbine. The cooled central blade utilizes uniform crystal temperature sensors (UCTS) with six sensors each on the blade pressure and suction surfaces distributed radially and from LE to TE. The post-processed UCTS quantified the maximum wall temperature seen at each installed sensor location. The test procedure consisted of establishing supercritical CO2 cooling flow temperature and flow rate and maintaining it throughout the test. The flow rate aims to match that for the actual in-service turbine blade design and is maintained through an orifice restriction to keep the pressure differential between internal cooling flow and external hot flow nearly constant. For the sCO2 flow path external to the blade, temperatures were ramped throughout the test via control of the test loop’s natural gas burner heater. The maximum temperature seen was 468°C and held constant for a duration of 10 minutes at which the blade metal temperature was predicted to be at its maximum before ramping down. For the turbine blade design for service inlet conditions, external flow path computational fluid dynamics (CFD) results and an internal cooling 1-D thermal and hydraulic flow network model using experimentally validated correlations served as thermal finite element (FE) boundary conditions to predict blade metal temperatures. These predicted temperatures were subsequently utilized in a structural FE model to predict blade life ratings dictated by Haynes 282 creep strength data, having a strong dependence on temperature. The boundary conditions experienced during testing are used in the same workflow and compared to the experimental results, with the goal of validating the analysis methodology and providing insight on the uncertainty in local metal temperature predictions.

20 FOSSIL-FUELED POWER PLANTS↗

Hydrogeological Interpretation of Candidate Origin Sites for Martian Meteorite ALH84001

Barlow (this meeting) has identified two potential source craters for the martian meteorite ALH84001. The craters are at 11.7 deg S, 243.3 deg W (Mare Tyrrhenum site) and 14.0 deg S, 343.5 deg W (Sinus Sabaeus site). As noted by Barlow, both craters lie in the heavily cratered terrain (HCT) and are adjacent to fluvial valleys, Here I explore the fluvial history of these areas based upon the surrounding valley morphology. The most prominent valley network at the Sabaeus site is Evros Vallis. This wide, flat-floored valley is approximately 600 km long with an average width of 2.5 km and a depth of 220 m. The eroded volume of the entire Evros network is approximately 6 x 10(exp 11) cc. This is typical for networks located in the heavily cratered terrain (e.g. Warrego and Parana Valles). Evros is also an isolated valley system. No similar networks are found in the surrounding terrain. Thus it is unlikely that Evros formed as a result of widespread rainfall. A localized water source, such as discharge of a hydrothermal system or localized melting of snowfall, seems more consistent Previous modeling has demonstrated that only hydrothermal systems associated with high permeability subsurfaces can discharge sufficient water to form a valley network. The bulk of the discharge from such systems is consequently low temperature, slightly heated water Precipitation of calcium carbonate by low temperature fluids is consistent with most interpretations of the geochemistry of ALH84001. Available imagery at the Tyrrhenum site is of lesser quality. While eroded units of the HCT are nearby, there are no comparable well developed valley networks at this site. Erosion is instead manifested predominantly as gullies on slopes. This style of erosion suggests that water was not present at this site for the length of time as at the more integrated Sabaeus site. The superposition of fluidized ejecta blankets suggests however that ground water or ground ice was still present at this locality fairly recently, if the chosen crater is indeed correct.

Gulick, Virginia C.↗

Multi-scale computational model of three-dimensional hemodynamics within a deformable full-body arterial network

In this article, we present a computational multi-scale model of fully three-dimensional and unsteady hemodynamics within the primary large arteries in the human. Computed tomography image data from two different patients were used to reconstruct a nearly complete network of the major arteries from head to foot. A linearized coupled-momentum method for fluid–structure-interaction was used to describe vessel wall deformability and a multi-domain method for outflow boundary condition specification was used to account for the distal circulation. We demonstrated that physiologically realistic results can be obtained from the model by comparing simulated quantities such as regional blood flow, pressure and flow waveforms, and pulse wave velocities to known values in the literature. We also simulated the impact of age-related arterial stiffening on wave propagation phenomena by progressively increasing the stiffness of the central arteries and found that the predicted effects on pressure amplification and pulse wave velocity are in agreement with findings in the clinical literature. This work demonstrates the feasibility of three-dimensional techniques for simulating hemodynamics in a full-body compliant arterial network.

97 MATHEMATICS AND COMPUTING↗

Pore-scale and Continuum Simulations of Solute Transport Micromodel Benchmark Experiments

Four sets of micromodel nonreactive solute transport experiments were conducted with flow velocity, grain diameter, pore-aspect ratio, and flow focusing heterogeneity as the variables. The data sets were offered to pore-scale modeling groups to test their simulators. Each set consisted of two learning experiments, for which all results was made available, and a challenge experiment, for which only the experimental description and base input parameters were provided. The experimental results showed a nonlinear dependence of the dispersion coefficient on the Peclet number, a negligible effect of the pore-aspect ratio on transverse mixing, and considerably enhanced mixing due to flow focusing. Five pore-scale models and one continuum-scale model were used to simulate the experiments. Of the pore-scale models, two used a pore-network (PN) method, two others are based on a lattice-Boltzmann (LB) approach, and one employed a computational fluid dynamics (CFD) technique. The learning experiments were used by the PN models to modify the standard perfect mixing approach in pore bodies into approaches to simulate the observed incomplete mixing. The LB and CFD models used these experiments to appropriately discretize the grid representations. The continuum model use published non-linear relations between transverse dispersion coefficients and Peclet numbers to compute the required dispersivity input values. Comparisons between experimental and numerical results for the four challenge experiments show that all pore-scale models were all able to satisfactorily simulate the experiments. The continuum model underestimated the required dispersivity values and, resulting in less dispersion. The PN models were able to complete the simulations in a few minutes, whereas the direct models needed up to several days on supercomputers to resolve the more complex problems.

benchmark↗

Resonant amplification of enzymatic chemical oscillations by oscillating flow

Using theory and simulation, we analyzed the resonant amplification of chemical oscillations that occur due to externally imposed oscillatory fluid flows. The chemical reactions are promoted by two enzyme-coated patches located sequentially on the inner surface of a pipe that transports the enclosed chemical solution. In the case of diffusion-limited systems, the period of oscillations in chemical reaction networks is determined by the rate of the chemical transport, which is diffusive in nature and, therefore, can be effectively accelerated by the imposed fluid flows. In this work, we first identify the natural frequencies of the chemical oscillations in the unperturbed reaction–diffusion system and, then, use the frequencies as a forcing input to drive the system to resonance. We demonstrate that flow-induced resonance can be used to amplify the amplitude of the chemical oscillations and to synchronize their frequency to the external forcing. In particular, we show that even 10% perturbations in the flow velocities can double the amplitude of the resulting chemical oscillations. Particularly, effective control can be achieved for the two-step chemical reactions where during the first half-period, the fluid flow accelerates the chemical flux toward the second catalytic patch, while during the second half-period, the flow amplifies the flux to the first patch. The results can provide design rules for regulating the dynamics of coupled reaction–diffusion processes and can facilitate the development of chemical reaction networks that act as chemical clocks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine Learning Assisted Safety Modeling and Analysis of Advanced Reactors

With the advances in computational power and numerical methods, analysts can now rely on first-principle simulations to predict ultra-fine details in a variety of applications. Advances in machine learning (ML) have produced algorithms that can now learn high-level abstractions via hierarchical models. This project aims to leverage advances in ML techniques and the available high-resolution simulation data to develop a novel modeling and simulation (M\&S) methodology for reactor safety analysis. While application-agnostic ML techniques are available, complex physics constraints need to be incorporated into ML techniques to build ML-based closures for computationally efficient predictive simulations. This project intends to develop a physics-guided data-driven multi-scale methodology for M\&S of advanced reactors. The project focuses on thermal fluid (T/F) phenomena, which play major roles in advanced reactor safety. Specifically, we propose a data-driven coarse-mesh turbulence model based on local flow features for the transient analysis of thermal mixing and stratification in a sodium-cooled fast reactor (SFR). The model has a coarse-mesh setup to ensure computational efficiency, while it is trained by fine-mesh computational fluid dynamics (CFD) data with Reynolds-averaged Navier-Stokes (RANS) turbulence model to ensure accuracy. Three different neural networks are developed and tested for loss-of-flow transients in the hot pool of SFR, i.e. the densely connected convolutional neural network (DCNN), long-short-term-memory network based on proper orthogonal decomposition (POD-LSTM), and the DCNN informed by LSTM (DCNN-LSTM). The performances of these three neural networks are evaluated based on baseline models. The DCNN-LSTM model has been chosen for further hyperparameter optimization. Furthermore, based on a simplified two-dimensional case, uncertainty quantification (UQ) of the developed ML-based closure are investigated with three methods, i.e. Monte Carlo dropout, deep ensemble, and Bayesian neural network. The developed ML-based turbulent viscosity closure relation based on deep ensemble is then integrated into the system analysis module SAM and serves as a term in the conservation equations. Such a SAM-ML based procedure guarantees that the obtained results are consistent with the physical constraints of the thermal-fluid system. The SAM-ML simulation on the same loss-of-flow transient showed comparable accuracy with the CFD simulation but with a much coarser mesh setup. Last but not least, the ML-based closure improvement with the support of higher-fidelity data from large eddy simulation (LES) is discussed. As a first step towards this direction, a baseline LES simulation is performed to obtain comparable data with RANS results. Based on the early results, future investigation on further improving the ML-based closure is discussed. We believe the developed approach that combines scientific machine learning with nuclear system analysis code can benefit the advanced reactor community as more accurate safety analyses will better characterize reactor safety margins and reduce licensing efforts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development of flow network analysis code for block type VHTR core by linear theory method

VHTR (Very High Temperature Reactor) is high-efficiency nuclear reactor which is capable of generating hydrogen with high temperature of coolant. PMR (Prismatic Modular Reactor) type reactor consists of hexagonal prismatic fuel blocks and reflector blocks. The flow paths in the prismatic VHTR core consist of coolant holes, bypass gaps and cross gaps. Complicated flow paths are formed in the core since the coolant holes and bypass gap are connected by the cross gap. Distributed coolant was mixed in the core through the cross gap so that the flow characteristics could not be modeled as a simple parallel pipe system. It requires lot of effort and takes very long time to analyze the core flow with CFD analysis. Hence, it is important to develop the code for VHTR core flow which can predict the core flow distribution fast and accurate. In this study, steady state flow network analysis code is developed using flow network algorithm. Developed flow network analysis code was named as FLASH code and it was validated with the experimental data and CFD simulation results. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Sequential geophysical and flow inversion to characterize fracture networks in subsurface systems

Subsurface applications, including geothermal, geological carbon sequestration, and oil and gas, typically involve maximizing either the extraction of energy or the storage of fluids. Fractures form the main pathways for flow in these systems, and locating these fractures is critical for predicting flow. However, fracture characterization is a highly uncertain process, and data from multiple sources, such as flow and geophysical are needed to reduce this uncertainty. We present a nonintrusive, sequential inversion framework for integrating data from geophysical and flow sources to constrain fracture networks in the subsurface. In this framework, we first estimate bounds on the statistics for the fracture orientations using microseismic data. These bounds are estimated through a combination of a focal mechanism (physics-based approach) and clustering analysis (statistical approach) of seismic data. Then, the fracture lengths are constrained using flow data. In conclusion, the efficacy of this inversion is demonstrated through a representative example.

97 MATHEMATICS AND COMPUTING↗

High Temperature, High Pressure Devices for Zonal Isolation in Geothermal Wells

The U.S. Department of Energy is leading the development of alternative energy sources that will ensure the long-term energy independence of our nation. One key renewable resource being advanced is geothermal energy which offers an environmentally benign, reliable source of energy for the nation. To utilize this resource, water will be introduced into wells 3 to 10 km deep to create a geothermal reservoir. This approach is known as an Enhanced Geothermal System (EGS). The high temperatures and pressures at these depths have become a limiting factor in the development of this energy source. For example, reliable zonal isolation for high-temperature applications at high differential pressures is needed to conduct mini-fracs and other stress state diagnostics. Zonal isolation is essential for many EGS reservoir development activities. To date, the capability has not been sufficiently demonstrated to isolate sections of the wellbore to: 1) enable stimulation; and 2) seal off unwanted flow regions in unknown EGS completion schemes and high-temperature (>200°C) environments. In addition, packers and other zonal isolation tools are required to eliminate fluid loss, to help identify and mitigate short circuiting of flow from injectors to producers, and to target individual fractures or fracture networks for testing and validating reservoir models. General-purpose open-hole packers do not exist for geothermal environments, with the primary barrier being the poor stability of elastomeric seals at high temperature above 175°C. Experimental packer systems have been developed for geothermal environments but they currently only operate at low pressure, they are not retrievable, and they are not commercially available. The development of the high-temperature, high-pressure (HTHP) zonal isolation device would provide the geothermal community with the capability to conduct mini-fracs, eliminate fluid loss, to help identify and mitigate short circuiting of flow from injectors to producers, and to target individual fractures or fracture networks for testing and validating reservoir models. The goal of this program, therefore, was to develop and demonstrate reliable high-temperature high-pressure zonal isolation devices that are compatible with the high-temperature downhole environment anticipated for Enhanced Geothermal Systems. Over the course of this 3 year program, CTD designed and demonstrated a high temperature high pressure zonal isolation device. CTD has utilized its expertise in high-temperature materials, shape memory composites and foams, and deployment mechanisms for spacecraft and oil tools to develop a new class of HTHP zonal isolation devices for use in EGS. Specifically, the objective was to demonstrate the ability to isolate sections of an EGS wellbore to: 1) enable stimulation; and 2) seal off unwanted flow regions in unknown EGS completion schemes and high-temperature (>200°C) environments at temperatures upwards of 300°C. A concept for a high-temperature, high-pressure zonal isolation device was designed for use in creating circulation paths in environments where temperatures are in excess of 300°C and differential pressures of up to 10,000 psi are required to stimulate solid granite. This new zonal isolation device distributes the high-pressure differential through a significant length of high temperature, shape memory polymer composite material. The development of differential pressure is initiated by a thermally triggered actuation of a ‘structural seal’ that fills the cross section of the bore with the shape memory material. After this preliminary seal is made, differential pressure begins to build through the material as the seal forces the flow to pass through pores in the shape memory material. The shape memory material is compressed into place to react the force of the internal

58 GEOSCIENCES↗

Application of a convolutional neural network in permeability prediction: A case study in the Jacksonburg-Stringtown oil field, West Virginia, USA

Permeability is a critical parameter for understanding subsurface fluid flow behavior, managing reservoirs, enhancing hydrocarbon recovery, and sequestering carbon dioxide. In general, permeability is measured in the laboratory based on subsurface core samples, calculated from well logs or estimated from well tests. However, laboratory measurements and well tests are expensive, time-consuming, and usually limited to a few core samples or wells in a hydrocarbon field or carbon storage site. Machine-learning techniques are good options for generating a rapid, robust, and cost-effective permeability prediction model because of their strengths to recognize the potential interrelationships between input and output variables. Convolutional neural networks (CNN), as a good pattern recognition algorithm, are widely used in image processing, natural language processing, and speech recognition, but are rarely used with regression problems and even less often in reservoir characterization. We have developed a CNN regression model to estimate the permeability in the Jacksonburg-Stringtown oil field, West Virginia, which is a potential carbon storage site and enhanced oil recovery operations field. We also evaluate the concept of the geologic feature image, which is converted from geophysical well logs. Five variables, including two commonly available conventional well logs (the gamma rays [GRs] and bulk density) and three well-log-derived variables (the slopes of the GR and bulk density curves, and shale content), are used to generate a geologic feature image. The CNN treats the geologic feature image as the input and the permeability as the desired output. In addition, the permeability predicted using traditional backpropagation artificial neural networks, which are optimized by genetic algorithms and particle swarm optimization, is compared with the permeability estimated using our CNN. Our results indicate that the CNN regression model provides more accurate permeability predictions than the traditional neural network.

Geochemistry & Geophysics↗

On the Versatility of Rheoreversible, Stimuli-responsive Hydraulic-Fracturing Fluids for Enhanced Geothermal Systems: Effect of Reservoir pH

Abstract The primary challenge for the feasibility of enhanced geothermal systems (EGS) is to cost-effectively create high-permeability reservoirs inside deep crystalline bedrock. Although fracturing fluids are commonly used for oil/gas, standard fracturing methods are not developed or proven for EGS temperatures and pressures. Furthermore, the environmental impacts of currently used fracturing methods are only recently being determined. These authors recently reported an environmentally benign, CO2-activated, rheoreversible fracturing fluid that enhances permeability through fracturing due to in situ volume expansion and gel formation. The potential of this novel fracturing fluid is evaluated in this work towards its application at geothermal sites under different pH conditions. Laboratory-scale fracturing experiments using Coso Geothermal rock cores under different pH environments were performed followed by X-ray microtomography characterization. The results demonstrate that CO2-reactive aqueous solutions of environmentally amenable polyallylamine (PAA) consistently and reproducibly creates/propagates fracture networks through highly impermeable crystalline rock from Coso EGS sites at considerably lower effective stress as compared to conventional fracturing fluids. In addition, permeability was significantly enhanced in a wide range of formation-water pH values. This effective, and environmentally-friendly fracturing fluid technology represents a potential alternative to conventional fracturing fluids.

CO2-reactive↗

On the Versatility of Rheoreversible, Stimuli-responsive Hydraulic-Fracturing Fluids for Enhanced Geothermal Systems: Effect of Reservoir pH

Abstract The primary challenge for the feasibility of enhanced geothermal systems (EGS) is to cost-effectively create high-permeability reservoirs inside deep crystalline bedrock. Although fracturing fluids are commonly used for oil/gas, standard fracturing methods are not developed or proven for EGS temperatures and pressures. Furthermore, the environmental impacts of currently used fracturing methods are only recently being determined. These authors recently reported an environmentally benign, CO2-activated, rheoreversible fracturing fluid that enhances permeability through fracturing due to in situ volume expansion and gel formation. The potential of this novel fracturing fluid is evaluated in this work towards its application at geothermal sites under different pH conditions. Laboratory-scale fracturing experiments using Coso Geothermal rock cores under different pH environments were performed followed by X-ray microtomography characterization. The results demonstrate that CO2-reactive aqueous solutions of environmentally amenable polyallylamine (PAA) consistently and reproducibly creates/propagates fracture networks through highly impermeable crystalline rock from Coso EGS sites at considerably lower effective stress as compared to conventional fracturing fluids. In addition, permeability was significantly enhanced in a wide range of formation-water pH values. This effective, and environmentally-friendly fracturing fluid technology represents a potential alternative to conventional fracturing fluids.

geothermal, reservoirs, crystalline, fracturing, o↗

Single and Multi-Node Modeling of Direct, Submerged, and Self-Pressurization of A Cryogenic Propellant Tank Using Nodal Tools

The pressurization of cryogenic propellant tanks will always be an important process so long as cryogenic liquids are being considered as fuel sources or used for other in-space applications. Pressure control of the tank ullage is necessary for managing propellant flowrates to an engine or a receiver tank, and modeling of the process is used to predict the pressurant requirements and the amount of propellant boiloff. Direct ullage pressurization is the more traditional approach to tank pressurization, as the physics are straight-forward, and ample test (flight) data have been collected and analyzed over the past several decades. Submerged injection pressurization is an alternate method for tank pressurization and has been shown to reduce pressurant requirements, subcool the propellant, and reduce the risk of ullage collapse. Additionally, the pressurant gas entering the ullage is usually much colder when using the submerged pressurization approach, resulting in reduced propellant boiloff. These benefits are at the expense of vaporizing a small percentage of the propellent. Both tank pressurization methods are viable options for current and future space missions, and it is important to have the capability of analyzing the tank ullage conditions for both approaches. Our previous work has demonstrated the development of a Generalized Fluid System Simulation Program (GFSSP) model, which contains a thermodynamic equilibrium heat and mass transfer subroutine capable of effectively analyzing both direct and submerged pressurization systems [1-2]. This subroutine has most recently been enhanced to include the non-equilibrium effect of pressurant dissolution into the propellant. To date the ullage has always been represented as a single node, and although the simulated single-node temperatures have good comparison with the volume-averaged ullage temperatures computed from test data, the physics of the thermal stratification in the ullage were never captured, and adjustment factors in the model were required. The purpose of this paper is to introduce the development of a multi-node ullage model using GFSSP and to discuss the improvements of the simulated ullage temperature distribution and its resulting effects on ullage heat transfer processes. Test data from the Cryogenic Propellant Storage and Transfer Engineering Developmental Unit (CPST EDU) was used for model validation. For additional comparison, a Thermal Desktop (TD) model was also developed to analyze the CPST EDU direct ullage pressurization tests using both a single node and multi-node approach. The model includes the direct pressurant line, vent line, fill/drain line, and a TD FloCAD Compartment. The TD FloCAD Compartment is employed to represent the liquid and ullage as single volumes inside the tank, to include a liquid/vapor interface, and to generate network level objects such as lumps (analogous to nodes in GFSSP), paths, and ties between the fluid and thermal elements. An established heat load on the model tank was leveraged from a pre-existing higher-fidelity model correlated to CPST EDU test data.

pressurization↗

Investigation of Stress Field and Fracture Development During Shale Maturation Using Analog Rock Systems

The emergence of hydrocarbons within shale as a major recoverable resource has sparked interest in fluid transport through these tight mudstones. Recent studies suggest the importance to recovery of microfracture networks that connect localized zones with large organic content to the inorganic matrix. This paper presents a joint modeling and experimental study to examine the onset, formation, and evolution of microfracture networks as shale matures. Both the stress field and fractures are simulated and imaged. A novel laboratory-scale, phase-field fracture propagation model was developed to characterize the material failure mechanisms that play a significant role during the shale maturation process. The numerical model developed consists of coupled solid deformation, pore pressure, and fracture propagation mechanisms. Benchmark tests were conducted to validate model accuracy. Laboratory-grade gelatins with varying Young’s modulus were used as scaled-rock analogs in two dimensional Hele-Shaw cell apparatus. Yeast within the gelatin generates gas in a fashion analogous to hydrocarbon formation as shale matures. These setups allow study and visualization of host rock elastic-brittle fracture and fracture network propagation mechanisms. The experimental setup was fitted to utilize photoelasticity principles coupled with birefringence properties of gelatin to explore visually the stress field of the gelatin as the fracture network developed. Stress optics image analysis and Linear Elastic Fracture Mechanics (LEFM) principles for crack propagation were used to monitor fracture growth for each gelatin type. Observed and simulated responses suggest gas diffusion within and deformation of the gelatin matrix as predominant mechanisms for energy dissipation depending on gelatin strength. LEFM, an experimental estimation of principal stress development with fracture growth, at different stages was determined for each gelatin rheology. Here, the interplay of gas diffusion and material deformation determines the resulting frequency and pattern of fractures. Results correlate with Young’s modulus. Experimental and computed stress fields reveal that fractures resulting from internal gas generation are similar to, but not identical to, type 1 opening mode.

58 GEOSCIENCES↗

Testing and Analytical Modeling for Purging Process of a Cryogenic Line

The purging operations for cryogenic main propulsion systems of upper stage are usually carried out for the following cases: 1) Purging of the Fill/Drain line after completion of propellant loading. This operation allows the removal of residual propellant mass; and 2) Purging of the Feed/Drain line if the mission is scrubbed. The lines would be purged by connections to a ground high-pressure gas storage source. The flowrate of purge gas should be regulated such that the pressure in the line will not exceed the required maximum allowable value. Exceeding the maximum allowable pressure may lead to structural damage in the line. To gain confidence in analytical models of the purge process, a test series was conducted. The test article, a 20-cm incline line, was filled with liquid hydrogen and then purged with gaseous helium (GHe). The influences of GHe flowrates and initial temperatures were evaluated. The Generalized Fluid System Simulation Program, an in-house general-purpose computer program for flow network analysis, was utilized to model and simulate the testing. The test procedures, modeling descriptions, and the results will be presented in the final paper.

Hedayat, A.↗

Testing and Analytical Modeling for Purging Process of a Cryogenic Line

The purging operations for cryogenic main propulsion systems of upper stage are usually carried out for the following cases: 1) Purging of the Fill/Drain line after completion of propellant loading. This operation allows the removal of residual propellant mass; and 2) Purging of the Feed/Drain line if the mission is scrubbed. The lines would be purged by connections to a ground high-pressure gas storage source. The flow-rate of purge gas should be regulated such that the pressure in the line will not exceed the required maximum allowable value. Exceeding the maximum allowable pressure may lead to structural damage in the line. To gain confidence in analytical models of the purge process, a test series was conducted. The test article, a 20-cm incline line, was filled with liquid hydrogen and then purged with gaseous helium (GHe). The influences of GHe flow-rates and initial temperatures were evaluated. The Generalized Fluid System Simulation Program, an in-house general-purpose computer program for flow network analysis, was utilized to model and simulate the testing. The test procedures, modeling descriptions, and the results will be presented in the final paper.

Hedayat, A.↗