Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Porous flow”

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 505 records · Page 28

Physics-Informed Neural Network Method for Forward and Backward Advection-Dispersion Equations

Advection-dispersion equations (ADEs) are commonly used to describe transport phenomena in porous media. Even though mature discretization-based numerical methods for ADEs exist, some challenges still remain, especially when it comes to solving advection-dominated forward ADEs and diffusion-dominated backward ADEs. The latter problem usually arises in the source identification context and leads to numerically unstable grid-based solutions that require a form of regularization or should be treated as an inverse problem that is computationally more expensive because it requires solving the forward problem multiple times. In this study, we propose a discretization-free approach based on the physics-informed neural network (PINN) method for solving coupled ADE and Darcy flow equations with space-dependent hydraulic conductivity. In this approach, the hydraulic conductivity, hydraulic head, and concentration fields are approximated with deep neural networks (DNNs). We assume that the conductivity field is given by its values on a grid, and we use these values to train the conductivity DNN. The head and concentration DNNs are trained by minimizing the residuals of the flow equation and ADE and using the initial and boundary conditions as additional constraints. The PINN method is applied to one- and two-dimensional forward ADE problems, where its performance for various P\'{e}clet numbers ($Pe$) is compared with the analytical and numerical solutions. We find that the PINN method is accurate with errors of less than 1\% and outperforms some conventional discretization-based methods for $Pe$ larger that 100. Next, we demonstrate that the PINN method remains accurate for the backward ADEs, with the relative errors in most cases staying under 5\% compared to the reference concentration field. Finally, we show that when available, the concentration measurements can be easily incorporated in the PINN method and significantly improve (by more than 50\% in the considered cases) the accuracy of the PINN solution of the backward ADE.

He, Qizhi↗

Real-time streaming tomographic reconstruction with on-demand data capturing and 3D zooming to regions of interest

Complex dynamic tomographic experiments at brilliant X-ray light sources require real-time feedback on the sample changes with respect to environmental conditions, selecting representative regions of interest for high-resolution scanning, and on-demand data saving mechanisms for storing only relevant projections acquired by fast area detectors and reducing data volumes. Here the implementation details of a 3D real-time imaging monitoring instrument, with zooming to a volume of interest with easy-to-use visualization via ImageJ , a tool familiar to most beamline users, is presented. The instrument relies on optimized data flow between the detector and processing machines and is implemented on commodity computers. The instrument has been developed at beamline 2-BM of the Advanced Photon Source, where the automatic lens changing mechanism for zooming is implemented with an Optique Peter microscope. Performance tests demonstrate the ability to process more than 3 GB of projection data per second and generate real-time 3D zooming with different magnification. These new capabilities are essential for new APS Upgrade instruments such as the projection microscope under development at beamline 32-ID. The efficacy of the proposed instrument was demonstrated during an in situ tomographic experiment on ice and gas hydrate formation in porous samples.

3D zooming↗

Synthesis of results for Brine Availability Test in Salt (BATS) DECOVALEX-2023 Task E

In the 2023 phase of the international collaborative DECOVALEX modeling project, Task E focused on understanding thermal, hydrological, and mechanical (THM) processes related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and non-porous in the far field. Investigation of the short-term, near-field behavior is important for radioactive waste disposal because this early period strongly controls the amount of inflowing brine. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment.The Task was divided into steps. Step 0 included matching unheated brine inflow data from boreholes at the Waste Isolation Pilot Plant (WIPP) and matching temperature observations during a Brine Availability Test in Salt (BATS) heater test. Step 1 included validation of models against a thermo-poroelastic analytical solution, and two-phase flow around an excavation. Finally, Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the same BATS heater test.There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. Task E was a learning experience for the teams involved, and feedback from the modeling teams has led to changes in follow-on BATS experiments at WIPP. The primary Task E lessons learned were the impact of hydrologic initialization methods (wetting up vs. drying down), the difference between confined and unconfined thermal expansion, and the large changes in permeability associated with heating and cooling.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

MODEL DEVELOPMENT AND VERIFICATION OF A 1,350°C AIR RECEIVER TEST SYSTEM

