SEARCH · Search NASA

Results for “Three dimensional computational fluid dynamics model”

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.

SourceResearch connectionFrom yearThrough year

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

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS2003

An Active Smart Material Control System for F/A-18 Buffet Alleviation

The vertical tail buffet problem of fighter aircraft occurs at high angles of attack when the vortical flow breaks down ahead of the vertical tails resulting in unsteady and unbalanced pressure loads on the vertical tails. The buffet loads imposed upon the vertical tails resulted in a premature fatigue failure of the tails, and consequently limits the performance and super maneuverability of twin-tail fighter aircraft. An active smart material control system using distributed piezoelectric actuators has been developed for buffet alleviation and is presented. The inboard and outboard surfaces of the vertical tail are equipped with piezoelectric actuators to control the buffet responses in the first bending and torsion modes. The electrodynamics of the piezoelectric actuators are expressed with a three-dimensional finite-element model. A single-input-single-output controller is designed to drive the active piezoelectric actuators. High-fidelity multidisciplinary analysis modules for the fluid dynamics, structure dynamics, electrodynamics of the piezoelectric actuators, control law, fluid structure interfacing, and grid motion are integrated into a multidisciplinary computing environment that controls the temporal synchronization of the analysis modules. At 30 degree angle of attack, RMS values of tip acceleration are reduced by as much as 12%. The peak values of the power spectral density of tail-tip acceleration are reduced by as much as 22% in the first bending mode and by as much as 82% in the first torsion mode. The actively controlled piezoelectric actuators were also effective in adding damping at wide range of angles of attack.

Sheta, Essam F.↗

NASA NTRS2014

Computational Aerodynamic Analysis of Three-Dimensional Ice Shapes on a NACA 23012 Airfoil

The present study identifies a process for performing computational fluid dynamic calculations of the flow over full three-dimensional (3D) representations of complex ice shapes deposited on aircraft surfaces. Rime and glaze icing geometries formed on a NACA23012 airfoil were obtained during testing in the NASA Glenn Research Center's Icing Research Tunnel (IRT). The ice shape geometries were scanned as a cloud of data points using a 3D laser scanner. The data point clouds were meshed using Geomagic software to create highly accurate models of the ice surface. The surface data was imported into Pointwise grid generation software to create the CFD surface and volume grids. It was determined that generating grids in Pointwise for complex 3D icing geometries was possible using various techniques that depended on the ice shape. Computations of the flow fields over these ice shapes were performed using the NASA National Combustion Code (NCC). Results for a rime ice shape for angle of attack conditions ranging from 0 to 10 degrees and for freestream Mach numbers of 0.10 and 0.18 are presented. For validation of the computational results, comparisons were made to test results from rapid-prototype models of the selected ice accretion shapes, obtained from a separate study in a subsonic wind tunnel at the University of Illinois at Urbana-Champaign. The computational and experimental results were compared for values of pressure coefficient and lift. Initial results show fairly good agreement for rime ice accretion simulations across the range of conditions examined. The glaze ice results are promising but require some further examination.

Aerodynamics↗

DOE OSTI2023

Coupling Thermal-Hydraulics with Reactor Pressure Vessel Fracture Models

Because of the importance of maintaining the integrity of the reactor pressure vessel (RPV) in light-water reactors, significant effort has been devoted over the last several decades to understand the potential failure mechanisms, characterize how material properties evolve due to environmental conditions, and develop tools to assess the probability of fracture in aged RPVs during potential transient events. One important aspect of this problem that has not yet been fully addressed is the effects of nonuniform cooling on the RPV inner wall due to the way low-temperature water injected into the primary coolant system is distributed within the RPV during potential loss-of-coolant accident transients. Colder temperatures are expected in the regions near the inlets, causing a “cold-plume” effect. Rapid cooling of the RPV wall increases the likelihood of fracture initiation, and localized cooling due to cold-plume effects could potentially exacerbate this problem. The Grizzly code can perform a three-dimensional simulation of the RPV thermomechanical response as well as the probabilistic assessment of fracture for a population of pre-existing flaws, which is necessary to account for spatial variations in coolant temperature. This report documents a first effort to couple Grizzly RPV fracture models with computational fluid dynamics models performed using the Cardinal code, which allows for a realistic representation of the temperature distribution. This report demonstrates this coupling approach on two test cases: a simplified proof-of-concept scenario and a realistic small-break loss-of-coolant accident scenario. In both cases, the effects of the spatially varying temperature are evident; although, because the small-break loss-of-coolant accident scenario was far from challenging the RPV’s margins, it is still inconclusive how important cold-plume effects are in realistic accident scenarios. However, this effort represents an important step toward performing realistic analyses to determine to what extent spatially nonuniform cooling might affect RPV integrity under more aggressive accident scenarios.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2002

