Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dependability modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,693 records · Page 94

An architecture and model for cognitive engineering simulation analysis - Application to advanced aviation automation

The process of designing crew stations for large-scale, complex automated systems is made difficult because of the flexibility of roles that the crew can assume, and by the rapid rate at which system designs become fixed. Modern cockpit automation frequently involves multiple layers of control and display technology in which human operators must exercise equipment in augmented, supervisory, and fully automated control modes. In this context, we maintain that effective human-centered design is dependent on adequate models of human/system performance in which representations of the equipment, the human operator(s), and the mission tasks are available to designers for manipulation and modification. The joint Army-NASA Aircrew/Aircraft Integration (A3I) Program, with its attendant Man-machine Integration Design and Analysis System (MIDAS), was initiated to meet this challenge. MIDAS provides designers with a test bed for analyzing human-system integration in an environment in which both cognitive human function and 'intelligent' machine function are described in similar terms. This distributed object-oriented simulation system, its architecture and assumptions, and our experiences from its application in advanced aviation crew stations are described.

Corker, Kevin M.↗

High-Resolution Mapping of Lunar Crustal Magnetic Fields: Correlations with Albedo Markings of the Reiner Gamma Class

During the last eight months of the Lunar Prospector mission (December 1999-July 1999), the spacecraft was placed in a relatively low-altitude (15-30-km perapsis), near-polar orbit that allowed high-resolution mapping of crustal magnetic fields. We report here initial studies of the correlation of locally strong magnetic anomalies with unusual, swirl-like albedo markings of the Reiner Gamma class. Based on this correlation, which is known from earlier studies of Apollo subsatellite magnetometer data, it has been proposed that the swirls represent regions whose higher albedos have been preserved via deflection of the solar-wind ion bombardment by strong crustal fields. This model in turn depends on the hypothesis that solar-wind implanted H is at least one component of the process that optically matures exposed silicate surfaces in the inner solar system . Specifically, it is hypothesized that implanted H acts as an effective reducing agent to enhance the rate of production of nanophase metallic Fe particles from preexisting silicates during micrometeoroid impacts. According to the model, the curvilinear shapes of these albedo markings are caused, at least in part, by the geometry of ion deflections in a magnetic field. The improved resolution and coverage of the Prospector data allow more detailed mapping of the fields, especially on the lunar farside. This permits a more quantitative test of whether all albedo markings of this class are associated with strong local magnetic fields.Only if the latter condition is met can the solar-wind deflection hypothesis he valid. The basic procedure for mapping crustal magnetic fields using Lunar Prospector magnetometer data follows that developed for analysis of Apollo subsatellite magnetometer data. The specific mapping steps are (1) selection of mission time intervals suitable for mapping crustal fields; these are limited essentially either to times when the Moon is in a lobe of the geomagnetic tail or to times when the Moon is in the solar wind but the spacecraft is in the lunar wake; the data are transformed to a radial, east, and north coordinate system with measurements given as a function of spacecraft latitude, longitude, and altitude; (2) visual editing of individual orbit segments selected for minimal external field disturbances; (3) minimization of remaining low-frequency external fields for individual orbit data segments by quadratic detrending; and (4) two-dimensional filtering of individual orbit segments to produce a vector field map along the slightly curved surface defined by the spacecraft altitude; maps of the three field components (radial, east, and north), the field magnitude, and the spacecraft altitude are constructed. For data obtained at low to middle latitudes, the horizontal resolution of the field maps is limited by the orbit-track separation (about 30 km at the equator). Maps of the field magnitude have been constructed within limited selenographic regions based mainly on data acquired in March, April, and May of 1999. This was a time period when the orbit plane was nearly aligned with the Sun-Moon line so that field mapping was possible at times when the Moon was in the solar wind as well as when the Moon was in the geomagnetic tail. Most of the coverage is across the lunar farside. However, a shows an example of a field map produced from solar-wind wake data for a region including Reiner Gamm on western Oceanus Procellarum (location: 58.5W, 7.5N). The contour interval is 3 nT and the mean spacecraft altitude is 18 km to within the accuracy allowed by the resolution of the map (30 km or about 1 deg.); strong magnetic anomalies correlate closely with swirl locations. Individual orbit profiles (whose resolution along the orbit track is comparable to the spacecraft altitude of 18 km) also demonstrate a good correlation of field magnitude with surface albedo. In order to investigate the correlation of magnetic fields with the location of swirl features, we have reexamined available lunar imagery (Lunar Orbiter, Apollo, and Clementine) to identify and map swirl locations within regions where swirls have previously been mapped. In these images, swirls were distinguished from other high-albedo features such as crater rays by their curvilinear shapes and increased visibility in forward-scattered light. Digital maps of swirls identified by all available imagery were then superposed on maps of the field magnitude at the spacecraft altitude. Based upon analysis of these composite magnetic/geologic maps, we draw the preliminary conclusion that swirl features are associated with magnetic anomalies revealed by Lunar Prospector. Detailed maps of these swirl features are currently being constructed for the magnetically strong regions antipodal to the Imbrium, Serenitatis, and Crisium Basins. Additional information contained in the original,

