Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex fluids”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,513 records · Page 84

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↗

Environmental assessment of the proposed nonelectric application of geothermal resources at Desert Hot Springs, California

The paper presents an environmental analysis performed in evaluating various proposed geothermal demonstration projects at Desert Hot Springs. These are categorized in two ways: (1) indirect, or (2) direct uses. Among the former are greenhouses, industrial complexes, and car washes. The latter include aquaculture, a cascaded agribusiness system, and a mobile home park. Major categories of environmental impact covered are: (1) site, (2) construction of projects, and (3) the use of the geothermal source. Attention is also given to the disposal of the geothermal fluid after use. Finally, it is concluded that there are no major problems forseen for each project, and future objectives are discussed.

Rosenberg, L.↗

CARS system for turbulent flame measurements

Simultaneous nitrogen number density and rotational-vibrational temperatures were measured in a turbulent diffusion flame with a Coherent Anti-Stokes Raman Scattering (CARS) instrument. The fuel jet was diluted with nitrogen (20 percent by volume) to allow temperature measurements across the entire jet mixing region. These measurements were compared with fluid dynamics computations. The CARS system incorporated a neodymium YAG laser, an intensified silicon photodiode array detector, and unique dynamic range enhancement methods. Theoretical calculations were based on a parabolic Navier-Stokes computer code. The comparison of these techniques will aid their development in the study of complex flowfields.

Antcliff, R. R.↗

Tethered elevator: A unique opportunity for space processing

The latest fluid dynamic and material science experiments in the microgravity environment have emphasized the importance of the residual gravity level and of the g-jitter on fluid physics phenomena. The tethered elevator presents the possibility of providing variable g-levels (both steady and g-jitter) around a very low steady g-level (that can be realized when the elevator is near the center of mass of the space station-tether complex). When positioning a variable periodic oscillation to the payload a clean g-jitter disturbance can be obtained that would not be otherwise obtainable by other systems. These two possibilities make the elevator a facility to help resolve a number of still open questions that are preventing wider utilization of the space environment in the microgravity area.

Monti, R.↗

Future requirements of wind tunnels for Computational Fluid Dynamics code verification

The role of experiment in the development of Computational Fluid Dynamics (CFD) for external aerodynamic flow-field prediction is discussed. Requirements for code verification from two sources pacing the development of CFD are described: (1) the need for developing adequate flow modeling, and (2) the need for establishing confidence in the use of CFD for predicting complex flows. The type of data needed to meet these needs differs in detail and scope and leads to definite wind tunnel requirements. Recent examples of testing to assess and develop turbulence models, and to verify code development are used to establish future requirements of wind tunnel testing for CFD verification.

Marvin, J. G.↗

Prediction of prompt NO(x) in hydrocarbon air flames

The gas turbine industry is directing particular attention to very low NOx combustors, whether for aircraft or land based CCGT systems. These low NOx combustors frequently use liquid fuels or natural gas burning under very lean premixed conditions with air or under rich-lean conditions, although only the first case is studied here. In land based systems, diluted steam or nitrogen are sometimes injected into the combustion chambers to reduce flame temperature. The NOx emissions from such systems are the product of three chemical mechanisms which are interrelated: the hydrocarbon prompt NO, the thermal NO (extended Zeldovich mechanism), and the nitrous oxide route to NO. Formation of NO2 from NO also occurs, as well as emission of carbon monoxide and unburnt hydrocarbons. When the fuel-oxidant proportion decreases towards leaner conditions, flame temperatures are lowered, resulting in the total NO being reduced and the thermal-NO contribution greatly diminished to the benefit of the remaining two mechanisms of NO formation. While knowledge of the elementary reactions and their chemical kinetics concerning methane and simple hydrocarbons combustion has existed for a number of years, its use for computer modeling is limited to simple flow dynamics configurations. Nevertheless, understanding of such combusting flows under a wide range of experimental conditions allows for analogies or speculations with more complex actual systems. Such understanding can be achieved by means of one dimensional laminar premixed flame modeling, with a full chemical mechanism which incorporates the three routes of NO formation. Complementary to this understanding is the modeling of the actual combustion system using a full description of the fluid dynamics coupled with a reduced chemical scheme, which is then compared against the first model. The objective of this investigation is to evaluate the relative importance of the three mechanisms of NO formation in lean premixed methane-air combustion with increasing pressure using the one dimensional plug flow package PREMIX, and to test the validity of a three dimensional model with a global chemical mechanism against the one dimensional model in the atmospheric pressure case. Methane is chosen because it is the only mechanism which is reasonably well known and is a good guide to the behavior of other hydrocarbons. The mixture ratio chosen is richer than that in lean gas turbines, but the combustion of this mixture with the low preheat gives realistic gas turbine final flame temperatures. Conditions of NO2 formation are also analyzed in the one dimensional model and results extrapolated to the case of gas turbines.

