Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Surface estimation”

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 271 records · Page 15

Physical Controls on the Hydrology of Perennially Ice-Covered Lakes, Taylor Valley, Antarctica (1996–2013)

The McMurdo Dry Valleys, Antarctica, are a polar desert populated with numerous closed-watershed, perennially ice-covered lakes primarily fed by glacial melt. Lake levels have varied by as much as 8 m since 1972 and are currently rising after a decade of decreasing. Precipitation falls as snow, so lake hydrology is dominated by energy available to melt glacier ice and to sublimate lake ice. To understand the energy and hydrologic controls on lake level changes and to explain the variability between neighboring lakes, only a few kilometers apart, we model the hydrology for the three largest lakes in Taylor Valley. Here, we apply a physically based hydrological model that includes a surface energy balance model to estimate glacial melt and lake sublimation to constrain mass fluxes to and from the lakes. Results show that lake levels are very sensitive to small changes in glacier albedo, air temperature, and wind speed. We were able to balance the hydrologic budget in two watersheds using meltwater inflow and sublimation loss from the ice-covered lake alone. A third watershed, closest to the coast, required additional inflow beyond model uncertainties. We hypothesize a shallow groundwater system within the active layer, fed by dispersed snow patches, contributes 23% of the inflow to this watershed. The lakes are out of equilibrium with the current climate. If the climate of our study period (1996–2013) persists into the future, the lakes will reach equilibrium starting in 2300, with levels 2–17 m higher, depending on the lake, relative to the 2020 level.

58 GEOSCIENCES↗

Wildfire aerosol deposition likely amplified a summertime Arctic phytoplankton bloom

Summertime wildfire activity is increasing in boreal forest and tundra ecosystems in the Northern Hemisphere. However, the impact of long range transport and deposition of wildfire aerosols on biogeochemical cycles in the Arctic Ocean is unknown. Here, we use satellite-based ocean color data, atmospheric modeling and back trajectory analysis to investigate the transport and fate of aerosols emitted from Siberian wildfires in summer 2014 and their potential impact on phytoplankton dynamics in the Arctic Ocean. We detect large phytoplankton blooms near the North Pole (up to 82°N in the eastern Eurasian Basin). Our analysis indicates that these blooms were induced by the northward plume transport and deposition of nutrient-bearing wildfire aerosols. We estimate that these highly stratified surface waters received large amounts of wildfire-derived nitrogen, which alleviated nutrient stress in the phytoplankton community and triggered an unusually large bloom event. Our findings suggest that changes in wildfire activity may strongly influence summertime productivity in the Arctic Ocean.

54 ENVIRONMENTAL SCIENCES↗

Impact of n > 1 neoclassical tearing modes on fast ions, plasma rotation, and the onset of disruptive tearing modes in DIII-D ITER baseline scenario

We present integrated TRANSP-kick analyses quantitatively characterizing how n > 1 magnetic islands redistribute fast ions (FIs) and modify momentum balance in low-torque DIII-D H-mode scenarios operated with the ITER normalized parameter set and shape (ITER baseline scenario, or IBS). In this plasma scenario, disruptive neoclassical tearing modes (NTMs) are most commonly seeded by nonlinear 3-wave coupling when the differential rotation between the q = 1 and q = 2 surfaces ($δf_{1, 2}$) approaches zero. Analysis of a DIII-D IBS database of shots unstable to the 2/1 NTM demonstrates that the flattening rotation profile is correlated with n > 1 NTM amplitude growth, with a strong drop in $δf_{1, 2}$ occurring when the radial magnetic perturbation amplitude at the rational surface ($δB$ r s ) is estimated to be about 20 G. The interpretive TRANSP-kick simulations show that the experimentally observed reduction in core rotation is due to the fast-ion kicks. Nonlinearity in momentum losses versus NTM amplitude is found to be correlated with overlap of FI island chains, occurring when $δB$ r s is about 35 G, in rough agreement with the experimental results. This agreement suggests the primary mechanism for the conductive momentum loss in the core is due to chaotic FI orbits and FI losses developing from overlapping resonant FI island chains.

ITER baseline scenario↗

