Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic 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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Ab-initio informed cluster dynamics simulation of self- and Xe diffusivity in uranium mononitride under irradiation

Uranium mononitride (UN) is one of the ceramic nuclear fuel alternatives to oxide fuels considered for light water reactor and advanced reactor designs, as it presents significant advantages such as high uranium density (better economics) and high thermal conductivity and melting point (increased safety). Self- and fission gas diffusivities need to be better understood, given that they influence key fuel performance phenomena like swelling and fission gas release. Recently, radiation enhanced diffusivity was investigated in UN by means of cluster dynamics simulations relying on empirical potential-based parameterizations, the reliability of which highly depends on the interatomic potential accuracy. Here, in this work, we refine this approach by determining, using ab-initio calculations, the properties of defect clusters containing vacancies, self-interstitials and Xe impurities. We also consider larger clusters than previous studies. The obtained dataset (formation enthalpies, entropies, and migration barriers) is used to parameterize a cluster dynamics model of mobile clusters, and to calculate the defect cluster concentrations under irradiation. This gives us access to the radiation enhanced self- and fission gas diffusivities. Although the resulting diffusivities are close to the values reported in the literature, we find important qualitative differences in the diffusion mechanisms. Capturing the correct mechanisms is crucial to properly describe the chemistry and fission rate dependence of the model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Coregistered positron emission particle tracking (PEPT) and X-ray computed tomography (CT) for engineering flow measurements

Increasingly, fully 3D experimental measurements of flow in complex engineering geometries are required to validate computational fluid dynamics models that support and inform reactor design and licensing. One barrier to such measurements is the complexity of typical reactor components and subsequent lack of optical access in these systems. To overcome this, the deployment of coregistered positron emission particle tracking (PEPT) and X-ray computed tomography (CT) is explored for flow measurement in reactor thermal hydraulic components and model (scaled) systems. Through this methodology, fully 3D flow information (via PEPT) and detailed internal geometry (via CT) are captured in opaque systems such as pipes, rod bundles, packed beds, etc. The reconstructed flow field and geometry can then be overlain to reveal detailed flow features around internal structures within a given test section. This is enabled through the use of a combined preclinical PET/CT scanner with overlapping PET and CT fields of view. Such measurements are useful for characterizing flow inside such intricate nuclear thermal hydraulic components as core geometries and heat exchangers, among others, and providing valuable 3D validation data for CFD models. In this work, basic tests of this 3D flow/geometry mapping are presented, and the implications of such measurements are discussed. Further, preliminary measurements are made with both point sources and flow in a simple pipe flow geometry to evaluate the capabilities of this technique. PEPT and CT features are coregistered with up to 0.1 mm precision, and pipe flow mean velocity and Reynolds stresses are reconstructed with similar accuracy to previous PEPT demonstrations. The utility of PEPT/CT is shown herein, and suggestions for future measurements are made.

3D flow measurement↗
DOE OSTI2026

Cosmic Overture Echoed in a Stellar Final Act: Implications of Nanohertz Gravitational Waves for Core Collapse Supernova Neutrinos

The growing evidence for nanohertz gravitational waves, from NANOGrav and other observations, may be pointing to a cosmological first-order phase transition at temperatures of 𝒪⁢(10−100)⁢MeV. Such an interpretation requires beyond the Standard Model dynamics in this energy range. If so, it may well be the case that galactic core-collapse supernova explosions would produce the key components related to the first-order phase transition, leaving detectable imprints on the spectrum of neutrinos emitted in the initial few seconds of the collapse. This might provide further evidence in support of the early universe interpretation of the nanohertz gravitational wave signal. The scenario proposed here is also suggestive of a low-mass seesaw mechanism to explain neutrino masses. We outline the prospects for future observations of Galactic supernovae to uncover the signals of this scenario, with further confirmation from future pulsar timing array measurements of primordial nanohertz gravitational waves.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

An integrated supply chain analysis for cobalt and rare earth elements under global electrification and constrained resources

