Search NASA⌕ Search

SEARCH · Search NASA

Results for “Grid Simulation”

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,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Near-Wall Models in Large Eddy Simulations of Flow Behind a Backward-Facing Step

Accurate large eddy simulation (LES) of a wall-bounded flow generally requires a near-wall resolution comparable to that in direct numerical simulation (DNS). As much as 50% of the total grid points and computational costs are expended in the near-wall regions in a typical simulation. This limits LES to fairly low Reynolds numbers on current computers. To perform practical flow applications at realistically high Reynolds numbers, such as flow over an airfoil, it is desirable to replace very thin, near-wall regions in the LES with easily and inexpensively computed wall models to specify the near-wall boundary conditions.

Cabot, W.↗
DOE OSTI2016

Extended-Term Dynamic Simulations with High Penetrations of Photovoltaic Generation.

The uncontrolled intermittent availability of renewable energy sources makes integration of such devices into today's grid a challenge. Thus, it is imperative that dynamic simulation tools used to analyze power system performance are able to support systems with high amounts of photovoltaic (PV) generation. Additionally, simulation durations expanding beyond minutes into hours must be supported. This report aims to identify the path forward for dynamic simulation tools to accom- modate these needs by characterizing the properties of power systems (with high PV penetration), analyzing how these properties affect dynamic simulation software, and offering solutions for po- tential problems. We present a study of fixed time step, explicit numerical integration schemes that may be more suitable for these goals, based on identified requirements for simulating high PV penetration systems. We also present the alternative of variable time step integration. To help determine the characteristics of systems with high PV generation, we performed small signal sta- bility studies and time domain simulations of two representative systems. Along with feedback from stakeholders and vendors, we identify the current gaps in power system modeling including fast and slow dynamics and propose a new simulation framework to improve our ability to model and simulate longer-term dynamics.

14 SOLAR ENERGY↗
DOE OSTI2013

An Efficient Algorithm for Mapping Imaging Data to 3D Unstructured Grids in Computational Biomechanics

Geometries for organ scale and multiscale simulations of organ function are now routinely derived from imaging data. However, medical images may also contain spatially heterogeneous information other than geometry that are relevant to such simulations either as initial conditions or in the form of model parameters. In this manuscript, we present an algorithm for the efficient and robust mapping of such data to imaging based unstructured polyhedral grids in parallel. We then illustrate the application of our mapping algorithm to three different mapping problems: 1) the mapping of MRI diffusion tensor data to an unstuctured ventricular grid; 2) the mapping of serial cyro-section histology data to an unstructured mouse brain grid; and 3) the mapping of CT-derived volumetric strain data to an unstructured multiscale lung grid. Execution times and parallel performance are reported for each case.

Einstein, Daniel R.↗
DOE OSTI2016

Three-dimensional simulations of low foot and high foot implosion experiments on the National Ignition Facility

In order to achieve the several hundred Gbar stagnation pressures necessary for inertial confinement fusion ignition, implosion experiments on the National Ignition Facility (NIF) [E. I. Moses et al., Phys. Plasmas 16, 041006 (2009)] require the compression of deuterium-tritium fuel layers by a convergence ratio as high as forty. Such high convergence implosions are subject to degradation by a range of perturbations, including the growth of small-scale defects due to hydrodynamic instabilities, as well as longer scale modulations due to radiation flux asymmetries in the enclosing hohlraum. Due to the broad range of scales involved, and also the genuinely three-dimensional (3D) character of the flow, accurately modeling NIF implosions remains at the edge of current simulation capabilities. This paper describes the current state of progress of 3D capsule-only simulations of NIF implosions aimed at accurately describing the performance of specific NIF experiments. Current simulations include the effects of hohlraum radiation asymmetries, capsule surface defects, the capsule support tent and fill tube, and use a grid resolution shown to be converged in companion two-dimensional simulations. The results of detailed simulations of low foot implosions from the National Ignition Campaign are contrasted against results for more recent high foot implosions. While the simulations suggest that low foot performance was dominated by ablation front instability growth, especially the defect seeded by the capsule support tent, high foot implosions appear to be dominated by hohlraum flux asymmetries, although the support tent still plays a significant role. For both implosion types, the simulations show reasonable, though not perfect, agreement with the data and suggest that a reliable predictive capability is developing to guide future implosions toward ignition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2017

