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 163 records · Page 9

TChem v2.0 - A Software Toolkit for the Analysis of Complex Kinetic Models

TChem is an open source software library for solving complex computational chemistry problems and analyzing detailed chemical kinetic models. The software provides support for: complex kinetic models for gas-phase and surface chemistry; thermodynamic properties based on NASA polynomials; species production/consumption rates; stable time integrator for solving stiff time ordinary differential equations; and, reactor models such as homogenous gas-phase ignition (with analytical Jacobian matrices), continuously stirred tank reactor, plug-flow reactor. This toolkit builds upon earlier versions that were written in C and featured tools for gas-phase chemistry only. The current version of the software was completely refactored in C++, uses an object-oriented programming model, and adopts Kokkos as its portability layer to make it ready for the next generation computing architectures i.e., multi/many core computing platforms with GPU accelerators. We have expanded the range of kinetic models to include surface chemistry and have added examples pertaining to Continuously Stirred Tank Reactors (CSTR) and Plug Flow Reactor (PFR) models to complement the homogenous ignition examples present in the earlier versions. To exploit the massive parallelism available from modern computing platforms, the current software interface is designed to evaluate samples in parallel, which enables large scale parametric studies, e.g. for sensitivity analysis and model calibration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stratospheric aerosol forcing for CMIP7 – Part 1: optical properties for pre-industrial, historical, and scenario simulations

Stratospheric aerosols, most of which originate from explosive volcanic sulfur emissions into the stratosphere, are a key natural driver of climate variability. They are thus a forcing provided by the Coupled Model Intercomparison Project (CMIP) Climate Forcings Task Team to climate modelling groups participating in phase 7 of CMIP. For the historical period, we provide two datasets covering 1750–2023: (i) a volcanic upper tropospheric-stratospheric sulfur emission dataset, documented in a companion paper; and (ii) a stratospheric sulfate aerosol optical property dataset, documented in this paper. For the satellite era (from 1979 onwards), stratospheric aerosol optical properties are derived from the Global Space-based Stratospheric Aerosol Climatology (GloSSAC) dataset. For the pre-satellite era (1750–1978), optical properties are derived from our volcanic SO 2 emission dataset using a new version of the reduced-complexity volcanic aerosol model Easy Volcanic Aerosol (Height) (EVA_H v2). A background, non-volcanic stratospheric aerosol climatology is derived from the 1998–2001 period with a trend over 1850–1978 accounting for increasing anthropogenic aerosols. A monthly stratospheric aerosol climatology is derived from the 1850–2021 average for both pre-industrial and scenario (future) simulations, with a 9-year ramp over 2022–2030 for scenario simulations to ensure a smooth transition from the historical period. CMIP7 simulations should use version 2.2.1 of our pre-industrial and historical datasets, and version 2.2.2 of our scenario dataset. Our methodology to produce historical aerosol optical properties significantly differs from CMIP6 for the pre-satellite era, and the resulting forcings in turn largely differ. In particular, the CMIP6 dataset was mostly based on the sparse and uncertain pyrheliometer record, which resulted in strongly underrepresented emissions from small-to-moderate magnitude eruptions. The resulting bias is addressed in CMIP7, which is entirely emission-derived in the pre-satellite era and uses more recent ice-core-based volcanic sulfur emission inventories than CMIP6. Our approach results in an overall larger volcanic aerosol forcing for CMIP7, with the 1850–2014 mean mid-visible global mean stratospheric aerosol optical depth (SAOD) in CMIP7 (0.0138) being 29 % higher than in CMIP6 (0.0107). The pre-industrial mean of the same variable is 26 % higher in CMIP7 (0.0135, derived from the historical 1850–2021) than CMIP6 (0.0107, derived from the historical 1850–2014 mean). Using a reduced-complexity climate model, we simulate a global mean surface temperature that is 0.07 °C colder for 1850–1900 when using the CMIP7 dataset instead of CMIP6, whereas 2000–2014 is 0.03 °C warmer in CMIP7. Our dataset also exhibits lower forcing for 1960–1980, resulting in temperatures 0.06 °C warmer when averaged over 1960–1990, a period for which CMIP6 climate models exhibit a cold bias. Given the large uncertainties characterizing the dataset, in particular for the pre-satellite era, we advise against treating the CMIP7 or CMIP6 dataset as uniquely superior for any specific year and highlight the need for further evaluation. We conclude the study by discussing sources of uncertainty for the dataset, future research avenues to improve it, as well as requirements to operationalize the production of the dataset, i.e. extend it and update it on an annual basis instead of every 5–7 years following CMIP cycles.

