SEARCH · Search NASA

Results for “MATHEMATICS - ITERATION”

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 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Enhancing occupant behavior representation for interoperability between building information modeling and building energy modeling

Building Performance Simulation (BPS) has been adopted as an essential tool for designing, operating, and retrofitting buildings to optimize energy efficiency throughout the building life cycle. The Green Building XML (gbXML) schema facilitates seamless data exchange between Building Information Modeling (BIM) and Building Energy Modeling (BEM) software tools. However, limited occupant behavior (OB) representation in BIM often leads to inconsistent and inaccurate energy simulation in BEM software. This paper presents 154 systematic enhancements to the existing occupant behavior XML (obXML) schema v1.3.4, initially developed for standardizing OB representation for BEM, to address existing limitations and improve interoperability with BIM models. The enhancements encompass improved integration with BIM models through extended building representations and system operations, expanded support for advanced OB models with additional environmental parameters and mathematical capabilities, and implementation of a standardized model documentation framework. To facilitate seamless data transformation between gbXML and obXML schemas, we developed a publicly available gb-obXML Schema Converter. Three case studies demonstrate the enhanced schema’s capabilities: representation of building information using a two-story office building model, documentation of a window operation behavior model, and validation of the schema converter’s functionality. The enhanced obXML schema v1.4 enables sophisticated modeling of occupant-building interactions while maintaining consistency with industry-standard BIM schemas. The standardized documentation framework facilitates reproducibility and knowledge sharing in the OB research community, while the schema converter automates the integration of building information into OB simulation workflows. These enhancements establish a foundation for more accurate building performance simulation by supporting sophisticated representation of occupant behavior within the BIM-to-BEM simulation workflows.

Chung, Jihoon↗

DOE OSTI2025

Comparative studies of cross-phase dynamics in turbulent momentum flux and particle flux at the tokamak edge

Turbulent transport events, including turbulent transport flux of momentum (i.e., turbulent momentum flux or Reynolds stress) and turbulent transport flux of particle (i.e., turbulent particle flux), have important effects on the confinement performance of magnetic confinement fusion devices. Poloidal Reynolds stress is the ensemble average of the product of radial velocity fluctuations and poloidal velocity fluctuations, i.e., $\langle {\widetilde{v}}_{r}{\widetilde{v}}_{\theta }\rangle$. Turbulent particle flux is the ensemble average of the product of radial velocity fluctuations and density fluctuations, i.e., $\langle \widetilde{n}{\widetilde{v}}_{r}\rangle$. Changes in either amplitude of fluctuations or cross phase between fluctuations can cause changes in turbulent transport. In this paper, cross-phase dynamics in the Reynolds stress and turbulent particle flux at the tokamak edge are studied in detail. Reynolds stress and turbulent particle flux are, respectively, written as the product of fluctuation amplitudes and an average cross-phase factor. The mathematical expressions of the average cross-phase factors are derived. The average cross-phase factors and the power spectra of cross phase are obtained using experimental measurement data. It is found that the cross-phase dynamics in Reynolds stress and particle flux are very different. Reynolds stress is found to be more sensitive to cross phase than particle flux is. In the strong $E\times B$ shear layer, spatial slips of cross phase lead to the obvious radial gradient of Reynolds stress. In the no/weak $E\times B$ shear region, the cross phase in Reynolds stress tends to lock. Here, phase locking refers to that the power spectra of phase tend to distribute around a fixed phase which does not change with radial position, while phase slip means that the power spectra of cross phase tend to distribute around a phase that varies with radial position. Phase slip or locking mainly describes the central phase weighted by the power spectra, while the phase scattering mainly describes the dispersion of the power spectrum distribution of the phase. The increased scattering of cross phase, which indicates the power spectra distribution of the phase is more dispersed, contributes to the decreased Reynolds stress for higher collisionality. The cross phase in particle flux tends to lock in both strong and no/weak shear regions. The degree of scattering of cross phase in the particle flux does not change obviously as collisionality increases. For higher collisionality, it is the increased density fluctuation amplitude rather than cross-phase dynamics that leads to the increased particle flux. The underlying physical mechanism that causes Reynolds stress and particle flux to exhibit different phase dynamics is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Generalization error guaranteed auto-encoder-based nonlinear model reduction for operator learning

