Search NASA⌕ Search

SEARCH · Search NASA

Results for “Latent control”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Latent-space time evolution of non-intrusive reduced-order models using Gaussian process emulation

Non-intrusive reduced-order models (ROMs) have recently generated considerable interest for constructing computationally efficient counterparts of nonlinear dynamical systems emerging from various domain sciences. They provide a low-dimensional emulation framework for systems that may be intrinsically high-dimensional. This is accomplished by utilizing a construction algorithm that is purely data-driven. It is no surprise, therefore, that the algorithmic advances of machine learning have led to non-intrusive ROMs with greater accuracy and computational gains. However, in bypassing the utilization of an equation-based evolution, it is often seen that the interpretability of the ROM framework suffers. This becomes more problematic when black-box deep learning methods are used which are notorious for lacking robustness outside the physical regime of the observed data. In this article, we propose the use of a novel latent-space interpolation algorithm based on Gaussian process regression. Notably, this reduced-order evolution of the system is parameterized by control parameters to allow for interpolation in space. The use of this procedure also allows for a continuous interpretation of time which allows for temporal interpolation. The latter aspect provides information, with quantified uncertainty, about full-state evolution at a finer resolution than that utilized for training the ROMs. This research assesses the viability of this algorithm for an advection-dominated system given by the inviscid shallow water equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Tailoring the Interfacial Interactions of Porous Polymer Membranes to Accelerate Atomic Layer Deposition: The Latent Path to Antifouling Membranes

Atomic layer deposition (ALD) is a powerful strategy to engineer hybrid organic-inorganic membranes with emergent functionalities. The combination of atomic-level thickness control, a wide materials palette, and unprecedented conformality allows the physiochemical properties (e.g., hydrophilicity) of mesoporous polymer membranes to be precisely tuned. The nucleation of ALD material growth on polymer surfaces relies on Lewis acid-base interactions and remains an overlooked motif with tremendous potential to accelerate ALD nucleation and growth. Strategies to enhance these interactions could enable desirable properties such as antifouling behavior to be imparted on inert polymer surfaces that lack the necessary functional groups for ALD nucleation. Here, in this study, we demonstrate that the reactivity of polyacrylonitrile (PAN) membranes toward ALD metal oxide (MO) precursors with Lewis acid characteristics is enhanced by introducing strong Lewis base functional groups (amidoxime: Am) on the PAN backbone (Am-PAN). The resulting Lewis acid-base interactions accelerate the MO nucleation in Am-PAN and reduce the number of deposition cycles required to achieve hydrophilicity compared with the untreated PAN membrane. Unveiling the reaction mechanism, in situ Fourier transform infrared (FTIR) spectroscopy measurements established enhanced interaction dynamics between the ALD MO precursors and the Am-PAN membrane, unlike the PAN membrane. Spectroscopic ellipsometry and thermogravimetric analysis measurements revealed higher MO loadings in Am-PAN membranes compared to PAN membranes for the same number of ALD cycles. We found that strong Lewis acid-base interactions accelerated the ALD for a range of materials including Al 2 O 3 , TiO 2 , SnO 2 , and ZnO. More broadly, our work demonstrates that tailoring metal-precursor-polymer interactions is a powerful strategy to accelerate and modulate the ALD. We used this design strategy to fabricate Al 2 O 3 -Am-PAN hybrid membranes that showed 2-fold higher antifouling capability compared to pristine PAN membranes prepared with an equivalent number of Al 2 O 3 ALD cycles. Our approach expands the scope of design options for fouling-resistant porous hybrid inorganic-organic membranes and may ultimately reduce the operational costs of water treatment.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Learning and Fast Adaptation for Grid Emergency Control via Deep Meta Reinforcement Learning

