Search NASA⌕ Search

SEARCH · Search NASA

Results for “process simulation”

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,711 records · Page 95

Isotherm Modeling and Techno-Economic Analysis of Contactor Technologies for New Tetraamine-Appended MOF for NGCC Applications

Recently, a family of tetraamine-functionalized metal-organic frameworks (MOF) has been reported as promising sorbent materials for capturing CO2 from flue gas conditions relevant to natural gas combined cycle (NGCC) applications. The main advantages of these materials are their two-step cooperative CO2 adsorption, which gives rise to unusual two step-shaped CO2 adsorption profiles and their high thermal stability. This work presents the modelling of the two-transition isotherm of the tetraamine-appended MOF, N,N'-bis(3-aminopropyl)-1,4-diaminobutane (3-4-3)-appended Mg2(dobpdc), and the techno-economic analysis (TEA) of carbon capture processes utilizing this sorbent. Due to the unusual isotherm shapes of the experimental CO2 adsorption data for tetraamine-appended Mg2(dobpdc) and the strong nonlinearity of CO2 loading with respect to temperature and pressure, we tested two different models which use logistic functions for representing the different isotherm behaviors in the different pressure ranges. The first model uses the quadratic isotherm model in the low-pressure region, the Langmuir isotherm model in the middle pressure range, and the dual site Langmuir isotherm model in the high-pressure range. To model the transition between regions we used the arctangent functions independent of temperature and the thermal effect was accounted by using the Clausius-Clapeyron relation]. The second model is an extension of the weighted dual-site Langmuir isotherm model presented by Hughes et al.. In this extended model, the dual-site Langmuir isotherm is employed in the three transition regions, using temperature-dependent logistic functions to activate or deactivate the isotherm model in the low, middle, and high-pressure ranges. Both models fit the experimental data quite well with root mean squared errors (RMSE’s) of 0.41 and 0.17 for model 1 and model 2, respectively. Since model 2 resulted in a lower RMSE, it was leveraged for the development of the gas/solid contactor models used by the TEA. Specifically, two different contactor models, an axial-flow fixed bed and moving bed contactor, were developed as part of this work. These models are dynamic, pressure-driven, and consist of mass, energy, and momentum conservation equations. A kinetic model was also developed by performing parameter estimation using experimental fixed bed breakthrough data. These models are then used to simulate CO2 capture processes from the flue gas generated from a ~600 gross MW NGCC power plant. A cost model was developed which considers the capital cost of the reactors and the significant operating costs such as steam and electricity. Using NETL’s Framework for the Optimization and Quantification of Uncertainty of Uncertainty and Surrogates tool (FOQUS), which has the capability of linking models built using numerous modelling platforms with derivative-free optimization solvers, a techno-economic optimization of the carbon capture processes was performed which minimizes the cost of capture.

Caballero, Daison↗

FULL RANGE TUNE SCAN STUDIES USING GRAPHICS PROCESSING UNITS WITH CUDA IN EIC BEAM-BEAM SIMULATIONS

The hadron beam in the Electron-Ion Collider (EIC) suffers high order betatron and synchro-betatron resonances. In this paper, we present a weak-strong full range (0.0 ~ 0.5) fractional tune scan with a step size as small as 0.001. Multiple Graphics Processing Units (GPUs) are used to speed up the simulation. A code parallelized with MPI and CUDA is implemented. The good tune region from weak-strong scan is further checked by the self-consistent strong-strong simulation. This study provides beam dynamics guidance in choosing proper working points for the future EIC.

43 PARTICLE ACCELERATORS↗

Pseudo-dynamic Aspen Plus{sup R} Model of the Defense Waste Processing Facility - 18475

