Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parallel application”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

High-power ultraviolet (UV) and vacuum ultraviolet (VUV) lamps with micro-cavity plasma arrays

A plasma lamp includes plates that are approximately parallel, with at least one array of microcavities formed in a surface of at least one plate. When desirable, the plates are separated a fixed distance by spacers with at least one spacer being placed near the plate's edge to form a hermetic seal therewith. A gas makes contact with the microcavity array. Electrodes capable of delivering a time-varying voltage are located on the surface of each plate. At least one electrode is located on an exterior surface of at least one interior plate. Optionally, protective windows may be placed over the electrodes. The application of the time-varying voltage interacts with the gas to form a glow discharge plasma in the microcavities and the fixed volume between the plates (when present). The glow discharge plasma efficiently and uniformly emits UV/VUV radiation over the entire surface of the lamp.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Numerical modeling of chemical looping oxidative dehydrogenation of ethane in parallel packed beds

Chemical looping oxidative dehydrogenation (CL-ODH) of ethane has the potential to be a highly efficient alternative to steam cracking for ethylene production. Accurate reactor modeling is of critical importance to efficiently scale up and optimize this new technology. This study reports a one-dimensional, heterogeneous packed bed model to simulate the CL-ODH of ethane to ethylene with a Na 2 MoO 4 -promoted CaTi 0.1 Mn 0.9 O 3 redox catalyst. Here, the overall reaction kinetics was well-described by coupling the gas-phase steam cracking of ethane with the reduction kinetics of the redox catalyst by H 2 and C 2 H 4 . The impact of H 2 on the formation rate of CO 2 byproduct from C 2 H 4 conversion was also thoroughly investigated to validate the applicability of the kinetic model under operational environments. The temperature variation within the different CL-ODH steps and the temperature distribution along the bed were also carefully considered. The accuracy of the model was validated by experiments conducted in a large lab-scale packed bed reactor (200 g catalyst loading), with an average deviation of 2.8% in terms of ethane conversion and ethylene yield. The model was subsequently used to optimize the operating parameters of the CL-ODH reactor, indicating that up to 63.7% single-pass C 2 +olefin yield can be achieved with the current redox catalyst bed whereas further optimization of the redox catalyst to inhibit C 2 H 4 activation can result in 69.4% single-pass C 2 +yield while maintaining low CO 2 selectivity.

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

Monofluorination of Naphthyls Promotes the Cofacial π–π Stacking and Increases the Electron Mobility of Non-Planar Zinc(II) Complexes of Di(naphthylethynyl)azadipyrromethene

The homoleptic zinc(II) complex of [2,8-di(1-naphthylethynyl) 3,7-diphenyl 1,9-(4-hexylphenyl)azadipyrromethene (ZnL2) 2 ] is a promising non-planar non-fullerene acceptor for organic photovoltaic applications, but it has a relatively low electron mobility that may limit its performance. Here, we explored the fluorination of peripheral aryl groups to increase intermolecular cofacial π–π stacking interactions, which are desirable for electron transport. Here complexes with fluorine on the distal phenyls [Zn(1F-L2) 2 ], on the naphthyls [Zn(2F-L2) 2 ], and on both [Zn(3F-L2) 2 ] were synthesized and characterized. All three complexes had similar optical and electrochemical properties. The crystal packing structure of Zn(2F-L2) 2 and Zn(3F-L2) 2 revealed cofacial parallel-displaced π–π stacking between the fluorinated 1-naphthylethynyl groups. Such a cofacial orientation was not observed in Zn(L2) 2 crystals, suggesting that fluorination of the naphthyl groups promotes the cofacial π–π stacking orientation. The hole mobility increased from 1.0 × 10 –4 cm 2 V –1 s –1 for Zn(L2) 2 to 0.8–1.0 × 10 –3 cm 2 V –1 s –1 for the fluorinated complexes. Fluorination on the naphthyl groups increased the electron mobility from 4.2 × 10 –5 cm 2 V –1 s –1 for Zn(L2) 2 and Zn(1F-L2) 2 to 2.0 × 10 –4 cm 2 V –1 s –1 for Zn(2F-L2) 2 and Zn(3F-L2) 2 , consistent with cofacial π–π stacking being favorable for electron transport. The three complexes were tested in OPVs using regioregular poly(3-hexylthiophene) (P3HT) as the p-type material, and the best power conversion efficiencies were 5.2, 5.4, and 5.8% for Zn(2F-L2) 2 , Zn(1F-L2) 2 , and Zn(3F-L2) 2 , respectively, compared to 5.5% for Zn(L2) 2 . The fluorination combination found in Zn(3F-L2) 2 resulted in the best device performance. This study points to a viable strategy to increase the electron mobility and performance of non-planar zinc(II) complexes of azadipyrromethene.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Promotion of CaCO 3 Nucleation by Carboxyl- and Amine-Terminated Peptoid Nanotubes