As power systems are undergoing a significant transformation with more uncertainties, less inertia and closer to operation limits, there is increasing risk of large outages. Thus, there is an imperative need to enhance grid emergency control to maintain system reliability and security. Towards this end, great progress has been made in developing deep reinforcement learning (DRL) based grid control solutions in recent years. However, existing DRL-based solutions have two main limitations: 1) they cannot handle well with a wide range of grid operation conditions, system parameters, and contingencies; 2) they generally lack the ability to fast adapt to new grid operation conditions, system parameters, and contingencies, limiting their applicability for real-world applications. Here, in this paper, we mitigate these limitations by developing a novel deep meta-reinforcement learning (DMRL) algorithm. The DMRL combines the meta strategy optimization together with DRL, and trains policies modulated by a latent space that can quickly adapt to new scenarios. We test the developed DMRL algorithm on the IEEE 300-bus system. We demonstrate fast adaptation of the meta-trained DRL polices with latent variables to new operating conditions and scenarios using the proposed method, which achieves superior performance compared to the state-of-the-art DRL and model predictive control (MPC) methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

On the predictability of turbulent fluxes from land: PLUMBER2 MIP experimental description and preliminary results

Accurate representation of the turbulent exchange of carbon, water, and heat between the land surface and the atmosphere is critical for modelling global energy, water, and carbon cycles in both future climate projections and weather forecasts. Evaluation of models' ability to do this is performed in a wide range of simulation environments, often without explicit consideration of the degree of observational constraint or uncertainty and typically without quantification of benchmark performance expectations. We describe a Model Intercomparison Project (MIP) that attempts to resolve these shortcomings, comparing the surface turbulent heat flux predictions of around 20 different land models provided with in situ meteorological forcing evaluated with measured surface fluxes using quality-controlled data from 170 eddy-covariance-based flux tower sites. Predictions from seven out-of-sample empirical models are used to quantify the information available to land models in their forcing data and so the potential for land model performance improvement. Sites with unusual behaviour, complicated processes, poor data quality, or uncommon flux magnitude are more difficult to predict for both mechanistic and empirical models, providing a means of fairer assessment of land model performance. When examining observational uncertainty, model performance does not appear to improve in low-turbulence periods or with energy-balance-corrected flux tower data, and indeed some results raise questions about whether the energy balance correction process itself is appropriate. In all cases the results are broadly consistent, with simple out-of-sample empirical models, including linear regression, comfortably outperforming mechanistic land models. In all but two cases, latent heat flux and net ecosystem exchange of CO 2 are better predicted by land models than sensible heat flux, despite it seeming to have fewer physical controlling processes. Land models that are implemented in Earth system models also appear to perform notably better than stand-alone ecosystem (including demographic) models, at least in terms of the fluxes examined here. The approach we outline enables isolation of the locations and conditions under which model developers can know that a land model can improve, allowing information pathways and discrete parameterisations in models to be identified and targeted for future model development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

On the validity and limitations of 1D model for heat and mass transfer performance evaluation in a multilayer binder-free desiccant dehumidifier: isothermal dehumidification with internal cooling

Efficient humidity control is essential for maintaining indoor thermal comfort, yet conventional vapor-compression-based dehumidifiers are energy-intensive. Employing separate sensible and latent cooling through desiccant-coated heat exchangers (DCHEs) combined with evaporative coolers offers energy savings of up to 80 % compared to conventional systems. However, the dehumidification performance of DCHEs remains limited due to the use of polymer binders for coating desiccant materials onto heat exchange surfaces. In our previous study, we developed a multilayer fixed-bed binder-free desiccant dehumidifier (MFBDD) that demonstrated high dehumidification capacity and low pressure drop compared to rotary desiccant wheels. Nevertheless, its potential for further enhancement through internal cooling and the use of step-shaped adsorption isotherms has not been explored. In this study, a physics-based one-dimensional (1D) transient model is developed and validated to capture the coupled heat and mass transfer processes in the MFBDD and extended to simulate internal cooling using a high-capacity composite metal–organic framework, MIL-101/GO-6 (water uptake ≈1.6 g/g within 35–47 % RH). The model enables detailed analysis of local air and bed temperature dynamics and quantifies how internal cooling affects the dehumidification performance under a wide range of operating conditions. Results show that integrating internal cooling and using MIL-101/GO-6 enhance mass adsorbed, moisture removal capacity, and dehumidification effectiveness by 50 %–99 % compared with the M.S. Gel baseline. The study further reveals that achieving near-isothermal operation requires simultaneous enhancement of the convective heat transfer coefficient and heat exchange surface area. In conclusion, this work provides the first detailed physical insight into the interplay between internal cooling and step-shaped isotherms in a binder-free desiccant device and establishes a validated modeling framework for scaling up and system-level performance evaluation of next-generation energy-efficient dehumidification systems.

