Search NASA⌕ Search

SEARCH · Search NASA

Results for “Microstructure simulation”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Focused Ion Beam Tomography of Alloy 617 Corroded in Molten Chloride Salt

Materials qualification of reactor structural materials is a critical step in rapid implementation of advanced nuclear reactor technologies, particularly to assess the corrosion performance in these designs. Accelerated qualification of reactor structural materials requires incorporating powerful computational toolsets, such as phase field modelling in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, to predict the evolution of structural materials due to corrosion. Accordingly, computational toolsets will require experimental data generated at appropriate length scales to validate accuracy. Focused ion beam (FIB) provides a high degree of control over manipulation of materials for analytical purposes, including capturing data on the evolution in the microstructure and elemental composition of materials at the mesoscale, an appropriate length scale for phase field modelling of intergranular diffusion phenomena using the MOOSE framework. For instance, the FEI Helios G4 UX dual beam plasma FIB microscope at the Irradiated Materials Characterization Laboratory (IMCL) is capable of backscatter diffraction (EBSD) and energy-dispersive x-ray spectroscopy (EDS) documenting the evolution in the microstructure and elemental composition, respectively. The Helios can perform EDS and EBSD three-dimensionally (3D) using tomography, which is then combined using different software packages to visualize 3D volumes correlating elemental composition to microstructural data. The purpose of this investigation was to develop a streamlined characterization and data processing workflow for 3D tomography studies on the FEI Helios G4 plasma FIB. The investigation is segmented into three parts: 1) Optimizing the data collection workflow, 2) identifying appropriate data processing and visualization software (i.e. DREAM.3D, MIPAR, and VGStudioMax), and 3) establishing an infrastructure for public release. The optimization of the data collection workflow is in collaboration with members of the U220 department to setup formal training on the tomography operation of the G4, through ThermoFisher Scientific, and exploring DREAM.3D, MIPAR, and VGStudioMax data processing/visualization software packages. VGStudioMax currently demonstrates the most promise for future use. Optimization of the data collection and processing workflow is still ongoing. A collaboration with INL High Performance Computing (HPC) established an open-source license for expediting the public release of FIB tomography datasets through HPC. FIB tomography data generated by the G4 will provide comprehensive data for validating 3D phase field mesoscale modelling tools within the MOOSE framework for accelerated qualification of reactor structural materials.

Copeland-Johnson, Trishelle↗
DOE OSTI2024

Effects of Rolling Reduction on Critical Current Density and Microstructure of Bi-2212 Wires

Bi-2212 Rutherford cables have been fabricated into flat racetrack coils and canted-cosine-theta dipole magnets. The performance gap between the magnets made with Rutherford cables and the “short-sample-limit” is about 30%. To better understand the influence of Rutherford cable processing on the strand performance, we studied three Bi-2212 wires with filament architectures of 37 × 18 and 55 × 18 and diameters of 0.8 and 1.0 mm. To simulate the deformation caused by cabling process, the three wires were rolled with thickness reductions ranging from 10% to 30%. The aspect ratios of rolled strands are between 1.29 and 2.05. The low aspect-ratio wire is also an interesting form for fabricating solenoid coils with higher packing density. The round and rolled strands were heat-treated under 50 bar and with maximum heat treatment temperatures of 885.5 °C and 890.5 °C. The rolling deformation reduced filament size uniformity, resulting in filament merging in fully heat-treated wires. It was found that rolling reduction reduced wire critical current density (JE) by 16 to 18%, but the JE decrease saturated at 15 to 20% of the thickness reduction. It is believed that the reduced JE results from the filament merging caused by rolling and non-uniform shrinking during overpressure heat treatment.

Bi-2212 wire↗
DOE OSTI2025

Strongly nonlinear wave propagation in elasto-plastic metamaterials: Low-order dynamic modeling