The Defense Waste Processing Facility (DWPF) at the Savannah River Site converts legacy high-level waste from the Cold War into borosilicate glass for permanent disposal in a geological repository. The waste sludge is chemically adjusted by the addition of acid, mixed with frit, and concentrated in the DWPF Chemical Process Cell (CPC) to ensure that the final glass has the desired properties. The resulting mixture is fed into the melter where the waste is melted and poured into stainless steel storage canisters. We present an Aspen Plus{sup R} model of the DWPF CPC, which is being developed to aid chemical processing decisions pertaining to the current flowsheet and to aid decisions regarding optimizing future processing. Given a waste feed composition with acid and frit additions, the model will predict the sludge and off-gas compositions throughout the processing, as well as the final melter feed composition. This allows different processing scenarios to be investigated easily, reducing the reliance on experimental testing, and hence saving time and money. The model has multiple potential applications, including calculating the acid and frit additions required to give acceptable glass properties, off-gas permitting, tracking mercury steam-stripped from the waste in the off-gas system, tracking flammable gases in the offgas system, investigating the effect of reducing air purges in the condensers, and general process optimization. The waste is processed in the CPC in batches, requiring dynamic simulation to track the time-dependent operations. Accurate modeling of the processing also requires calculations of phase and electrolyte chemical equilibria and parallel rate-dependent reactions. Due to existing technical limitations with dynamic simulation of the complex chemical system, we have adopted a pseudo-dynamic approach in which time intervals are treated as consecutive steady-state operations in Aspen Plus{sup R}. Thus, the model consists of ∼100 replicas of the processing vessel and its off-gas system, with the sludge output by one replica (representing one time interval) becoming the feed to the next (representing the next time interval). By examining the results from these replicas in sequence, one can see the time evolution of the model. This general approach can be used wherever simultaneous calculation of phase and electrolyte chemical equilibria in a dynamically evolving system is needed. In addition to describing the general structure and chemistry of the model, we present results from tests of the model against experimental data. These data are from laboratory-scale process simulations of the CPC with a nonradioactive waste sludge simulant. These tests include the model's predictions of the offgas composition, the sludge pH, and the final sludge composition. We will also discuss planned future developments of the model. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NEPTUNE ON TIPTOES: DYNAMICAL HISTORIES THAT PRESERVE THE COLD CLASSICAL KUIPER BELT

The current dynamical structure of the Kuiper Belt was shaped by the orbital evolution of the giant planets, especially Neptune, during the era following planet formation when the giant planets may have undergone planet-planet scattering and/or planetesimal-driven migration. Numerical simulations of this process, while reproducing many properties of the Belt, fail to generate the high inclinations and eccentricities observed for some objects while maintaining the observed dynamically 'cold' population. We present the first of a three-part parameter study of how different dynamical histories of Neptune sculpt the planetesimal disk. Here we identify which dynamical histories allow an in situ planetesimal disk to remain dynamically cold, becoming today's cold Kuiper Belt population. We find that if Neptune undergoes a period of elevated eccentricity and/or inclination, it secularly excites the eccentricities and inclinations of the planetesimal disk. We demonstrate that there are several well-defined regimes for this secular excitation, depending on the relative timescales of Neptune's migration, the damping of Neptune's orbital inclination and/or eccentricity, and the secular evolution of the planetesimals. We model this secular excitation analytically in each regime, allowing for a thorough exploration of parameter space. Neptune's eccentricity and inclination can remain high for a limited amount of time without disrupting the cold classical belt. In the regime of slow damping and slow migration, if Neptune is located (for example) at 20 AU, then its eccentricity must stay below 0.18 and its inclination below 6 Degree-Sign .

79 ASTRONOMY AND ASTROPHYSICS↗

SPECTRAL INDEX AS A FUNCTION OF MASS ACCRETION RATE IN BLACK HOLE SOURCES: MONTE CARLO SIMULATIONS AND AN ANALYTICAL DESCRIPTION

We present herein a theoretical study of correlations between spectral indexes of X-ray emergent spectra and mass accretion rate ( m-dot ) in black hole (BH) sources, which provide a definitive signature for BHs. It has been firmly established, using the Rossi X-ray Timing Explorer (RXTE) in numerous BH observations during hard-soft state spectral evolution, that the photon index of X-ray spectra increases when m-dot increases and, moreover, the index saturates at high values of m-dot . In this paper, we present theoretical arguments that the observationally established index saturation effect versus mass accretion rate is a signature of the bulk (converging) flow onto the BH. Also, we demonstrate that the index saturation value depends on the plasma temperature of converging flow. We self-consistently calculate the Compton cloud (CC) plasma temperature as a function of mass accretion rate using the energy balance between energy dissipation and Compton cooling. We explain the observable phenomenon, index- m-dot correlations using a Monte Carlo simulation of radiative processes in the innermost part (CC) of a BH source and we account for the Comptonization processes in the presence of thermal and bulk motions, as basic types of plasma motion. We show that, when m-dot increases, BH sources evolve to high and very soft states (HSS and VSS, respectively), in which the strong blackbody(BB)-like and steep power-law components are formed in the resulting X-ray spectrum. The simultaneous detections of these two components strongly depends on sensitivity of high-energy instruments, given that the relative contribution of the hard power-law tail in the resulting VSS spectrum can be very low, which is why, to date RXTE observations of the VSS X-ray spectrum have been characterized by the presence of the strong BB-like component only. We also predict specific patterns for high-energy e-fold (cutoff) energy (E{sub fold}) evolution with m-dot for thermal and dynamical (bulk) Comptonization cases. For the former case, E{sub fold} monotonically decreases with m-dot , in the latter case, the E{sub fold} decrease is followed by its increase at high values of m-dot . The observational evolution of E{sub fold} versus m-dot can be another test for the presence of a converging flow effect in the formation of the resulting spectra in the close vicinity of BHs.

