SEARCH · Search NASA

Results for “Mass Transfer Kinetics”

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 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Understanding the Reaction of Nuclear Graphite with Molecular Oxygen

A thorough understanding of oxidation is important when considering the health and integrity of graphite components in graphite reactors. For the next generation of graphite reactors, HTGRs specifically, an unlikely air ingress has been deemed significant enough to have made its way into the licensing applications of many international licensing bodies. While a substantial body of literature exists on nuclear graphite oxidation in the presence of molecular oxygen and significant efforts have been made to characterize oxidation kinetics of various grades, the value of existing information is somewhat limited. Often, multiple competing processes, including reaction kinetics, mass transfer, and microstructural evolution, are lumped together into a single rate expression that limits the ability to translate this information to different conditions. This article reviews the reaction of graphite with molecular oxygen in terms of the reaction kinetics, gas transport, and microstructural evolution of graphite. It also presents the foundations of a model for the graphite-molecular oxygen reaction system that is kinetically independent of graphite grade, and is capable of describing both the bulk and local oxidation rates under a wide range of conditions applicable to air-ingress.

J. Kane, Joshua↗

DOE OSTI2025

Kinetically Dormant Ni‐Rich Layered Cathode During High‐Voltage Operation

Abstract The degradation of Ni‐rich cathodes during long‐term operation at high voltage has garnered significant attention from both academia and industry. Despite many post‐mortem qualitative structural analyses, precise quantification of their individual and coupling contributions to the overall capacity degradation remains challenging. Here, by leveraging multiscale synchrotron X‐ray probes, electron microscopy, and post‐galvanostatic intermittent titration technique, the thermodynamically irreversible and kinetically reversible capacity loss is successfully deconvoluted in a polycrystalline LiNi 0.83 Mn 0.1 Co 0.07 O 2 cathode during long‐term charge/discharge cycling in full cell configuration. Contradicting the dramatic capacity loss, the layered structure remains highly alive even after 1000 cycles at 4.6 V while undergoing a three‐order of magnitude reduction in the mass transfer kinetics, leading to almost fully recoverable capacity under kinetic‐free conditions. Such kinetic dormant behavior after cycling is not simply ascribed to poor chemical diffusion by reconstructed cathode surface but highly synchronizes with the lattice strain evolution stemming from the structural heterogeneity between deeply delithiated layered and degraded rock‐salt phases at high voltage. These findings deepen the degradation mechanism of high‐voltage cathodes to achieve long‐cycling and fast‐charging performance.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Thermal stability and coalescence dynamics of exsolved metal nanoparticles at charged perovskite surfaces

Exsolution reactions enable the synthesis of oxide-supported metal nanoparticles, which are desirable as catalysts in green energy conversion technologies. It is crucial to precisely tailor the nanoparticle characteristics to optimize the catalysts’ functionality, and to maintain the catalytic performance under operation conditions. We use chemical (co)-doping to modify the defect chemistry of exsolution-active perovskite oxides and examine its influence on the mass transfer kinetics of Ni dopants towards the oxide surface and on the subsequent coalescence behavior of the exsolved nanoparticles during a continuous thermal reduction treatment. Nanoparticles that exsolve at the surface of the acceptor-type fast-oxygen-ion-conductor SrTi 0.95 Ni 0.05 O 3–δ (STNi) show a high surface mobility leading to a very low thermal stability compared to nanoparticles that exsolve at the surface of donor-type SrTi 0.9 Nb 0.05 Ni 0.05 O 3–δ (STNNi). Our analysis indicates that the low thermal stability of exsolved nanoparticles at the acceptor-doped perovskite surface is linked to a high oxygen vacancy concentration at the nanoparticle-oxide interface. For catalysts that require fast oxygen exchange kinetics, exsolution synthesis routes in dry hydrogen conditions may hence lead to accelerated degradation, while humid reaction conditions may mitigate this failure mechanism.

25 ENERGY STORAGE↗

DOE OSTI2024

Revealing the role of redox reaction selectivity and mass transfer in current–voltage predictions for ensembles of photocatalysts

