SEARCH · Search NASA

Results for “gas turbine engine modeling”

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 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Bayesian learning of orthogonal embeddings for multi-fidelity Gaussian Processes

Uncertainty propagation in complex engineering systems often poses significant computational challenges related to modeling and quantifying probability distributions of model outputs, as those emerge as the result of various sources of uncertainty that are inherent in the system under investigation. Gaussian Processes regression (GPs) is a robust meta-modeling technique that allows for fast model prediction and exploration of response surfaces. Multi-fidelity variations of GPs further leverage information from cheap and low fidelity model simulations in order to improve their predictive performance on the high fidelity model. In order to cope with the high volume of data required to train GPs in high dimensional design spaces, a common practice is to introduce latent design variables that are typically projections of the original input space to a lower dimensional subspace, and therefore substitute the problem of learning the initial high dimensional mapping, with that of training a GP on a low dimensional space. Here in this paper, we present a Bayesian approach to identify optimal transformations that map the input points to low dimensional latent variables. The \projection" mapping consists of an orthonormal matrix that is considered a priori unknown and needs to be inferred jointly with the GP parameters, conditioned on the available training data. The proposed Bayesian inference scheme relies on a two-step iterative algorithm that samples from the marginal posteriors of the GP parameters and the projection matrix respectively, both using Markov Chain Monte Carlo (MCMC) sampling. In order to take into account the orthogonality constraints imposed on the orthonormal projection matrix, a Geodesic Monte Carlo sampling algorithm is employed, that is suitable for exploiting probability measures on manifolds. We extend the proposed framework to multi-fidelity models using GPs including the scenarios of training multiple outputs together. We validate our framework on three synthetic problems with a known lower-dimensional subspace. The benefits of our proposed framework, are illustrated on the computationally challenging aerodynamic optimization of a last-stage blade for an industrial gas turbine, where we study the effect of an 85-dimensional shape parameterization of a three-dimensional airfoil on two output quantities of interest, specifically on the aerodynamic efficiency and the degree of reaction

42 ENGINEERING↗

DOE OSTI2021

Carbon and Hydrocarbon Particle Seeding in Air-Breathing Rotating Detonation Engine

Within the power generation community, the rotating detonation engine (RDE) is only growing in popularity with its increased performance, simple mechanism, and operation. Although significant testing is underway to characterize the RDE for integration with conventional gas turbines, this entire system is still at a relatively low technology readiness level. In the midst of RDE research, there is an initiative to understand solid particle seeding effects in the detonation performance. Under investigation at the University of Central Florida is a Department of Energy (DOE) 15.24 cm (6 in.) RDE, with a solid particle seeder in parallel with its H2 and air flow lines. Previous work on this system involved carbon particle detonation; however, the tested particles were taken one step further to include more sustainable, greener hydrocarbon particles. Testing of powdered sugar, peanut flour, and cornstarch, along with previous carbon black tests have shown not only successful detonability, but a noticeable effect on the detonation wave dynamics. Side-by-side with a particle burning model being developed, an operational map can be determined for the hydrocarbon particles particularly, which can be tuned with the local flow conditions to achieve peak operability while replacing fuels with sustainable alternatives that could even be grown.

Engineering↗

DOE OSTI2024

Advanced Compressors for CO2-Based Power Cycles and Energy

Pumped thermal energy storage (PTES) is a cost-effective method for grid-scale, long-duration electrical energy storage (LDES). PTES uses a heat pump cycle to transfer thermal energy from a low temperature reservoir (LTR) to a high temperature reservoir (HTR), later using a heat engine cycle to reverse the process. A key component of the PTES system, the heat pump compressor, represents a significant driver to the cost, performance and operating characteristics of the PTES system. For grid-scale charging (>100 MW), traditional compressor scaling charts indicate that the operating conditions needed would be best served by a multistage axial compressor. While frame gas turbine compressors at these power ratings exist and operate at higher pressure ratios than needed for the CO2 PTES system, the inlet pressure and fluid density of this application exceed experience values. A conceptual design of a large-scale CO2 axial compressor was completed, including mean-line estimates of the compressor performance at full power conditions. The results of the conceptual design were used to refine the PTES cycle design, and updated operating conditions provided for further aero design optimization. A roughly one third scale of the first three stages of the 100 MW compressor was designed, built, and tested at relevant PTES charging conditions. Results of this test program are compared against high-resolution computational fluid dynamics models.

