Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermodynamic separation”

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.

464 records · Page 26

Tuning Thermal Stability through Dopant Size in Chemically Doped DPP–Thiophene Polymers

Molecular doping of conjugated polymers (CPs) is a key strategy for improving the performance of organic electronics devices, particularly thermoelectrics. Doped donor–acceptor (D–A) conjugated polymers, characterized by a tunable energy gap between the Fermi level and the transport band, show great promise in achieving high electrical conductivity (σ) while preserving a favorable Seebeck coefficient (S). Despite the promising performance enhancement of chemically doped D–A polymers, their thermal stability remains largely underexplored, a crucial consideration for the long-term operation of organic thermoelectric devices. In this study, we investigated the dopant size-dependent thermal stability of a diketopyrrolopyrrole-thiophene (DPP-T) D–A copolymer, utilizing two p-dopants: 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F 4 TCNQ) and Mo(tfd-CO 2 Me) 3 . Temperature-dependent UV–vis–NIR spectroscopy revealed that DPP-T/F 4 TCNQ is more prone to dedoping under a high temperature thermal stress than DPP-T/Mo(tfd-CO 2 Me) 3 . Although the F 4 TCNQ doped polymer shows higher initial in-plane conductivity than its Mo(tfd-CO 2 Me) 3 counterpart, it undergoes a conductivity loss of more than an order of magnitude after annealing at 120 °C for 30 min. In contrast, the in-plane conductivity of DPP-T/Mo(tfd-CO 2 Me) 3 remains stable under the same thermal conditions. Thermogravimetric analysis ruled out dopant sublimation as a primary contributor to dedoping, leading us to attribute the conductivity loss in F 4 TCNQ-doped DPP-T to dopant phase separation and migration. This observation was further confirmed by X-ray scattering studies and nanoscale infrared microscopy and spectroscopy studies. This work could provide further insights into the thermal stability of doped conjugated polymers and suggests that incorporating bulkier dopants is an effective strategy to enhance the thermal robustness of doped DPP-type systems.

Conjugated polymers↗

Hierarchical multi-time-scale predictive thermal management and fuel optimization for heavy-duty compression ignition engines

For heavy-duty diesel engines, NO X emissions reduction is strongly constrained by fuel efficiency. This paper presents a hierarchical model predictive controller (H-MPC) for coordinated control of tailpipe NO X emissions and fuel consumption. The H-MPC uses the separation of slow and fast dynamics that exist in the engine and its aftertreatment system. The controller is synthesized with an architecture in which a high-level MPC uses a longer prediction horizon compared to the low-level predictive controller which tracks the high-level controller command and manages the thermal dynamics of the aftertreatment system. Engine load preview enables the high-level controller to estimate the desired catalyst temperature ahead of time and addresses the selective catalytic reduction (SCR) slow thermal dynamics. Calculated by the high-level controller, the intake manifold pressure, and the start of injection (SOI) crank angle is used as reference trajectories in the low-level controller that regulates fast dynamical behaviors such as engine out NO X emissions. Hardware-in-the-loop (HIL) validation of this integrated H-MPC on a rapid prototype controller shows that when the SCR catalyst temperature is above light-off temperature (warmed-up condition), the engine operation is shifted to operate with the best fuel economy since the warmed-up SCR can efficiently reduce the engine-out NO X emissions. Results indicate that up to 0.8% benefit in cycle averaged BSFC along with a 13% reduction in tailpipe NO X compared to a stock engine calibration can be achieved with the coordinated engine and aftertreatment system through H-MPC.

Engineering↗

Modeling and Optimization of a Rotating Packed Bed Contactor with a Tetraamine-Appended Metal−Organic Framework for CO 2 Capture

