SEARCH · Engineering Papers

Results for “coupled”

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 685 records · Page 38

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Matter-wave stripe solitons induced by helicoidal spin–orbit coupling

Highlights: • Four kinds of exact stripe soliton solutions in Bose–Einstein condensates with helicoidal spin–orbit coupling (SOC) are represented. • The influence of helicoidal SOC on the structure of stripe solitons are described in detail. • The corresponding spin polarization is expressed on the Bloch sphere. • The formation mechanism of stripe solitons and the interaction properties of moving stripe solitons are discussed. We present the exact analytical solution of stripe solitons on different backgrounds in Bose–Einstein condensates with helicoidal spin–orbit coupling (SOC). They are divided into two categories, one is stationary stripe soliton on zero background or period-wave background, the other is moving stripe soliton on zero background or plane-wave background. The period of stripe is determined by both strength and frequency of helicoidal SOC, and the oscillation amplitude of stripe is adjustable in the case of fixed total density. The corresponding spin polarization is expressed on Bloch sphere. For stationary soliton, the spin density vector take values in three directions, which is different from the Rashba–Dresslhaus coupling scheme. The spin density vectors display more complicated spatial distribution for moving stripe solitons. Meanwhile, the formation mechanism of stripe soliton and the interaction properties of moving stripe solitons are discussed. Our work enhances the understanding of stripe solitons, and is helpful for the observation of stripe phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Adiabatic quantum decoherence in many non-interacting subsystems induced by the coupling with a common boson bath

Highlights: • System–environment quantum correlation: a main solid state NMR decoherence channel. • Non-separable system–environment model yields realistic spin decoherence rates. • New open quantum system approach explains irreversible decay of refocused NMR echoes. • Adiabatic quantum decoherence is inherently irreversible and eigen-selective. This work addresses adiabatic quantum decoherence of many-body spin systems coupled with a boson field in the framework of open quantum systems theory. We generalize the traditional spin-boson model by considering a system–environment interaction Hamiltonian that represents a partition of non-interacting subsystems and highlights the collective correlation that appears exclusively due to the coupling with a common environment. Remarkably, this simple, exactly solvable model encompasses relevant aspects of a many-body open quantum system and features the subtle quantum effects that arise when the size scales up to a macroscopic level. We derive an analytical expression for the time dependence of the density matrix elements (in the preferred basis) without assuming coarse-graining. The resulting decoherence function is eigen-selective and is a complex exponential whose exponent has a real part that introduces a decay similar to that in the spin-boson model. On the contrary, the imaginary part depends on the quantum numbers and geometry of the whole partition and does not reflect the system temperature. Motivated by decoherence in solid-state NMR, and in search of realistic numerical estimations, we apply the theoretical results to a partition of dipole-coupled spin pairs in contact with a common phonon bath, using typical parameters of hydrated salts. The proposal allows estimating the decoherence time scale in terms of the system physical constants: sound velocity and eigenvalue distribution width. As a significant novelty, the decoherence function phase depends on the eigenvalue distribution throughout the sample. It plays the leading role, overshadowing the mechanism associated with the bath thermal state. Finally, we apply the formalism to describe decoherence in the “magic echo” NMR reversal experiment. We find that the system–environment correlation explains the origin of irreversibility, and both the decoherence rate value and its dependence on the dipolar frequency, are remarkably similar to the experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Simulation and power quality analysis of a Loose-Coupled bipolar DC microgrid in an office building

With distributed generation and battery storage technologies thriving in microgrids, the use of direct current (DC) microgrids in the building sector offers multiple advantages in energy efficiency and power quality compared with alternating current (AC) systems. This study developed a new concept of a loose-coupled bipolar DC building power system. In this work, the concept was used to design a real-world office building in Shenzhen, China. A power system model was developed to study the stability and control of the DC power system and to verify DC power quality. The design and modeling of the DC power control system is discussed in detail. The study developed a few common fault scenarios in DC building microgrids that were simulated in the MATLAB-Simulink environment to validate the design of a loose-coupled bipolar DC system. The results indicate that the proposed loose-coupled bipolar DC system schema, when implemented with proper control algorithms, can achieve good fault-tolerant performance with reliable power quality, even during disruptive system events.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2022

Coupling between DNA replication, segregation, and the onset of constriction in Escherichia coli

