Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation 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 1,369 records · Page 76

Exploring Causal Physical Mechanisms via Non-Gaussian Linear Models and Deep Kernel Learning: Applications for Ferroelectric Domain Structures

Rapid emergence of multimodal imaging in scanning probe, electron, and optical microscopies has brought forth the challenge of understanding the information contained in these complex data sets, targeting the intrinsic correlations between different channels, and further exploring the underpinning causal physical mechanisms. Here, we develop such an analysis framework for Piezoresponse Force Microscopy. We argue that under certain conditions, we can bootstrap experimental observations with the prior knowledge of materials structure to get information on certain nonobserved properties, and demonstrate linear causal analysis for PFM observables. We further demonstrate that the strength of individual causal links between complex descriptors can be ascertained using the deep kernel learning (DKL) model. In this DKL analysis, we use the prior information on domain structure within the image to predict the physical properties. This analysis demonstrates the correlative relationships between morphology, piezoresponse, elastic property, etc., at nanoscale. The prediction of morphology and other physical parameters illustrates a mutual interaction between surface condition and physical properties in ferroelectric materials. Overall, this analysis is universal and can be extended to explore the correlative relationships of other multichannel data sets, and allow for high-fidelity reconstruction of underpinning functionalities and physical mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermal performance of FRSI/graphite epoxy materials for the Orbiter TPS application

Radiant tests have been conducted to evaluate the performance of the Flexible Reusable Surface Insulation (FRSI) when bonded to Graphite/Epoxy (G/E) honeycomb sandwich structures. The results from this study indicate reasonable agreement with the baseline FRSI properties that existed prior to this test program. An updated Thermal Math Model (TMM) prediction technique was developed, with the use of the well-defined boundary conditions, which appears to be valid in analyzing the complex heat transfer associated with the honeycomb structure. The updated TM analysis of flight data resulted in good agreement of predicted temperatures with flight data for bondline responses of Orbiter FRSI/GE TPS applications.

Norman, I.↗

CFD validation for aerodynamic flows - Challenge for the '90's

The process of computational fluid dynamics (CFD) validation is described from two perspectives, numerical error validation and physical model validation. Errors associated with time and space discretization, grid refinement, numerical dissipation, and level of solution convergence, are all possible sources of numerical error. The second aspect of CFD code validation is associated with the governing equations and the physical models, e.g., chemistry and turbulence models. This type of validation requires comparisons with carefully conducted experiments and is the primary aspect of the validation process discussed in this paper. Examples of validation are shown for selected Reynolds-averaged Navier-Stokes codes for conditions from transonic through hypersonic. Since code applications are becoming more complex, it no longer suffices to use data from surface or integral quantities alone to provide the required validation. Flow field surveys and experimental boundary condition measurements are emerging as critical data that must be obtained for CFD code validation. Progress in the overall level of validation is shown to be improving. However, future challenges remain, some of which are enumerated.

Marvin, Joseph G.↗

Solution adaptive Cartesian grid methods for aerodynamic flows with complex geometries

Cartesian methods for computational fluid dynamics are considered to offer an accurate and robust approach for the simulation of aerodynamic flows around geometrically complex bodies. A part of the ongoing research in this domain is reviewed with the aim of providing insight into the fundamental challenges faced by the practitioners of this approach, and a guide to further research. The integration schemes used in Cartesian solvers are similar to those used in other approaches. The emphasis is on the geometric algorithms, surface modeling and boundary conditions required to design a successful Cartesian mesh scheme.

Aftosmis, Michael J.↗

Improving the Accuracy of Satellite Sea Surface Temperature Measurements by Explicitly Accounting for the Bulk-Skin Temperature Difference