Aubry, Thomas J. [Univ. of Exeter, Devon (United K↗

Measurements of surface energy budgets in the rocky mountains of Colorado

The surface energy balance of a mountain valley and a mountain peak were continuously monitored for periods of 126 and 34 days, respectively. The data yielded temporally resolved, high-quality measurements of the radiative, soil, and total atmospheric turbulent heat fluxes. Comparative data for the two sites reveal the effects of terrain setting and surface cover in modulating the heat balance. A procedure is suggested for accommodating local advective affects due to variable terrain and surface cover to allow partitioning of the total turbulent flux into values for sensible and latent heat. This procedure relies on a stability-dependent coefficient of turbulent transfer which is locally calibrated to systematically minimize advection-related errors. Scaling functions for the transfer coeffient are based on the effective blackbody temperature of the surface. Consequently, the procedure offers an opportunity for developing reliable bulk representations of heat transfer at the surface for use with remote sensing data and with models of the soil and surface cover in complex terrain.

Sheaffer, John D.↗

Variable Complexity Structural Optimization of Shells

Structural designers today face both opportunities and challenges in a vast array of available analysis and optimization programs. Some programs such as NASTRAN, are very general, permitting the designer to model any structure, to any degree of accuracy, but often at a higher computational cost. Additionally, such general procedures often do not allow easy implementation of all constraints of interest to the designer. Other programs, based on algebraic expressions used by designers one generation ago, have limited applicability for general structures with modem materials. However, when applicable, they provide easy understanding of design decisions trade-off. Finally, designers can also use specialized programs suitable for designing efficiently a subset of structural problems. For example, PASCO and PANDA2 are panel design codes, which calculate response and estimate failure much more efficiently than general-purpose codes, but are narrowly applicable in terms of geometry and loading. Therefore, the problem of optimizing structures based on simultaneous use of several models and computer programs is a subject of considerable interest. The problem of using several levels of models in optimization has been dubbed variable complexity modeling. Work under NASA grant NAG1-1808 has been concerned with the development of variable complexity modeling strategies with special emphasis on response surface techniques. In addition several modeling issues for the design of shells of revolution were studied.

Haftka, Raphael T.↗

Variable Complexity Structural Optimization of Shells

Structural designers today face both opportunities and challenges in a vast array of available analysis and optimization programs. Some programs such as NASTRAN, are very general, permitting the designer to model any structure, to any degree of accuracy, but often at a higher computational cost. Additionally, such general procedures often do not allow easy implementation of all constraints of interest to the designer. Other programs, based on algebraic expressions used by designers one generation ago, have limited applicability for general structures with modem materials. However, when applicable, they provide easy understanding of design decisions trade-off. Finally, designers can also use specialized programs suitable for designing efficiently a subset of structural problems. For example, PASCO and PANDA2 are panel design codes, which calculate response and estimate failure much more efficiently than general-purpose codes, but are narrowly applicable in terms of geometry and loading. Therefore, the problem of optimizing structures based on simultaneous use of several models and computer programs is a subject of considerable interest. The problem of using several levels of models in optimization has been dubbed variable complexity modeling. Work under NASA grant NAG1-2110 has been concerned with the development of variable complexity modeling strategies with special emphasis on response surface techniques. In addition, several modeling issues for the design of shells of revolution were studied.

Haftka, Raphael T.↗

Correlations Between Sea-Surface Salinity Tendencies and Freshwater Fluxes in the Pacific Ocean

Temporal changes in sea-surface salinity (SSS) from 21 years of a high resolution model integration of the Pacific Ocean are correlated with the freshwater flux that was used to force the integration. The correlations are calculated on a 1 x10 grid, and on a monthly scale to assess the possibility of deducing evaporation minus precipitation (E-P) fields from the salinity measurements to be taken by the upcoming Aquarius/SAC-D mission. Correlations between the monthly mean E-P fields and monthly mean SSS temporal tendencies are mainly zonally-oriented, and are highest where the local precipitation is relatively high. Nonseasonal (deviations from the monthly mean) correlations are highest along mid-latitude storm tracks and are relatively small in the tropics. The response of the model's surface salinity to surface forcing is very complex, and retrievals of freshwater fluxes from SSS measurements alone will require consideration of other processes, including horizontal advection and vertical mixing, rather than a simple balance between the two.

Li, Zhen↗

Copper migration and surface oxidation of CuxBi2Se3 in ambient pressure environments

Abstract Chemical modifications such as intercalation can be used to modify surface properties or to further functionalize the surface states of topological insulators (TIs). Using ambient pressure x-ray photoelectron spectroscopy, we report copper migration in C u x B i 2 S e 3 , which occurs on a timescale of hours to days after initial surface cleaving. The increase in near-surface copper proceeds along with the oxidation of the sample surface and large changes in the selenium content. These complex changes are further modeled with core-level spectroscopy simulations, which suggest a composition gradient near the surface which develops with oxygen exposure. Our results shed light on a new phenomenon that must be considered for intercalated TIs—and intercalated materials in general—that surface chemical composition can change when specimens are exposed to ambient conditions.

36 MATERIALS SCIENCE↗

Cathodic Protection Modeling for Hanford Underground Double-Shell Tank Farms

Hanford stores millions of gallons of radioactive and chemically hazardous waste from the production of weapon materials in tank farms consisting of underground carbon-steel storage tanks surrounded by reinforced concrete. Six of these Hanford tank farms use double-shell storage tanks (DSTs). The DST farms were constructed from 1968 to 1986 with a planned 40–50 year design life, so some are already operating beyond their initial life expectancy. Ultrasonic testing (UT) has indicated significant thinning on the bottom of the secondary (outer) liner of these tanks, believed to arise from groundwater intrusion driving concrete side corrosion. There is no direct access to the steel/concrete interface between the tank and the concrete pad, making it difficult to apply a chemical-based mitigation strategy or to conduct repairs, but cathodic protection (CP) is a possible method to inhibit further concrete-side corrosion. Hanford already uses CP to protect below grade steel piping within the tank farms and connected to the tanks, but this system was not designed to protect the tank bottoms. CP design must account for the structures surrounding the DSTs, including the steel reinforcing bars (rebar) within the concrete pad and vault, various process lines, and the existing CP system. In this study, finite element analysis (FEA) modeling was carried out to simulate CP protection of 1) a single tank and CP anode to develop options for modeling the rebar and to compare to a simpler circuit model and 2) the entire Hanford AN tank farm as a representative example consisting of seven tanks, associated piping, and both existing and new CP anodes. Both circuit and FEA models predict that significant protective current could be delivered to the bottoms of the tanks with the addition of tank-protection anodes below the depth of the tanks. Simulations with only the existing pipe-protection anodes active confirmed that only a very small current to the tank bottoms is predicted under present conditions. Multiple simplified representations of the dome and wall rebar were tested to reduce the computational complexity of the tank-farm simulations, resulting in modeling the rebar as edge elements with a prescribed effective circumference that matches the real rebar surface area. The geometry of the rebar is also simplified into horizontal hoops around the tank walls and radial rebar over the dome with increased effective circumference to retain the target surface area. This simplification was found to greatly reduce the complexity and solution time of the models without large changes in current distributions, especially to the tank bottom. A range of values were tested for model parameters such as soil and concrete resistivities and polarization resistance to investigate their impact on the current and electric potential distributions. Depending on the parameters used, FEA simulations predict some risk of overprotection, particularly on the piping system; since overprotection can also lead to surface damage associated with hydrogen gas generation at the interface (e.g. hydrogen embrittlement or damage to coatings), this needs to be considered when refining the design of the new CP system. Comparison between the FEA models and the circuit model representation demonstrated that the circuit model could not match the predicted FEA current distribution, even when using the exact same surface areas. This discrepancy appeared to be at least partly attributable to the impact of the relative positions of the tank components and anodes to each other and to the ground surface. The FEA model accounts for the relative positions since it solves the governing equations in three dimensions, but the circuit model cannot account for the positioning. In particular, the circuit model underpredicts the current to the tank bottom and overpredicts the current to the dome compared to FEA for the baseline geometry. The FEA models omitted the electrically isolated rebar in the bottom concrete slab. However, a circuit based stray current model estimated that only 2.1% of the total current through the slab would stray into the rebar, corresponding to ~0.21 A for a target current density of 2 mA/ft2 to the tank bottom. The estimated corrosion driven by this amount of stray current is predicted to yield a lifetime of >400 years for the minimum rebar diameter, assuming an acceptable cross-section area loss of 10%.