Using ARM Data to Retrieve Entrainment Rates in Stratocumulus Cloud Systems (Final Report)

A stratocumulus-topped boundary layer (STBL) is a shallow convecting layer in which both the updraft and the downdraft branches of the convective circulation are saturated (cloudy) in the upper part of the boundary layer. The convective circulation is vigorous, yet the cloud-top entrainment velocity, $w_e$, is small and difficult to measure. Differences between currently used parameterizations of $w_e$ in GCMs can lead to STBLs that differ by as much as a factor of two in climatologically important properties such as liquid water path and boundary layer depth. We have developed a new method to retrieve entrainment rates in STBLs which uses vertical velocity retrievals from surface-based cloud radar, in conjunction with other surface-based measurements. These new estimates will allow a more comprehensive evaluation of parameterizations of entrainment rate for STBLs. The new method will also provide ground-based estimates of entrainment in STBLs over long time periods. We use retrievals of vertical velocity and liquid water content profiles to estimate the turbulent liquid water flux, which is then used in combination with observed cloud-top jumps of temperature, water vapor, and calculated radiative flux profiles to estimate the entrainment rate. We present results for three days for which observational data was not significantly contaminated by drizzle. The resulting average entrainment rates are 0.58 cm/s, 1.51 cm/s, and 0.31 cm/s.

54 ENVIRONMENTAL SCIENCES↗

Enhancements Of Interior Flux In Magnetic Bodies & Exterior Surface Fields On Conducting Objects Of Various Shapes

The magnetic flux induced in magnetic bodies of several shapes are examined. In addition the surface magnetic field intensity on the exterior surface of conducting objects is estimated since it is the driver for narrow slot aperture inductive ports of entry. In this report lower frequencies are the focus, but we also include the magnetic field intensity of a conducting cylinder near the first resonance and a reference to the high frequency limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The SDSS-Gaia View of the Color–Magnitude Relation for Blue Horizontal-branch Stars

We present an updated sample of blue horizontal-branch (BHB) stars selected from the photometric and spectroscopic data from Sloan Digital Sky Survey and its associated project Sloan Extension for Galactic Understanding and Exploration (SEGUE). With these data, we selected candidates for A-type stars in the color–color space and then a mixture modeling technique was implemented in order to distinguish between BHB and main-sequence/blue-straggler stars based on their surface gravity values (log g) estimated by the SEGUE Stellar Parameter Pipeline. Our robust approach allows us to attribute individual probabilities of each star truly being in the BHB stage. Hence, our method is advantageous in comparison to previous SEGUE BHB selections that adopted simple log g cuts. We also revisit the color–magnitude relation for these stars and propose two calibrations, based on updated distances for Galactic globular clusters, to estimate absolute magnitudes with (g – r) 0 and (u – r) 0 colors.

79 ASTRONOMY AND ASTROPHYSICS↗

Advanced Instrumentation for Metal Additive Manufacturing

Laser powder bed fusion (LPBF) is the most widely used process for metal additive manufacturing (AM), particularly where complex geometries provide performance advantages unattainable with traditional manufacturing techniques. However, LPBF is highly sensitive to innate variability in both the powder spreading and fusion steps, often leading to defects such as pores that are difficult to detect yet significantly impair component mechanical properties and fatigue life. This thesis presents a range of novel instruments enabling both precise assessment of powder layer characteristics and in-situ thermal metrology of metal AM to advance the quality control of LPBF. First, leveraging a custom X-ray microscope and a radiation-transport model developed through this work, transmission X-ray imaging is used to study spreading of thin metal powder layers. Effective layer depth is directly mapped at a process-relevant size scale, surpassing optical techniques that can only estimate local deposition from layer surface topography. Layer packing density and quality are shown to be influenced by powder flowability and particle size relative to nominal powder layer thickness. Layer quality is additionally connected to the geometry of the spreading implement and its velocity. This technique and its presented findings enable pairing feedstocks with spreading strategies that create layers with consistent packing density and uniformity. Second, a twofold approach is employed to optically interrogate the laser fusion step of LPBF for observing signatures of defect formation. Aperture division multiplexing is conceptualized, providing for simultaneous laser delivery and high-fidelity infrared (IR) process monitoring through a common optic. In-situ microscopy at 50 μm spatial resolution and at mid-wave IR wavelengths is proven readily achievable with the first purpose-built optic of this type. Next, a bespoke imaging spectrometer, along with a temperature-emissivity separation technique, is used to retrieve accurate process temperatures over a 1000 K range. Data from these instruments are correlated to porosity as fine as 4.3 μm in two LPBF test artifacts, as verified using computed tomography (CT), establishing the viability of robust optically-based component qualification.