As global electrification becomes more prominent, demand is projected to increase for electric vehicles (EVs), their components (EV batteries, rare earth traction motors), and their raw materials (cobalt, rare earths). To analyze the supply chain constraints associated with these components and raw materials, a system dynamics model was developed to simulate the interactions between the rare earth and cobalt supply chains. The results showed that under current understanding of raw material resources, production of the projected volume of EVs could be inadequate. To meet the aggressive EV volume sale projections set by the International Energy Agency, EV manufacturers will need to implement a large portion of traction motors without rare earths and utilize non-cobalt batteries as much as possible. On the supply side, a much faster ramp-up than what being modeled would be needed to meet much higher demand than current market forecasts.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Molecular Mass Growth Processes to Polycyclic Aromatic Hydrocarbons through Radical–Radical Reactions Exploiting Photoionization Reflectron Time-of-Flight Mass Spectrometry

Polycyclic aromatic hydrocarbons (PAHs) represent critical building blocks in molecular mass growth processes to carbonaceous nanoparticles, referred to as interstellar and circumstellar grains along with soot particles in astrophysical environments and combustion systems, respectively. Recent advancements on elucidating elementary steps to PAHs have utilized reactions of aromatic radicals, resonantly stabilized free radicals, and aliphatic radicals with closed shell hydrocarbons. However, the role of radical–radical reactions (RRRs) leading to PAHs has remained largely unexplored on the molecular level due to preceding experimental challenges in producing sufficiently high number densities of radical reactants for isomer-selective detection of products from bimolecular and termolecular reactions. This Account offers the latest developments in our knowledge on the mechanisms and pathways to PAHs via RRRs probed in a chemical microreactor at temperatures as high as 1600 K. Product preservation in a molecular beam coupled with synchrotron vacuum ultraviolet photoionization reflectron time-of-flight mass spectrometry and photoelectron photoion coincidence spectroscopy enabled isomer-selective detection of PAHs of up to three rings by their photoionization efficiency curves, which were fit with a linear combination of reference curves for identification. Experiments were combined with computational fluid dynamics modeling of the physicochemical processes in the microreactor, as well as high-level electronic structure calculations to reveal the reaction pathways of each system. Six distinct reaction mechanisms were discovered in this work: propargyl addition─benzannulation (PABA), methyl addition─ring expansion (MARE), cyclopentadienyl addition─naphthylization (CPAN), fulvenallenyl addition─cyclization─aromatization (FACA), benzyl addition─aromatization (BAA), and phenyl addition─pentacyclization (PAP). By systematically varying the number of carbon atoms in the radical reactants, molecular mass growth processes involving reactions between radicals with odd numbers of carbon atoms access aromatics carrying one, two, or three six-membered rings, whereas reactions between even- and odd-carbon-numbered radicals produce aromatics combining five- and six-membered rings. Our investigations reveal unconventional cycloadditions on excited state triplet surfaces, additions of radicals to low spin density carbon-centered radicals, spiroaromatic and fulvene-type intermediates, and highly strained bicyclic reaction intermediates, challenging current perceptions of PAH molecular mass growth processes. All of the listed mechanisms, except for FACA, feature endoergic reactions or barriers which lie above the separated reactants and therefore might be central to circumstellar environments of carbon-rich stars and planetary nebulae as their descendants, but they play no role in the gas phase of cold molecular clouds where temperatures as low as 10 K dominate. Altogether, this work provides detailed reaction mechanisms of PAH growth processes, advancing our knowledge of the chemistry of carbonaceous matter in the universe.

Addition reactions↗
DOE OSTI2025

Mechanical Roles of Polysaccharide Assembly and Interactions in Plant Cell Walls

Plants synthesize polysaccharide-based primary cell walls that possess unique microstructures and mechanical properties to accommodate plant growth and provide protection. Here, it remains challenging to assess the role of polysaccharide organization and interactions in the mechanical behavior of primary cell walls owing to their complex microstructure and highly nonlinear mechanical responses. Employing a coarse-grained molecular dynamics model developed for onion epidermal walls, this work explores the conditions under which polysaccharide assembly and interactions might play a significant role in primary cell wall mechanics. Cellulose–cellulose adhesion plays a dominant role in the wall load-bearing capacity, but when cellulose–cellulose adhesion was disrupted computationally, cellulose–xyloglucan adhesion could influence the wall load-bearing capacity. Contrary to the common concept that xyloglucans mechanically tether well-separated cellulose microfibrils, xyloglucans functioned in this case as interfibrillar adhesives capable of transmitting tensile forces between cellulose microfibrils. Our findings may inform design criteria of new materials inspired by plant cell walls.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Dynamic Process Intensification via Data-Driven Dynamic Optimization: Concept and Application to Ternary Distillation