Model Validation for Propulsion - On the TFNS and LES Subgrid Models for a Bluff Body Stabilized Flame

This paper summarizes the reacting results of simulating a bluff body stabilized flame experiment of Volvo Validation Rig using a releasable edition of the National Combustion Code (NCC). The turbulence models selected to investigate the configuration are the sub-grid scaled kinetic energy coupled large eddy simulation (K-LES) and the time-filtered Navier-Stokes (TFNS) simulation. The turbulence chemistry interaction used is linear eddy mixing (LEM).

Flame↗
DOE OSTI2022

Large-Eddy Simulation of a wind turbine using a Filtered Actuator Line Model

When dealing with multirotor devices such as quadcopters or wind farms, the cost of blade-resolved large-eddy simulation (LES) becomes prohibitive. Combining LES with a family of lower-fidelity models, called actuator line models (ALMs), has grown in popularity in the past decade. ALM replaces full blade resolution with an array of actuator points or lines parameterized by aerodynamic lift/drag polar plots along the blades. Body forces computed based on these actuator points are then projected onto the LES flow mesh, mimicking the effect of rotating blades on the flow. However, the optimal projection radius and the associated LES grid size is often too restrictive for multirotor simulations. Recently, a new tip-correction-based filtered ALM (F-ALM) was proposed by Martinez-Tossas and Meneveau (2019), which allows coarser-than-optimal grids by avoiding the associated overprediction of thrust. In this work, F-ALM is implemented into a high-order, in-house LES code to simulate National Renewable Energy Laboratory Phase VI wind turbine. It is then followed by a comparison between the baseline ALM and the newly implemented F-ALM in terms of instantaneous and time-averaged flow fields and blade loads, revealing the advantages of F-ALM in preventing the overprediction of power on coarse grids. Finally, this encourages accurate and affordable simulations of multirotor devices in the future.

17 WIND ENERGY↗
DOE OSTI2024

Simulation study of ion beam used to produce Mo-99

An 820 mA CW positive ion source is being developed to produce Mo-99 using the fusion of deuterium and tritium ion beams on a rotating target to produce neutrons for use in the production of radiopharmaceuticals. The ion source consists of an RF plasma source, a multi-aperture extractor, and 300 kV accelerating column. This paper will describe a simulation study of the beam through the extractor grid and the accelerator to the target. The uniformity of beam distribution on the target is an important aspect of the simulation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

A Novel State-Constrained Primary Control for Grid-Forming Inverters

The problem of operating an AC microgrid with purely inverter-based generating units is considered, and a novel vector control for the constrained grid-forming (GFM) inverters is proposed. The design aims to achieve the primary control objectives of the GFM inverters: voltage establishment, synchronization, and voltage/frequency regulation. It is novel that the proposed control also ensures all operational variables constrained both real-time and for all the time except for low voltage during the period of grid forming. The constraints include frequency, output AC voltage, input DC power, and active/reactive power injections to the grid. Here, in this paper, outcomes of a detailed simulation are presented to demonstrated that the above design objectives are met and that the overall performance is superior due to its unique ability of simultaneously controlling the whole vectors of the inverter state and output.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2015

HYDRODYNAMIC SIMULATIONS OF H ENTRAINMENT AT THE TOP OF He-SHELL FLASH CONVECTION

