Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random fields”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Full 3D Kinetic Modeling and Quantification of Positive Streamer Evolution in an Azimuthally Swept Pin-to-Plane Wedge Geometry

Cathode-directed streamer evolution in near atmospheric air is modeled in 3D pin-to-plane geometries using a 3D kinetic Particle-In-Cell (PIC) code that simulates particle-particle collisions via the Direct Simulation Monte Carlo (DSMC) method. Due to the computational challenges associated with a complete 360° volumetric domain, a practical alternative was achieved using a wedge domain and a range of azimuthal angles was explored (5°, 15°, 30°, and 45°) to study possible effects on the streamer growth and propagation due to the finite wedge angle. A DC voltage of 6 kV is administered to a hemispherical anode of radius 100 μm, with a planar cathode held at ground potential, generating an over-volted state with an electric field of 4 MV/m across a 1500 μm gap. The domain is seeded with an initial ion and electron density of 10 18 m -3 at 1 eV temperature confined to a spherical region of radius 100 μm centered at the tip of the anode. The air chemistry model [1] includes standard Townsend breakdown mechanisms (electron-neutral elastic, excitation, ionization, attachment, and detachment collision chemistry and secondary electron emission) as well as streamer mechanisms (photoionization and ion-neutral collisions) via tracking excited state neutrals which can then either quench via collisions or spontaneously emit a photon based on specific Einstein-A coefficients [2, 3]. In this work, positive streamer dynamics are formally quantified for each wedge angle in terms of electron velocity and density as temporal functions of coordinates r, Φ, and z. Applying a random plasma seed for each simulation, particles of interest are tracked with near femtosecond temporal resolution out to 1.4 ns and spatially binned. This process is repeated six times and results are averaged. Prior 2D studies have shown that the reduced electric field, E/n, can significantly impact streamer evolution [4]. We extend the analysis to 3D wedge geometries, to limit computational costs, and examine the wedge angle’s effect on streamer branching, propagation, and velocity. Results indicate that the smallest wedge angle that produced an acceptably converged solution is 30°. The potential effects that a mesh, when under-resolved with respect to the Debye length, can impart on streamer dynamics and numerical heating were not investigated, and we explicitly state here that the smallest cell size was approximately 10 times the minimum λ D in the streamer channel at late times. This constraint on cell size was the result of computational limitations on total mesh count.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Nanocomposites Based on Coil-Comb Diblock Copolymers

The polymer chain architecture is an important factor determining the phase behavior of nanoparticle (NP) assembly in polymer matrices. Block copolymers (BCPs) containing a random copolymer (RCP) block present a convenient variation on traditional BCPs to tune the interaction parameters between the polymer blocks and the nanofillers as well as to evaluate the effect of the chain architecture on the NP arrangements within BCP microdomains. Here, we synthesized BCPs with a coil polystyrene (PS) block and a comb RCP block through reversible addition-fragmentation chain transfer polymerization. The RCP block consists of methyl- and lauryl acrylates, the latter of which confers a long-chain alkyl moiety to favorably interact with alkyl-passivated NPs. BCPs showing lamellar, cylindrical, and mixed morphologies were obtained by varying the volume fractions of the RCP block (f RCP ). In comparison to coil-coil BCP, the coil-comb BCPs show highly asymmetric phase behavior with respect to fRCP, where lamellar morphologies were observed at fRCP from 0.31 to 0.51. NPs in the size of 4-5 nm were successfully incorporated in the RCP block of the BCPs with periodicities of 30-60 nm. An order-to-order phase transition from lamellae to PS cylinders was observed after the addition of only 1-2 vol % 5 nm NPs into the BCP with a periodicity of 25 nm and f RCP of 0.51. Self-consistent field theory-density functional theory simulations qualitatively described the observed morphologies and phase transitions in the nanocomposites. Here, the current study presents a platform to fabricate nanocomposites with NP assemblies in coil-comb BCPs that contain a random copolymer block and provides insight into how polymer chain architectures can affect the phase behavior of BCPs and nanocomposites.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multi-scale investigation of short-range order and dislocation glide in MoNbTi and TaNbTi multi-principal element alloys

