Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity Matrix”

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 109 records · Page 6

Analysis of semivolatile organics in liquid radioactive residue using sorptive stir bar and solvent back-extraction

Current methods for semivolatiles analysis in radioactive samples can produce large volumes of radioactive solvent residue. A method utilizing stir-bar sorptive extraction has been explored in this work for its applicability to radioactive waste samples. This low solvent analytical method may accelerate remediation, minimize hazardous solvent waste, and reduce exposure risk to workers. Organic compounds (polyaromatic hydrocarbons, chlorinated aromatics, and phenolic compounds) were chosen as surrogates for common Liquid Waste System (LWS) contaminants at Savannah River Site (Aiken, SC). Stir-bar extraction parameters (extraction time, matrix modification, and effective pH range) and solvent back extraction parameters (solvent type, volume, and extraction time) were optimized experimentally for the chosen compounds. Affinity of the stir-bar extraction polymer to radionuclides Cs-137 and Am-241 was observed to determine radionuclide concentration effects. The stir-bar method achieved mean recovery of 100 ± 0.7% (1σ), relative to 114 ± 7% using solvent extraction, while reducing weekly method hands-on time by 93.4% and solvent volume consumption by 99.3%. Sensitivity was improved by 378% in simulated tank waste and 278% in real-world LWS matrix, relative to solvent extraction. This work has produced a safe and optimized method for the low solvent analysis of organics in legacy radioactive tank waste by stir-bar sorptive extraction.

GC-MS↗

Response Surface Methodology As a New Approach for Finding Optimal MALDI Matrix Spraying Parameters for Mass Spectrometry Imaging

Automated spraying devices have become ubiquitous in laboratories employing matrixassisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI), in part because they permit control of a number of matrix application parameters that can easily be reproduced for intra- and interlaboratory studies. Determining the optimal parameters for MALDI matrix application, such as temperature, flow rate, spraying velocity, number of spraying cycles, and solvent composition for matrix application, is critical for obtaining high-quality MALDI-MSI data. However, there are no established approaches for optimizing these multiple parameters simultaneously. Instead optimization is performed iteratively (i.e., one parameter at a time), which is time consuming and can lead to overall non-optimal settings. In this report, we demonstrate the use a novel experimental design and the response surface methodology to optimize five parameters of MALDI matrix application using a robotic sprayer. Thirtytwo combinations of MALDI matrix spraying conditions were tested, which allowed us to elucidate relationships between each of the application parameters as determined by MALDI-MS (specifically, using a 15 Tesla Fourier transform ion cyclotron resonance mass spectrometer). As such, we were able to determine the optimal automated spraying parameters that minimized signal de-localization and enabled high MALDI sensitivity. We envision this optimization strategy can be utilized for other MALDI-MSI applications, tissue types, and matrix application approaches.

60 APPLIED LIFE SCIENCES↗

Sorption‐Induced Fiber Optic Plasmonic Gas Sensing via Small Grazing Angle of Incidence

Abstract Sensing technologies based on plasmonic nanomaterials are of interest for various chemical, biological, environmental, and medical applications. In this work, an incorporation strategy of colloidal plasmonic nanoparticles (pNPs) in microporous polymer for realizing distinct sorption‐induced plasmonic sensing is reported. This approach is demonstrated by introducing tin‐doped indium oxide pNPs into a polymer of intrinsic microporosity (PIM‐1). The composite film (pNPs‐polymer) provides distinct and tunable optical features on the fiber optic (FO) platform that can be used as a signal transducer for gas sensing (e.g., CO 2 ) under atmospheric conditions. The resulting pNPs‐polymer composite demonstrates high sensitivity response on FO in the evanescent field configuration, provided by the dramatic response of modes above the total‐internal‐reflection angle. Furthermore, by varying the pNPs content in the polymer matrix, the optical behavior of the pNPs‐polymer composite film can be tuned to affect the operational wavelength by over several hundred nanometers and the sensitivity of the sensor in the near‐infrared range. It is also shown that the pNPs‐polymer composite film exhibits remarkable stability over a period of more than 10 months by mitigating the physical aging issue of the polymer.

Kim, Ki‐Joong↗

A proposed method for addressing large unphysical uncertainties in mubar