Penny, Ryan↗

Realistic Cost to Execute Practical Quantum Circuits using Direct Clifford+T Lattice Surgery Compilation

We report a resource estimation pipeline that explicitly compiles quantum circuits expressed using the Clifford+T gate set into a surface code lattice surgery instruction set. The cadence of magic state requests from the compiled circuit enables the optimization of magic state distillation and storage requirements in a post-hoc analysis. To compile logical circuits into lattice surgery operations, we build upon the open-source Lattice Surgery Compiler. The revised compiler operates in two stages: the first translates logical gates into an abstract, layout-independent instruction set; the second compiles these into local lattice surgery instructions that are allocated to hardware tiles according to a specified resource layout. The second stage retains logical parallelism while avoiding resource contention in the fault-tolerant layer, aiding realism. Additionally, users can specify dedicated tiles at which magic states are replenished, enabling resource costs from the logical computation to be considered independently from magic state distillation and storage. We demonstrate the applicability of our pipeline to large practical quantum circuits by providing resource estimates for the ground state estimation of molecules. Finally, we find that variable magic state consumption rates in real circuits can cause the resource costs of magic state storage to dominate unless production is varied to suit.

97 MATHEMATICS AND COMPUTING↗

Scaling of Transition State Vibrational Frequencies and Application of d -Band Theory to the Brønsted–Evans–Polanyi Relationship on Surfaces

Semiempirical energy relations provide a means of estimating thermodynamic properties. Specifically, linear scaling relationships (LSRs) and Brønsted–Evans–Polanyi (BEP) relationships correlate adsorption energies between adsorbates across surfaces and reaction energies with activation barriers, respectively. Although vibrational scaling relations (VSRs) exist between adsorbates at identical sites, scaling between vibrational frequencies of adsorbed local minima and transition states is lacking. Here, we present density functional theory calculations for AH X (A = C, N, O) diffusions on transition metal surfaces and reveal linear scaling between frequencies of local minima and the transition state between those minima. Using d-band theory and linear muffin tin orbital theory (LMTO), we derive the slopes of these transition state vibrational scaling relations (TSVSRs) and, in so doing, provide a rigorous theory extending the original BEP relations developed for solution chemistry to surface chemistry. Furthermore, with a single reference DFT calculation, we predict the slopes and quantify uncertainty in the predictions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A comprehensive model of Drosophila epithelium reveals the role of embryo geometry and cell topology in mechanical responses

In order to understand morphogenesis, it is necessary to know the material properties or forces shaping the living tissue. In spite of this need, very few in vivo measurements are currently available. Here, using the early Drosophila embryo as a model, we describe a novel cantilever-based technique which allows for the simultaneous quantification of applied force and tissue displacement in a living embryo. By analyzing data from a series of experiments in which embryonic epithelium is subjected to developmentally relevant perturbations, we conclude that the response to applied force is adiabatic and is dominated by elastic forces and geometric constraints, or system size effects. Crucially, computational modeling of the experimental data indicated that the apical surface of the epithelium must be softer than the basal surface, a result which we confirmed experimentally. Further, we used the combination of experimental data and comprehensive computational model to estimate the elastic modulus of the apical surface and set a lower bound on the elastic modulus of the basal surface. More generally, our investigations revealed important general features that we believe should be more widely addressed when quantitatively modeling tissue mechanics in any system. Specifically, different compartments of the same cell can have very different mechanical properties; when they do, they can contribute differently to different mechanical stimuli and cannot be merely averaged together. Additionally, tissue geometry can play a substantial role in mechanical response, and cannot be neglected.

59 BASIC BIOLOGICAL SCIENCES↗