Nonlinear elastic metamaterials are known to support a variety of dynamic phenomena that enhance our capacity to manipulate elastic waves. Since these properties stem from complex, subwavelength geometry, full-scale dynamic simulations are often prohibitively expensive at scales of interest. Prior studies have therefore utilized low-order effective medium models, such as discrete mass-spring lattices, to capture essential properties in the long-wavelength limit. While models of this type have been successfully implemented for a wide variety of nonlinear elastic systems, they have predominantly considered dynamics depending only on the instantaneous kinematics of the lattice, neglecting history-dependent effects, such as wear and plasticity. Here, to address this limitation, the present study develops a lattice-based modeling framework for nonlinear elastic metamaterials undergoing plastic deformation. Due to the history- and rate-dependent nature of plasticity, the framework generally yields a system of differential-algebraic equations whose computational cost is significantly greater than an elastic system of comparable size. We demonstrate the method using several models inspired by classical lattice dynamics and continuum plasticity theory and explore means to obtain empirical plasticity models for general geometries, thereby gaining insight into the influence of microstructural plasticity on effective material performance, which can be used to improve the design of nonlinear mechanical metamaterials.

Dynamic simulation↗
DOE OSTI2025

Phase-field modeling of radiation-induced composition redistribution: An application to additively manufactured austenitic Fe–Cr–Ni

Multicomponent alloys undergoing irradiation damage develop radiation-induced composition redistribution at point defect sinks such as grain boundaries (GBs) and dislocations. Such redistribution results in undesired changes to their mechanical behavior and corrosion resistance. Additively manufactured alloys proposed for future nuclear applications are expected to demonstrate a distinct response to irradiation owing to their unique microstructure with as-solidified dislocation density and chemical microsegregation. To capture the composition redistribution in such systems, we develop a mesoscale model with coupled evolution of atomic and point defect components in the presence of dislocation density, dislocation heterogeneity, and thermodynamic interactions at the GB. The model is parameterized for an FCC Fe–Cr–Ni alloy as a representative system for austenitic stainless steels, and simulations are performed in 1D and 2D as a function of irradiation temperature, dose, dislocation density, and grain size. Radiation-induced segregation (RIS) characterized by Cr depletion and Ni enrichment is predicted at both the GB and the dislocation cell wall, with RIS being lower in magnitude but wider at the cell wall. Strongly biased absorption of self-interstitials by dislocations is found to suppress Ni enrichment but slightly enhance Cr depletion under certain conditions. Thermodynamic segregation at the GB is predicted to be narrower and opposite in sign to RIS for both Cr and Ni. Importantly, non-monotonic segregation is found to occur when both thermodynamic and RIS mechanisms are considered, providing a novel physical interpretation of experimental observations. The model is expected to serve as a key tool in accelerated qualification of irradiated materials.

additively manufactured microstructure↗
DOE OSTI2025

Fatigue life prediction of powder bed fused–laser beam AlSi10Mg: Incorporating critical defects via crystal plasticity modelling

The current study provides a microstructurally-based computational framework to predict the fatigue life of additive manufactured (AM), i.e., powder bed fused–laser beam (PBF-LB), AlSi10Mg specimens using the crystal plasticity finite element method (CPFEM). The fractography analysis, electron backscatter diffraction (EBSD), uniaxial and cyclic responses, and fatigue life of specimens were used to inform the computational framework. CPFE simulation was used to compute fatigue indicator parameters (FIPs) as fatigue driving forces. A new fatigue criterion is introduced based on FIPs, which was calibrated using experimental fatigue data. The proposed fatigue measure was evaluated versus the specimens with critical defects of various sizes and locations subjected to different stress amplitudes. The results show that the developed framework can capture the fatigue life of samples with different critical defect locations and sizes along with different stress amplitudes for both high-cycle fatigue (HCF) and very high-cycle fatigue (VHCF) regimes.

Additive manufacturing↗
DOE OSTI2026

Correlation of coercivity, microstructure, and surface defects in HDDR-processed Nd-Fe-B powders for bonded magnet applications