Valerie Dupont↗

A Class of High-Resolution Explicit and Implicit Shock-Capturing Methods

The development of shock-capturing finite difference methods for hyperbolic conservation laws has been a rapidly growing area for the last decade. Many of the fundamental concepts, state-of-the-art developments and applications to fluid dynamics problems can only be found in meeting proceedings, scientific journals and internal reports. This paper attempts to give a unified and generalized formulation of a class of high-resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock waves, perfect gases, equilibrium real gases and nonequilibrium flow computations. These numerical methods are formulated for the purpose of ease and efficient implementation into a practical computer code. The various constructions of high-resolution shock-capturing methods fall nicely into the present framework and a computer code can be implemented with the various methods as separate modules. Included is a systematic overview of the basic design principle of the various related numerical methods. Special emphasis will be on the construction of the basic nonlinear, spatially second and third-order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and flux-vector splitting approaches. Generalization of these methods to efficiently include real gases and large systems of nonequilibrium flows will be discussed. Some perbolic conservation laws to problems containing stiff source terms and terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for one-, two- and three-dimensional gas-dynamics problems. The use of the Lax-Friedrichs numerical flux to obtain high-resolution shock-capturing schemes is generalized. This method can be extended to nonlinear systems of equations without the use of Riemann solvers or flux-vector splitting approaches and thus provides a large savings for multidimensional, equilibrium real gases and nonequilibrium flow computations.

Yee, H. C.↗

Numerical Modeling of Three-Dimensional Fluid Flow with Phase Change

We present a numerical method to compute phase change dynamics of three-dimensional deformable bubbles. The full Navier-Stokes and energy equations are solved for both phases by a front tracking/finite difference technique. The fluid boundary is explicitly tracked by discrete points that are connected by triangular elements to form a front that is used to keep the stratification of material properties sharp and to calculate the interfacial source terms. Two simulations are presented to show robustness of the method in handling complex phase boundaries. In the first case, growth of a vapor bubble in zero gravity is studied where large volume increase of the bubble is managed by adaptively increasing the front resolution. In the second case, growth of a bubble under high gravity is studied where indentation at the rear of the bubble results in a region of large curvature which challenges the front tracking in three dimensions.

Esmaeeli, Asghar↗

Runtime and Pressurization Analyses of Propellant Tanks

