SEARCH · Engineering Papers

Results for “Composite models”

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 631 records · Page 35

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Uncertainty in determining carbon dioxide removal potential of biochar

A quantitative and systematic assessment of uncertainty in life-cycle assessment is critical to informing sustainable development of carbon dioxide removal (CDR) technologies. Biochar is the most commonly sold form of CDR to date and it can be used in applications ranging from concrete to agricultural soil amendments. Previous analyses of biochar rely on modeled or estimated life-cycle data and suggest a cradle-to-gate range of 0.20–1.3 kg CO 2 net removal per kg of biomass feedstock, with the range reported driven by differences in energy consumption, pyrolysis temperature, and feedstock sourcing. Herein, we quantify the distribution of CDR possible for biochar production with a compositional life-cycle inventory model paired with scenario-aware Monte Carlo simulation in a 'best practice' (incorporating lower transportation distances, high pyrolysis temperatures, high energy efficiency, recapture of energy for drying and pyrolysis energy requirements, and co-generation of heat and electricity) and 'poor practice' (higher transportation distances, lower pyrolysis temperatures, low energy efficiency, natural gas for energy requirements, and no energy recovery) scenarios. In the best-practice scenario, cradle-to-gate CDR (which is representative of the upper limit of removal across the entire life cycle) is highly certain, with a median removal of 1.4 kg of CO 2 e/kg biomass and results in net removal across the entire distribution. In contrast, the poor-practice scenario results in median net emissions of 0.090 kg CO 2 e/kg biomass. Whether this scenario emits (66% likelihood) or removes (34% likelihood) carbon dioxide is highly uncertain. The emission intensity of energy inputs to the pyrolysis process and whether the bio-oil co-product is used as a chemical feedstock or combusted are critical factors impacting the net carbon dioxide emissions of biochar production, together responsible for 98% of the difference between the best- and poor-practice scenarios.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Discovering tungsten-based composites as plasma facing materials for future high-duty cycle nuclear fusion reactors

Abstract Despite of excellent thermal properties and high sputtering resistance, pure tungsten cannot fully satisfy the requirements for plasma facing materials in future high-duty cycle nuclear fusion reactions due to the coupled extreme environments, including the high thermal loads, plasma exposure, and radiation damage. Here, we demonstrated that tungsten-based composite materials fabricated using spark-plasma sintering (SPS) present promising solutions to these challenges. Through the examination of two model systems, i.e., tungsten-zirconium composite for producing porous tungsten near the surface and dispersoid-strengthened tungsten, we discussed both the strengths and limitations of the SPS-fabricated materials. Our findings point towards the need for future studies aimed at optimizing the SPS process to achieve desired microstructures and effective control of oxygen impurities in the tungsten-based composite materials.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2025

Integrated modeling of RF-induced tungsten erosion at ICRH antenna structures in the WEST tokamak *