Hood, L. L.↗

Removing Numerical Pathologies in a Turbulence Parameterization Through Convergence Testing

Abstract Discretized numerical models of the atmosphere are usually intended to faithfully represent an underlying set of continuous equations, but this necessary condition is violated sometimes by subtle pathologies that have crept into the discretized equations. Such pathologies can introduce undesirable artifacts, such as sawtooth noise, into the model solutions. The presence of these pathologies can be detected by numerical convergence testing. This study employs convergence testing to verify the discretization of the Cloud Layers Unified By Binormals (CLUBB) model of clouds and turbulence. That convergence testing identifies two aspects of CLUBB's equation set that contribute to undesirable noise in the solutions. First, numerical limiters (i.e., clipping) used by CLUBB introduce discontinuities or slope discontinuities in model fields. Second, nonlinear artificial diffusion employed for improving numerical stability can introduce unintended small‐scale features into the solution of the model equations. Smoothing the limiters and using linear artificial diffusion reduces the noise and restores the expected first‐order convergence in CLUBB's solutions. These model reformulations enhance our confidence in the trustworthiness of solutions from CLUBB by eliminating the unphysical oscillations in high‐resolution simulations. The improvements in the results at coarser, near‐operational grid spacing and timestep are also seen in cumulus cloud and dry turbulence tests. In addition, convergence testing is proven to be a valuable tool for detecting pathologies, including unintended discontinuities and grid dependence, in the model equation set.

58 GEOSCIENCES↗

Time Dilated Bundt Cake Analysis of PV Output [Poster]

We present a novel method for modeling time-dependent statistics in the power signal generated by a photovoltaic (PV) system. Our white-box machine learning method is interpretable and auditable, based on principles of multiperiodic basis functions and convex optimization. Our proposed method of time dilating the daily signal to remove night time values results in a novel representation of PV power signals, evocative of a ‘Bundt cake’. The proposed model describes the marginal distribution of power output as a function of date and time. The resulting probabilistic model of a PV system can be used to perform a variety of tasks, and here, we demonstrate the application of clear sky detection.

14 SOLAR ENERGY↗

Quantifying the Effect of Pore‐Size Dependent Wettability on Relative Permeability Using Capillary Bundle Model

Abstract Relative permeability is a key parameter for characterizing the multiphase flow dynamics in porous media at macroscopic scale while it can be significantly impacted by wettability. Recently, it has been reported in microfluidic experiments that wettability is dependent on the pore size (Van Rooijen et al., 2022). To investigate the effect of pore‐size‐dependent wettability on relative permeability, we propose a theoretical framework informed by digital core samples to quantify the deviation of relative permeability curves due to wettability change. We find that the significance of impact is highly dependent on two factors: (i) the function between contact angle and pore size (ii) overall pore size distribution. Under linear function, this impact can be significant for tight porous media with a maximum deviation of 1,000%.

Yu, Siqin↗

A priori analysis of a power-law sub-filter scale mixing model for transported PDF modeling of high Karlovitz turbulent premixed flames