d'Entremont, Anna [Savannah River National Laborat↗

The role of dissolved pyrogenic carbon from biochar in the sorption of As(V) in biogenic iron (oxyhydr)oxides

Water contamination by arsenic (As) affects millions of people around the world, making techniques to immobilize or remove this contaminant a pressing societal need. Biochar and iron (oxyhydr)oxides [in particular, biogenic iron (oxyhydr)oxides (BIOS)] offer the possibility of stabilizing As in remediation systems. However, little is known about the potential antagonism in As sorption generated by the dissolved organic carbon (DOC) from biochar, or whether DOC affects how As(V) interacts with BIOS. For this reason, our objectives were to evaluate the i) As(V) sorption potential in BIOS when there is presence of DOC from pyrolyzed biochars at different temperatures; and ii) identify whether the presence of DOC alters the surface complexes formed by As(V) sorbed in the BIOS. We conducted As(V) sorption experiments with BIOS at circumneutral pH conditions and in the presence of DOC from sugarcane (Saccharum officinarum) straw biochar at pyrolyzed 350 (BC350) and 750 °C (BC750). The As(V) content was quantified by inductively coupled plasma mass spectrometry, and the BIOS structure and As(V) sorption mechanisms were investigated by X-ray absorption spectroscopy. In addition, the organic moieties comprising the DOC from biochars were investigated by attenuated total reflectance Fourier transform infrared spectroscopy. The addition of DOC did not change the biomineral structure or As(V) oxidation state. The presence of DOC, however, reduced by 25 % the sorption of As(V), with BC350 being responsible for the greatest reduction in As(V) sorption capacity. Structural modeling revealed As(V) predominantly formed binuclear bidentate surface complexes on BIOS. The presence of DOC did not change the binding mechanism of As(V) in BIOS, suggesting that the reduction of As(V) sorption to BIOS was due to site blocking. In conclusion, our results bring insights into the fate of As(V) in surface waters and provide a basis for understanding the competitive sorption of As(V) in environments with biochar application.

