Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical modelling”

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 73 records · Page 4

Numerical Modeling & Optimization of the iProTech Pitching Inertial Pump (PIP) Wave Energy Converter (WEC) (CRADA Final Report)

This project represents a continuation of the collaboration between iProTech and NLR to simulate, optimize and design the iProTech Pitching Inertial Pump (PIP) device. The objectives of this TEAMER project are twofold: 1. Refining the physical characteristics of the existing iProTech PIP WEC-Sim model to enhance the model’s fidelity and include controllable components. Key model enhancements target the inclusion of Coulomb friction, the introduction of a controllable bypass valve, and the replacement of traditional check valves with advanced motorized ones. 2. Exploring traditional and advanced control algorithms. From traditional methods like latching control to cutting-edge reinforcement learning (RL) algorithms, the goal is to ensure the PIP device's adaptability and optimal performance across a range of ocean conditions. NLR is tasked with augmenting the WEC-Sim model and implementing the control algorithms, culminating in performance comparison analyses. iProTech will update their existing 3D models, advise on model improvements, and determine crucial system metrics. WEC-Sim, developed in MATLAB/SIMULINK with Simscape Multibody, is the main piece of software that will be used in this project. Coupled with the MATLAB RL Toolbox, it offers a robust platform for in-depth simulation and optimization of the iProTech PIP device. Building on previous work to explore the PIP design space and optimize its geometry, mass distribution, center of gravity and other key parameters, this project aims to refine iProTech’s existing numerical models and develop effective control algorithms that can seamlessly integrate into their future hardware testing campaigns.

16 TIDAL AND WAVE POWER↗

Ducted Fuel Injection vs. Free-Spray Injection: A Study of Mixing and Entrainment Effects Using Numerical Modeling

Diesel engines are an important technology for transportation of both people and goods. However, historically they have suffered a significant downside of high soot and nitrogen oxides (NO x ) emissions. Recently, ducted fuel injection (DFI) has been demonstrated to attenuate soot formation in compression-ignition engines and combustion vessels by 50% to 100%. This allows for diesel engines to be run at low-NO x emissions that would have otherwise produced significantly more soot due to the soot/NO x tradeoff. Currently the root causes of this soot attenuation are not well understood. To be able to better optimize DFI for use across a variety of engines and conditions, it is important to understand clearly how it works. This study expands on the current understanding of DFI by using numerical modeling under nonreacting conditions to provide insights about the roles of entrainment and mixing that would have been much more challenging to obtain experimentally. We report this study found that DFI enhances charge gas entrainment upstream of the duct and blocks entrainment inside of the duct. Mixing is enhanced by the duct, which results in lower peak equivalence ratios at the exit of the duct.

42 ENGINEERING↗

Numerical Model Development of a Variable-Geometry Attenuator Wave Energy Converter: Preprint

Because the wave energy industry is still in its infancy, an optimal design for wave energy converters (WECs) has yet to be established; more work is needed to explore various cost reduction pathways. The primary cost-reduction pathway considered for this work is the optimization of the geometric profile on an attenuator WEC to maximize power production while, at the same time, minimizing capital expenditures through the use of variable-geometry modules. In this investigation, the variable geometry modules consist of inflatable bags placed on either side of a base central steel cylinder that would be inflated in low moderate sea states to maximize power capture and then deflated in moderate-extreme sea states to minimize wave loading. The numerical model and simulation of the attenuator WEC were developed and completed using WEC-Sim, which is an open-source code that is appropriate for use in evaluating the dynamic response of the different WEC models in operational seas. The power production estimates were obtained from the Wave Energy Prize (WEP) sea states, which are representative of U. S. deployment sites, to calculate the average climate capture width that is used in the WEP ACE calculation. Preliminary capital expenditure costs were obtained assuming the base central steel cylinder mass was equal to the fluid displaced mass, minus the mass of the variable-geometry bags. The additional weight required to offset the additional buoyancy from the variable-geometry bags was assumed to come from the addition of seawater ballast. The variable-geometry attenuator model was found to have a similar power capture efficiency as a fixed-body model, but is expected to have a lower characteristic capital expenditure given its more streamlined profile, which demonstrates that variable-geometry modules may provide a realistic cost-reduction pathway to help design a more cost-competitive WEC.

