Search NASA⌕ Search

SEARCH · Search NASA

Results for “Enthalpy method”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A comprehensive modeling of falling film evaporators subject to vapor flow, pass arrangements, and refrigerants

Improving the heat transfer performance of falling film evaporators is a crucial step for improving the energy efficiency of the heat pump or refrigeration systems. This study conducts a numerical investigation for the practical-scale falling film evaporator based on the epsilon–NTU method with an updated heat transfer correlation. The algorithm was validated with lab-scale and real-scale falling film evaporator experimental results, and the prediction reaches a mean absolute deviation of 12.5%. Here, the parametric study encompasses eight refrigerants: R-134a, R-410A, R-600a, R-717, R-1270, R-152a, R-1234yf, and R-1234ze(E). The results indicate that vaporization enthalpy of refrigerant is a key property in selecting an appropriate working fluid because it helps minimize severe heat transfer degradation caused by dry-out. Additionally, the vapor-flow–induced heat transfer degradation can be predicted using the critical Weber number. Furthermore, the trade-off between extending the tube length and increasing the number of tubes for heat transfer improvement is discussed. Finally, different two-pass arrangements show deviations of less than 4 %.

Dry-out↗
DOE OSTI2023

Thermal Conductivity of a Novel Salt-Hydrogel Complex as a Phase Change Material for Building Thermal Management

Use of electricity for space heating, ventilation and air conditioning in buildings may be better managed through temporary energy storage. Compared with electrochemical storage, the levelized cost of thermal energy storage systems can be lower depending on geographical location and thermophysical properties [1]. In any phase change material (PCM), the rate of absorption and release of thermal energy increases with higher thermal conductivity [2]. However, leading PCMs like paraffin waxes have thermal conductivities on the order of 0.1 W/mK. Glauber’s salt, a salt hydrate, possesses a relatively high thermal conductivity but suffers from supercooling and phase segregation. Research has shown that polymer-based composites with salt hydrates can prevent phase segregation and limit supercooling over hundreds of cycles [3] though the thermal conductivity of the composite may be reduced leading to low energy storage rates. Here, we report thermal conductivity measurements on polymer hydrogels which are stable in solutions of Glauber’s salt. To measure thermal conductivity of soft gels, we use a 3-omega method [4] where heat flows bidirectionally through a substrate below and through a gel placed on top. This suits measurements of both gels and liquids unlike transient plane heat source or transient hot wire methods. We report measurements of thermal conductivity of salt-hydrogel complexes as a function of temperature, cross-linking and salt concentration. In addition, we present modeling of thermal conductivity of sodium sulfate solutions with available theories. We finally employ polarized optical microscopy to observe inside the transparent hydrogel networks, showing that precipitating crystals may influence thermal conductivity measurements for high salt concentrations. Our thermal conductivity results combined with enthalpy of phase change, mass density and specific heat are essential to accurately design thermal storage systems for energy-efficient buildings.

composite, thermal conductivity↗
DOE OSTI2023

Study of unripe and inferior banana flours pre-gelatinized by four different physical methods

This study aimed to prepare the pre-gelatinized banana flours and compare the effects of four physical treatment methods (autoclaving, microwave, ultrasound, and heat-moisture) on the digestive and structural characteristics of unripe and inferior banana flours. After the four physical treatments, the resistant starch (RS) content values of unripe and inferior banana flours were decreased from 96.85% (RS2) to 28.99–48.37% (RS2 + RS3), while C∞ and k values were increased from 5.90% and 0.039 min −1 to 56.22–74.58% and 0.040–0.059 min −1 , respectively. The gelatinization enthalpy (ΔHg) and I 1047/1022 ratio (short-range ordered crystalline structures) were decreased from 15.19 J/g and 1.0139 to 12.01–13.72 J/g, 0.9275–0.9811, respectively. The relative crystallinity decreased from 36.25% to 21.69–26.30%, and the XRD patterns of ultrasound (UT) and heat-moisture (HMT) treatment flours maintained the C-type, but those samples pre-gelatinized by autoclave (AT) and microwave (MT) treatment were changed to C + V-type, and heat-moisture (HMT) treatment was changed to A-type. The surface of pre-gelatinized samples was rough, and MT and HMT showed large amorphous holes. The above changes in structure further confirmed the results of digestibility. According to the experimental results, UT was more suitable for processing unripe and inferior banana flours as UT had a higher RS content and thermal gelatinization temperatures, a lower degree and rate of hydrolysis, and a more crystalline structure. The study can provide a theoretical basis for developing and utilizing unripe and inferior banana flours.

