SEARCH · Search NASA

Results for “frictional force”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2022

In-situ, nanoscale fracture toughness measurements for improved mechanical interfaces

In this project, we demonstrated stable nanoscale fracture in single-crystal silicon using an in-situ wedge-loaded double cantilever beam (DCB) specimen. The fracture toughness KIC was calculated directly from instrumented measurement of force and displacement via finite element analysis with frictional corrections. Measurements on multiple test specimens were used to show K IC = 0.72 ± 0.07 MPa m 1/2 on {111} planes and observe the crack-growth resistance curve in <500 nm increments. The exquisite stability of crack growth, instrumented measurement of material response, and direct visual access to observe nanoscale fracture processes in an ideally brittle material differentiate this approach from prior DCB methods.

36 MATERIALS SCIENCE↗

DOE OSTI2022

UAE6 - Wind Tunnel Tests Data - UAE6 - Sequence E - Raw Data

Sequence E: Yaw Releases (P) This test sequence used a downwind, rigid turbine with a 3.4° cone angle. The wind speeds ranged from 7 m/s to 17 m/s. Initial yaw angles of ±90° were achieved. The blade tip pitch was 3°. The rotor rotated at 72 RPM. Blade and probe pressure measurements were collected. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to –99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. The turbine was positioned at each specified yaw angle. A 30- second data set was collected with the yaw brake engaged. These points used a letter for the last digit beginning with A and proceeding through the alphabet as repeat points were needed. These points were plotted and compared to ascertain the functionality of the instrumentation. The files ending in numbers beginning with 0 and increasing with each repetition represent yaw release points. Once the fixed-position yaw test was complete, the yaw drive was engaged to hold the turbine. The yaw brake was released, and the yaw drive was disabled allowing the turbine to yaw freely. The yaw drive was disabled about 5 seconds into the campaign, and the turbine was allowed to yaw freely for the rest of the 30-second duration. When the yaw drive is disabled, no torque is applied to the motor, and the inertia and yawing force of the nacelle and rotor overcomes the friction and inertia of the yaw gear and motor. The yaw releases were repeated five times at each condition. Four additional campaigns were collected to determine the turbine’s natural yaw error by releasing the brake at 0° yaw error. These file names use the E designation, followed by two digits for wind speed, followed by XXXX, followed by the repetition digit.

17 WIND ENERGY↗

DOE OSTI2024

Torsional force microscopy of van der Waals moirés and atomic lattices

In a stack of atomically thin van der Waals layers, introducing interlayer twist creates a moiré superlattice whose period is a function of twist angle. Changes in that twist angle of even hundredths of a degree can dramatically transform the system’s electronic properties. Setting a precise and uniform twist angle for a stack remains difficult; hence, determining that twist angle and mapping its spatial variation is very important. Techniques have emerged to do this by imaging the moiré, but most of these require sophisticated infrastructure, time-consuming sample preparation beyond stack synthesis, or both. In this work, we show that torsional force microscopy (TFM), a scanning probe technique sensitive to dynamic friction, can reveal surface and shallow subsurface structure of van der Waals stacks on multiple length scales: the moirés formed between bi-layers of graphene and between graphene and hexagonal boron nitride (hBN) and also the atomic crystal lattices of graphene and hBN. In TFM, torsional motion of an Atomic Force Microscope (AFM) cantilever is monitored as it is actively driven at a torsional resonance while a feedback loop maintains contact at a set force with the sample surface. TFM works at room temperature in air, with no need for an electrical bias between the tip and the sample, making it applicable to a wide array of samples. It should enable determination of precise structural information including twist angles and strain in moiré superlattices and crystallographic orientation of van der Waals flakes to support predictable moiré heterostructure fabrication.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2023

Rheology of concentrated fiber suspensions with a load-dependent friction coefficient