We present the first three-dimensional, fully compressible gas-dynamics simulations in 4π geometry of He-shell flash convection with proton-rich fuel entrainment at the upper boundary. This work is motivated by the insufficiently understood observed consequences of the H-ingestion flash in post-asymptotic giant branch (post-AGB) stars (Sakurai's object) and metal-poor AGB stars. Our investigation is focused on the entrainment process at the top convection boundary and on the subsequent advection of H-rich material into deeper layers, and we therefore ignore the burning of the proton-rich fuel in this study. We find that for our deep convection zone, coherent convective motions of near global scale appear to dominate the flow. At the top boundary convective shear flows are stable against Kelvin-Helmholtz instabilities. However, such shear instabilities are induced by the boundary-layer separation in large-scale, opposing flows. This links the global nature of thick shell convection with the entrainment process. We establish the quantitative dependence of the entrainment rate on grid resolution. With our numerical technique, simulations with 1024{sup 3} cells or more are required to reach a numerical fidelity appropriate for this problem. However, only the result from the 1536{sup 3} simulation provides a clear indication that we approach convergence with regard to the entrainment rate. Our results demonstrate that our method, which is described in detail, can provide quantitative results related to entrainment and convective boundary mixing in deep stellar interior environments with very stiff convective boundaries. For the representative case we study in detail, we find an entrainment rate of 4.38 ± 1.48 × 10{sup –13} M {sub ☉} s{sup –1}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Modeling Framework for Bulk Electric Grid Impacts from HEMP E1 and E3 Effects (Tasks 3.1 Final Report)

This report presents a framework to evaluate the impact of a high-altitude electromagnetic pulse (HEMP) event on a bulk electric power grid. This report limits itself to modeling the impact of EMP E1 and E3 components. The co-simulation of E1 and E3 is presented in detail, and the focus of the paper is on the framework rather than actual results. This approach is highly conservative as E1 and E3 are not maximized with the same event characteristics and may only slightly overlap. The actual results shown in this report are based on a synthetic grid with synthetic data and a limited exemplary EMP model. The framework presented can be leveraged and used to analyze the impact of other threat scenarios, both manmade and natural disasters. This report d escribes a Monte-Carlo based methodology to probabilistically quantify the transient response of the power grid to a HEMP event. The approach uses multiple fundamental steps to characterize the system response to HEMP events, focused on the E1 and E3 components of the event. 1) Obtain component failure data related to HEMP events testing of components and creating component failure models. Use the component failure model to create component failure conditional probability density function (PDF) that is a function of the HEMP induced terminal voltage. 2) Model HEMP scenarios and calculate the E1 coupled voltage profiles seen by all system components. Model the same HEMP scenarios and calculate the transformer reactive power consumption profiles due to E3. 3) Sample each component failure PDF to determine which grid components will fail, due to the E1 voltage spike, for each scenario. 4) Perform dynamic simulations that incorporate the predicted component failures from E1 and reactive power consumption at each transformer affected by E3. These simulations allow for secondary transients to affect the relays/protection remaining in service which can lead to cascading outages. 5) Identify the locations and amount of load lost for each scenario through grid dynamic simulation. This can be an indication of the immediate grid impacts from a HEMP event. In addition, perform more detailed analysis to determine critical nodes and system trends. 6) To help realize the longer-term impacts, a security constrained alternating current optimal power flow (ACOPF) is run to maximize critical load served. This report describes a modeling framework to assess the systemic grid impacts due to a HEMP event. This stochastic simulation framework generates a large amount of data for each Monte Carlo replication, including HEMP location and characteristics, relay and component failures, E3 GIC profiles, cascading dynamics including voltage and frequency over time, and final system state. This data can then be analyzed to identify trends, e.g., unique system behavior modes or critical components whose failure is more likely to cause serious systemic effects. The proposed analysis process is demonstrated on a representative system. In order to draw realistic conclusions of the impact of a HEMP event on the grid, a significant amount of work remains with respect to modeling the impact on various grid components.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

FY 2021 Isolated Grids and Grid-Connected Turbine Reference Systems; Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL)