Zhang, Siwei↗
DOE OSTI2023

Structurally Constrained Evolutionary Algorithm for the Discovery and Design of Metastable Phases

Metastable materials are abundant in nature and technology, showcasing remarkable properties that inspire innovative materials design. However, traditional crystal structure prediction methods, which rely solely on energetic factors to determine a structure’s fitness, are not suitable for predicting the vast number of potentially synthesizable phases that represent a local minimum corresponding to a state in thermodynamic equilibrium. Here, we present a new approach for the prediction of metastable phases with specific structural features, and interface this method with the XTALOPT evolutionary algorithm. Our method relies on structural features that include the local crystalline order (e.g., the coordination number or chemical environment), and symmetry (e.g., Bravais lattice and space group) to filter the breeding pool of an evolutionary crystal structure search. The effectiveness of this approach is benchmarked on three known metastable systems: XeN 8 , with a two-dimensional polymeric nitrogen sublattice, brookite TiO 2 , and a high pressure BaH 4 phase that was recently characterized. Additionally, a newly predicted metastable melaminate salt, P1¯WC 3 N 6 , was found to possess an energy that is lower than two phases proposed in a recent computational study. Here, the method presented here could help in identifying the structures of compounds that have already been synthesized, and developing new synthesis targets with desired properties.

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

Gas‐Phase Ion‐Molecule Interactions in a Collision Reaction Cell with ICP‐MS / MS : Investigations with CO 2 as the Reaction Gas

Carbon dioxide (CO 2 ) was used as a reaction gas to investigate the gas‐phase ion‐molecule interactions using the Agilent 8900 ICP‐MS/MS. A solution containing forty‐five elements representative of the periodic table was used to supply the ions to react with CO 2 in the collision/reaction cell (CRC). The only significant product ions formed were monoxides. The general reactivity was shown to be consistent with density functional theory (DFT)‐predicted reaction enthalpies, such that all predicted exothermic reactions produced product ions at levels of at least 1% of the unreacted ion. Most endothermic reactions observed had sufficient kinetic energy in excess of the reaction enthalpies. Our results suggest that reaction enthalpy is a reasonable predictor of reactivity with CO 2 on the timescales of the interactions in non‐thermal ICP‐MS/MS systems. The ease and rapidity of data collection with the ICP‐MS/MS and DFT calculations using the NWChem suite has value given the scarcity of thermochemical data of CO 2 reactions in the literature. These studies are especially useful for the identification of targeted reaction chemistries to be leveraged for analytical method development, such as for the inline separation of isobaric interferences from analytes of interest.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The structure of CaO–MgO–Al 2 O 3 –SiO 2 melts and glasses doped with FeO X –NiO

Neutron and x-ray diffraction measurements have been performed on CaO–MgO–Al 2 O 3 –SiO 2 (CMAS) glasses doped with NiO–Fe X O at room temperature, along with x-ray measurements on aerodynamically levitated liquids at ≥2000 K. The disordered structures have been modeled using empirical potential structure refinement to investigate the relation between the aluminosilicate network and the modifying cations. The SiO 4 and AlO 4 tetrahedra are found to have wider Si–O and Al–O bond distance distributions in the glass, and the first Ca–O n coordination shell is highly distorted, redistributing different populations of long and short bonds between the liquid and the glass. The addition of Fe and Ni at low aluminosilicate content increases the number of free oxygens not bonded to AlO 4 or SiO 4 . Mg–O and Fe–O are both found to be predominantly fourfold and fivefold in the liquid and glassy states. Despite these low coordination numbers, their bond angle distributions indicate that they are predominantly in nontetrahedral-type geometries, with ferrous and ferric iron possessing similar coordination environments. The Ca–O and Mg–O average coordination numbers and enthalpies of solution are consistent with their higher reactivity within relatively acidic aluminosilicate melts.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Uncertainty Quantification and Error Propagation in the Enthalpy and Entropy of Surface Reactions Arising from a Single DFT Functional

