Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase variation”

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 703 records · Page 39

MESA Isochrones and Stellar Tracks (MIST). III. The White Dwarf Cooling Sequence

We present a substantial update to the MESA Isochrones and Stellar Tracks (MIST) library, extending the MIST model grids and isochrones down the white dwarf (WD) cooling sequence with realistic physics for WD cooling timescales. This work provides a large grid of MESA models for carbon–oxygen core WDs with hydrogen atmospheres (spectral type DA/DC), descended from full prior stellar evolution calculations. The model tracks, isochrones, and WD cooling timescale contours are available on the MIST project website and at doi:10.5281/zenodo.15242047. Our WD models provide a very large, publicly available grid with detailed physics for WD cooling timescales: realistic interior and envelope compositions, with element diffusion and heavy-element sedimentation, nuclear burning at the base of the WD hydrogen envelope, core crystallization, and C/O phase separation. As a large grid of open-source stellar evolution models, these WD models provide both out-of-the-box model tracks for comparison with observations and a framework for building further WD models to investigate variations in WD physics.

Astronomy and AstroPhysics↗

Dependence of the impurity transport on the dominant turbulent regime in ELM-y H-mode discharges on the DIII-D tokamak

Laser blow-off injections of aluminum and tungsten have been performed on the DIII-D tokamak to investigate the variation of impurity transport in a set of dedicated ion and electron heating scans with a fixed value of the external torque. The particle transport is quantified via the Bayesian inference method, which, constrained by a combination of a charge exchange recombination spectroscopy, soft X-ray measurements, and VUV spectroscopy provides a detailed uncertainty quantification of transport coefficients. Contrasting discharge phases with a dominant electron and ion heating reveal 3-fold drop in the impurity confinement time and order of magnitude increase in midradius impurity diffusion, when additional electron heating is applied. Further, the calculated stationary aluminum density profiles reverse from peaked in electron heated to hollow in the ion heated case, following a similar trend as electron and carbon density. Comparable values of a core diffusion have been observed for W and Al ions, while differences in the propagation dynamics of these impurities are attributed to pedestal and edge transport. Modeling of the core transport with non-linear gyrokinetics code CGYRO, significantly underpredicts the magnitude of the variation in Al transport. Diffusion increases 3-times steeper with additional electron heat flux and 10-times lower diffusion is observed in ion heated case than predicted by the modeling. Here, the CGYRO model quantitatively matches the increase of the Al diffusion when approaching the linear threshold for the transition from the ion temperature gradient (ITG) to trapped electron mode (TEM).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multitask graph neural networks for elastoplastic response prediction in dual-phase polycrystals

Microstructure-sensitive prediction of elastoplastic response remains a recurring bottleneck in multiscale damage and fatigue modeling, where large ensembles of statistically distinct polycrystals are required to quantify variability and extreme-value behavior. In this work, we develop a multitask graph neural network (GNN) surrogate that maps dual-phase ferrite–martensite polycrystal microstructures to Statistical Volume Element (SVE)-level elastoplastic Quantities of Interest (QoIs). Each SVE is represented as a grain-adjacency graph, with node features encoding phase, geometry, and crystallographic orientation, and edge features encoding relative misorientation. A message-passing graph convolution generates node embeddings, which are pooled into a graph representation and passed to a multitask regression head that jointly predicts 10 scalar QoIs and vector-valued stress–strain responses in orthogonal loading directions across multiple martensite volume fractions and SVE sizes. Results show high accuracy for scalar QoIs and strong agreement for full stress–strain trajectories, with population envelopes reproducing both median behavior and finite-SVE variability across compositions and partition scales. A unified model trained on pooled volume-fraction data preserves most within-regime accuracy relative to regime-specific models while also capturing the broader cross-regime variation reflected in the pooled test set. Distributional comparisons further demonstrate that the surrogate preserves heterogeneity under SVE partitioning, enabling statistically consistent block-wise random-field construction for mesoscale analyses. Overall, the proposed grain-graph surrogate provides a practical pathway to accelerate ensemble-based studies of SVE-level constitutive variability in dual-phase polycrystals.

