SEARCH · Search NASA

Results for “COMSOL finite element model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Reduction of Methane Leaks through Corrosion Mitigation Pre-treatments for Pipelines with Field Applied Coatings

Corrosion of buried, coated steel pipelines transporting natural gas is a significant source of methane emissions, from pipeline venting required for maintenance and repairs and from pipeline leaks and incidents. Corrosion of steel under field applied coatings is an important safety concern for the pipeline industry. This project investigates the application of a field applied alloy over girth welds to mitigate external corrosion of buried coated steel pipelines. Various metallic coating options were considered, which were required to meet several criteria: (1) it must resist corrosion under open-circuit or mild cathodic protection conditions, (2) it must protect the substrate steel, and (3) it must not negatively affect the adhesion of the field coating. Finite element models and lab testing were performed of alloy coating compositions to identify promising alloy types underneath disbonded coatings. Polarization curves of coating alloys were generated to provide the boundary conditions for the COMSOL model to compute potential and current distributions around coated areas. Sacrificial and corrosion-resistant metal alloy coatings were evaluated and optimized using corrosion modeling and laboratory electrochemical testing, where aluminum alloy 5356 (5% Mg) and steel alloy B9 (9% Cr) were selected. Corrosion test coupons were designed and fabricated using thermal spray aluminum 5356 and welded B9 steel overlays on API 5L grade X42 line pipe steel. The corrosion test coupons, with simulated pipe coating damage, were tested in a laboratory soil box and a field pipeline site in Texas for 3-months. Corrosion test coupons were then tested for 6-months at field pipeline sites in Texas and Tennessee to quantify corrosion rates and performance of the aluminum and steel alloys under polyethylene tape and 2-part epoxy coatings, various coating holidays, and with and without cathodic protection.

03 NATURAL GAS↗

NASA NTRS2014

Mean Flow Augmented Acoustics in Rocket Systems