Accounting for parametric uncertainty in models is essential for quantifying the models’ predictive ability. Recently, approaches have been introduced to estimate parametric uncertainty in kinetic models while accounting for correlations among energy parameters. However, correlations have been estimated indirectly and correlations in entropies have not been accounted for. For surface-catalyzed microkinetic models of >C2 (more than two carbon-containing) molecules, which consist of thousands of reaction steps and intermediate surface species, first-principles density functional theory (DFT) is costly, and thus, estimation of thermochemistry and reaction barriers requires surrogate methods of DFT, such as group additivity and Brønsted–Evans–Polanyi relationships, respectively. For such parametrization, model uncertainty is unclear. This work develops a framework to overcome these gaps using group additivity and a single DFT functional. We estimate correlations in parameters of kinetic models and quantify uncertainty for thermochemistry, reaction barriers, reaction paths, and ultimately reaction rates, accounting also for the contribution of entropic uncertainty. Furthermore, the approach is illustrated on propane combustion and ethane oxidative dehydrogenation reactions.

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

Optimal Design and Techno-Economic Analysis of 3D-Printed, Intensified Packings for Absorbers and Strippers in Solvent-Based CO 2 Capture

A potential technology for the CO 2 absorption process is utilizing intensified structured packing with embedded cooling/heating channels for continuous heat exchange, which can overcome limitations of discrete methods, such as discrete intercooling and centralized reboilers, to aid in reducing energy consumption and decreasing costs. This work investigates the modeling of intensified packing (IP) for the stripper tower, extending on previous work for the absorber, which distributes heat internally within the column, improving the thermodynamics for the solvent regeneration process. The model includes submodels for steam turbine extraction to produce steam at various qualities as well as a surrogate model for calculating steam enthalpy. A cost model for a plant-scale absorption capture process was developed, allowing for the design of the plant to be optimized, subject to minimizing capture cost using two different power plant flue gas sources. In this optimization, the placement of IP in both towers is optimized to balance the trade-off between enhanced heat transfer and reduced mass transfer volume. For natural gas combined cycle flue gas, the standard process configuration had a minimum cost of $\$$65.40/tonne CO 2 , and considering IP, the minimum capture cost is reduced to $\$$62.73/tonne, with utilization in the stripper column, which reduces yearly costs by up to $\$$2.67 MM/yr. Cooling the absorber through IP, or intercoolers, was only found to be beneficial at higher capture rates, with IP in both towers having a cost of capture of $\$$68.08/tonne at 99.9% capture, a reduction of $\$$12.64/tonne when using only intercoolers at the same capture rate. When capturing from pulverized-coal power plants, the minimum cost of capture when using IP in both towers is $\$$44.18/tonne (at 97% capture), while the standard configuration with and without intercoolers was $\$$45.69 and $\$$47.22 per tonne, respectively. This results in a reduction in yearly costs of $\$$16.98 MM/yr from the base-case configuration. At this higher CO 2 concentration, cooling in the absorber from the IP becomes extremely beneficial, reducing energy consumption by up to 6%.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2026

Numerical implementation of the MOOSE subchannel module (SCM) algorithm

The subchannel module of moose (SCM), previously referred to as Pronghorn-SC, is a subchannel code designed to resolve single-phase flow fields and calculate the relevant flow variables, in nuclear reactor fuel-pin assemblies. The assemblies it models are: water-cooled, with bare fuel-pins in a square ducted, quadrilateral lattice and liquid-metal/water-cooled, with bare/wire-wrapped fuel pins in a hexagonal ducted, triangular lattice. Previous publications have presented the development, validation and verification of SCM. This work presents an comparative overview of the different solvers implemented within SCM. This includes the solver algorithm description, and a performance comparison. The two test cases chosen to demonstrate the solver performance were taken from the PSBT enthalpy mixing benchmark and the ORNL-19 pin benchmark.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Heterostructural Alloy Phase Diagram for (Cd 1-x Zn x ) 3 As 2