Crystal plasticity↗

Multidecadal Fluctuations in the Observed ENSO‐Tropical Cyclone Teleconnection

Abstract El Niño‐Southern Oscillation (ENSO) is a skillful predictor for seasonal tropical cyclone (TC) activity in most TC basins. This study examines recent changes in the observed ENSO‐TC teleconnection strength, as measured by ENSO modulation of hurricane frequency. We find that the ENSO‐North Atlantic TC teleconnection fluctuated over time, with the strongest relationship occurring from the 1980s to the mid‐2000s. In the western and eastern North Pacific, the ENSO‐TC teleconnection has strengthened in recent decades. Periods with a strong ENSO‐TC teleconnection are associated with more favorable environmental conditions for TCs, with higher values of genesis potential indices. Positive phases of the Atlantic Multidecadal Oscillation coincided with periods of strong ENSO‐TC teleconnections in the Atlantic and North Pacific basins. A weaker Atlantic ENSO‐TC relationship was associated with negative phases of the Pacific Decadal Oscillation and the North Atlantic Oscillation. This research reveals climate conditions that modulate ENSO's utility for seasonal TC prediction. Plain Language Summary El Niño‐Southern Oscillation (ENSO) is a useful predictor for seasonal tropical cyclone (TC) activity in many basins. Here we found that the strength of the ENSO‐TC teleconnection, represented as the correlation between ENSO and the number of hurricanes and accumulated cyclone energy, has changed in the historical record. The ENSO‐TC teleconnection in the North Atlantic fluctuated over time, with a weak relationship during the 1960s and 1970s and a strong relationship during the 1980s to mid‐2000s. Meanwhile, the ENSO‐TC teleconnection strengthened in the North Pacific in recent decades, with strong teleconnections after the 1980s in the western North Pacific and after the 2000s in the eastern North Pacific. Periods of strong ENSO‐TC teleconnections are associated with more favorable environmental conditions for TCs, including higher values of genesis potential indices and higher mid‐tropospheric humidity, as well as positive phases of the Atlantic Multidecadal Oscillation. Additionally, the negative phase of the Pacific Decadal Oscillation leads to strong/weak ENSO‐TC teleconnections in the eastern North Pacific and North Atlantic, respectively. Furthermore, a negative North Atlantic Oscillation is associated with a weak ENSO‐North Atlantic TC teleconnection. This research highlights variations in ENSO's effectiveness for seasonal TC prediction. Key Points The observed impact of ENSO on tropical cyclone (TC) activity exhibits multidecadal fluctuations The ENSO‐TC teleconnection was strong in the Atlantic from the 1980s to mid‐2000s and strengthened over the North Pacific in recent decades The ENSO‐TC teleconnection is stronger in the Atlantic and North Pacific basins during a positive Atlantic Multidecadal Oscillation

ENSO↗

Temperature dependence of London penetration depth anisotropy in superconductors with anisotropic order parameters

Here, we study the effects of anisotropic order parameters on the temperature dependence of London penetration depth anisotropy γ λ ( T ) . After MgB 2 , this dependence is commonly attributed to distinct gaps on multiband Fermi surfaces in superconductors. We have found, however, that the anisotropy parameter may depend on temperature also in one-band materials with anisotropic order parameters Δ ( T , k F ) ; a few such examples are given. We have also found that for different order parameters, the temperature dependence of Δ ( T ) / Δ ( 0 ) can be represented with good accuracy by the interpolation suggested by Einzel, which simplifies considerably the evaluation of γ λ ( T ) . Of particular interest are mixed order parameters of two symmetries for which γ λ ( T ) may go through a maximum for a certain relative weight of two phases. Also, for this case we find that the ratio Δ max ( 0 ) / T c may exceed substantially the weak-coupling limit of 1.76. It, however, does not imply strong coupling; rather, it is due to significantly anisotropic angular variation of Δ .

36 MATERIALS SCIENCE↗

Testing for pair density wave order in La 1.875 Ba 0.125 CuO 4