Process intensification is a design philosophy aimed at making chemical processes safer and more efficient. Its implementation often results in significant modifications to the design and structure of the process, with several conventional unit operations occurring in the same physical device. Traditionally, process intensification has focused on steady-state operation. In our previous works, we introduced dynamic process intensification (DPI) as a new intensification paradigm based on operational changes for conventional or intensified units. DPI is predicated on switching operation between two auxiliary steady states selected via a steady-state optimization calculation that ensures that the system generates, on average and over time, the same products as in nominal steady-state operation, but with favorable economics. This paper extends the DPI concept and introduces a novel dynamic optimization-based DPI strategy (Do-DPI) that involves imposing a true cyclic operation rather than switching between two discrete states. We discuss its implementation using surrogate dynamic models learned via system identification. Here, an extensive case study concerning a ternary distillation column separating a canonical hydrocarbon mixture shows that Do-DPI can reduce energy use by more than 4% relative to steady-state operation, with no significant deviations in product quality and production rate.

42 ENGINEERING↗
DOE OSTI2021

Signatures of Coherent Phonon Transport in Ultralow Thermal Conductivity Two-Dimensional Ruddlesden–Popper Phase Perovskites

An emerging class of methylammonium lead iodide (MAPbI 3 )-based Ruddlesden-Popper (RP) phase perovskites, BA 2 MA n-1 Pb n I 3n+1 (n = 1-7), exhibit enhanced stability to environmental conditions relative to MAPbI 3 , yet still degrade at elevated temperatures. Here, we experimentally determine the thermal conductivities of these layered RP phases for n = 1-6, where n defines the number of repeated perovskite octahedra per layer. We measure thermal conductivities of 0.37 ± 0.13/0.12, 0.17 ± 0.08/0.07, 0.21 ± 0.05/0.04, and 0.19 ± 0.04/0.03 W/m∙K in thin films of n = 1-4 and 0.08 ± 0.06/0.04, 0.06 ± 0.04/0.03, 0.06 ± 0.03/0.03, and 0.08 ± 0.07/0.04 W/m∙K in single crystals of n = 3-6. With the exception of n = 1, these thermal conductivities are lower than the range of 0.34-0.50 W/m∙K reported for single-crystal MAPbI 3 . Reduced-order lattice dynamics modeling suggests that the initially decreasing trend of thermal conductivity in similarly oriented perovskites with increasing n may result from the transport properties of coherent phonons, emergent from the superstructure, that do not scatter at the interfaces of organic butylammonium chains and perovskite octahedra. Reduced group velocity of coherent phonons in n = 3-6, a consequence of band flattening in the phonon dispersion, is primarily responsible for their ultralow thermal conductivities. Similar effects on thermal conductivity have been experimentally demonstrated in deposited superlattices, but never in naturally defined materials such as RP phases. GIWAXS measurements reveal that higher n RP phase thin films are less orientationally controlled and therefore possess apparently elevated thermal conductivities relative to single crystals of the same n.

2D perovskites↗
DOE OSTI2025

Storage-Induced Collapse of Lignin Macromolecular Structure and Its Impacts on the Biorefinery