Using Atmospheric Inverse Modelling of Methane Budgets with Copernicus Land Water and Wetness Data to Detect Land Use-Related Emissions

Climate change mitigation requires countries to report their annual greenhouse gas (GHG) emissions and sinks, including those from land use, land use change, and forestry (LULUCF). In Finland, the LULUCF sector plays a crucial role in achieving net-zero GHG emissions, as the sector is expected to be a net sink. However, accurate estimates of LULUCF-related GHG emissions, such as methane (CH 4 ), remain challenging. We estimated LULUCF-related CH 4 emissions in Finland in 2013–2020 by combining national land cover and remote-sensed surface wetness data with CH 4 emissions estimated by an inversion model. According to our inversion model, most of Finland’s CH 4 emissions were attributed to natural sources such as open pristine peatlands. However, our research indicated that forests with thin tree cover surrounding open peatlands may also be a significant source of CH 4 . Unlike open pristine peatlands and pristine peatlands with thin tree cover, surrounding transient forests are included in the Finnish GHG inventory if they meet the criteria used for forest land. The current Finnish national GHG inventory may therefore underestimate CH 4 emissions from forested organic soils surrounding open peatlands, although more precise methods and data are needed to verify this. Given the potential impact on net GHG emissions, CH 4 emissions from transitional forests on organic soils should be further investigated. Furthermore, the results demonstrate the potential of combining atmospheric inversion modelling of GHGs with diverse data sources and highlight the need for methods to more easily combine atmospheric inversions with national GHG inventories.

54 ENVIRONMENTAL SCIENCES↗

A Multiple Filter Based Neural Network Approach to the Extrapolation of Adsorption Energies on Metal Surfaces for Catalysis Applications

Computational catalyst discovery involves the development of microkinetic reactor models based on estimated parameters determined from density functional theory (DFT). For complex surface chemistries, the number of reaction intermediates can be very large, and the cost of calculating the adsorption energies by DFT for all surface intermediates even for one active site model can become prohibitive. In this paper, we have identified appropriate descriptors and machine learning models that can be used to predict a significant part of these adsorption energies given data on the rest of them. Furthermore, our investigations also included the case when the species data used to train the predictive model are of different size relative to the species the model tries to predict—this is an extrapolation in the data space which is typically difficult with regular machine learning models. Due to the relative size of the available data sets, we have attempted to extrapolate from the larger species to the smaller ones in the current work. Here, we have developed a neural network based predictive model that combines an established additive atomic contribution based model with the concepts of a convolutional neural network that, when extrapolating, achieves a statistically significant improvement over the previous models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Direct Demonstration of Unified Brønsted-Evans-Polanyi Relationships for Proton-Coupled Electron Transfer Reactions on Transition Metal Surfaces

Brønsted-Evans-Polanyi (BEP) relationships, which relate elementary reaction barriers to reaction thermodynamics, have long been discussed in electrochemical science, but demonstration of their existence in heterogeneous electrocatalysis, across many different catalyst surfaces and voltages, is lacking. In this contribution, the BEP model is demonstrated to describe both multiple catalyst surfaces and variable voltages in electrochemical environments. For proton-coupled electron transfer (PCET) to surface nitrogen (N*) and nitric oxide (NO*), reaction energies and activation barriers are calculated using Density Functional Theory (DFT) on a parallel plate capacitor model at three different potentials and for nine different transition metal surfaces. Linear BEP relations that describe all potentials and catalyst surfaces are obtained for these elementary reactions, and the slopes of the correlations are shown to be directly related to the fractional coordinate of the transition states (FCTS) of the reactions. The results, which are explained in terms of Marcus Theory, prove a direct equivalence between unified BEP coefficients, describing both variable catalyst surfaces and voltages, and electrochemical symmetry factors and provide a straightforward means of estimating this quantify for PCET reactions on transition metal surfaces. The resulting relationships, in turn, could lead to predictions of electrocatalytic reactivity trends of enhanced accuracy and efficiency.

Electrochemistry↗

Ecosystem-Level Energy and Water Budgets Are Resilient to Canopy Mortality in Sparse Semiarid Biomes