Many physical processes in science and engineering are naturally represented by operators between infinite-dimensional function spaces. The problem of operator learning, in this context, seeks to extract these physical processes from empirical data, which is challenging due to the infinite or high dimensionality of data. An integral component in addressing this challenge is model reduction, which reduces both the data dimensionality and problem size. In this paper, we utilize low-dimensional nonlinear structures in model reduction by investigating Auto-Encoder-based Neural Network (AENet). AENet first learns the latent variables of the input data and then learns the transformation from these latent variables to corresponding output data. Our numerical experiments validate the ability of AENet to accurately learn the solution operator of nonlinear partial differential equations. Furthermore, we establish a mathematical and statistical estimation theory that analyzes the generalization error of AENet. Finally, our theoretical framework shows that the sample complexity of training AENet is intricately tied to the intrinsic dimension of the modeled process, while also demonstrating the robustness of AENet to noise.

Auto-encoder↗

DOE OSTI2025

Elucidating grain boundary energy minimization mechanisms in textured Ca-doped alumina with inclination-dependent Monte Carlo Potts simulations

The grain growth behavior of textured Ca-doped alumina is compared to Monte Carlo Potts (MCP) simulations to investigate the effect of anisotropic grain boundary (GB) properties on local boundary migration. Experimentally, the growth of textured Ca-doped alumina results in highly elongated grains. The relative GB energy distribution is measured using the thermal groove method before and after heat treating at 1600°C, finding that high energy GBs are eliminated during grain growth. No significant difference in the GB energy distributions is found between the long and short axes of the elongated grains, suggesting that anisotropic mobility may be responsible for the grain shape. However, MCP simulations with anisotropic mobility as a function of plane inclination do not result in grains with distinct morphologies, regardless of the degree of anisotropy introduced. The final grain shape after grain growth of textured Ca-doped alumina resembles that of the MCP simulations using an anisotropic GB energy as a cosine function of plane inclination. Several energy functions are tested and only those that mathematically impose a torque (second derivative of energy with respect to the plane inclination angle) result in elongated grains. Although area reduction is the dominant energy minimization mechanism, these results suggest that local GB migration is affected by anisotropic GB energy and torque and alternative mechanisms like GB replacement and reorientation.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Localized material compression to correct distortion in wire arc additive manufacturing

Wire Arc Additive Manufacturing (WAAM) is an advanced manufacturing technology which utilizes welding systems to generate three dimensional geometries in a layer-by-layer fashion. Distortion or warping of a print substrate and WAAM components due to thermally induced residual stresses is an ongoing challenge limiting the widespread adoption of WAAM technologies for producing components. In this manuscript, a novel approach is described to address thermal distortion in deposited components by applying lateral compressions along the length of the deposited material. To demonstrate this method, a series of single-track walls were printed and compressed at evenly spaced intervals using a modified hydraulic cutter tool. The jaws of the tool were modified to compress material rather than to shear it. A mathematical model was developed to relate the curvature of the deposited material to the volume of compression required to eliminate this distortion. Validation of this model was performed using 3D scan data to compare the change in wall curvature induced by compression to the volume of the applied compressions. Substrate deflection was also compared against a control wall, and implementation of wall compression reduced maximum deflections by 93% across a series of four depositions and subsequent compressions. Wall cross sections were also analyzed to determine the impact of compression on material hardness and grain structure. The results demonstrate that successively placed lateral compressions can effectively control and potentially eliminate bending distortion in printed parts. This methodology can be further developed to form a robust model for correction of thermally-induced distortion in WAAM components.

Additive manufacturing↗

DOE OSTI2024

A dynamic volumetric heat source model for laser additive manufacturing