A potential contactor technology for sorbent-based CO 2 capture is the rotating packed bed that contains separate sections for continuous adsorption and desorption. A heat exchanger can be embedded to remove heat in the adsorption section and add heat in the desorption section. In this work, we develop a two-dimensional (2D) model of a rotating packed bed for use in CO 2 capture applications. Mass and energy balances for the model are developed based on a Ljungström-type air preheater, which accounts for the counter-current axial flow of gas phases in separate sections of the bed and the rotation of a solid sorbent, which cycles between adsorption and desorption sections. The sorbent used for this analysis is the tetraamine-appended metal−organic framework Mg 2 (dobpdc)(3−4− 3), chosen for its stability and affinity for CO 2 at low partial pressures, such as those from a natural gas power plant source. An optimization problem is solved that considers the trade-off between maximizing the productivity of the bed and minimizing energy consumption. Maximum productivity and minimum energy are found to be 8.53 kg/h/m 3 and 3.84 MJ/kg, respectively, when these objectives are optimized independently. It is observed that the flue gas pressure and bed rotational speed are the desired operating variables to vary for model-based design of experiments to reduce uncertainty in parameter estimation, as these two variables yielded the most information content based on the Fisher information matrix.

20 FOSSIL-FUELED POWER PLANTS↗

Atomistic origin of the entropy of melting from inelastic neutron scattering and machine learned molecular dynamics

The latent heat, L, is central to melting, but its atomic origin remains elusive. It is proportional to the entropy of fusion, ΔS fus = L/T m (T m is the melting temperature), which depends on changes of atom configurations, atom vibrations, and thermal electron excitations. Here, we combine inelastic neutron scattering and machine-learned molecular dynamics to separate ΔS fus into these components for Ge, Si, Bi, Sn, Pb, and Li. When the vibrational entropy of melting, ΔS vib , is zero, ΔS fus ≃ 1.2 kB per atom. This result provides a baseline for ΔSconfig and nearly coincides with “Richard’s Rule” of melting. The ΔSfus deviates from this value for most elements, however, and we show that this deviation originates with extra ΔS vib and extra ΔS config . These two components are correlated for positive and negative deviations from Richard’s rule – the extra ΔS config is consistently ~ 80% of ΔS vib . Our results, interpreted with potential energy landscape theory, imply a correlation between the change in the number of basins and the change in the inverse of their curvature for the melting of pure elements.

36 MATERIALS SCIENCE↗

Reaction path model of the formation of abiotic immiscible hydrocarbon fluids in subducted carbonated serpentinites, Lanzo Massif (Western Italian Alps)

Fluids generated from subducted slabs participate in the cycling of deep carbon in the crust and upper mantle. In these fluids, aqueous carbon species vary in oxidation state between +IV and -IV depending on whether the fluids are oxidizing or reducing, respectively. Most studies of subduction-zone fluids have focused on oxidized carbon species. However, recent studies of natural samples have demonstrated the occurrence of deep, reducing fluids, generated in both subducted oceanic upper mantle and crustal rocks. CH 4 -H 2 -rich fluid inclusions in subducted carbonated serpentinites have demonstrated the existence of abiotic, immiscible, hydrocarbon fluids at upper mantle conditions. To investigate the formation of such immiscible hydrocarbon fluids during the evolution of subducted carbonated serpentinites, we used equilibrium constants from the Deep Earth Water model to carry out predictive chemical mass transfer modeling to simulate the alteration reactions. A novel feature of the models was the inclusion of an immiscible hydrocarbon fluid containing six components (CH 4,f , C 2 H 6,f , C 3 H 8,f , isoC 4 H 10,f , CO 2,f , H 2,f ). This feature enabled prediction of the formation of a separate immiscible fluid in equilibrium with aqueous species and minerals. We developed a predictive reaction path model of invasive H 2,f reacting with carbonated serpentinites and interstitial aqueous fluids for comparison with the natural samples from the Lanzo Massif, western Italian Alps. Over a range of temperatures and pressures, immiscible hydrocarbon fluids formed in association with altered mineral assemblages. Reaction progress caused the transformation of carbonated serpentinites and the formation of clinopyroxene, brucite, graphite, and hydrocarbon fluids, along with changes of pH, logfO 2 , and aqueous species. CH 4,f was the most abundant hydrocarbon species in all the models. The overall results at 2.0 GPa and 400 to 450 °C were consistent with the natural samples from the Lanzo Massif. Interestingly, large amounts of H 2 O formed due to oxidation of H 2 . More hydrocarbons and H 2 O formed in models with lower fluid/rock mass ratios or with more reactant H 2 . Models at different pressure and temperature conditions showed similar results with some variation in the relative stabilities of aragonite, graphite and olivine solid solution, and associated differences in mineral sequences, hydrocarbon fluids, values of aqueous species, and the final log fO 2 and pH. As a result, our models strongly support the laboratory and field evidence that reduction of carbonated serpentinites by infiltrating H 2 fluids can cause the formation of immiscible, abiotic hydrocarbon fluids in subduction zones.