For individuals, businesses, and communities focused on building resilient electrical grid infrastructure, wind energy can provide an affordable, accessible, and compatible distributed energy resource option that also enhances the capabilities of local grid operations. The Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL) project is a multi-year distributed wind research effort, driven through a partnership between four Department of Energy National Laboratories and industry to develop and improve the planning, design, and operation of wind-centered microgrids to complement solar, energy storage, and other distributed energy resources for grid-tied and isolated operation (U.S. Department of Energy, 2021). This report documents the application of methods developed through the initial three years of the MIRACL project to two real-world distributed wind reference systems. Specifically, the methods demonstrated in this report include 1) a market valuation framework to comprehensively value the services distributed wind can provide and 2) a resilience framework that enables stakeholders to characterize distribution system resilience and compare grid investment decisions from a resilience perspective. Additional methods mentioned in this report include distributed hybrid system design methods for grid resilience, advanced control co-simulation platforms, and power hardware-in-the-loop (PHIL) models. Preliminary results from these additional methods are presented in this report and will be demonstrated and/or applied to the reference systems in the coming year. The purpose of applying these methods to reference systems is to drive technology transfer of the theories, methodologies, and technologies developed under the MIRACL project and increase the number of referenceable case studies available to stakeholders interested in additional value-added capabilities of wind systems beyond bulk energy supply (i.e. kilowatt-hours).

17 WIND ENERGY↗
NASA NTRS1994

High performance real-time flight simulation at NASA Langley

In order to meet the stringent time-critical requirements for real-time man-in-the-loop flight simulation, computer processing operations must be deterministic and be completed in as short a time as possible. This includes simulation mathematical model computational and data input/output to the simulators. In 1986, in response to increased demands for flight simulation performance, personnel at NASA's Langley Research Center (LaRC), working with the contractor, developed extensions to a standard input/output system to provide for high bandwidth, low latency data acquisition and distribution. The Computer Automated Measurement and Control technology (IEEE standard 595) was extended to meet the performance requirements for real-time simulation. This technology extension increased the effective bandwidth by a factor of ten and increased the performance of modules necessary for simulator communications. This technology is being used by more than 80 leading technological developers in the United States, Canada, and Europe. Included among the commercial applications of this technology are nuclear process control, power grid analysis, process monitoring, real-time simulation, and radar data acquisition. Personnel at LaRC have completed the development of the use of supercomputers for simulation mathematical model computational to support real-time flight simulation. This includes the development of a real-time operating system and the development of specialized software and hardware for the CAMAC simulator network. This work, coupled with the use of an open systems software architecture, has advanced the state of the art in real time flight simulation. The data acquisition technology innovation and experience with recent developments in this technology are described.

Cleveland, Jeff I., II↗
DOE OSTI2021

Implementation of a 9-point stencil in SOLPS-ITER and implications for Alcator C-Mod divertor plasma simulations

The SOLPS-ITER code suite is used worldwide for plasma edge modeling, the interpretation of experiments, as well as for the design of the ITER divertor. The numerical scheme of the plasma solver of the code, B2.5, is based on the assumption of perfectly field-aligned grids, while in practice grids are often strongly distorted to match divertor target shapes. Neglecting these grid distortion leads to qualitatively and quantitatively incorrect results for fluid neutral simulations, and may affect results in cold (detached) divertors even when using kinetic neutral simulations. In this contribution, we present the first results of a newly implemented 9-point stencil in B2.5 to properly handle misaligned grids. The new scheme is then applied to fluid neutral simulations of a well-diagnosed and previously modeled Alcator C-Mod discharge. Results are compared with the original 5-point scheme neglecting grid distortion effects, as well as with simulations including a full kinetic neutral model. We conclude that the 9-point stencil is essential to correctly model the transport of fluid neutrals on distorted grids, and to capture the effects of divertor closure on the fluid neutral behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Real-Time Testbed for Studying Cyberattacks and Defense in DER-integrated Smart Inverter Systems

In this paper, we propose a Hardware-in-the-Loop (HIL) simulation testbed suitable for the implementation and testing of realistic cyberattacks on grid-tied smart inverter systems integrated with Distributed Energy Resources (DER) that use the Distributed Network Protocol-3 (DNP3) protocol for communications between grid components. Specifically, our testbed combines a Real-Time Digital Simulator (RTDS) NovaCor device, outfitted with GNETx2 network interface cards, a gridtied DER topology implemented via the RTDS software package RSCAD, and a custom virtual network that emulates a man in the middle attacker. The Man-in-the-Middle (MITM) attacker captures DNP3 traffic and falsifies telemetry data in DNP3 packets to trigger unwarranted commands from a DNP3 controller that exploit smart inverter grid support functions. We choose DNP3 and implement grid support functions according to the IEEE Std. 1547-2018 mandated for the interconnection and interoperability of DER power systems with associated power components. Furthermore, we develop a protocol payload agnostic attack detection framework that leverages the round-trip time (RTT) anomalies between DNP3 requests and responses and can detect the presence of attacks without having to analyze the payload’s contents, while balancing trade-offs between false alarm counts, missed detections, and time to detection. To facilitate further research, we publicly release benign and attack network traffic exchanged between various sensors, controllers, and actuators in our grid-tied inverter testbed.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRSDate not supplied

