Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Vertical Motion Simulator”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Trajectory Simulation for Air Traffic Management Employing a Multirotor Urban Air Mobility Aircraft Model

This paper describes a point-mass model of a multirotor electric vertical takeoff and landing aircraft designed to carry up to five passengers and a pilot. The model employs thrust, thrust vector angle and bank angle as controls, which are computed in part using the acceleration, heading-angle rate, and flight-path-angle rate generated by control systems in response to speed, heading angle, and flight-path angle commands. The equations of motion are integrated forward in time using the controls needed for following the desired vertical climb/descent and speed profiles, and horizontal path and cruise speed in the presence of wind to generate the flight trajectory.

trajectory simulation, urban air mobility aircraft↗
NASA NTRS1975

A study of moving base simulation motion cues utilizing washout technique

The present study was conducted to derive data on non-visual motion thresholds utilizing washout technique, and to develop specific threshold values for use as washout parameters. It describes the results of acceleration detection studies carried out using the NASA MSFC General Purpose Simulator which provides six degree-of-freedom cab motion. A series of experiments was performed to test the hypothesis that constant velocity visual cues might suppress non-visual deceleration sensitivity. The psychophysical method employed was the forced-choice procedure which theoretically yields a pure sensitivity measure. The MSFC General Purpose Simulator was programmed to provide velocity ramps in three axes - pitch, fore-aft translation, and vertical translation. Data were collected with and without a constant velocity visual input. Comparison of the current results with classical data suggested significantly greater sensitivity of the human observer to fore-aft accelerations than has generally been reported, sensitivity tending to depend on motion direction and sign of velocity change.

Kirkpatrick, M., III↗
NASA NTRS1987

Three-dimensional analysis of liquid oxygen sloshing in Space Shuttle external tank

A three-dimensional numerical simulation has been performed in order to investigate the hydrodynamics within the liquid oxygen tank of the Space Shuttle external tank after liftoff. The results indicate that the LOX surface undergoes very high vertical accelerations (up to 5 g) and that splashing almost certainly occurs. Although the number of slosh baffles is not found to affect the surface, it is noted that the number of baffles may affect the structural motion of the tank. Surface accelerations are similar to those observed in previous two-dimensional simulations.

Kannapel, M. D.↗
DOE OSTI2024

Investigating the interactions between wakes and floating wind turbines using FAST.Farm

As floating offshore wind progresses to commercial maturity, wake and array effects across a farm of floating offshore wind turbines (FOWTs) will become increasingly important. While wakes of land-based and bottom-fixed offshore wind turbines have been extensively studied, only recently has this topic become relevant for floating turbines. This work presents an investigation of the mutual interaction between the motions of floating wind turbines and wakes using FAST.Farm. While FAST.Farm has been extensively validated across a wide range of conditions, it has never been validated for FOWT applications. Hence, in the first part of this work, we validate FAST.Farm by comparing simulations of a single FOWT against high-fidelity results from large-eddy simulations available in the literature. The validation is based on wake meandering, mean wake deflection, and velocity deficit at different downstream locations. This validation showed that the original axisymmetric (polar) wake model of FAST.Farm overpredicts the vertical wake deflection induced by shaft tilt and floater pitch, while the new curled wake model is capable of properly capturing the vertical wake deflection. In the second part, we use FAST.Farm to analyze a small three-unit array of FOWTs with a spacing of 7 diameters across a wide range of environmental conditions. The same National Renewable Energy Laboratory 5 MW reference wind turbine atop the OC4-DeepCwind semisubmersible is adopted for the three FOWTs and for the validation against high-fidelity simulations. To assess the effect of the floating substructure, we compare the power production, tower-base moments, and blade-root moments obtained for the floating turbines with the results obtained in a fixed-bottom configuration. The main differences introduced by the floating substructure are the motions induced by the waves, the change in the natural frequencies of the tower caused by differences in the boundary condition at its base, and the larger vertical deflection of the wake deficit due to the mean pitch of the platform. The impact of these differences, as well as other minor effects, are analyzed in detail.

16 TIDAL AND WAVE POWER↗
NASA NTRS2011

Long-Term Evolution of the Aerosol Debris Cloud Produced by the 2009 Impact on Jupiter