We report charge order is commonly believed to compete with superconducting order. An intertwined form of a superconducting wave function, known as pair density wave (PDW) order, has been proposed; however, direct evidence, theoretical or experimental, that it forms the ground state of any cuprate superconductor is lacking. As a test case, we consider La 2-x Ba x CuO 4 with x = 1/8, where charge and spin stripe orders within the CuO 2 planes compete with three-dimensional superconducting order. We report measurements of the superconducting critical current perpendicular to the planes in the presence of an in-plane magnetic field. The variation of the critical current with orientation of the field is inconsistent with a theoretical prediction specific to the PDW model. It appears, instead, that the orientation dependence of the critical-current density might be determined by a minority phase of d-wave superconductivity that is present as a consequence of doped-charge inhomogeneity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Radial-velocity Search for Binary RR Lyrae Variables

We report 272 radial velocities for 19 RR Lyrae variables. For most of the stars we have radial velocities for the complete pulsation cycle. These data are used to determine robust center-of-mass radial velocities that have been compared to values from the literature in a search for evidence of binary systems. Center-of-mass velocities were determined for each star using Fourier Series and template fits to the radial velocities. Our center-of-mass velocities have uncertainties from ±0.16 km s{sup −1} to ±2.5 km s{sup −1}, with a mean uncertainty of ±0.92 km s{sup −1}. We combined our center-of-mass velocities with values from the literature to look for deviations from the mean center-of-mass velocity of each star. Fifteen RR Lyrae show no evidence of binary motion (BK And, CI And, Z CVn, DM Cyg, BK Dra, RR Gem, XX Hya, SZ Leo, BX Leo, TT Lyn, CN Lyr, TU Per, U Tri, RV UMa, and AV Vir). In most cases this conclusion is reached due to the sporadic sampling of the center-of-mass velocities over time. Three RR Lyrae show suspicious variation in the center-of-mass velocities that may indicate binary motion but do not prove it (SS Leo, ST Leo, and AO Peg). TU UMa was observed by us near a predicted periastron passage (at 0.14 in orbital phase) but the absence of additional center-of-mass velocities near periastron makes the binary detection, based on radial velocities alone, uncertain. Two stars in our sample show Hγ emission in phases 0.9–1.0: SS Leo and TU UMa.

47 OTHER INSTRUMENTATION↗

Development of Risk Mitigation Guidance for Hydrogen Sensor Placement Indoors and Outdoors

Guidance on Sensor Placement remains one of the top priorities for the safe deployment of hydrogen and fuel cell equipment in the commercial marketplace. Building on the success of Phase 1 work reported at ICHS2019 and published in IJHE, this paper discusses the consecutive steps to further develop and validate such guidance for mechanically ventilated enclosures. The key step included a more in-depth analysis of sensitivity to variation of physical parameters in a small enclosure, and finally, expansion of the developed approach to confined spaces in an outdoor environment.

codes and standards↗

Unprecedented non-hysteretic superelasticity of [001]-oriented NiCoFeGa single crystals

Superelasticity associated with the martensitic transformation has found a broad range of engineering applications. However, the intrinsic hysteresis and temperature sensitivity of the first-order phase transformation significantly hinder the usage of smart metallic components in many critical areas. In this paper, we report a large superelasticity up to 15.2% strain in [001]-oriented NiCoFeGa single crystals, exhibiting non-hysteretic mechanical responses, a small temperature dependence and high-energy-storage capability and cyclic stability over a wide temperature and composition range. In situ synchrotron X-ray diffraction measurements show that the superelasticity is correlated with a stress-induced continuous variation of lattice parameter accompanied by structural fluctuation. Neutron diffraction and electron microscopy observations reveal an unprecedented microstructure consisting of atomic-level entanglement of ordered and disordered crystal structures, which can be manipulated to tune the superelasticity. The discovery of the large elasticity related to the entangled structure paves the way for exploiting elastic strain engineering and development of related functional materials. NiCoFeGa single crystals exhibit large non-hysteretic superelasticity over broad temperature and composition ranges. It is attributed to the continuous phase transition with applied stress, which is related to the fluctuation of entangled ordered and disordered crystal structures.

36 MATERIALS SCIENCE↗

Mass enhancement in 3 d and s - p perovskites from symmetry breaking