Heat and mass transfer↗
DOE OSTI2025

Meteorological Variables and Energy Fluxes at the Pumphouse Site, Crested Butte, CO 2017-2019

This data contains output from the pumphouse eddy covariance tower that includes shortwave radiation, longwave radiation, net radiation, air temperature, relative humidity, as well as sensible, latent, and ground heat fluxes. Also included is calculated evapotranspiration from the latent heat flux and the latent heat of vaporization. All data are on a daily timestep and displayed in Mountain Time. The data has been processed, and Quality Assurance / Quality Control (QA/QC) was done, but any daily gaps in the data have not been filled in. This research was funded by the Department of Energy and performed as part of the Watershed Function Scientific Focus Area. This research aimed to constrain evapotranspiration in a high-elevation catchment.The dataset includes one comma-separated values (CSV) data file (EddyCovariance_MeteorlogicalVariables_CrestedButtePumphouse.csv). Additionally, three metadata CSV files are included: (1) location metadata file (locations.csv), which contains location metadata and coordinates; (2) a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and (3) a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type.This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Novel Representation of Leaf Phenology Improves Simulation of Amazonian Evergreen Forest Photosynthesis in a Land Surface Model

Leaf phenology in the humid tropics largely regulates the seasonality of forest carbon and water exchange. However, it is inadequately represented in most global land surface models due to limited understanding of its controls. Based on intensive field studies at four Amazonian evergreen forests, we propose a novel, quantitative representation of tropical forest leaf phenology, which links multiple environmental variables with the seasonality of new leaf production and old leaf litterfall. The new phenology simulates higher rates of leaf turnover (new leaves replacing old leaves) in dry seasons with more sunlight, which is then implemented in ORCHIDEE, together with recent findings of ontogeny-associated photosynthetic capacity, and is evaluated against ground-based measurements of leaf phenology (canopy leaf area index and litterfall), eddy covariance fluxes (photosynthesis and latent heat), and carbon allocations from field observations. Results show the periodical cycles of solar radiation and vapor pressure deficit are the two most important environmental variables that are empirically related to new leaf production and old leaf abscission in tropical evergreen forests. The model with new representation of leaf phenology captures the seasonality of canopy photosynthesis at three out of four sites, as well as the seasonality of litterfall, latent heat, and light use efficiency of photosynthesis at all tested sites, and improves the seasonality of carbon allocations to leaves, roots, and sapwoods. This study advances understanding of the environmental controls on tropical leaf phenology and offers an improved modeling tool for gridded simulations of interannual CO 2 and water fluxes in the tropics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Modular Hydronic Room Conditioning System (CRADA NFE-24-10120 Final Report)

This project presents the development and validation of high-fidelity numerical models for fin-tube heat exchangers to enable accurate performance prediction and informed design optimization. The modeling framework integrates detailed geometric specifications, thermophysical property data, and system-level constraints to simulate the heat exchanger’s behavior under a range of operating conditions. The model is calibrated using real-world product specifications and validated against experimental data collected from controlled cooling and heating tests. In cooling mode, the model captures the overall trends in capacity and outlet air temperature but tends to underpredict latent effects, especially at lower air flow rates. In heating mode, the simulation consistently overestimates both the thermal capacity and outlet air temperature, indicating the need for refinement in air-side heat transfer assumptions. Despite these deviations, the model provides a solid foundation for optimization, allowing key design variables to be tuned within physical and performance-based constraints. This research advances the ability to simulate, validate, and optimize fin-tube heat exchanger designs with greater confidence and efficiency.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Modular Hydronic Room Conditioning System

