Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Displacement”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Deep Learning–Assisted Multiobjective Optimization of Geological CO 2 Storage Performance under Geomechanical Risks

In geological CO 2 storage, designing the optimal well control strategy for CO 2 injection to maximize CO 2 storage while minimizing the associated geomechanical risks is not trivial. This challenge arises due to pressure buildup, CO 2 plume migration, the highly nonlinear nature of geomechanical responses to rock-fluid interaction, and the high computational cost associated with coupled flow and geomechanics simulations. In this paper, we introduce a novel optimization framework to address these challenges. The optimization problem is formulated as follows: maximize total CO 2 storage while minimizing geomechanical risks by adjusting the injection schedules within bounded constraints. The geomechanical risks are primarily driven by injection-induced pressure build-up, which is characterized by ground displacement and the induced microseismicity. We used the Fourier neural operator (FNO)-based deep learning model to construct surrogate models, replacing the time-consuming coupled flow and geomechanics simulations for evaluating the aforementioned objective functions. The developed surrogate models have been incorporated into a multiobjective optimization framework through a genetic algorithm to reduce the computational burden. The proposed optimization framework reduces the computational cost from approximately 2,400 hours, when using objective function evaluations based on physics-based simulations, to around 20 minutes. A set of Pareto-optimal solutions of the proposed workflow yields nontrivial optimal decisions, reducing the microseismicity potential and the vertical displacement. This Pareto front highlights the optimal trade-offs between CO 2 storage amount, safety, and ground displacement, emphasizing the need for careful optimization and management of injection strategies to achieve a balanced outcome. The novelty of this work is twofold. First, we demonstrate the importance of incorporating the minimization of the geomechanical risks as objective functions into the CO 2 storage optimization workflow to mitigate the potential risk of induced microseismicity and ground displacement. Second, we leverage the FNO-based surrogate models to optimize a real-field CO 2 storage operation.

42 ENGINEERING↗
DOE OSTI2025

Analysis of temperature effects on self-powered neutron detectors in a gamma field

Self-powered neutron detectors (SPND) are small, robust, in-core radiation sensors that have the potential to operate in high temperature environments (above 600°C). However, previous experiments indicated temperature related effects in the form of displacement currents. Given the recent irradiation at The Ohio State University Research Reactor using a delayed-response rhodium-based SPND alongside a prompt-response tantalum-based SPND, it was concluded that only low-energy electrons from gamma interactions were the cause of these temperature effects. A series of heated gamma irradiations with temperatures up to 750°C of SPNDs were performed at the Advanced Test Reactor Gamma Facility in order to better characterize these displacement currents. The results from the experiments demonstrated a strong relation between the displacement currents that occurred during temperature ramps to thermoluminescence (TL) glow peaks—the SPND with Al2O3 insulation has one glow peak within 750°C and the SPND with MgO has three. A demonstration was performed to calculate the trap energy levels of the displacement currents through TL glow peak analysis. Unfortunately, the methods employed provided inconsistent results, requiring an updated model based on ongoing irradiation instead of post-irradiation measurements.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
DOE OSTI2020

Molecular-Scale Considerations of Enhanced Oil Recovery in Shale

With only less than 10% recovery, the primary production of hydrocarbon from shale reservoirs has redefined the energy equation in the world. Similar to conventional reservoirs, Enhanced Oil Recovery (EOR) techniques could be devised to enhance the current recovery factors. However, shale reservoirs possess unique characteristics that significantly affect the fluid properties. Therefore, we are adopting a molecular simulation approach that is well-suited to account for these effects to evaluate the performance of three different gases, methane, carbon dioxide and nitrogen, to recover the hydrocarbons from rough pore surfaces. Our hydrocarbon systems consists of either a single component (decane) or more than one component (decane and pentane). We simulated cases where concurrent and countercurrent displacement is studied. For concurrent displacement (injected fluids displace hydrocarbons towards the production region), we found that nitrogen and methane yielded similar recovery; however nitrogen exhibited a faster breakthrough. On the other hand, carbon dioxide was more effective in extracting the hydrocarbons when sufficient pressure was maintained. For countercurrent displacement (gases are injected and hydrocarbons are produced from the same direction), methane was found to be more effective, followed by carbon dioxide and nitrogen. In all cases, confinement reduced the recovery factor of all gases. This work provides insights to devise strategies to improve the current recovery factors observed in shale reservoirs.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2020