79 ASTRONOMY AND ASTROPHYSICS↗

THE SDSS-III BARYON OSCILLATION SPECTROSCOPIC SURVEY: THE QUASAR LUMINOSITY FUNCTION FROM DATA RELEASE NINE

We present a new measurement of the optical quasar luminosity function (QLF), using data from the Sloan Digital Sky Survey-III: Baryon Oscillation Spectroscopic Survey (SDSS-III: BOSS). From the SDSS-III Data Release Nine, a uniform sample of 22,301 i {approx}< 21.8 quasars are selected over an area of 2236 deg{sup 2}, with confirmed spectroscopic redshifts between 2.2 < z < 3.5, filling in a key part of the luminosity-redshift plane for optical quasar studies. The completeness of the survey is derived through simulated quasar photometry, and this completeness estimate is checked using a sample of quasars selected by their photometric variability within the BOSS footprint. We investigate the level of systematics associated with our quasar sample using the simulations, in the process generating color-redshift relations and a new quasar K-correction. We probe the faint end of the QLF to M{sub i} (z = 2.2) Almost-Equal-To -24.5 and see a clear break in the QLF at all redshifts up to z = 3.5. A log-linear relation (in log {Phi}* - M*) for a luminosity evolution and density evolution model is found to adequately describe our data within the range 2.2 < z < 3.5; across this interval the break luminosity increases by a factor of {approx}2.6 while {Phi}* declines by a factor of {approx}8. At z {approx}< 2.2 our data are reasonably well fit by a pure luminosity evolution model, and only a weak signature of ''AGN downsizing'' is seen, in line with recent studies of the hard X-ray luminosity function. We compare our measured QLF to a number of theoretical models and find that models making a variety of assumptions about quasar triggering and halo occupation can fit our data over a wide range of redshifts and luminosities.

79 ASTRONOMY AND ASTROPHYSICS↗

THE CENTRAL DARK MATTER DISTRIBUTION OF NGC 2976

We study the mass distribution in the late-type dwarf galaxy NGC 2976 through stellar kinematics obtained with the Visible Integral Field Replicable Unit Spectrograph Prototype and anisotropic Jeans models as a test of cosmological simulations and baryonic processes that putatively alter small-scale structure. Previous measurements of the H{alpha} emission-line kinematics have determined that the dark matter halo of NGC 2976 is most consistent with a cored density profile. We find that the stellar kinematics are best fit with a cuspy halo. Cored dark matter halo fits are only consistent with the stellar kinematics if the stellar mass-to-light ratio is significantly larger than that derived from stellar population synthesis, while the best-fitting cuspy model has no such conflict. The inferred mass distribution from a harmonic decomposition of the gaseous kinematics is inconsistent with that of the stellar kinematics. This difference is likely due to the gas disk not meeting the assumptions that underlie the analysis such as no pressure support, a constant kinematic axis, and planar orbits. By relaxing some of these assumptions, in particular the form of the kinematic axis with radius, the gas-derived solution can be made consistent with the stellar kinematic models. A strong kinematic twist in the gas of NGC 2976's center suggests caution, and we advance the mass model based on the stellar kinematics as more reliable. The analysis of this first galaxy shows promising evidence that dark matter halos in late-type dwarfs may in fact be more consistent with cuspy dark matter distributions than earlier work has claimed.

79 ASTRONOMY AND ASTROPHYSICS↗

Four-Dimensional Dose Reconstruction for Scanned Proton Therapy Using Liver 4DCT-MRI