This project presents the development and validation of high-fidelity numerical models for fin-tube heat exchangers to enable accurate performance prediction and informed design optimization. The modeling framework integrates detailed geometric specifications, thermophysical property data, and system-level constraints to simulate the heat exchanger’s behavior under a range of operating conditions. The model is calibrated using real-world product specifications and validated against experimental data collected from controlled cooling and heating tests. In cooling mode, the model captures the overall trends in capacity and outlet air temperature but tends to underpredict latent effects, especially at lower air flow rates. In heating mode, the simulation consistently overestimates both the thermal capacity and outlet air temperature, indicating the need for refinement in air-side heat transfer assumptions. Despite these deviations, the model provides a solid foundation for optimization, allowing key design variables to be tuned within physical and performance-based constraints. This research advances the ability to simulate, validate, and optimize fin tube heat exchanger designs with greater confidence and efficiency.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Latent Twins

Over the past decade, scientific machine learning has transformed the development of mathematical and computational frameworks for analyzing, modeling, and predicting complex systems. From inverse problems to numerical partial differential equations (PDEs), dynamical systems, and model reduction, these advances have pushed the boundaries of what can be simulated. Yet they have often progressed in parallel, with representation learning and algorithmic solution methods evolving largely as separate pipelines. With Latent Twins, we propose a unifying mathematical framework that creates a hidden surrogate in latent space for the underlying equations. Whereas digital twins mirror physical systems in the digital world, Latent Twins mirror mathematical systems in a learned latent space governed by operators. Through this lens, classical modeling, inversion, model reduction, and operator approximation all emerge as special cases of a single principle. We establish the fundamental approximation properties of Latent Twins for both ordinary differential equations (ODEs) and PDEs and demonstrate the framework across three representative settings: (i) canonical ODEs, capturing diverse dynamical regimes; (ii) a PDE benchmark using the shallow-water equations, contrasting Latent Twin simulations with deep operator network and forecasts with a four-dimensional variational method baseline; and (iii) a challenging real-data geopotential reanalysis dataset, reconstructing and forecasting from sparse, noisy observations. Latent Twins provide a compact, interpretable surrogate for solution operators that evaluate across arbitrary time gaps in a single-shot, while remaining compatible with scientific pipelines such as assimilation, control, and uncertainty quantification. Looking forward, this framework offers scalable, theory-grounded surrogates that bridge data-driven representation learning and classical scientific modeling across disciplines.

Latent Twins↗
DOE OSTI2022

Optimization of WAAM Process to Produce AUSC Components with Increased Service Life

Additive manufacturing has the potential to revolutionize industrial hardware and unlock efficiency gains through the fabrication of geometries and architectures not possible by conventional processing. Wire Arc Additive manufacturing (WAAM) process is a class of directed energy deposition process enabling higher build rate and allowing custom wire feedstock allowing spatial variation of microstructure. The larger spot size and lower speed creates a larger melt pool, reducing residual stress and often time creating directional/columnar microstructure for Ni-superalloy. The use of flexible platform provides freedom in deposition strategy, which can accommodate complex substrates, including feature addition onto existing structures, non-flat layers, and repair methods. However, the certification of the final component needs to match the strength requirement. Hence, the quality of the component produced is stringently monitored to avoid buildup of residual stress, cracks, porosity and to reduce detrimental segregated phases commonly observed during alloy solidification. Simulation plays a huge role in predicting the melt pool dimension and can be used to optimize the process parameter. Similar development is also required to perform physics-based modeling of microstructural development in WAAM that can predict the microstructural features during solidification and can be used to optimize the process more effectively. To move toward this goal, Raytheon Technologies Research Center together with Siemens worked to create a set of computational tools to control the process parameters, enable on-line measurements and acquisition with feedback to the optimized process parameters, and eventually track material evolution through each step of the additive process. Computational fluid dynamics is used for accurate prediction and calibration of the thermal field during WAAM process. and phase field models for microstructure evolution as a function of processing parameters to establish a connection between additive parameters and the final microstructure. Here we report cellular automata (CA) model development to predict the dendritic microstructure evolution with surface and bulk nuclei for single track and multiple layers. The CA model was developed to account for secondary element addition and predict segregation, local melting, and latent heat release as well as prediction of Euler angles from orientation information and validated against experiments. This framework was utilized to tailor spatially-varying composition in a part by appropriately controlling the microstructure evolution during the additive process. Functionally graded Haynes 282 alloy with high Cr content at the surface was tested for oxidation and mechanical properties. An advanced physics-based reaction-diffusion model predicting the simultaneous creation of chromium oxide and alumina is developed and validated to extend life expectancy of the WAAM manufactured high temperature part. A machine-learning data-driven framework establishing the process-structure relationship from a dataset of real microstructure images and corresponding process history data has been developed and implemented in the NX Siemens design system. The digital twin configuration along with the tool path generation enabled prediction of WAAM component buildup time and the techno-economic analysis provided a favorable option for all 4 cases with 15-40% cost reduction.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Dynamic Thermal Performance Analysis of PCM Products Used for Energy Efficiency and Internal Climate Control in Buildings