Here we investigate the effects of fiber aspect ratio, roughness, flexibility, and volume fraction on the rheology of concentrated suspensions in a steady shear flow using direct numerical simulations. We model the fibers as inextensible continuous flexible slender bodies with the Euler-Bernoulli beam equation governing their dynamics suspended in an incompressible Newtonian fluid. The fiber dynamics and fluid flow coupling is achieved using the immersed boundary method. In addition, the fiber surface roughness might lead to interfiber contacts, resulting in normal and tangential forces between the fibers, which follow Coulomb’s law of friction. The surface roughness is modeled as hemispherical protrusions on the fiber surfaces. Their deformation results in a normal load-dependent friction coefficient. Our simulations accurately predict the experimentally observed shear thinning in fiber suspensions. Furthermore, we find that the suspension viscosity η increases with increasing the volume fraction, roughness, fiber rigidity, and aspect ratio. The increase in η is the macroscopic manifestation of a similar increase in the microscopic contact contribution to the total stress with these parameters. In addition, we observe positive and negative first N 1 and second N 2 normal stress differences, respectively, with |N 2 | < |N 1 |, in agreement with previous experiments. Last, we propose a modified Maron-Pierce law to quantify the reduction in the jamming volume fraction by increasing the fiber aspect ratio and roughness. Our results and analysis establish the use of fiber surface tribology to tune the suspension flow behavior.

42 ENGINEERING↗

DOE OSTI2023

A constitutive model for sheared dense fiber suspensions

Here we propose a constitutive model to predict the viscosity of fiber suspensions, which undergoes shear thinning, at various volume fractions, aspect ratios, and shear stresses/rates. We calibrate the model using the data from direct numerical simulation and prove the accuracy by predicting experimental measurements from the literature. We use a friction coefficient decreasing with the normal load between the fibers to quantitatively reproduce the experimentally observed shear thinning in fiber suspensions. In this model, the effective normal contact force, which is directly proportional to the bulk shear stress, determines the effective friction coefficient. A rise in the shear stress reduces the effective friction coefficient in the suspension. As a result, the jamming volume fraction increases with the shear stress, resulting in a shear thinning in the suspension viscosity. Moreover, we extend the model to quantify the effects of fiber volume fraction and aspect ratio in the suspension. We calibrate this model using the data from numerical simulations for the rate-controlled shear flow. Once calibrated, we show that the model can be used to predict the relative viscosity for different volume fractions, shear stresses, and aspect ratios. The model predictions are in excellent agreement with the available experimental measurements from the literature. The findings of this study can potentially be used to tune the fiber size and volume fraction for designing the suspension rheology in various applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Stress- and Time-Dependent Formation of Self-Lubricating In Situ Carbon (SLIC) Films on Catalytically-Active Noble Alloys

Abstract Low shear strength (30 MPa) organic films were grown in situ on Pt 0.9 Au 0.1 surfaces via cyclic sliding contact in dry N 2 with trace concentrations of ambient hydrocarbons. We present a systematic investigation of the stress- and time-dependent film formation. Steady-state friction coefficients were found to be as low as µ ~ 0.015 and inversely proportional to contact pressure, revealing non-Amontonian behavior. Above a Hertzian contact pressure of ~500 MPa, shear strength dropped, indicating an activated process. Raman spectroscopy identified non-uniformity in areal coverage and relative order with contact pressure. Regions of steady-state low-friction behavior exhibited spectra similar to DLC coatings. Atomic force microscopy was used to study the formation and growth of films at the nanoscale. Stress- and time-dependent measurements suggested a sublinear increase of film volume with time, and a transition from growth to wear at a Hertzian contact pressure of ~1.2 GPa.

Jones, Morgan R.↗

DOE OSTI2025

Challenges in Continuous In-Field Critical Current Testing of High-Temperature Superconducting Tapes: Thermal and Mechanical Perspectives