The focus of this research was to determine whether the accuracy of satellite measurements of sea surface temperature (SST) could be improved by explicitly accounting for the complex temperature gradients at the surface of the ocean associated with the cool skin and diurnal warm layers. To achieve this goal, work was performed in two different major areas. The first centered on the development and deployment of low-cost infrared radiometers to enable the direct validation of satellite measurements of skin temperature. The second involved a modeling and data analysis effort whereby modeled near-surface temperature profiles were integrated into the retrieval of bulk SST estimates from existing satellite data. Under the first work area, two different seagoing infrared radiometers were designed and fabricated and the first of these was deployed on research ships during two major experiments. Analyses of these data contributed significantly to the Ph.D. thesis of one graduate student and these results are currently being converted into a journal publication. The results of the second portion of work demonstrated that, with presently available models and heat flux estimates, accuracy improvements in SST retrievals associated with better physical treatment of the near-surface layer were partially balanced by uncertainties in the models and extra required input data. While no significant accuracy improvement was observed in this experiment, the results are very encouraging for future applications where improved models and coincident environmental data will be available. These results are included in a manuscript undergoing final review with the Journal of Atmospheric and Oceanic Technology.

Wick, Gary A.↗

Analysis of Summer-Time Ozone and Precursor Species in the Southeast United States

Ozone (O3) is a greenhouse gas and toxic pollutant which plays a major role in air quality and atmospheric chemistry. The understanding and ability to model the horizontal and vertical structure of O3 mixing ratios is difficult due to the complex formation/destruction processes and transport pathways that cause large variability of O3. The Environmental Protection Agency has National Ambient Air Quality Standards for O3 set at 75 ppb with future standards proposed to be as low as 65 ppb. These lower values emphasize the need to better understand/simulate the transport processes, emission sources, and chemical processes controlling precursor species (e.g., NOx, VOCs, and CO) which influence O3 mixing ratios. The uncertainty of these controlling variables is particularly large in the southeast United States (US) which is a region impacted by multiple different emission sources of precursor species (anthropogenic and biogenic) and transport processes resulting in complex spatio-temporal O3 patterns. During this work we will evaluate O3 and precursor species in the southeast US applying models, ground-based and airborne in situ data, and lidar observations. In the summer of 2013, the UAH O3 Differential Absorption Lidar (DIAL) (part of the Tropospheric Ozone Lidar Network (TOLNet)) measured vertical O3 profiles from the surface up to approximately 12 km. During this period, the lidar observed numerous periods of dynamic temporal and vertical O3 structures. In order to determine the sources/processes impacting these O3 mixing ratios we will apply the CTM GEOS-Chem (v9-02) at a 0.25 deg x 0.3125 deg resolution. Using in situ ground-based (e.g., SEARCH Network, CASTNET), airborne (e.g., NOAA WP-3D - SENEX 2013, DC-8 - SEAC4RS), and TOLNet lidar data we will first evaluate the model to determine the capability of GEOS-Chem to simulate the spatio-temporal variability of O3 in the southeast US. Secondly, we will perform model sensitivity studies in order to quantify which emission sources (e.g., anthropogenic, biogenic, lighting, wildfire) and transport processes (e.g., stratospheric, long-range, local scale) are contributing to these TOLNet-observed dynamic O3 patterns. Results from the evaluation of the model and the study of sources/processes impacting observed O3 mixing ratios will be presented.

Optical Radar↗

The Interplay Between Radiation Pressure and the Photoelectric Instability in Optically Thin Disks of Gas and Dust

In optically thin disks, dust grains are photoelectrically stripped of electrons by starlight, heating nearby gas and possibly creating a dust clumping instability-the photoelectric instability (PeI)-that significantly alters global disk structure. In the current work, we use the Pencil Code to perform the first numerical models of the PeI that include stellar radiation pressure on dust grains in order to explore the parameter regime in which the instability operates. In some models with low gas and dust surface densities, we see a variety of dust structures, including sharp concentric rings. In the most gas- and dust-rich models, nonaxisymmetric clumps, arcs, and spiral arms emerge that represent dust surface density enhancements of factors of ∼5-20. In one high gas surface density model, we include a large, low-order gas viscosity and find that it observably smooths the structures that form in the gas and dust, suggesting that resolved images of a given disk may be useful for deriving constraints on the effective viscosity of its gas. Our models show that radiation pressure does not preclude the formation of complex structure from the PeI, but the qualitative manifestation of the PeI depends strongly on the parameters of the system. The PeI may provide an explanation for unusual disk morphologies, such as the moving blobs of the AU Mic disk, the asymmetric dust distribution of the 49Ceti disk, and the rings and arcs found in the HD 141569A disk.

photoelectric instability (PeI)↗