Calculation and Correlation of the Unsteady Flowfield in a High Pressure Turbine

Forced vibrations in turbomachinery components can cause blades to crack or fail due to high-cycle fatigue. Such forced response problems will become more pronounced in newer engines with higher pressure ratios and smaller axial gap between blade rows. An accurate numerical prediction of the unsteady aerodynamics phenomena that cause resonant forced vibrations is increasingly important to designers. Validation of the computational fluid dynamics (CFD) codes used to model the unsteady aerodynamic excitations is necessary before these codes can be used with confidence. Recently published benchmark data, including unsteady pressures and vibratory strains, for a high-pressure turbine stage makes such code validation possible. In the present work, a three dimensional, unsteady, multi blade-row, Reynolds-Averaged Navier Stokes code is applied to a turbine stage that was recently tested in a short duration test facility. Two configurations with three operating conditions corresponding to modes 2, 3, and 4 crossings on the Campbell diagram are analyzed. Unsteady pressures on the rotor surface are compared with data.

Bakhle, Milind A.↗

NASA NTRS2008

Scalability of Localized Arc Filament Plasma Actuators

Temporal flow control of a jet has been widely studied in the past to enhance jet mixing or reduce jet noise. Most of this research, however, has been done using small diameter low Reynolds number jets that often have little resemblance to the much larger jets common in real world applications because the flow actuators available lacked either the power or bandwidth to sufficiently impact these larger higher energy jets. The Localized Arc Filament Plasma Actuators (LAFPA), developed at the Ohio State University (OSU), have demonstrated the ability to impact a small high speed jet in experiments conducted at OSU and the power to perturb a larger high Reynolds number jet in experiments conducted at the NASA Glenn Research Center. However, the response measured in the large-scale experiments was significantly reduced for the same number of actuators compared to the jet response found in the small-scale experiments. A computational study has been initiated to simulate the LAFPA system with additional actuators on a large-scale jet to determine the number of actuators required to achieve the same desired response for a given jet diameter. Central to this computational study is a model for the LAFPA that both accurately represents the physics of the actuator and can be implemented into a computational fluid dynamics solver. One possible model, based on pressure waves created by the rapid localized heating that occurs at the actuator, is investigated using simplified axisymmetric simulations. The results of these simulations will be used to determine the validity of the model before more realistic and time consuming three-dimensional simulations are conducted to ultimately determine the scalability of the LAFPA system.

Brown, Clifford A.↗

NASA NTRS2009

Data Parallel Line Relaxation (DPLR) Code User Manual: Acadia - Version 4.01.1

Data-Parallel Line Relaxation (DPLR) code is a computational fluid dynamic (CFD) solver that was developed at NASA Ames Research Center to help mission support teams generate high-value predictive solutions for hypersonic flow field problems. The DPLR Code Package is an MPI-based, parallel, full three-dimensional Navier-Stokes CFD solver with generalized models for finite-rate reaction kinetics, thermal and chemical non-equilibrium, accurate high-temperature transport coefficients, and ionized flow physics incorporated into the code. DPLR also includes a large selection of generalized realistic surface boundary conditions and links to enable loose coupling with external thermal protection system (TPS) material response and shock layer radiation codes.

Wright, Michael J.↗

DOE OSTI2014

Oxygen transport properties estimation by DSMC-CT simulations