AGC-4 Experiment Irradiation Monitoring Data Qualification Final Report

The Graphite Technology Development Program ran a series of six experiments to quantify the effects of irradiation on nuclear grade graphite. This report focuses on the fourth experiment, Advanced Graphite Creep 4 (AGC 4). The Advanced Reactor Development (ARD) Technology Development Office (TDO) Program for research and development activities require documentation of qualified monitoring data to design and license the first high-temperature reactor nuclear plant. Qualified data meets the requirements for use as described in the experiment planning and quality assurance documents. Failed data do not meet the requirements and provide no useable information. Trend data may not meet all requirements but still provide some useable information. Use of Trend data requires assessment of how any deficiencies affect a particular use of the data. AGC-4 began with Advanced Test Reactor (ATR) Cycle 157D on May 30, 2015. After irradiating the graphite for two cycles, the capsule was removed from the reactor after ATR Cycle 158A, which ended on January 2, 2016, due to interference with another experiment. Irradiation was resumed with Cycle 162A on October 7, 2017 after the interfering experiment was removed from the reactor. AGC-4 irradiation completed by the end of Cycle 166B on January 10, 2020. Between Cycles 162A and 166B, AGC-4 capsule was removed from the reactor core during two Powered Axial Locator Mechanism (PALM) Cycles 163A and 165A to avoid overheating. All thermocouples (TCs) have functioned throughout the AGC-4 experiment. A total of 9,256,184 out of all 11,167,213 temperature records (or 82.9%) are Qualified for use by the ARD Program and 1,911,029 missing values are Failed. Argon, helium, and total gas flow data were within expected ranges, except only eight out-of-range values occurred during outages. A total of 21,240,627 out of all 22,288,5103 gas flow rates (or 92.8%) are Qualified for use by the ARD Program and 1,644,476 are Failed records mostly due to missing values. Discharge gas line moisture values were consistently low during full ATR power, except for 162B, when moisture content was increasing to more than 200 ppmv by the end of this cycle. During outages, moisture content reached as high as 700 ppmv. 1,089,941 out of a total of 1,107,331 moisture values (or 98.4%) are within the measurement range of the instrument and are Qualified for use by the ARD Program and 17,390 missing moisture values are Failed. Graphite creep specimens were subjected to one of three loads: 393, 491, or 589 lbf. For a brief period during Cycle 157D between 12:19 on June 2, 2015, and 08:23 on June 11, 2015, the load cells were wired incorrectly, resulting in missing stack load data. Missing stack loads were estimated from measured ram pressures using regression equations developed from the existing data from Cycle 157D. Estimated stack loads during this period are considered to be an accurate representation of actual load applied to the stacks. These loads deviate slightly from the planned loads. This deviation does not prevent the data from being Qualified for use but must be taken into account when analyzing the effect of load on creep. 6,095,728 out of a total of 6,403,985 moisture values (or 95.2%) are within the measurement range of the instrument and are Qualified for use by the ARD Program and 308,257 missing load values are Failed. Stack displacement increased consistently throughout the eight cycles, with total displacement reached highest value of 2.4 in by the end of irradiation. During ATR outages, a set of pneumatic rams raised the stacks of graphite creep specimens to ensure the specimens were not stuck within the test train. This stack raising was performed seven times throughout irradiation. All stacks were raised successfully each time. 4,744,974 out of a total of 6,094,513 displacement values (or 77.9%) are within the measurement range of the instrument and are Qualified for use by the ARD Program and 1,349,539 displacement values are Failed mostly due to missing values. Analyses were conducted on correlations between TCs to look for trends and step changes that might indicate instrument degradation or failure. Correlation analysis was used to identify instances when TCs form short circuits, referred to as virtual junctions, which result in TCs reporting temperatures from some location in the capsule other than the location where they were intended to read. No evidence of virtual junctions was found. Analyses were also conducted on control charts of temperature differences of between two TCs, which are expected to behave consistently throughout the entire irradiation period. Upward or downward trend over time indicates at least one TC in the pair was drifted. Examining control charts for all 66 possible pairs out of twelve TCs installed in the AGC-4 capsules reveal no clear drift failures occurred, except unstable behavior of two TCs in Zone 3, TC-7 and TC-8, over irradiation time. In conclus

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Application of Dynamical System Scaling to Bubble Dynamics