Mubar uncertainties in Section MF 34 MT 2 of ENDF/B-VIII.0 data are sometimes too large. Physical limitations of the bounds of mubar limit the maximum uncertainty to be < 1.0 and usually $\ll$ 1.0. Two ad hoc methods are proposed to artificially constrain the randomly sampled values of mubar to allowable values. The application of NJOY’s mubar covariance matrix to the “sandwich rule” is also discussed and a consistent way to apply sensitivity coefficients is described. Mubar uncertainties for the Jezebel k eff were found to be 156 pcm – a value comparable to other evaluated cross section uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Development of NDE/NDT Tools for High-Volume & High-Speed Inspection of CFRP Structures in Automotive Manufacturing

Main advantages of the air-coupled ultrasound testing (ACUT) and electromagnetic testing (EMT) techniques for NDE of CFRP composites were non-contact sensing, scalability for high-speed inspection, cost-effectiveness, and non-hazardous operation. Despite these advantages, no systems that would satisfy the project requirements were commercially available. Hence, one of the major efforts of the Michigan State University (MSU) team at the initial stage of the project was to close this technological gap by developing, optimizing, and validating array sensors that would provide sufficient sensitivity, spatial coverage, and resolution for robust defect detection. Optimization of the ACUT and EMT sensor designs was performed using experimentally validated finite element models. Initial experiments using array probes were conducted on relatively flat CFRP samples. In parallel, the MSU team designed and assembled a portable platform with two robotic arms. The robots were equipped with newly designed sensors that enabled high-speed NDE of curved CFRP parts. Presently, the developed robotic platform can be used as a demo/template NDE system, which is easily adaptable to manufacturing environments and in-line NDE. The ACUT NDE system developed by the MSU team used a high-power 4-channel pulser receiver for parallel data acquisition. The array probes were designed by stacking commercially available ACUT transducers, which operated in the frequency range between 100 kHz and 500 kHz. MSU optimized the excitation procedure and developed wave focusing cones so as to reduce the crosstalk between the transducers and to provide higher pulse repletion frequency (PRF). The through-transmission (TT) and single-side access (SSA) inspection modes were successfully implemented. In the TT-ACUT, structural defects in CFRP were detected by passing ultrasonic waves through the test part. Hence, the ACUT transmitters and receivers needed to be placed on the opposite sides of the test part. In the SSA-ACUT, guided waves (GW) were excited in the test part using the transmitters and were sensed by the receivers from the same side. Multi-channel TT-ACUT and SSA-ACUT provided high-speed NDE, and were successfully validated on CFRP test samples with interlaminar delaminations and other embedded defects The EM techniques developed by the MSU team included: 1) eddy current testing (ECT), 2) capacitive imaging (CI) and hybrid dual-mode imaging. In ECT, structural damage was detected in CFRP using coils sensor arrays. In ECT, the excitation magnetic field is generated by passing an alternating current through a coil, which is placed above the test sample. The excitation field penetrates the conductive sample and induces the eddy currents in its transect. In turn, the eddy currents generate the reaction field, which affects the total field sensed by a coil. Hence, the presence of structural flaws will alter the eddy current flow and the picked-up signal. ECT is mostly sensitive to local changes of the electric conductivity of the test sample, and CFRPs are mostly conductive in the direction of carbon fibers. Hence, ECT was well suited for the detection of fiber damage/fiber irregularities. The MSU team developed printed circuit boards (PCB) with coil sensor arrays optimized for NDE of CFRP. Unlike most commercial probes designed for ECT of metallic structures, the MSU array probes were designed for operation in [1-10] MHz frequency range, which was optimal for low-conductive CFRP. Multiple sensing topologies (coil groups excitation/sensing arrangements) were implemented and successfully validated. Capacitive Imaging (CI) technique developed by MSU was complementary to ECT. In contrast to ECT, which was sensitive to local changes of the electrical conductivity, the CI was sensitive to local changes of the dielectric constant. Therefore, CI could provide information about matrix damage/matrix irregularities in CFRP. The MSU CI sensor arrays were made of multiple circular or rectangular open-plate capacitors printed on PCB. Sensors of this type are not commercially available. In addition to ECT and CI, the MSU team developed a hybrid (dual-mode) inductive/capacitive measurement technique that synergistically combined the benefits of inductive and capacitive sensing for rapid NDE of fiber reinforced polymer (FRP) composite structures. Fiber damage and fiber irregularities in FRPs were detected by configuring hybrid sensors as coil sensors. Similarly, matrix damage, matrix irregularities and interlaminar delaminations were detected by configuring hybrid sensors as capacitive sensors. ECT and CI were performed sequentially by means of electronic switching. Hence, eliminating the need for mounting two separate sensor arrays on the probe. Portable robotic platform was developed by MSU for multi-technique high-speed NDE of CFRP test parts. The platform had two 6-axis robots, which enabled inspection of curved parts in approximately a 6×6×6 ft 3 active scan area. On the software side, the MSU team integrated scripts for NDE hardware control with scripts for robot motion control. MSU also implemented automated path planning for the robots, reconstruction of part’s surfaces via stereovision, 3D rendering of inspection data, and image processing algorithms for enhanced defect detection. Automotive composite parts manufactured by Plasan Composites from Phase I were used to validate the ACUT and EMT techniques on representative testbeds. Among those parts were three X-braces for a Dodge Viper, one composite calibration plaque with known defects at known locations, and four other test sections, including sections from a front splitter, a corner section from a composite hood, and a high-pressure RTM panel made using non crimp fabric. Other test samples included CFRP and GFRP calibration plates with fiber/matrix defects fabricated at MSU/CVRC.