Multi-element unstructured CFD has been utilized at NASA SSC to carry out analyses of propellant tank systems in different modes of operation. The three regimes of interest at SSC include (a) tank chill down (b) tank pressurization and (c) runtime propellant draw-down and purge. While tank chill down is an important event that is best addressed with long time-scale heat transfer calculations, CFD can play a critical role in the tank pressurization and runtime modes of operation. In these situations, problems with contamination of the propellant by inclusion of the pressurant gas from the ullage causes a deterioration of the quality of the propellant delivered to the test article. CFD can be used to help quantify the mixing and propellant degradation. During tank pressurization under some circumstances, rapid mixing of relatively warm pressurant gas with cryogenic propellant can lead to rapid densification of the gas and loss of pressure in the tank. This phenomenon can cause serious problems during testing because of the resulting decrease in propellant flow rate. With proper physical models implemented, CFD can model the coupling between the propellant and pressurant including heat transfer and phase change effects and accurately capture the complex physics in the evolving flowfields. This holds the promise of allowing the specification of operational conditions and procedures that could minimize the undesirable mixing and heat transfer inherent in propellant tank operation. It should be noted that traditional CFD modeling is inadequate for such simulations because the fluids in the tank are in a range of different sub-critical and supercritical states and elaborate phase change and mixing rules have to be developed to accurately model the interaction between the ullage gas and the propellant. We show a typical run-time simulation of a spherical propellant tank, containing RP-1 in this case, being pressurized with room-temperature nitrogen at 540 R. Nitrogen, shown in blue on the right-hand side of the figures, enters the tank from the diffuser at the top of the figures and impinges on the RP-1, shown in red, while the propellant is being continuously drained at the rate of 1050 lbs/sec through a pipe at the bottom of the tank. The sequence of frames in Figure 1 shows the resultant velocity fields and mixing between nitrogen and RP-1 in a cross-section of the tank at different times. A vortex is seen to form in the incoming nitrogen stream that tends to entrain propellant, mixing it with the pressurant gas. The RP-1 mass fraction contours in Figure 1 are also indicative of the level of mixing and contamination of the propellant. The simulation is used to track the propagation of the pure propellant front as it is drawn toward the exit with the evolution of the mixing processes in the tank. The CFD simulation modeled a total of 10 seconds of run time. As is seen from Figure 1d, after 5.65 seconds the propellant front is nearing the drain pipe, especially near the center of the tank. Behind this pure propellant front is a mixed fluid of compromised quality that would require the test to end when it reaches the exit pipe. Such unsteady simulations provide an estimate of the time that a high-quality propellant supply to the test article can be guaranteed at the modeled mass flow rate. In the final paper, we will discuss simulations of the LOX and propellant tanks at NASA SSC being pressurized by an inert ullage. Detailed comparisons will be made between the CFD simulations and lower order models as well as with test data. Conditions leading to cryo collapse in the tank will also be identified.

Field, Robert E.↗

Training Knowledge Bots for Physics-Based Simulations Using Artificial Neural Networks

Millions of complex physics-based simulations are required for design of an aerospace vehicle. These simulations are usually performed by highly trained and skilled analysts, who execute, monitor, and steer each simulation. Analysts rely heavily on their broad experience that may have taken 20-30 years to accumulate. In addition, the simulation software is complex in nature, requiring significant computational resources. Simulations of system of systems become even more complex and are beyond human capacity to effectively learn their behavior. IBM has developed machines that can learn and compete successfully with a chess grandmaster and most successful jeopardy contestants. These machines are capable of learning some complex problems much faster than humans can learn. In this paper, we propose using artificial neural network to train knowledge bots to identify the idiosyncrasies of simulation software and recognize patterns that can lead to successful simulations. We examine the use of knowledge bots for applications of computational fluid dynamics (CFD), trajectory analysis, commercial finite-element analysis software, and slosh propellant dynamics. We will show that machine learning algorithms can be used to learn the idiosyncrasies of computational simulations and identify regions of instability without including any additional information about their mathematical form or applied discretization approaches.

Samareh, Jamshid A.↗

NASA,FAA,ONERA Swept-Wing Icing and Aerodynamics: Summary of Research and Current Status

NASA, FAA, ONERA, and other partner organizations have embarked on a significant, collaborative research effort to address the technical challenges associated with icing on large scale, three-dimensional swept wings. These are extremely complex phenomena important to the design, certification and safe operation of small and large transport aircraft. There is increasing demand to balance trade-offs in aircraft efficiency, cost and noise that tend to compete directly with allowable performance degradations over an increasing range of icing conditions. Computational fluid dynamics codes have reached a level of maturity that they are being proposed by manufacturers for use in certification of aircraft for flight in icing. However, sufficient high-quality data to evaluate their performance on iced swept wings are not currently available in the public domain and significant knowledge gaps remain.

Aircraft Icing↗

Initial Low-Reynolds Number Iced Aerodynamic Performance for CRM Wing

NASA, FAA, ONERA, and other partner organizations have embarked on a significant, collaborative research effort to address the technical challenges associated with icing on large scale, three-dimensional swept wings. These are extremely complex phenomena important to the design, certification and safe operation of small and large transport aircraft. There is increasing demand to balance trade-offs in aircraft efficiency, cost and noise that tend to compete directly with allowable performance degradations over an increasing range of icing conditions. Computational fluid dynamics codes have reached a level of maturity that they are being proposed by manufacturers for use in certification of aircraft for flight in icing. However, sufficient high-quality data to evaluate their performance on iced swept wings are not currently available in the public domain and significant knowledge gaps remain.