Melt pool scale models of laser powder bed fusion (LPBF) offer insights into the process-structure-property relationships in additive manufacturing (AM). These models often neglect physical phenomena such as vapor cavity formation and fluid mechanics to reduce computational demands. Instead, volumetric heat source models are used to represent the effects that these phenomena have on the predicted melt pool dimensions. Generally, the dimensions and effective absorption of the volumetric heat source are calibrated to reproduce melt pool dimensions observed in metallographic cross sections taken from single-track experiments on bare plate. However, the transient nature of LPBF often deviates the melt pool dimensions from the assumed steady-state conditions of single-track experiments, motivating the need for a volumetric heat source model that more generally considers the dynamic relationship between melt pool shape and laser-material interactions. Here, we introduce a two-parameter volumetric heat source model that integrates several existing models into a generalized mathematical expression, providing independent control over the radial heat distribution via the parameter k and the volumetric shape of the heat source via the parameter m. This parameterization enables the calibration of melt pool shape predictions through simultaneous adjustment of these parameters, while keeping the radial heat source dimensions consistent with the experimental spot size (D4σ) and constraining the heat source depth and absorption to physically derived expressions for cavities. Consequently, the proposed volumetric heat source model adapts to changes in the local melt pool conditions due to scanning strategy and part geometry by dynamically adjusting the heat source depth and absorption. We demonstrate the capabilities of the proposed model through comparisons with a collection of experiments from the Additive Manufacturing Benchmark (AMBench).

36 MATERIALS SCIENCE↗

DOE OSTI2025

Catalytic hydrogenation of HMF to BHMF over copper catalysts

2,5-Bis(hydroxymethyl)furan (BHMF) is a bio-derived building block for polyester production, obtained via the hydrogenation of 5-hydroxymethylfurfural (HMF). First-principles thermodynamic equilibrium calculations indicate that this reaction is not thermodynamically limited under relevant conditions (e.g., 100 °C and high H 2 partial pressure). In this work, crude HMF was employed as the feedstock for BHMF synthesis. Initially, acidic impurities and humins were removed from unrefined HMF through filtration using a packed bed of γ-alumina. A comprehensive study of the filtration process is presented, including filtration kinetics, breakthrough curve analysis, and mathematical modeling. The purified HMF was subsequently hydrogenated over a 10 wt% CuZrO 2 catalyst, using ethanol as the reaction solvent. Batch reactions were first performed for collection of kinetic data to guide the transition to continuous flow operation. Kinetic data was collected in a fixed bed reactor at varying contact time, time on stream, temperature, and HMF concentration. This data was used to develop a kinetic model for HMF hydrogenation. Maximum BHMF production rates were achieved at 130 °C, accompanied by minor formation of byproducts from BHMF ring-opening reactions. The BHMF selectivity was 100 % at 100 °C although with lower reaction rates. Furthermore, catalyst stability tests revealed a loss of up to 50 % in catalytic activity within the first 24 h, likely due to the adsorption of HMF-derived oligomers that are not easily removed by filtration.

Crude HMF filtration↗

DOE OSTI2025

CFD modeling of turbulent air flow in self-heated gyroid TPMS structures: Thermal-hydraulic performance and validation

The application of mathematically derived geometries, such as triply periodic minimal surface (TPMS) lattices, has garnered significant interest across various fields, including the nuclear sector, due to their superior thermal-hydraulic characteristics for heat transfer compared to traditional plain or finned tubes. Here, this study validates a computational fluid dynamics (CFD) model, evaluates different turbulence models and CFD model settings, and performs uncertainty quantification to provide a comprehensive analysis. Despite extensive research on CFD modeling of TPMS lattices, such as gyroid and diamond geometries, there is a notable lack of publicly available literature providing comprehensive details on numerical analysis aspects, including convergence and methodological best practices. This study embarks on a benchmark analysis of a gyroid geometry to evaluate its thermal-hydraulic performance under turbulent flow conditions and scrutinize various CFD model configurations. The main contributions of this work include validating the CFD model, assessing and comparing different turbulence models, and enhancing pressure drop and temperature prediction capabilities. The results aim to support the development of methodologies needed to benchmark and enhance numerical analysis techniques for TPMS lattices. This work seeks to complement the existing body of knowledge, support the development of TPMS reactor concepts, and improve best practices for CFD modeling of TPMS lattices, ultimately advancing methodologies to support future applications in this domain.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Analysis of the weighted shifted boundary method for the Poisson and Stokes problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. Accuracy is maintained by properly shifting the location and values of the boundary conditions. This avoids integration over cut cells and the associated implementation issues. Recently, the Weighted SBM (WSBM) was proposed for the Navier-Stokes equations with free surfaces and the Stokes flow with moving boundaries. The attribute “weighted” in the name WSBM stems from the fact that its variational form is weighted with the elemental volume fraction of active fluid. The motivation for the development of the WSBM was the preservation of the volume of active fluid to a higher degree of accuracy, which in turn resulted in improved stability and robustness characteristics in moving-boundary, time-dependent simulations. In this article, we present the numerical analysis of the WSBM formulations for the Poisson and Stokes problems. We give mathematical conditions under which the bilinear forms defining the discrete variational formulations are uniformly coercive (Poisson problem) or inf-sup stable (Stokes problem). By these results, stability and optimal convergence is proven in the natural norm; L2-error estimates can also be derived.