We present a study of the long-term evolution of the cloud of aerosols produced in the atmosphere of Jupiter by the impact of an object on 19 July 2009. The work is based on images obtained during 5 months from the impact to 31 December 2009 taken in visible continuum wavelengths and from 20 July 2009 to 28 May 2010 taken in near-infrared deep hydrogen-methane absorption bands at 2.1-2.3 micron. The impact cloud expanded zonally from approximately 5000 km (July 19) to 225,000 km (29 October, about 180 deg in longitude), remaining meridionally localized within a latitude band from 53.5 deg S to 61.5 deg S planetographic latitude. During the first two months after its formation the site showed heterogeneous structure with 500-1000 km sized embedded spots. Later the reflectivity of the debris field became more homogeneous due to clump mergers. The cloud was mainly dispersed in longitude by the dominant zonal winds and their meridional shear, during the initial stages, localized motions may have been induced by thermal perturbation caused by the impact's energy deposition. The tracking of individual spots within the impact cloud shows that the westward jet at 56.5 deg S latitude increases its eastward velocity with altitude above the tropopause by 5- 10 m/s. The corresponding vertical wind shear is low, about 1 m/s per scale height in agreement with previous thermal wind estimations. We found evidence for discrete localized meridional motions with speeds of 1-2 m/s. Two numerical models are used to simulate the observed cloud dispersion. One is a pure advection of the aerosols by the winds and their shears. The other uses the EPIC code, a nonlinear calculation of the evolution of the potential vorticity field generated by a heat pulse that simulates the impact. Both models reproduce the observed global structure of the cloud and the dominant zonal dispersion of the aerosols, but not the details of the cloud morphology. The reflectivity of the impact cloud decreased exponentially with a characteristic timescale of 15 days; we can explain this behavior with a radiative transfer model of the cloud optical depth coupled to an advection model of the cloud dispersion by the wind shears. The expected sedimentation time in the stratosphere (altitude levels 5-100 mbar) for the small aerosol particles forming the cloud is 45-200 days, thus aerosols were removed vertically over the long term following their zonal dispersion. No evidence of the cloud was detected 10 months after the impact.

Sanchez-Lavega, A.↗
DOE OSTI2015

Turbulent kinetics of a large wind farm and their impact in the neutral boundary layer

High-resolution large-eddy simulation of the flow over a large wind farm (64 wind turbines) is performed using the HIGRAD/FIRETEC-WindBlade model, which is a high-performance computing wind turbine–atmosphere interaction model that uses the Lagrangian actuator line method to represent rotating turbine blades. These high-resolution large-eddy simulation results are used to parameterize the thrust and power coefficients that contain information about turbine interference effects within the wind farm. Those coefficients are then incorporated into the WRF (Weather Research and Forecasting) model in order to evaluate interference effects in larger-scale models. In the high-resolution WindBlade wind farm simulation, insufficient distance between turbines creates the interference between turbines, including significant vertical variations in momentum and turbulent intensity. The characteristics of the wake are further investigated by analyzing the distribution of the vorticity and turbulent intensity. Quadrant analysis in the turbine and post-turbine areas reveals that the ejection motion induced by the presence of the wind turbines is dominant compared to that in the other quadrants, indicating that the sweep motion is increased at the location where strong wake recovery occurs. Regional-scale WRF simulations reveal that although the turbulent mixing induced by the wind farm is partly diffused to the upper region, there is no significant change in the boundary layer depth. The velocity deficit does not appear to be very sensitive to the local distribution of turbine coefficients. However, differences of about 5% on parameterized turbulent kinetic energy were found depending on the turbine coefficient distribution. Furthermore, turbine coefficients that consider interference in the wind farm should be used in wind farm parameterization for larger-scale models to better describe sub-grid scale turbulent processes.

17 WIND ENERGY↗
DOE OSTI2017

Frontogenesis and frontal arrest of a dense filament in the oceanic surface boundary layer