36 MATERIALS SCIENCE↗

Effects of dose rate on the microstructure and deuterium retention in γ-LiAlO 2 pellets

This report presents experimental findings obtained from November 2024 to September 2025. Defect accumulation and microstructural evolution during ion irradiation at elevated temperatures are governed by two competing processes: defect production, driven by dose rate, and defect recovery, controlled by defect diffusion, interaction, and annihilation. At a given dose, the resulting microstructural evolution depends on both the dose rate and irradiation temperature. As a continuation of our FY24 tritium science project, which investigated temperature effects at a fixed dose rate, this study focuses on dose-rate effects at a fixed temperature to provide deeper insights into the irradiated microstructure and compositional changes in γ-LiAlO 2 pellets. The study aims to elucidate the impact of dose rate on microstructure, precipitate morphology, deuterium retention, and lithium volatility in γ-LiAlO 2 pellets under sequential 120 keV He + and 80 keV D 2 + ion irradiation. Three dose rates of 7.3×10 4 , 2.9×10 4 and 6.8×10 5 dpa/s were applied to achieve the same total ion fluence of 2×10 17 (He + +D + )/cm 2 at 500 °C, corresponding to a maximum combined dose of 7.55 dpa at ~255 nm. The irradiated pellets were subsequently characterized using scanning transmission electron microscopy (STEM) and time-of-flight secondary ion mass spectrometry (ToF-SIMS). The microstructural response of γ-LiAlO 2 to ion irradiation was found to be strongly dose-rate dependent. At medium and high dose rates, irradiation produced a surface amorphized layer and an underlying crystalline layer containing LiAl 5 O 8 precipitates, with implanted gases accumulating and forming blisters at the crystalline-amorphous interface. It remains to be investigated whether the amorphized layer was produced by He + ion irradiation prior to D 2 + ion irradiation. In contrast, at low dose rates, the material remained crystalline, with cavities, likely gas-filled, distributed around precipitates, within the γ-LiAlO 2 matrix, and along grain boundaries. While precipitate morphology exhibited anisotropy, their size showed little sensitivity to dose rate in the applied range of this study. This result, however, does not rule out the possibility that further reductions in dose rate could influence precipitate size. High dose-rate irradiation enhanced protonium–deuterium isotopic exchange and lithium depletion in the amorphized region. Collectively, the results show that dose rate governs amorphization, gas redistribution, isotopic exchange, and lithium depletion, providing important insights into the mechanisms underlying structural evolution in γ-LiAlO 2 pellets under reactor-relevant irradiation conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Optimization of gold germanium (Au 0.17 Ge 0.83 ) thin films for high sensitivity resistance thermometry