Escherichia coli cell cycle features two critical cell-cycle checkpoints: initiation of replication and the onset of constriction. While the initiation of DNA replication has been extensively studied, it is less clear what triggers the onset of constriction and when exactly it occurs during the cell cycle. Here, using high-throughput fluorescence microscopy in microfluidic devices, we determine the timing for the onset of constriction relative to the replication cycle in different growth rates. Our single-cell data and modeling indicate that the initiation of constriction is coupled to replication-related processes in slow growth conditions. Furthermore, our data suggest that this coupling involves the mid-cell chromosome blocking the onset of constriction via some form of nucleoid occlusion occurring independently of SlmA and the Ter linkage proteins. This work highlights the coupling between replication and division cycles and brings up a new nucleoid mediated control mechanism in E. coli.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

A continuum framework for coupled solid deformation–fluid flow through anisotropic elastoplastic porous media

We present a continuum framework for coupled solid deformation–fluid flow in anisotropic elastoplastic porous media. A thermodynamic formulation of the coupled processes gives rise to an anisotropic Biot tensor that is a function of drained elastic tangent moduli tensor of the solid skeleton and the intrinsic bulk modulus of the solid constituent. Two effective stress measures emerge from the formulation, namely, σ', which is energy-conjugate to the elastic strain, and σ'', which is energy-conjugate to the plastic strain. For the special case of transverse isotropy that is commonly encountered in natural rocks, the Biot tensor can be expressed in terms of its normal and tangential components to the bedding plane, along with a microstructure tensor. Apart from its thermodynamic consistency, an advantage of this new formulation is that standard mixed finite element formulation can be employed to discretize the domain and solve initial boundary-value problems. We conduct plane strain simulations of coupled solid deformation–fluid flow in a transversely isotropic porous medium to demonstrate the impacts of material anisotropy, stress history, and the Biot tensor on the system response.

42 ENGINEERING↗

DOE OSTI2022

Enriched immersed boundary method (EIBM) for interface-coupled multi-physics and applications to convective conjugate heat transfer

It is challenging to develop numerical methods that simultaneously maintain accuracy of resolving boundary conditions and mesh flexibility to handle the interface in interface-coupled multi-physics problems involving large property discontinuity. On the one hand, boundary-fitted methods possess high accuracy in capturing the interfacial phenomenon but involve complicated volumetric mesh generation, mesh-motion, and even re-meshing procedures. On the other hand, immersed boundary methods (IBM) provide mesh flexibility but sometimes suffer from inferior interface representations, leading to poor enforcement of boundary conditions. This paper presents an enriched immersed boundary method (EIBM) to overcome this challenge and demonstrates its efficacy in conjugate heat transfer, a representative example in interface-coupled multi-physics systems that have implications for many industrial processes. Further, the core technique of the method is to enhance IBM’s accuracy of the fluid–solid interface by enriching the degrees of freedom of the cut elements to enforce temperature and flux compatibilities and resolve all the physical unknowns on the background mesh to simplify the volumetric mesh generation. We implement the EIBM under the framework of a variational multiscale formulation for coupled Navier–Stokes and thermodynamics equations. The enriched DoFs enable better enforcement of temperature and flux compatibilities with large conductivity ratios across the fluid–solid interface. At the same time, the immersed nature of the proposed method still attains mesh flexibility. The EIBM’s accuracy is thoroughly evaluated through a set of examples, ranging from benchmark problems with analytical solutions to real-world cooling processes of a moving metallic structure with complex geometry.

42 ENGINEERING↗

DOE OSTI2026

Flame dynamics and kinetic coupling of ammonia and dimethyl-ether in non-premixed cool and warm flames at elevated pressure