The use of surfaces to promote heterogeneous nucleation of CaCO 3 has been extensively studied to better understand natural processes of biomineralization, to develop bioinspired approaches for synthesizing composite materials with diverse functionalities, and as a route toward accelerating CO 2 mineralization. Peptoids have emerged as a material of particular interest, because they offer tunable side chain chemistry while maintaining a constant long-range self-assembled architecture. Here, in this work, we investigate the ability of films of self-assembled -NH 2 and −COOH terminated peptoid nanotubes to promote the nucleation of CaCO 3 . Using a combination of atomic force, scanning electron, and optical microscopy, we find that interwoven films of nanotubes formed from peptoids with -NH 2 side chains promote the preferential formation of vaterite with the (001) plane parallel to the film and exhibit a lower interfacial energy than that reported for inorganic surfaces. In contrast, the film of −COOH-functionalized peptoid nanotubes is a much stronger promoter, giving a particle density 5 orders of magnitude larger than for the -NH 2 terminated peptoid nanotubes and enabling complete coverage of individual peptoid tubes with CaCO 3 nanoparticles that results in a composite material with a 50-fold increase in Young’s modulus. Thus, interwoven films of peptoid nanotubes provide a useful platform for investigating the impact of specific chemical groups on nucleation and offer potential application to accelerated mineralization of CO 2 for removal from the environment or modulation of mineralogy and properties of carbonate-based building materials.

carbon nanotubes↗
DOE OSTI2024

Conceptual design of a Doppler spectrometer for 102 m/s cross-field flows in tokamak divertors

The cross-field transport in the scrape-off-layers (SOLs) and divertors in tokamaks is of a similar size to the poloidal component of the parallel flow, thereby significantly impacting the plasma transport there. However, its direct observation has been challenging because the drift velocity (102–103 m/s) is significantly below the detection limit of conventional diagnostics. To realize cross-field ion flow measurement, a variety of systematic uncertainties in the system must be narrowed down. Here, we develop a conceptual design of the Doppler spectrometry that enables us to measure the impurity flows with 102-m/s accuracy based on an in situ wavelength-calibration technique developed in the astrophysics field, the iodine-cell method. We discuss its properties and applicability. In particular, the scaling relation between wavelength accuracy and various spectroscopic parameters is newly presented, which suggests the high importance of the wavelength resolution of the system. Based on transport simulations for the JT-60SA divertor, the feasibility of the system is assessed.

Fujii, Keisuke↗
DOE OSTI2023

Exciton Fine Structure in 2D Perovskites: The Out‐of‐Plane Excitonic State

2D Ruddlesden‐Popper metal‐halide perovskites feature particularly strong excitonic effects, making them a fascinating playground for studying exciton physics. A complete understanding of the properties of this quasi‐particle is crucial to fully exploit the tremendous potential of 2D perovskites (2DP) in light emission applications. Despite intense investigations, some of the exciton properties remain elusive to date, for example, the energy‐ordering of the exciton states within the so‐called fine structure manifold. Using optical spectroscopy, it demonstrates that in the archetypical 2DP (PEA) 2 PbI 4 , in contradiction to theoretical predictions, the energy of the bright out‐of‐plane exciton state is higher than that of two in‐plane states. Having elucidated the order of exciton fine structure, it determines the g ‐factor of the dark exciton transition, together with the values of the electron and hole g ‐factors in the direction parallel to the c ‐axis of the crystal. In this way, it provides for the first time, a complete picture of the exciton fine structure in (PEA) 2 PbI 4 2DP.

2D perovskites↗
DOE OSTI2020

2020 IEEE PES Innovative Smart Grid Technologies Europe (ISGT-Europe)

