Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent scheme”

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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Hypervelocity impact facility for simulating materials exposure to impact by space debris

The Space Power Institute at Auburn University has constructed an electromagnetically driven particle accelerator for simulating the effects of space debris on the materials for use in advanced spacecraft. The facility consists of a capacitively driven accelerator section, a drift tube and a specimen impact chamber. The drift tube is sufficiently long that all electrical activity has ceased prior to impact in the specimen chamber. The impact chamber is large enough to allow a wide range of specimen geometries, ranging from small coupons to active portions of advanced spacecraft. The electric drive for the accelerator consists of a 67 kJ, 50 k capacitor bank arranged in a low inductance configuration. The bank is discharged through an aluminum armature/plastic ablator plate/projectile load in roughly 1.2 microsec. The evaporation of the ablaitor plate produces an expanding gas slug, mostly H2, traveling at a velocity of some 60 km/sec. Because of the pressure and local density, the expanding gas cloud accelerates projectiles due to plasma drag. To date, we have utilized projectiles consisting of 100 micron SiC, 100 and 400 micron Al2O3, 100 and 145 micron olivines. Since many particles are accelerated in a given experiment, there is a range of velocities for each shot as well as some particle breakup. Advanced diagnostics techniques allow determination of impact coordinates, velocity, and approximate size for as many as 50 individual impacts in a given experiment. We routinely measure velocities in the range 1-15 km/sec. We have used this facility to study a variety of impact generated phenomena on coated surfaces, both paint and plastic, thermal blanket material, solar cell arrays, and optical materials such as glass and quartz lenses. The operating characteristics of the gun, the advanced diagnostic scheme, and the results of studies of crater morphology are described in detail. Projectile residue analysis, as a function of impact velocity for the materials listed above, is also discussed. Wherever possible, these results are compared to those obtained by LDEF investigators and future experiments suggested which could help to explain unique features associated with LDEF impacts.

Rose, M. Frank↗
NASA NTRS2001

Ocean Color Optical Property Data Derived from OCTS and POLDER: A Comparison Study

We describe our efforts in studying and comparing the ocean color data derived from the Japanese Ocean Color and Temperature Scanner (OCTS) and the French Polarization and Directionality of the Earth's Reflectances (POLDER). OCTS and POLDER were both on board Japan's Sun-synchronous Advanced Earth Observing Satellite (ADEOS-1) from August 1996 to June 1997, collecting about 10 months of global ocean color data. This provides a unique opportunity for developing methods and strategies for the merging of ocean color data from multiple ocean color sensors. In this paper, we describe our approach in developing consistent data processing algorithms for both OCTS and POLDER and using a common in situ data set to vicariously calibrate the two sensors. Therefore, the OCTS and POLDER-measured radiances are effectively bridged through common in situ measurements. With this approach in processing data from two different sensors, the only differences in the derived products from OCTS and POLDER are the differences inherited from the instrument characteristics. Results show that there are no obvious bias differences between the OCTS and POLDER-derived ocean color products, whereas the differences due to noise, which stem from variations in sensor characteristics, are difficult to correct. It is possible, however, to reduce noise differences with some data averaging schemes. The ocean color data from OCTS and POLDER can therefore be compared and merged in the sense that there is no significant bias between two.

Wang, Menghua↗
NASA NTRS2006

In-Space Radiator Shape Optimization using Genetic Algorithms

Future space exploration missions will require the development of more advanced in-space radiators. These radiators should be highly efficient and lightweight, deployable heat rejection systems. Typical radiators for in-space heat mitigation commonly comprise a substantial portion of the total vehicle mass. A small mass savings of even 5-10% can greatly improve vehicle performance. The objective of this paper is to present the development of detailed tools for the analysis and design of in-space radiators using evolutionary computation techniques. The optimality criterion is defined as a two-dimensional radiator with a shape demonstrating the smallest mass for the greatest overall heat transfer, thus the end result is a set of highly functional radiator designs. This cross-disciplinary work combines topology optimization and thermal analysis design by means of a genetic algorithm The proposed design tool consists of the following steps; design parameterization based on the exterior boundary of the radiator, objective function definition (mass minimization and heat loss maximization), objective function evaluation via finite element analysis (thermal radiation analysis) and optimization based on evolutionary algorithms. The radiator design problem is defined as follows: the input force is a driving temperature and the output reaction is heat loss. Appropriate modeling of the space environment is added to capture its effect on the radiator. The design parameters chosen for this radiator shape optimization problem fall into two classes, variable height along the width of the radiator and a spline curve defining the -material boundary of the radiator. The implementation of multiple design parameter schemes allows the user to have more confidence in the radiator optimization tool upon demonstration of convergence between the two design parameter schemes. This tool easily allows the user to manipulate the driving temperature regions thus permitting detailed design of in-space radiators for unique situations. Preliminary results indicate an optimized shape following that of the temperature distribution regions in the "cooler" portions of the radiator. The results closely follow the expected radiator shape.