Developing advanced low-temperature combustion engines with ammonia-biofuel blends requires a comprehensive understanding of low-temperature flame dynamics and kinetic interactions between ammonia and oxygenated fuels at elevated pressures. This work aims to study the dynamics and kinetics of non-premixed Dimethyl Ether (DME)/Ammonia (NH 3 ) cool and warm flames, and their reignition to hot flames. A counterflow burner is employed to establish DME/NH 3 cool/warm flames at pressures up to 5 atm. The extinction limits of cool flame and the reignition limits of warm flame to hot flame are measured by varying NH 3 concentrations and compared to simulations to quantitatively examine the effects on DME/NH 3 flames. It is found that NH 3 inhibits low-temperature DME oxidation and results in lower cool flame extinction limits. Warm flames in the presence of NH 3 are observed for the first time, revealing a non-monotonic effect of NH 3 addition: a small amount of NH 3 presence enhances warm flame chemistry and promotes reignition to hot flames, while a high NH 3 concentration weakens the warm flame. This trend is further explained by 0-D PSR kinetic simulations and 1-D S-curve flame dynamic calculations. Three flame transition regimes between cool flames (CF), warm flames (WF), and hot flames (HF) by different levels of NH 3 additions at a specific strain rate are identified, namely WFHF reignition, WF-CF transition, and WF extinction. Reaction sensitivity analyses of OH at low temperatures show that NH3 inhibits DME oxidation through OH consumption via H-abstraction and the kinetic couplings of RO 2 /NH 2 , RO 2 /NO x , R/NO x , and O 2 QOOH/NO x further suppress the low-temperature branching. At intermediate-temperatures, NH 2 /NO x /HO 2 coupling promotes warm flames via the pathway NH 2 → H 2 NO → HNO → NO by converting O 2 → HO 2 → OH. At even higher NH₃ concentrations, radical termination reactions of NH 2 + NO/NO 2 and excessive OH consumption via H-abstraction inhibit the flame. The insights into the kinetic coupling between NH 3 and low-temperature chemistry at elevated pressure and its impact on the dynamics of cool-warm-hot flame transitions will contribute to advancing combustion technologies with reduced emissions and improved energy-efficiency.

42 ENGINEERING↗

DOE OSTI2021

Improved adhesion of polyurethane-based nanocomposite coatings to tin surface through silane coupling agents

Conformal coatings with higher adhesion strength is critical for mitigating the growth of whiskers from the pure tin (Sn) or Sn-rich surfaces of the Pb-free electronics. Silane coupling agents have the ability to form a durable bond between polymeric coatings and metallic materials. The main goal of this paper is to investigate the effect of silane agents on the adhesion behavior of polyurethane (PU), polyurethane acrylate (PUA) and both coatings filled with silicone rubber (SR) nanoparticles on the tin solder pads. It was found that an isocyanate-based silane provided excellent adhesion for PU and PU + SR coatings and an acrylic-based silane for PUA and PUA + SR coatings on the tin surface. This was due to the reaction bond formation between the silane coupling agent and the coating. Hence, the selection of silane coupling agents with the organofunctional group similar to the coating material is crucial to improve the adhesion properties of the coating.

36 MATERIALS SCIENCE↗

DOE OSTI2021

TOUGH-RFPA: Coupled thermal-hydraulic-mechanical Rock Failure Process Analysis with application to deep geothermal wells

This paper presents a coupled thermal-hydraulic-mechanical (THM) model for simulating failure processes in geological media. The coupling between stress/strain and heat transfer describing how thermal expansion affects stress is based on the theory of thermo-elasticity. Additionally, the coupling between stress/strain and fluid flow describing how pore-pressure affects stress is governed by Biot's consolidation theory. Stress/strain-dependent porosity/permeability/capillary-pressure defines the influence of stress on fluid flow and thermal fields. The THM model is implemented in the solver module that links the multiphase, multicomponent fluid flow simulator TOUGH2 with the geomechanical simulator RFPA (Rock Failure Process Analysis), which is based on the concept that heterogeneity leads to non-linearity and progressive failure behavior. Capability of TOUGH-RFPA simulator was verified by two benchmark problems related to hydraulic-mechanical and thermal-mechanical response in a hollow cylinder. Then the applicability of TOUGH-RFPA was demonstrated for modeling cooling-induced damage around a very deep geothermal well with mesoscopic heterogeneous properties, which showed complex damage evolution and fracturing during well cooling and fluid pressurization.

42 ENGINEERING↗

DOE OSTI2020

Efficient bi-directional coupling of 3D computational fluid-particle dynamics and 1D Multiple Path Particle Dosimetry lung models for multiscale modeling of aerosol dosimetry