Purpose: Four-dimensional computed tomography-magnetic resonance imaging (4DCT-MRI) is an image-processing technique for simulating many 4DCT data sets from a static reference CT and motions extracted from 4DMRI studies performed using either volunteers or patients. In this work, different motion extraction approaches were tested using 6 liver cases, and a detailed comparison between 4DCT-MRI and 4DCT was performed. Methods and Materials: 4DCT-MRI has been generated using 2 approaches. The first approach used motion extracted from 4DMRI as being “most similar” to that of 4DCT from the same patient (subject-specific), and the second approach used the most similar motion obtained from a motion library derived from 4DMRI liver studies of 13 healthy volunteers (population-based). The resulting 4DCT-MRI and 4DCTs were compared using scanned proton 4D dose calculations (4DDC). Results: Dosimetric analysis showed that 93% ± 8% of points inside the clinical target volume (CTV) agreed between 4DCT and subject-specific 4DCT-MRI (gamma analysis: 3%/3 mm). The population-based approach however showed lower dosimetric agreement with only 79% ± 14% points in the CTV reaching the 3%/3 mm criteria. Conclusions: 4D CT-MRI extends the capabilities of motion modeling for dose calculations by accounting for realistic and variable motion patterns, which can be directly employed in clinical research studies. We have found that the subject-specific liver modeling appears more accurate than the population-based approach. The former is particularly interesting for clinical applications, such as improved target delineation and 4D dose reconstruction for patient-specific QA to allow for inter- and/or intra-fractional plan corrections.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Reduction of CO 2 Emissions Through Lightweight Body Panels (Project Final Report)

Lightweight construction is an integral part of Volkswagen’s overall strategy of reducing CO 2 emissions. Due to its low cost, steel is the most commonly used material for automotive exterior body panels today. Unfortunately, steel has a high density, resulting in a relatively low specific strength. Glass fiber based sheet molding compounds (SMC) provide high properties in combination with lower density. The high specific strength of SMC offers an enormous lightweight potential. To unlock the full potential of SMC materials in combination with cutting edge manufacturing processes, Volkswagen Group of America worked together with IACMI (Institute for Advanced Composites Manufacturing Innovation) and the academic partners: University of Tennessee Knoxville, Purdue University and Michigan State University; as well as the industry partners: Ashland, IDI, Owens Corning and Continental Structural Plastics. Leveraging the expertise of all project partners, reaching over the entire supply chain the project demonstrated the potential of these materials. This report will highlight the major steps in the development process on the way to technology readiness for SMC using the example of the Volkswagen Atlas Liftgate. Over the 36-month period of the project, the work focused on three R&D areas: material development, design and simulation, and development of the manufacturing process. Material selection included alternative fibers and resin systems, and accounted for material availability, properties, and cost-efficiency. The work undertaken in the field of design and simulation has pushed the envelope of short fiber reinforced thermoset molding compound process simulation. Design iterations were virtually tested, while the final design was used to validate the simulation software against physical parts. Manufacturing development used cutting-edge technology, while experts along the supply chain were working together to ensure the best possible results. In the final stage of the project, liftgates were molded, trimmed, bonded, painted and assembled before exhaustive testing. The result is an e-coat (electrophoretic dip coating) capable Class-A SMC liftgate, which is ready for high-volume production, and can be used as a technology demonstrator. The prototypes manufactured in the scope of this work have exhibited a mass reduction for the Volkswagen Atlas liftgate of up to 35% compared to the series production steel version, without a degradation of the functionalities.

36 MATERIALS SCIENCE↗

Metal additive manufacturing simulation across length, time, and computing scales

Metal additive manufacturing (AM) offers a unique opportunity for production of advanced materials and complex geometries. However, variability in microstructure and properties challenges conventional approaches to design, process optimization, qualification, and materials selection. Modeling and simulation can improve understanding of AM processing and materials, but also poses major challenges for existing computational methods. Simultaneously, modern scientific computing hardware has become increasingly complex, most notably with the adoption of hybrid architectures such as Graphical Processing Units (GPUs). If appropriately utilized, emerging computational capabilities provide an opportunity to reveal new insight into AM processing and the resulting material structure and properties. In this review we describe the computational AM landscape, identify critical gaps, and highlight opportunities to impact the development and application of AM. First, the requirements and challenges of representative AM problem statements will be defined. Here, these problems range from scientific studies to industrial applications and are designed to capture the breadth of challenges facing the AM community. Next, the current state of AM modeling and simulation is evaluated, broken down by enabling hardware and software, process simulation, microstructure simulation, and property simulation. Each section describes the diversity of simulation approaches and associated trade-offs in physical fidelity and computational expense. Each area is then assessed based on their suitability and readiness for current and developing computational architectures. Lastly, the greatest opportunities for future research and application are highlighted, including gaps in modeling capabilities, opportunities for near-term application, and key scientific challenges.