Photocatalysts are conceptually simple reaction units where nanoscale semiconductors integrated with catalysts drive a pair of redox reactions on illumination. However, the proximity of reaction sites performing cathodic and anodic reactions poses dire challenges to realize large light-to-fuel conversion efficiencies. In this study, a powerful, yet straightforward, equivalent-circuit detail-balance modeling framework is developed and applied to evaluate the performance of photocatalytic systems featuring multiple light absorbers. Specifically, low bandgap iridium-doped strontium titanate is modeled as a Z-scheme photocatalyst to achieve desirable hydrogen evolution and iron-based redox shuttle oxidation reactions. Our model has unique capabilities to simulate competing redox reactions and address mass-transfer limitations. In a significant departure from state-of-the-art circuit models, our study develops tools to perform load-line analyses by incorporating a net electrochemical load curve that includes both desired and competing redox reactions. Consequently, reaction selectivity is predicted from equivalent circuit models for photocatalytic and photoelectrochemical systems. Our investigation into ensembles comprised of multiple, semi-transparent light absorbers reveals their potential to outperform a single, optically thick light absorber, particularly when operated under mass-transfer-limited conditions. However, this outcome hinges on minimizing mass-transfer rates of select redox species to prevent undesired reactions of hydrogen oxidation and/or redox shuttle reduction. Our findings demonstrate that reaction selectivity can be achieved by tuning asymmetry in redox species mass-transfer even with perfectly symmetric electrocatalytic charge-transfer coefficients. The influences of various kinetic, mass-transfer, and thermodynamic parameters are explored to offer crucial insights for synthesis of the next-generation of photocatalysts and selective coatings, and reactor designs.

25 ENERGY STORAGE↗

DOE OSTI2025

Physics-Informed Neural Network (PINN) Prediction of Mixed Mass-Heat-Crystallization Limited Methane Hydrate Formation and Dissociation in Micro-Confinement

The creation and use of Physics-Informed Neural Networks (PINNs) for simulating the dynamics of methane hydrate formation and dissociation will be presented. The PINN framework's main benefit is its capacity to impose physical consistency with only a partial comprehension of the governing equations. This makes the algorithm especially useful for systems with little experimental evidence or a lack of theoretical knowledge. A strong basis for forecasting methane hydrate behavior over the verified operating ranges of 30.0-80.9 bar pressure and 1.0-4.0 K sub-cooling conditions is provided by the combination of conductive heat transfer equations and mixed mass-transfer–crystallization kinetics. PINNs were more accurate at predicting the mixed mass-heat-crystallization limited kinetics than conventional Artificial Neural Networks (ANNs), demonstrating remarkable predictive accuracy for methane hydrate production over the ANN model. The efficiency of incorporating physical limitations from first principles into machine learning frameworks for methane hydrate crystallizations is reinforced by these findings. For hydrate-related applications in energy generation, carbon sequestration, and climate modelling, our study establishes PINNs as a computational tool that is both scalable and efficient. The proven capacity to close the gap between conventional physics-based simulations and solely data-driven models creates new opportunities for expedited hydrate research and practical applications.

Hartman, Ryan L [NYU Tandon School of Engineering]↗

DOE OSTI2025

Multi-species electrolyte migration kinetics with directed flow

Here, this paper investigates an electrochemical process designed to control mixing kinetics in liquid electrolytes. The process utilizes an electrochemical cell with anode and cathode chambers connected by a channel without a membrane. This configuration depletes multiple electroactive species in the cathode chamber while selectively replenishing them in the anode chamber. A key application is the removal of impurity species from the cathode chamber electrolyte and the selective replenishment of primary electroactive species in the anode chamber. Aqueous cell experiments were conducted to validate the theoretical model, and simulations were performed for molten salt electrolytes. This approach minimizes waste by reducing the need for additional electrolyte supplies and extends the utilization of anion species, contributing to environmentally sustainable electrochemical materials processing.

Electrochemical Process↗

DOE OSTI2026

Intrinsic kinetics of water-inhibited ultra-lean methane oxidation over PtPd-Mg/ θ -Al 2 O 3 catalyst

Here, this study develops intrinsic methane oxidation kinetics for ultra-lean methane conditions in the presence of water over a highly active and stable PtPd–Mg/θ-Al 2 O 3 catalyst. Comprehensive laboratory experiments were conducted over a wide range of methane concentrations (150–1200 ppm CH 4 ), water contents (1–5% H 2 O), and industrially relevant space velocities (80,000 ≤ GHSV ≤ 160,000 h -1 ). These systematic experiments informed a two-dimensional, axisymmetric, multiscale reactor model that was used to develop and validate methane oxidation kinetics under practically relevant conditions, including non-isothermal operation and high conversion regimes. Combined experimental and modeling results revealed significant intraparticle diffusion resistance and transport-induced reaction exotherm at elevated temperatures, which limited catalyst utilization despite high intrinsic activity. These transport effects were explicitly incorporated into the reactor model, enabling accurate estimation of intrinsic kinetic parameters without reliance on conventional effectiveness-factor corrections. The resulting kinetic model successfully captured both kinetically controlled and mass-transfer-limited regimes and reliably predicted CH 4 conversion across broad ranges of temperature, methane concentration, and water content