Ice Formation↗

hydro-hifem v.1.0

SAND2021-3180 O The hydro-hifem 1.0 software is used to simulate transient physical processes in large geologic environments consisting of complex fractures and man-made infrastructures with affordable computational cost and mesh refinement. The Earth models with such complex model features can be designed by considering the reduction of insignificant dimensions of fractures/infrastructures and embracing their corresponding hierarchical material properties, and their computational meshes are generated by using Cubit. The finite element solutions are automatically distributed across nodes of parallel cluster for efficient computation of transient fluid flow or heat conduction in geologic media. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Beskardes, GungorDidem↗

General overview of multiphysics modeling and simulation at CEA/DES/IRESNE supporting safety, operating and development of nuclear reactors and facilities

This paper presents a general overview of the innovative multiphysics research/development activities carried on by the CEA/DES/IRESNE institute. Eight main different multiphysics fields addressed by IRESNE have been listed (among others: fuel fabrication, behavior of fuels, of reactors, of corium, radionuclide transfer to environment,...) in relation with the involved coupled physics: neutron transport, fluid mechanics, electro-magnetism, heat transfer, solid mechanics and (physico-)chemistry. The major mission of IRESNE institute is to study current and future nuclear systems integrated into a low carbon energy system. To fulfill this mission, IRESNE's activities rely on the development and the implementation of multiphysics modeling, advanced coupling scheme, data transfer techniques and dedicated multiphysics experimentations. Indeed, such multiphysics simulation approach of nuclear systems is today required to deal with the complex intrinsic features of studied phenomena. This approach aims at taking into account physical phenomena interactions, at each interest modeling scale, in order to reach an accurate and predictive behavior representation of a system component or of the whole system. The final objective is to contribute to the development of digital reactor, core and fuel at CEA, by progressing on applied mathematics, high-performance computing and physical mechanism understanding.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Microseismic at HFTS2: A Story of Three Stimulated Wells

A microseismic data set was acquired at Hydraulic Fracturing Test Site 2 (HFTS2) in the Delaware Basin (Texas, USA) by the HFTS2 consortium as a part of data-collection program designed to help maximize the efficiency of hydraulic fracturing, understand the impacts of nearby depleted zones, and increase the production of hydrocarbons. Advanced methods of microseismic monitoring applied at HFTS2 include 1) native synchronization of five monitor arrays, 2) real-time location of microseismic events, 3) joint event hypocenter/ velocity-model inversion that locates microseismic events simultaneously with constructing azimuthally anisotropic velocity models, 4) high-resolution reservoir imaging, utilizing microseismic events as downhole sources. HFTS2 microseismic data set illustrated the influence of geologic complexity, variations in treatment designs, and the presence of parent wells on the efficiency of hydraulic fracturing. Three Wolfcamp wells, Boxwood 1H, 2H, and 4H, were stimulated sequentially, exhibiting very different patterns of microseismicity. When the Boxwood 1H well is stimulated, its toe stages quickly reacted to the presence of depleted Bitterroot section to the east of the Boxwood pad, extending microseismic events over 3,000 ft to the east. This behavior changes dramatically after the treatment passes the Bitterroot section at the toe, and the heel stages in 1H showed tight microseismic grouping around the stimulated well, indicative of treatment of virgin reservoir rock. The Boxwood 2H demonstrates this same pattern of behavior between the toe and heel stages of the well, suggesting that the depleted Bitterroot section remains easier to refracture than adjacent non-depleted rock. Microseismic activity associated with Boxwood 4H treatment is distinctly different from those in Boxwood wells 1H and 2H. The events generated during stimulation of 4H are hardly influenced by the depleted Bitterroot section and fill the space between the previously stimulated 1H and 2H wells. Altogether, microseismic patterns suggest that treatment fluids tend to find the weakest passages through geologic formations, and the relative well positions and treatment sequence have a significant impact on fracture geometry, as indicated by microseismic event locations.

58 GEOSCIENCES↗

[?]Nonlinear Issues in the Aerothermochemistry of Gases and Materials and the Associated Physics and Dynamics of Interfaces