58 GEOSCIENCES↗

Chirality Transfer to the Magnetic Sublattice in the Hybrid Perovskite (R)-/(S)-3-Fluoropyrrolidinium Copper(II) Chloride

Incorporating chiral organic cations into organic–inorganic hybrid materials has been shown to enable the inorganic sublattice to display chiroptical properties. We report a new two-dimensional magnetic (S = 1/2) chiral metal halide perovskite, (R)- and (S)-(C 4 H 9 FN) 2 CuCl 4 (where (C 4 H 9 FN) + is 3-fluoropyrrolidinium), which consists of Cu–Cl inorganic layers separated by (C 4 H 9 FN) + organic cations. The presence of the chiral (C 4 H 9 FN) + organic cation induces the formation of chiral magnetic order, even though the inorganic sublattice itself is nearly structurally centrosymmetric. We also report the racemic variant, containing an equal amount of (R)- and (S)- cations, which shows no evidence of chiral magnetic order. When the magnetic susceptibility is measured perpendicular to the inorganic Cu–Cl layer propagation direction, an antiferromagnetic phase transition at Néel temperature T N = 2.23 K is observed in both the chiral and racemic materials, and the existence of the magnetic phase transition is supported by specific heat capacity measurements. Field-induced magnetic chirality is observed through the existence of a second-order magnetoelectric effect in the chiral variant, while no magnetoelectric signal is observed for the racemic material, indicating the absence of magnetic chirality. Our findings demonstrate that materials exhibiting chiral magnetic order can be created through the incorporation of a chiral cation into an organic–inorganic hybrid magnetic material, potentially allowing for the design of tailored materials that combine chiral magnetism with other desirable optical and electronic properties that come from structural chirality.

Cations↗

Insight into pseudo-film-boiling in supercritical water under simplified interfacial conditions

The focus of this paper is the physics of supercritical water film boiling through the simulation of natural convection of supercritical water in a rectangular cavity. It has been widely believed that above a fluid's critical point, only one phase exists. However, recent research has indicated the existence of distinct gas-like and liquid-like regions separated by the Widom line, which is the locus of the maxima of the specific heat. Along the Widom line, density decreases sixfold, viscosity drops by a factor of two, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film-boiling, denoted here as pseudo-two-phase flow dynamics. When the cavity is heated from the bottom, a pseudo-film of gas-like fluid is formed at the heated wall, and the flow complexity evolves, characterized by the formation of multiple vortices. Results demonstrate the ability of modern conjugate heat transfer computational fluid dynamics to capture the pseudo two-phase regime, with the formation of a pseudo-film at the heated wall. Further, the boiling curve is studied, showing an increase in heat flux with an increase in the temperature differential between the wall and the pseudo-critical temperature. Furthermore, results also show a previously uncharacterized instability, where the convection cells in the cavity oscillate and can flip directions under certain conditions.

Computational fluid dynamics↗

Laser Measurements and Modeling of Shock Tunnel Freestream Velocity and Multispecies Thermal Nonequilibrium

Coherent anti-Stokes Raman scattering (CARS) and nitric oxide molecular tagging velocimetry (NO-MTV) were used to characterize the freestream in Sandia’s Hypersonic Shock Tunnel (HST) using a burst-mode laser operated at 100 kHz. Experiments were performed at nominal freestream velocities of 3 and 4 km/s using both air and N 2 test gas. The CARS diagnostic provides nonequilibrium characterization of the flow by measuring vibrational and rotational temperatures of N 2 and O 2 , which are compared to NO temperatures from separate laser absorption experiments. Simultaneous, collinear freestream velocities were measured using NO MTV along with pitot pressure measurements. Furthermore, this extensive freestream dataset is compared to nonequilibrium CFD that is capable of modeling species-specific, vibrational temperatures throughout the nozzle expansion. Significant nonequilibrium between vibrational and rotational temperatures was measured at each flow condition. N 2 exhibits the most nonequilibrium, followed by O 2 and NO. The CFD model captures this trend, although it consistently overpredicts N 2 and O 2 vibrational temperatures. At 3 km/s, the modeled NO nonequilibrium is underpredicted, whereas it is overpredicted at 4 km/s. Good agreement is seen between CFD and the velocity and rotational temperature measurements. Experiments with water added to the test gas yielded no discernable difference in vibrational relaxation.