Approximate domain boundaries↗

DOE OSTI2025

Degrees of rate control and AutoDiff-driven direct sensitivity analysis in heterogeneous catalysis

Despite the wide application and benefits of the degree of rate control (DRC) analysis, several details remain argued, particularly about the conservation of DRCs at transient (TR) and steady-state (SS) conditions, especially for complex reaction networks. This work argues that previous proofs about the conservation properties of DRCs have been incomplete, and we provide new mathematical proofs at TR and SS conditions. In addition, we use both analytical (automatic differentiation) and numerical (finite difference) approaches to compute DRCs for the case study of ethane hydrogenolysis (EH) over Pt(111). This work confirms that at both TR and SS conditions, the sum of all DRCs, i.e., sum of the degrees of kinetic (DKRC) and thermodynamic rate control (DTRC), is conserved at zero. At SS conditions, the sum of DKRC is conserved at 1 while the sum of DTRC is conserved at −1. In corroboration of previous works, we show that the DTRC for any adsorbate at SS is equal to the product of the species coverage and a constant. In contrast, at TR conditions, the individual sums of both DTRC and DKRC are not conserved and can be any real number, with potential implications for the novel field of dynamic catalysis. Finally, we show that the conventional finite difference (FD) approach, only useful at SS, is prone to inaccuracy and very sensitive to the value of the differential change applied. The optimal differential value also varies significantly with system and rate definition. Consequently, we describe and illustrate in this work the application of the automatic differentiation (AD) approach for the more accurate determination of DRCs at both TR and SS conditions.

Automatic differentiation↗

DOE OSTI2025

Implementation and evaluation of multi-dual mode counter-current chromatography in the CUP Modeler software

Counter-current chromatography (CCC) is a separation technique that utilizes immiscible solvent pairs as stationary and mobile phases, which imparts numerous benefits compared to solid-liquid chromatography including the ability to treat either the more-dense or less-dense solvent layer as the mobile phase. Multi-dual mode (MDM) is a CCC elution mode capable of improving the separation of closely eluting compounds by alternating upper- and lower-layer solvent flows in opposing directions within the same separation. While some effort has been made to model MDM, implementation of these models in experimental design has yet to be widely adopted. Accordingly, we further developed our previously published cell utilized partitioning (CUP) model to include MDM predictions with CCC and packaged the full suite of CUP modeling capabilities into a user-friendly, open-source tool called the CUP Modeler. The mathematical model for MDM CCC was derived and validated with experimental separation of ethyl guaiacol (EG) and ethyl phenol (EP), two compounds that co-elute in our previously demonstrated reductive catalytic fractionation (RCF) lignin monomer isolation method. The developed MDM model provided insights into the effect of multiple operating parameters - including stationary phase retention, flow rate, column efficiency, feed concentration ratio, selectivity factor, and solute distribution ratios - on the separation yields, productivity, and purities. Our model agreed with prevailing understanding of MDM but also revealed new insights including that the ideal distribution ratios for co-eluting solutes to be separated by MDM is between 1.1 and 1.5, with the lower value ideally close to 1.25. Overall, this work provides fundamental insights for MDM process design and enables broader adoption of general liquid-liquid chromatography with a new, open-source user-friendly interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Constrained or unconstrained? Neural-network-based equation discovery from data