Refractory multi-principal element alloys (RMPEAs) are promising materials for high-temperature structural applications. Here, we investigate the role of short-range ordering (SRO) on dislocation glide in the MoNbTi and TaNbTi RMPEAs using a multi-scale modeling approach. Monte carlo/molecular dynamics simulations with a moment tensor potential show that MoNbTi exhibits a much greater degree of SRO than TaNbTi and the local composition has a direct effect on the unstable stacking fault energies (USFEs). From mesoscale phase-field dislocation dynamics simulations, we find that increasing SRO leads to higher mean USFEs and stress required for dislocation glide. The gliding dislocations experience significant hardening due to pinning and depinning caused by random compositional fluctuations, with higher SRO decreasing the degree of USFE dispersion and hence, amount of hardening. Finally, we show how the morphology of an expanding dislocation loop is affected by the applied stress.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Two-Level Sketching Alternating Anderson Acceleration for Complex Physics Applications

We present a novel two-level sketching extension of the Alternating Anderson–Picard (AAP) method for accelerating fixed-point iterations in challenging single- and multiphysics simulations governed by discretized PDEs. Our approach combines a static, physics-based projection that reduces the least-squares (LS) problem to the most informative field (e.g., via Schur-complement insight) with a dynamic, algebraic sketching stage driven by a backward stability analysis under Lipschitz continuity. We introduce inexpensive estimators for stability thresholds and cache-aware randomized selection strategies to balance computational cost against memory access overhead. The resulting algorithm solves reduced LS systems in place, minimizes memory footprints, and seamlessly alternates between low-cost Picard updates and Anderson mixing. Implemented in Julia, our two-level sketching AAP achieves up to 50% time-to-solution reductions compared to standard Anderson acceleration—without degrading convergence rates—on benchmark problems including Stokes, 𝑝-Laplacian, bidomain, and Navier–Stokes formulations at varying problem sizes. These results demonstrate the method’s robustness, scalability, and potential for integration into high-performance scientific computing frameworks. Our implementation is available open source in the AAP.jl library.

Barnafi, Nicolas [University of Chile, Santiago]↗
DOE OSTI2025

Transcripts and genomic intervals associated with variation in metabolite abundance in maize leaves under field conditions

Abstract Plants exhibit extensive environment-dependent intraspecific metabolic variation, which likely plays a role in determining variation in whole plant phenotypes. However, much of the work seeking to use natural variation to link genes and transcript’s impacts on plant metabolism has employed data from controlled environments. Here, we generated and analyzed data on the variation in the abundance of 26 metabolites across 660 maize inbred lines under field conditions. We employ these data and previously published transcript and whole plant phenotype data reported for the same field experiment to identify both genomic intervals (through genome-wide association studies (GWAS)) and transcripts (using both transcriptome-wide association studies (TWAS) and an explainable artificial intelligence (AI) approach based on random forest (RF)) associated with variation in metabolite abundance. Both genome-wide association and random forest-based methods identified substantial numbers of significant associations including genes with plausible links to the metabolites they are associated with. In contrast, the transcriptome-wide association identified only six significant associations. In three cases, genetic markers associated with metabolic variation in our study colocalized with markers linked to variation in non-metabolic traits scored in the same experiment. We speculate that the poor performance of transcriptome-wide association studies in identifying transcript-metabolite associations may reflect a high prevalence of non-linear interactions between transcripts and metabolites and/or a bias towards rare transcripts playing a large role in determining intraspecific metabolic variation.

Mathivanan, Ramesh Kanna↗
DOE OSTI2022

Objective Phenotyping of Root System Architecture Using Image Augmentation and Machine Learning in Alfalfa (Medicago sativa L.)

Active breeding programs specifically for root system architecture (RSA) phenotypes remain rare; however, breeding for branch and taproot types in the perennial crop alfalfa is ongoing. Phenotyping in this and other crops for active RSA breeding has mostly used visual scoring of specific traits or subjective classification into different root types. While image-based methods have been developed, translation to applied breeding is limited. This research is aimed at developing and comparing image-based RSA phenotyping methods using machine and deep learning algorithms for objective classification of 617 root images from mature alfalfa plants collected from the field to support the ongoing breeding efforts. Our results show that unsupervised machine learning tends to incorrectly classify roots into a normal distribution with most lines predicted as the intermediate root type. Encouragingly, random forest and TensorFlow-based neural networks can classify the root types into branch-type, taproot-type, and an intermediate taproot-branch type with 86% accuracy. With image augmentation, the prediction accuracy was improved to 97%. Coupling the predicted root type with its prediction probability will give breeders a confidence level for better decisions to advance the best and exclude the worst lines from their breeding program. This machine and deep learning approach enables accurate classification of the RSA phenotypes for genomic breeding of climate-resilient alfalfa.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Experimental Investigation of Nanosecond and Subnanosecond Pulsed DBD in Atmospheric Air: Fast Imaging and Spectroscopy