50 EE - Wind and Water Power Program - Water (EE-4↗

Numerical Model Development of a Variable-Geometry Attenuator Wave Energy Converter

Because the wave energy industry is still in its infancy, an optimal design for wave energy converters (WECs) has yet to be established; more work is needed to explore various cost-reduction pathways. The primary cost-reduction pathway considered for this work is the optimization of the geometric profile on an attenuator WEC to maximize power production while, at the same time, minimizing capital expenditures through the use of variable-geometry modules. In this investigation, the variable-geometry modules consist of inflatable bags placed on either side of a base central steel cylinder that would be inflated in low-moderate sea states to maximize power capture and then deflated in moderate-extreme sea states to minimize wave loading. The numerical model and simulation of the attenuator WEC were developed and completed using WEC-Sim, which is an open-source code that is appropriate for use in evaluating the dynamic response of the different WEC models in operational seas. The power production estimates were obtained from the Wave Energy Prize (WEP) sea states, which are representative of U. S. deployment sites, to calculate the average climate capture width that is used in the WEP ACE calculation. Preliminary capital expenditure costs were obtained assuming the base central steel cylinder mass was equal to the fluid displaced mass, minus the mass of the variable-geometry bags. The additional weight required to offset the additional buoyancy from the variable-geometry bags was assumed to come from the addition of seawater ballast. The variable-geometry attenuator model was found to have a similar power capture efficiency as a fixed-body model, but is expected to have a lower characteristic capital expenditure given its more streamlined profile, which demonstrates that variable-geometry modules may provide a realistic cost-reduction pathway to help design a more cost-competitive WEC.

50 EE - Wind and Water Power Program - Water (EE-4↗

Experiments and numerical modeling of surrogate droplet combustion for a 93 octane gasoline calibration fuel

Gasoline is the world’s most prolific liquid transportation fuel for ground transportation vehicles. Yet, a fundamental understanding of its burning characteristics is still lacking. Complicating this matter is the highly multicomponent nature of gasoline (comprised of hundreds of chemical species) and the seasonable variability of its composition due to variations in the crude oil feedstocks used in refining. The need to bring order to this complex problem has led to development of calibration fuels (Tier II in presented study) which on average will provide repeatable metrics that replicate performance of refined gasoline for various engine designs. This study presents results of experiments using the configuration of an isolated droplet to study the combustion dynamics of calibration fuel and three different surrogates (S4, S5, S6) for the calibration fuel. Droplet burning is under microgravity conditions, where gas transport is determined entirely by fuel evaporation in order to promote spherical symmetry in the burning process. The initial droplet diameters ranged between 0.55 mm and 0.60 mm. The droplets were spark-ignited and their burning histories were recorded by digital video cameras to obtain quantitative data. For each fuel, the evolutions of droplet, flame, and soot-shell diameters were measured to calculate the burn rates, flame-droplet diameter ratios (FSR) and soot-droplet diameter ratios. Additionally, a detailed numerical modeling is presented. This simulation includes multicomponent phase equilibrium, unsteady gas and liquid transport, radiative heat transfer, internal liquid diffusion and detailed combustion chemistry of different fuels during their burning process. Experiment results show that the calibration fuel and surrogates have a close burn rates and FSR, and S5 presents a relatively similar results with calibration fuel. Predictions from the numerical simulation using the surrogate are compared with experiment measurements, which showed a good agreement with the experiment results, including the droplet burning rate and FSR.

Guo, Songtao↗

Numerical Modeling of Anisotropic Elastic-Wave Sensitivity Propagation for Optimal Design of Time-Lapse Seismic Surveys

Reliable subsurface time-lapse seismic monitoring is crucial for many geophysical applications, such as enhanced geothermal system characterization, geologic carbon utilization and storage, and conventional and unconventional oil/gas reservoir characterization, etc. In this work, we develop an elastic-wave sensitivity propagation method for optimal design of cost-effective time-lapse seismic surveys considering the fact that most of subsurface geologic layers and fractured reservoirs are anisotropic instead of isotropic. For anisotropic media, we define monitoring criteria using qP- and qS-wave sensitivity energies after decomposing qP- and qS-wave components from the total elastic-wave sensitivity wavefield using a hybrid time- and frequency-domain approach. Geophones should therefore be placed at locations with significant qP- and qS-wave sensitivity energies for cost-effective time-lapse seismic monitoring in an anisotropic geology setting. Our numerical modeling results for a modified anisotropic Hess model demonstrate that, compared with the isotropic case, subsurface anisotropy changes the spatial distributions of elastic-wave sensitivity energies. Consequently, it is necessary to consider subsurface anisotropies when designing the spatial distribution of geophones for cost-effective time-lapse seismic monitoring. This finding suggests that it is essential to use our new anisotropic elastic-wave sensitivity modeling method for optimal design of time-lapse seismic surveys to reliably monitor the changes in subsurface reservoirs, fracture zones or target monitoring regions.

97 MATHEMATICS AND COMPUTING↗

A simple and effective five-equation two-phase numerical model for liquid-vapor phase transition in cavitating flows

Numerical difficulties, notably the non-monotonic behavior of the Wood speed of sound and the volume fraction positivity, associated with the reduced five-equation two-phase flow model of Kapila et al. (2001) [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, 2001. Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10), 3002–3024] have been resolved in the past through the introduction of a frozen speed of sound and an algebraic approach for mechanical relaxation afforded by a pressure non-equilibrium six-equation model proposed by [R. Saurel, F. Petitpas, R.A. Berry, 2009. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixture, J. Comput. Phys. 228, 1678–1712]. By contrast, it is explored and demonstrated in this work that these difficulties can in fact still be resolved within the numerical scheme for solving the reduced five-equation model by numerically replacing the Wood speed of sound for the estimates of wave speeds in the approximate Riemann solver HLLC with the monotonic mixture speed of sound for a transport five-equation model. For shock interface (artificial mixture separating pure or nearly pure fluids) interaction problems, with the apparent advantage of monotonic behavior of the speed of sound in the interface, the effect of the numerical replacement is also confined to the interface. Differences other than the behavior of the speed of sound within the interface in the solutions due to the replacement diminish with increasing resolution when reasonable solution can be obtained with Wood speed of sound. For cavitating/expansion problems in physical fluid mixture, it is pointed out and explained why acoustics in the numerical solutions still propagate at the Wood speed of sound (therefore consistent with the reduced five-equation model) even though in some cases a much higher speed of sound like the numerical replacement above for solving the reduced five-equation model or the frozen speed of sound for solving a six-equation model is used for the estimates of wave speeds in the HLLC scheme. A variant of the five-equation two-phase flow model by Saurel et al. (2008) [R. Saurel, F. Petitpas, R. Abgrall, 2008. Modelling phase transition in metastable liquids: application to cavitating and flashing flows, J. Fluid Mech. 607, 313–350] is then constructed for liquid-vapor phase transition in cavitating flows. The relaxation toward thermo-chemical equilibrium during phase transition is achieved by solving a simple system of algebraic equations for the equilibrium state variables for better efficiency, following Pelanti and Shyue (2014) [M. Pelanti, K.-M. Shyue, 2014. A mixture-energy-consistent six-equation two-phase numerical model for fluid with interfaces, cavitation and evaporation waves. J. Comput. Phys. 259, 331–357]. Therefore, the current model retains both the simplicity afforded by the five-equation model and the efficiency of the algebraic relaxation solver. An alternative algebraic approach for handling the non-conservative term (the so-called K∇ · u term) in the reduced five-equation model for mechanical equilibrium of a liquid-vapor mixture is also explored by enforcing the thermal equilibrium at the same time. Finally, numerical results of sample tests in both one and two dimensions in the literature as well as that in three dimensions demonstrate the effectiveness and ability of the proposed model to simulate cavitating flows. An interesting mechanism of shock generation by acoustics in water due to phase transition is then found by the numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Droplet burning of n-butyl acetate/n-heptane mixtures in the standard atmosphere: experiments and numerical modeling

Blending petroleum fuels with oxygenates is a common approach to stem the depletion of crude oil while also mitigating the impact of their combustion on the environment. It has recently been considered that n-butyl acetate (BA, C6H12O2 , boiling point of 399K) could be a viable oxygenate additive to diesel fuel. In this application it is important to determine the influence of the fractional amount of BA on burning of the mixture. This presentation considers this problem from an experimental and computational approach using n-heptane (C7H14 , boiling point of 372K) as an essentially single component surrogate for diesel fuel. The burning configuration considered was ostensibly spherical symmetry as promoted by burning droplets under conditions where forced and natural convection effects were minimized. In this configuration the droplet and flame are concentric spheres and soot aggregates are trapped in a shell structure between the droplet and flame. Measurements were made of the droplet, flame and soot shell diameters through video imaging of the burning process. The experimental results were used to validate a detailed numerical model of the spherically symmetric droplet burning process that incorporated a consistent combustion kinetic mechanism for BA/heptane mixtures, comprised of 402 species and 16,872 reactions inclusive of soot chemistry. A model for soot formation was also included in the simulation to predict formation of the soot shell. The simulations agreed well with measured droplet and flame diameters. It was found both experimentally and computationally that adding BA to heptane had a minimal effect on the mixture burning rate, while the flame was positioned closer to the droplet surface and the simulated soot volume fraction decreased. The implications of these results are discussed for using BA as a potential additive to diesel fuel.

09 BIOMASS FUELS↗

Numerical Modeling of CO2 Sequestration within a Five-Spot Well Pattern in the Morrow B Sandstone of the Farnsworth Hydrocarbon Field: Comparison of the TOUGHREACT, STOMP-EOR, and GEM Simulators

The objectives of this study were (1) to assess the fate and impact of CO2 injected into the Morrow B Sandstone in the Farnsworth Unit (FWU) through numerical non-isothermal reactive transport modeling, and (2) to compare the performance of three major reactive solute transport simulators, TOUGHREACT, STOMP-EOR, and GEM, under the same input conditions. The models were based on a quarter of a five-spot well pattern where CO2 was injected on a water-alternating-gas schedule for the first 25 years of the 1000 year simulation. The reservoir pore fluid consisted of water with or without petroleum. The results of the models have numerous broad similarities, such as the pattern of reservoir cooling caused by the injected fluids, a large initial pH drop followed by gradual pH neutralization, the long-term persistence of an immiscible CO2 gas phase, the continuous dissolution of calcite, very small decreases in porosity, and the increasing importance over time of carbonate mineral CO2 sequestration. The models differed in their predicted fluid pressure evolutions; amounts of mineral precipitation and dissolution; and distribution of CO2 among immiscible gas, petroleum, formation water, and carbonate minerals. The results of the study show the usefulness of numerical simulations in identifying broad patterns of behavior associated with CO2 injection, but also point to significant uncertainties in the numerical values of many model output parameters.

STOMP↗

Numerical modeling of ion transport and adsorption in porous media: A pore-scale study for capacitive deionization desalination

In this report a pore-scale model is presented to simulate the dynamic ion transport and adsorption processes in porous electrodes used for capacitive deionization (CDI). The Stokes equation governing water flow is solved using the lattice Boltzmann method and Nernst-Planck equation describing ion transport is solved using the finite volume method. The ion adsorption process is considered at the surface of carbon electrodes. After validation against analytical solutions and published results, the model is used to study the coupled water flow, ion transport and adsorption in both two-dimensional and three-dimensional porous CDI electrodes at the pore scale. The effect of electrode microstructure, electrical potential and flow velocity on the adsorption processes is quantitatively investigated, and the relative importance of various parameters is determined. The presented model can be a powerful numerical tool to quantitatively analyze ion transport and adsorption in porous electrodes, and may provide useful information for the design and optimization of CDI electrodes and operating conditions for desired desalination efficiency, water throughput and cost.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Numerical modeling of hydrogen mixing in a direct-injection engine fueled with gaseous hydrogen

Hydrogen is considered as one of the most promising options to achieve effective decarbonization of the energy and transportation sectors. As such, it has recently been receiving increasing attention because of its promising potential as an energy carrier for advanced energy and propulsion systems. With a focus on internal combustion engines, direct injection (DI) of gaseous hydrogen during the compression stroke offers great potential for high engine efficiency and specific power while reducing the risk of backfiring and pre-ignition issues. Therefore, many experimental and numerical efforts have recently been dedicated to understanding the physical and chemical behaviors of hydrogen in engine during mixing and combustion. This study focuses on computational fluid dynamics (CFD) modeling of the hydrogen DI process in a hydrogen optical research engine. Under the conditions studied, gaseous hydrogen is injected into the combustion chamber via a centrally located single-hole injector at a pressure of 100 bar. Two configurations, namely low-and high-tumble, are investigated to understand the impact of different in-cylinder flow patterns on the fuel-air mixture preparation. Simulations are carried out using the commercial CFD software CONVERGE. Here, the in-cylinder turbulence is modeled with an unsteady Reynolds-averaged Navier-Stokes (URANS) formulation closed by the renormalization group (RNG) k-ε model. Several numerical methods and model constants, including but not limited to turbulent Schmidt number, are evaluated. The numerical results are systematically compared against experimental measurements of velocity and hydrogen concentration fields on the vertical center plane to assess the performance of the CFD model, unveil the physics of hydrogen mixing, and establish best practices for modeling hydrogen DI under relatively high injection pressure conditions.

33 ADVANCED PROPULSION SYSTEMS↗

Future Response of Coastal Wetlands to Environmental Stresses: Algorithm Comparison of Numerical Models

Coastal wetlands are a critical component of the earth system that strongly influence the global water and biogeochemical cycles. They are also likely important sentinel of climate change. Because eco-geomorphological processes have long been recognized to be important for coastal wetland survival under accelerated sea-level rise (SLR), many eco-geomorphology models have been developed to assess the impact of climate change on coastal wetlands. Although these models differ substantially in complexity and numerical methods, few studies have investigated the algorithm-level uncertainties in these models. In this study, we developed a multiple-algorithm model framework of coastal wetlands that represents coastal hydrodynamics (such as water level, significant wave height and bottom shear stress) and four eco-geomorphological processes: mineral accretion, organic matter accretion, storm surge erosion and landward migration. We validated the model at three representative coastal wetland sites (Venice Lagoon, Plum Island Estuary and Hunter Estuary) for hydrodynamics, mineral accretion and organic matter accretion. Through model-data comparison, we showed that the model can well capture the dynamics of hydrodynamical and eco-geomorphological conditions in the study sites. Importantly, analysis of the multiple-algorithm simulations suggests that differences in the process representation of mineral and organic matter accretion may contribute to the recent contradicting predictions of coastal wetland evolution under accelerated SLR.

54 ENVIRONMENTAL SCIENCES↗

TEAMER Technical Support of Ramboll's for Numerical Modeling of WECs to Support OES Task 10: Cooperative Research and Development (Final Report)

National Technology & Engineering Solutions of Sandia, LLC (NTESS) in collaboration with the National Renewable Energy Laboratory (NREL), and with guidance from Ramboll, will perform fluid dynamics simulations to support of The Ocean Energy Systems (OES) Energy Technology Collaboration Program Task 10 Wave Energy Converters (WEC) Modelling Verification and Validation effort. Specific numerical simulations include the fixed device wave impingement studies to benchmark the performance of simulation techniques against physical testing results.

16 TIDAL AND WAVE POWER↗

Numerical simulation of the hot-tail runaway electron production mechanism using CQL3D and comparison with Smith–Verwichte analytical model

Abstract The hot-tail mechanism of runaway electron (RE) production (Harvey et al 2000 Phys. Plasmas 7 4590) is the primary source of RE in the case of rapidly cooling tokamak plasma. Quantifying this mechanism is very important as it can provide most of the post-thermal-quench (TQ) current, or a seed current for the secondary source of RE through the avalanche mechanism. An analytic model which omits pitch-angle scattering is often used in literature for estimating the hot-tail RE density (Smith and Verwichte 2008 Phys. Plasmas 15 072502). In the present study, we use the CQL3D bounce-averaged Fokker–Planck code (Harvey and McCoy 1992 Proc. IAEA Technical Committee Meeting on Advances in Simulation and Modeling of Thermonuclear Plasmas p 527) to test the limits of validity of the model. In particular, we examine the cases of Z = 1 and Z = 18 ions, for sets of different initial temperature, density, electric field and the characteristic time of temperature decay. We show that for Z = 1 plasma, the ratio of RE density computed by CQL3D to that estimated from the model is within 0.6–6.0 in studied cases. For the Z = 18 case, this factor is systematically a much smaller number, typically 0.02–0.6. We suggest a simple correction to the model that narrows down the range of this ratio to 0.3–3.8 in all of the cases, including Z = 1 and Z = 18 plasmas.

Physics↗

Numerical Model of IProTech PIP WEC Device

iProTech PIP wave energy converter (WEC) is a slack moored, single hull device with no moving parts in the water, joints or bearings. This submission includes data of the simulation, reports, and code for the iProTech PIP (WEC) project. The organization of the data included in the provided archive is detailed below and in the data description of the archive. The data teamer-iprotech-nrel folder includes and explains matlab and python code developed to hydrodynamically model the PIP WEC device in WEC-Sim. The subfolders cover the following steps: 1) report: explanatory information on device geometry 2) pip_mesher: python code to generate mesh panels from device profile data 3) wec-sim_models: matlab code to run WEC-Sim The data uploaded is a snapshot as of 11/02/2121 of code residing in a Github repository administered by David Ogden of NREL.