Benchmarking and Exploring Parameter Space of the 2-Phase Bubble Tracking Model for Liquid Mercury Target Simulation

High intensity proton pulses strike the Spallation Neutron Source (SNS)’s mercury target to provide bright neutron beams. These strikes deposit extensive energy into the mercury and its steel vessel. Prediction of the resultant loading on the target is difficult when helium gas is intentionally injected into the mercury to reduce the loading and to mitigate the pitting damage on the vessel. A 2-phase material model that incorporates the Rayleigh-Plesset (R-P) model is expected to address this complex multi-physics dynamics problem by including the bubble dynamics in the liquid mercury. We present a study comparing the measured target strains in the SNS target station with the simulation results of the solid mechanics simulation framework. We investigate a wide range of various physical model parameters, including the number of bubble families, bubble size distribution, viscosity, surface tension, etc. to understand their impact on simulation accuracy. Our initial findings reveal that using 8-10 bubble families in the model renders a simulation strain envelope that covers the experimental ones. Further optimization studies are planned to predict the strain response more accurately.

Lin, Lianshan↗

Modeling Surface-Enhanced Raman Scattering of Au-Pyrazine and Au-Pyrazine-Au Nanorod Dimer Systems with the TD-DFTB Method

This study investigates Raman scattering from a pyrazine molecule adsorbed on gold nanorod surfaces using the time-dependent density functional tight-binding (TD-DFTB) method. We analyze surface enhanced Raman scattering (SERS) of the molecule for two configurations: one where the molecule is adsorbed on a single gold nanorod forming Au-pyrazine complex, and another where it forms a nano-junction Au-pyrazine-Au dimer system. These two configurations offer distinct chemical environments for the molecule and different local field enhancements, with the dimer systems generating enhanced hotspot regions at the junctions. We present results for nanojunction, Au-pyrazine-Au dimer structures, for variable nanogap sizes that lead to shift in the plasmon energy. The study identifies contributions to the SERS enhancements from chemical, electromagnetic, and resonance charge transfer mechanisms. Our results show that the chemical mechanism provides enhancement factors in the range of 10 2 –10 3 , while the electromagnetic enhancements are on the order of 10 3 –10 5 for monomer systems and 10 6 –10 8 for dimer structures. We present results for a nanorod with a length of 5.3 nm and a width of 0.52 nm, consisting of 121 gold atoms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluating impacts of snow, surface water, soil and vegetation on empirical vegetation and snow indices for the Utqiaġvik tundra ecosystem in Alaska with the LVS3 model

Satellite observations for the Arctic and boreal region may contain information of vegetation, soil, snow, snowmelt, and/or other surface water bodies. We investigated the impacts of vegetation, soil, snow and surface water on empirical vegetation/snow indices on a tundra ecosystem area located around Utqiaġvik (formerly Barrow) of Alaska with the Moderate Resolution Imaging Spectrometer (MODIS) images in 2001–2014. Empirical vegetation indices, such as normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), the index of near infrared of vegetation (NIRv), and modified EVI (EVI2), have been used to monitor vegetation. Normalized difference snow index (NDSI) has been widely applied to monitor snow. The vegetation cover fraction (VGCF), the soil cover fraction (SOILCF), the snow cover fraction (SNOWCF), the surface water body cover fraction (WaterBodyCF), the fractional absorption of photosynthetically active radiation (PAR) by vegetation chlorophyll (fAPARchl), the fractional absorption of PAR by non-chlorophyll components of the vegetation (fAPARnon-chl), and the fractional absorption of PAR by the entire canopy (fAPARcanopy) are retrieved with the MODIS images and a coupled Leaf-Vegetation-Soil-Snow-Surface water body radiative transfer model, LVS3. The vegetation indices (NDVI, EVI, EVI2 and NIRv) differ from VGCF, fAPARchl, fAPARnon-chl, and fAPARcanopy. In addition to vegetation, we find that soil, snow and surface water also have impacts on vegetation indices NDVI, EVI (EVI2), and NIRv. Presence of snow makes lower the observed values of NDVI, EVI2 and NIRv. After snowmelt is gone, the vegetation indices (NDVI, EVI, EVI2 and NIRv) linearly decrease with SOILCF and WaterBodyCF, and WaterBodyCF has stronger impacts on these vegetation indices than SOILCF. The relationship between EVI and snow is complicated. NDSI non-linearly increases with SNOWCF, but linearly increases with sum of SNOWCF and WaterBodyCF (sum = 0.5893 × NDSI +0.4342, R2 = 0.976). NDSI linearly decreases with VGCF, and the relationship between NDSI and SOILCF is complex. Retrievals of VGCF, fAPARchl, fAPARnon-chl and fAPARcanopy with the LVS3 model provide alternatives for vegetation monitoring and ecological modeling.