The evolution of upper ocean currents involves a set of complex, poorly understood interactions between submesoscale turbulence (e.g. density fronts and filaments and coherent vortices) and smaller-scale boundary-layer turbulence. Here we simulate the lifecycle of a cold (dense) filament undergoing frontogenesis in the presence of turbulence generated by surface stress and/or buoyancy loss. This phenomenon is examined in large-eddy simulations with resolved turbulent motions in large horizontal domains using ${\sim}10^{10}$ grid points. Steady winds are oriented in directions perpendicular or parallel to the filament axis. Due to turbulent vertical momentum mixing, cold filaments generate a potent two-celled secondary circulation in the cross-filament plane that is frontogenetic, sharpens the cross-filament buoyancy and horizontal velocity gradients and blocks Ekman buoyancy flux across the cold filament core towards the warm filament edge. Within less than a day, the frontogenesis is arrested at a small width, ${\approx}100~\text{m}$ , primarily by an enhancement of the turbulence through a small submesoscale, horizontal shear instability of the sharpened filament, followed by a subsequent slow decay of the filament by further turbulent mixing. The boundary-layer turbulence is inhomogeneous and non-stationary in relation to the evolving submesoscale currents and density stratification. The occurrence of frontogenesis and arrest are qualitatively similar with varying stress direction or with convective cooling, but the detailed evolution and flow structure differ among the cases. Thus submesoscale filament frontogenesis caused by boundary-layer turbulence, frontal arrest by frontal instability and frontal decay by forward energy cascade, and turbulent mixing are generic processes in the upper ocean.

Mechanics↗
DOE OSTI2024

Numerical validation of scaling laws for stratified turbulence

Recent theoretical progress using multiscale asymptotic analysis has revealed various possible regimes of stratified turbulence. Notably, buoyancy transport can either be dominated by advection or diffusion, depending on the effective Péclet number of the flow. Two types of asymptotic models have been proposed, which yield measurably different predictions for the characteristic vertical velocity and length scale of the turbulent eddies in both diffusive and non-diffusive regimes. The first, termed a ‘single-scale model’, is designed to describe flow structures having large horizontal and small vertical scales, while the second, termed a ‘multiscale model’, additionally incorporates flow features with small horizontal scales, and reduces to the single-scale model in their absence. By comparing predicted vertical velocity scaling laws with direct numerical simulation data, we show that the multiscale model correctly captures the properties of strongly stratified turbulence within regions dominated by small-scale isotropic motions, whose volume fraction decreases as the stratification increases. Meanwhile its single-scale reduction accurately describes the more orderly, layer-like, quiescent flow outside those regions.

Mechanics↗
DOE OSTI2015

Epicyclic helical channels for parametric resonance ionization cooling

Proposed next-generation muon colliders will require major technical advances to achieve rapid muon beam cooling requirements. Parametric-resonance Ionization Cooling (PIC) is proposed as the final 6D cooling stage of a high-luminosity muon collider. In PIC, a half-integer parametric resonance causes strong focusing of a muon beam at appropriately placed energy absorbers while ionization cooling limits the beam’s angular spread. Combining muon ionization cooling with parametric resonant dynamics in this way should then allow much smaller final transverse muon beam sizes than conventional ionization cooling alone. One of the PIC challenges is compensation of beam aberrations over a sufficiently wide parameter range while maintaining the dynamical stability with correlated behavior of the horizontal and vertical betatron motion and dispersion. We explore use of a coupling resonance to reduce the dimensionality of the problem and to shift the dynamics away from non-linear resonances. PIC simulations are presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Bayesian Cloud Property Retrievals from ARM Active and Passive Measurements

The optimum use of the continuous measurements of thermodynamics, radiation, aerosols, clouds and precipitation from the DOE Atmospheric Radiation Measurement (ARM) program is key to achieve the DOE Atmospheric System Research (ASR)’s objectives. One of the key mission requirements is to retrieve cloud and precipitation properties, as well as vertical motion parameters, along the vertical cross- section defined by the profiling active sensors. Such retrievals are challenging to perform continuously in the entire spectrum of cloud and precipitation conditions due to the large natural microphysical and dynamical variability, the often-limited information content in the measurements, and the lack of proper characterization of measurement quality and uncertainty. Today, the acquisition of new remote and in-situ sensors by the ARM program creates opportunities to address the microphysical retrieval problem by exploiting new, more robust retrieval techniques and integrating various scattered advancements in both sensor techniques and retrieval algorithms. During this project, we constructed a robust Bayesian Markov chain Monte Carlo (MCMC) cloud property retrieval algorithm that includes a state of the art radar forward model. Our MCMC-based retrieval produces both the best estimate of height-resolved cloud and precipitation properties in the radar profile, as well as an estimate of the in-cloud vertical motion and turbulence. In addition, the MCMC algorithm automatically produces robust and flexible estimates of retrieval uncertainty. We tested the algorithm on several synthetic cloud profiles obtained from large eddy simulation (LES) models with bin-resolved microphysics.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2010