Dynamical System Scaling (DSS) provides a useful method for analyzing, categorizing and scaling time-dependent processes. A key feature of DSS is the temporal displacement rate, D, which relates the natural process time to the reference clock time. Because of its property of being invariant under a two-parameter affine transformation, it provides an underlying basis for process scaling. Application of DSS to bubble dynamics serves to elucidate the physical significance of the temporal displacement rate and its role in scaling a variety of bubble dynamics processes. This paper shows that the temporal displacement rate consists of the sum of dimensionless groups that govern the bubble dynamics. Preserving the temporal displacement rate for a prototype and a scaled model results in similitude of the time-dependent normalized bubble size, growth rates and interfacial accelerations for different fluid conditions. For vapor bubble growth in a superheated liquid, it is shown that inertial controlled bubble growth occurs when D=0 and thermally controlled bubble growth occurs when D=1 .

dynamical system scaling, SMR, bubble dynamics, DS↗
DOE OSTI2023

Machine-Learning-aided Approach for Predicting the Thermal Expansion Behaviors in Advanced Test Reactor Capsules (NURETH-20 full paper)

Instrumented experiments at test reactors are essential to deploying new advanced reactor systems. Designing new experiments and generating data on specific conditions require both time and cost investment. A high-fidelity model of the experiment environment can be created using finite element analysis software to support the actual experiments, but computation time is still a concern in applying outcomes to real-time usage (e.g., a digital twin). This research proposes a machine-learning-aided approach to temperature and displacement predictions, based on the thickness of the outer gas gap on the experimental capsule used for the in-pile demonstration of a novel thermal conductivity probe in the Advanced Test Reactor. The capsule consisted of U10Zr fuel, a rodlet, sodium, and inner and outer capsules. There were gas gaps between the fuel and rodlet and between the inner and outer capsule. The learning data consisted of an experimental capsule’s radial distributions of temperature and displacement, as obtained from Abaqus and the physical features. For the first step, temperature was predicted using three positional parameters. Then the displacement was predicted using six different positional parameters. Each physical feature was normalized to be both nondimensional and standardized. The temperature and displacement predictions showed good agreement in all cases involving interpolation and extrapolation. Also, data similarity enhancement increased the similarity between training and target data increasing the predictive accuracy of machine-learning models. In some cases of extrapolation, the accuracy of the machine-learning model showed limited performance, but still data similarity enhancement improved the accuracy.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Role of electronic energy loss on defect production and interface stability: Comparison between ceramic materials and high-entropy alloys

High-entropy alloys (HEAs) and some complex alloys exhibit desirable properties and significant structural stability in harsh environments, including possible applications in advanced reactors. Energetic ion irradiation is often used as a surrogate for neutron irradiation; however, the impact of ion electronic energy deposition and dissipation is often neglected. Moreover, differences in recoil energy spectrum and density of cascade events on damage evolution must also be considered. In many chemically complex alloys, the mean free path of electrons is reduced significantly, thus their decreased thermal conductivity and slow dissipation of localized radiation energy can have noticeable effects on displacement cascade evolution that is greatly different from metals with high thermal conductivity. In this work, nanocrystalline HEAs of Ni 20 Fe 20 Co 20 Cr 20 Cu 20 and nonequiatomic (NiFeCoCr) 97 Cu 3 , both having much lower room-temperature thermal conductivity than pure Ni or Fe, are chosen as model HEAs to reveal the role that electronic energy loss during ion irradiation has in complex alloys. The response of nanocrystalline HEAs is investigated under irradiation at room temperature using MeV Ni and Au ions that have different ratios of electronic energy to damage energy, which is the energy dissipated in displacing atoms. Different from previously reported amorphization of nanocrystalline SiC, experimental results on these HEAs show that, similar to the process in nanocrystalline oxide materials, both inelastic thermal spikes via electron–phonon coupling and elastic thermal spikes via collisions among atomic nuclei contribute to the overall grain growth. The growth follows a power law dependence with the total deposited ion energy, and the derived value of the power-exponent suggests that the irradiation-induced instability at and near grain boundaries leads to local rapid atomic rearrangements and consequently grain growth. The high power-exponent value can be attributed to the sluggish diffusion and delayed defect evolution arising from the chemical complexity intrinsic to HEAs. Here, this work calls attention to quantified fundamental understanding of radiation damage processes beyond that of simplified displacement events, especially in simulating neutron environments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Simultaneous 10 kHz three-dimensional CH 2 O and tomographic PIV measurements in a lifted partially-premixed jet flame