Lignin plays a vital role in the economics of biorefineries, serving as a source of process energy and a feedstock for sustainable fuels and chemical production. While understanding lignin’s chemical composition is crucial, emerging evidence suggests that a more comprehensive understanding of its macromolecular structure is critical to explaining its complex behavior in the biorefinery. This study investigated the partial collapse of the lignin network in corn stover feedstock after harvest and storage as a result of the microbial digestion of hemicellulose. Fluorescence microscopy was used to detect the collapse of lignin in terms of lignin’s inter-molecular interaction and the re-orientation of lignin’s chromophores, by the changes in lignin’s fluorescence lifetime, anisotropy, and the number of effective emitters. With minimal sample perturbation, our in-situ microscopic results revealed lignin's coil-globule transition phenomena, which was only previously predicted by molecular dynamics modeling extracted lignin in solvent. This collapse of lignin macromolecular structure was confirmed by results from NMR, IR, Raman, and powder X-ray diffraction. We also investigated the impact of this storage-induced collapse on the downstream biorefinery processes. Our study revealed that the two major approaches for lignin valorization in the lignin-first biorefinery model, namely monomer extraction and milled wood lignin extraction, were negatively impacted by the lignin collapse. As changes during storage are a source of feedstock variability, our study highlights the importance of understanding the effect of feedstock handling on biorefinery operations and economics.

09 BIOMASS FUELS↗
DOE OSTI2020

Permeability Hysteresis From Microchannels Opening During Dissolution/Reprecipitation Cycle

Abstract Permeability is a critical parameter for geological resources characterization. Its evolution with respect to porosity is particularly interesting and many research initiatives focus on deriving such relationships, to understand some hydraulic impacts of microstructure alteration. Permeability evolution from chemical reactions, for instance, can become complex as flow channels may open during rock dissolution. In this contribution, we show that this phenomenon can lead to irregular porosity‐permeability curves and permeability hysteresis after reprecipitation. Current approaches describing permeability as a simple function of porosity can therefore not capture this behavior, and we advocate instead the use of dynamic modeling for such scenarios. We demonstrate our approach by modeling a dissolution/precipitation cycle for a unit cell pore channel and quantify the process at larger scale on three different rock samples, whose microstructures are reconstructed from segmented micro‐Computerised Tomography scans.

Lesueur, Martin↗
DOE OSTI2021

Sensitivity of Grass Fires Burning in Marginal Conditions to Atmospheric Turbulence

Abstract Atmospheric forcing and interactions between the fire and atmosphere are primary drivers of wildland fire behavior. The atmosphere is known to be a chaotic system that, although deterministic, is very sensitive to small perturbations to initial conditions. We assume that as a result of the tight coupling between fire and atmosphere; wildland fire behavior, in turn, should also be sensitive to perturbations in atmospheric initial conditions. Observations suggest that low intensity prescribed fire, in particular, is susceptible to small perturbations in the wind field, which can significantly alter fire spread. Here, we employ a computational fluid dynamics model of coupled fire‐atmosphere interactions to answer the question: How sensitive is fire behavior to small variations in atmospheric turbulence? We perform ensemble simulations of fires in homogenous grass fuels. The only difference between ensemble members is the state of the turbulent atmosphere provided to the model throughout the simulation. The atmospheric state is a function of the initial conditions applied at the start of the simulation and boundary conditions applied throughout the simulation. We find a wide range of outcomes, with area burned ranging from 2,212 to 11,236 m 2 (>400% change), driven primarily by sensitivity to initial conditions, with nonnegligible contributions from boundary condition variability during the initial 30 s of simulation. Our results highlight the need for ensemble simulations, especially when considering fire behavior in marginal burning conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Hyporheic Exchange in Sand Dunes Under a Freely Deforming River Water Surface

Bedform-driven hyporheic exchange is conditioned by the head gradients at the sediment-water interface. Local exchange phenomena between the surface and the subsurface are often driven by the dynamic forces related to the velocity distributions around sediment waves. In open channels, the static forces are represented by the water depth, yet most computational fluid dynamics models use a rigid-lid approximation. Here, we investigated whether and when the deformation of the river’s free-surface influences bedform-driven hyporheic exchange. This was done through simulations of coupled open channel and hyporheic flows with the air-water interface modeled either as a free-surface or a rigid-lid across increasing subcritical Froude numbers from ~0.05 to ~0.95. The normalized hyporheic flux was higher when considering free-surface deformation across most of the range of Froude numbers considered. When the Froude number was larger than 0.6 and smaller than 0.85, both hydrostatic and nonhydrostatic-driven fluxes increased significantly compared to the rigid-lid approach and the total flux was about five times that predicted with the rigid-lid. These results indicate that bedform-driven hyporheic fluxes are underestimated in most studies that use a rigid-lid assumption in subcritical flows.