In some d-electron oxides, the measured effective mass m* exptl has long been known to be significantly larger than the model effective mass m* model deduced from mean-field band theory, i.e., m* exptl =βm* model , where β>1 is the “mass-enhancement” or “mass-renormalization” factor. Previous applications of density functional theory (DFT), based on a symmetry-restricted structure with the smallest number of possible magnetic, orbital, and structural degrees of freedom, missed such mass enhancement. This fact has been taken as evidence of strong electronic correlation, often described via the symmetry-restricted dynamic mean-field approach of the many-body theory, being the exclusive enabling physics. This paper uses instead a static density functional approach that does not restrict positional or spin degrees of freedom (symmetry-broken structures). This approach analyzes the contributions of different symmetry-broken modalities to mass enhancement for a few nominally highly correlated d-electron perovskites as well as the nominally uncorrelated, closed-shell s-p bonding perovskites. It shows that the energy-lowering symmetry-broken spin effects (e.g., nonzero local moment in the paramagnetic phase) and structural effects (e.g., atomic displacement) as described in mean-field DFT already manifest mass enhancement for both electrons and holes in a range of d-electron perovskites SrVO 3 , SrTiO 3 , BaTiO 3 , and LaMnO 3 , as well as p-electron perovskites CsPbI 3 and SrBiO 3 , including both metals (SrVO 3 ) and insulators (the rest). This is revealed only when enlarged unit cells of the same parent global symmetry, which are large enough to allow for symmetry-breaking distortions and concomitant variations in spin order, are explored for their ability to lower the total energy. Positional symmetry breaking that leads to mass enhancement includes octahedral rotation in halide perovskites such as CsPbI 3 , Jahn-Teller-like Q${+}\atop{2}$ distortion in LaMnO 3 , and bond disproportionation in SrBiO 3 , while magnetic symmetry breaking resulting in mass enhancement includes the formation of a distribution of local moments in SrVO3 that averages to zero in the paramagnetic phase. Not all symmetry breaking leads to significant mass enhancement, e.g., the rather small octahedral rotations in the nearly perfectly cubic SrTiO 3 cause negligible mass enhancement, as do the paraelectric displacements in BaTiO 3 . In principle, by ergodicity, the two descriptions, i.e., the symmetry-restricted dynamic approach with a single time-fluctuating magnetic moment and the symmetry-broken mean-field approach with a static distribution of spatially fluctuated local moments, are not mutually exclusive but are a choice of representation and consequently, a choice of computational efficiency. In approximate implementations, the symmetry-broken mean-field approach appears to remove much of what was strong correlation in dynamically correlated symmetry-restricted solutions, leaving smaller (“weak”) residual correlation with respect to the exact solution.

36 MATERIALS SCIENCE↗

Assessing free tropospheric quasi-equilibrium for different GCM resolutions using a cloud-resolving model simulation of tropical convection

Abstract This study examines the free-tropospheric quasi-equilibrium at different global climate model (GCM) resolutions using the simulation of tropical convection by a cloud-resolving model during the Tropical Western Pacific International Cloud Experiment. The simulated dynamic and thermodynamic fields within the model domain are averaged over subdomains of different sizes equivalent to different GCM resolutions. These coarse-grained fields are then used to compute CAPE and its change with time, and their relationships with simulated convection. Results show that CAPE change with time is controlled predominantly by variations of thermodynamic properties in the planetary boundary layer for all subdomain sizes ranging from 64 to 4 km. Lag correlation analysis shows that CAPE generation by the free-tropospheric dynamical advection (dCAPE ls ) leads convective precipitation but is in phase with convective mass flux at 600 mb and 500 mb vertical velocity for all subdomain sizes. However, the correlation coefficients and regression slopes decrease as the subdomain size decreases for subdomain sizes smaller than 16 km. This is probably due to increased randomness of convection and more scale-dependence of the relationships when the subdomain size reaches the grey zone. By examining the sensitivity of the relationships of convection with dCAPE ls to temporal scales in different subdomain size, it shows that the quasi-equilibrium between dCAPE ls and convection holds well for timescales of 30 min or longer at all subdomain sizes. These results suggest that the free tropospheric quasi-equilibrium assumption may still be useable even for GCM resolutions in the grey zone.