Coupling DSMC simulations with classical trajectories calculations is emerging as a powerful tool to improve predictive capabilities of computational rarefied gas dynamics. The considerable increase of computational effort outlined in the early application of the method (Koura,1997) can be compensated by running simulations on massively parallel computers. In particular, GPU acceleration has been found quite effective in reducing computing time (Ferrigni,2012; Norman et al.,2013) of DSMC-CT simulations. The aim of the present work is to study rarefied Oxygen flows by modeling binary collisions through an accurate potential energy surface, obtained by molecular beams scattering (Aquilanti, et al.,1999). The accuracy of the method is assessed by calculating molecular Oxygen shear viscosity and heat conductivity following three different DSMC-CT simulation methods. In the first one, transport properties are obtained from DSMC-CT simulations of spontaneous fluctuation of an equilibrium state (Bruno et al, Phys. Fluids, 23, 093104, 2011). In the second method, the collision trajectory calculation is incorporated in a Monte Carlo integration procedure to evaluate the Taxman’s expressions for the transport properties of polyatomic gases (Taxman,1959). In the third, non-equilibrium zero and one-dimensional rarefied gas dynamic simulations are adopted and the transport properties are computed from the non-equilibrium fluxes of momentum and energy. The three methods provide close values of the transport properties, their estimated statistical error not exceeding 3%. The experimental values are slightly underestimated, the percentage deviation being, again, few percent.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2003

Numerical Models of Human Circulatory System under Altered Gravity: Brain Circulation

A computational fluid dynamics (CFD) approach is presented to model the blood flow through the human circulatory system under altered gravity conditions. Models required for CFD simulation relevant to major hemodynamic issues are introduced such as non-Newtonian flow models governed by red blood cells, a model for arterial wall motion due to fluid-wall interactions, a vascular bed model for outflow boundary conditions, and a model for auto-regulation mechanism. The three-dimensional unsteady incompressible Navier-Stokes equations coupled with these models are solved iteratively using the pseudocompressibility method and dual time stepping. Moving wall boundary conditions from the first-order fluid-wall interaction model are used to study the influence of arterial wall distensibility on flow patterns and wall shear stresses during the heart pulse. A vascular bed modeling utilizing the analogy with electric circuits is coupled with an auto-regulation algorithm for multiple outflow boundaries. For the treatment of complex geometry, a chimera overset grid technique is adopted to obtain connectivity between arterial branches. For code validation, computed results are compared with experimental data for steady and unsteady non-Newtonian flows. Good agreement is obtained for both cases. In sin-type Gravity Benchmark Problems, gravity source terms are added to the Navier-Stokes equations to study the effect of gravitational variation on the human circulatory system. This computational approach is then applied to localized blood flows through a realistic carotid bifurcation and two Circle of Willis models, one using an idealized geometry and the other model using an anatomical data set. A three- dimensional anatomical Circle of Willis configuration is reconstructed from human-specific magnetic resonance images using an image segmentation method. The blood flow through these Circle of Willis models is simulated to provide means for studying gravitational effects on the brain circulation under auto-regulation.

Kim, Chang Sung↗

NASA NTRS1996

Numerical Simulation of Flow in a Whirling Annular Seal and Comparison With Experiments

The turbulent flow field inside a whirling annular seal was simulated by using SCISEAL, a three-dimensional computational fluid dynamics code. The rotor center described a circular synchronous whirl. A rotating frame transformation was used to make the problem quasi-steady. The flow field at an axial Reynolds number of 24,000 and a Taylor number of 6600 was simulated. The standard kappa-epsilon model with wall functions and the low-Reynolds-number model were used to treat turbulence. An experimentally measured velocity field was used at the inlet boundary. Numerical predictions of the velocities and the stator wall pressures compared well with experimental data. Both turbulence models yielded nearly the same results. The capability of the SCISEAL code to analyze this complex flow field was demonstrated; the isotropic turbulence models performed adequately on this nonisotropic turbulence flow.

Athavale, M. M.↗

NASA NTRS2001

Fluid Aspects of Solar Wind Disturbances Driven by Coronal Mass Ejections