Heat and mass transfer limitations↗

DOE OSTI2024

Comparative study of machine learning techniques for post-combustion carbon capture systems

Computational analysis of countercurrent flows in packed absorption columns, often used in solvent-based post-combustion carbon capture systems (CCSs), is challenging. Typically, computational fluid dynamics (CFD) approaches are used to simulate the interactions between a solvent, gas, and column's packing geometry while accounting for the thermodynamics, kinetics, heat, and mass transfer effects of the absorption process. These simulations can then be used explain a column's hydrodynamic characteristics and evaluate its CO 2 -capture efficiency. However, these approaches are computationally expensive, making it difficult to evaluate numerous designs and operating conditions to improve efficiency at industrial scales. In this work, we comprehensively explore the application of statistical ML methods, convolutional neural networks (CNNs), and graph neural networks (GNNs) to aid and accelerate the scale-up and design optimization of solvent-based post-combustion CCSs. We apply these methods to CFD datasets of countercurrent flows in absorption columns with structured packings characterized by several geometric parameters. We train models to use these parameters, inlet velocity conditions, and other model-specific representations of the column to estimate key determinants of CO 2 -capture efficiency without having to simulate additional CFD datasets. We also evaluate the impact of different input types on the accuracy and generalizability of each model. We discuss the strengths and limitations of each approach to further elucidate the role of CNNs, GNNs, and other machine learning approaches for CO 2 -capture property prediction and design optimization.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Integration of the fundamental knowledge on solvent-packing interactions into the multiscale framework for column scale design and optimization

The interfacial area, also known as the effective mass transfer area, is a key factor for determining the mass transfer for carbon dioxide (CO 2 ) capture via the chemical absorption process in a packed column, and thus the overall capture efficiency of the packed column. Most of the widely used empirical and semi-empirical models for interfacial area were derived indirectly through absorption mass transfer with simplifications based on fast chemical kinetics. This report presents the comprehensive unique multiscale approach to develop a surrogate model for effective mass transfer area in structured packed columns that accounts local hydrodynamics as well as variation in physical properties, and changes in solid surface characteristics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Accretion onto a Supermassive Black Hole Binary before Merger

While supermassive binary black holes (SMBBHs) inspiral toward merger they may also accrete matter from a surrounding disk. To study the dynamics of this system requires simultaneously describing the evolving spacetime and the magnetized plasma. We present the first relativistic calculation simulating two equal-mass, nonspinning black holes as they inspiral from a 20 M (G = c = 1) initial separation almost to merger. Our results imply important observational consequences: for instance, the accretion rate $\dot{M}$ onto the black holes first decreases and then plateaus, dropping by only a factor of ~3 despite the rapid inspiral. An estimated bolometric light curve follows the same profile, suggesting some merging SMBBHs may be significantly luminous past the predicted circumbinary disk decoupling. The minidisks are nonstandard: Reynolds, not Maxwell, stresses dominate, and they oscillate between two states. In one part of the cycle, "sloshing" streams transfer mass between minidisks, carrying kinetic energy at a rate sometimes as high as the peak minidisk bolometric luminosity. We also discover that episodic accretion drives time-varying minidisk tilts. These complex dynamics all contribute to unique cyclical behavior in the light curves of late-time inspiraling SMBBHs. The poloidal magnetic flux on the black holes is roughly constant at a dimensionless level $\phi$ ~ 2–3, but doubles just before merger; for significant black hole spin, this flux predicts powerful jets with variability driven by binary dynamics, another potentially unique electromagnetic signature. This simulation is the first to employ our multipatch infrastructure PATCHWORK MHD, decreasing the computational expense to ~3% of conventional single-grid methods' cost.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Opportunities and challenges in process modeling and simulation of electrochemical systems