additive manufacturing↗

The impact of ARM on climate modeling

Climate models are among humanity’s most ambitious and elaborate creations. They are designed to simulate the interactions of the atmosphere, ocean, land surface, and cryosphere on time scales far beyond the limits of deterministic predictability and including the effects of time-dependent external forcings. The processes involved include radiative transfer, fluid dynamics, microphysics, and some aspects of geochemistry, biology, and ecology. The models explicitly simulate processes on spatial scales ranging from the circumference of Earth down to 100 km or smaller and implicitly include the effects of processes on even smaller scales down to a micron or so. In addition, the atmospheric component of a climate model can be called an atmospheric global circulation model (AGCM).

58 GEOSCIENCES↗

Open‐Source Anaerobic Digestion Modeling Platform, Anaerobic Digestion Model No. 1 Fast (ADM1F)

An open‐source modeling platform, called Anaerobic Digestion Model No. 1 Fast (ADM1F), is introduced to achieve fast and numerically stable simulations of anaerobic digestion processes. ADM1F is compatible with an iPython interface to facilitate model configuration, simulation, data analysis, and visualization. Faster simulations and more stable results are accomplished by implementing an advanced open‐source library of numerical methods called Portable Extensive Toolkit for Scientific Computation (PETSc) to solve the ADM1 system of equations. Leveraging PETSc, ADM1F can consistently complete a steady‐state simulation under 0.2 s, over 99% faster than a benchmark ADM1 model implemented with MATLAB while achieving agreement of model outputs within 1% of those obtained with the benchmark model. For dynamic simulations, however, ADM1F has a computational speed advantage only when the influent characteristics update more frequently than every 4 h. The ability of ADM1F to be useful as a tool to study anaerobic digestion systems is demonstrated through two example implementations of ADM1F: (1) a two‐phase co‐digestion scenario evaluating the impact of the organic loading rate and the substrate composition on reactor performance and stability, and (2) a conventional digester scenario assessing the effectiveness of recovery strategies after disruptions that led to instability. These examples demonstrate how the high simulation speed and the convenience of the iPython interface allow ADM1F to complete complex analyses within minutes, much faster than computational strategies currently reported in the literature.

anaerobic co-digestion↗

Computer-aided Veress needle guidance using endoscopic optical coherence tomography and convolutional neural networks

During laparoscopic surgery, the Veress needle is commonly used in pneumoperitoneum establishment. Precise placement of the Veress needle is still a challenge for the surgeon. In this study, a computer-aided endoscopic optical coherence tomography (OCT) system was developed to effectively and safely guide Veress needle insertion. This endoscopic system was tested by imaging subcutaneous fat, muscle, abdominal space, and the small intestine from swine samples to simulate the surgical process, including the situation with small intestine injury. Each tissue layer was visualized in OCT images with unique features and subsequently used to develop a system for automatic localization of the Veress needle tip by identifying tissue layers (or spaces) and estimating the needle-to-tissue distance. We used convolutional neural networks (CNNs) in automatic tissue classification and distance estimation. In conclusion, the average testing accuracy in tissue classification was 98.53 ± 0.39%, and the average testing relative error in distance estimation reached 4.42 ± 0.56% (36.09 ± 4.92 μm).

59 BASIC BIOLOGICAL SCIENCES↗

A sparse Gaussian sigmoid basis function approximation of hyperspectral data for detection of solids

We define a new characterization of emissivity and reflectance curves for compositional exploitation of hyperspectral data. Here, our method decomposes each spectrum into a sparse set of Gaussian sigmoid components using penalized regression. Detection is based on the combination of Gaussian sigmoid components unique to a target material. Focusing on the presence of spectral upslopes and downslopes rather than spectral correlations makes detection more robust to both target variation and spectral variability from atmosphere and background encountered during the collection process. We present simulation studies that demonstrate the potential to reduce false positive rates without compromising sensitivity. Characterization of long–wave infrared (LWIR) experimental data validates our method using minerals of different particle sizes, measurement angles, and collection conditions.

97 MATHEMATICS AND COMPUTING↗

Strain Field Around Individual Dislocations Controls Failure