54 ENVIRONMENTAL SCIENCES↗

Analysis of the performance of different culvert boundary conditions in 2D shallow flow models

Culverts allow roads to safely traverse small streams or drainage ditches, and their proper design is critical to ensure a safe and reliable transportation network. A correct modelization of these hydraulic structures becomes crucial in the assessment of flood footprints or discharge peak estimation in a risk evaluation plan. The question of how to include culverts comes up frequently when assembling a hydraulic model that requires the presence of as many singular elements as possible. In this work, three different culvert integrations with the surface domain are studied and compared in the context of a 2D shallow water (SW) model. All of them are based on the Federal Highway Administration (FHWA) formulation for the culvert discharge estimation but differ in complexity and in the interaction with the numerical model for surface flow, some of them as internal boundary conditions. Several steady and unsteady validation test cases are presented and the numerical results are compared with the predictions from HEC-RAS 1D and HY-8 software. Furthermore, the culvert area, shape and their sensitivity to the 2D computational mesh is also analyzed.

54 ENVIRONMENTAL SCIENCES↗

Modeling of surface mechanical behaviors of soft elastic solids: theory and examples

Surfaces of soft solids can have significant surface stress, extensional modulus and bending stiffness. Previous theoretical studies have usually examined cases in which both the surface stress and bending stiffness are constant, assuming small deformation. In this work we consider a general formulation in which the surface can support large deformation and carry both surface stresses and surface bending moments. We show that the large deformation theory can be reduced to the classical linear theory (Shuttleworth equation). We obtain exact solutions for problems of an inflated cylindrical shell and bending of a plate with a finite thickness. Our analysis illustrates the different manners in which surface stiffening and surface bending stabilize these structures. We discuss how the complex surface constitutive behaviors affect the stress field of the bulk. Our calculation provides insights into effects of strain-dependent surface stress and surface bending in the large deformation regime, and can be used as a model to implement surface finite elements to study large deformation of complex structures.