Aerodynamics↗

Mechanisms enabling reconfigurability and long-term retention in vanadium oxide electrochemical memory

Phase coexistence in nanoscale electrochemical random-access memory (ECRAM) has recently been demonstrated to enable both information storage and extraordinary reconfigurability. These proof-of-principle demonstrations have left the mechanistic details of such a process unresolved. Particularly, the mechanisms that stabilize the multiple phases, and the underlying processes behind sustained memory retention, remain unclear, and are necessary to design such devices. Here we report microscale ECRAM devices composed of V⁢O𝑥, which enables us to directly probe the active region in an operando fashion using optical techniques. Using Raman mapping, we show the phase coexistence driven by the electrochemical injection of O vacancies to be spatially uniform (i.e., with no filaments). The stability was observed to be unusually long, with 1% loss over 14 years in ambient conditions. First-principles calculations of the oxygen vacancy formation energies in V⁢O 𝑥 further support the thermodynamic coexistence of multiple V⁢O 𝑥 phases and clarify the origin of the observed long-term retention in the ECRAM devices. Further, we demonstrate single devices that can be voltage programmed to exhibit synaptic, neuronal, and reconfigurable logic gate functionalities. Furthermore, we not only uncover the phase coexistence mechanism that may help device design, but also demonstrate the circuit-level applications of reconfigurability.

Electrical conductivity↗

Understanding Differences in Water Adsorption Isotherms: Structural Variations, Force Fields, and Monte Carlo Simulation Approaches

Accurate prediction of water adsorption in micro- and mesoporous materials with hydrophobic pores is essential for the design and characterization of advanced adsorbent materials for separation and energy applications. Here, we assess the reproducibility and consistency of water adsorption isotherms in two microporous all-silica MFI zeolite structures (MFI-K and MFI-O) using two different zeolite force fields and three simulation approaches: grand canonical Monte Carlo (GCMC), Gibbs ensemble Monte Carlo (GEMC), and transition matrix Monte Carlo (TMMC). We demonstrate that consistent treatment of the bulk fluid phase in GCMC and TMMC simulations is critical for reconciling isotherms across methods, and we construct simulation-based equations of state for the TIP4P water model to enable rigorous fugacity-to-pressure conversions. Large shifts in the isotherms are observed for two zeolite force fields developed using different parametrization strategies, with the GCS force field representing implicitly a defect-containing all-silica zeolite, whereas the TraPPE-zeo force field accurately represents an essentially defect-free all-silica zeolite. While water in the van Koningsveld structure of MFI exhibits a first-order phase transition and condensation-like step for adsorption near room temperature, water in the Olson structure of MFI displays continuous adsorption, attributed to differences in the adsorption free energy landscapes. Structural analysis reveals that small geometric variations, particularly Si–O–Si bond angles near the strongest adsorption sites, lead to these substantial differences in adsorption behavior. Furthermore, our results highlight the sensitivity of simulated water adsorption isotherms in hydrophobic frameworks to seemingly small differences in the framework structures, force field parametrization, and simulation approaches.

36 MATERIALS SCIENCE↗

Toward Hydrogen Isotope Separations through Strong Hydrogen Adsorption at Open Copper(I) Sites in an Ultramicroporous Metal-Organic Framework