Oscillatory motion in solid rocket motors and liquid engines has long been a subject of concern. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process and gas dynamics. The customary approach to modeling acoustic waves inside a rocket chamber is to apply the classical inhomogeneous wave equation to the combustion gas. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while the acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The converging section of a rocket nozzle, where gradients in pressure, density, and velocity become large, is a notable region where this approach is not applicable. The expulsion of unsteady energy through the nozzle of a rocket is identified as the predominate source of acoustic damping for most rocket systems. An accurate model of the acoustic behavior within this region where acoustic modes are influenced by the presence of a steady mean flow is required for reliable stability predictions. Recently, an approach to address nozzle damping with mean flow effects was implemented by French [1]. This new approach extends the work originated by Sigman and Zinn [2] by solving the acoustic velocity potential equation (AVPE) formulated by perturbing the Euler equations [3]. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, (del squared)(psi) − (lambda/c)(exp 2)(psi) − M(dot)[M(dot)(del)(del(psi))] − 2(lambda(M/c) + (M(dot)del(M))(dot)del(psi)−2(lambda)(psi)[M(dot)del(1/c)]=0 (1) with M as the Mach vector, c as the speed of sound, and lambda as the complex eigenvalue. French apply the finite volume method to solve the steady flow field within the combustion chamber and nozzle with inviscid walls. The complex eigenvalues and eigenvector are determined with the use of the ARPACK eigensolver. The present study employs the COMSOL Multphysics framework to solve the coupled eigenvalue problem using the finite element approach. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluated the gas flow inside of a solid rocket motor using St. Robert's law modeling solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are used by the coefficient form of the mathematics module to determine the eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. Pertinent results from these analyses are the complex valued eigenvalue and eigenvectors. Comparisons are made to the French results to evaluate the modeling approach. A comparison of the French results with that of the present analysis is displayed in figures 1 and 2, respectively. The graphic shows the first tangential eigenvector's real (a) and imaginary (b) values.

Fischbach, Sean R.↗

DOE OSTI2023

Cryogenic System Preliminary Design for a 0.5m-Long, Conduction-Cooled Nb3Sn Undulator Magnet Prototype

There are several NbTi superconducting undulator (SCU) magnets currently in operation at the Advanced Photon Source (APS) at Argonne National Laboratory (ANL). The development on Nb3Sn-based superconducting undulator magnets at APS is underway due to the potential to further enhance the performance of the SCUs. Superconducting undulator magnets need to keep temperature gradients minimized in order to retain thermal and operating current margin. We have designed the thermal links for efficient heat conduction using 3D finite element analysis (FEA) simulation in COMSOL Multiphysics software, which was later used for the calculation of the temperature distribution across a 0.5 m long, conduction-cooled, Nb3Sn undulator magnet prototype that includes both conductive heat transfer and radiative heating components. We have modelled the evolution of the thermal properties of the magnet winding as well as other cold parts during cool-down from ambient temperature with the operation of a SHI RDK-415D cryocooler to examine the estimated time that is needed for the cooldown as well as the baseline temperature we could achieve. A key result was that a maximum coil delta T = 0.11 K along the designed thermal links and a temperature range of 3.42 K-3.53 K for the winding were predicted at steady state. A pair of G10 support rods for carrying the magnet as well as a pair of current leads made of metals were also designed and optimized for the testing system, targeting for an operation temperature of 4 K and a coil current of 950 A DC.

Zhang, Danlu↗

DOE OSTI2025

Comprehensive model for evaluating voltage losses and performance improvements in thin-film photovoltaic devices

Progress of state-of-the-art and next-generation thin-film photovoltaic devices is often stymied by open-circuit voltage (𝑉 oc ) that is significantly lower than theoretical and practical limits. Yet, effectively diagnosing the primary sources of voltage loss remains challenging. Herein, a sequence of device-level characterization techniques and simulations are employed to identify and rank loss mechanisms. For the research-based Cd⁡(Se,Te) device under study, most of the loss was at the front semiconductor heterointerface due to a clifflike conduction-band offset that lowered the recombination activation energy. Additional losses due to band tails were quantified by photoluminescence analysis. The latter provided the absorption coefficient and activation energy reduction associated with band tails as inputs to device models. Simulations showed that alleviating front-interface issues would improve 𝑉 oc , but it would then be limited by bulk recombination. Further improvement of the bulk would then lead to back-contact limitations. Reducing band tails is beneficial in any circumstance. In conclusion, this analysis provides guidance for reaching toward the radiative 𝑉 oc limit.

14 SOLAR ENERGY↗

DOE OSTI2023

Modeling of the effectiveness of novel edge seal designs for fast, low‐Cap‐Ex manufacturing

Abstract In the manufacturing of photovoltaic (PV) modules, the lamination process can take up to 20 min to complete. In this work, new lamination processes are being developed, and have been prototyped, which hope to be able to cut this time down to as little as 30 s. This could provide significant savings in the cost of lamination equipment, floor space, and energy. PV modules are expected to have a lifespan exceeding 20 to 30 years. For moisture‐sensitive PV technologies, the edge seal between the two layers of glass can be the weakest point of its reliability. There is an inherent challenge when evaluating edge seal materials due to their low permeation rates. As part of Colorado State University's Photovoltaic Research and Development 2, work at the National Renewable Energy Laboratory has developed models to evaluate edge seal configurations in glass‐glass PV modules. Here, this new manufacturing process is evaluated for long‐term moisture durability. Different edge seal design options within glass–glass PV modules are explored. Most of these designs are targeting a superstrate on glass configuration, e.g CdTe, but some designs could be used on conventional crystalline Si cells. Using COMSOL finite element simulation software, we investigated the edge seal and interlayer design configurations containing silicone perimeter edge adhesive, desiccated polyisobutylene‐based edge seal, air, and polyolefin while integrating climate conditions equivalent to a hot and humid environment such as Miami, Florida. We found optimized configurations that will allow the module to prevent moisture ingress over 50 years minimizing the amount of time and material used while utilizing polymers that are easily dispensed.

14 SOLAR ENERGY↗

DOE OSTI2024

Berkeley Lab Finite Element Framework (BELFEM) v0.1

The software program, referred to as BELFEM, is a specialized finite element code designed for the magnetodynamic modeling of high-temperature superconducting (HTS) tapes. It incorporates novel mixed finite element formulations, particularly the h-ϕ-formulation with thin-shell simplification, to efficiently simulate larger geometries. This methodology is extremely promising for predicting the electrodynamic performance of HTS tapes used in superconducting cables and magnets, offering the benefit of reduced computational cost. Compared to similar technologies like COMSOL Multiphysics and GetDP, BELFEM's performance benchmarking indicates superior efficiency in its thin-shell implementation. In the future, it will also support features like thermal coupling and inter-tape current sharing, enhancing its utility in research and development, particularly in nuclear fusion applications. The intent is to develop BELFEM as a robust and efficient tool for the HTS community, contributing to the analysis and design of superconducting cables and magnets.

Messe, Christian↗

NASA NTRS2016

Benchmarking of Computational Models for NDE and SHM of Composites

Ultrasonic wave phenomena constitute the leading physical mechanism for nondestructive evaluation (NDE) and structural health monitoring (SHM) of solid composite materials such as carbon-fiber-reinforced polymer (CFRP) laminates. Computational models of ultrasonic guided-wave excitation, propagation, scattering, and detection in quasi-isotropic laminates can be extremely valuable in designing practically realizable NDE and SHM hardware and software with desired accuracy, reliability, efficiency, and coverage. This paper presents comparisons of guided-wave simulations for CFRP composites implemented using three different simulation codes: two commercial finite-element analysis packages, COMSOL and ABAQUS, and a custom code implementing the Elastodynamic Finite Integration Technique (EFIT). Comparisons are also made to experimental laser Doppler vibrometry data and theoretical dispersion curves.

composite materials↗

DOE OSTI2025

Simulations of the fault current limiting operation of a long-length REBCO CORC ® superconducting cable cooled by helium gas

Conductor-on-round-core (CORC®) cables composed of rare-earth-barium-copper-oxide high-temperature superconducting (HTS) tapes are of great interest for power transmission applications due to their many advantages such as high power density, light weight, and low loss. Closed circulation loops of cryogenic helium gas can be used to cool HTS cables down to low temperatures to significantly improve their current-carrying capacity. Coupled circuit-electromagnetic-thermal finite element simulations implemented in the COMSOL Multiphysics package were developed, validated, and then used for simulating the fault current limiting (FCL) performance and the cooling processes of an 8-layer CORC® cable cooled with a flow of cryogenic helium gas. In the simulations, the temperature dependence of the electrical and thermal properties of all component materials is implemented for improved accuracy. To overcome computational challenges caused by the considerable difference in geometrical scales (i.e. few-µm-thick HTS layers versus 10 m-long HTS cable), the model is divided into two separate simulations. The first simulation is performed on the transverse cross-section of the cable to calculate the electric field, heating power and temperature rise in each component of a CORC® cable during FCL operation. The heating power calculated in the first simulation is transferred to the second model to simulate the cooling of a 10 m-long cable after the fault is cleared. The effect of the helium gas flow rate on the cooling process is also investigated to develop strategic approaches for optimizing cooling systems for HTS cables with FCL capability. The simulations indicated that a 40 ms fault with a voltage drop of 20 V m −1 along the cable can result in a temperature increase from 60 K to about 165 K inside the cable, and it takes about 500 s to cool the cable back to nearly 60 K with a flow of cold helium gas at a rate of 5 g s −1 .

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Numerical and Analytical Modeling of the Effect of Cracks on the Self-Inductance of a COTS YJ-41003-TC Toroid

COTS inductors and transformers often contain partial cracks whose effect on inductance, a key performance parameter, have not been carefully studied. In this report, the effects of both partial and complete cracks on the self-inductance of a 100 turn square cross section COTS YJ-41003-TC toroid comprised of J Material was comprehensively investigated using both analytically derived closed form expressions and 3D computational techniques employing commercial codes. Both partial (half-penny) and complete (air gap) cracks of 10 and 25 μm were investigated. The crack is defined as the physical distance between two faces of the toroid's magnetic core, such that the surface normal of either face is along the Φ-direction, in alignment with the B-field. For the purposes of validation, two different approaches were incorporated for both the analytical and numerical models. The two analytical methods are comprised of a first principles approach based on the physics of electromagnetics, as well as linear circuit theory. The former directly utilizes the integral form of Maxwell's equations while the latter exploits the interchangeable relationship between electric and magnetic circuits. Validation within the computational scheme is realized through a code-to-code comparison between commercial solvers, COMSOL Multiphysics and CST, with the former employing the Finite Element Method (FEM) and the latter the Finite Difference Time Domain (FDTD) technique. Sound agreement between all four methods (ie., two analytical and two numerical) is observed, with results indicating that only a perturbation in self-inductance occurs for the half-penny cracks, while a substantial reduction takes place for the case of complete cracks. It is important to note that even though a static μ r is applied, representing the linear region of the BH curve (based on manufacturer specifications), the complete crack results still place a lower conservative bound on the inductance. This follows from the fact that even in the case of a half-penny crack, if the magnetic core portion of the crack approaches saturation, the crack begins to behave like an air gap, or complete crack. When an air gap is introduced into a magnetic core, a substantial reduction in inductance can occur due to the significant difference in permeabilities between the two mediums (ie., μ core >> μ air ). The once intact bulk magnetic core of the toroid essentially begins to behave like an air core.

36 MATERIALS SCIENCE↗

NASA NTRSDate not supplied

Assessment of Acoustic Behavior for Perforate-Over-Large-Cell Liners

This study explores the effects of increasing the cell size for large-cell acoustic liners. Tests are conducted in the NASA Langley Grazing Flow Impedance Tube to evaluate liners with increasingly larger cell dimensions (up to 2” x 4” cross-section). Conventional impedance eduction confirms that liners with 2” x 3” cells (or larger) must be evaluated using nonlocally reacting assumptions. In addition, due to the sound propagation within cavities, the liners must be modeled using higher fidelity techniques. Thus, grazing flow duct and acoustic liner are modeled simultaneously using finite element methods. The facesheet is modeled using a transfer impedance, while the rest of the domain is modeled using the convected Helmholtz equation. The acoustic pressures predicted are shown to compare favorably with those measured in the Grazing Flow Impedance Tube.

large diameter↗

DOE OSTI2022

Finite Element Simulation of the Acoustic Pressure Inside a Beverage Container for Non-Thermal, Ultrasound-based Pasteurization

The purpose of this effort is to investigate whether large acoustic pressure waves can be transmitted inside beverage containers to enable pasteurization. Acoustic waves are known to induce large nonlinear compressive forces and shock waves in fluids, suggesting that compression waves may be capable of damaging bacteria inside beverage containers without appreciably increasingly the temperature or altering the freshness and flavor of the beverage contents. Although a combined process such as thermosonication (e.g., sonication with heating) is likely more efficient, it is instructive to compute the acoustic pressure field distribution inside the beverage container. The COMSOL simulations used two and three-dimensional models of beverage containers placed in a water bath to compute the acoustic pressure field. A limitation of these COMSOL models is that they cannot determine the bacterial lysis efficiency, rather the models provide an indirect metric of bacterial lysis based on the magnitude of the pressure field and its distribution.

42 ENGINEERING↗

DOE OSTI2025

Finite element modeling of electropolishing of radio frequency dipole Nb crab cavity in hydrofluoric-sulfuric acid electrolyte

The superior performance of niobium superconducting radio frequency (SRF) cavities is achieved through meticulous surface treatment, notably via chemical electropolishing, ensuring exceptionally smooth surfaces. While this technique has been extensively employed for cylindrically symmetric structures, addressing more intricate geometries poses significant challenges in achieving uniform polishing and controlled material removal, especially when moving away from retractable cathodes. This presents a multifaceted electrochemical, thermal, and fluid dynamics conundrum. A prime example is the 197 MHz radio frequency dipole (RFD) crabbing cavity proposed for the Electron Ion Collider (EIC) project, exemplifying such complex structures. Our groundbreaking work integrates the localized oxide thickness variation, considering its impact on current distribution and Joule heating, within the framework of multi-physics modeling using the COMSOL® simulation suite. This approach was applied to a comprehensive model of the RFD cavity, allowing us to investigate optimal external cooling water flow conditions necessary for achieving desirable outcomes. In conclusion, this illustrates the increasing utility of such multi-physics codes to greatly aid in designing solutions to challenging processing requirements for increasingly complex accelerator cavities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

Iodine Effects on Medical Isotope Producer Delay & Guard Bed Temperature and Performance in Support of NA-23 Program

In support of NA-23’s objective to accelerate the establishment of a reliable commercial United States (U.S.) domestic supply of 99Mo produced without the use of highly enriched uranium (HEU), the Pacific Northwest National Laboratory (PNNL) was tasked with modeling the performance of activated carbon delay and guard beds that may be typical for non-HEU 99Mo production processes. The specific tasks accomplished by PNNL included calculations of long-term effects of iodine on delay bed performance and development of a finite element simulation to better understand delay bed performance under various scenarios specifically with radioxenon and radioiodine decay heating of the delay beds. The main issue of interest in this study is understanding the effect of iodine on the abatement performance of the beds if iodine is not sufficiently removed from gases prior to reaching the delay beds. Finite element simulations of delay bed performance were developed using COMSOL Multiphysics® Software. These simulations are based on partial differential equations that describe xenon and iodine adsorption, heat transfer and decay of radioxenon and radioiodine isotopes. Heat transfer in the activated carbon beds is simulated with a heat source tied to the decay of radioactive xenon and iodine species. These simulations were used to investigate delay bed temperature and hold up time using radioxenon and radioiodine activities and flowrates expected from a non-HEU production process. A bounding case scenario was simulated, where all the iodine was released into the delay beds due to insufficient guarding of the delay beds. In addition, preliminary investigations of the guard beds were conducted.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2022

Charge carrier motion and effect of fixed oxide charge in a microstructured silicon radiation detector

Signal formation in a microstructured semiconductor neutron detector is more complex than in planar diode geometry. Three-dimensional microstructures are laterally smaller than the ionization cloud length, and the electric fields may be weak enough to exhibit plasma time effects. This work is the first detailed treatment of charge carrier motion in these complex semiconductor devices to replicate the time profile and signal magnitude. Simulations were performed using COMSOL Multiphysics to investigate various parameters that affect the propagation of the charge cloud. It was observed that the size of the simulated three-dimensional structure had an impact on the induced current pulse, indicating the importance of simulation geometry optimization to accurately simulate charge cloud expansion. COMSOL Multiphysics was used to replicate accurate charge creation profiles using energy deposition information imported from radiation transport codes. A detailed simulation methodology is presented to benchmark preamplifier event pulses along with complexities in modeling the charge carrier motion along the etched microstructured trenches with Si–SiO 2 boundary conditions, including fixed oxide charge and interface trapping.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS2012

Finite Element Analysis of Three Methods for Microwave Heating of Planetary Surfaces

In-Situ Resource Utilization will be Ground Breaking technology for sustained exploration of space. Volatiles are present in planetary regolith, but water by far has the most potential for effective utilization. The presence of water at the lunar poles and Mars opens the possibility of using the hydrogen for propellant on missions beyond Earth orbit. Likewise, the oxygen could be used for in-space propulsion for lunar ascent/descent and for space tugs from low lunar orbit to low Earth orbit. Water is also an effective radiation shielding material as well as a valuable expendable (water and oxygen) required for habitation in space. Because of the strong function of water vapor pressure with temperature, heating regolith effectively liberates water vapor by sublimation. Microwave energy will penetrate soil and heat from within, much more efficiently than heating from the surface with radiant heat. This is especially true under vacuum conditions since the heat transfer rate is very low. The depth of microwave penetration is a strong function of the microwave frequency and to a lesser extent on regolith dielectric properties. New methods for delivery of microwaves into lunar and planetary surfaces is being prototyped with laboratory experiments and modeled with COMSOL MultiPhysics. Recent results are discussed.

Ethridge, Edwin↗

DOE OSTI2025

Design workflow of a symmetric traveling wave antenna for fast ion production on DD tokamaks

Initial computational plasma physics scoping and a finite element method antenna modeling design workflow for a symmetric center-fed high-field side high harmonic fast wave traveling wave array (TWA) antenna are reported here. The TWA is designed to generate a test population of fast deuterium ions in an existing D–D tokamak by heating neutral beam deuterium ions, accelerating them from 80 keV to several hundred keV. The resulting fast particles are tailored to mimic key reactor energetic particle parameters with regards to exciting Alfven eigenmode instabilities, allowing for a D–D tokamak like DIII-D or ASDEX-U to replicate reactor-relevant conditions experimentally. Initial scenario scoping for high single-pass absorption as well as good preferential fast ion damping relative to electron damping was completed using the ray-tracing/Fokker–Planck codes GENRAY and CQL3D. Python RF network analysis packages were used to create a custom TWA optimization tool to inform a COMSOL flat antenna design, and Petra-M was used to study cold plasma effects. The TWA produced by this workflow has several novel features when compared to previous TWA studies, including symmetric center feeding, and passive end straps for image current cancellation for reduced impurity production. We show here that the antenna design workflow can readily produce TWA antennas optimized for reflection coefficient, image current cancellation, and launched power spectrum shape; and that a population of fast ions can be generated in the correct region of parameter space, warranting future more detailed studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Simulation and Modeling of Optical Properties of U, Th, Pb, and Co Nanoparticles of Interest to Nuclear Security Using Finite Element Analysis

In this work, the optical characteristics of uranium (U), lead (Pb), cobalt (Co), and thorium (Th) nanoparticles are fashioned and simulated employing the finite element analysis (FEA) approach concerning multiple particle sizes. Applying finite element analysis, it was found that the simulated absorption peaks of electronic excitations of nuclear nanoparticles are red-shifted from 365 nm to 555 nm for U; from 355 nm to 550 nm for Pb; from 415 nm to 610 nm for Co; and from 350 nm to 540 nm for Th, comparing expanding particle sizes from 60 nm to 100 nm (except for Co, which varied from 70 nm to 100 nm). The FEA-simulated optical band gap energies and far-field radiation patterns were also obtained for nuclear materials. The simulation approach in this research enables the prediction of optical properties and design of nuclear materials before manufacture for nuclear security applications.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Modeling Isothermal Reduction of Iron Ore Pellet Using Finite Element Analysis Method: Experiments & Validation

Iron ore pellet reduction experiments were performed with pure hydrogen (H2) and mixtures with carbon monoxide (CO) at different ratios. For direct reduction processes that switch dynamically between reformed natural gas and hydrogen as the reductant, it is important to understand the effects of the transition on the oxide reduction kinetics to optimize the residence time of iron ore pellets in a shaft reactor. Hence, the reduction rates were studied by varying experimental parameters such as the temperature (800, 850 & 900 °C), reactant gas flow rate (100, 150 & 200 cm3/min), pellet size and composition of the reactant gas mixture. The rate of reduction was observed to increase with an increase in temperature and reactant gas flow rate, but it decreased with an increase in pellet size. SEM greyscale analysis was performed to analyze the porosity and phase composition of partially reduced pellets. The porosity of the pellets was observed to increase from 0.3 for unreacted pellet to 0.42 for a completely reduced pellet. Energy-dispersive X-ray spectroscopy (EDAX) analysis was performed to identify the phases observed in the SEM images. The fraction of iron phase was observed to increase from the shell region of the pellet to the core region with an increase in the degree of reduction. A 2D-axisymmetric numerical model was developed on COMSOL Multiphysics, and it was validated using the conversion (X) vs. time curves obtained from each experiment. The model was able to accurately predict the total time needed for the complete conversion of a single iron ore pellet for multiple experiments. Effects of changes in the porosity and tortuosity of the pellet on the model were also studied and the rate of reduction was observed to be sensitive to changes in both porosity and tortuosity. The SEM analysis and the model results show that tortuosity is higher for pellets reduced with H2 than for pellets reduced with H2-CO gas mixtures.

08 HYDROGEN↗