Alloying the topological semimetal Cd 3 As 2 with Zn 3 As 2 provides a potential route for controlling the electronic properties. We predict the alloy phase diagram from first-principles calculations, considering that both end members have a crystal structure derived from the antifluorite lattice, but with different arrangements of the unoccupied cation sites. To overcome the limitations of the regular solution approximation and to include short-range order effects, we perform Monte Carlo simulations, parameterize the temperature dependence of the mixing enthalpy ΔH m , and perform thermodynamic integration of the free energy. The resulting phase diagram exhibits features that are unique to heterostructural alloy systems and provides computational predictions of solubility limits and composition ranges that are stable against spinodal decomposition.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Computing the solubility of argon and xenon in molten sodium chloride and potassium chloride salts

Molten salt reactors (MSRs) offer significant advancements in nuclear reactor safety and efficiency by operating at higher temperatures and lower pressures compared to traditional reactors. A critical aspect of MSR operation involves understanding the solubility of fission byproducts, particularly noble gases, in the molten salts used. This study employs molecular dynamics (MD) simulations to compute Henry’s law constants and enthalpies of solvation for argon and xenon in molten sodium chloride (NaCl) and potassium chloride (KCl). We developed a new pairwise potential for the noble gas and salt interactions based on first principles calculations. We then used this potential to calculate Henry’s law constants of the two gases in the molten salts, which were modeled using both a rigid ion model (RIM) and a polarizable ion model (PIM). The solubility calculations, performed using the Widom insertion method, show qualitative agreement with limited experimental data, highlighting the temperature dependence and greater solubility of both gases in KCl compared to NaCl. Additionally, free volume analysis elucidated the role of available space within the molten salts in governing solubility trends. Our findings suggest that PIM trajectories provide more reliable predictions for noble gas solubility than RIM due to their accurate density representation. Furthermore, these results enhance understanding of gas solubility in MSR environments, and the methods can be readily extended to other systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Gas-phase ion-molecule interactions in a collision reaction cell with triple quadrupole-inductively coupled plasma mass spectrometry: Investigations with N 2 O as the reaction gas

Nitrous oxide (N 2 O) was used as a reaction gas to investigate the gas phase ion-molecule interactions using the Agilent 8900 QQQ-ICP-MS. A multi-element standard containing 45 elements with masses ranging from 9 to 208 u was measured in the presence and absence of N 2 O. The main product ion species observed were oxides and nitrides. Comparison of the N 2 O reaction results with similar measurements conducted with O 2 revealed that N 2 O was more effective at forming oxides in general: the elements Cd and Pb were shown to produce oxides with N 2 O where the reaction did not occur with O 2 . Nitrous oxide was also shown to produce a significant amount of nitride species in a few cases. The general reactivity was shown to be consistent with density functional theory (DFT)-predicted reaction enthalpies, such that all predicted exothermic reactions produced product ions at levels at least 1% of the unreacted ion. Our results show that reaction enthalpy is a reasonable predictor of reactivity with N2O on the timescales of the interactions in non-thermal ICP-MS/MS systems. Our work demonstrates the utility of two relatively new platforms (commercial elemental ICP-MS/MS and EMSL Arrows interface to the NWChem program suite), which allows for the study of a large number of elements within a short period. While DFT with the basis sets utilized here is not the most accurate computational method, it is also not computationally expensive and is shown to be suitable for predicting gas phase reactivity in the QQQ-ICP-MS for the majority of ions studied. Here, the ease and rapidity of data collection and DFT calculations has the potential to be very impactful for the identification of targeted reaction chemistries to be leveraged for analytical method development, such as for the inline separation of isobaric interferences from analytes of interest.

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

Dynamic Thermal Performance Analysis of PCM Products Used for Energy Efficiency and Internal Climate Control in Buildings