This paper introduces STRIPE (Simulated Transport of RF Impurity Production and Emission), an advanced modeling framework developed to analyze material erosion and the global transport of eroded impurities originating from radio-frequency (RF) antenna structures in magnetic confinement fusion devices. STRIPE integrates multiple physics modules: SolEdge3x for scrape-off-layer plasma profiles, COMSOL for 3D RF rectified sheath potentials, RustBCA for erosion yields and surface interactions, and global impurity transport for 3D ion energy-angle distributions and impurity transport. The framework is applied to an ion cyclotron RF-heated L-mode discharge (#57877) in the WEST tokamak, where it predicts a thirty-fold increase in gross tungsten erosion at antenna limiters during the transition from ohmic to ICRH operation. Additionally, under ICRH conditions, a tenfold enhancement in erosion is observed when comparing RF sheath effects to purely thermal sheath conditions. High-charge-state oxygen ions ($\mathrm{O}$ 6+ and above) are identified as the dominant contributors to tungsten sputtering. To validate the model, a synthetic diagnostic tool based on inverse photon efficiency (S/XB coefficients) from the ColRadPy collisional-radiative model enables direct comparison with spectroscopic measurements. Model predictions using a plasma composition of 1% oxygen and 99% deuterium show good agreement with observed W − I (400.9 nm) emission for discharge #57877, supporting the accuracy of the STRIPE framework. This study focuses specifically on gross erosion calculations to demonstrate STRIPE’s capabilities. Future extensions of this work will incorporate net erosion, re-deposition, self-sputtering effects, and whole-device modeling of sputtered tungsten impurity transport. STRIPE is also being applied to other RF-heated linear and toroidal devices, offering valuable insights for antenna design, impurity control, and performance optimization in next-generation fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Interseasonal Connections between the Timing of the Stratospheric Final Warming and Arctic Sea Ice

Connections across seasons in atmospheric circulation and sea ice have long been sought to advance seasonal prediction. This study presents a link between the springtime stratosphere and Arctic sea ice in summer through autumn. The polar stratospheric vortex dominates the winter stratosphere before breaking down each spring, which is called the stratospheric final warming, as solar radiation returns to the pole. Interannual variability of this breakdown is dynamically driven, leading to different springtime tropospheric and surface circulation patterns. To examine the different impacts of delayed and early final warmings, a multimodel composite was generated from selected CMIP5 models. Additionally, regressions were performed on JRA-55 against an index of springtime polar vortex strength. In both the multimodel composites and reanalysis regressions, significant anomalies in sea ice thickness persist several months following an anomalous timing of the final warming. A later final warming or stronger springtime polar stratospheric vortex leads to negative sea ice thickness anomalies in the East Siberian Sea and positive anomalies in the Beaufort Sea in comparison with an earlier final warming or weaker polar vortex. The spring polar stratospheric vortex is related to spring polar surface circulation patterns. The winds associated with this pattern induce anomalous sea ice motion, moving ice from the East Siberian Sea toward the Beaufort Sea. Reduced sea ice in the East Siberian Sea is linked to anomalous warmth over this region in autumn. Our results suggest that the timing of the stratospheric final warming exerts an influence on the tropospheric circulation and sea ice through autumn, which has implications for seasonal climate prediction.

Kelleher, Michael↗

DOE OSTI2024

Modelling the ejection of primary aerosols during the fast pyrolysis of biomass anisotropic particles

A model for the fast pyrolysis of anisotropic biomass particles is presented which considers bubbling dynamics within the liquid intermediate phase (metaplast) and aerosol ejection from this phase. The model employs the population balance equation and the method of moments to estimate the production rate and resultant size distribution of aerosol ejections, incorporating a detailed CRECK reaction mechanism, and considers the effect of anisotropic biomass microstructure on the intraparticle transport of mass and energy. Here, this study investigates the impact of particle size, heating rate (heat transfer coefficient), and lignocellulosic composition on aerosol ejection. The model predicts that, at high heating rates (convective heat transfer coefficient of 359 W/m 2 .K), aerosols can contribute over 20% to the heavy fraction yield in bio-oil for small particles (1 mm diameter, 4 mm length). The model can predict aerosol size distribution and surface area, indicating an average size of 20 μm for bubbles and 5 μm for aerosols during increased bubble production and aerosol ejection rates. These findings are consistent with prior experimental results and provide essential information for future modeling of extra-particle reactions of the aerosols as they progress through the reactor.

09 BIOMASS FUELS↗

DOE OSTI2023

Evaluation of Composite Materials for Heliostat Cost Reduction

This is an evaluation of composite materials for heliostat cost reduction, including the motivation, reference heliostat model, characteristic wind loads, selection of pultruded structural composites, and technoeconomic analysis using equivalent deformation.

composite materials↗

DOE OSTI2023

A comparative study of solute trapping in Fe-(33–45 at%) Cu alloys manufactured by laser directed energy deposition

The high cooling rates in laser directed energy deposition (DED-LB) of alloys lead to substantial amounts of solute trapping as solute atoms cannot diffuse away from the solid/liquid interface before it advances. In some concentrated alloys, this results in supersaturated solid phases that form nanoscale hierarchical microstructures when the solute atoms precipitate out during reheating from subsequent laser passes. We choose the Iron-Copper (Fe-Cu) binary alloy as model system as it is chemically homogeneous in the liquid phase and has negligible solid solubility at room temperature. Two alloys with nominal compositions in atomic (at.) %, Fe67Cu33 and Fe55Cu45, were manufactured using DED-LB. Scanning transmission electron microscopy (STEM), energy dispersive spectroscopy (EDS) and wavelength dispersive spectroscopy (WDS) were used to characterize the nanostructures and heterogeneous chemical compositions. A non-equilibrium solute partitioning model was used to compute the supersaturated chemistries of the constituent phases and validated with experimentally measured compositions. The measured phase compositions of the two alloys were very similar, at roughly 12 at.% Cu and 4 at.% Fe in the α(bcc)-Fe and ε(fcc)-Cu phases respectively, despite having different processing parameters and mechanical behavior. This indicates that the total thermal history, that depends on both the processing parameters and laser scan pattern, plays a stronger role on the final microstructure evolution than just the initial quantity of trapped solute. In addition, we find that current non-equilibrium solute partitioning models applied on the continuum scale fall short of predicting accurate quantitative phase compositions in concentrated alloys, although the qualitative trends are captured correctly.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Unleashing the Potential of Fast Charging Batteries: Leveraging Anion Redox Chemistry in Ni- and Co-Free Cathodes

Designing Li-ion battery cathodes free from critical raw materials such as Co and Ni has a huge technological and societal impact. Though anion redox-based Li-rich oxide cathodes allow designing Co and Ni free cathode compositions, the Li-rich oxides witnessed voltage fade, voltage hysteresis, and irreversible oxygen release despite their high capacity. Conversely, anion redox through highly covalent chalcogenides (S/Se) is emerging due to the improved covalency between metal d and ligand p bands. Here, we investigate the tuning of multi-chalcogen (S/Se) p-band and redox-active metal d-band in a model Li-rich chalcogen composition Li 1.13 Ti 0.57 Fe 0.3 S 2-y Se y (y = 0 - 1) through in-depth electrochemical, X-ray spectroscopy, and DFT-based electronic structure investigations. Introducing the appropriate amount of Se p band character in anion redox sulfides increases interlayer distance and metal - ligand covalency without modifying the original crystal structure, promoting significant electrochemical reversibility through mixed anionic (Se 2- /Se n- , S 2- /S n- , wherein n<2) and cationic (Fe 2+ /Fe 3+ ) redox reactions. Here we show the detailed Fe, S, and Se redox contributions during Li insertion/extraction through X-ray Absorption (XAS) and Hard X-ray Photoemission Spectroscopy (HAXPES) measurements. The orbital tuning approach improves rate capability for more than 10 C charge-discharge rate, exhibiting more than 50% of its original capacity obtained at C/20 rate. The buffer cation in the lattice (Ti 4+ ) remains electrochemically inactive even after significant Se p-band introduction in the sulfide framework. Overall, this work takes advantage of multi-anion redox chemistry to uncover practically demanding fast charging-discharging characteristics in intercalation cathodes. The obtained knowledge of this design can be extended to other oxide and chalcogen cathodes for high performance Li-ion batteries.

25 ENERGY STORAGE↗

DOE OSTI2021

Cement stress and microstructure evolution during curing in semi-rigid high-pressure environments

The process of cement slurry curing into a solid has a direct effect on the initial state of stress in a composite concrete and steel structure. Failure models for critical infrastructure and wellbores could benefit from measurements of this initial stress because it greatly influences structural performance. Here, we present measurements of the initial stress state (total stress, effective stress, and pore pressure) of multiple cement formulations cured under a high-pressure of 40 MPa within a semi-rigid steel pipe. Our results show that cements cured in this restrained environment exhibit dense microstructures with low permeability (nano-Darcy), develop an anisotropic stress state greater than zero and less than the curing pressure, and undergo a staged curing process where hydration and mechanical behaviors evolve from a slurry into a pre-stressed porous medium. Both drained and undrained scenarios were investigated, revealing a feedback from the semi-rigid boundary towards improved cement performance.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Parameter Study of the Running-In Process for the Generic Pebble Bed Reactor (GPBR200)

The run-in period of a pebble bed reactor is complex and difficult to model given significant heterogeneity in core composition, power, and temperature. While it is understood that the initial composition of the core should eventually result in the same equilibrium core composition, the approach to equilibrium can vary significantly depending on factors such as start-up fuel enrichment and power ramp rate. To explore this, a high-fidelity model of the General Pebble Bed Reactor was used to vary power ramp schemes and start-up core compositions. It was found that both the initial core composition and power ramp rate had a significant impact on the flow rate of pebbles during early time steps, with higher ramp rates and low enrichment resulting in non-physical flow rates. Power ramp rate alone was found to dictated maximum pebble power peaking observed during the run-in process, with higher ramp rates resulting in greater peak pebble powers. Start-up fuel enrichment and power ramp were found to both impact total fuel consumption, although the impact of start-up fuel enrichment was generally secondary to ramp rate.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Clumpy structures within the turbulent primordial cloud

ABSTRACT The primordial clouds in the mini-haloes hatch the first generation stars of the Universe, which play a crucial role in cosmic evolution. In this paper, we investigate how turbulence impacts the structure of primordial star-forming clouds. Previous cosmological simulations of the first star formation predicted a typical mass of around $\mathrm{ 100 \, M_\odot }$. This conflicts with recent observations of extremely metal-poor stars, suggesting a lower mass scale of about $\mathrm{25 \, M_\odot }$. The discrepancy may arise from unresolved turbulence in the star-forming cloud, driven by primordial gas accretion during mini-halo formation in the previous simulations. To quantitatively examine the turbulence effect on the primordial cloud formation, we employ the adaptive mesh refinement code Enzo to model the gas cloud with primordial composition, including artificially driven turbulence on the cloud scale and relevant gas physics. This artificially driven turbulence utilizes a stochastic forcing model to mimic the unresolved turbulence inside mini-haloes. Our results show that the turbulence with high Mach number and compressional mode effectively fragments the cloud into several clumps, each with dense cores of $\mathrm{22.7 - 174.9 \, M_\odot }$ that undergo Jeans instability to form stars. Fragmentation caused by intense and compressive turbulence prevents a runaway collapse of the cloud. The self-bound clumps with smaller masses in the turbulent primordial clouds suggest a possible pathway to decrease the theoretical mass scale of the first stars, further reconciling the mass discrepancy between simulations and observations.

Astronomy & Astrophysics↗

DOE OSTI2020

High operating temperature heat transfer fluids for solar thermal power generation. Final report

This is the final technical report for the DOE MURI project on developing high-operating temperature heat transfer fluids for solar thermal power generation. This project aimed to evaluate liquid metals as a heat-transport fluid and potential energy storage medium for concentrated solar powered (CSP) systems. The project consisted of five phases, starting with fluid selection and screen and ending with potential CSP integration. The initial fluid selection arrived at heavy liquid metals as candidate fluids such as lead bismuth eutectic (LBE) as the prime liquid-alloy candidate to be studied, given the conditions above. Fifty-seven liquid compositions were initially identified with liquidus temperatures lower than 200°C, composed of elements that provide no driving force to the reduction of Cr 2 O 3 or Al 2 O 3 passivation layers in pipe material which is required for corrosion prevention. LBE was identified as the primary candidate fluid after these compositions were compared against computational thermodynamic modeling to screen all of the selected systems that would be compatible with our design objectives. Twenty-eight of these compositions were also investigated with the rapid characterization methods; the accuracies of the liquidus temperature and the specific heats predicted by the computations were better than 3 percent and 9 percent, respectively. We conducted numerus corrosion tests of various Fe-Cr-Al based steels in LBE and its variants at temperate as high as 800°C. It was found that minor additions to the liquid did not provide sufficient corrosion protection to warrant their use in real applications. Furthermore, we investigated the detailed mechanism of corrosion between heavy liquid metal and the pipe material at these very high temperatures. A liquid metal loop was constructed to operate at these high temperatures. The loop construction and operation lead to valuable lessons learned that will be needed for industrial scale loops. Furthermore, we learned that electromagnetic pumps at these extreme temperatures will have to be reengineered and redesigned to be used for CSP applications.

14 SOLAR ENERGY↗

DOE OSTI2021

A Multiscale Modeling Approach to Cryocompressed Hydrogen Storage Pressure Vessels – Part II: Constitutive Modeling and Finite Element Analysis

Subjected to thermomechanical cycling, composite cryogenic hydrogen (H2) storage pressure vessels experience high stresses in the carbon fiber (CF)/epoxy overwrap which can lead to vessel failure due to a combination of degradation mechanisms such as matrix cracking, fiber/matrix debonding, delamination, and fiber rupture. The present paper is the second part of a set of two articles on this subject that addresses analysis and design of these pressure vessels. Predictive finite element (FE) modeling capabilities have been used to support a material acceptance process that aims at providing guidance to evaluate specialty resins, vessel liner options, and CF composites through thermomechanical testing. We have applied a multiscale modeling approach recently developed (B.N. Nguyen et al., International Journal of Hydrogen Energy, 2019, https://doi.org/10.1016/j.ijhydene.2019.09.200) to perform three-dimensional (3D) FE analysis and design of a cryo-compressed H2 storage pressure vessel for its operation in a large temperature range from room to cryogenic temperatures. The developed approach termed the micro-meso-macro approach determines the laminar stresses in different layers of the vessel in addition to constituent (i.e., fiber and matrix) stresses that are important to the design of the composite overwrap for its layup and material combination to reduce the risk of vessel failure.

Nguyen, Ba Nghiep↗

DOE OSTI2024

Structure and Reactivity of Binuclear Cu Active Sites in Cu-CHA Zeolites for Stoichiometric Partial Methane Oxidation to Methanol

Aluminosilicate zeolites exchanged with copper ions facilitate partial methane oxidation (PMO) to methanol in stoichiometric oxidation and reduction cycles, yet the identities of active Cu sites and details of the reaction mechanism remain debated. Here, we use the high-symmetry chabazite (CHA) zeolite framework as a model support to probe the relationship between bulk composition, Cu speciation, and response to various oxidizing and reducing treatments. Density functional theory and first-principles thermodynamics combined with statistical models reveal that Cu speciation and composition depend strongly on Al configuration and external gas conditions. Cu-CHA samples were synthesized to survey broad regions of Si/Al and Cu/Al composition space and framework Al proximity. Characterization by in situ X-ray absorption and UV–visible spectroscopy during exposure to different oxidation conditions reveal that the extent of Cu oxidation is sensitive to activation conditions and thus that both kinetic and thermodynamic factors influence Cu oxidizability in a given material. Similar characterizations during CO reduction reveal that CO titrates Cu 2+ in amounts suggesting the presence of both O- and O 2 -bridged species. In contrast, CH 4 and autoreduction (He) treatments reduce similar but smaller numbers of Cu sites than CO, implicating O 2 -bridged Cu dimers as a potential common intermediate in the former reduction pathways. Furthermore, a systematic increase in methanol yields (per Cu) in stoichiometric PMO cycles increase with the fraction of binuclear O x -bridged Cu sites suggests these species as active sites, as depicted in an updated PMO reaction mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be ${12.39}_{-0.98}^{+1.30}$ km and ${2.072}_{-0.066}^{+0.067}$ M ⊙ respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is ${\mathrm{log}}_{10}(T\,[{\rm{K}}])={5.99}_{-0.06}^{+0.05}$ for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Furthermore, our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters (within the prior support) and can thus be computed from the posterior samples.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Fixing Amdahl's Law within the Limits of Accelerated Systems: FALLACY

Closeout report for FALLACY project. The performance of Data Model Convergence Initiative (DMC) applications on parallel machines is far below the limit set by Amdahl’s law. Whether the machine is based on many-core, GPUs, FPGAs, or a heterogeneous combination, usually the most significant bottleneck is accessing data from the memory system. Aligning with DMC’s HW/architecture thrust, this project developed a set of memory-centric tools called ‘MemGaze’ that inform the HW/SW stack about an application’s memory behavior, including data access latency and diagnosing poor data layout and data composition. Our approach uses architectural modeling and analysis of workload data accesses.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Towards stacking fault energy engineering in FCC high entropy alloys

Stacking Fault Energy (SFE) is an intrinsic alloy property that governs much of the plastic deformation mechanisms observed in fcc alloys. While SFE has been recognized for many years as a key intrinsic mechanical property, its inference via experimental observations or prediction using, for example, computationally intensive first-principles methods is challenging. This difficulty precludes the explicit use of SFE as an alloy design parameter. In this work, we combine DFT calculations (with necessary configurational averaging), machine-learning (ML) and physics-based models to predict the SFE in the fcc CoCrFeMnNiV-Al high-entropy alloy space. The best-performing ML model is capable of accurately predicting the SFE of arbitrary compositions within this 7-element system. Finally, this efficient model along with a recently developed model to estimate intrinsic strength of fcc HEAs is used to explore the strength–SFE Pareto front, predicting new-candidate alloys with particularly interesting mechanical behavior.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Design, Modeling, and Testing of a Novel Inductor for Electric Vehicles: Iron Nitride Soft Magnetic Composites

New technology for electric vehicles (EVs) must meet the requirements of higher energy usage, lower costs, and more sustainable source materials. One promising material for EV power system components is iron nitride (IN) soft magnetic composites (SMCs) because of their competitive magnetic properties and high abundance of the source materials. As part of an ongoing program at Sandia National Laboratories, this project focused on using computer modeling to optimize the prototyping process for an iron nitride SMC toroidal inductor to reach a target inductance of 600 μH. Four inductors with different combinations of wiring (26 AWG and 20 AWG) and vol% loading of iron nitride (65 vol% and 50 vol%) were fabricated at Cal Poly and characterized using an LCR meter. These inductors were also modeled using COMSOL Multiphysics™ with the Magnetic Fields module. The inductance data from the experiment and the model show that the 65 vol% IN prototypes and models agree with about 8% difference, while the 50 vol% IN samples show about a 9% difference between the prototype and the model. These results suggest that the model can predict inductance with both accuracy and precision with low confidence for the given sample size of four. An additional parameter of AC resistance is studied but the AC resistance results from the inductors and from the model generally do not agree closely, suggesting that the current model used in the project does not fully capture the mechanisms behind AC resistance of the inductor. With the focus of the project on inductance, the percent difference results of less than 9% across the four inductors that were tested increases confidence in the model’s predictive capabilities for inductance only. Using the inductance results from both the model and experiment, the final suggested inductor design is a 65 vol% core with 150 windings of 20 AWG wire that is 8 cm across and 1.5 cm tall to reach the inductance goal of 600 μH based on analysis using the optimized COMSOLTM model.

33 ADVANCED PROPULSION SYSTEMS↗