High-temperature superconductors (HTS) are essential for ultra-high-field applications requiring exceptional current-carrying capacity under extreme conditions. However, systematic characterization of critical current in long-length conductors remains challenging due to complex thermal, electromag netic, and mechanical interactions during continuous testing. This study reports the development of a continuous in-field magnetization testing system for position-dependent critical current measurement in HTS tapes at 20 K under 7.5 T fields applied normal to the tape plane, enabling identification of performance-limiting regions that could compromise magnet stability. Here, the system addresses two fundamental challenges inherent to cryogenic reel to-reel testing. First, thermal management requires continuous cooling of a moving conductor to 20 K, achieved through liquid nitrogen precooling combined with a 100 W@20 K Gifford McMahon cryocooler. Second, screening currents in high fields generate Lorentz forces that induce twisting, bowing, and potential delamination. To mitigate these risks, we propose mechanical reinforcement and active current density suppression strategies. Numerical simulations using the stream function formulation reveal four primary failure modes: frictional heating at guide interfaces, unstable equilibria causing deformation, transverse current-induced stresses at guide transitions, and unsupported forces in vertical spans. Our mitigation strategies include PTFE coated guides to minimize friction, spring-loaded stabilization mechanisms to maintain tape alignment, controlled pre-heating using the liquid nitrogen thermal jacket to suppress critical current at stress points, and optimized guide positioning to minimize force accumulation. The experimental system is nearing completion, with testing planned to commence within two months. Preliminary validation at 65 K under 0.5 T demonstrates strong correlation between simulation-predicted mechanical instabilities and observed critical current variations during conductor tran sitions through the measurement region. These findings establish a robust foundation for quality assurance protocols essential to next-generation superconducting magnet applications.