PCMs are attractive for the future generation of buildings, where energy efficiency targets and thermal comfort expectations are increasingly prioritized. Experimental analysis of local thermal processes in these dynamic components and whole-building energy consumption predictions are essential for the proper implementation of PCMs in buildings. This paper discusses the experimental analysis of the thermophysical characteristics of both a latent heat storage material (PCM) and a product containing this PCM. The prototype product under investigation is a panelized PCM technology containing inorganic, salt-hydrate-based PCM. The thermal analysis includes studies of melting and freezing temperatures, enthalpy changes during phase change processes, nucleation intensity, sub-cooling effects, and PCM stability. The PCM’s stability is also investigated, as is the ability of PCM products to control local temperatures and peak load transmission times. Two inorganic PCM formulations based on calcium chloride hexahydrate (CaCl 2 .6H 2 O) were prepared and tested in laboratory conditions. Material-scale testing results were compared with outcomes from the system-scale analysis, using both laboratory test methods as well as field exposure in test huts. This work demonstrates that PCM technologies used in buildings can effectively control both the magnitude of thermal storage capacity as well as the time of the peak thermal load. It was found that commonly used material-scale testing methods may not always be beneficial in assessing the dynamic thermal performance characteristics of building technologies containing PCMs.

42 ENGINEERING↗
DOE OSTI2023

Latent code-based fusion: A Volterra neural network approach

We propose a deep structure encoder using Volterra Neural Networks (VNNs) to seek a latent representation of multi-modal data whose features are jointly captured by a union of subspaces. The so-called self-representation embedding of the latent codes leads to a simplified fusion which is driven by a similarly constructed decoding. The Volterra Filter architecture achieved reduction in parameter complexity is primarily due to controlled non-linearities being introduced by the higher-order convolutions in lieu of generalized activation functions. Experimental results on two different datasets have shown a significant improvement in the clustering performance for VNNs auto-encoder over conventional Convolutional Neural Networks (CNNs) auto-encoder. In addition, we also show that the proposed approach demonstrates a much-improved sample complexity over CNN-based auto-encoder with a robust classification performance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Activity Time-Use Patterns and Experienced Satisfaction: Travel Behavioral Dynamics from the 2010–2013 American Time-Use Survey Data