Compressor Heat↗

DOE OSTI2025

Novel Design and Fabrication of a High Frequency Transient Heat Flux Sensor for Use in an RDE

Rotating detonation engine (RDE) combustion systems have been a topic of interest in the pressure gain combustion community for their benefits over traditional gas turbine engine combustors. However, cooling requirements for these engines are significantly higher and less predictable than non-detonating engines. To understand the high-speed heat transfer dynamics inside an RDE, a novel, high-frequency heat flux gage is presented. This study aims to design a robust, single-sided sensor that can withstand the high temperature and harsh environment of an RDE for extended durations. Sensor bench testing is performed using a hot plate as a heat source, and the sensor response is compared to a finite-element analysis (FEA) model. The sensor response is then tested inside a water-cooled RDE and the wall heat flux is compared to calorimetry data.

rotating detonation engines↗

DOE OSTI2025

Novel Design and Fabrication of a High Frequency Transient Heat Flux Sensor for Use in an RDE

rotating detonation engines↗

DOE OSTI2020

Thermal power plant upgrade via a rotating detonation combustor and retrofitted turbine with optimized endwalls

Over the past decade, pressure gain combustion research has promised over 10 percentage-points of increase in power plant thermal efficiency. Alas, to realize such potential gain, one must effectively couple the turbine with the detonation combustor, whose exhaust conditions differ substantially from current state of the art gas turbines. This paper presents a modeling approach that enables a superior thermodynamic cycle with a rotating detonation combustor and a retrofitted gas turbine by including a diffuser downstream of the combustor and by contouring the turbine endwall, while preserving the airfoil geometry. We propose a multi-step optimization strategy, parametrizing the endwall geometry with a few control points, without altering the airfoil geometry, and with the stage turbine efficiency as objective function. In a first step the turbine performance is assessed with steady inlet conditions by solving the steady three-dimensional Reynolds-Averaged Navier-Stokes equations. In a second step, the inlet conditions are unsteady, as predicted from a detonation combustor and a diffuser, and three-dimensional full unsteady simulations are performed with an unsteady Reynolds-Averaged Navier-Stokes solver. By altering the vane endwall, the steady optimization yielded an efficiency increase of 12% relative to the baseline, while the unsteady optimization resulted in 21% increase compared to the datum turbine. Finally, a full engine analysis demonstrated the superiority of pressure gain combustion which included a realistic thermodynamic cycle of the combustor, the diffuser and the optimized turbine components.

42 ENGINEERING↗

DOE OSTI2021

MACHINE LEARNING BASED CHEMICAL EXPLOSIVE MODE ANALYSIS (ML-CEMA)

The software consists of a Machine Learning based Chemical Explosive Mode Analysis (ML-CEMA) tool for advanced computational flame diagnostics. CEMA, originally based on eigen-analysis of the local thermochemical system, is capable of identifying reaction fronts and limit phenomena such as auto-ignition and extinction in practical combustion systems, such as internal combustion engines and gas turbine combustors. However, the original CEMA is computationally expensive for large reaction mechanisms that are typically needed to describe fuel chemistry of practical large-hydrocarbon fuels. This novel ML-CEMA tool employs a ML technique to accelerate the eigen-analysis of the basic CEMA approach by orders of magnitude, thus making it suitable for practical fuels. In ML-CEMA,zero-dimensional (0D) reactors and one-dimensional (1D) premixed flames are first used to generate a large number of data points for neural network based ML training. The trained ML model is then used to perform CEMA prediction. This ML-CEMA tool has been demonstrated in canonical 0D and 1D configurations as well as highly-transient three-dimensional spray flames exhibiting multi-mode turbulent combustion, showing promising results. ML-CEMA, as a standalone tool, can be used for computationally-efficient diagnostics of massive datasets generated from both experiments and simulations. For example, based on spatially resolved measurements of a small set of reactive scalars(such as temperature, hydroxyl radical and formaldehyde), ML-CEMA can effectively identify flame fronts and rare events. ML-CEMA also provides a robust online or offline flame feature detection tool. When used for on-the-fly simulations, ML-CEMA further enables zone-adaptive combustion modeling, in which the predicted eigenvalue is used as a robust mode indicator for judicious assignment of locally-valid combustion models. This ML-CEMA based zone-adaptive model can lead to substantial computational cost savings when used for large-scale simulations of multi-mode combustion systems. Third Party Code Web Page to Download Code Web Page Location of Third Party License