Hull, Patrick V.↗
DOE OSTI2025

A segmented approach to modeling building height: Delineating high-rise and low-rise buildings for enhanced height estimation

Understanding building height is imperative to the overall study of energy efficiency, population distribution, urban morphologies, emergency response, among others. Currently, existing approaches for modeling building height at scale are hindered by two pervasive issues. First, there is no consistent approach to quantify what a high-rise building is at a macro scale, leaving researchers unable to accurately compare results across geographies and domains. Second, high-rise buildings represent a small fraction of the built environment, implying data imbalance challenges that negatively affect current approaches. This is a problem of practical relevance since information on high-rise buildings is important for studies on urban heat islands, population dynamics, and pollution dispersion. Here, we introduce a novel approach to map building height which first identifies two distinct distributions within the built environment, with one being composed of low-rise buildings and one composed of high-rise buildings. We then develop an ensemble scheme where discrete specialist models are trained for each subset of low-rise buildings and high-rise buildings to infer building height from morphology features. For experiments mapping heights of 4.85 million buildings in Japan, we show an increase of 34 % in accuracy within 3m error when compared to the current state-of-the-art when modeling high-rise buildings, which based on KNN experimentation we define as any building > 12m . Our findings show that such an ensemble framework outperforms the current state-of-the-art approaches, which is especially relevant in relation to inferring height for high-rise buildings, a prominent issue of existing approaches for mapping the built environment.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Coupled cluster and dislocation dynamics modeling of microstructure evolution in irradiated materials

We develop here a coupled cluster and dislocation dynamics framework to study the microstructure evolution of irradiated materials. The framework not only accounts for the three dimensional diffusion of radiation-generated clusters, but also their interaction with dislocation networks and the resultant climb motion of discrete dislocations within finite crystals. The framework is solved with a superposition solution scheme, and is applied to investigate the evolution of the irradiation-induced dislocation loops in zirconium (Zr), considering the effects of various bias factors including the diffusion anisotropy difference (DAD) of interstitials and interstitial clusters, the dislocation bias of defects to discrete dislocation segments, and the production bias of defects from the radiation cascade. We find that the DAD is the most critical factor influencing the kinetics of the loop evolution in Zr, while the recombination/interaction of mobile defects can induce a strong spatial dependence of the loop evolution together with the DAD. Here, the method is also adopted to study the evolution of interstitial $\langle$a$\rangle$ and vacancy $\langle$c$\rangle$ dislocation loop ensembles consistent with the microstructure observed during irradiation-induced growth of Zr. Our findings not only reveal the spatial dependence of the size and ellipticity of the dislocation loops, but also suggest a limit on the anisotropy factor of interstitials to reproduce the co-growth of $\langle$a$\rangle$ and $\langle$c$\rangle$ loops in zirconium, in good agreement with experimental observations and other simulation results.

Bias factors↗
DOE OSTI2023

Electron Scattering from 1-Methyl-5-Nitroimidazole: Cross-Sections for Modeling Electron Transport through Potential Radiosensitizers

In this study, we present a complete set of electron scattering cross-sections from 1-Methyl-5-Nitroimidazole (1M5NI) molecules for impact energies ranging from 0.1 to 1000 eV. This information is relevant to evaluate the potential role of 1M5NI as a molecular radiosensitizers. The total electron scattering cross-sections (TCS) that we previously measured with a magnetically confined electron transmission apparatus were considered as the reference values for the present analysis. Elastic scattering cross-sections were calculated by means of two different schemes: The Schwinger multichannel (SMC) method for the lower energies (below 15 eV) and the independent atom model-based screening-corrected additivity rule with interferences (IAM-SCARI) for higher energies (above 15 eV). The latter was also applied to calculate the total ionization cross-sections, which were complemented with experimental values of the induced cationic fragmentation by electron impact. Double differential ionization cross-sections were measured with a reaction microscope multi-particle coincidence spectrometer. Using a momentum imaging spectrometer, direct measurements of the anion fragment yields and kinetic energies by the dissociative electron attachment are also presented. Cross-sections for the other inelastic channels were derived with a self-consistent procedure by sampling their values at a given energy to ensure that the sum of the cross-sections of all the scattering processes available at that energy coincides with the corresponding TCS. This cross-section data set is ready to be used for modelling electron-induced radiation damage at the molecular level to biologically relevant media containing 1M5NI as a potential radiosensitizer. Nonetheless, a proper evaluation of its radiosensitizing effects would require further radiobiological experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1989