54 ENVIRONMENTAL SCIENCES↗

Local residual stresses and microstructure within recrystallizing grains in iron

The intragranular strains and the microstructure within recrystallizing grains in single phase pure iron have been systematically studied in 3D using synchrotron micro-diffraction. In contradiction to common knowledge we observe that residual elastic strains of up to 1–2 × 10 –3 are present in the recrystallizing grains. The local strain variations within individual grains can be similar in magnitude to those between different grains. The effects of processing parameters, including rolling reduction, recrystallized volume fraction, and cooling rate, as well as microstructural parameters, including grain size, orientation and internal misorientation, on the development of local residual elastic strains and microstructure have been quantified. It is suggested that the difference in local defect density between recrystallized grains and the surrounding deformed matrix, as well as the redistribution of medium to long-range residual stresses in the deformed grains, are the main reasons for the development of such strains in the recrystallized grains. The local residual stress levels are calculated and their possible impact on the mechanical properties and recrystallization processes are discussed.

36 MATERIALS SCIENCE↗

Toward Energy-Efficient HPC: Insights from Power Profiling a Cloud-Resolving Earth System Model

Power is a fundamental constraint as supercomputing advances to exascale. Efficient operation within strict power budgets requires application-aware power management based on a detailed understanding of application-level power behavior. This work analyzes the Energy Exascale Earth System Model (E3SM) atmosphere component, SCREAM, on Perlmutter (NERSC) and Frontier (OLCF). We characterize power variation across inputs, concurrency levels, and power caps, evaluate the energy impact of code optimizations, and attribute energy within the code using a newly developed GPU energy model. Results show that SCREAM’s peak power remains stable during its core execution phase and decreases gradually as concurrency increases. Power capping experiments reveal a performance–energy "sweet spot". On Perlmutter, limiting GPU power to 50% of thermal design power (TDP) achieves up to 15% energy savings with a 7% performance penalty. On Frontier, a 40% TDP cap yields up to 10% energy savings with less than 10% performance loss. Code optimizations reduce SCREAM energy by shortening run time without increasing power. Modeling reveals a critical insight: data movement accounts for approximately 70% of SCREAM’s GPU energy. This fundamentally shifts the optimization focus from FLOPS to data transfer reduction for this class of applications, offering the most impactful strategy for improving energy efficiency. This work establishes a foundation for practical, application-aware power management at exascale.