Xu, Chao↗

DOE OSTI2021

Role of Inlet Boundary Conditions on Fuel-Air Mixing at Supercritical Conditions

Advanced engine design and alternative fuels present the possibility of fuel injection at purely supercritical conditions in diesel engines and gas turbines. The complex interactions that govern this phenomenon still need significant research, particularly the boundary conditions for fuel injection are critical for accurate simulation. However, the flow inside the injector itself is often omitted to reduce the computational efforts, and thus, velocity, mass flux, or total pressure is specified at the injector exit (or domain inlet), often with simplified velocity profiles and turbulence levels. This simplified inlet boundary treatment has minimal effects on results for conventional fuel injection conditions, however, the validity of this approach at supercritical conditions has not been assessed. Comprehensive real-gas and binary fluid mixing models have been implemented for computational fluid dynamics (CFD) analysis of fuel-air mixing at supercritical conditions. The model is verified using prior CFD results from the literature. The model is used to investigate the effects of the shape of axial velocity and mass fraction profiles at the inlet boundary with the goal to improve the comparison of predictions to experimental data. Results show that the boundary conditions have a significant effect on the predictions, and none of the cases match precisely with experimental data. Furthermore, the study reveals that the physical location of the inlet boundary might be difficult to infer correctly from the experiments and highlights the need for high-quality, repeatable measurements at supercritical conditions to support the development of relevant high-fidelity models for fuel-air mixing.

02 PETROLEUM↗

DOE OSTI2021

Impact on Cycle Efficiency of Small CHP Plants from Increasing Firing Temperature Enabled by AM of Turbine Blades

Combined Heat and Power (CHP) systems are gaining popularity because of their potential for high overall thermodynamic efficiency and increased need for distributed power generation. Most CHP systems include a gas turbine for electricity generation and a heat recovery steam generator (HRSG) for steam generation. The steam can be used to power a steam turbine for additional electricity generation or to drive rotating equipment, for space heating, for absorption chillers, etc. A CHP cycle configuration is often driven by the ‘principal’ utility for the facility that it services; steam or electricity. In either case, performance improvements in the gas turbine have the potential to increase the steam and power output from the cycle, which is a direct result of improvements in the gas turbine efficiency and power output. One possible opportunity to improve an existing gas turbine’s performance is to increase the firing temperature with improved turbine cooling and increased compression ratio.<br>In this study, the impact on CHP cycle performance from increasing the turbine firing temperature by 100 °C and improving the turbine blade cooling for a 6-MW scale gas turbine is estimated using an analytical cooled gas turbine model and a steam cycle model. A sensitivity analysis was performed to understand the impact of increasing the internal cooling effectiveness, thermal barrier coating performance and blade material upgrades on gas turbine and CHP cycle efficiency. The impacts of turbine blade cooling improvements were studied for three common CHP cycle configurations identified from the literature. Various definitions for CHP cycle efficiency from the literature are used in the comparisons. The results showed that a 100 °C increase in firing temperature can increase the gas turbine efficiency by 1 percentage point without improving cooling effectiveness and add 2 additional percentage points in efficiency by using advanced turbine blades with higher internal cooling efficiency. Studied engine upgrades showed potential for increasing the CHP cycle efficiency by 3 percentage points while increasing the steam generation rate by 8%.<br>

Uysal, Selcuk Can↗

DOE OSTI2023

Surface Roughness Effects from Additive Manufacturing in High Efficiency Gas Turbine Combustion Systems