Bandwidth efficient coding for fading channels - Code construction and performance analysis

The authors apply a general method of bounding the event error probability of trellis-coded modulation schemes to fading channels and use the effective length and the minimum-squared-product distance to replace the minimum-free-squared-Euclidean distance as code design parameters for Rayleigh and Rician fading channels with a substantial multipath component. They present 8-PSK trellis codes specifically constructed for fading channels that outperform equivalent codes designed for the additive white Gaussian noise channel when v is greater than or equal to 5. For quasiregular trellis codes there exists an efficient algorithm for evaluating event error probability, and numerical results on Pe which demonstrate the importance of the effective length as a code design parameter for fading channels with or without side information have been obtained. This is consistent with the case for binary signaling, where the Hamming distance remains the best code design parameter for fading channels. The authors show that the use of Reed-Solomon block codes with expanded signal sets becomes interesting only for large values of E(s)/N(0), where they begin to outperform trellis codes.

Schlegel, Christian↗
NASA NTRS2000

Autonomous Performance Monitoring System: Monitoring and Self-Tuning (MAST)

Maintaining the long-term performance of software onboard a spacecraft can be a major factor in the cost of operations. In particular, the task of controlling and maintaining a future mission of distributed spacecraft will undoubtedly pose a great challenge, since the complexity of multiple spacecraft flying in formation grows rapidly as the number of spacecraft in the formation increases. Eventually, new approaches will be required in developing viable control systems that can handle the complexity of the data and that are flexible, reliable and efficient. In this paper we propose a methodology that aims to maintain the accuracy of flight software, while reducing the computational complexity of software tuning tasks. The proposed Monitoring and Self-Tuning (MAST) method consists of two parts: a flight software monitoring algorithm and a tuning algorithm. The dependency on the software being monitored is mostly contained in the monitoring process, while the tuning process is a generic algorithm independent of the detailed knowledge on the software. This architecture will enable MAST to be applicable to different onboard software controlling various dynamics of the spacecraft, such as attitude self-calibration, and formation control. An advantage of MAST over conventional techniques such as filter or batch least square is that the tuning algorithm uses machine learning approach to handle uncertainty in the problem domain, resulting in reducing over all computational complexity. The underlying concept of this technique is a reinforcement learning scheme based on cumulative probability generated by the historical performance of the system. The success of MAST will depend heavily on the reinforcement scheme used in the tuning algorithm, which guarantees the tuning solutions exist.

Peterson, Chariya↗
NASA NTRS2001

X-Ray Probes of Cosmic Star Formation History

We discuss the imprints left by a cosmological evolution of the star formation rate (SFR) on the evolution of X-ray luminosities Lx of normal galaxies, using the scheme earlier proposed by us, wherein the evolution of LX of a galaxy is driven by the evolution of its X-ray binary population. As indicated in our earlier work, the profile of Lx with redshift can both serve as a diagnostic probe of the SFR profile and constrain evolutionary models for X-ray binaries. We report here the first calculation of the expected evolution of X-ray luminosities of galaxies, updating our work by using a suite of more recently developed SFR profiles that span the currently plausible range. The first Chandra deep imaging results on Lx evolution are beginning to probe the SFR profile of bright spiral galaxies; the early results are consistent with predictions based on current SFR models. Using these new SFR profiles, the resolution of the "birthrate problem" of low-mass X-ray binaries and recycled, millisecond pulsars in terms of an evolving global SFR is more complete. We discuss the possible impact of the variations in the SFR profile of individual galaxies and galaxy types.

Ghosh, Pranab↗
NASA NTRS2017

Status and Plans for Reanalysis at NASA/GMAO