Climate-driven woody vegetation mortality is a defining feature of semiarid biomes that drives fundamental changes in ecosystem structure. However, the observed impacts of woody mortality on ecosystem-scale energy and water budgets and the responses of surviving vegetation are highly variable among studies in water-limited environments. A previous girdling manipulation experiment in a piñon-juniper woodland suggested that although ecosystem-scale evapotranspiration was not altered by large-scale piñon mortality, soil water content decreased and the surviving juniper experienced greater water stress than juniper in an undisturbed woodland. Here we experimentally explored to what extent mortality-induced changes in energy balance components can explain these results. We compared energy fluxes measured above two adjacent piñon-juniper woodlands where piñon girdling was implemented at one site and the other subsequently experienced large-scale natural piñon mortality. We found that the mortality-induced decrease in canopy area was not sufficient to alter surface reflectance, roughness, and partitioning between energy budget components at both sites. A radiative transfer model estimated that because of the sparse premortality canopy, surface reflectance is more sensitive to a large increase in understory leaf area than further loss of crown area. Increased water stress in the remaining juniper following both mortality events can be explained by an increase in radiation on the ground that promoted higher soil temperature and evaporation. We found similar responses of ecosystem and tree-level functions to both girdling and natural mortality. This suggests that girdling is an appropriate approach to explore the impact of tree mortality on ecosystem structure, function, and energy balance.

54 ENVIRONMENTAL SCIENCES↗

SNF Interim Storage Canister Corrosion and Surface Environment Investigations (FY2020 Status Report)

This progress report describes work performed during FY20 at Sandia National Laboratories (SNL) to assess the localized corrosion performance of container/cask materials used in the interim storage of spent nuclear fuel (SNF). Of particular concern is stress corrosion cracking (SCC), by which a through-wall crack could potentially form in a canister outer wall over time intervals that are shorter than possible dry storage times. Work in FY20 further defined our understanding of the potential chemical and physical environment present on canister surfaces, evaluated the relationship between the environment and the resultant corrosion that occurs, and initiated crack growth rate testing under relevant environmental conditions. In FY20, work to define dry storage canister surface environments included several tasks. First, collection of dust deposition specimens from independent spent fuel storage installation (ISFSI) site locations helped to establish a more complete understanding of the potential chemical environment formed on the canister. Second, the predicted evolution of canister surface relative humidity RH) values was estimated using ISFSI site weather data and the horizontal canister thermal model used by the SNL probabilistic SCC model. These calculations determined that for typical ISFSI weather conditions, seasalt deliquescence to produce MgCl2-rich brines could occur in less than 20 years at the coolest locations on the canister surface, and, even after nearly 300 years, conditions for NaCl deliquescence (75% RH) are not reached. This work illustrates the importance of understanding the stability of MgCl2-rich brines on the heated canister surface, and the potential impact of brine composition on corrosion processes, including pitting and stress corrosion cracking. In an additional study, the description of the canister surface environment was refined in order to define more realistic corrosion testing environments including diurnal cycles, soluble salt chemistries, and inert mineral particles. The potential impacts of these phenomena on canister corrosion are being evaluated experimentally. Finally, work over the past few years to evaluate the stability of magnesium chloride brines continued in FY20. MgCl2 degassing experiments were carried out, confirming that MgCl2 brines slowly degas HCl on heated surfaces, converting to less deliquescent magnesium hydroxychloride phases and potentially leading to brine dryout.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Planetary Boundary-Layer Height (PBLHT) Value-Added Product: Remote-Sensing Retrievals