Granular Flow and Dynamics of Lunar Simulants in Excavating Implements

The exploration of the lunar surface will rely on properly designed excavation equipment for surface preparations and for collection of lunar regolith in In-Situ Resource Utilization (ISRU) processes. Performance efficiency, i.e minimizing loading forces while maximizing material collection, and mass and volume reductions are major design goals. The NASA Glenn Research Center has embarked on an experimental program to determine the flow characteristics and dynamic forces produced by excavation operations using various excavator bucket designs. A new large scale soil bin facility, 2.27 m x 5.94 m x 0.76 m (nominally 8 ft. x 20 ft. x 27 in.) in size, capable of accommodating moderately large test implements was used for the simulations of lunar operations. The soil bin is filled with GRC-3simulant (a mixture of industrial sands and silt with a particle size distribution and the bulk mechanical (shear) strength representative of an average of lunar regolith from different regions) and uses motorized horizontal rails and a vertical actuator to drive the implement through the lunar simulant soil. A six-axis load cell and encoders provide well resolved measurements of the three dimensional forces and torques and motion of the bucket. In addition, simultaneous video allows for the analysis of the flow behavior and structure formation of the regolith during excavation. The data may be useful in anchoring soil mechanic models and to provide engineering data for design consideration.

Agui, Juan H.↗
DOE OSTI2012

TEMPORAL EVOLUTION OF VELOCITY AND MAGNETIC FIELD IN AND AROUND UMBRAL DOTS

We study the temporal evolution of umbral dots (UDs) using measurements from the CRISP imaging spectropolarimeter at the Swedish 1 m Solar Telescope. Scans of the magnetically sensitive 630 nm iron lines were performed under stable atmospheric conditions for 71 minutes with a cadence of 63 s. These observations allow us to investigate the magnetic field and velocity in and around UDs at a resolution approaching 0.''13. From the analysis of 339 UDs, we draw the following conclusions: (1) UDs show clear hints of upflows, as predicted by magnetohydrodynamic simulations. By contrast, we could not find systematic downflow signals. Only in very deep layers, we detect localized downflows around UDs, but they do not persist in time. (2) We confirm that UDs exhibit weaker and more inclined fields than their surroundings, as reported previously. However, UDs that have strong fields above 2000 G or are in the decay phase show enhanced and more vertical fields. (3) There are enhanced fields at the migration front of UDs detached from penumbral grains, as if their motion were impeded by the ambient field. (4) Long-lived UDs travel longer distances with slower proper motions. Our results appear to confirm some aspects of recent numerical simulations of magnetoconvection in the umbra (e.g., the existence of upflows in UDs), but not others (e.g., the systematic weakening of the magnetic field at the position of UDs).

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRSDate not supplied

Developing Urban Air Mobility Vehicle Models to Support Air Traffic Management Concept Development

To support Urban Air Mobility (UAM) research efforts at NASA, the Airspace Target Generator (ATG) software used in the FutureFlight Central (FFC) air traffic control tower simulator is undergoing updates to support physics-based UAM vehicle models. A process was developed to integrate UAM vertical takeoff and landing (VTOL) aircraft into the fixed-wing ATG modeling environment without significant change to the underlying equations of motion and vehicle model database. The VTOL aircraft models were converted from a six degrees-of-freedom (6-DOF) representation into a four degrees-of-freedom (4-DOF) representation for integration within ATG. Three vehicle designs from the NASA Revolutionary Vertical-Lift Technologies (RVLT) project were selected: a lift-plus-cruise (LPC) aircraft model and quadrotor, electric-powered (QEP) 1-seater and 6-seater models. With the LPC model comprised of a nonlinear force and moment build-up, and the QEP models comprised of linearized stability derivatives, two separate processes were developed to convert the lift, drag, and propulsion characteristics of each model into the ATG model database. Key aircraft performance characteristics including climb, cruise, and descent performance were preserved during the conversion process. Because ATG simulates fixed-wing aircraft through ground taxi and takeoff to approach and landing, acceleration command algorithms were developed to model the vertical takeoff and vertical landing phase of UAM operations. A strategy was then developed to transition the aircraft model to- and from- the new control mode.