Reanalysis plays a critical role in GMAOs goal to enhance NASA's program of Earth observations, providing vital data sets for climate research and the development of future missions. As the breadth of NASAs observations expands to include multiple components of the Earth system, so does the need to assimilate observations from currently uncoupled components of the system in a more physically consistent manner. GMAOs most recent reanalysis of the satellite era, MERRA-2, has completed the period 1980-present, and is now running as a continuing global climate analysis with two- to three-week latency. MERRA-2 assimilates meteorological and aerosol observations as a weakly coupled assimilation system as a first step toward GMAOs longer term goal of developing an integrated Earth system analysis (IESA) capability that will couple assimilation systems for the atmosphere, ocean, land and chemistry. The GMAO strategy is to progress incrementally toward an IESA through an evolving combination of coupled systems and offline component reanalyses driven by, for example, MERRA-2 atmospheric forcing. Most recently, the GMAO has implemented a weakly coupled assimilation scheme for analyzing ocean skin temperature within the existing atmospheric analysis. The scheme uses background fields from a near-surface ocean diurnal layer model to assimilate surface-sensitive radiances plus in-situ observations along with all other observations in the atmospheric assimilation system. In addition, MERRA-2-driven simulations of the ocean (plus sea ice) and atmospheric chemistry (for the EOS period) are currently underway, as is the development of a coupled atmosphere-ocean assimilation system. This talk will describe the status of these ongoing efforts and the planned steps toward an IESA capability for climate research.

MERRA-2↗
DOE OSTI2024

Optimization of ray-tracing simulations to confirm performance of the GP-SANS instrument at the High-Flux Isotope Reactor

The CG-2 beamline at the High Flux Isotope Reactor (HFIR) exhibits a notable discrepancy between observed count rates and the count rates we would expect based on a Monte-Carlo neutron ray-trace simulation. These simulations consistently predict count rates approximately five times greater than those observed in four separate experimental runs involving different instrument configurations. This discrepancy suggests that certain factors are causing losses in measurements that are not adequately accounted for in the simulation, in particular guide reflectivity or misalignment. To investigate these discrepancies, a high-dimensional simulation parameter approach is applied in order to understand the losses. Region of Interest (ROI) groups along the instrument are assigned to different surfaces of the guide components within the simulation. This allows the parameters of those guide components to be varied as a group to minimize the complexity of the search space. The result is an optimization of simulation parameters using an iterative scheme that aims to minimize the difference between experimentally measured count rates and simulated count rates across all tested collimator combinations. This proposed methodology holds the potential to reveal previously unrecognized sources of intensity loss in the CG-2 beamline at HFIR and improve the accuracy of simulations, leading to enhanced understanding and performance of the beamline for various scientific applications.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1995

Evaluation of Turbulence Models for Unsteady Flows of an Oscillating Airfoil

Unsteady flowfields of a two-dimensional oscillating airfoil are calculated using an implicit, finite-difference, Navier Stokes numerical scheme. Five widely used turbulence models are used with the numerical scheme to assess the accuracy and suitability of the models for simulating the retreating blade stall of helicopter rotor in forward flight. Three unsteady flow conditions corresponding to an essentially attached flow, light-stall, and deep-stall cases of an oscillating NACA 0015 wing experiment were chosen as test cases for computations. Results of unsteady airloads hysteresis curves, harmonics of unsteady pressures, and instantaneous flowfield patterns are presented. Some effects of grid density, time-step size, and numerical dissipation on the unsteady solutions relevant to the evaluation of turbulence models are examined. Comparison of unsteady airloads with experimental data show that all models tested are deficient in some sense and no single model predicts airloads consistently and in agreement with experiment for the three flow regimes. The chief findings are that the simple algebraic model based on the renormalization group theory (RNG) offers some improvement over the Baldwin Lomax model in all flow regimes with nearly same computational cost. The one-equation models provide significant improvement over the algebraic and the half-equation models but have their own limitations. The Baldwin-Barth model overpredicts separation and underpredicts reattachment. In contrast, the Spalart-Allmaras model underpredicts separation and overpredicts reattachment.

Srinivasan, G. R.↗
NASA NTRS2004

Interaction and Synergism of Microbial Fuel Cell Bacteria within Methanogenesis