Accurate modeling of mixing in large-eddy simulation (LES) / transported probability density function (PDF) modeling of turbulent combustion remains an outstanding issue. The issue is particularly salient in turbulent premixed combustion under extreme conditions such as high-Karlovitz number Ka. Here, the present study addresses this issue by conducting an a priori analysis of a power-law scaling based mixing timescale model for the transported PDF model. A recently produced DNS dataset of a high-Ka turbulent jet flame is used for the analysis. A power-law scaling is observed for a scaling factor used to model the sub-filter scale mixing timescale in this high-Ka turbulent premixed DNS flame when the LES filter size is much greater than the characteristic thermal thickness of a laminar premixed flame. The sensitivity of the observed power-law scaling to the different viewpoints (local or global) and to the different scalars for the data analysis is examined and the dependence of the model parameters on the dimensionless numbers Ka and Re (the Reynolds number) is investigated. Different model formulations for the mixing timescale are then constructed and assessed in the DNS flame. The proposed model is found to be able to reproduce the mixing timescale informed by the high-Ka DNS flame significantly better than a previous model.

42 ENGINEERING↗

A comparison of viscoelastic damping models

Modern finite element methods (FEM's) enable the precise modeling of mass and stiffness properties in what were in the past overwhelmingly large and complex structures. These models allow the accurate determination of natural frequencies and mode shapes. However, adequate methods for modeling highly damped and high frequency dependent structures did not exist until recently. The most commonly used method, Modal Strain Energy, does not correctly predict complex mode shapes since it is based on the assumption that the mode shapes of a structure are real. Recently, many techniques have been developed which allow the modeling of frequency dependent damping properties of materials in a finite element compatible form. Two of these methods, the Golla-Hughes-McTavish method and the Lesieutre-Mingori method, model the frequency dependent effects by adding coordinates to the existing system thus maintaining the linearity of the model. The third model, proposed by Bagley and Torvik, is based on the Fractional Calculus method and requires fewer empirical parameters to model the frequency dependence at the expense of linearity of the governing equations. This work examines the Modal Strain Energy, Golla-Hughes-McTavish and Bagley and Torvik models and compares them to determine the plausibility of using them for modeling viscoelastic damping in large structures.

Slater, Joseph C.↗

Redshift-space streaming velocity effects on the Lyman- α forest baryon acoustic oscillation scale

The baryon acoustic oscillation (BAO) scale acts as a standard ruler for measuring cosmological distances and has therefore emerged as a leading probe of cosmic expansion history. However, any physical effect that alters the length of the ruler can lead to a bias in our determination of distance and expansion rate. One of these physical effects is the streaming velocity, the relative velocity between baryons and dark matter in the early Universe, which couples to the BAO scale due to their common origin in acoustic waves at recombination. In this work, we investigate the impact of streaming velocity on the BAO feature of the Lyman- α forest auto-power spectrum, one of the main tracers being used by the recently commissioned Dark Energy Spectroscopic Instrument (DESI). To do this, we develop a perturbative model for Lyman- α flux fluctuations which is complete to second order for a certain set of fields, and applicable to any redshift-space tracer of structure since it is based only on symmetry considerations. We find that there are 8 biasing coefficients through second order. We find streaming velocity-induced shifts in the BAO scale of 0.081%–0.149% (transverse direction) and 0.053%–0.058% (radial direction), depending on the model for the biasing coefficients used. These are smaller than, but not negligible compared to, the DESI Lyman- α BAO error budget, which is 0.46% on the overall scale. Overall, the sensitivity of these results to our choice of bias parameters underscores the need for future work to measure the higher-order biasing coefficients from simulations, especially for future experiments beyond DESI.

79 ASTRONOMY AND ASTROPHYSICS↗

The Size Evolution of Passive Galaxies: Observations From the Wide-Field Camera 3 Early Release Science Program

We present the size evolution of passively evolving galaxies at z approximately 2 identified in Wide-Field Camera 3 imaging from the Early Release Science program. Our sample was constructed using an analog to the passive BzK galaxy selection criterion, which isolates galaxies with little or no ongoing star formation at z greater than approximately 1.5. We identify 30 galaxies in approximately 40 arcmin(sup 2) to H less than 25 mag. By fitting the 10-band Hubble Space Telescope photometry from 0.22 micrometers less than approximately lambda (sub obs) 1.6 micrometers with stellar population synthesis models, we simultaneously determine photometric redshift, stellar mass, and a bevy of other population parameters. Based on the six galaxies with published spectroscopic redshifts, we estimate a typical redshift uncertainty of approximately 0.033(1+z).We determine effective radii from Sersic profile fits to the H-band image using an empirical point-spread function. By supplementing our data with published samples, we propose a mass-dependent size evolution model for passively evolving galaxies, where the most massive galaxies (M(sub ∗) approximately 10(sup 11) solar mass) undergo the strongest evolution from z approximately 2 to the present. Parameterizing the size evolution as (1 + z)(sup - alpha), we find a tentative scaling of alpha approximately equals (−0.6 plus or minus 0.7) + (0.9 plus or minus 0.4) log(M(sub ∗)/10(sup 9 solar mass), where the relatively large uncertainties reflect the poor sampling in stellar mass due to the low numbers of highredshift systems. We discuss the implications of this result for the redshift evolution of the M(sub ∗)-R(sub e) relation for red galaxies.

fundamental parameters↗

Preliminary Result of Neoclassical Impurity Transport Affected by ICRH with a 4-D Fokker-Planck Code

To understand the role of ion cyclotron resonance heating (ICRH) in mitigating tungsten accumulation, we evaluate the impact of ICRH-induced poloidal asymmetries on neoclassical impurity transport using a newly developed 4-D Fokker–Planck code (FP4D). Unlike conventional approaches that separate RF heating and transport modelling, FP4D solves both fast-ion distribution formation and impurity transport self-consistently without bounce-averaging. Benchmarking against the NEO code confirms the accuracy of the neoclassical physics in FP4D. Using FP4D, we analyze anisotropic and non-Maxwellian minority distributions under ICRH and their effects on the equilibrium potential and tungsten transport. Results show that anisotropy induced poloidal potential variations can reduce inward tungsten flux, depending on the modelling approach. Our findings highlight the importance of resolving the full velocity and poloidal structure of minority species in predicting impurity transport under strong RF heating conditions.

Jeong, Yunho [Hanyang Univ., Seoul (Korea, Republi↗

Comparative Study of Advanced Turbulence Models for Turbomachinery

A computational study has been undertaken to study the performance of advanced phenomenological turbulence models coded in a modular form to describe incompressible turbulent flow behavior in two dimensional/axisymmetric and three dimensional complex geometry. The models include a variety of two equation models (single and multi-scale k-epsilon models with different near wall treatments) and second moment algebraic and full Reynolds stress closure models. These models were systematically assessed to evaluate their performance in complex flows with rotation, curvature and separation. The models are coded as self contained modules that can be interfaced with a number of flow solvers. These modules are stand alone satellite programs that come with their own formulation, finite-volume discretization scheme, solver and boundary condition implementation. They will take as input (from any generic Navier-Stokes solver) the velocity field, grid (structured H-type grid) and computational domain specification (boundary conditions), and will deliver, depending on the model used, turbulent viscosity, or the components of the Reynolds stress tensor. There are separate 2D/axisymmetric and/or 3D decks for each module considered. The modules are tested using Rocketdyn's proprietary code REACT. The code utilizes an efficient solution procedure to solve Navier-Stokes equations in a non-orthogonal body-fitted coordinate system. The differential equations are discretized over a finite-volume grid using a non-staggered variable arrangement and an efficient solution procedure based on the SIMPLE algorithm for the velocity-pressure coupling is used. The modules developed have been interfaced and tested using finite-volume, pressure-correction CFD solvers which are widely used in the CFD community. Other solvers can also be used to test these modules since they are independently structured with their own discretization scheme and solver methodology. Many of these modules have been independently tested by Professor C.P. Chen and his group at the University of Alabama at Huntsville (UAH) by interfacing them with own flow solver (MAST).

Hadid, Ali H.↗

Probabilistic Failure Criterion of SiC/SiC Composites Under Multiaxial Loading

Owing to its excellent mechanical properties and stability under high temperature and neutron irradiation conditions, SiC/SiC composites have emerged as a promising material for light water reactors (LWRs) in the development of accident-tolerant fuel (ATF) systems. Structural integrity and retention of hermeticity are two crucial requirements for SiC/SiC claddings during normal operations, and both of them are closely related to the proportional limit stress (PLS) of the material. Understanding the behavior of SiC/SiC composites under multiaxial stress states and developing a probabilistic approach for evaluating the structural vulnerability are of paramount importance for reliability-based analysis and design of SiC/SiC composite claddings. So far, there has been very limited effort towards experimental and analytical investigations of probabilistic failure of SiC/SiC claddings. This critical knowledge gap motivates this research. A probabilistic failure criterion for SiC/SiC composites under multi-axial loading is developed, and this criterion is incorporated into reliability analysis of the structural integrity of SiC/SiC fuel cladding. The research consists of two parts: 1) experimental investigation of multiaxial failure behavior of SiC/SiC composites, and 2) theoretical modeling of time-dependent probabilistic failure of SiC/SiC cladding. In the experimental investigation, the PLS is determined through the examination of stress-strain response, the acoustic emission measurement, as well as the X-ray computed tomography. The theoretical framework is derived by combin- ing the finite weakest-link statistical model and the subcritical damage growth model. This theoretical model captures the time-dependent failure mechanism of the material, which has a major consequence for predicting the lifetime distribution of the cladding. Meanwhile, the model also predicts that the failure statistics of the cladding depends strongly on the cladding length. The results of the multiaxial experiments reveal the level of statistical variation of the PLS of SiC/SiC materials under different stress states. The theoretical model provides a robust analytical tool for extrapolation of small-scale laboratory test results to the behavior of full-scale claddings. These findings establish a scientific foundation for the development of reliability-based design of SiC/SiC fuel claddings, which will play an essential role in improving the structural safety and integrity of LWRs.

42 ENGINEERING↗

Contribution to polar albedo from a mesospheric aerosol layer

An examination is made of the impact of a layer of particulate matter, assumed to be ice crystals, on the albedo of the polar region. The model is time dependent, includes the growth of the layer, and incorporates the diffuse nature of radiation reflected from the surface and atmosphere. Although the magnitude of the effect is about an order of magnitude less than previous results, the impact is one of heating instead of cooling. It is also shown that ignoring the diffuse nature of the radiation reflected from the underlying earth-atmosphere system, as has been done in many previous simple models, can result in overestimation of the climatological impact of aerosols in sign and magnitude by a factor of up to 4-6.

Hummel, J. R.↗

Effect of Accurate Modeling of Converter Interfaced Generation (CIG) on a Practical Bulk Power System

With increased penetration of converter interfaced generation (CIG) in power systems, it becomes imperative to adequately model them in both steady-state and dynamic simulations. There are several modeling approaches and simulation platforms for CIG. In addition, it is noted that suitability of modeling approaches and models is also dependent on the system conditions. However, their suitability has yet to be systematically analyzed, and this becomes a challenge in the industry without a systemic method. To tackle these challenges, in this paper, we first develop a methodology for identifying the maximum penetration of aggregated rooftop sources that can be represented as negative loads in both the steady state and dynamic analysis. Then we verify the suitability of positive-sequence performance-based model and detailed electromagnetic transient (EMT) model for utility-scale PV plants under different fault conditions. Based on the results, we perform dynamic simulations on a large practical system with both positive-sequence transient stability(TS) and EMT-TS hybrid simulation. Our study shows that positive-sequence performance based PV plant model is more suitable for study cases with faults relatively far away from the point of connection, and that positive-sequence models and TS simulation approach can produce comparable results with reference to EMT-TS hybrid simulation for normally-cleared single-line-to-ground faults, but the inadequacy of positive-sequence based model and simulation approach becomes obvious for prolonged-tripping of unbalanced faults in proximity of the PV plants.

Huang, Qiuhua↗

Comparative Analysis of Machine Learning Models for Day-Ahead Photovoltaic Power Production Forecasting

A main challenge for integrating the intermittent photovoltaic (PV) power generation remains the accuracy of day-ahead forecasts and the establishment of robust performing methods. The purpose of this work is to address these technological challenges by evaluating the day-ahead PV production forecasting performance of different machine learning models under different supervised learning regimes and minimal input features. Specifically, the day-ahead forecasting capability of Bayesian neural network (BNN), support vector regression (SVR), and regression tree (RT) models was investigated by employing the same dataset for training and performance verification, thus enabling a valid comparison. The training regime analysis demonstrated that the performance of the investigated models was strongly dependent on the timeframe of the train set, training data sequence, and application of irradiance condition filters. Furthermore, accurate results were obtained utilizing only the measured power output and other calculated parameters for training. Consequently, useful information is provided for establishing a robust day-ahead forecasting methodology that utilizes calculated input parameters and an optimal supervised learning approach. Finally, the obtained results demonstrated that the optimally constructed BNN outperformed all other machine learning models achieving forecasting accuracies lower than 5%.

14 SOLAR ENERGY↗

Computation of Sound Generated by Viscous Flow Over a Circular Cylinder

The Lighthill acoustic analogy approach combined with Reynolds-averaged Navier Stokes is used to predict the sound generated by unsteady viscous flow past a circular cylinder assuming a correlation length of 10 cylinder diameters. The two-dimensional unsteady flow field is computed using two Navier-Stokes codes at a low Mach number over a range of Reynolds numbers from 100 to 5 million. Both laminar flow as well as turbulent flow with a variety of eddy viscosity turbulence models are employed. Mean drag and Strouhal number are examined, and trends similar to experiments are observed. Computing the noise within the Reynolds number regime where transition to turbulence occurs near the separation point is problematic: laminar flow exhibits chaotic behavior and turbulent flow exhibits strong dependence on the turbulence model employed. Comparisons of far-field noise with experiment at a Reynolds number of 90,000, therefore, vary significantly, depending on the turbulence model. At a high Reynolds number outside this regime, three different turbulence models yield self-consistent results.

Cox, Jared S.↗

Computation of Vortex Shedding and Radiated Sound for a Circular Cylinder

The Lighthill acoustic analogy approach combined with Reynolds-averaged Navier Stokes is used to predict the sound generated by unsteady viscous flow past a circular cylinder assuming a correlation length of ten cylinder diameters. The two- dimensional unsteady ow field is computed using two Navier-Stokes codes at a low Mach number over a range of Reynolds numbers from 100 to 5 million. Both laminar ow as well as turbulent ow with a variety of eddy viscosity turbulence models are employed. Mean drag and Strouhal number are examined, and trends similar to experiments are observed. Computing the noise within the Reynolds number regime where transition to turbulence occurs near the separation point is problematic: laminar flow exhibits chaotic behavior and turbulent ow exhibits strong dependence on the turbulence model employed. Comparisons of far-field noise with experiment at a Reynolds number of 90,000, therefore, vary significantly, depending on the turbulence model. At a high Reynolds number outside this regime, three different turbulence models yield self-consistent results.

Cox, Jared S.↗

What Is the Right Level of Activation of a High-Spin {FeNO} 7 Complex to Enable Direct N–N Coupling? Mechanistic Insight into Flavodiiron NO Reductases

Flavodiiron nitric oxide reductases (FNORs), found in pathogenic bacteria, are capable of reducing nitric oxide (NO) to nitrous oxide (N 2 O) to detoxify NO released by the human immune system. Previously, we reported the first FNOR model system that mediates direct NO reduction (Dong, H. T.; et al. J. Am. Chem. Soc. 2018, 140, 13429-13440), but no intermediate of the reaction could be characterized. Here, we present a new set of model complexes that, depending on the ligand substitution, can either mediate direct NO reduction or stabilize a highly activated high-spin (hs) {FeNO} 7 complex, the first intermediate of the reaction. The precursors, [{Fe II (MPA-(RPhO) 2 )} 2 ] (1, R = H and 2, R = t Bu, Me), were prepared first and fully characterized. Complex 1 (without steric protection) directly reduces NO to N 2 O almost quantitatively, which constitutes only the second example of this reaction in model systems. Contrarily, the reaction of sterically protected 2 with NO forms the stable mononitrosyl complex 3, which shows one of the lowest N-O stretching frequencies (1689 cm -1 ) observed so far for a mononuclear hs-{FeNO} 7 complex. Here this study confirms that an N-O stretch & LE;1700 cm -1 represents the appropriate level of activation of the FeNO unit to enable direct NO reduction. The higher activation level of these hs-{FeNO} 7 complexes required for NO reduction compared to those formed in FNORs emphasizes the importance of hydrogen bonding residues in the active sites of FNORs to activate the bound NO ligands for direct N-N coupling and N2O formation. The implications of these results for FNORs are further discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