Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “combined data integration”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

System Modeling Frameworks for Wind Turbines and Plants: Review and Requirements Specifications

System modeling frameworks for wind turbines and plants are used by research groups and industry to design wind energy systems that take into account key trade-offs across performance, cost, and reliability at both the turbine and plant level. The frameworks are exercised using a variety of multi-disciplinary design, analysis and optimization (MDAO) methods. To improve inter-operability and foster collaboration, this report proposes a classification system for the frameworks along dimensions of model fidelity and scope. The classification system is first motivated with reviews the state-of-the-art in the development of software frameworks for integrated wind turbine and plant simulation. Within each major wind turbine and power plant subsystem, a matrix is developed for the disciplines used and the fidelity levels with which each discipline can be modeled. The existing frameworks are then classified according to the matrix. Next, an ontology is proposed that will allow for standardizing how data is transferred between the most common discipline-fidelity combinations used in the frameworks. A common representation of data creates the ability to 1) share system descriptions and analysis results, supporting more transparent benchmarks and comparison, and 2) integrate models together into workflows within and across organizations for improving the efficiency and performance of wind turbine and power plant design processes. Ultimately, this integration leads to better overall wind energy system designs with high performance and low costs.

17 WIND ENERGY↗
DOE OSTI2025

The Analysis Description Language Ecosystem: Latest developments and physics applications

We present latest developments in Analysis Description Language (ADL), a declarative domain-specific language describing the physics algorithm of a HEP data analysis decoupled from software frameworks. Analyses written in ADL can be integrated into any framework for various tasks. ADL is a multipurpose construct with uses ranging from analysis design to preservation, reinterpretation, queries, visualisation, combination, etc. The most advanced infrastructure to execute ADL on events is the CutLang runtime interpreter. Recent technical developments include an automated interface with different data types, generation of the abstract syntax tree, a visualization tool that that auto-converts analysis flows to graphs, incorporation of trained machine learning models and a Jupyter-based plotting tool. We also report physics implications including a large scale LHC analysis implementation and validation effort for beyond the standard model reinterpretation purposes and studies with ATLAS and CMS open data.