Chen, Siwei [Princeton Plasma Physics Laboratory (↗

DOE OSTI2026

Towards High-Speed friction stir welding of 25 mm Thick AA2139-T8: tool innovation and process development

Joining thick plates (≥ 12 mm) of aluminum (Al) alloys is challenging due to high tool forces, uneven material flow, and non-uniform heat distribution through the material’s thickness. Recent advancements in tool design and welding parameters for friction stir welding (FSW) in thick plates butt joining have encountered a developmental plateau, highlighting the need for innovative approaches to overcome existing limitations. Here, this study focuses on the systematic development of single-pass and double-pass FSW processes for high-strength aluminum alloy AA2139-T8, to improve joint efficiency and enable high-speed welding capabilities. Experimental evidence is presented for a novel tool design with opposing pin threads, enabling high-speed (178 mm/min) single-pass friction stir butt welding of 25 mm thick AA2139-T8. A series of FSW trials was conducted both in air and with a trailing water spray, using steel backing plates (BPs) to investigate the impact of quenching and cooling rates on process response and joint performance. A joint efficiency of 83% was attained using the novel tool features and effectively controlling process forces and thermal boundary conditions.

AA2139-T8↗

DOE OSTI2021

Microstructural evolution in Cu-Nb processed via friction consolidation

Immiscible alloys, whether in well-mixed or layered forms, are of increasing interest based on their novel structural and functional properties, such as enhanced thermal stability against grain growth or radiation-induced defect trapping at the interfaces. To address the need for new approaches to tailor microstructures, the microstructural development of an immiscible Cu-Nb alloy processed via friction consolidation of elemental powders is investigated. Friction consolidation is a solid-phase processing technique that imparts severe plastic strain into a deforming volume resulting in elevated temperatures below the melting temperature of the alloy. Two distinct processing pathways were chosen to understand the effect of thermomechanical conditions on the final microstructure. The microstructure was characterized using scanning electron microscopy, scanning transmission electron microscopy, and X-ray diffraction techniques. Path 1 exhibited larger strain, strain rate, and temperature as compared with path 2. In path 1, agglomerated Nb particles were present in the recrystallized ultrafine-grained Cu matrix, while in path 2 extremely fine and dispersed Nb particles were present in a highly deformed Cu matrix. In both pathways, supersaturation of Cu in Nb lattices was noted, but not vice-versa. The asymmetry in mixing is explained based on deformation-based, thermodynamic and kinetic factors. These findings provide a pathway for creation of novel tailored microstructures and improved properties in any number of binary immiscible alloy systems.

Cu-Nb alloys, Microstructural refinement, Friction↗

DOE OSTI2026

C‐Cracking of Brittle‐Material Spheres From Eccentric Hertzian Contact

Contact of brittle-material spheres can occur while they are manufactured, handled, or used in operation as roller elements or in deformably processed composite material. Statistically, the directional contact between two spheres will nearly always be oblique or eccentric in an unconfined and mobile sphere population. Although analyses of oblique contact exist, the portrayal of its maximum (tensile) first principal stress ( S 1 ) field is lacking. This information is needed for improved judgment of the prospect of crack initiation. Given these, the S 1 due to eccentric contact and subsequent crack initiation was analyzed using finite element analysis (FEA) and corroborative demonstration of c-crack creation. The FEA shows that an asymmetric S 1 field is created about the contact patch and is caused by superimposed shear intrinsic to eccentric contact. A crack will initiate if the maximum S 1 exceeds the material's tensile strength, and its field will cause the crack to propagate and arrest to form a C-shaped crack, or c-crack. C-crack initiation trends with lower applied forces, with greater amount of contact eccentricity, higher coefficient of friction, and smaller spherical radius. Finally, these metrics deserve recognition, especially when c-cracked brittle-material spheres are used in a system whose functionality or reliability is potentially compromised by the c-cracks.

Hertzian contact↗

DOE OSTI2024

Optimizing Deep Geothermal Drilling for Energy Sustainability in the Appalachian Basin

This study investigates the geological and geomechanical characteristics of the MIP 1S geothermal well in the Appalachian Basin to optimize drilling and address the wellbore stability issues encountered. Data from well logs, sidewall core analysis, and injection tests were used to derive elastic and rock strength properties, as well as stress and pore pressure profiles. A robust 1D-geomechanical model was developed and validated, correlating strongly with wellbore instability observations. This revealed significant wellbore breakout, widening the diameter from 12 ¼ inches to over 16 inches. Advanced technologies like Cerebro Force™ In-Bit Sensing were used to monitor drilling performance with high accuracy. This technology tracks critical metrics such as bit acceleration, vibration in the x, y, and z directions, Gyro RPM, stick-slip indicators, and bending on the bit. Cerebro Force™ readings identified hole drag caused by poor hole conditions, including friction between the drill string and wellbore walls and the presence of cuttings or debris. This led to higher torque and weight on bit (WOB) readings at the surface compared to downhole measurements, affecting drilling efficiency and wellbore stability. Optimal drilling parameters for future deep geothermal wells were determined based on these findings.

Environmental Sciences & Ecology↗

DOE OSTI2020

Literature Review on the Theory and Measurement of Dry Sliding Friction of Metals

Understanding the friction behavior between two sliding bodies can inform the design of machines, processing of materials, and simulation of dynamic processes. Kinetic friction is a complex phenomenon that depends on a multitude of factors such as sliding velocity, normal force, contact area, surface roughness, material properties, lubrication conditions, and thermal effects. This literature review covers the major known effects of sliding velocity, normal load, and surface roughness on the measured kinetic friction coefficient in the context of microscopic friction phenomena. Classic macroscale friction models are reviewed to illustrate approaches for simulating friction behavior. Prominent experimental friction setups within the literature are discussed with respect to achievable velocity and pressure regimes. The background information gathered here will be used to inform experimental procedures and modeling strategies of the exploratory research (ER) project titled “Measurement of Dynamic Friction via Kolsky Bar” (20200418ER).

36 MATERIALS SCIENCE↗

DOE OSTI2025

Influence of Initial fabric and water-wetting on particle fracture and force chain evolution in natural silica sand

Particle fracture has a significant influence on the engineering behavior of granular materials. However, the combined influence of the initial fabric and wetting on particle fracture and force chain evolution in granular materials, particularly in sands, remains insufficiently explored in the current literature. Here, in this paper, one-dimensional (1D) confined compression experiments were conducted on dry and wet specimens, and particle-level fracture was captured using three-dimensional (3D) in-situ Synchrotron Micro-Computed Tomography (SMT). A quantitative assessment of particle fracture revealed that wet specimens, which have disturbed fabric, exhibited a higher fracture percentage in comparison to dry specimens. Complementary 3D finite element (FE) simulations were performed using dry and wet (only solid fabric was considered and the influence of pore water is not considered) sand assemblies to assess interparticle contact forces and particle-scale stress distributions within the specimens. Force chain analysis was subsequently conducted, encompassing quantification of interparticle forces, characterization of force network structures, and monitoring the dynamic evolution of force chains under different strain levels. The results show that the specimens with disturbed fabric led to a more dynamic and less persistent force network, more fabric instability, and thus more reorganization of force chain structures. In addition, the frequent rearrangement of the force network in the presence of water (with reduced inter-particle contact friction) likely exacerbates localized stress concentrations, promoting failure in previously unengaged or weakly connected particles. The results reported in this paper offer a new insight into how the initial fabric and wetting cause different fracture behavior. The findings can also pave the way for more in-depth future investigations into the mechanics governing particle fracture in granular assemblies.

Finite element analysis↗

DOE OSTI2021

Emerging contact force heterogeneity in ordered soft granular media

Under external perturbations, inter-particle forces in disordered granular media are well known to form a heterogeneous distribution with filamentary patterns. Better understanding these forces and the distribution is important for predicting the collective behavior of granular media, the media second only to water as the most manipulated material in global industry. However, studies in this regard so far have been largely confined to granular media exhibiting only geometric heterogeneity, leaving the dimension of mechanical heterogeneity a rather uncharted area. In this work, through a FEM contact mechanics model, we show that a heterogeneous inter-particle force distribution can also emerge from the dimension of mechanical heterogeneity alone. Specifically, we numerically study inter-particle forces in packing of mechanically heterogeneous disks arranged over either a square or a hexagonal lattice and under quasi-static isotropic compression. Our results show that, at the system scale, a hexagonal packing exhibit a more heterogeneous inter-particle force distribution than a square packing does; At the particle scale, for both packing lattices, preliminary analysis shows the consistent coexistence of outliers (i.e., softer disks sustaining larger forces while stiffer disks sustaining smaller forces) in comparison to their homogeneous counterparts, which implies the existence of nonlocal effect. Further analysis on the portion of outliers and on spatial contact force correlations suggest that the hexagonal packing shows more pronounced nonlocal effect over the square packing under small mechanical heterogeneity. However, such trend is reversed when assemblies becomes more mechanically heterogeneous. Lastly, we confirm that, in the absence of particle reorganization events, contact friction merely plays the role of packing stabilization while its variation has little effect on inter-particle forces and their distribution.

36 MATERIALS SCIENCE↗

DOE OSTI2025

pyTCR: A tropical cyclone rainfall model for python

pyTCR is a climatology software package developed in the Python programming language. It integrates the capabilities of several legacy physical models and increases computational efficiency to allow rapid estimation of tropical cyclone (TC) rainfall consistent with the large-scale environment. Specifically, pyTCR implements a horizontally distributed and vertically integrated model [Zhu et al., 2013] for simulating rainfall driven by TCs. Along storm tracks, rainfall is estimated by computing the cross-boundary-layer, upward water vapor transport caused by different mechanisms including frictional convergence, vortex stretching, large-scale baroclinic effect (i.e., wind shear), topographic forcing, and radiative cooling [Lu et al., 2018]. The package provides essential functionalities for modeling and interpreting spatio-temporal TC rainfall data. pyTCR requires a limited number of model input parameters, making it a convenient and useful tool for analyzing rainfall mechanisms driven by TCs. To sample rare (most intense) rainfall events that are often of great societal interest, pyTCR adapts and leverages outputs from a statistical-dynamical TC downscaling model [Lin et al., 2023] capable of rapidly generating a large number of synthetic TCs given a certain climate. As a result, pyTCR significantly reduces computational effort and improves the efficiency in capturing extreme TC rainfall events at the tail of the distributions from limited datasets. Furthermore, the TC downscaling model is forced entirely by large-scale environmental conditions from reanalysis data or coupled General Circulation Models (GCMs), simplifying the projection of TC-induced rainfall and wind speed under future climate using pyTCR. Finally, pyTCR can be coupled with hydrological and wind models to assess risks associated with independent and compound events (e.g., storm surges and freshwater flooding).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Effects of Urbanization in China on the East Asian Summer Monsoon as Revealed by Two Global Climate Models

Abstract The East Asian summer monsoon (EASM), vulnerable to anthropogenic activities, exerts profound impacts on precipitation in East Asia. In recent decades, China has been experiencing rapid urbanization. It has not yet reached a consensus on its influences on EASM. Here, we introduce satellite‐derived urban impervious areas of China during 1985–2017 into two global climate models with contrasting biases of climatological EASM simulation. Both models consistently show that urbanization weakens EASM by increasing surface friction drag. The increased surface roughness reduces climatological southerly winds in the concentrated urbanization area. As a result, moisture convergence occurs south of the concentrated urbanization area with ascending motion, while moisture divergence with descending motion appears to the north, resulting in the precipitation anomalies of “southern flood and northern drought.” The dominant role of urban dynamic forcing in regulating EASM is confirmed by reducing urban roughness to the level of croplands in sensitivity experiments. As urban dynamic forcing is largely mitigated, the urbanization impact mainly driven by urban thermal forcing conversely strengthens EASM due to enhanced land‐sea thermal contrast via urbanization‐induced warming. Spatial patterns and magnitudes of the urbanization‐induced EASM weakening vary between the two models due to disparities in their simulations of urban thermodynamic forcing and atmospheric adjustments. These results underscore the robust weakening effect of urbanization on EASM by urban friction drag and document uncertainties in the related spatial pattern and magnitude resulting from model discrepancies.

Meteorology & Atmospheric Sciences↗

DOE OSTI2020

Microscopic Theory of Long-Time Center-of-Mass Self-Diffusion and Anomalous Transport in Ring Polymer Liquids

We construct a microscopic theory at the level of segment-scale correlated space–time intermolecular forces for the long-time center-of-mass (CM) diffusion constant and intermediate-time non-Fickian transport in dense solutions and melts of ring polymers. The approach combines ideas of polymer, colloid, and liquid-state statistical mechanics to quantify how the multifractal intra-ring conformational structure and inter-ring packing correlations determine dynamic caging constraints and time-dependent friction. Breakdown of Rouse theory is predicted to occur due to length scale-dependent temporal correlation of forces exerted on pairs of tagged ring segments from surrounding polymers. At large enough degrees of polymerization (N), a stronger scaling of the CM diffusion constant (D ∝ N –2 ) is predicted, with a crossover N D proportional to the product of the system-specific macromolecular volume fraction and dimensionless compressibility. In analogy with caging effects in glass-forming fluids, the theory does appear to begin to fail at sufficiently high N/N D for the center-of-mass diffusivity, likely due to another crossover to an even slower activated transport regime. However, use of N/N D with the theoretically predicted N D , in conjunction with dynamic blob scaling ideas for local Rouse friction, collapses semidilute and concentrated solution simulation data onto a master curve. Based on the same physical ideas employed to predict the diffusion constant, a generalized Langevin equation description is formulated for intermediate-time CM transport. It predicts two subdiffusive regimes that emerge due to the self-similar nature of internal ring structure. Analytic and numerical predictions for the apparent non-Fickian exponent as a function of time, N/N D and dimensionless compressibility, the evolution of the maximum degree of subdiffusive motion with N/N D , and the time scale for recovering Fickian diffusion are made, all of which are in good accordance with melt simulations. Furthermore, the present work sets the stage to address activated dynamics and glass formation on the macromolecular scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Analytical modeling and simulation of electrical contact resistance for elastic rough electrode surface contact including frictional temperature rise

An improved electrical contact resistance (ECR) model for elastic rough electrode contact is proposed, incorporating the effects of asperity interactions and temperature rise by frictional and joule heating. The analytical simulation results show that the ECR decreases steeply at the beginning of the contact between Al and Cu. However, it becomes stabilized after reaching a specific contact force. It is also found that the longer elapsed sliding contact time, the higher ECR due to the increase in electrical resistivity of electrode materials by the frictional temperature rise at the interface. The effects of surface roughness parameters on ECR are studied through the 3 2 full-factorial design-of-experiment analysis. Based on the two representative roughness parameters, i.e., root-mean-square (rms) roughness and asperity radius, their individual and coupled effects on the saturated ECR are examined. The saturated ECR increases with the rms roughness for a rough machined surface condition, but it is hardly affected by the asperity radius. On the other hand, the saturated ECR increases with both the rms roughness and the asperity radius under a smooth thin film surface condition.

42 ENGINEERING↗