Qingyuan Zhang↗

Thermal emission from particulate surfaces: A comparison of scattering models with measured spectra

Emissivity spectra of particulate mineral samples are highly dependent on particle size when that size is comparable to the wavelength of light emitted (5-50 micrometers for the midinfrared). Proper geologic interpretation of data from planetary infrared spectrometers will require that these particle size effects be well understood. To address this issue, samples of quartz powders were produced with narrow, well-characterized particle size distributions. Mean particle diameters in these samples ranged from 15 to 227 micrometers. Emission spectra of these powders allow the first detailed comparison of the complex spectral variations with particle size observed in laboratory data with the predictions of radiative transfer models. Four such models are considered here. Hapke's relectance theory (converted to emissivity via Kirchoff's law) is the first model tested. Hapke's more recently published emission theory is also employed. The third model, the 'Mie/Conel' model, uses Mie single scattering with a two-stream approximation for multiple scattering. This model, like the first, is a converted reflec- tance model. Mie scattering assumes particles are both spherical and well separated, which is not true for the quartz powders, but includes diffraction effects. The fourth model uses the Mie solution for single scattering by spheres and inputs those results into the multiple scattering formalism of Hapke's emission theory. The results of the four models are considered in relation to the values of the optical constants n and k. We have grouped these as class 1 (k large), class 2 (k moderate, n is approximately 2), class 3 (k small, n is approximately 2), and class 4 (k small, n is approximately 1). In general, the Mie/Hapke hybrid model does best at predicting variations with grain size. In particular, it predicts changes of the correct pattern, although incorrect magnitude, for class 1 bands, where large increases in emissivity with decreasing grain size are observed. This model also does an excellent job on moderate (class 2) and very weak and intraband (class 3) regions, and correctly predicts the emission maximum and its invariance with grain size near the Christiansen frequency (class 4). The Mie/Hapke hybrid model also has the fewest free parameters of the four models examined, while maintaining the most physical treatment of the radiative transfer. The Mie/Conel model performs as well as the Mie/Hapke hybrid model in strong bands (class 1) but does not accurately model the behavior of moderate (class 2) and very weak (class 3) bands.

Moersch, J. E.↗

TPSAS-NF1676L-33873-DND

NOAA's National Air Quality Forecasting Capability (NAQFC) model provides 48-hour forecast of surface ozone and other atmospheric constituents serving as the nation's key tool actively used by regulatory professionals in assessing air quality events and the issuance of public air quality alerts. NAQFC incorporates input meteorology from North American Mesoscale Model (NAM) at 12 km horizontal resolution that is coupled to the Community Multi-scale Air Quality Model (CMAQ) to predict emissions, tropospheric gas and aerosol chemistry, and related deposition and removal processes. The NAQFC forecast simulations contain 35 steps in altitude, with about one third of these in a typical afternoon daytime boundary layer (0-1.5 km). While several studies have been conducted comparing NAQFC predictions to surface observations, investigations comparing the vertical distribution to observations remain limited. Vertical comparisons are needed in order to fully assess forecast model behavior and atmospheric mixing conditions leading to the formation of ozone at the surface. The Tropospheric Ozone Lidar Network (TOLNet) instruments provide a unique observational capability in this regard, as ozone vertical measurements can be directly compared to models, such as NAQFC, and provide additional insight into atmospheric processes and enhance forecast capabilities. In this study, we utilize multiple TOLNet lidar observations obtained during the Ozone Water Land Environmental Transition Studies (OWLETS-2017 & OWLETS-2018) and the Long Island Sound Tropospheric Ozone Study (LISTOS-2018) to compare with NAQFC vertical predictions of ozone. These campaigns took place at different coastal locations along the US East coast with each campaign having multiple TOLNet lidars present, providing a unique opportunity to examine vertical structure of ozone in regions with complex water-land boundaries. Preliminary comparison results will be presented on ozone temporal and vertical differences between NAQFC and the TOLNet lidars along with example case studies to illustrate insights into boundary layer and surface processes of trace gases.