The development of predictive aerosol dosimetry models has been a major focus of environmental toxicology and pharmaceutical health research for decades. One-dimensional (1D) models successfully predict overall deposition averages but fail to accurately predict local deposition. Computational fluid-particle dynamics (CFPD) models provide site-specific predictions but at a computational cost that prohibits whole lung predictions. Thus, there is a need for developing multiscale strategies to provide a realistic subject-specific picture of the fate of inhaled aerosol in the lungs. CT-based 3D/CFPD models of the large airways were bidirectionally coupled with individualized 1D Navier-Stokes airflow and particle transport based upon the widely used Multiple Path Particle Dosimetry Model (MPPD). Distribution of airflows among lobes was adjusted by measured lobar volume changes observed in CT images between FRC and FRC + 1.5 L. Additionally, as a test of the effectiveness of the coupling procedures, deposition modeling of previous 1 µm aerosol exposure studies was performed. The complete coupled model was run for 3 breaths, with the computation-intense portion being the 3D CFD Lagrangian particle tracking calculation. The average deposition per breath was 11% in the combined multiscale model with site-specific doses available in the CFPD portion of the model and airway- or region-specific deposition available for the MPPD portion. In conclusion, the key methods developed in this study enable predictions of ventilation heterogeneities and aerosol deposition across the lungs that are not captured by 3D or 1D models alone. Overall, these methods can be used as the foundation for multi-scale modeling of the full respiratory system.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Tusas: A fully implicit parallel approach for coupled phase-field equations

In this study, we develop a fully-coupled, fully-implicit approach for phase-field modeling of solidification in metals and alloys. Predictive simulation of solidification in pure metals and metal alloys remains a significant challenge in the field of materials science, as microstructure formation during the solidification process plays a critical role in the properties and performance of the solid material. Our simulation approach consists of a finite element spatial discretization of the fully-coupled nonlinear system of partial differential equations at the microscale, which is treated implicitly in time with a preconditioned Jacobian-free Newton-Krylov method. The approach is algorithmically scalable as well as efficient due to an effective preconditioning strategy based on algebraic multigrid and block factorization. We implement this approach in the open-source Tusas framework, which is a general, flexible tool developed in C++ for solving coupled systems of nonlinear partial differential equations. The performance of our approach is analyzed in terms of algorithmic scalability and efficiency, while the computational performance of Tusas is presented in terms of parallel scalability and efficiency on emerging heterogeneous architectures. We demonstrate that modern algorithms, discretizations, and computational science, and heterogeneous hardware provide a robust route for predictive phase-field simulation of microstructure evolution during additive manufacturing.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗

DOE OSTI2022

Interfacial magnetic vortex formation in exchange-coupled hard-soft magnetic bilayers

The exchange coupling between a hard magnetic layer MnBi and a soft magnetic layer Co-Fe has been found to significantly improve the maximum energy product. In this work, the spin structure of exchange coupled MnBi:Co-Fe bilayers is experimentally investigated by X-ray magnetic circular dichroism (XMCD) and polarized neutron reflectometry (PNR). We find that the out-of-plane magnetization reversal process of the MnBi:Co-Fe bilayer structure involves formation of a curling-type twisting of the magnetization in the film plane at low or intermediate reversal fields. Additionally, micromagnetic simulations are further performed to provide a detailed view of the spins at the curling center. Reminiscent of chiral spin structures known as spin bobbers, this curling in the exchange-coupled hard-soft magnetic bilayers is a new type of skyrmionic spin structure and worth further investigation.

36 MATERIALS SCIENCE↗

DOE OSTI2022

A concurrent irradiation-mechanics multiscale coupling model

The vast majority of our current knowledge regarding the basic mechanisms controlling irradiation effect on mechanical properties is based almost entirely on the results of post-irradiation experiments or theoretical models. However, the concurrent effects of irradiation, mechanical stress, and thermal damage on the failure phenomena of materials and components remain largely unexplored due to its internal multiscale-multiphysics coupling nature. Here, we present here a concurrent irradiation-mechanics multiscale coupling model. The concurrent evolutions of nanoscale irradiation defect clusters, microscale dislocation configurations, and mechanical responses are well captured based on coupling cluster dynamics, discrete dislocation dynamics, and the finite element methods using an effective time marching scheme. Model predictions of defect densities and size are in general agreement with experimental observations. Irradiation hardening is shown to take place also in samples undergoing concurrent irradiation-mechanical loading, similar to samples tested post-irradiation. However, the occurrence of plastic flow localization and dislocation channel formation is not accompanied with apparent yield drop (softening) under concurrent irradiation-mechanical loading conditions, which is different from the post-irradiation case.

42 ENGINEERING↗

DOE OSTI2022

Double echo symmetry-based REDOR and RESPDOR pulse sequences for proton detected measurements of heteronuclear dipolar coupling constants