urban air mobility↗
DOE OSTI2015

Coupling between flexural modes in free vibration of single-walled carbon nanotubes

The nonlinear thermal vibration behavior of a single-walled carbon nanotube (SWCNT) is investigated by molecular dynamics simulation and a nonlinear, nonplanar beam model. Whirling motion with energy transfer between flexural motions is found in the free vibration of the SWCNT excited by the thermal motion of atoms where the geometric nonlinearity is significant. A nonlinear, nonplanar beam model considering the coupling in two vertical vibrational directions is presented to explain the whirling motion of the SWCNT. Energy in different vibrational modes is not equal even over a time scale of tens of nanoseconds, which is much larger than the period of fundamental natural vibration of the SWCNT at equilibrium state. The energy of different modes becomes equal when the time scale increases to the microsecond range.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Evaluating Large-Eddy Simulations of Boundary-Layer Cloud Response to Climate Change using MAGIC Observations

Over the cool waters off California, a persistent and extensive deck of low-lying clouds forms atop a shallow, turbulent layer of moist air extending 500-2000 meters above the ocean, above which the air is much warmer and drier. Its west edge varies with weather and season and is typically half way to Hawaii. This cloud deck and similar cloud decks elsewhere reflect sunlight back to space, cooling the underlying ocean and indeed the entire planet. They must be accurately represented by a climate model to correctly simulate our current climate or future global warming. Most climate models suggest these cloud decks will thin and break up in a warmer CO 2 -enriched climate, so extra sunlight would reach the ocean surface and amplify warming worldwide. But these models have difficulty accurately simulating the cloud thickness and extent because the turbulent cloud-forming air currents are much too small to be simulated by their coarse computational grid. In 2012-2013, the DOE Atmospheric Radiation Measurement MAGIC program instrumented a container ship transiting between California and Hawaii to observe clouds and rain, small particles (aerosols), and weather conditions. Data were obtained from 20 4-6 day cruises in each direction. This study tested whether a sophisticated ‘large-eddy simulation’ (LES) of the turbulent air motions that create the clouds in this region accurately predicts the variations in cloud properties seen across the different cruises. If so, that model is more believable for predicting cloud changes in this region in a future warmer climate. Our LES used a numerical mesh over a 6.4x6.4 km region with 50 m horizontal and 5 m vertical grid spacing to simulate the cloud-forming eddies and the sharp overlying transition between moist and dry air. To better compare with MAGIC observations, we configured our model to simulate a region that moved along with the ship. We used a weather forecast model to specify effects of the large-scale air flow. We compared our results with balloon ascents, a radar and other instruments pointed upward at the clouds, measurements of downward sunlight and infrared radiation, aerosols, and weather variables for all suitable cruise periods. Our LES was remarkably skillful in predicting the cloudiness variations measured by the ship. MAGIC also sampled mesoscale cells, regions of cloud thickening and thinning over distances of 10-50 km. An LES run over a larger region 100 km on a side shows these cells form spontaneously, enhanced by strong radiative heat loss at the cloud tops.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2010

Visual Performance Challenges to Low-Frequency Perturbations After Long-Duration Space Flight, and Countermeasure Development