Transient disturbances in the solar wind initiated by coronal eruptions have been modeled for many years, beginning with the self-similar analytical models of Parker and Simon and Axford. The first numerical computer code (one-dimensional, gas dynamic) to study disturbance propagation in the solar wind was developed in the late 1960s, and a variety of other codes ranging from simple one-dimensional gas dynamic codes through three-dimensional gas dynamic and magnetohydrodynamic codes have been developed in subsequent years. For the most part, these codes have been applied to the problem of disturbances driven by fast CMEs propagating into a structureless solar wind. Pizzo provided an excellent summary of the level of understanding achieved from such simulation studies through about 1984, and other reviews have subsequently become available. More recently, some attention has been focused on disturbances generated by slow CMEs, on disturbances driven by CMEs having high internal pressures, and disturbance propagation effects associated with a structured ambient solar wind. Our purpose here is to provide a brief tutorial on fluid aspects of solar wind disturbances derived from numerical gas dynamic simulations. For the most part we illustrate disturbance evolution by propagating idealized perturbations, mimicking different types of CMEs, into a structureless solar wind using a simple one-dimensional, adiabatic (except at shocks), gas dynamic code. The simulations begin outside the critical point where the solar wind becomes supersonic and thus do not address questions of how the CMEs themselves are initiated. Limited to one dimension (the radial direction), the simulation code predicts too strong an interaction between newly ejected solar material and the ambient wind because it neglects azimuthal and meridional motions of the plasma that help relieve pressure stresses. Moreover, the code ignores magnetic forces and thus also underestimates the speed with which pressure disturbances propagate in the wind.

Gosling, J. T.↗

DOE OSTI2024

Computational study on the impact of gasoline-ethanol blending on autoignition and soot/NO x emissions under low-load gasoline compression ignition conditions

Here, in the present work, computational fluid dynamics (CFD) simulations of a single-cylinder gasoline compression ignition (GCI) engine are performed to investigate the impact of gasoline-ethanol blending on autoignition, nitrogen oxide (NO x ), and soot emissions under low-load conditions. In order to represent the test gasoline (RD5-87), a four-component toluene primary reference fuel (TPRF)+ethanol (ETPRF) surrogate (with 10% ethanol by volume; E10) is employed. A three-dimensional (3D) engine CFD model employing finite-rate chemistry with a skeletal kinetic mechanism (including NO x sub-mechanism), adaptive mesh refinement (AMR), and hybrid method of moments (HMOM) is adopted to capture the in-cylinder combustion phenomena and soot/NO x emissions. The engine CFD model is validated against experimental data for three gasoline-ethanol blends: E10, E30 and E100, with varying ethanol content by volume. Model validation is carried out for a broad range of start-of-injection (SOI) timings (−21, −27, −36, and −45 crank angle degrees (°CA) after top-dead-center (aTDC)) with respect to in-cylinder pressure, heat release rate, combustion phasing, NO x and soot emissions. For relatively later injection timings (−21 and −27 °CA aTDC), E30 yields higher amount of soot than E10; while the trend reverses for early injection cases (−36 and −45 °CA aTDC ). On the other hand, E100 yields the lowest amount of soot among all fuels irrespective of SOI timing. Further, E10 shows a non-monotonic trend in soot emissions with SOI timing: SOI-36>SOI-45>SOI-21>SOI-27, while soot emissions from E30 exhibit monotonic decrease with advancing SOI timing. NO x emissions from various fuels follow a trend of E10>E30>E100. On the other hand, NO x emissions increase as SOI timing is advanced for all fuels, with an anomaly for E10 and E100 where NO x decreases when SOI is advanced beyond −36 °CA aTDC. Detailed analysis of the numerical results is performed to investigate the soot/NO x emission trends and elucidate the impact of chemical composition and physical properties on autoignition and emissions characteristics.

Computational fluid dynamics↗

DOE OSTI2022

An arbitrary Lagrange-Eulerian investigation of HRAM shallow jet pre-spurt formation and time sensitivities to impact plate dynamics

With dry bay fires persisting as a significant contributor to aircraft vulnerability despite chronicled developments in survivability technologies, an accurate fire prediction capability remains paramount for credible vulnerability assessments. Physics based modeling of the hydrodynamic ram (HRAM) fluid deposition process is a key component of such capability, wherein capturing the first instance of fluid spurt, referred to herein as shallow jet spurts, is a core focus. Such pre-spurts as they have been formerly identified have only been witnessed sporadically in HRAM spurt experiments. In this work, ALE3D, a first-principles multi-physics code was employed to model the shallow jet spurt phenomenon with spherical projectiles impacting water-filled tanks faced with aluminum panels, such that the underlying physics and sensitivities could be explored. Development and verification of the 2D axisymmetric model is described relative to trends observed in a prior experimental campaign. Results from the verified model suggest shallow jet spurts have at least a quadratic sensitivity to the fundamental vibrational mode of the impact plate across impact velocities 610–1829 m/s (2000–6000 ft/s). It is further explained how shallow jet spurts are arrested for impacts at the two extremes of plate rigidity. This research constitutes the first fluid-structure modeling of shallow jet spurts which future three-dimensional analyses will expound upon.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1991