Sekmen, Sezen [Kyungpook National Univ., Daegu (Ko↗
DOE OSTI2021

Search for long-lived particles using displaced jets in proton-proton collisions at $\sqrt{s}$ = 13 TeV

An inclusive search is presented for long-lived particles using displaced jets. The search uses a data sample collected with the CMS detector at the CERN LHC in 2017 and 2018, from proton-proton collisions at a center-of-mass energy of 13 TeV. The results of this search are combined with those of a previous search using a data sample collected with the CMS detector in 2016, yielding a total integrated luminosity of 132 fb - 1 . The analysis searches for the distinctive topology of displaced tracks and displaced vertices associated with a dijet system. For a simplified model, where pair-produced long-lived neutral particles decay into quark-antiquark pairs, pair production cross sections larger than 0.07 fb are excluded at 95% confidence level (C.L.) for long-lived particle masses larger than 500 GeV and mean proper decay lengths between 2 and 250 mm. For a model where the standard model-like Higgs boson decays to two long-lived scalar particles that each decays to a quark-antiquark pair, branching fractions larger than 1% are excluded at 95% C.L. for mean proper decay lengths between 1 mm and 340 mm. A group of supersymmetric models with pair-produced long-lived gluinos or top squarks decaying into various final-state topologies containing displaced jets is also tested. Gluino masses up to 2500 GeV and top squark masses up to 1600 GeV are excluded at 95% C.L. for mean proper decay lengths between 3 and 300 mm. The highest lower bounds on mass reach 2600 GeV for long-lived gluinos and 1800 GeV for long-lived top squarks. These are the most stringent limits to date on these models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

ExaWorks: Workflows for Exascale

Exascale computers will offer transformative capabilities to combine data-driven and learning-based approaches with traditional simulation applications to accelerate scientific discovery and insight. These software combinations and integrations, however, are difficult to achieve due to challenges of coordination and deployment of heterogeneous software components on diverse and massive platforms. We present the ExaWorks project, which can address many of these challenges: ExaWorks is leading a co-design process to create a workflow Software Development Toolkit (SDK) consisting of a wide range of workflow management tools that can be composed and interoperate through common interfaces. We describe the initial set of tools and interfaces supported by the SDK, efforts to make them easier to apply to complex science challenges, and examples of their application to exemplar cases. Furthermore, we discuss how our project is working with the workflows community, large computing facilities as well as HPC platform vendors to sustainably address the requirements of workflows at the exascale.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Fracture Network Quantification during CO2 Injection

This is the presentation prepared for the ARMA 2025 (59th US Rock Mechanics/Geomechanics Symposium) Conference held in Santa Fe, New Mexico, June 8-11, 2025. Accurate mapping and quantification of these networks are essential to ensure the integrity of CO2 storage reservoirs, understand and reduce potential leakage, and maintain long-term environmental safety. This study presents a novel machine learning-driven approach, integrated with geomechanical analysis, to quantify fracture networks and assess their spatial distribution during CO2 injection. This paper combines microseismic monitoring data with principles of hydraulic diffusivity and geomechanical analysis to characterize reservoir scale fracture network. The novelty of our approach lies in its capacity to assimilate time-dependent pressure data and microseismicity into a cohesive framework, which not only identifies microseismic triggering fronts but also tracks fracture distribution during active injection. Besides, leveraging image log data and analysis our approach also provides another angle of the insights to solidate the fracture networks understanding and geomechanical impacts. Key results from our study include the detection of over 100 distinct fracture clusters across the injection site, with fracture orientations strongly correlated with the prevailing in-situ stress field.

CO2 storage and sequestration↗
DOE OSTI2025

Landscaper v1

Understanding the inner workings of machine learning models through their loss landscapes offers crucial insights into model properties, optimization dynamics, and generalizability. However, accessing these insights has traditionally required specialized mathematical expertise, limiting broader adoption. Landscaper is an open-source Python package designed to bridge this gap. Landscaper seamlessly integrates a suite of multi-dimensional loss landscape analyses with cutting-edge topological data analysis (TDA) methods. This powerful combination makes both fundamental loss landscape analysis and advanced TDA techniques accessible to the broader scientific ML community, without requiring deep pre-existing mathematical knowledge. Landscaper offers three key functionalities: * Construction: Builds detailed loss landscape representations through versatile low and high-dimensional sampling techniques. * Quantification: Applies advanced metrics, including a novel topological data analysis (TDA) based smoothness metric, enabling new perspectives on model behavior. * Visualization: Offers intuitive tools to visualize and interpret loss landscapes, providing actionable insights beyond traditional performance metrics.

Weber, Gunther [Lawrence Berkeley National Laborat↗
DOE OSTI2025

Fracture Network Quantification during CO2 Injection

This is the conference paper accompanying an oral presentation at the ARMA 2025 (59th US Rock Mechanics/Geomechanics Symposium) Conference held in Santa Fe, New Mexico, June 8-11, 2025. Accurate mapping and quantification of these networks are essential to ensure the integrity of CO2 storage reservoirs, understand and reduce potential leakage, and maintain long-term environmental safety. This study presents a novel machine learning-driven approach, integrated with geomechanical analysis, to quantify fracture networks and assess their spatial distribution during CO2 injection. This paper combines microseismic monitoring data with principles of hydraulic diffusivity and geomechanical analysis to characterize reservoir scale fracture network. The novelty of our approach lies in its capacity to assimilate time-dependent pressure data and microseismicity into a cohesive framework, which not only identifies microseismic triggering fronts but also tracks fracture distribution during active injection. Besides, leveraging image log data and analysis our approach also provides another angle of the insights to solidate the fracture networks understanding and geomechanical impacts. Key results from our study include the detection of over 100 distinct fracture clusters across the injection site, with fracture orientations strongly correlated with the prevailing in-situ stress field.

CO2 storage and sequestration↗
DOE OSTI2020

HPC for the EEC: Industrial Trade Tools for the Aging Energy Infrastructure

The petrochemical and refining sectors are challenged with reliably delivering safer and cleaner energy to US consumers, while meeting an ever-growing global demand. The petrochemical and refining infrastructure in the USA is aging, and corrosion and damage mechanisms are constant threats to mechanical integrity, safety, and profitability. However, governments and industry stakeholders are reluctant to replace or upgrade the existing infrastructure due to the immense cost. To better understand and mitigate the risks of aging energy infrastructure, strong technical analysis capabilities, combined with optimized monitoring and decision making, is critical. Today, this is accomplished via complex simulations and data analysis. To this end, advanced High-Performance-Computing (HPC) software will be leveraged and integrated in easy-to-use industrial trade tools, to lower the barrier for new users, increase the ease of access for experienced users, and allow for smarter decisions to be made in the midstream and downstream energy sectors.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Stepwise Dynamic Calibration of a Hydromechanical Simulation Using Time-Lapse Vertical Seismic Profile

This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse Vertical Seismic Profiling (VSP) integration. The selected study site is the 13-10A injector well within the ongoing CO2-EOR operation of the Farnsworth Field Unit. The Time-lapse VSP dataset carries the combined effects of fluid substitution and mean effective stress changes, thereby providing a dataset amenable for the calibration of production and injection-induced stress changes. The concept is similar to calibrating a reservoir simulation model in that the process honor real field data to set up an inverse problem. The solution optimizes the independent and impactful geomechanical parameters that replicate the observed time-lapse seismic velocity changes. This stress calibration is enabled by 4D geomechanical modeling and the VSP Integration workflow. This calibration benefits from extensive geological, geophysical and geomechanical characterization through 3D seismic data, geophysical well logs, and core assessed as part of the 1D MEM conducted on the 13-10A subject well. These data are used to develop a site-specific rock physics model. The Biot Gassmann workflow combines rock physics and reservoir simulation outputs to determine the fluid substitution contribution to seismic velocity change. Additionally, modeled seismic velocity attributed to mean effective stress are determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from the ultrasonic seismic velocity measurements on the extracted Morrow B core. A penalty function is then formed between the modeled seismic velocities and the observed time-lapse VSP dataset. Four independent and impactful geomechanical parameters have been determined. These are the bulk modulus and shear modulus for zero porosity and the shear and compressional seismic velocity to mean effective stress derivatives. The dataset of numerous coupled hydromechanical- geomechanical simulation realizations is built by combining variations of the four stated geomechanical parameters. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer are used to converge on the optimal geomechanical parameters. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling.

02 PETROLEUM↗
DOE OSTI2023

Combination of inclusive top-quark pair production cross-section measurements using ATLAS and CMS data at $ \sqrt{s} $ = 7 and 8 TeV

A combination of measurements of the inclusive top-quark pair production cross-section performed by ATLAS and CMS in proton–proton collisions at centre-of-mass energies of 7 and 8 TeV at the LHC is presented. The cross-sections are obtained using top-quark pair decays with an opposite-charge electron–muon pair in the final state and with data corresponding to an integrated luminosity of about 5 fb$^{−1}$ at $ \sqrt{s} $ = 7 TeV and about 20 fb$^{−1}$ at $ \sqrt{s} $ = 8 TeV for each experiment. The combined cross-sections are determined to be 178.5 ± 4.7 pb at $ \sqrt{s} $ = 7 TeV and $ {243.3}_{-5.9}^{+6.0} $ pb at $ \sqrt{s} $ = 8 TeV with a correlation of 0.41, using a reference top-quark mass value of 172.5 GeV. The ratio of the combined cross-sections is determined to be R$_{8/7}$ = 1.363 ± 0.032. The combined measured cross-sections and their ratio agree well with theory calculations using several parton distribution function (PDF) sets. The values of the top-quark pole mass (with the strong coupling fixed at 0.118) and the strong coupling (with the top-quark pole mass fixed at 172.5 GeV) are extracted from the combined results by fitting a next-to-next-to-leading-order plus next-to-next-to-leading-log QCD prediction to the measurements. Using a version of the NNPDF3.1 PDF set containing no top-quark measurements, the results obtained are $ {m}_t^{\textrm{pole}}={173.4}_{-2.0}^{+1.8} $ GeV and $ {\alpha}_{\textrm{s}}\left({m}_Z\right)={0.1170}_{-0.0018}^{+0.0021} $.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Combined searches for the production of supersymmetric top quark partners in proton–proton collisions at $\sqrt{s} = 13\,\text {Te}\text {V} $

A combination of searches for top squark pair production using proton–proton collision data at a center-of-mass energy of 13$\,\text {Te}\text {V}$ at the CERN LHC, corresponding to an integrated luminosity of 137$\,\text {fb}^{-1}$ collected by the CMS experiment, is presented. Signatures with at least 2 jets and large missing transverse momentum are categorized into events with 0, 1, or 2 leptons. New results for regions of parameter space where the kinematical properties of top squark pair production and top quark pair production are very similar are presented. Depending on the model, the combined result excludes a top squark mass up to 1325$\,\text {Ge}\text {V}$ for a massless neutralino, and a neutralino mass up to 700$\,\text {Ge}\text {V}$ for a top squark mass of 1150$\,\text {Ge}\text {V}$. Top squarks with masses from 145 to 295$\,\text {Ge}\text {V}$, for neutralino masses from 0 to 100$\,\text {Ge}\text {V}$, with a mass difference between the top squark and the neutralino in a window of 30$\,\text {Ge}\text {V}$ around the mass of the top quark, are excluded for the first time with CMS data. The results of theses searches are also interpreted in an alternative signal model of dark matter production via a spin-0 mediator in association with a top quark pair. Upper limits are set on the cross section for mediator particle masses of up to 420$\,\text {Ge}\text {V}$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

From Micro-environments to Macroscopic Effects: How the Alkaline Hydrogen Evolution Reaction Drives Cu Cathodic Corrosion

Cathodic corrosion of copper (Cu) has posed a significant challenge for over a century, impeding various technological progresses such as electrochemical conversion of CO 2 (eCO 2 RR) into fuels and other value-added carbon products. Here, in this study, employing a combined Density Functional Theory (DFT) and kinetic Monte Carlo (kMC) simulation approach, we delve into the atomistic level mechanism driving this phenomenon in Cu. Our hypothesis posits the pivotal role of alkaline hydrogen evolution reaction (HER) in facilitating cathodic corrosion in Cu. We rigorously develop a pH-dependent hydroxide (OH) adsorption mechanism and calculate the equilibrium OH coverage (𝜃 OH ) at varying pH levels, the thermodynamic stability of subsurface oxygen (O sub ), as well as the Cu-vacancy mediated diffusion of subsurface oxygens (O sub ). Through comprehensive analysis, we establish correlation among various microenvironments, including oxygen diffusion in subsurface layers, pH-dependent OH adsorption, and Cu dissolution into the electrolyte as (Cu-OH) complexes. Furthermore, our investigation explores the correlation between surface coordination environment of active sites and cathodic corrosion of Cu. Finally, by integrating DFT-derived thermodynamic data into a kMC model, we successfully predict the formation of experimentally observed corrosion pits on Cu-surfaces. This combined approach not only advances our fundamental understanding of Cu cathodic corrosion but also offers insights crucial for developing effective corrosion mitigation strategies.

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

Impedance-Based Detection of NO 2 Using Ni-MOF-74: Influence of Competitive Gas Adsorption

Chemically robust, low-power sensors are needed for the direct electrical detection of toxic gases. Metal–organic frameworks (MOFs) offer exceptional chemical and structural tunability to meet this challenge, though further understanding is needed regarding how coadsorbed gases influence or interfere with the electrical response. To probe the influence of competitive gases on trace NO 2 detection in a simulated flue gas stream, a combined structure–property study integrating synchrotron powder diffraction and pair distribution function analyses was undertaken, to elucidate how structural changes associated with gas binding inside Ni-MOF-74 pores correlate with the electrical response from Ni-MOF-74-based sensors. Data were evaluated for 16 gas combinations of N 2 , NO 2 , SO 2 , CO 2 , and H 2 O at 50 °C. Fourier difference maps from a rigid-body Rietveld analysis showed that additional electron density localized around the Ni-MOF-74 lattice correlated with large decreases in Ni-MOF-74 film resistance of up to a factor of 6 × 10 3 , observed only when NO 2 was present. These changes in resistance were significantly amplified by the presence of competing gases, except for CO 2 . Without NO 2 , H 2 O rapidly (<120 s) produced small (1–3×) decreases in resistance, though this effect could be differentiated from the slower adsorption of NO 2 by the evaluation of the MOF’s capacitance. Furthermore, samples exposed to H 2 O displayed a significant shift in lattice parameters toward a larger lattice and more diffuse charge density in the MOF pore. Evaluating the Ni-MOF-74 impedance in real time, NO 2 adsorption was associated with two electrically distinct processes, the faster of which was inhibited by competitive adsorption of CO 2 . Together, this work points to the unique interaction of NO 2 and other specific gases (e.g., H 2 O, SO 2 ) with the MOF’s surface, leading to orders of magnitude decrease in MOF resistance and enhanced NO2 detection. Finally, understanding and leveraging these coadsorbed gases will further improve the gas detection properties of MOF materials.

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

LevSeq: Rapid Generation of Sequence-Function Data for Directed Evolution and Machine Learning

Sequence-function data provides valuable information about the protein functional landscape but is rarely obtained during directed evolution campaigns. Here, we present Long-read every variant Sequencing (LevSeq), a pipeline that combines a dual barcoding strategy with nanopore sequencing to rapidly generate sequence-function data for entire protein-coding genes. LevSeq integrates into existing protein engineering workflows and comes with open-source software for data analysis and visualization. The pipeline facilitates data-driven protein engineering by consolidating sequence-function data to inform directed evolution and provide the requisite data for machine learning-guided protein engineering (MLPE). LevSeq enables quality control of mutagenesis libraries prior to screening, which reduces time and resource costs. Simulation studies demonstrate LevSeq’s ability to accurately detect variants under various experimental conditions. Lastly, we show LevSeq’s utility in engineering protoglobins for new-to-nature chemistry. Widespread adoption of LevSeq and sharing of the data will enhance our understanding of protein sequence-function landscapes and empower data-driven directed evolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Exascale workflow applications and middleware: An ExaWorks retrospective

Exascale computers offer transformative capabilities to combine data-driven and learning-based approaches with traditional simulation applications to accelerate scientific discovery and insight. However, these software combinations and integrations are difficult to achieve due to the challenges of coordinating and deploying heterogeneous software components on diverse and massive platforms. Here, we present the ExaWorks project, which addresses many of these challenges. We developed a workflow Software Development Toolkit (SDK), a curated collection of workflow technologies that can be composed and interoperated through a common interface, engineered following current best practices, and specifically designed to work on HPC platforms. ExaWorks also developed PSI/J, a job management abstraction API, to simplify the construction of portable software components and applications that can be used over various HPC schedulers. The PSI/J API is a minimal interface for submitting and monitoring jobs and their execution state across multiple and commonly used HPC schedulers. We also describe several leading and innovative workflow examples of ExaWorks tools used on DOE leadership platforms. Furthermore, we discuss how our project is working with the workflow community, large computing facilities, and HPC platform vendors to address the requirements of workflows sustainably at the exascale.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Offshore Wind Farm Turbine and Energy Storage Optimization

Abstract This paper evaluates the technical and economic feasibility of repurposing decommissioned offshore oil and gas platforms as electrical substations for offshore wind projects in the U.S. Gulf of America, a region characterized by relatively low and highly variable wind speeds, extensive legacy offshore infrastructure, and exposure to merchant electricity markets. A unified techno-economic framework is developed using the Repurposing Offshore Infrastructure for Continued Energy (ROICE) Economic Model (REM) to integrate Gulfspecific wind resource assessment, commercial wind turbine performance, offshore infrastructure cost modeling, and wholesale electricity market exposure. Gulf wind speed data are vertically extrapolated to turbine hub height and combined with manufacturer power curves to compute annual energy production and capacity factors across a broad portfolio of commercial turbines, enabling identification of turbine designs best suited for low-wind offshore environments. Hourly electricity price data from the Midcontinent Independent System Operator (MISO) day-ahead market are incorporated to characterize revenue potential, price volatility, and the temporal alignment between wind generation and market conditions. In addition, a conceptual framework for offshore battery energy storage system (BESS) integration is developed to support future investigation of market-responsive energy shifting at repurposed platforms. Results from the turbine evaluation demonstrate that machines with lower cut-in wind speeds and earlier ‘rated-power’ characteristics significantly outperform larger, industry-standard offshore turbines for the same net power under Gulf wind conditions, underscoring the need for region-specific technology selection. Market analysis further reveals substantial price variability and limited intrinsic alignment between wind production and high-price periods, motivating consideration of operational flexibility mechanisms. While storage optimization is not implemented in this study, the REM framework establishes a transparent and replicable foundation for co-evaluating turbine selection, infrastructure constraints, and market exposure, providing a practical pathway for assessing the potential role of repurposed offshore platforms in enabling economically viable offshore wind development in the Gulf of America.

02 PETROLEUM↗
DOE OSTI2025

Evidence for the Dimuon Decay of the Higgs Boson in 𝑝⁢𝑝 Collisions with the ATLAS Detector

A search for the dimuon decay of the Higgs boson is presented based on 𝑝⁢𝑝 collision data recorded by ATLAS during Run 3 of the Large Hadron Collider, corresponding to an integrated luminosity of 165 fb −1 at $\sqrt{s}$ = 13.6 TeV. To enhance the sensitivity, the results are combined with those from Run 2. An excess of events over the background is observed with a significance of 3.4⁢𝜎 (2.5⁢𝜎 expected). The best-fit signal strength is 𝜇 = 1.4 ± 0.4. This result provides evidence for the 𝐻 → 𝜇⁢𝜇 decay with ATLAS data and offers a direct probe of the Higgs-boson Yukawa coupling to second-generation fermions.

Hadron colliders↗
DOE OSTI2025

Energy-resolved neutron imaging and diffraction including grain orientation mapping using event camera technology

Time-of-flight neutron diffraction and energy-resolved imaging each provide unique perspectives into material properties. Neutron diffraction is useful for assessing microstructural parameters such as phase composition, texture, and dislocation densities, though it typically provides averaged data over the sampled volume. Energy-resolved imaging, on the other hand, offers both spatial and spectral information by detecting Bragg edges and neutron absorption resonances, which enables detailed mapping of microstructure and isotopic composition. When combined, these techniques have the potential to enrich our understanding of material behavior across different scales, enhancing our understanding of complex materials. Traditionally, these modalities are conducted on separate instruments, which is time-consuming and poses challenges for data integration. Here, we report the integration of the LumaCam, an event-mode energy-resolved neutron imaging camera with the HIPPO time-of-flight diffractometer at LANSCE. This integration enables simultaneous diffraction and imaging across the full spectrum, with analysis optimized for diffraction and Bragg-edge imaging in the thermal range (0.45–10 Å) and resonance imaging in the epithermal range (0.5–3000 eV), facilitating comprehensive multi-modal analysis. We demonstrate its capabilities through case studies, including spatial mapping of grain orientations in a steel sample and accurate thickness estimations for irregular samples including a depleted uranium cylinder and a natural silver-containing mineral specimen. The combined setup enhances real-time sample alignment and provides comprehensive data for crystal structure, texture, and isotopic composition analysis. This approach opens new possibilities for advanced applications in nuclear engineering, archaeology, and materials science.

36 MATERIALS SCIENCE↗