36 MATERIALS SCIENCE↗

Transforming ESM Physical Parameterization Development Using Machine Learning Trained on Global Cloud-Resolving Models and Process Observations

ESMs robustly predict that 21st century greenhouse warming will slowly increase global mean precipitation, rapidly increase extreme precipitation, and increase subtropical drought. ESMs agree less about precipitation trends and extremes over particular land regions critical to human societies, e. g. in semi-arid regions such as California or the Sahel, or in wetter climates prone to monsoonal rainfall (e. g. southeast Asia) or to tropical cyclones and flooding from mesoscale convective systems (e. g. the southeastern U.S.) Deep convective parameterizations and poor representation of orography and complex vegetated land surfaces contribute to this inter-model spread; clouds, aerosols and sea-surface temperature biases are also key. Reducing regional precipitation projection uncertainty has enormous planning value for water supplies, land use, wildfire, hydropower, flood control, etc. IPCC-class ESMs are making painfully slow progress on this.

54 ENVIRONMENTAL SCIENCES↗

Insights Into the Hydrology of the Congo Peatlands Through Land Surface Modeling and Data Assimilation

The 16.8 million ha of peatlands in the Cuvette Centrale wetland complex in the Congo Basin is one of the largest peatland regions on Earth but still highly understudied. Understanding the hydrological functioning of these peatlands and the effects of external disturbances thereon remains a major challenge. Recent research suggested fundamental hydrological differences between the Congo peatlands and the well-studied Southeast Asian peatlands. The Congo peatlands have a doming gradient that is up to ten times smaller, and they are influenced by river hydrology to some extent. In this study, we explore the Congo peatland hydrology through land surface modeling and data assimilation. We build upon our recently developed tropical PEATCLSM module (Apers et al., 2022) that was parameterized based on data from Southeast Asian peatlands due to the lack of field data from other tropical peatland regions. In a first step, we derive Congo-specific peat hydraulic and discharge function parameters from a scalar parametrization of water level dynamics in the Congo peatlands, using observed water level data at two locations. These Congo-specific parameters differ considerably from the original literature-based parameters from Southeast Asian peatlands. In a second step, we apply our original and Congo-specific parameters in an assimilation scheme for L-band brightness temperature (Tb) data from the Soil Moisture and Ocean Salinity (SMOS) mission. The data assimilation results are used in two ways. First, the effect of these parameters on the simulated peatland hydrology and the observation-minus-forecast Tb residuals is evaluated. It is hypothesized that the new parameters reduce the previously reported modeling errors over the Congo peatlands and reduce the residuals in Tb as well. Second, we analyze the data assimilation diagnostics to learn about other model improvement possibilities. In preliminary results, we found long periods of temporally autocorrelated total water storage increments (difference of pre- and post-update) that coincided with anomalies in river stages measured upstream of the peatlands. Since PEATLCSM neglects possible river influence, this concurrence suggests that the typically used grid-based approach of land surface models should be combined with a river routing scheme over the Congo peatlands.

Sebastian Apers↗

Experimental and modeling investigation of binary liquid mixtures

Binary mixtures of liquids may be encountered in industrial or remote sensing scenarios and present challenges to positive identification compared to neat single-component liquids. Our investigation examines whether one can predict the optical properties of the mixture, i.e. its complex index of refraction, by assuming a linear superposition of the real and imaginary components of the index of refraction in proportion to the ratio of each constituent. To investigate this hypothesis various liquid mixtures were created using mass ratios. The mixtures were then characterized as to their complex index of refraction and used in numerical modeling calculations of thin liquid mixture films on surfaces and compared with composed mixtures using linear n and k synthetic mixtures where the n and k components of the complex index of refraction were combined in similar ratios. The comparison of modeling and experimental results is presented with recommendations for further investigation.

infrared (IR) spectroscopy, Liquid mixtures, compl↗

On the interaction of surface heating anomalies with zonally symmetric and asymmetric atmospheric flows