An evaluation of three two-dimensional computational fluid dynamics codes including low Reynolds numbers and transonic Mach numbers

Full-potential, Euler, and Navier-Stokes computational fluid dynamics (CFD) codes were evaluated for use in analyzing the flow field about airfoils sections operating at Mach numbers from 0.20 to 0.60 and Reynolds numbers from 500,000 to 2,000,000. The potential code (LBAUER) includes weakly coupled integral boundary layer equations for laminar and turbulent flow with simple transition and separation models. The Navier-Stokes code (ARC2D) uses the thin-layer formulation of the Reynolds-averaged equations with an algebraic turbulence model. The Euler code (ISES) includes strongly coupled integral boundary layer equations and advanced transition and separation calculations with the capability to model laminar separation bubbles and limited zones of turbulent separation. The best experiment/CFD correlation was obtained with the Euler code because its boundary layer equations model the physics of the flow better than the other two codes. An unusual reversal of boundary layer separation with increasing angle of attack, following initial shock formation on the upper surface of the airfoil, was found in the experiment data. This phenomenon was not predicted by the CFD codes evaluated.

Hicks, Raymond M.↗

NASA NTRS1992

Numerical procedures for three-dimensional computational surface thermochemistry

Models and equations for surface thermochemistry and near-surface thermophysics of aerodynamically-heated thermal protection materials are reviewed, with particular emphasis on computational boundary conditions for surface mass and energy transfer. The surface energy and mass balances, coupled with an appropriate ablation or surface catalysis model, provide complete thermochemical boundary conditions for a true multidisciplinary solution of the fully coupled fluid-dynamics/solid mechanics problem. Practical approximate solutions can be obtained by using a detailed model with full thermophysics for either the solid or fluid phase amd a semianalytic method for the other half of the problem. A significant increase in the state-of-the-art in aerothermal computational fluid dynamics is possible by uniting CFD methodology with surface thermochemistry boundary conditions and the heat-balance-integral method.

Milos, Frank S.↗

DOE OSTI2013

Flow and axial dispersion in a sinusoidal-walled tube: Effects of inertial and unsteady flows

Dispersion in porous media flows has been the subject of much experimental, theoretical and numerical study. Here we consider a wavy-walled tube (a three-dimensional tube with sinusoidally-varying diameter) as a simplified conceptualization of flow in porous media, where constrictions represent pore throats and expansions pore bodies. A theoretical model for effective (macroscopic) longitudinal dispersion in this system has been developed by volume averaging the microscale velocity field. Direct numerical simulation using computational fluid dynamics (CFD) methods was used to compute velocity fields by solving the Navier-Stokes equations, and also to numerically solve the volume averaging closure problem, for a range of Reynolds numbers (Re) spanning the low-Re to inertial flow regimes, including one simulation at Re = 449 for which unsteady flow was observed. Dispersion values were computed using both the volume averaging solution and a random walk particle tracking method, and results of the two methods were shown to be consistent. Our results are compared to experimental measurements of dispersion in porous media and to previous theoretical results for the low-Re, Stokes flow regime. In the steady inertial regime we observe an power-law increase in effective longitudinal dispersion (DL) with Re, consistent with previous results. This rapid rate of increase is caused by trapping of solute in expansions due to flow separation (eddies). For the unsteady case (Re = 449), the rate of increase of DL with Re was smaller than that observed at lower Re. Velocity fluctuations in this regime lead to increased rates of solute mass transfer between the core flow and separated flow regions, thus diminishing the amount of tailing caused by solute trapping in eddies and thereby reducing longitudinal dispersion.

CFD↗

NASA NTRS1993

Development of a computational testbed for numerical simulation of combustion instability