Dielectric barrier discharge (DBD), as an easy and simple way of generation of non-thermal plasma, has found a number of applications in variuos fields. However, development of a homogeneous, or uniform, DBD that would operate at atmospheric pressure conditions in atmospheric air, would open a number of new applications in various fields from thin film coatings to plasma medicine. Unfortunately, at atmospheric pressure, a uniform DBD can be easily transformed into a filamentary dielectric DBD; therefore some serious issues arise, such as gas heating due to strong discharges in the random microdischarge channel and non-uniform energy distribution, which adversely affect applications. These issues traditionally are solved by the use of an appropriate working gas composition, an alternating current driving frequency, lowering of gas pressure, etc. The transitions between discharge modes in the same experimental conditions have been thoroughly investigated in nitrogen, rare gases and their mixtures with air and other gases. In many cases (for example, in plasma medicine), these methods, especially those related to gas composition and pressure, may not be applied in a convenient manner. Recent advances in pulsed power technology permitted application of much faster voltage rise times (including the subnanosecond range) and short (few nanoseconds) pulses, and revealed that uniform DBD can, in fact, be generated in atmospheric air. Such discharges are in the great interest for many applications, however to date there still little understanding of the mechanisms of their operation and characteristics. Currently, there is no adequate model of the uniform dielectric barrier discharge development in atmospheric air. Although extensive studies have been performed on understanding of the nature of the pulsed DBD uniformity, until now there is still little understanding of the mechanism of the DBD transition from the filamentary mode to uniform mode. One of the reasons for this is that development of streamers, and later – filaments, occurs on the sub-nanosecond and nanosecond time scales, and therefore requires imaging and other diagnostic techniques with corresponding speed of registration. Recently developed technologies allow such studies. Here we demonstrate that DBD uniformity strongly depends on applied electric field in the discharge gap. More specifically, the discharge uniformity may be achieved in the case when two conditions are satisfied: (1) stong overvoltage in the discharge gap (provided by fast rise times), when anode-directed streamers are formed, and (2) short pulse duration that prevents discharge overheating due to rising conductivity (current) which leads to formation of filaments. We show that by controlling the applied (global) electric field in ns-pulsed DBD, it is possible to control the uniformity of the discharge. In addition, this offers better control of the discharge chemistry due to local changes of electric fields and therefore electron energy distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Supercapacitance and superinductance of TiN and NbTiN films in the vicinity of superconductor-to-insulator transition

We investigate the low-temperature complex impedance of disordered insulating thin TiN and NbTiN films in the frequency region 400 Hz–1 MHz in close proximity to the superconductor–insulator transition (SIT). The frequency, temperature, and magnetic field dependencies of the real and imaginary parts of the impedance indicate that in full accord with the theoretical predictions and earlier observations, the films acquire self-induced electronic granularity and become effectively random arrays of superconducting granules coupled via Josephson links. Accordingly, the inductive component of the response is due to superconducting droplets, while the capacitive component results from the effective Josephson junctions capacitances. The impedance crosses over from capacitive to inductive behavior as films go across the transition.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

An improved method to access initial states in relativistic heavy-ion collisions

Abstract Observables in heavy-ion collisions are generally categorized into centralities, which reflect an average over events within a range of impact parameters including a wide variety of initial-state configurations. A multiple binning method using spectator neutrons within each centrality has been shown to provide access to events with rare initial-state conditions. This work suggests an improvement in quantifying the difference between standard centrality and spectator neutron binning towards accessing the initial-state properties. A selection of events with higher initial-state density at a fixed participating nucleon number was observed to result in larger final-state particle production and smaller elliptic flow. The relative difference between observables in centrality and spectator binning shows reduced sensitivity for the observables dominated by impact parameter fluctuations in the initial state, such as triangular flow. This property renders the spectator binning method a good candidate for separating geometric contributions from random fluctuations in the initial state towards final-state observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Anomaly Identification of Synchronized Voltage Waveform for Situational Awareness of Low Inertia Systems