1 H{X} symmetry-based rotational echo double resonance pulse sequences (S-REDOR) and symmetry-based rotational echo saturation pulse double resonance (S-RESPDOR) solid-state NMR experiments have found widespread application for 1 H detected measurements of difference NMR spectra, dipolar coupling constants, and internuclear distances under conditions of fast magic angle spinning (MAS). In these experiments the supercycled $\text{R4}^2_1 (\text{SR4}^2_1)$ symmetry-based recoupling pulse sequence is typically applied to the 1 H spins to reintroduce heteronuclear dipolar couplings. However, the timing of $\text{SR4}^2_1$ and other symmetry-based pulse sequences must be precisely synchronized with the rotation of the sample, otherwise, the evolution of 1 H CSA and other interactions will not be properly refocused. For this reason, significant distortions are often observed in experimental dipolar dephasing difference curves obtained with S-REDOR or S-RESPDOR pulse sequences. Here we introduce a family of double echo (DE) S-REDOR/S-RESPDOR pulse sequences that function in an analogous manner to the recently introduced $t_1$-noise eliminated (TONE) family of dipolar heteronuclear multiple quantum coherence (D-HMQC) pulse sequences. Through numerical simulations and experiments the DE S-REDOR/S-RESPDOR sequences are shown to provide dephasing difference curves similar to those obtained with S-REDOR/S-RESPDOR. However, the DE sequences are more robust to the deviations of the MAS frequency from the ideal value that occurs during typical solid-state NMR experiments. In this work, the DE sequences are shown to provide more reliable 1 H detected dipolar dephasing difference curves for nuclei such as 15 N (with isotopic labelling), 183 W and 35 Cl. The double echo sequences are therefore recommended to be used in place of conventional S-REDOR/S-RESPDOR sequences for measurement of weak dipolar coupling constants and long-range distances.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Observation of large Rashba spin–orbit coupling at room temperature in compositionally engineered perovskite single crystals and application in high performance photodetectors

Indirect absorption extended below the direct transition edge and increase in carrier lifetime derived from Rashba spin–orbit coupling may advance the optoelectronic applications of metal halide perovskites. Spin-orbit coupling in halide perovskites is due to the presence of heavy elements in their structure. However, when these materials lack an inversion symmetry, for example by the application of strain, spin–orbit coupling becomes odd in the electron’s momentum giving rise to a splitting in the electronic energy bands. Here we report on the observation of a large Rashba splitting of 117 meV at room temperature, as predicted by relativistic first-principles calculations, in halide perovskite single crystals through a facile compositional engineering approach. Partial substitution of organic cations by rubidium in single crystals induces significant indirect absorption and dual peak photoluminescence as a result of a large Rashba splitting. We measured circularly polarized photoluminescence and magneto-photoluminescence in perovskite films printed by single crystals as well as magneto-electroluminescence and magneto-photocurrent in spin-LEDs based on perovskite single crystals. They indicated significant spin-momentum locking due to the large Rashba effect. A hybrid perovskite single crystal photodetector achieved record figures of merit, including detectivity of more than 1.3 × 10 18 Jones which represents a three orders of magnitude improvement compared to the to date record. Furthermore, these findings show that facile compositional engineering of perovskite single crystals holds great promise for further advancing the optoelectronic properties of existing materials.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Resistive AC-Coupled Silicon Detectors: principles of operation and first results from a combined analysis of beam test and laser data

This paper presents the principles of operation of Resistive AC-Coupled Silicon Detectors (RSDs) and measurements of the temporal and spatial resolutions using a combined analysis of laser and beam test data. RSDs are a new type of n-in-p silicon sensor based on the Low-Gain Avalanche Diode (LGAD) technology, where the n+ implant has been designed to be resistive, and the read-out is obtained via AC-coupling. The truly innovative feature of RSD is that the signal generated by an impinging particle is shared isotropically among multiple read-out pads without the need for floating electrodes or an external magnetic field. Careful tuning of the coupling oxide thickness and the n+ doping profile is at the basis of the successful functioning of this device. Several RSD matrices with different pad width-pitch geometries have been extensively tested with a laser setup in the Laboratory for Innovative Silicon Sensors in Torino, while a smaller set of devices have been tested at the Fermilab Test Beam Facility with a 120 GeV/c proton beam. The measured spatial resolution ranges between 2.5μm for 70–100 pad-pitch geometry and 17μm with 200–500 matrices, a factor of 10 better than what is achievable in binary read-out (bin size/12). Beam test data show a temporal resolution of ~40ps for 200 μm pitch devices, in line with the best performances of LGAD sensors at the same gain.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

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↗