Zhao, Zhengji [Lawrence Berkeley National Laborato↗

Doping Metallic Grain Boundaries to Control Atomic Structure and Damage Tolerance

Grain boundaries often act as sites for crack and void nucleation during plastic deformation of metallic materials. While it is known that grain boundary character and structure can greatly influence this damage nucleation process, the current level of control over such details is limited. The objective of this project was to obtain a fundamental understanding of how metallic grain boundary structure can be controlled through intelligent doping, with the idea of inducing planned amorphous grain boundary phases or complexions. The effect of amorphous complexion structure on dislocation accommodation mechanisms was studied, to improve the field’s understanding of damage nucleation at a promising type of interface. Different microstructural descriptors were studied, as grain boundaries can have large variations and complexity of local structure, and influence the mechanical damage resistance of these features was tested. Nanocrystalline systems that contain amorphous grain boundary complexions were prioritized, as these have more extreme variations in interfacial structure and damage tolerance, with an emphasis on isolating the importance of complexion population (type, thickness, etc.), network topology, local structure, and local chemistry. This research used a combination of computational, experimental, and characterization techniques to isolate and understand the importance of nanoscale grain boundary structure and interfacial chemistry. Amorphous grain boundary complexions were found to clearly increase a material’s resistance to mechanical damage, with the local distribution of structural short-range order within the complexions found to be an important descriptor for damage and the proposed focus of future work in this area.

36 MATERIALS SCIENCE↗

Neural Network Enhanced RKPM for Electrochemical-Mechanical Coupled Damage Modeling of Energy Storage Materials

Energy storage materials undergo significant charge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. For Li-ion batteries, for example, lithium moving between electrodes during charging and discharging process causes expansion and contraction of grains, and the strongly anisotropic and nonlinearly [Li]-dependent grain material properties can cause grains to expand into and contract away from each other, leading to chemo-mechanical cracking. In the first part of this work, a RKPM based computational framework for solving the coupled solid-phase lithium conservation with Fickian diffusion and the lithium concentration dependent anisotropic mechanical problem subjected to a highly nonlinear Butler-Volmer boundary condition is introduced. The choice of RKPM completeness conditions for lithium concentration and mechanical deformation fields, and the variational consistency condition for the domain integration of the coupled problem is first determined. In the second part of this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material, by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that is necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure.

damage modeling↗

Ultrasonic characterization of Ti-5Al-5Mo-5V-3Cr

Metastable β titanium alloys such as Ti-5Al-5Mo-5V-3Cr are of increasing interest due to their excellent corrosion resistance, high strength, and in particular, the ability to manipulate microstructure to a high degree to control the mechanical properties of the material. It is therefore important to consider how this manipulation of mechanical properties will affect the inspection of this material with ultrasonic testing both for the purposes of flaw detection and material characterization. The connection between material properties and ultrasonic characterization techniques will be discussed and then dissected in the context of Ti-5553. This metastable alloy was shown to exhibit remarkable variation in microstructure as well as ultrasonic velocity and attenuation demonstrating these techniques as potential candidates for characterizing Ti-5553 and Ti-5553 as a good candidate for aiding in future studies aimed at learning more about the role of secondary phases in ultrasonic wave propagation.

Sunderman, Ruth↗

Grafted nickel-promoter catalysts for dry reforming of methane identified through high-throughput experimentation

High-throughput synthesis of a series of monometallic and bimetallic catalysts (45 bimetallic and 50 monometallic samples) consisting of nickel and one of nine different metal promoters (B, Co, Cu, Fe, Mg, Mn, Sn, V and Zn) supported on one of six different metal oxides alumina, ceria, magnesia, silica and titania) is carried out via organometallic grafting using a robotic platform. The catalysts are evaluated for their activity and selectivity for the dry reforming of methane at a feed ratio of CH 4 :CO 2 of 1 at 650–800 °C in a parallel flow reactor system. The type of oxide support prevails over the type of additive for both catalyst activity and stability. On Al 2 O 3 and MgO, Fe was found to be the best promoter; on SiO 2 , Cu is the best promoter at 700 °C and higher, while on TiO 2 , Mn is found to enhance the conversion at 800 °C. On CeO 2 , all additives except Fe have beneficial effects. Twenty-five catalysts show > 90% methane conversion with ten catalysts showing > 95% conversion at 800 °C with the H 2 :CO ratios ranging from 0.8 to 1.2. Amongst the ten highest performers, NiFe/Al 2 O 3 and NiFe/MgO are more active than Ni/Al 2 O 3 and Ni/MgO, respectively and were stable over a period of 25 h at 800 °C. Characterization on the as-prepared samples reveals highly dispersed phase, while after reduction in H 2 , highly dispersed and reduced nickel particles up to 10 nm are formed. The particles do not increase in size under dry reforming reaction conditions at 800 °C. An increased hydrogen consumption observed during H 2 -TPR of the nickel particles is positively correlated with methane conversion for Al 2 O 3 -based catalysts. The resistance to deactivation by coking and variation in coke structure are investigated by spectroscopic and microscopic methods to identify the relationship between metal promoters, alloy formation, and type of surface carbon deposits. Carbon whiskers were observed on the ten selected spent samples and are preferentially deposited on Ni rather than on the promoters. Carbon nanotube formation and metal particle removal from support were not observed to cause deactivation while amorphous carbon formation was clearly linked to catalyst deactivation, as amorphous carbon could encapsulate nickel, either on the support or at the end of the carbon nanotube. Furthermore, the organometallic grafting technique is an efficient and suitable technique for synthesizing highly dispersed and homogeneous phases which lead to high conversion and high durability for dry reforming of methane.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