PCMs are attractive for the future generation of buildings, where energy efficiency targets and thermal comfort expectations are increasingly prioritized. Experimental analysis of local thermal processes in these dynamic components and whole-building energy consumption predictions are essential for the proper implementation of PCMs in buildings. This paper discusses the experimental analysis of the thermophysical characteristics of both a latent heat storage material (PCM) and a product containing this PCM. The prototype product under investigation is a panelized PCM technology containing inorganic, salt-hydrate-based PCM. The thermal analysis includes studies of melting and freezing temperatures, enthalpy changes during phase change processes, nucleation intensity, sub-cooling effects, and PCM stability. The PCM’s stability is also investigated, as is the ability of PCM products to control local temperatures and peak load transmission times. Two inorganic PCM formulations based on calcium chloride hexahydrate (CaCl 2 .6H 2 O) were prepared and tested in laboratory conditions. Material-scale testing results were compared with outcomes from the system-scale analysis, using both laboratory test methods as well as field exposure in test huts. This work demonstrates that PCM technologies used in buildings can effectively control both the magnitude of thermal storage capacity as well as the time of the peak thermal load. It was found that commonly used material-scale testing methods may not always be beneficial in assessing the dynamic thermal performance characteristics of building technologies containing PCMs.

42 ENGINEERING↗
DOE OSTI2023

Validation of Octane Hyperboosting Phenomenon in Prenol and Structurally Related Olefinic Alcohols

Hyperboosting is a recently discovered phenomenon in which the research octane number (RON) of a blend is higher than both the neat blendstock and the neat fuel it was blended into. RON is a measure of a fuel's resistance to knock, and knock is a cause of engine inefficiency. Blends which exhibit hyperboosting are blends in which an overall improvement in engine efficiency may be expected. The first discovery of hyperboosting came from blending experiments in which prenol was blended into several different base fuels. Here, ignition delay time (IDT) measurements taken using a commercially available constant volume combustion chamber called the Advanced Fuel Ignition Delay Analyzer (AFIDA) are presented. The data show that some prenol blends have longer IDTs (lower reactivity) than either neat prenol or the base fuel, providing further evidence of hyperboosting. Additionally, more blending data is presented in which the base fuel is varied, which allows for a better understanding of hyperboosting sensitivity to chemical classes. The data indicate that aromatics may inhibit, and branched alkanes may enhance the magnitude of hyperboosting observed. Enthalpy of vaporization estimates are also given for several molecules which are blended into a 4-component surrogate. These estimates are derived from Equation of State simulations and reveal that there is no observable correlation between hyperboosting and enthalpy of vaporization. Blending data for molecules which share structural similarities with prenol are also presented. Structure property relationships are suggested, in which the double bond motif of prenol seems to play an important role in hyperboosting. This may help to understand hyperboosting and its underlying mechanism. Lastly, blending curves of surrogate blends with prenol experienced hyperboosting under lean (Homogeneous Charged Compression Ignition-HCCI) operating conditions, which validates that hyperboosting is not an artifact of the octane test methods, but inherent to the properties of prenol.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Wood template-supported phase change material composites for durable and form-stable thermal energy storage in buildings

Here, to reduce and shift peak energy loads in buildings, phase change materials (PCMs) with high transition enthalpies and transition temperatures near human thermal comfort are desirable for thermal energy storage (TES). Traditional solid/liquid PCMs suffer from leakage during thermal cycling, requiring encapsulation that lowers heat storage capacity. Wood templates (WTs), with porous and hierarchical structures, provide natural encapsulation scaffolds for PCM containment. Polyethylene glycol (PEG) is a compatible PCM with WTs, but when infiltrated alone, ~30 wt% of PEG binds to wood with no detectable phase change, limiting TES efficiency. To address this, we developed a method to (1) improve the form stability of balsa and pine-based composites (BWT + PCMs and PWT + PCMs) and (2) reduce inactive PCM within the composites to ~10 wt%. BWT + PCMs exhibit transition properties upwards of 114.2 J/g at 25.4 degrees C, with no degradation after 1000 thermal cycles, and similar stiffness compared to raw balsa. Meanwhile, PWT + PCMs exhibit 25 % higher storage efficiency with the addition of poly(ethylene glycol) diacrylate compared to solely PEG-infiltrated PWT. PWT + PCM retains 84 % of raw pine's mechanical stiffness, sufficient for light-duty construction. Our shape-stabilized WT + PCM composites enrich the functionality of wood materials as both ideal TES material candidates and light-duty building construction applications.