High-speed, three-dimensional (3D) scalar-velocity field measurements are demonstrated in a lifted partially-premixed dimethyl-ether/air jet flame using simultaneous laser-induced fluorescence (LIF) of formaldehyde (CH 2 O) and tomographic particle image velocimetry (TPIV). The 3D LIF measurements are performed by raster scanning the laser beam from a 100 kHz pulse-burst laser across the probe volume using an acousto-optic deflector. The volumetric reconstruction of the LIF signal from ten parallel planes provides quasi-instantaneous 3D LIF measurements that are synchronized with 10 kHz TPIV measurements. Here, the temporally resolved CH 2 O-LIF and velocity field data are used to analyze Lagrangian particle trajectories and displacement speeds at the base of the lifted flame. The particle trajectories reveal flow structures that are difficult to observe in an Eulerian reference frame. Positive and negative displacement speeds are observed at the CH 2 O-LIF surfaces at the inner and outer regions of the jet flame with a maximum displacement speed of approximately eight times the laminar flame speed of a stoichiometric dimethyl-ether/air mixture. The observations in present study support formation of lean premixed flames near the flame base whose role in flame stabilization should be investigated in more detail.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Wettability and capillary effects: Dynamics of pinch-off in unconstricted straight capillary tubes

We study the interfacial evolution of immiscible two-phase flow within a capillary tube in the partial wetting regime using direct numerical simulation. We investigate the flow patterns resulting from the displacement of a more viscous fluid by a less viscous one under a wide range of wettability conditions. We find that beyond a wettability dependent critical capillary number, a uniform displacement by a less viscous fluid can transition into a growing finger that eventually breaks up into discrete blobs by a series of pinch-off events for both wetting and non-wetting contact angles. This study validates previous experimental observations of pinch-off for wetting contact angles and extends those to non-wetting contact angles. We find that the blob length increases with the capillary number. We observe that the time between consecutive pinch-off events decreases with the capillary number and is greater for more wetting conditions in the displaced phase. Here, we further show that the blob separation distance as a function of the difference between the inlet velocity and the contact line speed collapses into two monotonically decreasing curves for wetting and non-wetting contact angles. For the phase separation in the form of pinch-off, this work provides a quantitative study of the emerging length and time scales and their dependence on the wettability conditions, capillary effects, and viscous forces.

42 ENGINEERING↗
DOE OSTI2025

ARC-DPA 1-MeV Neutron Fluence Equivalence Model for GaAs

Here, the legacy data used in the production of the American Society for Testing and Materials (ASTM) E722 standard for gallium arsenide (GaAs) were fit with the athermal recombination corrected-displacements per atom (ARC-DPA) model for the development of a novel 1-MeV neutron equivalent fluence metric. Improving on previous work, the coefficients for the ARC-DPA model were optimized using a conjugate gradient method, and the uncertainties in the observed damage in different benchmark neutron fields were sampled to quantify their contribution to the damage metric. The optimized parameters of −0.474 ± 0.03 for b arcdpa and 0.016 ± 0.002 for c arcdpa provided good agreement and indicated that the solution is not sensitive to the corresponding fluence uncertainties. Compared to the current 1-MeV neutron fluence equivalence standard for displacement damage in GaAs, our results indicate that for a thermal reactor neutron environment, displacement damage has been underestimated by around 25% and that the saturation value for the damage efficiency had been higher than previously modeled.

1-MeV equivalent fluence↗
DOE OSTI2021

Sentinel-1 Input Data for PSInSAR Analysis