58 GEOSCIENCES↗
DOE OSTI2022

Characterization of Acoustic Emissions From Analogue Rocks Using Sparse Regression‐DMDc

Abstract Moisture loss in rock is known to generate acoustic emissions (AE). Phenomena that result in AE during drying are related to the movement of fluids through the pores and induced‐cracks that arise from differential mineral shrinkage, especially in clay‐bearing rock. AE from the movement of fluids occurs from the reconfiguration of fluid interfaces during drying, while AE from mineral shrinkage involves the debonding within or between minerals. Here, analogue rock samples were used to examine the differences in the AE signatures when one or both AE source‐types are present. An unsupervised sparse regression model, Dynamic Mode Decomposition with control, that extends Dynamic Mode Decomposition is used to characterize the AE signals recorded during the drying of porous analogue rock samples fabricated with ordinary Portland cement, with and without clay. This method can effectively and accurately reconstruct acoustic signals emitted from samples that only experience moisture loss without cracking. However, the method struggles to reconstruct signals from samples with intricate crack networks that formed during drying because AE generating mechanisms can emit contemporaneously, and the resulting waves propagate through drying‐induced cracks that can lead to multiple internal reflections. Thus, the differential reconstruction accuracy of time series generated by different underlying physical processes provides a robust filter for reducing large data catalogs. In general, both dynamics and sparse initiating events are learned directly from data and this method exposes a data hierarchy based on the complexity of the intrinsic dynamics.

58 GEOSCIENCES↗
DOE OSTI2023

Use of Shallow Ice Core Measurements to Evaluate and Constrain 1980–1990 Global Reanalyses of Ice Sheet Precipitation Rates

Abstract Sea‐level rise (SLR) projections by Earth System Models (ESMs) depend on ice sheet surface mass balances. Accurate, global atmosphere reanalyses would be ideal for providing equilibrated ice sheet model initial conditions in fully coupled ESM simulations. Here we present the first evaluation of 1980–1990 global reanalysis precipitation over Greenland and Antarctica that uses independent observations of net accumulation rates derived from shallow ice cores. Precipitation distributions from both the European Centre for Medium‐Range Weather Forecast's Reanalysis (ERA5) and the Modern‐Era Retrospective Analysis for Research and Applications (MERRA‐2) are highly correlated with contemporaneous co‐located net accumulation rates from Greenland ( r 2 > 0.95) and West Antarctica ( r 2 > 0.7). Three other commonly used reanalyses (WFDE5, CRUNCEP, and GSWP3) exhibit significantly weaker correlations on one or both ice sheets. Our findings imply that ESMs should use ERA5 or MERRA‐2 in data‐forced simulations to validate ice sheet model dynamics and precondition firn for SLR projections.

Schneider, Adam↗
DOE OSTI2022

Sugar loading is not required for phloem sap flow in maize plants

We report phloem transport of photoassimilates from leaves to non-photosynthetic organs, such as the root and shoot apices and reproductive organs, is crucial to plant growth and yield. For nearly 90 years, evidence has been generally consistent with the theory of a pressure-flow mechanism of phloem transport. Central to this hypothesis is the loading of osmolytes, principally sugars, into the phloem to generate the osmotic pressure that propels bulk flow. Here we used genetic and light manipulations to test whether sugar import into the phloem is required as the driving force for phloem sap flow. Using carbon-11 radiotracer, we show that a maize sucrose transporter1 (sut1) loss-of-function mutant has severely reduced export of carbon from photosynthetic leaves (only ~4% of the wild type level). Yet, the mutant remarkably maintains phloem pressure at ~100% and sap flow speeds at ~50–75% of those of wild type. Potassium (K+) abundance in the phloem was elevated in sut1 mutant leaves. Fluid dynamic modelling supports the conclusion that increased K+ loading compensated for decreased sucrose loading to maintain phloem pressure, and thereby maintained phloem transport via the pressure-flow mechanism. Furthermore, these results suggest that sap flow and transport of other phloem-mobile nutrients and signalling molecules could be regulated independently of sugar loading into the phloem, potentially influencing carbon–nutrient homoeostasis and the distribution of signalling molecules in plants encountering different environmental conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Genetic circuit design automation for yeast