Hydrogenation-Disproportionation-Desorption-Recombination (HDDR) is an effective method for producing high coercivity, anisotropic powder from either fresh Nd-Fe-B alloy or recycled Nd-Fe-B magnets for use in bonded magnets. We investigated the impact of surface defects on the coercivity of HDDR-processed Nd-Fe-B using scanning electron microscopy (SEM), magnetic measurements, and micromagnetic simulations. We observed that coercivity decreases as particle size reduces, with SEM revealing surface defects and the detachment of Nd 2 Fe 14 B grains and Nd-rich phases from the particle surface. Micromagnetic simulations indicate that demagnetization initiates at the particle surface, where these defects are most concentrated, leading to reduced coercivity. The reduction in squareness of demagnetization curve, knee point field and coercivity for smaller HDDR particles is attributed to an increased specific surface area, which exhibits reduced nucleation field and weak domain wall pinning field during magnetization reversal. By addressing the role of surface defects in coercivity degradation, this study provides insights for improving both new powder production and recycling strategies, ultimately leading to enhanced performance of bonded magnets and contributing to more sustainable practices in the rare earth supply chain. One potential strategy to enhance the performance of HDDR Nd-Fe-B materials involves reducing the fraction of fine particles ($<$35 µm) and promoting the formation of grain boundary phases on particle surfaces.