Our research and technology are focused on nonlinear issues in the aerothermochemistry of gases and materials and the associated physics and dynamics of interfaces. Our program is now organized to aggressively support the NASA Aeronautics Enterprise so as to: (a) develop a new generation of environmentally compatible, economic subsonic aircraft; (b) develop the technology base for an economically viable and environmentally compatible high-speed civil transport; (c) develop the technology options for new capabilities in high-performance aircraft; (d) develop hypersonic technologies for air-breathing flight; and (e) develop advanced concepts, understanding of physical phenomena, and theoretical, experimental, and computational tools for advanced aerospace systems. The implications from our research for aeronautical and aerospace technology have been both broad and deep. For example, using advanced computational techniques, we have determined exact solutions for the Schrodinger equation in electron-molecule scattering allowing us to evaluate atmospheric models important to reentry physics. We have also found a new class of exact solutions for the Navier Stokes equations. In experimental fluid dynamics, we have found explicit evidence of turbulence modification of droplet sizes in shock tube flow with condensation. We have developed a new diagnostic tool for the direct estimation of flow velocities at MHz sampling rates in quasi-one dimensional turbulent flow. This procedure suggests an unexpected confirmation of the possibility of 'natural' closure in Reynolds stresses with deep implications for the development of turbulent models. A transient increase is observed in both the spectral energy decay rate and the degree of chaotic complexity at the interface of a shock wave and a turbulent ionized gas. Even though the gas is apparently brought to rest by the shock wave, no evidence is found either of the expected relaminarization. A unique diamond-shaped nozzle has been designed for a detailed investigation of the effect of significant streamwise vorticity on the acoustic and IR characteristics of supersonic jets. Our results provide convincing evidence of the significant effect of vorticity on the far-field noise for the diamond jet as compared to the conventional round jets. We have found that the countercurrent shear layer mixes much more efficiently than conventional coflowing shear layers. We also developed the fluid thrust vectoring procedures which use counter flow to vector a jet. Our materials research has shown that the steep stress gradients at the fiber-matrix interface could be the primary cause of interface cracks after the processing of metallic and intermetallic matrix composites. New techniques have been evolved for: the microcharacterization of materials including microplastic strain and, point by point, the misorientation and plasticity for matrix composites; thermally induced stress measurements and load relaxation; the growth and characterization of metallic matrix composite interfaces; and for the growth of ferrite materials by pulsed laser deposition. The FAMU commitment to the HBCU Research Center also continues to be broad and deep.

Johnson, Joseph A., III↗

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗

Stochastic fracture generation and thermo-hydro-mechanical modeling in an equivalent continuum framework for enhanced geothermal systems

Enhanced geothermal systems (EGS) involve fracturing low permeability material to establish well connectivity and then injecting and circulating fluid into the fractured subsurface for geothermal power production. Changes in fracture aperture from contraction of the cooling matrix rock may alter network connectivity and risk thermal short-circuiting. Thermo-hydro-mechanical (THM) models are a useful tool to study these processes. However, as fracture networks are complex, and data may be limited, fracture networks in THM models are often stochastically generated. Given reliance on stochastic fracture networks and THM modeling to represent the subsurface and assess productivity of EGS, increased understanding of the influence of such statistically derived fracture networks on flow and heat transport in THM models is needed. Here, a new fracture process model is developed in the reactive transport code PFLOTRAN to stochastically generate fracture families and simulate changes in fracture aperture over time due to temperature changes of the rock matrix. Sixty-four different fracture networks ranging from well to poorly-connected, are modeled in PFLOTRAN with and without mechanical processes (THM vs TH). Results indicate that for well-connected fracture networks, thermal short-circuiting is less of a concern due to the abundance of available alternative flowpaths. For poorly-connected fracture networks, inclusion of mechanical processes showed steep thermal drawdown coincident with increase in fracture aperture along developing colder flowpaths, demonstrating the risk of thermal short-circuiting. Simulations with additional, larger fractures engineered to establish connectivity in a poorly-fractured subsurface, indicate that while stochastic variation of fracture orientation of the background network had limited influence, such variation in the engineered fractures significantly affected flow and heat transport.

Discrete fracture networks (DFN)↗