25 ENERGY STORAGE↗
DOE OSTI2020

Extracting meaningful standard enthalpies and entropies of activation for surface reactions from kinetic rates

While analyses based on the Arrhenius kinetic model have been successfully applied since its introduction, the Arrhenius model neglects to describe pre-exponential factor in a scientifically meaningful fashion. Since the 1930s, transition-state theory (TST), has met with success in interpreting the pre-exponential factor’s value, allowing a standard entropy of activation to be estimated. However, analyses based on TST’s assumptions have been applied inconsistently in the literature, particularly in corrosion science, leading to difficulty in comparison of standard entropy of activations from different studies. In this work, the foundational principles of TST and standard states are discussed and a standard method to apply TST in analyzing rates of surface reactions is recommended. When full details are not available, reacting species’ concentrations should be normalized to the concentration of active surface sites. For corrosion reactions, conversion relationships are given to convert from units of corrosion rate to surface reaction rate, consistent with TST. This method is dubbed a surface reactant equi-density approximation. Application of this standard to reported data results in adjustments of the standard entropy of activation between - 65 and +50 J/mol K and brings reported entropies into a narrower range of values.

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

Grain boundary segregation in BCC vanadium-based alloys: Quantum-accurate computed segregation spectra and targeted experimental validations

Grain boundaries are critically important to the material performance of fusion reactor materials such as vanadium, particularly mechanical properties and irradiation resistance. A key challenge to the design and control of grain boundaries in vanadium alloys is the lack of quantitative data on grain boundary segregation. In this study, we combine computational and experimental methods to address this gap. Furthermore, using a machine learning-accelerated quantum mechanics/molecular mechanics approach, we calculated the segregation spectra for 28 transition metal elements in polycrystalline vanadium, and validated these predictions experimentally for a subset of solutes that sample a range of segregation behavior, specifically zirconium, titanium, and tungsten, using analytical transmission electron microscopy. Furthermore, the agreement between experiment and theory highlights the predictive capability of our approach. Critically, this work provides a comprehensive database of quantum-accurate solute segregation enthalpies in vanadium, enabling the development of advanced alloys for fusion reactors applications.

Fusion materials↗
DOE OSTI2025

Isopotential Electron Titration of Ammonia Charge Transfer on Metal Catalysts

Electron transfer between adsorbates and surfaces determines the binding strength and reactivity of chemical moieties on materials designed for separations and catalysis. To quantify electron exchange, the extent of charge transfer resulting from ammonia adsorption on a Ru surface was measured by isopotential electron titration (IET) on a Ru catalytic condenser, where isopotential conditions were maintained between Ru and silicon separated by an insulating HfO 2 layer during gas phase ammonia adsorption. Charge transfer upon ammonia adsorption on a Ru catalytic condenser increased from 40 to 1200 nC/cm 2 at 75 and 225 °C, respectively. Charge transfer measurements provided a direct estimate of ammonia adsorption thermodynamics on Ru without knowing surface coverages a priori, revealing an adsorption enthalpy of −53 ± 10 kJ/mol and entropy of −61 ± 26 J/mol·K. Combining experimentally-measured charge transfer with kinetic Monte Carlo simulations informed adsorbate surface coverages determined that 0.058 electrons were transferred to the Ru surface for each molecular ammonia adsorption event (δ NH 3 = 0.058 ± 0.005 e – /NH 3 *), consistent with calculated Bader charges. The ability to measure the extent of charge transfer for adsorbed species provides a fundamental descriptor to understand existing and new chemically functional surfaces, providing a foundational method for the emerging field of thermochemical surface coulometry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