Timothy Berkoff↗

Mars Pathfinder Landing Site: Evidence for a Change in Wind Regime and Climate from Lander and Orbiter Data

Surface features related to the wind are observed in data from the Mars Pathfinder lander and from orbit by the Viking Orbiter and Mars Global Surveyor missions. Features seen from the surface include wind tails associated with small rocks, barchanoid duneforms, ripplelike patterns, and ventifact flutes cut into some rocks. Features seen from orbit include wind tails associated with impact craters, ridges inferred to be duneforms, and modified crater rims interpreted to have been eroded and mantled by windblown material. The orientations of these features show two prevailing directions, one inferred to represent winds from the northeast which is consistent with strongest winds predicted by a general circulation model to occur during the Martian northern winter under current conditions, and a second wind pattern oriented approx. 90 degrees to the first. This latter wind could be from the W-NW or from the E-SE and was responsible for cutting the ventifacts and modifying the crater rims. The two wind regimes could reflect a change in climate related to Mars' obliquity or some other, unknown factor. Regardless of the cause, the MPF area has been subjected to a complex pattern of winds and supply of small particles, in which the original surface formed by sedimentary processes from Tiu and Ares Vallis events has been modified by repeated burial and exhumation.

Greeley, R.↗

Peptoid-Loaded Microgels Self-Defensively Inhibit Staphylococcal Colonization of Titanium in a Model of Operating-Room Contamination

Bacterial contamination of an exposed implantable medical device by the atmosphere of an operating room (OR) is increasingly implicated as a cause of device-associated infection. Here, OR contamination is modeled in vitro using an aerosolizing system to spray small quantities of staphylococci onto titanium rods. Contaminated rods always manifest culturable bacteria. Self-assembly is used to create a self-defensive Ti surface that substantially enhances the rod's resistance to such contamination. Poly(acrylic acid) microgels are electrostatically deposited onto small Ti rods and subsequently loaded by complexation with a cationic antimicrobial peptoid (TM1). The microgels are visualized in situ by optical microscopy, and changes in microgel diameter indicate the loading state. Further, these measurements show that TM1 can be quickly loaded from low-ionic-strength buffer and subsequently remained sequestered within the microgels for up to 4 weeks when soaked in phosphate buffered saline. TM1-loaded microgel-modified Ti surfaces are contaminated with aerosolized staphylococci, and subsequent assays indicate few or no culturable bacteria. In the absence of nutrients to enable metabolism, this finding suggests that bacteria trigger local TM1 release by contact transfer. The modified surfaces exhibit good in vitro cytocompatibility as manifested by the adhesion, spreading, and metabolic activity of human fetal osteoblasts.

36 MATERIALS SCIENCE↗

Phosphorus recycling by mineral-catalyzed ribonucleotide cleavage on iron and manganese oxides

Abstract Phosphorus is an essential element influencing both food security via plant fertilization, and water pollution through excessive phosphorus use, yet the phosphorus cycle in ecosystems is poorly known. In particular, beyond adsorption, the role of iron and manganese oxides in catalyzing the abiotic dephosphorylation of biomolecules is debated. Here, we studied the reactions of ribonucleotides, containing different phosphate bonding, with goethite, hematite, and birnessite. We employed both high-resolution mass spectrometry of solution species and molecular modeling simulations of ribonucleotide-mineral complexes. Results disclose an up to fivefold preferential hydrolytic cleavage of a phosphoanhydride bond over a phosphoester bond, indicating that mineral-catalyzed reactions reflect the hierarchy reported for the activity of phosphatase enzymes. The fourfold higher catalytic reactivity of goethite and birnessite versus hematite is explained by mineral-specific binding rather than surface area differences. Corresponding simulated adsorbate conformations at the water–mineral interfaces are proposed. Overall, our findings provide new insights on the catalytic recycling of organic phosphorus species by mineral oxides.