Understanding material failure on a fundamental level is a key aspect in the design of robust structural materials, especially for metals and alloys capable to undergo plastic deformation. In the last decade, significant progress is made in quantifying the stresses associated with failure in both experiments and simulations. Nonetheless, the processes occurring on the most essential level of individual dislocations that govern semi-brittle and ductile fracture are still experimentally not accessible, limiting the failure prediction capabilities. Therefore, in the present work, a one-of-a-kind nanoscale fracture experiment is conducted on a single crystalline Cr bending beam in situ in the transmission electron microscope and for the first time quantify the transient strains around individual dislocations, as well as of the whole dislocation network during crack opening. The results reveal the importance of both pre-existing and newly emitted dislocations for crack-tip shielding via their intrinsic strain field and provide guidelines to design more damage tolerant materials.

36 MATERIALS SCIENCE↗

Process-driven roadmap for depositing super duplex stainless steel via wire Arc additive manufacturing

Here, this study systematically investigates the effects of shielding gas, bead spacing, weld mode, and travel speed on the phase balance, porosity, and hardness of wire arc additively manufactured (WAAM) super duplex stainless steel (ER2594). Robotic WAAM was employed to fabricate multilayer walls under systematically varied process conditions, followed by phase transformation simulations, X-ray computed tomography (XCT), electron backscatter diffraction (EBSD), and microhardness evaluation. Thermodynamic simulations predicted rapid cooling of AM process can suppress the potential formation of deleterious precipitates which was later validated via cross-sectional microstructure analyses of printed samples. XCT revealed porosity levels below 0.2% for all deposits, with reduced travel speed significantly lowering defect volume. Microstructural analyses revealed the evolution of various austenite precipitates, including grain boundary austenite (GBA), Widmanstätten austenite (WA), and intergranular austenite (IGA), sequentially upon cooling of the ferrite phase. Among all process parameters, weld transfer mode exhibited the strongest influence on phase balance; pulsed mode promoted higher ferrite retention (~ 36%) compared to RapidX mode. No consistent relationship between stabilized phase fraction and captured microhardness was observed. This work provides critical insights for optimizing WAAM parameters to control phase balance and mechanical performance, which is essential for producing high-integrity super duplex stainless-steel components for nuclear and marine applications.

Grain orientation↗

Advances in Multimodal Characterization of Structural Materials

The myriad detectors and instruments now available for materials characterization provide researchers with an ever-growing suite of tools to probe material behavior. Progress in the development of instrumentation and workflows that enable the collection, and leverage the potential, of various data modalities have provided novel insights into material behavior. Using data across multiple length scales, or performing complementary analyses of in situ and ex situ data, can help reveal a more complete picture of dynamic processes or material structure. However, the accurate combination, or fusion, of these disparate data modalities presents new challenges. Differences in resolution, as well as the varying length scales at which physical phenomena are exploited to generate these data, necessitate novel approaches to accurately interpret and combine these data. Furthermore, the papers within this special topic focus on the collection and fusion of multimodal data to better understand structural materials. From new frameworks and workflows for data segmentation and analysis, process monitoring, enhancing simulations, or interrogating mechanical response, these papers reveal the potential benefits of utilizing multimodal data.

36 MATERIALS SCIENCE↗

Permeability variations within mining-induced fractured rock mass and its influence on groundwater inrush

This paper is concerned with the evaluation of permeability of fractured rock mass due to the cover stress re-establishment, which is a major factor in controlling water and gas flow rate induced by mining operations in fractured rock. The case study considered in this paper is based on the results of observations of groundwater inrush and a decrease in water inflow from the fractured roof strata due to mining advancing in the Taiping Coalmine, Shandong Province, China. A conceptual model of an effective porous media was used to assess the permeability distribution in the fractured zone induced by coal mining. The cover stress re-establishment in gob fractured rock mass was evaluated using an empirical formula based on the surface subsidence. A simplified conceptual model of the fractured zone was used to evaluate the deformation of fractured zone along with the evaluation of changes in the rock permeability above the gob due to the cover stress re-establishment. These data were then used to calculate the water inflow rate into the panel. Predicted water inflow rates have been found to be in good agreement with those from monitoring data. This study improved the understanding of the mechanisms of the post-mining cover stress re-establishment on permeability change of the overburden fracture rock strata. These results can then be applied for numerical simulations of the process of overburden failure and consequent groundwater inrush due to coal mining.

54 ENVIRONMENTAL SCIENCES↗