Astronauts experience sensorimotor disturbances after long-duration space flight. After a water landing, crewmembers may need to egress the vehicle within a few minutes for safety and operational reasons in various sea state conditions. Exposure to even low-frequency motions induced by sea conditions surrounding a vessel can cause significant motor control problems affecting critical functions. The first objective of this study was to document human visual performance during simulated wave motion below 2.0 Hz. We examined the changes in accuracy and reaction time when subjects performed a visual target acquisition task in which the location of the target was offset vertically during horizontal rotation at an oscillating frequency of 0.8 Hz. The main finding was that both accuracy and reaction time varied as a function of target location, with greater performance decrements occurring when vertical targets were acquired at perturbing frequencies of 0.8 Hz in the horizontal plane. A second objective was to develop a countermeasure, base d on stochastic resonance (SR), to enhance sensorimotor capabilities with the aim of facilitating rapid adaptation to gravitational transitions after long-duration space flight. SR is a mechanism by which noise can enhance the response of neural systems to relevant sensory signals. Recent studies have shown that applying imperceptible stochastic electrical stimulation to the vestibular system (SVS) significantly improved balance and oculomotor responses. This study examined the effectiveness of SVS on improving balance performance. Subjects performed a standard balance task while bipolar SVS was applied to the vestibular system using constant current stimulation through electrodes placed over the mastoid process. The main finding of this study was that balance performance with the application of SR showed significant improvement in the range of 10%-25%. Ultimately an SR-based countermeasure might be fielded either as preflight training to enhance adaptability, or as a miniature patch-type stimulator worn post flight.

Mulavara, Ajitkumar↗
DOE OSTI2016

Rise of an argon bubble in liquid steel in the presence of a transverse magnetic field

The rise of gaseous bubbles in viscous liquids is a fundamental problem in fluid physics, and it is also a common phenomenon in many industrial applications such as materials processing, food processing, and fusion reactor cooling. In this work, the motion of a single argon gas bubble rising in quiescent liquid steel under an external magnetic field is studied numerically using a Volume-of-Fluid method. To mitigate spurious velocities normally generated during numerical simulation of multiphase flows with large density differences, an improved algorithm for surface tension modeling, originally proposed by Wang and Tong [“Deformation and oscillations of a single gas bubble rising in a narrow vertical tube,” Int. J. Therm. Sci. 47, 221–228 (2008)] is implemented, validated and used in the present computations. The governing equations are integrated by a second-order space and time accurate numerical scheme, and implemented on multiple Graphics Processing Units with high parallel efficiency. The motion and terminal velocities of the rising bubble under different magnetic fields are compared and a reduction in rise velocity is seen in cases with the magnetic field applied. The shape deformation and the path of the bubble are discussed. An elongation of the bubble along the field direction is seen, and the physics behind these phenomena is discussed. The wake structures behind the bubble are visualized and effects of the magnetic field on the wake structures are presented. A modified drag coefficient is obtained to include the additional resistance force caused by adding a transverse magnetic field.

42 ENGINEERING↗
DOE OSTI2018

Electronic spectrum and characterization of diabatic potential energy surfaces for thiophenol

The electronic spectrum of thiophenol was simulated by a normal-mode sampling approach combined with TDDFT in the Tamm–Dancoff approximation (TDA). The vertical excitation energies were compared with electronic structure calculations by completely renormalized equation-of-motion coupled cluster theory with single and double excitations and noniterative inclusion of connected triples (CR-EOM-CCSD(T)) and by multi-reference perturbation theory. The spectrum was computed both with and without solvation effects, and these spectra are compared to each other and to experiment. Using multireference-perturbation-theory adiabatic wave functions and model-space diabatization by the fourfold way, diabatic potential energy surfaces of the lowest three singlet states ( 1 ππ, 1 ππ * , and 1 n π σ * ) were constructed along the S–H stretching coordinate, the C–C–S–H torsion coordinate, and the v 16a and v 16b normal coordinates. The first two of these two are primary coordinates for the photodissociation, and the diabatic crossing seams of the three states were calculated and plotted as functions of the two coordinates. The other two coordinates are out-of-plane ring distortion modes studied to assess the extent of their role in coupling the states near the first conical intersection, and the v 16a mode was shown to be an important coupling mode there. The current study is the first step toward a detailed mechanistic analysis of the photoinduced S–H fission process of thiophenol, a test system to understand 1 n π σ * -mediated reactions, at the same time already providing a better understanding of the thiophenol electronic excitations by clarifying the assignment of the experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