Models of intermediate complexity have been used to study some aspects of the climatic effects of anomalous heating, but many aspects of the problem have yet to be explored thoroughly. The present study represents a preliminary investigated of the gaps in scientific understanding of the interaction of heating and atmospheric dynamics. The principle research tool is a model of intermediate complexity, including a time-dependent, nonlinear-two-layer quasi-geostrophic model of relatively high horizontal resolution which incorporates simple heating parameterizations. The model is used to examine systematically the interaction of heating arising from anomalies in surface temperature with zonally symmetric and zonally asymmetric flows characterized by different values of static stability and mean vertical wind shear.

Phillips, T. J.↗

Radiative force model performance for TOPEX/Poseidon precision orbit determination

The TOPEX/Poseidon spacecraft was launched to study the Earth's oceans. To maximize the benefit from the alimetric data collected, mission requirements dictate that TOPEX/Poseidon's orbit must be computed at extremely high accuracy. To meet these demands, a nonconservative force model which accounts for the satellite's complex geometry, attitude and surface properties has been developed. The 'box-wing' representation treats the spacecraft as the combination of flat plates arranged in the shape of a box and a connected solar array. Model performance and parameter sensitivities are discussed.

TRT-THEORETICAL↗

Heat Transfer from Radiatively Heated Material in a Low Reynolds Number Microgravity Environment

A mathematical model of the transient three-dimensional heat transfer between a slowly moving ambient gas stream and a thermally thick or thin flat surface heated by external radiation in a microgravity environment is presented. The problem is motivated in part by fire safety issues in spacecraft. The gas phase is represented by variable property convection-diffusion energy and mass conservation equations valid at low Reynolds numbers. The absence of gravity and low Reynolds number together permit the flow to be represented by a self-consistent velocity potential determined by the ambient velocity and the thermal expansion in the gas. The solid exchanges energy with the gas by conduction/convection and with the surroundings by surface absorption and re-emission of radiation. Heat conduction in the solid is assumed to be one dimensional at each point on the surface as a consequence of the limited times (of order of 10 seconds) of interest in these simulations. Despite the apparent simplicity of the model, the results show a complex thermally induced flow near the heated surface. The thermal exchange between the gas and solid produces an outward sourcelike flow upstream of the center of the irradiated area and a sinklike flow downstream. The responses of the temperature fields and the associated flows to changes in the intensity of the external radiation and the ambient velocity are discussed.

Yamashita, H.↗

Oxide– and Silicate–Water Interfaces and Their Roles in Technology and the Environment

Interfacial reactions drive all elemental cycling on Earth and play pivotal roles in human activities such as agriculture, water purification, energy production and storage, environmental contaminant remediation, and nuclear waste repository management. The onset of the 21st century marked the beginning of a more detailed understanding of mineral aqueous interfaces enabled by advances in techniques that use tunable high-flux focused ultrafast laser and X-ray sources to provide near-atomic measurement resolution, as well as by nanofabrication approaches that enable transmission electron microscopy in a liquid cell. This leap into atomic- and nanometer-scale measurements has uncovered scale-dependent phenomena whose reaction thermodynamics, kinetics, and pathways deviate from previous observations made on larger systems. A second key advance is new experimental evidence for what scientists hypothesized but could not test previously, namely, interfacial chemical reactions are frequently driven by “anomalies” or “non-idealities” such as defects, nanoconfinement, and other nontypical chemical structures. Third, progress in computational chemistry has yielded new insights that allow a move beyond simple schematics, leading to a molecular model of these complex interfaces. In combination with surface-sensitive measurements, we have gained knowledge of the interfacial structure and dynamics, including the underlying solid surface and the immediately adjacent water and aqueous ions, enabling a better definition of what constitutes the oxide– and silicate–water interfaces. This critical review discusses how science progresses from understanding ideal solid–water interfaces to more realistic systems, focusing on accomplishments in the last 20 years and identifying challenges and future opportunities for the community to address. Further, we anticipate that the next 20 years will focus on understanding and predicting dynamic transient and reactive structures over greater spatial and temporal ranges as well as systems of greater structural and chemical complexity. Closer collaborations of theoretical and experimental experts across disciplines will continue to be critical to achieving this great aspiration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