Cells can be programmed to monitor and react to their environment using genetic circuits. Design automation software maps a desired circuit function to a DNA sequence, a process that requires units of gene regulation (gates) that are simple to connect and behave predictably. This poses a challenge for eukaryotes due to their complex mechanisms of transcription and translation. To this end, we have developed gates for yeast (Saccharomyces cerevisiae) that are connected using RNA polymerase flux as the signal carrier and are insulated from each other and host regulation. They are based on minimal constitutive promoters (~120 base pairs), for which rules are developed to insert operators for DNA-binding proteins. Using this approach, in this work we constructed nine NOT/NOR gates with nearly identical response functions and 400-fold dynamic range. In circuits, they are transcriptionally insulated from each other by placing ribozymes downstream of terminators to block nuclear export of messenger RNAs resulting from RNA polymerase readthrough. Based on these gates, Cello 2.0 was used to build circuits with up to 11 regulatory proteins. A simple dynamic model predicts the circuit response over days. Genetic circuit design automation for eukaryotes simplifies the construction of regulatory networks as part of cellular engineering projects, whether it be to stage processes during bioproduction, serve as environmental sentinels or guide living therapeutics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Traceable and Scalable Food Balance Sheets from Agricultural Commodity Supply and Utilization Accounts (2010–2022)

Abstract The Food Balance Sheets (FBS), compiled by the Food and Agriculture Organization (FAO), serve as a cornerstone dataset for studies on agricultural development, food security, and dietary health, providing a broad overview of global and regional food systems. However, its limited transparency and scalability hinder its application in empirical analysis and multisector dynamic modeling. Here, we present a traceable Food Balance Sheets (T-FBS) dataset, developed from detailed Supply Utilization Accounts (SUA) using a novel Primary Commodity equivalent (PCe) aggregation approach. This framework enables the aggregation of commodity flows along supply chains while ensuring consistency and balance across multiple dimensions. The T-FBS dataset includes 57 PCe commodities across 195 regions for the period 2010–2022, consolidated from over 500 SUA products. While T-FBS closely aligns with FAO-FBS at aggregate levels for dietary energy and macronutrients, it identifies key uncertainties in other elements (e.g., feed, trade, stocks). By enhancing methodological transparency, traceability, and scalability, T-FBS strengthens the robustness of food system studies and fosters future research and collaboration within the open-source community.

agriculture↗
DOE OSTI2020

Decadal predictability of North Atlantic blocking and the NAO

Can multi-annual variations in the frequency of North Atlantic atmospheric blocking and mid-latitude circulation regimes be skilfully predicted? Recent advances in seasonal forecasting have shown that mid-latitude climate variability does exhibit significant predictability. However, atmospheric predictability has generally been found to be quite limited on multi-annual timescales. New decadal prediction experiments from NCAR are found to exhibit remarkable skill in reproducing the observed multi-annual variations of wintertime blocking frequency over the North Atlantic and of the North Atlantic Oscillation (NAO) itself. This is partly due to the large ensemble size that allows the predictable component of the atmospheric variability to emerge from the background chaotic component. The predictable atmospheric anomalies represent a forced response to oceanic low-frequency variability that strongly resembles the Atlantic Multi-decadal Variability (AMV), correctly reproduced in the decadal hindcasts thanks to realistic ocean initialization and ocean dynamics. The occurrence of blocking in certain areas of the Euro-Atlantic domain determines the concurrent circulation regime and the phase of known teleconnections, such as the NAO, consequently affecting the stormtrack and the frequency and intensity of extreme weather events. Therefore, skilfully predicting the decadal fluctuations of blocking frequency and the NAO may be used in statistical predictions of near-term climate anomalies, and it provides a strong indication that impactful climate anomalies may also be predictable with improved dynamical models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