Recent proliferation of distributed energy sources in distribution or sub-transmission systems necessitates close monitoring of these three-phase power grids which typically operate under unbalanced loading conditions. Unlike the transmission systems where the network equations are commonly based on the positive sequence component models, a detailed three phase model will have to be used in implementing network applications for these systems. In the specific case of the state estimator, where measurement and parameter errors may bias the solution, bad data and parameter error detection algorithms should also be incorporated. Implementing the state estimator and error detection algorithms for three-phase systems impose additional computational burden and modifications to the state estimation code. This paper proposes a practical solution to avoid these issues by using synchronized phasor measurements and modal decoupling. The previously developed parameter error detection algorithm based on the normalized Lagrange multipliers (NLM) test is applied to the measurements independently in each mode in parallel, not only saving CPU time but also avoiding new code development for a three-phase estimator. Different parameter error scenarios are created and tested to verify the effectiveness of the proposed error detection approach.

Khalili, Ramtin↗
DOE OSTI2023

Editorial: Advanced water splitting technologies development: Best practices and protocols

As the level of deployment and utilization of renewable energy sources, including wind and solar, continues to rise, large-scale, long-term energy storage technologies that could accommodate weekly and seasonal energy fluctuations will play a significant role in the overall deployment of renewable energies in the future. Harnessing and storing renewable energy resources via electrochemical, photoelectrochemical, or thermochemical processes by converting renewable energy into sustainable (energy storage) fuels have the potential to meet the long-term, terawatt scale energy storage challenge. Renewable hydrogen production is the cornerstone for sustainable fuel production and deep decarbonization of multiple sectors in our society. Cost-competitive clean hydrogen provides value to applications, such as 1) in the transportation sector for fuel cell vehicles, 2) in the electric grid sector for system stability and load balancing, and 3) in the industrial sector with metal refineries, cement production, and biomass upgrading (carbon-free fertilizer production). In addition, coupling clean renewable hydrogen with the carbon and nitrogen cycles enables known and well-established thermal-chemical processes to generate renewable hydrocarbon fuels and ammonia. The Advanced Water Splitting Technologies (AWST): low temperature electrolysis (LTE), high temperature electrolysis (HTE), photoelectrochemical (PEC) and solar thermo-chemical hydrogen (STCH) provide four unique and parallel approaches to produce low cost, low greenhouse gas (GHG) emission hydrogen at scale (Figure 1). Cost competitive clean hydrogen production using these four technologies is a current high priority focus for governments and industry. In June of 2022, the U.S. Department of Energy (DOE) launched the first in a series of Earthshot Initiatives. The Hydrogen Shot, “1 1 1” aims to reduce the cost of clean hydrogen by more than 80% to one dollar per one kilogram in 1 decade ($\$$1/kg H 2 ). The European Green Deal and the International Energy Agency (IEA) have implemented a strong focus on green hydrogen production for a clean and secure energy future.

benchmarking, low temperature electrolysis↗
DOE OSTI2022

ExaCA: A performance portable exascale cellular automata application for alloy solidification modeling