A synergistic hierarchy of analytical and computational fluid dynamic techniques is used to analyze three-dimensional combustion instabilities in liquid rocket engines. A mixed finite difference/spectral procedure is employed to study the effects of a distributed vaporization zone on standing and spinning instability modes within the chamber. Droplet atomization and vaporization are treated by a variety of classical models found in the literature. A multi-zone, linearized analytical solution is used to validate the accuracy of the numerical simulations at small amplitudes for a distributed vaporization region. This comparison indicates excellent amplitude and phase agreement under both stable and unstable operating conditions when amplitudes are small and proper grid resolution is used. As amplitudes get larger, expected nonlinearities are observed. The effect of liquid droplet temperature fluctuations was found to be of critical importance in driving the instabilities of the combustion chamber.

Grenda, Jeffrey↗

NASA NTRS1999

Design of a Low Speed Fan Stage for Noise Suppression

This report describes the design of a low tip speed, moderate pressure rise fan stage for demonstration of noise reduction concepts. The fan rotor is a fixed-pitch configuration delivering a design pressure ratio of 1.378 at a specific flow of 43.1 lbm/sec/sq ft. Four exit stator configurations were provided to demonstrate the effectiveness of circumferential and axial sweep in reducing rotor-stator interaction tone noise. The fan stage design was combined with an axisymmetric inlet, conical convergent nozzle, and nacelle to form a powered fan-nacelle subscale model. This model has a 22-inch cylindrical flow path and employs a rotor with a 0.30 hub-to-tip radius ratio. The design is fully compatible with an existing NASA force balance and rig drive system. The stage aerodynamic and structural design is described in detail. Three-dimensional (3-D) computational fluid dynamics (CFD) tools were used to define optimum airfoil sections for both the rotor and stators. A fan noise predictive system developed by Pratt & Whitney under contract to NASA was used to determine the acoustic characteristics of the various stator configurations. Parameters varied included rotor-to-stator spacing and vane leading edge sweep. The structural analysis of the rotor and stator are described herein. An integral blade and disk configuration was selected for the rotor. Analysis confirmed adequate low cycle fatigue life, vibratory endurance strength, and aeroelastic suitability. A unique load carrying stator arrangement was selected to minimize generation of tonal noise due to sources other than rotor-stator interaction. Analysis of all static structural components demonstrated adequate strength, fatigue life, and vibratory characteristics.

Dalton, W. N.↗

DOE OSTI2011

Progress on a New Integrated 3-D UCG Simulator and its Initial Application

A comprehensive simulator is being developed for underground coal gasification (UCG), with the capability to support site selection, design, hazard analyses, operations, and monitoring (Nitao et al., 2010). UCG is computationally challenging because it involves tightly-coupled multi-physical/chemical processes, with vastly different timescales. This new capability will predict cavity growth, product gas composition and rate, and the interaction with the host environment, accounting for site characteristics, injection gas composition and rate, and associated water-well extraction rates. Progress on the new simulator includes completion and system integration of a wall model, a rock spalling model, a cavity boundary tracking model, a one-dimensional cavity gas reactive transport model, a rudimentary rubble heat, mass, and reaction model, and coupling with a pre-existing hydrology simulator. An existing geomechanical simulator was enhanced to model cavity collapse and overburden subsidence. A commercial computational fluid dynamics (CFD) code is being evaluated to model cavity gas flow and combustion in two and three dimensions. Although the simulator is midway in its development, it was applied to modeling the Hoe Creek III field test (Stephens, 1981) conducted in the 1970s, in order to evaluate and demonstrate the simulator's basic capabilities, gain experience, and guide future development. Furthermore, it is consistent with our philosophy of incremental, spiral software development, which helps in identifying and resolving potential problems early in the process. The simulation accounts for two coal seams, two injection points, and air and oxygen phases. Approximate extent and shape of cavity growth showed reasonable agreement with interpreted field data. Product gas composition and carbon consumed could not be simultaneously matched for a given set of parameter values due to the rudimentary rubble model currently used, although they can be matched using separate parameter sets. This result is not surprising and confirms plans for a more sophisticated rubble model as our next step, as well as adding geomechanical collapse modeling and higher accuracy cavity gas reactive transport models. The results are very encouraging and demonstrate that our approach is sound.

01 COAL, LIGNITE, AND PEAT↗