Open-loop volumetric receivers work with air at atmospheric pressure and are suitable for single-cycle or multicycle power plants. A novel additively manufactured (AM) silicon carbide (SiC), ceramic matrix composite (CMC) volumetric receiver was developed that consisted of a lattice structure, which absorbs solar radiation and converts it into heat energy. Heat energy from the porous receiver was transferred to heated ambient air. The receiver acted as a convective heat exchanger, transferring heat to the fluid through convection. Engineering was facilitated to develop a high-temperature air receiver test bed at Sandia national Laboratories (SNL) capable of demonstrating a 50kWth open-loop volumetric SiC air receiver module developed by General Electric Aerospace Research (GE Aerospace). This paper presents the development of the high temperature AM SiC-CMC air receiver, the test bed, and the first testing operations of this device, which was experimentally demonstrated to achieve 1,350°C for over 3.5 hours of operation, based on input flux levels of 67–91 W/cm2 0.5-2.5 kg/s flow rates. From the test campaign results, based on the measured convective heat transfer to the airflow and bounded radiative-loss assumptions, the solar-thermal efficiency of HOTSSTAR test module was estimated to be approximately as high as 74% for representative test conditions.

14 SOLAR ENERGY↗

Surface ionization waves propagating over non-planar substrates: wavy surfaces, cut-pores and droplets

Atmospheric pressure plasmas intersecting with dielectric surfaces will often transition into surface ionization waves (SIWs). Several applications of these discharges are purposely configured to be SIWs. During propagation of an SIW over a dielectric surface, the plasma charges the surface while responding to changes in geometrical and electrical material properties. This is particularly important for non-planar surfaces where polarization of the dielectric results in local electric field enhancement. In this paper, we discuss results from computational investigations of negative and positive SIWs propagating over nonplanar dielectrics in three configurations—wavy surfaces, cuts through porous materials and water droplets on flat surfaces. Here we found that negative SIWs are particularly sensitive to the electric field enhancement that occurs at the crests of non-planar surfaces. The local increase in ionization rates by the electric field enhancement can result in the SIW detaching from the surface, which produces non-uniform plasma exposure of the surface. Positive SIWs tend to adhere to the surface to a greater degree. These trends indicate that treatment of pathogen containing droplets on surfaces may be best performed by positive SIWs. The same principles apply to the surfaces cut through pores. Buried pores with small openings to the SIW may be filled by plasma by either flow of plasma into the pore (large opening) or initiated by photoionization (small opening), depending on the size of the opening compared to the Debye length.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A hybrid machine-learning approach for analysis of methane hydrate formation dynamics in porous media with synchrotron CT imaging

Fast multi-phase processes in methane hydrate bearing samples pose a challenge for quantitative micro-computed tomography study and experiment steering due to complex tomographic data analysis involving time-consuming segmentation procedures. This is because of the sample's multi-scale structure, which changes over time, low contrast between solid and fluid materials, and the large amount of data acquired during dynamic processes. Here, a hybrid approach is proposed for the automatic segmentation of tomographic data from time-resolved imaging of methane gas-hydrate formation in sandy granular media, which includes a deep-learning 3D U-Net model. To prepare a training dataset for the 3D U-Net, a technique to automate data labeling based on sample-specific information about the mineral matrix immobility and occasional fluid movement in pores is proposed. Automatic segmentation allowed for studying properties of the hydrate growth in pores, as well as dynamic processes such as incremental flow and redistribution of pore brine. Results of the quantitative analysis showed that for typical gas-hydrate stability parameters (100 bar methane pressure, 7°C temperature) the rate of formation is slow (less than 1% per hour), after which the surface area of contact between brine and gas increases, resulting in faster formation (2.5% per hour). Hydrate growth reaches the saturation point after 11 h of the experiment. Finally, the efficacy of the proposed segmentation scheme in on-the-fly automatic data analysis and experiment steering with zooming to regions of interest is demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Alkaline Zinc Passivation Mechanism is Controlled by Hydroxide Concentration