Throughout many fields, practitioners often rely on differential equations to model systems. Yet, for many applications, the theoretical derivation of such equations and/or the accurate resolution of their solutions may be intractable. Instead, recently developed methods, including those based on parameter estimation, operator subset selection, and neural networks, allow for the data-driven discovery of both ordinary and partial differential equations (PDEs), on a spectrum of interpretability. The success of these strategies is often contingent upon the correct identification of representative equations from noisy observations of state variables and, as importantly and intertwined with that, the mathematical strategies utilized to enforce those equations. Specifically, the latter has been commonly addressed via unconstrained optimization strategies. Representing the PDE as a neural network, we propose to discover the PDE (or the associated operator) by solving a constrained optimization problem and using an intermediate state representation similar to a physics-informed neural network (PINN). The objective function of this constrained optimization problem promotes matching the data, while the constraints require that the discovered PDE is satisfied at a number of spatial collocation points. We present a penalty method and a widely used trust-region barrier method to solve this constrained optimization problem, and we compare these methods on numerical examples. Our results on several example problems demonstrate that the latter constrained method outperforms the penalty method, particularly for higher noise levels or fewer collocation points. This work motivates further exploration into using sophisticated constrained optimization methods in scientific machine learning, as opposed to their commonly used, penalty-method or unconstrained counterparts. For both of these methods, we solve these discovered neural network PDEs with classical methods, such as finite difference methods, as opposed to PINNs-type methods relying on automatic differentiation. Here, we briefly highlight how simultaneously fitting the data while discovering the PDE improves the robustness to noise and other small, yet crucial, implementation details.

Data-driven discovery↗

DOE OSTI2025

The latent variable proximal point algorithm for variational problems with inequality constraints

The latent variable proximal point (LVPP) algorithm is a framework for solving infinite-dimensional variational problems with pointwise inequality constraints. The algorithm is a saddle point reformulation of the Bregman proximal point algorithm. At the continuous level, the two formulations are equivalent, but the saddle point formulation is more amenable to discretization because it introduces a structure-preserving transformation between a latent function space and the feasible set. Working in this latent space is much more convenient for enforcing inequality constraints than the feasible set, as discretizations can employ general linear combinations of suitable basis functions, and nonlinear solvers can involve general additive updates. LVPP yields numerical methods with observed mesh-independence for obstacle problems, contact, fracture, plasticity, and others besides; in many cases, for the first time. The framework also extends to more complex constraints, providing means to enforce convexity in the Monge–Ampère equation and handling quasi-variational inequalities, where the underlying constraint depends implicitly on the unknown solution. Here, in this paper, we describe the LVPP algorithm in a general form and apply it to ten problems from across mathematics.

Inequality constraints↗

DOE OSTI2024

Deep operator network surrogate for phase-field modeling of metal grain growth during solidification

A deep operator network (DeepONet) has been constructed that generates accurate representations of phase-field model simulations for evolving two dimensional metal grain morphology growing from melt. These representations serve as lower resolution, computationally efficient stand-ins for quick parameter space exploration of solutions to the the Allen-Cahn equations that dictate the phase-field model simulations. The experimental target for the phase-field model is a uranium casting system cooling a 434 g uranium charge from a maximum temperature of 1400° C at an average rate of 30° C / min , traversing the crystallographic phases of the pure metal. Experimental parameters inform the phase-field model, whose higher resolution computational model solutions are used to train the DeepONet in a given parameter space with the aim of developing a faster, more efficient method for predicting the solidifying metal's microstructure at different potential experimental values. The final DeepONet generates high accuracy, lower resolution predictions with cumulative relative approximation error over all timesteps of less than 0.5%, while ensuring solutions remain within physically feasible ranges. Further, these relative error values are comparable with other state-of-the-art DeepONet models for microstructure evolution, while significantly reducing the amount of training data required. Training a convolutional neural network simultaneously with the DeepONet, enforcing realistic values at the complex metal grain boundaries, and mathematically encoding boundary conditions into the structure of the DeepONet improved prediction accuracy and computational efficiency over a standard DeepONet model.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Crop models: integrating systems from the molecular to global for agricultural productivity and sustainability

Mathematical models that simulate crop growth in response to environmental conditions and management practices are essential tools for exploring agriculture-based strategies to address food security and environmental sustainability challenges. Early applications of crop models focused on supporting farmers in making management decisions. Applications have since expanded to estimating future impacts on local and global food production from changing climates. Emerging applications of crop models aim to leverage how these models integrate plant processes across biological scales to identify engineering or breeding strategies that account for environmentally-responsive dynamics at field scales and for exploring solutions to improve sustainability. In this review, we highlight recent studies across these four broad application areas and highlight potential future directions for the crop modeling field.