A Large-Eddy Simulation Of Turbulent Flow Over A Backward Facing Step

Large-eddy simulation (LES) was used to study the turbulent flow over a backward facing step. The simulations examine the details of the turbulent flow's separated and reattaching regions created by the step geometry. An explicit high-order of accuracy/high resolution finite-difference computational fluid dynamics solver was used for both LES and Reynolds-averaged Navier-Stokes (RANS) simulations. The experimental data from the backward facing step experiment of Driver and Seegmiller was used for comparison. Data presented are derived from a RANS simulation using the Spalart-Allmaras turbulence model and LES simulations using the implict LES assumption on three different grids of varying resolution. The RANS solution agrees with past work. It does an adequate job of predicting the flowfield and reattachment point, but misses details such as the strength of the recirculating region and magnitudes of the Reynolds stresses. The synthetic eddy method was used to provide a turbulent inflow to the quasi-2D LES domains. For the LES, the coarse grid was of insufficient resolution and under-predicted the skin friction coefficient and over-predicted the normal stresses. The finer grid LES results accurately predict the skin friction coefficient and provide more accurate Reynolds stresses. Triple product velocity correlations are also computed from the LES and the finer grid results did a good job predicting the magnitudes of the curves and the shapes were reasonably represented. It is recommended that future work include investigating the turbulent inflow domain length, width of the quasi-2D domain, and additional grid refinement.

turbulence↗
DOE OSTI2013

Transforming the representation of the boundary layer and low clouds for high-resolution regional climate modeling: Final report

Stratocumulus and shallow cumulus clouds in subtropical oceanic regions (e.g., Southeast Pacific) cover thousands of square kilometers and play a key role in regulating global climate (e.g., Klein and Hartmann, 1993). Numerical modeling is an essential tool to study these clouds in regional and global systems, but the current generation of climate and weather models has difficulties in representing them in a realistic way (e.g., Siebesma et al., 2004; Stevens et al., 2007; Teixeira et al., 2011). While numerical models resolve the large-scale flow, subgrid-scale parameterizations are needed to estimate small-scale properties (e.g. boundary layer turbulence and convection, clouds, radiation), which have significant influence on the resolved scale due to the complex nonlinear nature of the atmosphere. To represent the contribution of these fine-scale processes to the resolved scale, climate models use various parameterizations, which are the main pieces in the model that contribute to the low clouds dynamics and therefore are the major sources of errors or approximations in their representation. In this project, we aim to 1) improve our understanding of the physical processes in thermal circulation and cloud formation, 2) examine the performance and sensitivity of various parameterizations in the regional weather model (Weather Research and Forecasting model; WRF), and 3) develop, implement, and evaluate the advanced boundary layer parameterization in the regional model to better represent stratocumulus, shallow cumulus, and their transition. Thus, this project includes three major corresponding studies. We find that the mean diurnal cycle is sensitive to model domain in ways that reveal the existence of different contributions originating from the Southeast Pacific land-masses. The experiments suggest that diurnal variations in circulations and thermal structures over this region are influenced by convection over the Peruvian sector of the Andes cordillera, while the mostly dry mountain-breeze circulations force an additional component that results in semi-diurnal variations near the coast. A series of numerical tests, however, reveal sensitivity of the simulations to the choice of vertical grid, limiting the possibility of solid quantitative statements on the amplitudes and phases of the diurnal and semidiurnal components across the domain. According to our experiments, the Mellor-Yamada-Nakanishi-Niino (MYNN) boundary layer scheme and the WSM6 microphysics scheme is the combination of schemes that performs best. For that combination, mean cloud cover, liquid water path, and cloud depth are fairly wellsimulated, while mean cloud top height remains too low in comparison to observations. Both microphysics and boundary layer schemes contribute to the spread in liquid water path and cloud depth, although the microphysics contribution is slightly more prominent. Boundary layer schemes are the primary contributors to cloud top height, degree of adiabaticity, and cloud cover. Cloud top height is closely related to surface fluxes and boundary layer structure. Thus, our study infers that an appropriate tuning of cloud top height would likely improve the low-cloud representation in the model. Finally, we show that entrainment governs the degree of adiabaticity, while boundary layer decoupling is a control on cloud cover. In the intercomparison study using WRF single-column model experiments, most parameterizations show a poor agreement of the vertical boundary layer structure when compared with large-eddy simulation models. We also implement a new Total-Energy/Mass- Flux boundary layer scheme into the WRF model and evaluate its ability to simulate both stratocumulus and shallow cumulus clouds. Result comparisons against large-eddy simulation show that this advanced parameterization based on the new Eddy-Diffusivity/Mass-Flux approach provides a better performance than other boundary layer parameterizations.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2011