Gold–germanium (Au x Ge 1–x ) solid solutions have been demonstrated as highly sensitive thin film thermometers for cryogenic applications. However, little is known regarding the performance of the films for thicknesses less than 100 nm. In response, here we report on the resistivity and temperature coefficient of resistance (TCR) for sputtered films with thicknesses ranging from 10 to 100 nm and annealed at temperatures from 22 to 200°C. The analysis is focused upon composition x=0.17, which demonstrates a strong temperature sensitivity over a broad range. The thinnest films are found to provide an enhancement in TCR, which approaches 20% K –1 at 10 K. Furthermore, reduced anneal temperatures are required to crystallize the Ge matrix and achieve a maximum TCR for films of reduced thickness. These features favor the application of ultra-thin films as high-sensitivity, on-device thermometers in micro- and nanolectromechanical systems.

59 BASIC BIOLOGICAL SCIENCES↗

Extraction of low-energy constants of single- and double-β decays from lattice QCD: A sensitivity analysis

Lattice quantum chromodynamics (LQCD) has the promise of constraining low-energy constants (LECs) of nuclear effective field theories (EFTs) from first-principles calculations that incorporate the dynamics of quarks and gluons. Given the Euclidean and finite-volume nature of LQCD outputs, complex mappings are developed in recent years to obtain the Minkowski and infinite-volume counterparts of LQCD observables. In particular, as LQCD is moving toward computing a set of important few-nucleon matrix elements at the physical values of the quark masses, it is important to investigate whether the anticipated precision of LQCD spectra and matrix elements will be sufficient to guarantee tighter constraints on the relevant LECs than those already obtained from phenomenology, considering the nontrivial mappings involved. With a focus on the leading-order LECs of the pionless EFT, L 1,A and $g$$^{NN}_{v}$, which parametrize, respectively, the strength of the isovector axial two-body current in a single-β decay (and other related processes such pp fusion), and of the isotensor contact two-body operator in the neutrinoless double-β decay within the light neutrino exchange scenario, the expected uncertainty on future extractions of L 1,A and $g$$^{NN}_{v}$ are examined using synthetic data at the physical values of the quark masses. It is observed that achieving small uncertainties in L 1,A will be challenging, and (sub)percent-level precision in the two-nucleon spectra and matrix elements is essential in reducing the uncertainty on this LEC compared to the existing constraints. On the other hand, the short-distance coupling of the neutrinoless double-β decay, $g$$^{NN}_{v}$, is shown to be less sensitive to uncertainties on both LQCD energies and the matrix element, and can likely be constrained with percent-level precision in the upcoming LQCD calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

High-precision monitoring of outgassing species in model extreme ultraviolet photoresists with a cavity ring-down spectrometer

Background Extreme ultraviolet (EUV) photoresists play a pivotal role in advancing nanopatterning technologies by balancing image quality and sensitivity. The outgassing behavior of photoresist thin films under EUV and deep ultraviolet (DUV) exposure reveals chemical details relevant to their performance. Aim Here, we focus on utilizing an analytical technique not previously used in photolithography, the tabletop cavity ring-down spectrometer, to investigate outgassing dynamics in EUV photoresists, enabling precise chemical identification and deeper insights into resist processing. Approach The spectrometer’s enhanced laser path length (∼ 20 km) and broadband absorption capabilities in the C–H overtone region allow for sensitive and temporally resolved detection of mixtures of outgassed species. Using a model resist comprising a polymer matrix with a photoacid generator and quencher, we analyzed the influence of time delays between exposure and post-exposure bake (PEB) as well as storage under varying environmental conditions. Results Suppression of isobutylene outgassing and thickness loss was observed with extended delays between exposure and PEB, potentially linked to water absorption and acid deactivation. The technique proved highly effective in distinguishing subtle chemical differences between processing stages. Conclusions Delay times and their environmental conditions, particularly humidity, reduce outgassing and thickness loss of photoresists during PEB, suggesting decreased acid-driven deprotection. This can potentially impact sensitivity, defectivity, and roughness of resist patterns, necessitating precise monitoring and control.

Extreme ultraviolet photoresists↗

Hierarchical Bayesian Thermonuclear Rate for the 7 Be(n,p) 7 Li Big Bang Nucleosynthesis Reaction