We studied Zn passivation and oxide growth in Zincate (Zn(OH) 4 2− ) in 4 and 8 M KOH solutions using an electrochemical quartz crystal microbalance (EQCM), building on our initial work at 1 M KOH where passivation was kinetically controlled. A porous passivating oxide spontaneously forms on Zn electrodes when KOH is above 4 M and saturated with zincate. However, passivation does not occur when bulk zincate concentration is decreased, resulting in continual Zn dissolution. EQCM data suggests that the passivation mechanism is strongly affected by pOH. Mass transport and kinetic processes in the 4 M KOH electrolytes couple and govern Zn passivation. At 8 M, KOH concentration shifts passivation to mass transport control. We explain this by the increased solubility of Zn(OH) 3 − with increasing pOH. The variation in the mechanism of passivation implications for how passivation is handled in Zn-alkaline batteries. The importance of controlling the mass transport increases with increased pOH, suggesting that electrode design, additives, and potential flow should be optimized. At lower pOHs, kinetics and mass transport must be balanced to manage passivation effectively. Additionally, the changing nature of the native oxide layer on the surface has implications for the evenness of deposition and dissolution on the Zn electrode.

Wittman, Reed M. [Oak Ridge National Laboratory (O↗

Combining Distributed Electrochemical Measurements with a Semi-Empirical Model to Identify Local Variations in Overpotentials for PEM Water Electrolyzers

Abstract A segmented cell and PCB approach was used to demonstrate an advanced technique for measuring and understanding the electrochemical behavior of PEM water electrolyzers. By employing a semi-empirical model, discrete modes of overpotential and their spatial distributions were identified in an operating cell. Distributed area-specific resistance (DASR measurements indicate that local dehydration of the polymer electrolyte membrane is the primary driver of increased polarization, behaving similarly to mass transport-limited currents at low flow rates due to ineffective gas-phase bubble removal. To improve quantification, the overall ohmic overpotential is divided into background (η Ohm,b ) and mass transport-coupled (η Ohm,mt ) components. Two types of anode diffusion media, a baseline PTL and a thin foil liquid/gas diffusion layer (LGDL), were examined. The results show significant variations in η Ohm,b and highlight that ineffective bubble removal leads to pronounced local dehydration. The effect of cell temperature on various overpotentials was also examined. Negligible differences in diffusion overpotential were observed when comparing PEWE performance using the PTL and LGDL, despite major differences in two-phase fluid transport within each diffusion media. This underscores the critical role of the catalyst layer's micro-porous structure in diffusion processes. The study highlights the complex nature and intricacies of mass transport limitations within PEWEs.

Roy, Anirban (ORCID:0000000163803386)↗

Mineral Scale Prevention on Electrically Conducting Membrane Distillation Membranes Using Induced Electrophoretic Mixing

The growth of mineral crystals on surfaces is a challenge across multiple industrial processes. Membrane-based desalination processes, in particular, are plagued by crystal growth (known as scaling), which restricts the flow of water through the membrane, can cause membrane wetting in membrane distillation, and can lead to the physical destruction of the membrane material. Scaling occurs when supersaturated conditions develop along the membrane surface due to the passage of water through the membrane, a process known as concentration polarization. To reduce scaling, concentration polarization is minimized by encouraging turbulent conditions and by reducing the amount of water recovered from the saline feed. In addition, antiscaling chemicals can be used to reduce the availability of cations. Here, we report on an energy-efficient electrophoretic mixing method capable of nearly eliminating CaSO4 and silicate scaling on electrically conducting membrane distillation (ECMD) membranes. The ECMD membrane material is composed of a percolating layer of carbon nanotubes deposited on porous polypropylene support and cross-linked by poly(vinyl alcohol). The application of low alternating potentials (2 Vpp,1Hz) had a dramatic impact on scale formation, with the impact highly dependent on the frequency of the applied signal, and in the case of silicate, on the pH of the solution.

54 ENVIRONMENTAL SCIENCES↗

Preactivated zeolite nanosheet plate-tiled membrane on porous PVDF film: Synthesis and study of proton-selective ion conduction

Preactivated MFI zeolite nanosheet plates (ZNPs) with large areas (~2.0 × 2.0 μm 2 ) and nanometer thicknesses (~60 nm) were prepared directly in liquid-dispersed state. The individual ZNP was a stack of 4-nm-thick single crystalline zeolite nanosheets (ZNs) interlinked by Si–O–Si bonds between neighboring ZN surfaces. The dispersed open-pore ZNPs allowed formulating suspensions for tiling ZNP membranes on polymer substrates without requiring post-coating activation. A pinhole-free ZNP-tiled (ZNPT) membrane has been achieved on macroporous PVDF film by a self-repairing vacuum-assisted filtration coating method. The resultant ZNPT layer was ~500 nm-thick consisting of about 7 ZNP layers with inter-ZNP width and inter-ZNP entrance porosity around 9 nm and 2%, respectively. The ZNPT-PVDF membrane exhibited high selectivity to proton transport over vanadyl ion and low resistances to proton conduction in aqueous solutions. The membrane also demonstrated to function as an efficient ion separator for the vanadium redox flow battery. The reported synthesis of preactivated ZNP suspension may overcome the major hurdle to realizing polymer-supported ZN membranes, which is the inefficiency of existing methods in obtaining readily dispersible open-pore ZNs. Finally, the ZNPT-PVDF membrane can be a more affordable and sustainable alternative to the Nafion-based ion separation membranes.

36 MATERIALS SCIENCE↗

Model-based optimization strategies for direct hydrogenation of carbon dioxide to dimethyl ether

Here, this study discusses model-based optimization strategies for CO 2 hydrogenation to dimethyl ether (DME) over aCuZnZr (CZZ) and ferrite (FER) mixed catalyst system in a packed-bed reactor configuration. A two-dimensional axisymmetric, nonisothermal packed-bed reactor model was developed using COMSOL Multiphysics 6.2 software. The model solves two-dimensional (radial and axial) heat and mass transport equations in the packed-bed and integrates intraparticle diffusion and heat transfer in a 1D approach. This powerful feature differs from a traditional porous media approach and takes into account any heat and mass transfer limitations that may exist. Analysis shows that the heat transfer limitations are negligible, but strong internal mass transfer limitations were observed at 10 ≤ WHSV ≤ 90 h -1 on the FER catalyst and at 240 °C. The optimum catalyst composition (i.e., mixing ratio) varies depending on the operating regime. The FER catalyst weight in the mixture can be as low as 5 wt.%, but the ideal composition de-pends on the internal mass transfer limitation and its relationship with the operating regime (i.e., weight hourly space velocity, temperature). A catalyst composition of 80 wt.% CZZ and 20 wt.% FER was suggested; this composition can provide high CO 2 conversion and DME production rates at a wide range of temperatures and flow rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Preparation of illite coated geomaterial microfluidic surfaces: Effect of salinity and heat treatment

Unique characteristics of clay (e.g., high surface area, deficiency in positive charges) make any clay abundant shales a good source of adsorbed gas. However, the presence of clay in porous media can also affect shale rock's producibility. The clay-water interaction can cause significant swelling and fines migration in the formation and impair the overall hydrocarbon recovery from shale. Geomaterial micromodels, developed by functionalizing traditional glass or PDMS surfaces with geomaterials (e.g., calcite, Quartz, clay), can represent the physicochemical properties of the natural porous media and have been used during the last few years to understand and evaluate the solid-fluid physicochemical interactions. Here, this study focuses on developing a clay-coated geomaterial surface and investigates the effect of base fluid's salinity on clay adsorption to the glass surface. Glass capillary tubes and straight channel borosilicate glass micromodels are coated with Illite clay to represent the pore-scale clay chemistry of Caney Shale, a Mississippian unconventional play in Southern Oklahoma, USA. 10 wt% of Illite clay slurries made with brines of four different salinities are used to coat the glass-capillary tubes to understand the effect of salinity in clay adsorption on the glass surface and overall coating quality. To achieve a stable coating and evaluate the impact of heat treatment on coating's stability, straight channel glass flow cells coated with Illite clay are heat-treated at low (25 °C) and high temperature (125 °C), and the results are compared with an untreated coated surface. Flooding tests were carried out with brines of different salinities on the coated surface to evaluate the stability of the coating. The experimental results indicated a strong relationship between the brine's salinity and the adsorption of clay particles on the glass surface. An increase in brine concentration resulted in improved adsorption of clay particles on the glass surface. Experiments involving heat treating the glass surface following the coating demonstrated significant improvement in the stability of the coating.

02 PETROLEUM↗

PREPARATION OF CLAY-COATED GLASS SURFACES FOR MICROFLUIDICS-BASED STUDY ON CLAY-FLUID INTERACTIONS

Unique characteristics of clay, even though make any clay abundant shales a good source of adsorbed gas, make them equally difficult to produce. The clay-water interaction can cause significant swelling and fines migration in the shale formation, both of which can significantly impair the overall hydrocarbon recovery from shale. Geomaterial micromodels, developed by functionalizing traditional glass or PDMS surface with geomaterials (e.g., calcite, quartz, clay) can represent the physico-chemical properties of the natural porous media and have been used during the last few years to understand and evaluate the solid-fluid physicochemical interactions. This study focuses on developing a clay-coated geomaterial surface and investigates the effect of base fluid's salinity on clay adsorption to the glass surface. Glass capillary tubes and straight channel borosilicate glass micromodels are coated with Illite clay minerals to represent the pore-scale clay chemistry of Caney Shale, a Mississippian unconventional play in Southern Oklahoma, USA. 10 wt.% of Illite clay slurries made with brines of four different salinities are used to coat the glass-capillary tubes to understand the effect of salinity in clay adsorption on the glass surface and overall coating quality. To achieve a stable coating and evaluate the impact of heat treatment on coating's stability, straight channel glass flow cells coated with Illite clay are heat-treated at low (25 °C) and high temperature (125 °C). Dynamic flooding tests were carried out with brines on the coated surface to evaluate the stability of the coating. The experiments' outcome indicated a strong relationship between the brine's salinity and the adsorption of clay particles on the glass surface. An increase in brine concentration resulted in improved adsorption of clay particles on the glass surface. Experiments involving heat treating the glass surface following the coating demonstrated significant improvement in the stability of the coating.

02 PETROLEUM↗

Solar Thermochemical Carbon Dioxide Splitting Using Ceria and Iron Aluminate Foam Devices and Simulation of a Plant System for Demonstration

An international research project has been undertaken to integrate a unique solar thermal processing reactor system with ceria and iron aluminate as active redox materials for CO2 splitting. Experimental investigations for CO2 splitting were conducted using a solar simulator and tube furnace at Niigata University, followed by demonstrations using a high-flux solar furnace (HFSF) at the National Renewable Energy Laboratory (NREL) in Golden, CO. Each experimental setup consisted of foam devices composed of reticulated porous ceramic (RPC). The RPC has a full ceria or iron aluminate body. It fabricated using the replica method and subjected to a two-step redox reaction, which iteratively separated a stream of CO2 into O2 and CO. Reactivity was evaluated using CO production per mass of the reactive material. The tubular furnace yielded a CO production of 6.41 mL/g at a reduction temperature of 1600degrees C, showcasing a higher CO production rate and total amount than those obtained from experiments conducted with solar simulators and solar furnace setups. For iron aluminate RPC, the productivity was measured as 3.57 mL/g using HFSF at a reduction temperature of 1450degrees C. These results are somewhat higher than those of the previous experiment at lower reduction temperatures of 1400degrees C-1500degrees C. Additionally, the production of CO in the case of ceria RPC was compared with the steady flow model simulation, which assumed chemical equilibrium at various levels of oxygen partial pressure during the reduction process. On the basis of these results, this study proposes a solar fuel system with an open receiver that uses a high-temperature heat transfer fluid.

carbon dioxide thermochemical splitting↗

Micro-continuum approach for mineral precipitation

Abstract Rates and extents of mineral precipitation in porous media are difficult to predict, in part because laboratory experiments are problematic. It is similarly challenging to implement numerical methods that model this process due to the need to dynamically evolve the interface of solid material. We developed a multiphase solver that implements a micro-continuum simulation approach based on the Darcy–Brinkman–Stokes equation to study mineral precipitation. We used the volume-of-fluid technique in sharp interface implementation to capture the propagation of the solid mineral surface. Additionally, we utilize an adaptive mesh refinement method to improve the resolution of near interface simulation domain dynamically. The developed solver was validated against both analytical solution and Arbitrary Lagrangian–Eulerian approach to ensure its accuracy on simulating the propagation of the solid interface. The precipitation of barite (BaSO 4 ) was chosen as a model system to test the solver using variety of simulation parameters: different geometrical constraints, flow conditions, reaction rate and ion diffusion. The growth of a single barite crystal was simulated to demonstrate the solver’s capability to capture the crystal face specific directional growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mobility of soil-biodegradable nanoplastics in unsaturated porous media affected by protein-corona

Soil-biodegradable plastic has been increasingly used as mulches in agriculture, which provides not only agronomical benefits but also in situ disposal and biodegradation options. However, soil-biodegradable plastic mulches inevitably fragment into micro- and nanoplastics during biodegradation, which can reside in soils or migrate into deep soils, where they may not degrade readily due to reduced microbial activity. To date, little is known about the transport of soil-biodegradable micro- and nanoplastics in soils. Here, in this study, we studied the transport of soil-biodegradable nanoplastics (∼200 nm) made of polybutylene adipate co-terephthalate (PBAT) in unsaturated sand (proxy for soil). Specifically, we studied the mobility of pristine and weathered PBAT nanoplastics in the absence and presence of proteins (positively charged lysozyme and negatively charged bovine serum albumin, pH = 7.7). We found that (1) both pristine and the weathered PBAT nanoplastics were mobile; (2) positively charged lysozyme formed protein-coronas around PBAT nanoplastics and inhibited the transport; and (3) decreased water saturation promoted the retention of PBAT nanoplastics via physical straining. These results suggest that soil-biodegradable nanoplastics fragmented from soil-biodegradable plastic mulches are mobile and may readily migrate into deep soil layers, but positively charged proteins and unsaturated flow would prevent such transport via formation of protein-corona and physical straining.

Yu, Yingxue↗

Mechanisms of Zn removal from water by amorphous geopolymer: Molecular-level insights from X-ray absorption spectroscopy, isotope fractionation, and surface complexation modeling

Porous geopolymers have attracted widespread attention as promising heavy metal adsorbents that can be synthesized from aluminosilicate solid wastes. However, the precise microstructural evidence for adsorbed heavy metals on geopolymers remains unclear due to the insensitivity of conventional characterization techniques on minerals with amorphous structure and surface disorder. Batch adsorption and column experiments coupled with X-ray absorption spectroscopy (XAS), Zn stable isotope, and surface complexation model (SCM) were employed to reveal the Zn removal mechanisms with coal fly ash porous geopolymer (CFAPG) at a molecular scale. The macroscopic kinetic and isothermal adsorption of Zn on CFAPG were well described by the pseudo-second-order model and Bi_Langmuir equation, respectively, indicating the presence of abundant heterogeneous active sites on the CFAPG surface. Further, two types of active sites on the CFAPG surface were identified by XAS coupled with Zn isotopes in batch experiments at pH <= 6.0: one is pH-dependent and associated with tetrahedral zinc coordi-nation, and the other is pH-insensitive and associated with octahedral zinc coordination; these sites were confirmed by the bidentate SCM as the variable charge site (surface complexation, >S-OH) and the permanent negative charge site (cation exchange, >X - ), respectively. Furthermore, the important contribution of surface co -precipitation besides surface complexation and cation exchange to the Zn adsorption on CFAPG was identified by XAS coupled with SCM in a flow-through column experiment at pH >6.0. These investigations provide a systemic understanding of the Zn adsorption mechanisms on CFAPG and an SCM reference for the application and prediction of geopolymers in heavy metal-contaminated water remediations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