Inverter-based resources (IBRs) such as photovoltaics (PVs), wind turbines, and battery energy storage systems (BESSs) are widely deployed in low-carbon power systems. However, these resources typically do not provide the inertia needed for grid stability, resulting in a low-inertia power system. IBRs and lack of inertia have been known to cause anomalies such as waveform distortions and wideband oscillations in power systems due to the limited inertia level, leading to increased generation trips and load shedding. Here, to achieve effective anomaly identification, this paper proposes a synchro-waveform-based algorithm utilizing real-time synchronized voltage waveform measurements from waveform measurement units (WMUs). In the proposed method, different physical characteristics, as well as statistical features, are extracted from synchronized voltage waveform measurements to filter anomalies. Then, the anomaly identification approach based on the random forest is developed and deployed into the FNET/GridEye system considering trade-offs among accuracy, computational burden, and deployment cost. Moreover, four WMUs are specially designed and deployed on Kauai Island to receive instantaneous synchronized voltage waveform measurements. To verify the performance of the proposed algorithm, different experiments are carried out with collected field test data. The result demonstrates that the performance of the proposed synchro-waveform-based anomaly categorization algorithm can accurately identify anomalies 95.35% of the time, which has comparable performance among benchmarking algorithms.

Situational awareness↗
DOE OSTI2022

How to estimate soil organic carbon stocks of agricultural fields? perspectives using ex-ante evaluation

Estimating soil organic carbon (SOC) stocks of agricultural fields has a range of important applications from development of sustainable management practices to monitoring carbon stocks. There are many estimation strategies with the potential for more reliable estimates of SOC stock and more efficient use of soil sampling and analysis resources, especially by leveraging readily available auxiliary information such as remote sensing. However, concrete guidance for strategy selection is lacking. This study narrows this gap with a comparison of strategies for estimating deep SOC stock (0–60 cm) in a prototypical field. Using high density SOC stock measurements and simulation, we built on past studies by 1) ex-ante evaluating a large number of strategy options, 2) using a Bayesian approach to quantify the uncertainty of the comparison, and 3) considering multiple Bayesian models to assess sensitivity to this modeling choice. We found that, using readily available auxiliary information, both balanced and stratified sampling offer substantial improvements over simple random sampling. The auxiliary information most important for this improvement is a Sentinel-2 SOC index = blue / (green × red), followed by the topographic wetness index. We found that these results are robust to the choice of mapping method, but that there is uncertainty in the magnitude of improvement. Here, we recommend future studies implement this Bayesian approach for simulated ex-ante evaluation of SOC stock estimation strategies across more fields to investigate the generalizability of these findings.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Multi-site evaluation of stratified and balanced sampling of soil organic carbon stocks in agricultural fields

Estimating soil organic carbon (SOC) stocks in agricultural fields is essential for environmental and agronomic research, management, and policy. Stratified sampling is a classic strategy for estimating mean soil properties, and has recently been codified in SOC monitoring protocols. However, for the specific task of estimating the SOC stock of an agricultural field, concrete guidance is needed for which covariates to stratify on and how much stratification can improve estimation efficiency. It is also unknown how stratified sampling of SOC stocks compares to modern alternatives, notably doubly balanced sampling. To address these gaps, we collected high-density (average of 7 samples ha -1 ) and deep (average of 75 cm) measurements of SOC stocks at eight commercial fields under maize-soybean production in two US Midwestern states. We combined these measurements with a Bayesian geostatistical model to evaluate stratified and balanced sampling strategies that use a set of readily-available geographic, topographic, spectroscopic, and soil survey data. We examined the number of samples needed to achieve a given level of SOC stock estimation accuracy. While stratified sampling using these variables enables an average sample size reduction of 17% (95% CI, 11% to 23%) compared to simple random sampling, doubly balanced sampling is consistently more efficient, reducing sample sizes by 32% (95% CI, 25% to 37%). The data most important to these efficiency gains are a remotely-sensed SOC index, SSURGO estimates of SOC stocks, and the topographic wetness index. We conclude that in order to meet the urgent challenge of climate change, SOC stocks in agricultural fields could be more efficiently estimated by taking advantage of this readily-available data, especially with doubly balanced sampling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Self-assembly of cocontinuous nanostructured copolymer templates with compositional and architectural dispersity. Final Report