Metal-organic frameworks with coordinatively unsaturated metal sites (open metal sites) capable of engaging in orbital interactions with pi-acidic gases are of interest for enabling ambient-temperature gas separations, such as hydrogen isotope separations. In view of the weakly pi-acidic nature of H2, we sought to strengthen pi-backbonding-mediated H2 adsorption through pore confinement effects. Toward that end, we synthesized and characterized the ultramicroporous metal-organic framework CuxZn5-xCl4-yHz(bbta)3 (CuIZn-MFU-4; H2bbta = 1H,5H-benzo(1,2-d:4,5-d')bistriazole), featuring pi-basic trigonal pyramidal CuI sites that reside within 7 A of one another at their closest. Gas adsorption measurements reveal an H2 adsorption enthalpy of -38 kJ/mol, exceeding that of the larger-pore analog (CuIZn-MFU-4l; -33 kJ/mol) and representing the strongest H2 adsorption yet achieved in a metal-organic framework. The stronger H2 adsorption in CuIZn-MFU-4 is attributed to a combination of pore confinement effects and the increased ..sigma..-accepting nature of the CuI sites caused by a more electron-withdrawing bbta2- linker, as supported by structural, spectroscopic, and computational evidence. With the strongest H2 adsorption, equilibrium isotope effects in CuIZn-MFU-4 lead to a D2/H2 selectivity (as estimated by ideal adsorbed solution theory) of 1.35 even at 298 K, approaching the values reported below 200 K for conventional porous materials.

08 HYDROGEN↗

Progress Report on Performance of A709 and G91 Steels in Sodium

Specimens in six different processing and heat treatment conditions of A709 H58776 were exposed to sodium at 550, 600, and 650°C, respectively for various exposure times. G91 base metal was tested in sodium at 550, 600, and 650°C, and G91 weldment at 550 and 600°C, respectively. Parallel thermal aging experiments were conducted on G91 and A709 to obtain thermal aging data for comparison with sodium exposure results to separate the thermal and sodium effects. The corrosion data obtained on A709 steel continue to show low corrosion rates, which indicate good compatibility of A709 with sodium when oxygen content is controlled. G91 has also shown acceptable corrosion rates over the temperature range of 550-650°C in sodium. Thermal aging or sodium exposures reduced the tensile strength, uniform elongation, and total elongation for A709 ESR, AOD specimens. These specimens consistently showed higher yield stress and tensile strength and lower uniform and total elongations after sodium exposure than after thermal aging under comparable testing conditions. It implies that there was an additional effect of sodium exposure that gave rise to an increase in tensile strength of A709. Dynamic strain aginginduced flow serrations were completely removed by either thermal aging or sodium exposure at 650°C, which imply that carbon or nitrogen in solution was removed from the solution and formed precipitates during aging or sodium exposures. The tensile data ruled out the possibility of decarburization under the sodium exposure conditions. The HOMO specimens behaved somewhat differently from the ESR and AOD specimens in sodium. While thermal aging and sodium exposures at 550 and 600°C decreased the yield stress and the ultimate tensile strength of G91, there was virtually no additional effect resulting from sodium exposures at these two temperatures. In contrast, sodium exposures at 650°C had a drastic effect on the yield stress and the ultimate tensile strength of G91. It is suggested that G91 experienced decarburization in the 650°C sodium environments. Carbon concentrations in sodium in the SMT-1 and SMT-2 loops were determined by a foil equilibration method. The estimated carbon concentration was in the range of 0.8-1.2 ppm in the SMT-1 loop and 0.3-0.7 ppm in the SMT-2 loop. The carbon activity in sodium in the SMT-2 loop was estimated to be 0.03-0.08 at 600°C, and 0.08-0.2 at 550°C. The carbon activity in sodium at 650°C in the SMT-1 loop was 0.04-0.07. Equilibrium simulation of the carburization – decarburization processes was conducted for A709 and G91 steels exposed in sodium environments at temperatures of 550-700°C. The carbon activity-concentration relationship for G91 was re-evaluated by considering four phases in G91, i.e. bcc ferrite, M 23 C 6 , NbC and VC carbides. It was found that the carburization-decarburization process in G91 steel was dictated by M 23 C 6 carbides at high carbon activities, while NbC carbides dominated the process at low carbon activities. Formation of NbC increases the decarburization resistance of G91 steel. It remains to be understood whether the decarburization resistance provided by MC carbides in G91 can be maintained during long-term operations of SFRs. The carbon concentration-activity relationship for A709 was calculated based on the equilibrium of fcc-austenite and M 23 C 6 carbide phase. The calculations showed that A709 decarburizes at 650°C in the SMT-1 loop environments, which is different from the experimental findings. The carbon concentration-activity relationship in A709 was further evaluated by considering four phases, including fcc-austenite, M 23 C 6 , NbC, and TiC carbide phases. TiC is the most stable carbides among the three carbide phases. The carburization-decarburization process in A709 is dictated by M 23 C 6 carbides at high carbon activities and by TiC at low carbon activities. Formation of TiC increases the decarburization resistance of A709 steel. Because of the complexity of the precipitation process in A709 during sodium exposures, detailed characterization of precipitates in sodium-exposed specimens and improved thermodynamic models are needed to understand the carburization-decarburization behavior of A709. Kinetic analysis of carbon transfer will be investigated for G91 and A709 in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Mixed Nd 4+/3+ and Cluster Magnetism in Hexagonal Perovskite 12R-Ba 4 NdMn 3 O 12−γ

Hexagonal perovskite oxides with 12R stacking host well-separated face-sharing octahedral metal trimers with short metal–metal distances, leading in some cases to large degrees of magnetic frustration and cluster magnetism. Introducing magnetic ions in proximity to these trimers can influence the degree of frustration and cluster magnetism, but the extent to which the magnetism can be tuned by varying neighboring metal cations remains an open question. In this work, we test the impacts of using Nd in proximity to Mn trimers in the hexagonal perovskite 12R-Ba 4 NdMn 3 O 12−γ . Given this stoichiometry, Nd should assume the 4+ oxidation state with a spin state of S = 1, which would be the first realization of Nd 4+ in an oxide environment. Through detailed bulk magnetic, X-ray absorption spectroscopic, and powder neutron diffraction (PND) measurements, we find that Nd 4+ is realized in this material, but there is also partial reduction to Nd 3+ , which is charge-balanced by O vacancies. Magnetometry measurements indicate an antiferromagnetic ordering temperature T N ≈ 16 K, and PND measurements reveal a surprising collinear antiferromagnetic structure with magnetic space group Pc2/m (no. 10.49 in BNS notation), which has not previously been seen in this class of materials. Our results represent a comprehensive analysis of the structural, electronic, and magnetic properties of 12R-Ba 4 NdMn 3 O 12−γ , showing the first observation of partial Nd 4+ in an oxide, and demonstrating that this structural class can host a broad range of magnetic structures which are not easily predicted based on compositional trends.

14 SOLAR ENERGY↗

Engineering and Application of a Thermostable MHETase for PET Depolymerization

Enzymatic hydrolysis of poly(ethylene terephthalate) (PET) releases mono(2-hydroxyethyl) terephthalate (MHET) as a major product, the accumulation of which can prolong reactor residence times and complicate downstream monomer separations. The use of a MHETase enzyme can enable MHET hydrolysis to the monomers, terephthalic acid and ethylene glycol, but industrial PETases typically operate at thermophilic temperatures and the well-known MHETase from Ideonella sakaiensis is a mesophilic enzyme, thus warranting the development of thermophilic MHETases. Here, we characterize thermostable MHET-active enzymes from a natural diversity screen by applying a hidden Markov model based on the previously reported, archaeal ferulic acid esterase, PET46. We identified enzymes with higher thermostability than PET46 and quantified their MHETase activity in reactions at 70 °C. The crystal structure of MHT077, the homologue with the highest MHETase activity and an apparent melting temperature (T m,app ) of 94.6 °C, informed site saturation mutagenesis in the active site and lid-domain interface. MHT077 exhibited a ∼100-fold slower unfolding rate at 65 °C than PET46, indicating substantially greater kinetic stability. In parallel, we applied evolution-informed design, a probabilistic model that leverages coevolutionary patterns in large multiple sequence alignments, to improve the activity and thermostability of five ferulic acid esterases. One design, EV-MHT043–5 was identified with a comparable thermostability (T m,app = 96.1 °C) and a 3-fold improvement in its MHETase activity relative to the wildtype enzyme, MHT043. Combination variants of beneficial mutations were screened and afforded a variant, MHT077 LFK , which reduced MHET accumulation in bioreactor experiments with postconsumer PET waste. Overall, this study expands the known MHET-hydrolyzing protein scaffolds available for enzymatic PET recycling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