Modeling the as-solidified grain structures that form during alloy processing is a critical component in understanding process-property relationships, particularly for additive manufacturing (AM) where grain structure is very sensitive to processing conditions. While cellular automata (CA)-based models have proven able to predict aspects of microstructure for several alloys and AM process conditions, long run times and large resource sets required limit the utility and the problem size to which existing CA models can be applied. As part of the ExaAM project, an initiative within the Exascale Computing Project (ECP) to develop, test, and optimize an exascale-capable coupled and self-consistent model of AM parts, we developed ExaCA (https://github.com/LLNL/ExaCA) for the liquid–solid phase transformation in the wake of AM melt pools. The CA-based code is parallelized using MPI and the Kokkos programming model, the latter enabling simulation on both CPUs and GPUs within a single-source implementation. Here, we detail the steps taken to transform a baseline, MPI-based CA code into one that is performant on CPUs and GPUs. Performance testing of ExaCA on Summit (a pre-exascale machine at Oak Ridge National Laboratory) was used to quantify CPU–GPU speedup comparing with equal numbers of nodes. Testing showed comparable CPU performance to the MPI-only CA code and a 5-20x speedup when running AM-based test problems using GPUs. The improved performance of CA through GPU utilization and the performance portable nature of ExaCA will enable accurate part-scale modeling by harnessing the power of current and future generations of high performance computing resources. Future work will include improving the strong scaling of ExaCA on GPUs by reducing load imbalance associated with the locality of the problem, and continuing performance optimization across exascale hardware.

36 MATERIALS SCIENCE↗
DOE OSTI2020

MARTINI-Compatible Coarse-Grained Model for the Mesoscale Simulation of Peptoids

Peptoids (poly-N-substituted glycines) are a class of synthetic polymers that are regioisomers of peptides (poly-C-substituted glycines), in which the point of side-chain connectivity is shifted from the backbone C to the N atom. Peptoids have found diverse applications as peptidomimetic drugs, protein mimetic polymers, surfactants, and catalysts. Computational modeling is valuable in the understanding and design of peptoid-based nanomaterials. In this work, we report the bottom-up parameterization of coarse-grained peptoid force fields based on the MARTINI peptide force field against all-atom peptoid simulation data. Our parameterization pipeline iteratively refits coarse-grained bonded interactions using iterative Boltzmann inversion and nonbonded interactions by matching the potential of mean force for chain extension. We assure good sampling of the amide bond cis/trans isomerizations in the all-atom simulation data using parallel bias metadynamics. We develop coarse-grained models for two representative peptoids—polysarcosine (poly(N-methyl glycine)) and poly(N-((4-bromophenyl)ethyl)glycine)—and show their structural and thermodynamic properties to be in excellent accord with all-atom calculations but up to 25-fold more efficient and compatible with MARTINI force fields. Here, this work establishes a new rigorously parameterized coarse-grained peptoid force field for the understanding and design of peptoid nanomaterials at length and time scales inaccessible to all-atom calculations.

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

Direct conversion of methane to aromatics and hydrogen via a heterogeneous trimetallic synergistic catalyst

Abstract Non-oxidative methane dehydro-aromatization reaction can co-produce hydrogen and benzene effectively on a molybdenum-zeolite based thermochemical catalyst, which is a very promising approach for natural-gas upgrading. However, the low methane conversion and aromatics selectivity and weak durability restrain the realistic application for industry. Here, a mechanism for enhancing catalysis activity on methane activation and carbon-carbon bond coupling has been found to promote conversion and selectivity simultaneously by adding platinum–bismuth alloy cluster to form a trimetallic catalyst on zeolite (Pt-Bi/Mo/ZSM-5). This bimetallic alloy cluster has synergistic interaction with molybdenum: the formed CH 3 * from Mo 2 C on the external surface of zeolite can efficiently move on for C-C coupling on the surface of Pt-Bi particle to produce C 2 compounds, which are the key intermediates of oligomerization. This pathway is parallel with the catalysis on Mo inside the cage. This catalyst demonstrated 18.7% methane conversion and 69.4% benzene selectivity at 710 °C. With 95% methane/5% nitrogen feedstock, it exhibited robust stability with slow deactivation rate of 9.3% after 2 h and instant recovery of 98.6% activity after regeneration in hydrogen. The enhanced catalytic activity is strongly associated with synergistic interaction with Mo and ligand effects of alloys by extensive mechanism studies and DFT calculation.

03 NATURAL GAS↗
DOE OSTI2020

Efficacy of the Cell Perturbation Method in Large-Eddy Simulations of Boundary Layer Flow over Complex Terrain

A challenge to simulating turbulent flow in multiscale atmospheric applications is the efficient generation of resolved turbulence motions over an area of interest. One approach is to apply small perturbations to flow variables near the inflow planes of turbulence-resolving simulation domains nested within larger mesoscale domains. While this approach has been examined in numerous idealized and simple terrain cases, its efficacy in complex terrain environments has not yet been fully explored. Here, we examine the benefits of the stochastic cell perturbation method (CPM) over real complex terrain using data from the 2017 Perdigão field campaign, conducted in an approximately 2-km wide valley situated between two nearly parallel ridges. Following a typical configuration for multiscale simulation using nested domains within the Weather Research and Forecasting (WRF) model to downscale from the mesoscale to a large-eddy simulation (LES), we apply the CPM on a domain with horizontal grid spacing of 150 m. At this resolution, spurious coherent structures are often observed under unstable atmospheric conditions with moderate mean wind speeds. Results from such an intermediate resolution grid are often nested down for finer, more detailed LES, where these spurious structures adversely affect the development of turbulence on the subsequent finer grid nest. We therefore examine the impacts of the CPM on the representation of turbulence within the nested LES domain under moderate mean flow conditions in three different stability regimes: weakly convective, strongly convective, and weakly stable. In addition, two different resolutions of the underlying terrain are used to explore the role of the complex topography itself in generating turbulent structures. We demonstrate that the CPM improves the representation of turbulence within the LES domain, relative to the use of high-resolution complex terrain alone. During the convective conditions, the CPM improves the rate at which smaller-scales of turbulence form, while also accelerating the attenuation of the spurious numerically generated roll structures near the inflow boundary. During stable conditions, the coarse mesh spacing of the intermediate LES domain used herein was insufficient to maintain resolved turbulence using CPM as the flow develops downstream, highlighting the need for yet higher resolution under even weakly stable conditions, and the importance of accurate representation of flow on intermediate LES grids.

17 WIND ENERGY↗
DOE OSTI2025

Effects of magnetic field assisted heat treatment on the microstructure and mechanical properties of Fe-0.63 %C alloy

This study investigates the influence of an applied magnetic field on the microstructural evolution and mechanical properties of hypoeutectoid steels subjected to heat treatment. Tensile tests and microstructural analysis were performed on samples processed under varying magnetic field strengths (0 T, 5 T, and 9 T) and different austenitization incubation times. The results indicate that the application of a magnetic field alters the fraction of proeutectoid ferrite phase without changing the cooling rates and heat treatment process. Additionally, pearlite microstructural features such as lamellar spacing and misorientation angles exhibit variations under different field strengths. While the pearlite nodule diameter remains largely unaffected, an increase in percentage elongation and strength is observed in the 5 T treated sample, attributed to changes in microstructural features with the magnetic field. Additionally, the percentage elongation is reduced in the samples heat treated with reduced austenitization incubation times. The study further demonstrates that low-angle misorientations increased in the samples taken parallel to the magnetic field direction, influencing the mechanical response. These findings suggest that applying a magnetic field during heat treatment provides an additional driving force for phase transformations, offering a manufacturing process for tailoring microstructures and optimizing mechanical properties. Moreover, integrating magnetic fields in heat treatment processes has potential benefits in energy efficiency.

High magnetic field↗
DOE OSTI2025

PETSc/TAO developments for GPU-based early exascale systems

The Portable Extensible Toolkit for Scientific Computation (PETSc) library provides scalable solvers for nonlinear time-dependent differential and algebraic equations and for numerical optimization via the Toolkit for Advanced Optimization (TAO). PETSc is used in dozens of scientific fields and is an important building block for many simulation codes. During the U.S. Department of Energy’s Exascale Computing Project, the PETSc team has made substantial efforts to enable efficient utilization of the massive fine-grain parallelism present within exascale compute nodes and to enable performance portability across exascale architectures. We recap some of the challenges that designers of numerical libraries face in such an endeavor, and then discuss the many developments we have made, which include the addition of new GPU backends, features supporting efficient on-device matrix assembly, better support for asynchronicity and GPU kernel concurrency, and new communication infrastructure. In conclusion, we evaluate the performance of these developments on some pre-exascale systems as well as the early exascale systems Frontier and Aurora, using compute kernel, communication layer, solver, and mini-application benchmark studies, and then close with a few observations drawn from our experiences on the tension between portable performance and other goals of numerical libraries.

Exascale Computing Project (ECP)↗
DOE OSTI2023

Adaptive, Active Learning, and Multifidelity Monte Carlo Methods in the MOOSE Stochastic Tools Module

MOOSE is an open-source computational platform for constructing multi-physics models and executing them in a massively parallel fashion. It has a stochastic tools module (STM) for forward/inverse uncertainty quantification (UQ) and surrogate modeling. This presentation details some recent developments to the STM with respect to the implementation of adaptive, active learning, and multifidelity Monte Carlo methods for forward UQ of computational models. Specifically, the adaptive Monte Carlo methods include Markov Chain Monte Carlo (MCMC)-driven algorithms like adaptive importance sampling and parallelized subset simulation for statistical QoI estimation, rare events analysis, and stochastic gradient-free optimization. The active learning methods include Gaussian Process (GP) surrogates and their training via Adam optimization, design of acquisition functions, and integration with samplers like Monte Carlo, adaptive importance, and parallelized subset simulation. These active learning methods are also designed to work in a batch mode, wherein, the required calls to the full computational model are executed in parallel whenever a user-specified batch size is met. The multifidelity methods in STM are broadly divided into two categories: hierarchical, where a defined hierarchy exists among the low-fidelity models, and peer, where all the low-fidelity models are treated equally. A GP surrogate is used to learn the differences between the low- and high-fidelity models in both multifidelity categories, and acquisition functions from the active learning classes are used to decide whether to rely on a low-fidelity model or call the expensive high-fidelity model. Alongside the software description and usage, applications are also presented to nuclear engineering computational models including a TRISO nuclear fuel particle, a reactor pressure vessel, and a heat-pipe microreactor.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Neutron transport methods for multiphysics heterogeneous reactor core simulation in Griffin

Griffin is a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). This work discloses the methods, algorithms, and implementation for simulating heterogeneous reactor dynamics models. Griffin utilizes a discontinuous finite-element method with discrete ordinates (DFEM-S ) to discretize the field variable of the multigroup neutron transport equation. Multiphysics feedback is handled using two-step tabulated cross-section methodology. Feedback quantities are evaluated using the MOOSE-MultiApp system to couple various engineering phenomena, such as heat conduction and thermal fluids. The multiphysics DFEM-S system is solved using fixed-point iteration with a fully asynchronous parallel sweeper, unstructured coarse-mesh finite difference acceleration, and a multi-timescale improved quasi-static method scheme. The implementation is applied to a multiphysics microreactor model, with two transients: one initiated by a single heat-pipe failure and another by control drum rotation. Importantly, these examples demonstrate the ability of Griffin to tractably solve the neutron transport equation considering seven independent variables and feedback.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A sensitivity analysis of twinning crystal plasticity finite element model using single crystal and poly crystal Zircaloy