In the past, clean combustion systems were characterized by high system development costs due to growing complexity and escalating manufacturing costs from conventional manufacturing processes. As a result, high efficiency concepts are difficult to design and more difficult to be cost-effectively manufactured. In recent years, with the introduction of additive manufacturing (AM) technologies, the rapid prototyping and mass production processes of clean combustion systems are promising to be significantly simplified with significant reduction in terms of engine manufacturing cost and engine energy cost. Compared with conventional manufacturing built parts, the AM process enabled a simpler design to be adopted for the nozzle, reducing the number of required braze and weld joints from twenty-five to just five. The resulting nozzle was 25% lighter and five times more durable and contributed to a 15% reduction in fuel burn in comparison with the previous model produced. Applying these improvements to a fleet of 5,000 turbofan engines at 150 kN thrust would result in fuel cost savings of $ \$ $6B annually (at $ \$ $5/gallon Jet-A fuel price, 700 gallon/hour fuel consumption rate, and 2,300 operational hours per year), and reduce CO2 emissions by more than 25 million tons per year. These advances support core Department of Energy (DOE) missions in energy efficiency, improving productivity, and environmental sustainability. Maximizing the benefit of these new capabilities will require high prediction capability of high speed turbulent flow with wall modeled Large Eddy Simulation (LES). GE Aviation maintains that advanced simulation technology and supercomputing is required in order to provide the appropriate boundary conditions to realize the maximum potential of AM and Ceramic Matrix Composite to reduce cooling flow, a first order penalty on the Brayton cycle. The algorithm developed using the ANSYS/Fluent software would provide the foundation for entirely new avenues of research and development with potential multi-billion dollar impact to the US economy, and significant reduction in carbon based emissions across the aerospace and power generation industries. The formulation of this iWLES (integral wall model for LES) is generic and allows to capture the changes in flow dynamics that have a significant impact on the wall bounded flow characteristics, such as swirler effective area, bulk swirl number, local pressure distribution, exit velocity profile, and turbulence kinetic energy profile. A periodic channel flow with rough flat plate is simulated using LES (Wall-Adapting Local Eddy-viscosity) model to verify the implementation of iWLES in the Fluent User Defined Function (UDF). As the flow fields are highly sensitive due to surface roughness of the wall bounded flows in the combustion systems, there is significant potential to conduct further LES studies focusing on the turbulence boundary layer interaction. Accurately capturing near wall physics will elucidate the impact of rough surfaces on combustor flow and aero-thermal interactions.

42 ENGINEERING↗

DOE OSTI2022

Stresses within rare-earth doped yttria-stabilized zirconia thermal barrier coatings from in-situ synchrotron X-ray diffraction at high temperatures

There is a growing interest for smart coatings that can be integrated into turbine engines for in-situ temperature measurements or health monitoring. The addition of rare-earth dopants into standard thermal barrier ceramic top coat materials is used to obtain luminescent coatings that enable spectral measurements, for real-time temperature or health monitoring. The thermomechanical performance and durability of such novel coating compositions in extreme environments still remains to be evaluated. Consequently, the ability to manufacture sensor coatings which present suitable thermal properties needs to be demonstrated. For this study, highly luminescent erbium and europium doped yttria-stabilized zirconia and state-of-the-art yttria-stabilized zirconia coatings manufactured by air plasma spray were characterized to determine the effects of the embedded rareearth dopants on coating internal strain and stress and to quantify and compare their high temperature response using synchrotron X-ray diffraction. In-situ depth-resolved strain measurements were performed at 15 μm intervals along the depth of the coatings to evaluate materials response at key locations, specifically at layer interfaces. In-plane stress was calculated for the coatings and a finite element model was implemented to supplement the results and enable further predictions. The results show that the sensor coatings that were manufactured in this work revealed only minor variations in the strain response of sensor coatings under a typical thermal cycle and especially at temperatures closer to that of gas turbine operating conditions, compared to state-of-the-art coatings. This work demonstrates the viability of manufacturing rare-earth doped yttria-stabilized zirconia coatings that provide beneficial spectroscopic monitoring capabilities while having minimal impact on the thermomechanical response of the thermal barrier coatings.

Air plasma spray↗

DOE OSTI2023

Distributed Energy Resource (DER) Reliability for Backup Electric Power Systems