Files used to perform the Persistent Scatterer InSAR analysis with SARPROZ. The data is sourced from ESAs Sentinel-1 project and covers Brady Hot Springs and Desert Peak geothermal areas. The original titles are included for the Sentinel-1 data. The naming guide is included as a link in this submission. The data contains SAR (Radar) data from the Sentinel 1A satellite between July 2017 and December 2019. The data is necessary to replicate our results, and can be used for further PSInSAR, DInSAR and other modern interferometric analyses to determine line of sight, vertical and, possibly, East-West horizontal displacement. These displacement analysis show the soil deformation in time (including average displacement velocity), which can be used to indicate subsoil phenomena: temperature changes, subsidence, uplift, and inform other analyses like seismicity and porosity.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Efficient, Compact, and Smooth Variable Propulsion Motor (Final Report)

In this project, a new architecture of highly efficient hydraulic motor was developed for the propulsion of off-highway vehicles. The motor uses an adjustable linkage driving a cam to vary the displacement of the piston, resulting in a Variable Displacement Linkage Motor (VDLM). The motor uses low friction rolling element bearings to significantly reduce mechanical friction, especially in the demanding low-speed high-torque conditions experienced by off-highway vehicles. The VDLM has high torque capabilities for its size due to the radial piston packaging and use of a multi-lobe cam. A VDLM is very smooth due to the ability to tune the torque ripple through the design of the cam profile. The project was divided into three periods. During the first period, a dynamic model was constructed of the motor to predict the performance of the motor and the vehicle. During the second period, a single-cylinder learning prototype was designed, built, and tested to validate the models constructed in the first period. In the third period, a multi-cylinder prototype motor was optimized, designed, fabricated, and tested. The motor demonstrated excellent mechanical efficiency (above 92.5% across the range of displacements), but the experimentally measure volumetric efficiency was lower than expected due to higher leakage rates created by the poor tolerance control on the prototype. To validate the dynamic models developed in the first period and better understand design trade-offs. In the third period a multi-cylinder concept demonstration prototype will be designed, fabricated, and tested. The final prototype will be tested on a motor dynamometer and will be utilized in hardware-in-the-loop testing to demonstrate its efficiency and performance impacts on the overall drive train. The experimental results were used in a drive train simulation of a compact track loader operating through a drive cycle. Using the VDLM in a hydrostatic circuit yielded 17.1% reduction in fuel consumption and 36.5% reduction in a series hybrid transmission.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Examination of Haines Jump in Microfluidic Experiments via Evolution Graphs and Interface Tracking

This work examines a type of rapid pore-filling event in multiphase flow through permeable media that is better known as Haines Jump. While existing microfluidic experiments on Haines Jump mostly seek to maintain quasi-steady states through very low bulk flow rates over long periods of time, this work explores the combined use of a highly structured microscale transport network, high-speed fluorescent microscopy, displacement front segmentation algorithms, and a tracking algorithm to build evolution graphs that track displacement fronts as they evolve through high-speed video recording. The resulting evolution graph allows the segmentation of a high-speed recording in both space and time, potentially facilitating topology-cognitive computation on the transport network. Occurrences of Haines Jump are identified in the microfluidic displacement experiments and their significance in bulk flow rates is qualitatively analyzed. The bulk flow rate has little effect on the significance of Haines Jump during merging and splitting, but large bulk flow rates may obscure small bursts at the narrowest part of the throat.

permeable media↗
DOE OSTI2021

In-situ observation of pre-, co- and post-seismic shear slip at 1.5 km depth

Understanding the initiation and arrest of earthquakes is one of the long-standing challenges of seismology. Here we report on direct observations of borehole displacement by a meter-sized shear rupture induced by pressurization of metamorphic rock at 1.5 km depth. We observed the acceleration of sliding, followed by fast co-seismic slip and transient afterslip. Total displacements were about 7, 5.5 and 9.5 micrometers, respectively for the observed pre-slip, co-seismic slip and afterslip. The observed pre-slip lasted about 0.4 seconds. Co-seismic slip was recorded by the 1 kHz displacement recording and a 12-component array of 3-C accelerometers sampled at 100 kHz. The observed afterslip is consistent with analytical models of arrest in a velocity-strengthening region and subsequent stress relaxation. < p >The observed slip vector agrees with the activation of a bedding plane within the phyllite, which is corroborated by relocated seismic events that were observed during the later stages of the injection experiment.