Effects of Spatial Resolution on the Simulated Dust Aerosol Lifecycle: Implications for Dust Event Magnitude and Timing in the NASA GEOS-5 AGCM

The NASA GEOS-5 atmospheric transport model simulates global aerosol distributions with an online aerosol module. GEOS-5 may be run at various horizontal spatial resolutions depending on the research application. For example, long integration climate simulations are typically run at 2 deg or 1 deg grid spacing, whereas aerosol reanalysis and forecasting applications may be performed at O.5 deg or 0.25 deg resolutions. In this study, we assess the implications of varying spatial resolution on the simulated aerosol fields, with a particular focus on dust. Dust emissions in GEOS-5 are calculated with one of two parameterizations, one based on the Goddard Chemistry, Aerosol, Radiation, and Transport (GO CART) model and another based on the Dust Entrainment and Deposition (DEAD) model. Emission fluxes are parameterized in terms of the surface wind speed, either the 10-m (GO CART) or friction (DEAD) wind speed. We consider how surface wind speeds and thus the dust emission rates are a function of the model spatial resolution. We find that spatial resolution has a significant effect on the magnitude of dust emissions, as higher resolution versions of the model have typically higher surface wind speeds. Utilizing space-borne observations from MISR, MODIS, and CALIOP, we find that simulated Aerosol Optical Thickness (AOT) distributions respond differently to spatial resolution over the African and Asian source regions, highlighting the need to regional dust emission tuning. When compared to ground-based observations from AERONET, we found improved timing of dust events with as spatial resolution was increased. In an attempt to improve the representation of the dust aerosol lifecycle at coarse resolutions, we found that incorporating the effects of sub-grid wind variability in a course resolution simulation led to improved agreement with observed AOT magnitudes, but did not impact the timing of simulated dust events.

Nowottnick, E.↗
DOE OSTI2019

Automated Multiscale Approach To Predict Self-Diffusion from a Potential Energy Field

For large-scale screening studies there is a need to estimate the diffusion of gas molecules in nanoporous materials more efficiently than (brute force) molecular dynamics. In particular for systems with low diffusion coefficients molecular dynamics can be prohibitively expensive. An alternative is to compute the hopping rates between adsorption sites using transition state theory. For large-scale screening this requires the automatic detection of the transition states between the adsorption sites along the different diffusion paths. Here an algorithm is presented that analyzes energy grids for the moving particles. It detects the energies at which diffusion paths are formed, together with their directions. This allows for easy identification of nondiffusive systems. For diffusive systems, it partitions the grid coordinates assigned to energy basins and transitions states, permitting a transition state theory based analysis of the diffusion. We test our method on CH 4 diffusion in zeolites, using a standard kinetic Monte Carlo simulation based on the output of our grid analysis. We find that it is accurate, fast, and rigorous without limitations to the geometries of the diffusion tunnels or transition states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