Liu, Xubo B. [Ames Laboratory (AMES), Ames, IA (Un↗
DOE OSTI2025

Stochastic parametric skeletal dosimetry model for humans: Pediatric and adult computational skeleton phantoms for internal bone marrow dosimetry

Currently, computational phantoms that simulate skeletal tissues are used in active red bone marrow (AM) internal dosimetry. Up-to-date reference computational phantoms recommended by the ICRP are based on the analysis of CT-images of cadavers. Such phantoms have significant disadvantages. One disadvantage is that the assessment of uncertainty due to the population variability of skeleton dimensions and microstructure results from the limited availability of autopsy material. Another disadvantage is the simplified modelling of cortical layer and bone microarchitecture. A method of stochastic parametric skeletal dosimetry modelling of the bone structures – SPSD modelling – has been developed as an alternative to the ICRP reference phantoms. In the framework of this approach, skeletal phantom parameters are evaluated based on extensively reviewed results of published measurements of real bones. The SPSD approach allows for the assessment of both population-average values and their variability. SPSD-phantoms of the skeleton are modelled in voxel representation. They consist of smaller phantoms of the bone sites – segments – described by simple geometric shapes with uniform microarchitecture parameters. Such segmentation makes it possible to account for non-homogeneous skeletal microarchitecture and to model the bone structure with the required voxel resolution to elaborate suitable skeletal phantoms. The current study presents the parameters of the SPSD skeletal phantoms for the following age-groups: newborn, 1-year-old, 5-year-old, 10-year-old, 15-year-old (male and female), and adult (male and female). This skeletal phantom can be used for dosimetry as an alternative to available reference phantoms for bone-seeking radionuclides. The above-mentioned age- and sex-specific skeletal phantoms are comprised of 289 unique segments. The characteristics of the SPSD phantoms do not contradict published data and are in good agreement with the measurement results of real bones.

Science & Technology - Other Topics↗
DOE OSTI2025

Multi-modal characterization of the B2 phase in the Ta-Re binary system

The energy and transportation industries demand materials that retain their mechanical property at high temperatures. Refractory complex concentrated alloys (RCCAs) with a BCC + B2 microstructure offer a potential solution, where maintaining the high temperature mechanical properties can be achieved by precipitation strengthening. This depends on the B2 phase in RCCAs being thermodynamically stable with a high solvus temperature. Recently, we predicted the high temperature stability of the B2 structure in the Ta-Re binary system, using density functional theory. Here, in this study, we provide experimental evidence for the existence of this phase for the first time, using a Ta 65 Re 35 alloy. Despite Ta-Re binary phase diagrams predicting a single-phase BCC microstructure for Ta 65 Re 35 , we show that a high Z nanoscale secondary phase appears after heat treatment at 1550 °C and 1100 °C. Scanning transmission electron microscopy (STEM) revealed that this phase has a cubic structure and is equiatomic TaRe though B2 superlattice reflections were absent in fast Fourier transforms (FFT) and diffraction patterns (DPs). DP simulations indicate that the B2 TaRe superlattice reflections are up to two orders of magnitude weaker than their fundamental reflections, making their detection challenging via electron microscopy. Neutron diffraction confirmed the second phase had a B2 structure. This study identified a previously unobserved high temperature stable B2 phase in the Ta-Re system, enabling the development of new high temperature BCC + B2 RCCAs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mesh Computing Remote Automatic Workflow

The software suite uses a microservice architecture using Docker and `docker-compose`. The microservices are as follows: 1. User interface. This interface is written in JavaScript using the Svelte framework. It exposes form elements and a 3D visualizer to prompt the user through the definition of microstructure parameters, and setting parameters for mesh generation and refinement. 2. Mesh generator. This is a container running the Python package for DREAM3D to generate a voxelized mesh that represents a microstructure defined by the user in the interface. 3. Cubit runner. This is a secure shell protocol tool that makes the submitting the DREAM mesh to an HPC instance and starts to run Cubit shell commands to smooth the grain boundaries with its `sculpt` library, applies user-defined boundary node sets, and bundles and returns the simulation-ready meshes and input files as a zipped directory.

Harris, BrennanKay↗
DOE OSTI2024

Controlling reversible phase transitions in rare-earth nickelates for novel memory devices

Resistive switching in correlated complex oxides is lucrative for emerging applications in neuromorphic computing, and densely scaled non-volatile memory. Electrical conductance of such complex oxides can be controllable switched across multiple orders of magnitude by either (a) electroforming a conduction channel (e.g., in tungsten oxide), or (b) inducing Mott-Hubbard transition (e.g., in rare-earth nickelates)– both via controlled migration of defects (such as oxygen vacancies) under applied bias. Nevertheless, the promise of such defect-driven electronic transitions are far from realized due to a lack of fundamental understanding of the atomic-scale processes that underlie migration and spatiotemporal evolution of oxygen vacancies over nano-to-mesoscopic length/timescales under applied electric field. In this project, we employ a synergistic integration of density functional theory (DFT) calculations, ab initio/classical molecular dynamics (AIMD/CMD) simulations, machine learning (ML), precision synthesis, and multi-modal X-ray imaging experiments to address this knowledge gap. Such an integrated approach offers to elucidate the correlations between subtle structural distortion and oxidation states; treat localized charge carriers; describe defect/ion transport in the presence of electric field; and, in turn, greatly advance the current understanding of microstructural evolution in complex oxides under applied bias. The fundamental knowledge gained from this work will enable precise control over hierarchical defect structures and unravel new routes to manipulate resistance states in complex oxides. This, in turn, will accelerate design of novel devices with desired set of neural functionalities, and high-speed densely-scaled resistive random access memory technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Supersonic hot jet ablative testing and analysis of boron nitride nanotube hybrid composites

Boron nitride nanotubes (BNNTs) are high-strength, high-modulus nanotubes with high thermal and oxidative stabilities. Two hybrid composites were prepared with satin weave carbon fiber (CF) and resole-type phenolic resin: one with surface layers of BNNTs and one with alternating interlayers of BNNTs. The samples were subjected to hot jet tests that simulate realistic high-pressure-temperature conditions to understand the behavior of BNNTs under high-pressure erosion. Adding BNNTs to CF/phenolic laminates enhanced the ablation resistance by reinforcing the char material and mitigated localized thermal damage. Hybrid laminates exhibited up to 14% lower weight loss, 55% increase in flexural modulus, higher thermal diffusivity, and improved char yield and microstructure compared to CF/phenolic samples. The surface layer hybrid had many surviving nanotubes reinforcing the char and crystalline oxide structures that could mitigate further oxygen diffusion. Further, various characterization methods were used to deduce possible mechanisms and their products, indicating that BNNTs could serve as growth templates for direct crystalline boron oxide formation. Overall, hybrid BNNT/CF/phenolic laminates displayed better ablation resistance and favorable microstructure evolution under high-pressure conditions.

42 ENGINEERING↗
DOE OSTI2025

Quantifying Impacts of Biomass Pelletization on Fast Pyrolysis Using a Single-Particle Reactor, X-ray Computed Tomography, and Computational Modeling

The pore structure and density of lignocellulosic feedstocks dictate intraparticle transport phenomena and thereby play an important role in thermochemical conversion processes such as fast pyrolysis for biofuel and biochemical production. Variations in microstructure are inherent from different biomass species and can be introduced by preprocessing techniques such as cutting and pelletization. Morphological changes also occur during conversion and lead to vastly different pore structures and behavior during pyrolysis, which impact required conversion times and product distributions. The current work presents a comprehensive comparison of fast pyrolysis of neat and pelletized pine feedstocks, which includes single-particle experiments, modeling, and 3D imaging by X-ray computed tomography (XCT). The particle-scale model included anisotropic heat and mass transport in a shrinking particle with pyrolysis reactions based on the CRECK mechanism with boundary conditions informed by reactor-scale simulations of the single-particle reactor. The models were validated by measurements of the temperature and mass loss from single-particle pyrolysis experiments of neat and pelletized pine. Quantitative analysis of XCT geometries revealed that pyrolytic conversion yielded chars with increased porosity and permeability compared to the unpyrolyzed materials, along with decreased tortuosity and anisotropy. Pelletization of the pine feedstock resulted in a much denser, less permeable material, which converted slower and produced more residual char after pyrolysis compared to neat pine. The results from particle modeling revealed that accounting for the dynamic and anisotropic heat and mass transport caused by differences in pore structure is critical to achieving agreement with experimental results. Overall, this study highlights the dramatic differences in conversion behavior imparted by pelletization and the importance of capturing microstructural attributes in computational models to guide the design and optimization of pyrolysis processes for specific biomass feedstocks.

09 BIOMASS FUELS↗
DOE OSTI2024

Transforming an Ionic Conductor into an Electronic Conductor via Crystallization: In Situ Evolution of Transference Numbers and Structure in (La,Sr)(Ga,Fe)O 3–x Perovskite Thin Films

Mixed-conducting perovskites are workhorse electrochemically active materials, but typical high-temperature processing compromises their catalytic activity and chemo-mechanical integrity. Low-temperature pulsed laser deposition of amorphous films plus mild thermal annealing is an emerging route to form homogeneous mixed conductors with exceptional catalytic activity, but little is known about the evolution of the oxide-ion transport and transference numbers during crystallization. Here the coupled evolution of ionic and electronic transport behavior and structure in room-temperature-grown amorphous (La,Sr)(Ga,Fe)O 3-x films as they crystallize is explored. In situ ac-impedance spectroscopy with and without blocking electrodes, simultaneous capturing-synchrotron-grazing-incidence X-ray diffraction, dc polarization, transmission electron microscopy, and molecular dynamics simulations to evaluate isothermal and non-isothermal crystallization effects and the role of grain boundaries on transference numbers is combined. Ionic conductivity increases by ≈2 orders of magnitude during crystallization, with even larger increases in electronic conductivity. Consequently, as crystallinity increases, LSGF transitions from a predominantly ionic conductor to a predominantly electronic conductor. The roles of evolving lattice structural order, microstructure, and defect chemistry are examined. Grain boundaries appear relatively nonblocking electronically but significantly blocking ionically. The results demonstrate that ionic transference numbers can be tailored over a wide range by tuning crystallinity and microstructure without having to change the cation composition.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A microstructural signature of the coesite-quartz transformation: New insights from high-pressure experiments and EBSD

Ultra-high pressure (UHP) metamorphism is difficult to identify in continental crust as few petrological barometers are suitable for dominantly felsic lithologies. In such cases, burial to extreme depths is commonly identified through the preservation of coesite, a high-pressure polymorph of SiO 2 that typically forms at depths exceeding ∼ 100 km (i.e., > 2 GPa pressure). Unfortunately, coesite readily transforms to quartz upon exhumation, meaning that UHP terranes may often be overlooked. While some studies have suggested that quartz may inherit an orientation signature indicative of former coesite, both the specific nature of this signature and the conditions favouring its development remain uncertain. Here, to address this problem, we combine electron backscatter diffraction analysis of natural and experimental samples to explore microstructural evolution across the coesite-quartz phase transformation. We demonstrate that neighbouring domains of quartz commonly feature an 84 ± 4° rotation of [c] axes around the pole of a common {m} plane. This orientation relationship is a product of epitaxy, whereby the {$11\bar{2}2$} Japan twin plane in quartz nucleates on the (010) plane in coesite. In supercell simulations, the nucleation of Japan twins can be explained by the energetically favourable alignment of quartz tetrahedra on parental coesite tetrahedra. Through experiments, we demonstrate that this signature emerges over a broad range of conditions, regardless of the availability of nucleation sites (e.g., grain boundaries) or the density of crystal lattice defects (e.g., dislocations). Overall, our work provides a quantitative and unambiguous tool for identifying UHP terranes from quartz in isolation.

Coesite↗
DOE OSTI2024

Benchmarking machine learning strategies for phase-field problems

Abstract We present a comprehensive benchmarking framework for evaluating machine-learning approaches applied to phase-field problems. This framework focuses on four key analysis areas crucial for assessing the performance of such approaches in a systematic and structured way. Firstly, interpolation tasks are examined to identify trends in prediction accuracy and accumulation of error over simulation time. Secondly, extrapolation tasks are also evaluated according to the same metrics. Thirdly, the relationship between model performance and data requirements is investigated to understand the impact on predictions and robustness of these approaches. Finally, systematic errors are analyzed to identify specific events or inadvertent rare events triggering high errors. Quantitative metrics evaluating the local and global description of the microstructure evolution, along with other scalar metrics representative of phase-field problems, are used across these four analysis areas. This benchmarking framework provides a path to evaluate the effectiveness and limitations of machine-learning strategies applied to phase-field problems, ultimately facilitating their practical application.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Coupled cluster and dislocation dynamics modeling of microstructure evolution in irradiated materials

We develop here a coupled cluster and dislocation dynamics framework to study the microstructure evolution of irradiated materials. The framework not only accounts for the three dimensional diffusion of radiation-generated clusters, but also their interaction with dislocation networks and the resultant climb motion of discrete dislocations within finite crystals. The framework is solved with a superposition solution scheme, and is applied to investigate the evolution of the irradiation-induced dislocation loops in zirconium (Zr), considering the effects of various bias factors including the diffusion anisotropy difference (DAD) of interstitials and interstitial clusters, the dislocation bias of defects to discrete dislocation segments, and the production bias of defects from the radiation cascade. We find that the DAD is the most critical factor influencing the kinetics of the loop evolution in Zr, while the recombination/interaction of mobile defects can induce a strong spatial dependence of the loop evolution together with the DAD. Here, the method is also adopted to study the evolution of interstitial $\langle$a$\rangle$ and vacancy $\langle$c$\rangle$ dislocation loop ensembles consistent with the microstructure observed during irradiation-induced growth of Zr. Our findings not only reveal the spatial dependence of the size and ellipticity of the dislocation loops, but also suggest a limit on the anisotropy factor of interstitials to reproduce the co-growth of $\langle$a$\rangle$ and $\langle$c$\rangle$ loops in zirconium, in good agreement with experimental observations and other simulation results.

Bias factors↗
DOE OSTI2026

Splat Quenching as a High Throughput Rapid Solidification Testing Method for Developing Additive Manufacturing Alloys

Splat quenching as a high throughput rapid solidification screening tool is evaluated. Samples were made using known alloys (SS316, IN625, Ti-5553) to evaluate and establish the microstructures for a variety of alloy systems as well as the predicted cooling rates. The samples are inductively heated and levitated prior to being struck between two platens that produce a high contact pressure and a thin sample resulting in cooling rates on the order of 10 6 to 10 7 K/s. Process parameters were evaluated with analytical models in addition to numerical simulations to provide an effect of process variables such as substrate material, melt superheat, and platen velocity on the resulting solidification. The splat thickness was found to be controlled by platen velocity, melt superheat, as well as the feedstock volume. The sample thickness is the key controlling factor for varying the average cooling rate experienced by the splat quenched sample. The splat quenching techniques can reach regions of rapid solidification space that meet and exceed laser and electron beam techniques across the sample. In conclusion, the results of the study provide a useful foundation in understanding the splat quenching technique and its potential as a low effort tool to explore rapid solidification effects on alloys.

Alloy Design↗
DOE OSTI2025

Dusty Gas Model for Solid Oxide Fuel Cell Fuel Electrode

This model applies the Dusty Gas Model simulate multi-component species transport in SOFC (solid oxide fuel cell) anodes which considers the pressure gradient across the fuel electrode. This studyhas been verified with the analytical solution for different fuel electrode thicknesses and with literature values. The model was developed using the VoronoiFVM platform in Julia which is a built in implicit and semi implicit solver that integrates electrochemical behavior, microstructural effects, and transient analysis for accurate prediction of species transport under varying conditions.

dusty gas model (DGM)↗