The popularity of crystal plasticity finite element method (CPFEM) models is increasing due to their ability to predict the mechanical response of crystalline materials such as metals and metal alloys more accurately than traditional continuum mechanics models. This is since the crystal plasticity models consider the effect of atomic structure, microstructural morphology, and properties of individual grains. These CPFEM models use a large number of material parameters in order to capture the mesoscale physics which comes with the downside of the tedious calibration process. In this paper, a CPFEM code was developed to include the twinning induced grain reorientation and subsequent crystallographic slip for HPC material. The developed code is incorporated in a large-scale, parallelized nonlinear solver WARP3D. Further, a sensitivity analysis with respect to 22 material parameters was then conducted using single crystal and polycrystal representative volume element (RVE) of Zircaloy material. Loading was applied along five different crystallographic orientations for single crystal RVE and along three directions namely, rolling (RD), transverse (TD), and normal (ND) direction for polycrystal RVE. Results obtained from the sensitivity analysis were used for the calibration of material parameters for Zircaloy. Finally, developed code along with calibrated material parameters was used to investigate the effect of the hydride phase formation in Zircaloy which is a typical case observed for nuclear applications. It was found that the volume fraction of the hydride phase has a significant impact on the mechanical properties of Zircaloy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multiscale spatial analysis of fracture nodes in two dimensions

Spatial arrangement of fractures as a function of scale is an important component of fracture quantification for inferential and predictive modeling. Available methods that analyze fracture spatial arrangement are based on one-dimensional spacing data; therefore, they are limited to semi-parallel fractures. Such methods cannot be applied to fracture networks in higher dimensions, particularly when fractures have different orientations. Here to characterize fracture arrangements in two dimensions, we propose using Ripley’s K-function, as a method of point pattern analysis, to quantify spatial arrangement of fracture nodes. Fracture nodes, such as barycenters, intersection points, and tips, are point-based representations of fracture locations and connectivity within the fracture network. We introduce formulations for isotropic as well as directional analyses of spatial arrangement. In addition, we derive formulations for edge correction in circular and rectangular study domains. Finally, we demonstrate applications of Ripley’s K-function on two natural fracture datasets. Our proposed method supports quantification and characterization of fracture spatial arrangements that allow practitioners to build representative models of fractures in the subsurface.

02 PETROLEUM↗