Biological hydrogen production from waste biomass has both terrestrial and Martian advanced life support applications. On earth, biological hydrogen production is being explored as a greenhouse neutral form of clean and efficient energy. In a permanently enclosed space habitat, carbon loop closure is required to reduce mission costs. Plants are grown to revitalize oxygen supply and are consumed by habitat inhabitants. Unharvested portions must then be recycled for reuse in the habitat. Several biological degradation techniques exist, but one process, biophotolysis, can be used to produce hydrogen from inedible plant biomass. This process is two-stage, with one stage using dark fermentation to convert plant wastes into organic acids. The second stage, photofermentation, uses photoheterotrophic purple non-sulfur bacteria with the addition of light to turn the organic acids into hydrogen and carbon dioxide. Such a system can prove useful as a co-generation scheme, providing some of the energy needed to power a larger primary carbon recovery system, such as composting. Since butyrate is expected as one of the major inputs into photofermentation, a characterization study was conducted with the bacterium Rhodobacter sphaeroides SCJ, a novel photoheterotrophic non-sulfur purple bacteria, to examine hydrogen production performance at 10 mM-100 mM butyrate concentrations. As butyrate levels increased, hydrogen production increased up to 25 mM, and then decreased and ceased by 100 mM. Additionally, lag phase increased with butyrate concentration, possibly indicating some product inhibition. Maximal substrate conversion efficiency was 8.0%; maximal light efficiency was 0.89%; and maximal hydrogen production rate was 7.7 Umol/mg/cdw/hr (173 ul/mg cdw/hr). These values were either consistent or lower than expected from literature.

Klaus, David↗
NASA NTRS2011

Three-Dimensional Modeling of the Solar Wind: From the Coronal Base to the Outer Heliosphere

We have developed a global fully three-dimensional magnetohydrodynamic solar wind model for the region that extends from the coronal base to 100 AU. The simulation domain consists of tree spherical shell subdomains with computational boundaries between them placed at 20 solar radii and 0.3 AU. The location of the first boundary ensures that the flow at the boundary is both supersonic and super-Alfvenic. A steady-state solution in the innermost (coronal) region is obtained by the time-relaxation method. The solution uses a tilted dipole model or solar magnetograms as the boundary condition at the coronal base and includes a flux of Alfven waves in the WKB approximation which provide additional acceleration for the coronal outflow in the open field regions. The intermediate region solution is constructed by the integration of steady-state equations along radius using a marching scheme. The outer region solution (0.3-100 AU) is obtained again by the time relaxation and takes into account turbulence transport and heating as well as heating, flow deceleration, and other effects due to the interstellar pickup protons treated as a separate fluid. We use the model to simulate the global steady-state structure of the solar wind from the coronal base to the heliospheric boundary and compare the results with Ulysses and Voyager observations.

Usmanov, A. V↗
NASA NTRS1993

Navier-Stokes analysis of turbine flowfield and external heat transfer

An explicit 2D Navier-Stokes code has been modified and used to analyze the aerodynamics and heat transfer of a transonic turbine cascade. This code is based on a four-stage Runge-Kutta scheme. An algebraic Reynolds stress model (ARSM) and two versions of low Reynolds number (LRN) two-equation turbulence models, Chien's (1982) LRN k-epsilon model and Coakley's (1983) LRN q-omega model, have been employed in the computations. The surface pressure distributions and wake profiles are predicted well by all the models. The k-epsilon model and the k-epsilon/ARSM model yield better predictions of heat transfer than the q-omega model. The k-epsilon/ARSM solutions show some significant, though not dramatic, differences in the predicted skin friction coefficients, heat transfer rates, and performance parameters, as compared to the k-epsilon model. The predicted semiwake width is consistent with the measurement and correlation.

Luo, J.↗
NASA NTRS2002

Molten-Metal Droplet Deposition on a Moving Substrate in Microgravity: Aiding the Development of Novel Technologies for Microelectronic Assembly