Cocontinuous nanostructured materials in which multiple domains of different materials simultaneously span three dimensional space offer opportunities to achieve combined properties not possible with a single homogeneous material. These architectures have importance in a broad range of energy-relevant technologies including batteries, supercapacitors, fuel cells, separation membranes, and catalysts. Achieving such structures in polymeric materials has been of long-standing interest in the field, due to both the inherently attractive properties of cocontinuous polymer morphologies as well as their ability to serve as templates for other functional nanostructured materials. Our work on this project has established that randomly-linked polymer architectures constructed from two immiscible polymer strands provide robust and highly tunable approaches to disordered cocontinuous nanostructures. In particular, we developed a detailed understanding of how the parameters (linker functionality and strand length, asymmetry, and dispersity) of randomly-linked networks controlled the breadth of the cocontinuous window over which both phases remain percolated. We further characterized how these nanostructures undergo orientation, while retaining cocontinuity, under mechanical deformation. We also compared their behavior to that of random multi-block polymers of linear architecture, which show similar propensity to form disordered nanostructures, albeit over narrower ranges of parameter space. Finally, we have explored the development of functional polymer nanostructures and composites based in part on the fundamental understanding obtained via this project.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Fast Algorithm for Scanning Transmission Electron Microscopy Imaging and 4D-STEM Diffraction Simulations

Scanning transmission electron microscopy (STEM) is an extremely versatile method for studying materials on the atomic scale. Many STEM experiments are supported or validated with electron scattering simulations. However, using the conventional multislice algorithm to perform these simulations can require extremely large calculation times, particularly for experiments with millions of probe positions as each probe position must be simulated independently. Recently, the plane-wave reciprocal-space interpolated scattering matrix (PRISM) algorithm was developed to reduce calculation times for large STEM simulations. Here, we introduce a new method for STEM simulation: partitioning of the STEM probe into “beamlets,” given by a natural neighbor interpolation of the parent beams. This idea is compatible with PRISM simulations and can lead to even larger improvements in simulation time, as well requiring significantly less computer random access memory (RAM). We have performed various simulations to demonstrate the advantages and disadvantages of partitioned PRISM STEM simulations. We find that this new algorithm is particularly useful for 4D-STEM simulations of large fields of view. We also provide a reference implementation of the multislice, PRISM, and partitioned PRISM algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Single skyrmion true random number generator using local dynamics and interaction between skyrmions

Magnetic skyrmions are of great interest to both fundamental research and applications in post-von-Neumann computing devices. The successful implementation of skyrmionic devices requires functionalities of skyrmions with effective controls. Here we show that the local dynamics of skyrmions, in contrast to the global dynamics of a skyrmion as a whole, can be introduced to provide effective functionalities for versatile computing. A single skyrmion interacting with local pinning centres under thermal effects can fluctuate in time and switch between a small-skyrmion and a large-skyrmion state, thereby serving as a robust true random number generator for probabilistic computing. Moreover, neighbouring skyrmions exhibit an anti-correlated coupling in their fluctuation dynamics. Both the switching probability and the dynamic coupling strength can be tuned by modifying the applied magnetic field and spin current. Our results could lead to progress in developing magnetic skyrmionic devices with high tunability and efficient controls.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Ferroelectric Al1 − xBxN thin films integrated on Si