Big Bang nucleosynthesis provides the earliest probe of standard model physics, at a time when the universe was less than 1000 seconds old. It determines the abundances of the lightest nuclides, which give rise to the subsequent history of the visible matter in the universe. This work derives new 7 Be(n,p) 7 Li thermonuclear reaction rates based on all available experimental information. This reaction sensitively impacts the primordial abundances of 7 Be and 7 Li during big bang nucleosynthesis. We critically evaluate all available data and disregard experimental results that are questionable. For the nuclear model, we adopt an incoherent sum of single-level, two-channel, R -matrix approximation expressions, which are implemented into a hierarchical Bayesian model, to analyze the remaining six data sets we deem most reliable. In the fitting of the data, we consistently model all known sources of uncertainty, including discrepant absolute normalizations of different data sets, and also take the variation of the neutron and proton channel radii into account, hence providing less biased estimates of the 7 Be(n,p) 7 Li thermonuclear rates. From the resulting posteriors, we extract R -matrix parameters ($E_r, γ^2_n, γ^2_p$) and derive excitation energies and partial and total widths. Our fit is sensitive to the contributions of the first three levels above the neutron threshold. Reaction rates were computed by integrating 10,000 samples of the reduced cross section. Our 7 Be(n,p) 7 Li thermonuclear rates have uncertainties between 1.5% and 2.0% at temperatures of ≤1 GK. Finally, we compare our rates to previous results and find that the 7 Be(n,p) 7 Li rates most commonly used in big bang simulations have uncertainties that are too optimistic.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Low-rank Sachdev-Ye-Kitaev models

Motivated by recent proposals of experimental realization of fast scramblers, we study a family of solvable variants of the ($\textit{q}$ = 4) Sachdev-Ye-Kitaev model in which the rank and eigenvalue distribution of the coupling matrix $J_{ ij,kl}$ are tuneable. When the rank is proportional to the number of fermions, the low temperature behavior is sensitive to the eigenvalue distribution. We obtain a complete classification of the possible non-Fermi liquid quantum phases. These include two previously studied phases whose fermion scaling dimension depends continuously on the rank; we show that they are maximally chaotic, but necessitate an extensively degenerate or negative semidefinite coupling matrix. In this work, more generic distributions give rise to “almost Fermi liquids” with a scaling dimension Δ = 1/2, but which differ from a genuine Fermi liquid in quasiparticle decay rate, quantum Lyapunov exponent, and/or specific heat.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A simple and practical wax-encapsulation method for air-sensitive XAS samples

To facilitate X-ray absorption spectroscopy (XAS) measurements of air-sensitive samples, we present a simple method in which materials are encased in common paraffin wax to protect them from air and moisture. We demonstrate the efficacy of this approach using a highly reducing, air- and moisture-sensitive uranium(III) complex, the tris(amide) U[N(SiMe 3 ) 2 ] 3 (1). When finely dispersed in a boron nitride matrix and subsequently encased in inert paraffin wax, samples of 1 remain stable with no visible or spectroscopic degradation after several days under ambient conditions. The viability of this method for XAS measurements was further evaluated across a series of uranium compounds, ranging from uranyl species to highly air- and moisture-sensitive molecular complexes, at the uranium L 3 -edge. Edge energy determinations were highly reproducible (±0.1 eV between replicates) and, where available, showed excellent agreement with literature values. In conclusion, this low-cost, effective, and versatile method offers a viable solution for XAS studies of air-sensitive compounds and materials.

36 MATERIALS SCIENCE↗

Three-nucleon lepton-number-violating potentials in chiral effective field theory and their matrix elements in light nuclei

Here, we derive the three-nucleon neutrinoless double-𝛽 decay potential in a Δ-full chiral effective field theory through next-to-next-to-next-to leading order in Weinberg's power counting. The matrix elements of the resulting operators are computed in light nuclei using variational Monte Carlo with wave functions constructed from the Norfolk family of nuclear interactions. We find that three-nucleon corrections induce a modest quenching of the total nuclear matrix elements. We discuss model dependencies and the potential impact of these corrections on the sensitivity of experimental programs to probe lepton number violating parameters. These results provide a benchmark for many-body methods capable of reaching heavier nuclei of experimental interest.