Electrochemical technologies have garnered intense interest in both academic research and industrial applications due to their potential to increase energy efficiency and reduce carbon footprint. However, electrochemical process fundamentals have been absent from the chemical process simulators available to millions of chemical engineers worldwide. To expedite process research and development of electrochemical technologies in the chemical industry, it is imperative that essential electrochemical process fundamentals be incorporated into process simulators to support modeling and simulation of electrochemical processes. Here, this study examines three process fundamentals key to the research and development of electrochemical processes: electrochemical reaction kinetics, electrolyte thermodynamics, and heat and mass transfer. It further illustrates application of these process fundamentals with a case study modeling an electrochemical process for the conversion of acrylonitrile to adiponitrile. The modeling example highlights the roles of applied voltage on reaction rates, and electron flow rates on the performance of electrochemical processes. It suggests that the applied voltage and electron flow rates are two unique concepts that should be included in simulators for electrochemical processes.

09 BIOMASS FUELS↗

DOE OSTI2024

Interpenetrating 3D Electrodes for High-Rate Alkaline Water Splitting

Mass transfer is critical for the reaction kinetics and efficiency of alkaline water splitting (AWS). For AWS to operate at high current densities (hundreds of mA/cm 2 ), the device architecture must ensure a large catalytic surface area, rapid ion diffusion, and minimal solution and charge transfer resistances. Effective electrodes should also facilitate gas bubble detachment and release. 3D-printed electrodes have shown promise, but stacking them increases the ion diffusion length and solution resistance. Here we demonstrate a new device architecture with interpenetrating gyroid electrodes, providing a large ion-accessible surface area and gas diffusion channels. This design significantly reduces the interelectrode distance, lowering ion diffusion length and solution resistance. Simulations show faster ion diffusion and higher current density in the interpenetrating configuration compared with separate electrodes. This improved performance, especially at low temperatures and high current densities, highlights a promising strategy for enhancing AWS and other electrochemical systems limited by slow ion diffusion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Bridging the Gap in Carbon Free Iron Making: How Hydrogen Affects the Reduction of Iron Ore between 900 and 1590 °C

Hydrogen-based reduction of iron ore for iron and steel production has emerged as a promising alternative to coal and natural gas. Unlike other hydrogen-based iron ore reduction studies, this research focuses on a wide temperature range across 900–1590 °C, encompassing reduction in solid, mixed, and liquid (slag) phases. For a 20 min exposure to hydrogen, the reduction degree increased monotonically from ∼35% at 900 °C to >90% at 1550 °C, except between 1100 °C and 1400 °C, where it stagnated ∼60%. This experimental work challenges the widely accepted notion that higher temperatures enhance the reduction process. Instead, it reveals an overlooked kinetic bottleneck, suggesting complex thermodynamic and mass transfer limitations influenced by phase transformations, diffusion barriers, and microstructural changes. Density functional theory-based molecular dynamics simulations indicate that oxygen diffusivity in BCC iron is 3.88 × 10 –5 cm 2 /s which is ∼5–10 times higher than that in FCC iron. This study reports that this stagnant reduction degree in the mixed solid–liquid phase is due to competition of multiple mechanisms, such as surface- and bulk-diffusion, pore collapse mechanisms, and crystallographic transitions.

08 HYDROGEN↗

DOE OSTI2025

Multiphase computational fluid dynamics modeling of reacting flows in absorption columns for carbon capture

First-principles derived computational fluid dynamics (CFD) simulations have been proposed as a fundamental tool for investigating solvent-based CO 2 absorption in packed columns due to their ability to accurately represent the underlying nonlinear, multiscale dynamics. Numerous studies have previously utilized such CFD simulations to investigate hydrodynamics of columns with structured and random packings by assessing the key hydrodynamic metrics such as the interfacial and wetted areas. While mapping such metrics for different conditions is essential to the optimization of absorption columns, it is not sufficient, as the CO 2 capture rate depends also on the coupled, nonlinear dynamics from the underlying chemical reaction kinetics, thermodynamics, and heat-transfer rates. In this work, we present detailed CFD simulation results augmented by incorporating the effects of interfacial physical mass transfer of CO 2 , heat release from chemical reaction kinetics, and thermophysical property variations from resulting temperature gradients. We demonstrate the applicability of the proposed approach in numerically assessing the performance of packed columns by evaluating key hydrodynamic quantities, CO 2 absorption rates, and temperature rise in a reference column with packings that are structurally similar to the Sulzer Mellapak™ 250.Y packing, for different solvent inflow velocities and temperatures. Predictions from simulation results are found to be consistent with the trends in experimental observations from the literature, suggesting that the predictive capabilities of the simulation framework can be leveraged to guide the future development of absorber-column designs and optimized process flowsheets.

Absorption columns↗

DOE OSTI2025