58 GEOSCIENCES↗

The role of surface forces in environment-enhanced cracking of brittle solids

Fracture initiation and propagation in brittle materials is promoted in surface-reactive (sorptive) environments, a phenomenon known as subcritical crack growth (SCG). Laboratory measured crack-propagation velocity vs. stress intensity factor relationships typically exhibit highly nonlinear, multi-stage characteristics that are sensitive to environmental factors such as adsorbate concentration and temperature. For practical purposes, empirical relationships (e.g., a power law) have been used to describe this complex phenomenon. However, how the overall SCG behavior emerges from the underlying fundamental processes near the crack tip, such as the interaction of the crack surfaces separated by only a few nanometers and mass transport within the nano-confined space, is still not well understood. This paper develops a mechanistic, surface-force-based fracture theory (SFFT) which integrates surface force models, fluid transport models, and linear elastic fracture mechanics to quantitatively explain the multi-stage characteristics of SCG in brittle solids. A numerical model is developed based on SFFT and solved through an implicit partitioned scheme for efficiency and modularity. The results are validated by Wiederhorn's data on crack propagation in soda-lime glasses at a wide range of relative humidity levels. We show that, for the first time, the entire range of an SCG curve can be captured by a single physics-based model. The predicted SCG curves reveal that the development of repulsive disjoining pressure behind the crack tip can be responsible for the reduced apparent fracture toughness in a sorptive environment. In conclusion, the shape of the SCG curve, and its changes with respect to the environment, is found to critically depend on the assumed transport models.

36 MATERIALS SCIENCE↗

First-Principles Examination of Multiple Criteria of Organic Solvent Oxidative Stability in Batteries

Oxidative instability of the liquid electrolyte at or near battery cathode oxide surfaces has significant detrimental effects on batteries. Organic solvent molecules are often the fuel and precursors of such degradation processes, releasing electrons and protons that react with cathode oxides and electrolyte anions. These reactions contribute to cathode–electrolyte interphase (CEI) film formation, transition-metal ion dissolution, and phase transformation of the surface regions of the cathode. Here we apply density functional theory calculations to examine four criteria of oxidative stability (oxidation potential, hydrogen removal energies, and initial reactivity on two types of oxide facets) using four different solvent/additive molecules (ethylene carbonate, fluoroethylene carbonate, 1,3-dioxolane, and dimethyl ether). The ranking of molecular stability differs with each criterion. Surprisingly, the all-oxygen-terminated basal planes of layered oxides exhibit lower reaction barriers than spinel surface facets with exposed transition-metal cations, especially for ether solvents; the calculations also suggest basal planes contribute to the dissolution of transition-metal ions. Further, the structure–degradation relation complexity underscores the challenge of understanding the function of the CEI but also offers a guide to future degradation-mitigation strategies including facet engineering. Our predictions and models help establish a framework for future studies relevant to high-voltage conditions.

25 ENERGY STORAGE↗

Characterizing Surface Ice-Philicity Using Molecular Simulations and Enhanced Sampling

The formation of ice, which plays an important role in diverse contexts ranging from cryopreservation to atmospheric science, is often mediated by solid surfaces. Although surfaces that interact favorably with ice (relative to liquid water) can facilitate ice formation by lowering nucleation barriers, the molecular characteristics that confer ice-philicity to a surface are complex and incompletely understood. To address this challenge, here we introduce a robust and computationally efficient method for characterizing surface ice-philicity that combines molecular simulations and enhanced sampling techniques to quantify the free energetic cost of increasing surface–ice contact at the expense of surface–water contact. Using this method to characterize the ice-philicity of a family of model surfaces that are lattice matched with ice but vary in their polarity, we find that the nonpolar surfaces are moderately ice-phobic, whereas the polar surfaces are highly ice-philic. In contrast, for surfaces that display no complementarity to the ice lattice, we find that ice-philicity is independent of surface polarity and that both nonpolar and polar surfaces are moderately ice-phobic. Furthermore, our work thus provides a prescription for quantitatively characterizing surface ice-philicity and sheds light on how ice-philicity is influenced by lattice matching and polarity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