Schoenball, Martin↗
DOE OSTI2021

Floating platform effects on power generation in spar and semisubmersible wind turbines

The design and financing of commercial-scale floating offshore wind projects require a better understanding of how power generation differs between newer floating turbines and well-established fixed-bottom turbines. In floating turbines, platform mobility causes additional rotor motion that can change the time-averaged power generation. In this work, OpenFAST simulations examine the power generated by the National Renewable Energy Laboratory's 5-MW reference turbine mounted on the OC3-UMaine spar and OC4-DeepCWind semisubmersible floating platforms, subjected to extreme irregular waves and below-rated turbulent inflow wind from large-eddy simulations of a neutral atmospheric boundary layer. For these below-rated conditions, average power generation in floating turbines is most affected by two types of turbine displacements: an average rotor pitch angle that reduces power, caused by platform pitch; and rotor motion upwind-downwind that increases power, caused by platform surge and pitch. The relative balance between these two effects determines whether a floating platform causes power gains or losses compared to a fixed-bottom turbine; for example, the spar creates modest (3.1%–4.5%) power gains, whereas the semisubmersible creates insignificant (0.1%–0.2%) power gains for the simulated conditions. Furthermore, platform surge and pitch motions must be analyzed concurrently to fully capture power generation in floating turbines, which is not yet universal practice. Finally, a simple analytical model for predicting average power in floating turbines under below-rated wind speeds is proposed, incorporating effects from both the time-averaged pitch displacement and the dynamic upwind-downwind displacements.

17 WIND ENERGY↗
DOE OSTI2024

Embedded symmetric positive semi-definite machine-learned elements for reduced-order modeling in finite-element simulations with application to threaded fasteners

Here, we present a machine-learning strategy for finite element analysis of solid mechanics wherein we replace complex portions of a computational domain with a data-driven surrogate. In the proposed strategy, we decompose a computational domain into an “outer” coarse-scale domain that we resolve using a finite element method (FEM) and an “inner” fine-scale domain. We then develop a machine-learned (ML) model for the impact of the inner domain on the outer domain. In essence, for solid mechanics, our machine-learned surrogate performs static condensation of the inner domain degrees of freedom. This is achieved by learning the map from displacements on the inner-outer domain interface boundary to forces contributed by the inner domain to the outer domain on the same interface boundary. We consider two such mappings, one that directly maps from displacements to forces without constraints, and one that maps from displacements to forces by virtue of learning a symmetric positive semi-definite (SPSD) stiffness matrix. We demonstrate, in a simplified setting, that learning an SPSD stiffness matrix results in a coarse-scale problem that is well-posed with a unique solution. We present numerical experiments on several exemplars, ranging from finite deformations of a cube to finite deformations with contact of a fastener-bushing geometry. We demonstrate that enforcing an SPSD stiffness matrix drastically improves the robustness and accuracy of FEM–ML coupled simulations, and that the resulting methods can accurately characterize out-of-sample loading configurations with significant speedups over the standard FEM simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Numerical investigation of multiphase flow through self-affine rough fractures

Multiphase flow through fractures has great significance in subsurface energy recovery and gas storage applications. Different fracture and flow properties affect flow through a fracture which is difficult to control in laboratory experiments. Here, we perform lattice Boltzmann simulations in an ensemble of synthetically generated fractures. Drainage simulations are performed at different capillary numbers, wettability, and viscosity ratios. We track the invading front and quantify breakthrough saturations and show that roughness and wettability have a strong effect on fluid invasion through a complex fracture. Invading a more viscous fluid results in more stable displacement regardless of the capillary number while at very low capillary numbers, fluid migration is dependent on the inherent structure of the fracture. We develop a fluid displacement phase diagram in a single rough fracture and compare our results from that in the literature. Finally, we extend the phase diagrams across multiple fractures and demonstrate the importance of natural fracture features of roughness and wettability in identifying stable versus unstable displacement regimes during multiphase flow through rough fractures. Our work presents an end-to-end numerical pathway for testing on experimental data and expanding numerical data sets for testing combinations of different physical phenomenon and make valuable predictions on fluid flow through rough fractures.

02 PETROLEUM↗