Scanning Electrochemical Microscopy for Kinetic Investigations in Viscous Deep Eutectic Solvents: Identifying Practical Approach Curves and Deviations from Electron Transfer Models

Determining heterogeneous electrochemical electron transfer (ET) kinetics in electrolytes with a wide range of physical properties is of great interest for achieving high-performance redox flow batteries. Among such electrolytes, concentrated hydrogen-bonded electrolytes (CoHBEs), including deep eutectic solvents (DESs), have recently garnered significant attention. Unfortunately, traditional Tafel analysis using macroelectrodes often encounters issues with mass transfer limitations in CoHBEs with high viscosities, thereby restricting kinetic analysis to a narrow potential window. Here, in this work, we introduce a methodology for evaluating ET kinetics in viscous DES using the scanning electrochemical microscopy (SECM). We first determined practical solutions to SECM tip positioning in ethaline DES, which yield pseudopositive feedback responses. Lattice Boltzmann method (LBM) simulations helped us rationalize the impact of the fluid and concentration fields, as well as tip geometry, tip approach velocity v, and the solvent viscosity ηs, on the shape of the approach curves. In addition to successfully recreating approach curves over a variety of conditions, we found that approaching a conductor ensured a practical point where the normalized tip response (Ni T = 2) converged at L = 0.7 within ∼10% error regardless of tip velocity. With positioning capabilities at hand, we investigated the kinetics of Fe 3+ /Fe 2+ redox couple in aqueous and the ethaline media. The experimental kinetic results were interpreted using the Butler–Volmer (BV) and Marcus–Hush–Chidsey (MHC) models. For ethaline, a nonideal kinetic behavior was observed, potentially attributed to solvent dynamics within DESs or to the interplay of chloride anions in the charge transfer process.

electrodes↗

DOE OSTI2026

Impact of Small-Alkane Solvents on Polyolefin Hydrogenolysis over a Ruthenium Catalyst

Selective catalytic hydrogenolysis of polyolefins is a promising route to convert plastic waste into valuable liquid products, such as lubricants, waxes, and surfactants. However, the high viscosity of polymer melts imposes mass transfer limitations on this reaction. Solvents can mitigate these challenges, but their effects on reaction kinetics and product selectivity remain underexplored. Here, we systematically explore the effects of small n-alkanes and cycloalkanes on the hydrogenolysis of polyethylene and polypropylene over a Ru/TiO 2 catalyst. Using kinetic measurements and isotopic labeling, we show that n-octane at high mass fractions alters the mechanism from direct hydrogenation to solvent-mediated hydrogen transfer, reducing the rate of C–C bond cleavage. Longer alkanes further inhibit reactivity due to stronger surface binding. 1,4-Dimethylcyclohexane suppresses methane formation, favoring heavier products, while decalin likely forms surface-bound aromatics that poison the catalyst. Overall, alkane solvents modulate product selectivity and reduce the yield of methane byproduct, allowing for ∼35–40% selectivity to valuable C 20 -C 30 alkane products. This work highlights the complex impact of polymer–alkane mixtures on hydrogenolysis kinetics relevant to the design of commercial-scale plastic waste valorization processes.

Catalysts↗

DOE OSTI2026

Effects of Dissolution Regimes on Flow Channelization and Solute Transport in 3D Fracture Networks: Insights From Graph‐Based Reactive Transport Modeling

We investigate how mineral dissolution reshapes flow pathways and solute transport in three‐dimensional discrete fracture networks using a computationally efficient graph‐based reactive transport model. The DFNs are inspired by field‐site observations of fractured carbonate and represent realistic connectivity and structural heterogeneity. Flow is simulated with the Reynolds equation, and dissolution follows first‐order kinetics with diffusive limitations captured through an effective mass‐transfer coefficient. By systematically varying two key dimensionless parameters, the effective Damköhler number (Da), governing reaction versus advection rates, and a transport parameter (Da), analogous to the Thiele modulus, distinct flow channelization regimes emerge: mildly channelized at low G, highly channelized at intermediate Da, and extreme wormhole formation at high Da and low G. Eulerian and Lagrangian analyses, including breakthrough curves, particle tortuosity, dispersivity, and flow channeling indicators quantitatively characterize the progression of dissolution‐driven network restructuring. Across all regimes, initial fracture heterogeneity persists. The results underscore how the interplay between this initial structure, advection, reaction, and diffusion critically shapes subsurface flow pathways, with implications for applications ranging from groundwater remediation to enhanced geothermal systems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