Hospitals, emergency services, military bases, ports, airports, industries, commercial facilities, and others rely on backup power systems to provide electricity for their critical loads during grid outages. The purpose of this report it to provide accurate reliability information on commonly deployed distributed energy resources (DERs) to improve quantitative estimates for the reliability of these backup power systems during a grid outage. A backup power system consists of DERs, an electric distribution system with its associated switches and other devices, and mechanisms to control and manage the flow of electricity. Too often, facilities and campuses fail to properly quantify the reliability of their backup power systems. DERs are assumed to be 100% reliable, with the only concern being the availability of fuel. Such assumptions can lead to gross errors in the backup system's reliability estimates, particularly for long-duration outages. This report provides a set of estimates for reliability of emergency diesel generators (EDGs), natural gas prime generators and combined heat and power (CHP) prime movers, solar photovoltaics (PV), wind turbines, and Li-ion battery energy storage systems (BESS). The estimates are derived from empirical data when available and supplemented by modeling results when needed. These reliability estimates are for the DERs ability to provide power during a grid outage, ranging from an hour to 2 weeks.

14 SOLAR ENERGY↗

DOE OSTI2021

Data Driven Fault Detection of Premixer Centerbody Degradation in a Swirl Combustor

This paper introduces a data-driven framework for combustor-focused, performance-based condition monitoring of gas turbines. Commercial condition monitoring systems typically generate huge amounts of data that make efficient onboard monitoring challenging. This paper focuses on quantifying combustor component degradation, using premixer centerbody degradation in a swirl stabilized combustor as a case study. The input for these analyses is acoustic pressure measurements acquired at various locations on the combustor. The diagnosis methodology is based on a classification framework and consists of 3 steps: 1) Data curation, 2) Feature Engineering, and 3) Diagnosis. Data curation ensures good quality of the data that is passed through the algorithm. Feature engineering deals with the extraction of the most informative features, from the most informative sensors, that can accurately capture the introduced fault. To perform diagnosis, the classification model is trained using experimentally acquired data and is then tested on a separate data set. The framework was able to achieve high classification accuracy (>99%) for training size as low as 30% of the total recorded observations. The low number of features required to achieve this accuracy suggests high potential for integration into existing onboard condition monitoring systems.

data driven methods, fault detection, swirl flames↗

DOE OSTI2023

FWP FEAA149: “Next Generation Environmental Barrier Coatings”

Environmental barrier coatings (EBCs) are required coatings for utilization of SiC/SiC ceramic matrix composite (CMC) components in gas turbines, where the EBC represents the life-limiting factor for such components. EBC/CMC systems have shown success in aero-engine applications with increased turbine inlet temperatures and improved efficiencies, which are achieved through higher temperature stability, lower density, and decreased reliance on cooling air compared to traditional superalloys. While industrial gas turbines (IGTs) do not currently utilize SiC/SiC CMCs, the current shift towards low-carbon or carbon-free fuel sources for power generation could result in a need for EBC/CMC components with higher temperature capabilities. In this work, three tasks were outlined to improve understanding of EBC lifetimes to encourage use in IGTs with carbon-free fuel such as hydrogen: 1. Define the bond coating oxidation kinetics and EBC failure criteria, 2. Measure thermal expansion coefficients of each layered material, and 3. Perform advanced characterization and modeling to assess EBC lifetimes. Cyclic steam oxidation tests were conducted on various EBC/Si/SiC chemistries and EBC/SiC architectures to define substrate oxidation kinetics and EBC failure modes. An open-source code was developed to quantify the undulating thermally grown oxide thickness with thousands of measurements from specimen cross-section images. Bond coating oxidation kinetics were determined and used to develop a kinetic and thermodynamic model for predicting EBC lifetimes. High-temperature Raman spectroscopy was utilized for determining the SiO 2 thermally grown oxide phase transformation as the life-limiting feature for EBCs. Model efforts supported the claim that the SiO 2 phase transformation causes elevated stress during thermal cycling with associated cracking that decreases the adhesion strength of the EBC, eventually resulting in coating spallation. The finite element model subroutine will be made publicly available upon internal review. Further development of an EBC lifetime model for IGTs involves definition of a critical SiO 2 thickness for EBC spallation and must also consider both environmental (gas velocity, pressure, etc.) and specimen (EBC dopants, layer architectures, etc.) effects into predicted bond coating oxidation kinetics for long-lifetime components.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Emissions and Flame Structure Assessment of Aeroengine Micromixing Injectors for Lean Direct Injection of Hydrogen and Hydrogen/Natural Gas Blends