The planetary boundary layer (PBL) is fundamental to numerous atmospheric processes, including aerosol mixing and transport, cloud evolution, and precipitation formation. A critical parameter in these studies is the PBL height (PBLHT). This vertical depth is essential for characterizing PBL structures in numerical simulations and serves as a primary metric for estimating flux exchanges between the Earth’s surface and the atmosphere. Radiosonde (SONDE) observations provide high-vertical-resolution measurements of temperature and moisture profiles and are widely used to estimate PBLHT (Liu and Liang 2010, Seidel et al. 2010). The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility’s PBLHT value-added product (VAP) for radiosonde measurements, known as PBLHTSONDE, applies three commonly used methods—the Heffter (1980) method, the Liu and Liang (2010) method, and the bulk Richardson number approach (Seibert et al. 2000)—to derive PBLHT. The PBLHTSONDE VAP operates routinely at ARM observatories and mobile facilities, with data available from the ARM Data Center shortly after sounding observations are collected (Sivaraman et al. 2013). However, radiosonde observations are limited by their low temporal resolution. Most stations launch soundings only twice daily, which constrains the ability to investigate and characterize the temporal evolution of the PBL using radiosonde data alone. The use of continuous remote-sensing observations provides high temporal resolution of PBLHT estimates. These observations include aerosol lidars (Dang et al. 2019, Su et al. 2020), Doppler lidar (DL; Tucker et al. 2009, Krishnamurthy et al. 2021), and water vapor and/or temperature lidars and radiometers (Turner et al. 2014). These observations provide valuable data on the PBL’s thermodynamic properties (e.g., water vapor and/or temperature lidars and radiometers), dynamic properties (e.g., DL), and distribution of tracer substances (e.g., aerosol lidars), all of which can be used to estimate PBLHT. ARM developed PBLHT estimates from the micropulse lidar (MPL; PBLHTMPL), Doppler lidar (PBLHTDL), and combined Raman lidar (RL)/atmospheric emitted radiance interferometer (AERI) thermodynamic profiles (PBLHTTHERMO). Each estimate captures different physical characteristics of the boundary layer—aerosol tracers, vertical velocity turbulence, and thermodynamic structure—and exhibits distinct strengths and limitations depending on the PBL regime and time of day. In addition, the ARM ceilometer (CEIL) provides three potential PBLHT candidates derived from the vendor's built-in algorithm. Building on these individual retrievals, ARM developed the PBLHTBEML VAP, which combines the four remote-sensing-based estimates with ancillary meteorological variables using the machine learning approach of Zhang et al. (2025) to produce a best-estimate PBLHT at 10-minute resolution.

54 ENVIRONMENTAL SCIENCES↗

A computational fluid dynamics model to estimate local quantities in firebrand char oxidation

Firebrand burning is a complex phenomenon that is influenced by several parameters which are difficult to fully explore experimentally. Computational fluid dynamics models capable of predicting local quantities are essential for accurate prediction of char oxidation in firebrands. This article presents a computational fluid dynamics model to estimate firebrand mass loss, diameter change, and surface temperature during char oxidation. The model was validated using previously conducted wind tunnel experiments. These experiments were conducted for firebrands of two different aspect ratios, which were arranged in three different configurations (single, horizontal array, and vertical array), and for four different wind speeds (0.5, 1, 1.5, and 2 m/s). The computational fluid dynamics results were compared with a previous 1 D model. In all the test cases, the computational fluid dynamics model predicted the physical phenomena with significantly improved accuracy compared to a 1 D model. The char oxidation model presented in this article can be coupled with other models to study firebrand generation and trajectory, biomass pyrolysis, fluidized bed reactors, and coal combustion.

Engineering↗

CO adsorption on MnO(100): Experimental benchmarks compared to DFT

In this study, CO adsorption on the MnO(100) surface was studied using temperature programmed desorption (TPD) and density functional theory (DFT). TPD results show that CO is weakly-bound on MnO(100), with an experimental adsorption energy of -35.6 ± 2.1 kJ/mol at terrace sites in the low coverage limit. PBE simulations suggest that CO adsorption causes an implausible (2 × 2) surface reconstruction. PBE+U simulations show no signs of surface reconstruction, and provide an accurate estimate of the adsorption energy (-36.4 kJ/mol) when combined with the DFT-D3 method with Becke-Jonson damping to correct for van der Waals interactions. This simulation also shows that CO adsorbs C-down onto the Mn 2+ terrace site in a tilted geometry, which is also observed experimentally and computationally on the similarly-structured NiO(100) transition metal oxide surface. TPD results for large doses show a plateauing of the coverage at about 0.4 monolayers of CO at 85 K, with a defect coverage equivalent to 0.08 monolayers. Adsorption associated with defect sites is indicated by a high-temperature desorption tail which is not satisfactorily explained by DFT simulations of simple step or oxygen vacancy defects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