16 TIDAL AND WAVE POWER↗

Bayesian merging of numerical modeling and remote sensing for saltwater intrusion quantification in the Vietnamese Mekong Delta

Saltwater intrusion has become one of the most concerning issues in the Vietnamese Mekong Delta (VMD) due to its increasing impacts on agriculture and food security of Vietnam. Reliable estimation of salinity plays a crucial role to mitigate the impacts of saltwater intrusion. Here, this study developed a hybrid technique that merges satellite imagery with numerical simulations to improve the estimation of salinity in the VMD. The salinity derived from Landsat images and by numerical simulations was fused using the Bayesian inference technique. The results indicate that our technique significantly reduces the uncertainties and improves the accuracy of salinity estimates. The Nash–Sutcliffe coefficient is 0.74, which is much higher than that of numerical simulation (0.63) and Landsat estimation (0.6). The correlation coefficient between the ensemble and measured salinity is relatively high (0.88). The variance of the ensemble salinity errors (5.0 ppt 2 ) is lower than that of Landsat estimation (10.4 ppt 2 ) and numerical simulations (9.6 ppt 2 ). The proposed approach shows a great potential to combine multiple data sources of a variable of interest to improve its accuracy and reliability wherever these data are available.

54 ENVIRONMENTAL SCIENCES↗