Central to the theory of activity-based approaches is the role of time-use and its impacts on various decision processes related to individual activity-travel patterns, thereby determining ones’ daily well-being. This research aims at quantifying and investigating the underlying correlation of time-use and happiness using the repeated cross-section data from the 2010-2013 American Time Use Survey Well-Being module. The survey specifically records reported levels of activity satisfaction as measured through – happiness, meaningfulness, tiredness, sadness and stress experienced from activity participation. Using these scores, we first construct a composite degree of activity satisfaction that a person experiences using latent class cluster analysis. A series of regression models are then estimated to jointly analyze episodic duration and satisfaction separately for each of years from 2010-2013, by controlling for socioeconomics, household lifecycle stages (29-groups based on age, gender, marital status and presence of children), activity characteristics and global life and health assessment measures. Overall, model estimation results reveal a strong negative correlation between episodic duration and happiness from individual activity engagement patterns. This correlation pattern is temporally stable for the years 2010-2013, with minor variations. Results also suggest that household lifecycle stage has a major impact on activity time-use and satisfaction behavior. Further, single parent and younger couples (25 years to 45 years) with children are more likely to engage in longer and less happy episodes due to family commitments. The model findings also show a strong correlation between health conditions and activity happiness.

Ravulaparthy, Srinath↗
DOE OSTI2022

Patterned crystal growth and heat wave generation in hydrogels

Abstract The crystallization of metastable liquid phase change materials releases stored energy as latent heat upon nucleation and may therefore provide a triggerable means of activating downstream processes that respond to changes in temperature. In this work, we describe a strategy for controlling the fast, exothermic crystallization of sodium acetate from a metastable aqueous solution into trihydrate crystals within a polyacrylamide hydrogel whose polymerization state has been patterned using photomasks. A comprehensive experimental study of crystal shapes, crystal growth front velocities and evolving thermal profiles showed that rapid growth of long needle-like crystals through unpolymerized solutions produced peak temperatures of up to 45˚C, while slower-crystallizing polymerized solutions produced polycrystalline composites and peaked at 30˚C due to lower rates of heat release relative to dissipation in these regions. This temperature difference in the propagating heat waves, which we describe using a proposed analytical model, enables the use of this strategy to selectively activate thermoresponsive processes in predefined areas.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Biophysical impacts of Earth greening largely controlled by aerodynamic resistance

Satellite observations show widespread increasing trends of leaf area index (LAI), known as the Earth greening. However, the biophysical impacts of this greening on land surface temperature (LST) remain unclear. Here, we quantify the biophysical impacts of Earth greening on LST from 2000 to 2014 and disentangle the contributions of different factors using a physically based attribution model. We find that 93% of the global vegetated area shows negative sensitivity of LST to LAI increase at the annual scale, especially for semiarid woody vegetation. Further considering the LAI trends (P ≤ 0.1), 30% of the global vegetated area is cooled by these trends and 5% is warmed. Aerodynamic resistance is the dominant factor in controlling Earth greening’s biophysical impacts: The increase in LAI produces a decrease in aerodynamic resistance, thereby favoring increased turbulent heat transfer between the land and the atmosphere, especially latent heat flux.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Long–short-term memory encoder–decoder with regularized hidden dynamics for fault detection in industrial processes

The ability of recurrent neural networks (RNN) to model nonlinear dynamics of high dimensional process data has enabled data-driven RNN-based fault detection algorithms. Previous studies have focused on detecting faults by identifying the discrepancies in data distribution between the faulty and normal data, as reflected in prediction errors generated by RNN models. However, in industrial processes, variations in data distribution can also result from changes in normal control setpoints and compensatory control adjustments in response to disturbances, making it hard to differentiate between normal and faulty conditions. This paper proposes a fault detection method utilizing a long short-term memory (LSTM) encoder–decoder structure with regularized hidden dynamics and reversible instance normalization (RevIN) to compactly represent high-dimensional measurements for effective monitoring. During training, the hidden states of the model are regularized to form a low-dimensional latent space representation of the original multivariate time series data. As a result, the prediction errors of the latent states can be used to monitor the abnormal dynamic variations, while the reconstruction errors of the measured variables are used to monitor the abnormal static variations. Furthermore, the proposed indices can reflect operating conditions, even when the distribution of test data changes, which helps distinguish faults from normal adjustments and disturbances that controllers can settle. Here, data from numerical simulation and the Tennessee Eastman process are used to illustrate the effectiveness of the proposed fault detection method.

42 ENGINEERING↗