Piao, Ximin [Univ. of Illinois at Urbana-Champaign↗

DOE OSTI2024

Stress intensity factor models using mechanics-guided decomposition and symbolic regression

The finite element method can be used to compute accurate stress intensity factors (SIFs) for cracks with complex geometries and boundary conditions. In contrast, handbook solutions act as surrogate SIF models that provide significantly faster evaluation times. However, the development of conventional surrogate SIF models relies on manual development based on low-order parameterizations. This limits surrogate model accuracy and generalizability. Here, in this paper, we develop a framework for the automated development of mechanics-guided handbook SIF solutions by using interpretable machine learning via genetic programming for symbolic regression (GPSR). Formalizing the mechanics-based approach of Raju and Newman, SIF training data is decomposed into multiple subsets. This decomposition enables parallel GPSR model development of subfunctions, each of which accounts for specific geometrical corrections with respect to a known analytical model. Using this mechanics-based approach with GPSR allows for equations to be learned with improved accuracy and reduced complexity relative to the Raju Newman equations while maintaining the inherent interpretability of mathematical expressions. In this paper, we present equations that match the complexity of the Raju Newman equations while having reduced error, as well as equations with similar errors and reduced complexity.

42 ENGINEERING↗

DOE OSTI2024

An OpenStreetMaps based tool to study the energy demand and emissions impact of electrification of medium and heavy-duty freight trucks

In this paper, we present the mathematical formulation of an OpenStreetMaps (OSM) based tool that compares the costs and emissions of long-haul medium and heavy-duty (M&HD) electric and diesel freight trucks, and determines the spatial distribution of added energy demand due to M&HD EVs. The optimization utilizes a combination of information on routes from OSM, utility rate design data across the United States, and freight volume data, to determine these values. In order to deal with the computational complexity of this problem, we formulate the problem as a convex optimization problem that is scalable to a large geographic area. In our analysis, we further evaluate various scenarios of utility rate design (energy charges) and EV penetration rate across different geographic regions and their impact on the operating cost and emissions of the freight trucks. Our approach determines the net emissions reduction benefits of freight electrification by considering the primary energy source in different regions. Such analysis will provide insights to policy makers in designing utility rates for electric vehicle supply equipment (EVSE) operators depending upon the specific geographic region and to electric utilities in deciding infrastructure upgrades based on the spatial distribution of the added energy demand of M&HD EVs. To showcase the results, a case study for the U.S. state of Texas is conducted.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2024

Relative permeabilities for two-phase flow through wellbore cement fractures

Multiple fluids are likely to exist in fractures and flow paths associated with leaky wellbores, including liquids (e.g., crude oil) and gases (e.g., gas exsolved from liquid). These fluids occupy and move through different portions of the pore spaces within the fractures depending on many factors, including fluid properties, fracture size, and the amount of the different fluids. Upward leakage of any phase, through the fracture, can contaminate water-bearing formations, create hazardous surface conditions, and compromise the functionality of the wellbore. Early signs of wellbore leaks may be expressed by anomalous pressure behavior at surface monitoring points on cavern storage wells. These pressure anomalies are difficult to interpret, necessitating knowledge of the factors that affect the multiphase flow in fractures and porous media. These parameters are critical to modeling multiphase flow in fractures. This insight can guide further diagnosis and maximize leak remediation. Here, our study focuses on the relationship of the liquid–gas relative permeabilities for representative variable-aperture wellbore cement fracture. To obtain the relative permeability of each phase, two-phase flow tests were conducted where both fluids were flowing simultaneously through a fractured wellbore cement specimen under a range of factors, namely (1) aperture size, (2) capillary numbers, and (3) viscosity ratio. The flow experiments were conducted under a range of confining stresses and flow velocities, using nitrogen gas and silicone oils (of different viscosities) in a specially designed pressure vessel. The sum of gas and oil relative permeabilities were found to be less than one under all conditions, which indicates that the presence of one phase affects the permeability of the other phase, and vice versa. Since the gas phase flow conditions include a significant inertial flow component in addition to viscous flow, the inertial flow coefficients at different saturation states are presented. The factors affecting the relationship between the relative permeabilities are discussed in detail. A new mathematical model for estimating the relative permeability of wellbore cement fracture is presented and experimentally validated.

58 GEOSCIENCES↗