Microstructure and CPO evolution of dynamically recrystallized olivine during complex deformation conditions: a full-field numerical modeling approach

The rheological properties of mantle rocks are strongly dependent on their crystallographic preferred orientation (CPO). Olivine CPO, defined by the orientation of seismically fast [100] axes parallel to flow direction, is also thought to be a dominant contributor to seismic anisotropy in the Earth's upper mantle. However, the amount of deformation needed to overprint a new CPO on a pre-existing fabric and the impact of the inherited CPOs on the transient microstructure evolution, remain unknown. This study employs a full-field numerical approach (VPFFT-ELLE) to explore the dynamic recrystallization and microstructural evolution of olivine polycrystalline aggregates under complex deformation conditions. We test four combinations of successive pure shear and simple shear boundary conditions. Findings indicate that inherited CPOs influence subsequent deformation in a manner dependent on the kinematic relationship between successive stages. In all cases, a minor strain increment (ε ∼0.3–0.6) is sufficient to erase the previous microstructure and CPO. However, when deformation conditions change dramatically (e.g., stretching direction changes orthogonally), the intensity of the new CPO developed is significantly lower and strain distribution are specially altered. During a transient strain stage, pre-existing microstructures undergo extensive reworking, especially when deformation conditions are changed dramatically, such as switching from simple shear to a pure shear condition with a parallel shortening direction relative to the stretching direction. We estimate the significance of these results in interpreting observations of seismic velocity anisotropy, concluding that P-wave seismic anisotropy is significantly and positively correlated with the evolution of olivine CPO with deformation history. This research underscores the transient nature of microstructural rearrangement in olivine aggregates and the necessity for caution in interpreting seismic anisotropy in regions with complex deformation histories, as inherited CPOs can influence current fabric development and induce deviation to the present deformation conditions.

36 MATERIALS SCIENCE↗