Abstract The present study focuses on whether adapting low emissions liquid fuel injection concepts can be applied to gaseous mixtures. The goal is to develop an injector with specific geometric characteristics that will achieve the lowest emissions. This is conducted by evaluating emissions and flame diagnostics on a set of test hardware that has been developed for combustion of pure hydrogen and hydrogen/methane mixtures. The current test hardware adapts lean direct injection concepts from aeroengine injectors to a ground-based engine based on target flows associated with a Solar Turbines Industrial Engine gas turbine. A statistical test plan was developed to assess these geometric characteristics using a three level Box Behnken design. The test plan investigates the effect of injector flow splits, air swirl, fuel swirl, pressure drop, preheat temperature, fuel composition, and the adiabatic flame temperature on the flame structure and emissions performance. A total of 16 injectors were considered for a total of 688 test points for emission and imaging data. Emissions data were reported on a mass basis in ng/J and chemical reactor network analysis were conducted to further interpret emissions behavior. Images were also captured to analyze the flame structure features of interest including reaction volume, flame length, flame spread, and flame width. With these results, a predictive model for emissions was made to conclude the most optimal injector configuration is to have the highest level of air split (more volume of air through the injector plate versus the injector) and air swirl angle for both methane and hydrogen gas mixtures. The predictive emissions level for all of these cases were below 10 ng/J.

08 HYDROGEN↗

DOE OSTI2025

Phase-field predictions of the influence of cooling rates during AM on the Evolution of Microstructures in Nickel-Based Single Crystal Superalloys

Additive manufacturing of single crystals made of Ni-based superalloys offers major cost savings for gas turbine engines with the inclusion of internal cooling channels. However, the lack of understanding of the effect of transient thermal conditions on solidification grain structure during additive manufacturing hinders the potential for process control to maintain the single crystal quality. The use of high-fidelity simulations through high performance computing to predict the evolution of the solidification microstructure will enhance the abilities to tailor the microstructures through process optimization. Phase field simulations are used to determine the effect local thermal conditions and defects on the stability of the solidification morphology, specifically with respect to the onset of columnar-to-equiaxed transition that results in the loss of the single crystal. The results are expected to be instrumental for developing future surrogate models to speed up the integration of design and manufacturing of turbine blades under the harsh in-service conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Large Eddy Simulation of Externally Induced Ingress about an Axial Seal by Stator Vanes

Turbine inlet temperatures in advanced gas turbines could be as high as 2000 °C. To prevent ingress of this hot gas into the wheelspace between the stator and rotor disks, whose metals can only handle temperatures up to 850 °C, rim seals and sealing flows are used. This study examines the abilities of large eddy simulation (LES) based on the WALE subgrid model and Reynolds-averaged Navier–Stokes (RANS) based on the SST model in predicting ingress in a rotor–stator configuration with vanes but no blades, a configuration with experimental data for validation. Results were obtained for an operating condition, where the ratio of the external Reynolds number to the rotational Reynolds number is 0.538. At this operating condition, both LES and RANS were found to correctly predict the coefficient of pressure, C p , located downstream of the vanes and upstream of the seal, but only LES was able to correctly predict the sealing effectiveness. This shows C p by itself is inadequate in quantifying externally induced ingress. RANS was unable to predict the sealing effectiveness because it significantly under predicted the pressure drop in the hot gas path along the axial direction, especially about the seal region. This affected the pressure difference across the seal in the radial direction, which ultimately drives ingress.

42 ENGINEERING↗

DOE OSTI2025

Simple Heat Transfer Model for Film Cooling Applications

This report describes the development of a simple engineering model for film cooling. This model is used to derive a relationship between local wall temperature variations and key cooling performance parameters like local heat transfer coefficients and film effectiveness. This relation and method new and different from previously published models. The scope of this report includes the derivation of regression model equations for a flat plate with and without film cooling. The model equation for a flat plate without film cooling can be used to estimate local heat transfer coefficients using surface temperatures measured from infrared thermography. The model equation for the flat plate with film cooling can be used to estimate film cooling effectiveness, $η_f$, and heat transfer augmentation from the film cooling jet(s).

42 ENGINEERING↗