Driven by advancements in microelectronics manufacturing, this research investigates the oblique (non-axisymmetric) impact of liquid-metal droplets on flat substrates. The problem of interest is relevant to the development of the novel technology of on-demand dispension (printing) of microscopic solder deposits for the surface mounting of microelectronic devices. The technology, known as solder jetting, features on-demand deposition of miniature solder droplets (30 to 120 microns in diameter) in very fine, very accurate patterns using techniques analogous to those developed for the ink-jet printing industry. Despite its promise, severe limitations exist currently with regards to the throughput rates of the technology; some of these limitations are largely due to the lack of the capability for reliable prediction of solder bump positioning and shapes, especially under ballistic deposition conditions where the droplet impact phenomena are inherently three-dimensional. The study consists of a theoretical and an experimental component. The theoretical work uses a finite element formulation to simulate numerically the non-axisymmetric (3-D) fluid mechanics and heat transfer phenomena of a liquid solder droplet impacting at an angle alpha on a flat substrate. The work focuses on the pre-solidification regime. The modeling of the most challenging fluid mechanics part of the process has been completed successfully. It is based upon the full laminar Navier-Stokes equations employing a Lagrangian frame of reference. Due to the large droplet deformation, the surface (skin) as well as the volumetric mesh have to be regenerated during the calculations in order to maintain the high accuracy of the numerical scheme. The pressure and velocity fields are then interpolated on the newly created mesh. The numerical predictions are being tested against experiments, for cases where wetting phenomena are not important. For the impact parameters used in the example shown (We = 2.38, Fr = 16300, Re = 157), the droplet rolls along the substrate, but its shape remains practically axisymmetric for all impact angles within the range from 0 to 60 deg. Interestingly, the substrate/droplet contact area during the recoiling phase of the impact is not a monotonically decreasing function of time. The experimental component of the research tests the numerical predictions and provides necessary input data (contact angles) for the theoretical model. The experiments are performed in microgravity (2.2s drop tower of the NASA GRC) in order to allow for the use of mm-size solder droplets, which make feasible the performance of accurate measurements, while maintaining similitude of the relevant fluid dynamic groups (Re, Fr, We, Ste). Preliminary oblique impact experiments have been performed using water droplets in normal gravity.

Megaridis, C. M.↗
DOE OSTI2020

A Systematic Approach to Light Load Calculation for DG Un-Intentional Islanding

With the influx of distributed generation (DG) penetration, power utilities have adopted active anti-islanding protection methods to prevent an unintentional islanding condition. One commonly used protection scheme is implemented through the use of a transfer trip from the upline connected substation to the DG facility and is required when the capacity of connected DGs exceed a certain threshold of local minimum load. The threshold established by standard IEEE 1547.2-2009 states that a 3:1 minimum load to generation ratio is acceptable to ensure the DG will not sustain an unintentional island. However, the term “minimum load” is not precisely defined in any standard, and in practice, a light load value is chosen by engineers through manual methods that can be inconsistent, not representative and change drastically based on various cases. To solve this problem, this paper proposes a new light load calculation method creating consistency and better accuracy; the method has been implemented in. NET C# application and is used within Dominion Energy.

Chen, Le↗
DOE OSTI2023

Pressure Transient Analyses and Poroelastic Modeling of Hydraulic Fracture Dilation for Multiple Injections at the Devine Fracture Pilot Site

Our team has conducted electromagnetic (EM) surveys for the past six years to monitor hydraulic-fracture behavior at the Devine Fracture Pilot Site (DFPS). The sub-horizontal orientation of a shallow hydraulic fracture at the DFPS provides uniform access to the fracture area for interrogation and data collection. Ahmadian et al. (2023) suggested a possible correlation between spatiotemporal changes in the flow rate, bottomhole pressure (BHP), and the observed surface recorded electric field at the DFPS. In this paper, we present the development of poroelastic forward models and pressure transient analyses (PTAs) to support the development of a multiphysics inverse model for these EM surveys. First, we conducted PTAs of the shut-in periods after six injections out of 10 to determine the fracture closure pressure (FCP) or the overburden pressure used in a poroelastic fracture reopening model. Second, we developed a finite-element poroelastic model throughout five injection cycles to include the effect of the cumulative injected volumes due to the previous injections on current fracture dilation in the presence of highly permeable unpropped and propped zones adjacent to the cohesive layer that models fracture reopening. Fracture reopening in this poroelastic model is based on a calibrated traction-separation response using the bottomhole pressure collected in two injection campaigns in 2020 and 2022. We used the outcomes of a previous simulation study of the primary hydraulic-fracturing stimulation to define the dimension of an unpropped fracture zone ahead of the propped fracture area. The PTAs led to FCPs consistent with those obtained using the injection data collected at the DFPS in 2020. Further, these analyses showed that at later injections, the fracture closure occurred at a later time with respect to the shut-in time, inferring the effect of cumulative injected volumes in previous injections. The simulation results show that considering the propped and unpropped fracture zones improves our poroelastic model in predicting the injection-well BHP. The numerical simulation results demonstrate a significant excess pore pressure near the fracture because of the preceding formation loadings by the previous injections. The obtained fracture dilation area and fluid pressure distribution provide a basis to improve the development of a multiphysics inverse model. Furthermore, in an iteratively coupled scheme, this pressure distribution can be introduced into EM models to render a holistic view of the causative mechanisms for the surface signal anomalies.

02 PETROLEUM↗