Chambers-Wall, Graham [Washington University, St. ↗

Electric-field control of spin dynamics during magnetic phase transitions

Controlling magnetization dynamics is imperative for developing ultrafast spintronics and tunable microwave devices. However, the previous research has demonstrated limited electric-field modulation of the effective magnetic damping, a parameter that governs the magnetization dynamics. Here, we propose an approach to manipulate the damping by using the large damping enhancement induced by the two-magnon scattering and a nonlocal spin relaxation process in which spin currents are resonantly transported from antiferromagnetic domains to ferromagnetic matrix in a mixed-phased metallic alloy FeRh. This damping enhancement in FeRh is sensitive to its fraction of antiferromagnetic and ferromagnetic phases, which can be dynamically tuned by electric fields through a strain-mediated magnetoelectric coupling. In a heterostructure of FeRh and piezoelectric PMN-PT, we demonstrated a more than 120% modulation of the effective damping by electric fields during the antiferromagnetic-to-ferromagnetic phase transition. Our results demonstrate an efficient approach to controlling the magnetization dynamics, thus enabling low-power tunable electronics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evaluation of Extensional Stiffness Matrix [A]. of Braided Composite Tubes

A novel method for analyzing tubular textile composites was developed by obtaining a curved repeatable unit cell (RUC) from the transformation of a flat RUC. Periodic boundary conditions were imposed on the tubular RUC. The stiffness matrix of the composite was then evaluated using transversely isotropic yarn and isotropic matrix. It was found that the stiffness coefficients A11, A12 of the composite are more sensitive to a change in the radius of the composite compared to A22 and A66.

Thandaga Nagaraju, Hematnh↗

Impacts of Crystalline Host Rock on Repository Barrier Materials at 250 °C: Hydrothermal Co-Alteration of Wyoming Bentonite and Steel in the Presence of Grimsel Granodiorite

Direct disposal of dual-purpose canisters (DPC) has been proposed to streamline the disposal of spent nuclear fuel. However, there are scenarios where direct disposal of DPCs may result in temperatures in excess of the specified upper temperature limits for some engineered barrier system (EBS) materials, which may cause alteration within EBS materials dependent on local conditions such as host rock composition, chemistry of the saturating groundwaters, and interactions between barrier materials themselves. Here we report the results of hydrothermal experiments reacting EBS materials—bentonite buffer and steel—with an analogue crystalline host rock and groundwater at 250 °C. Experiment series explored the effect of reaction time on the final products and the effects of the mineral and fluid reactants on different steel types. Post-mortem X-ray diffraction, electron microprobe, and scanning electron microscopy analyses showed characteristic alteration of both bentonite and steel, including the formation of secondary zeolite and calcium silicate hydrate minerals within the bentonite matrix and the formation of iron-bearing clays and metal oxides at the steel surfaces. Swelling clays in the bentonite matrix were not quantitatively altered to non-swelling clay species by the hydrothermal conditions. The combined results of the solution chemistry over time and post-mortem mineralogy suggest that EBS alteration is more sensitive to initial groundwater chemistry than the presence of host rock, where limited potassium concentration in the solution prohibits conversion of the smectite minerals in the bentonite matrix to non-swelling clay species.

54 ENVIRONMENTAL SCIENCES↗

Effects of Strain Rate and Temperature on the Mechanical Properties of Simulated Silica Ionogels

Ionogels are hybrid materials formed by impregnating the pore space of a solid matrix with a conducting ionic liquid. By combining the properties of both component materials, ionogels can act as self-supporting electrolytes in Li batteries. In this study, molecular dynamics simulations are used to investigate the dependence of mechanical properties of silica ionogels on solid fraction, temperature, and pore width. Comparisons are made with corresponding aerogels. We find that the solid matrix fraction increases the moduli and strength of the ionogel. This varies nonlinearly with temperature and strain rate, according to the contribution of the viscous ionic liquid to resisting deformation. Additionally, owing to the temperature and strain sensitivity of the ionic liquid viscosity, the mechanical properties approach a linear mixing law at high temperature and low strain rates. The median pore width of the solid matrix plays a complex role, with its influence varying qualitatively with deformation mode. Narrower pores increase the relevant elastic modulus under shear and uniaxial compression but reduce the modulus obtained under uniaxial tension. Conversely, shear and tensile strength are increased by narrowing the pore width. All of these pore size effects become more pronounced as the silica fraction increases. Pore size effects, similar to the effects of temperature and strain rate, are linked to the ease of fluid redistribution within the pore space during deformation-induced changes in the geometry of the pores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