Ferroelectric Al0.93B0.07N thin films are prepared (100) Si substrates. The necessary c-axis out-of-plane orientation to observe macroscopic ferroelectric switching was achieved by implementing an initial Ar/N2 plasma treatment, followed by a thin layer of AlN to initiate the desired texture and a 150 nm W layer. The plasma treatment facilitates crystallinity enhancement of the AlN template layer, allowing for subsequent growth of highly oriented W and Al0.93B0.07N layers. The W layer exhibits random in-plane orientation and exclusive (110) out-of-plane orientation with a rocking curve width of 1.4°. When grown on these W surfaces, 175 nm thick Al0.93B0.07N films exhibit random in-plane orientation and exclusive (001) texture with rocking curve full-width-half-max values of 1.6° and RMS roughness values less than 1 nm. Polarization hysteresis measurements show robust hysteresis with coercive field values of 5.4 MV/cm and remanent polarization values of 136 μC/cm2. XPS depth profile analysis suggests that the plasma treatment converts the existing native oxide to a nitrogen rich oxynitride with approximate composition Si3O0.5N3.67. Cross-sectional TEM reveals that the oxynitride interlayer is amorphous and ∼3.4 nm thick, more than double the native oxide thickness measured by multiwavelength ellipsometry, implying that (oxy)nitride growth continues after conversion of the native oxide. This new family of ferroelectric wurtzites is interesting from an integration perspective given their chemical compatibility with mainstream semiconductors. Developing synthesis routes that promote needed texture while preserving compatible processing windows is an important step toward practical integration.

Physics↗
DOE OSTI2026

Single-shot imaging with randomized structured illumination at a free electron laser

Stroboscopic nanoscale imaging with free electron laser light is revolutionizing our understanding of fast dynamics in heterogeneous systems. The short wavelength of X-ray and extreme ultraviolet radiation makes it possible to achieve nanoscale resolution, while resonance with atomic transitions gives access to electronic and magnetic degrees of freedom. Here, we report on our implementation of a recently developed imaging method, randomized probe imaging, at a free electron laser. The advantage of randomized probe imaging over existing methods is its compatibility with extended and strongly scattering samples. Our implementation delivers robust single-shot reconstructions at up to a full-pitch resolution of 400 nm over a field of view with a 40 µm diameter. We also demonstrate single-shot imaging of magnetic domain structures using circular dichroism at resonance, paving the way to future time-resolved studies of magnetic dynamics, shock physics, and the dynamics of collective electronic phases.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Analyzing the Validity of Brazilian Testing Using Digital Image Correlation and Numerical Simulation Techniques

Characterizing the mechanical behavior of rocks plays a crucial role to optimize the fracturing process in unconventional reservoirs. However, due to the intrinsic anisotropy and heterogeneity in unconventional resources, fracture process prediction remains the most significant challenge for sustainable and economic hydrocarbon production. During the deformation tracking under compression, deploying conventional methods (strain gauge, extensometer, etc.) is insufficient to measure the deformation since the physical attachment of the device is restricted to the size of the sample, monitoring limited point-wise deformation, producing difficulties in data retrieval, and a tendency to lose track in failure points, etc. Where conventional methods are limited, the application of digital image correlation (DIC) provides detailed and additional information of strain evolution and fracture patterns under loading. DIC is an image-based optical method that records an object with a camera and monitors the random contrast speckle pattern painted on the facing surface of the specimen. To overcome the existing limitations, this paper presents numerical modeling of Brazilian disc tests under quasi-static conditions to understand the full-field deformation behaviors and finally, it is validated by DIC. As the direct tensile test has limitations in sample preparation and test execution, the Brazilian testing principle is commonly used to evaluate indirectly the tensile strength of rocks. The two-dimensional numerical model was built to predict the stress distribution and full-field deformation on Brazilian disc under compression based on the assumptions of a homogenous, isotropic and linear elastic material. The uniaxial compression test was conducted using the DIC technique to determine the elastic properties of Spider Berea sandstone, which were used as inputs for the simulation model. The model was verified by the analytical solution and compared with the digital image correlation. The numerical simulation results showed that the solutions matched reasonably with the analytical solutions where the maximum deviation of stress distribution was obtained as 14.59%. The strain evolution (normal and shear strains) and displacements along the central horizontal and vertical planes were investigated in three distinguishable percentages of peak loads (20%, 40%, and 90%) to understand the deformation behaviors in rock. The simulation results demonstrated that the strain evolution contours consistently matched with DIC generated contours with a reasonable agreement. The changes in displacement along the central horizontal and vertical planes showed that numerical simulation and DIC generated experimental results were repeatable and matched closely. In terms of validation, Brazilian testing to measure the indirect tensile strength of rocks is still an issue of debate. The numerical model of fracture propagation supported by digital image correlation from this study can be used to explain the fracturing process in the homogeneous material and can be extended to non-homogeneous cases by incorporating heterogeneity, which is essential for rock mechanics field applications.

58 GEOSCIENCES↗