Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase Field Model”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Microstructure-Level Modeling of Stage 3 Fission Gas Release in UO 2 Fuel

A new phase-field model of fission gas bubble evolution was developed and applied to gain an improved understanding of the microstructure-level processes leading to fission gas release from nuclear fuel and to inform improvements to the fission gas release model used in BISON. The phase-field model simulates the evolution of multiple fuel grains and intergranular bubbles simultaneously, tracks the local concentration of vacancies and gas atoms, and includes the effect of hydrostatic gas pressure on the surrounding fuel matrix. This represents the most realistic microstructure-level model of fission gas bubbles developed to date. The model was used to simulate the growth of grain face and grain edge bubbles in a hexagonally periodic 3D grain structure, which contains multiple grain boundaries and triple junctions. The coverage of grain boundaries and triple junctions by fission gas bubbles was calculated as a function of time. Preliminary results suggest that grain edge bubbles do not release a significant amount of gas prior to gas release by grain face bubbles, which supports the assumptions of the fission gas release model currently used in BISON. However, the simulations should be improved by including a physically realistic number of grain edge bubbles in the initial conditions before definitive conclusions can be drawn.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2012

Progress on Generic Phase-Field Method Development

In this report, we summarize our current collaborative efforts, involving three national laboratories: Idaho National Laboratory (INL), Pacific Northwest National Laboratory (PNNL) and Los Alamos National Laboratory (LANL), to develop a computational framework for homogenous and heterogenous nucleation mechanisms into the generic phase-field model. During the studies, the Fe-Cr system was chosen as a model system due to its simplicity and availability of reliable thermodynamic and kinetic data, as well as the range of applications of low-chromium ferritic steels in nuclear reactors. For homogenous nucleation, the relevant parameters determined from atomistic studies were used directly to determine the energy functional and parameters in the phase-field model. Interfacial energy, critical nucleus size, nucleation rate, and coarsening kinetics were systematically examined in two- and three- dimensional models. For the heterogeneous nucleation mechanism, we studied the nucleation and growth behavior of chromium precipitates due to the presence of dislocations. The results demonstrate that both nucleation schemes can be introduced to a phase-field modeling algorithm with the desired accuracy and computational efficiency.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Lack of robustness and accuracy of many numerical schemes for phase-field simulations

In this paper, we study the stability, accuracy and convergence behavior of various numerical schemes for phase-field modeling through a simple ODE model. Both theoretical analysis and numerical experiments are carried out on this ODE model to demonstrate the limitation of most numerical schemes that have been used in practice. One main conclusion is that the first-order fully implicit scheme is the only robust algorithm for phase-field simulations while all other schemes (that have been analyzed) may have convergence issue if the time step size is not exceedingly small. More specifically, by rigorous analysis in most cases, we have the following conclusions: (i) The first-order fully implicit scheme converges to the correct steady state solution for all time step sizes. In the case of multiple solutions, one of the solution branches always converges to the correct steady state solution. (ii) The first-order convex splitting scheme, which is equivalent to the first-order fully implicit scheme with a different time scaling, always converges to the correct steady state solution but may seriously lack numerical accuracy for transient solutions. (iii) For the second-order fully implicit and convex splitting schemes, for any time step size $δt > 0$, there exists an initial condition $u_0$ , with $|u_0| > 1$, such that the numerical solution converges to the wrong steady state solution. (iv) For $|u_0|$ ≤ 1, all second-order schemes studied in this paper converge to the correct steady state solution although severe numerical oscillations occur for most of them if the time step size is not sufficiently small. (v) An unconditionally energy-stable scheme (such as the modified Crank–Nicolson scheme) is not necessarily better than a conditionally energy-stable scheme (such as the Crank–Nicolson scheme). Finally, most, if not all, of the above conclusions are expected to be true for more general Allen–Cahn and other phase-field models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Microstructure development during rapid alloy solidification

Abstract Solidification processing of structural alloys can take place over an extremely wide range of solid–liquid interface velocities spanning six orders of magnitude, from the low-velocity constitutional supercooling limit of microns/s to the high-velocity absolute stability limit of m/s. In between these two limits, the solid–liquid interface is morphologically unstable and typically forms cellular-dendritic microstructures, but also other microstructures that remain elusive. Rapid developments in additive manufacturing have renewed the interest in modeling the high-velocity range, where approximate analytical theories provide limited predictions. In this article, we discuss recent advances in phase-field modeling of rapid solidification of metallic alloys, including a brief description of state-of-the-art experiments used for model validation. We describe how phase-field models can cope with the dual challenge of carrying out simulations on experimentally relevant length- and time scales and incorporating nonequilibrium effects at the solid–liquid interface that become dominant at rapid rates. We present selected results, illustrating how phase-field simulations have yielded unprecedented insights into high-velocity interface dynamics, shedding new light on both the absolute stability limit and the formation of banded microstructures that are a hallmark of rapid alloy solidification near this limit. We also discuss state-of-the-art experiments used to validate those insights. Graphical abstract

36 MATERIALS SCIENCE↗
DOE OSTI2012

Predicting Thermal Conductivity Evolution of Polycrystalline Materials Under Irradiation Using Multiscale Approach

A multiscale methodology was developed to predict the evolution of thermal conductivity of polycrystalline fuel under irradiation. In the mesoscale level, phase field model was used to predict the evolution of gas bubble microstructure. Generation of gas atoms and vacancies were taken into consideration. In the macroscopic scale, a statistical continuum mechanics model was applied to predict the anisotropic thermal conductivity evolution during irradiation. Microstructure predicted by phase field model was fed into statistical continuum mechanics model to predict properties and behavior. Influence of irradiation intensity, exposition time and morphology were investigated. This approach provides a deep understanding on microstructure evolution and property prediction from a basic scientific viewpoint.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Preferential Paths of Air-water Two-phase Flow in Porous Structures with Special Consideration of Channel Thickness Effects

Accurate understanding and predicting the flow paths of immiscible two-phase flow in rocky porous structures are of critical importance for the evaluation of oil or gas recovery and prediction of rock slides caused by gas-liquid flow. A 2D phase field model was established for compressible air-water two-phase flow in heterogenous porous structures. The dynamic characteristics of air-water two-phase interface and preferential paths in porous structures were simulated. The factors affecting the path selection of two-phase flow in porous structures were analyzed. Transparent physical models of complex porous structures were prepared using 3D printing technology. Tracer dye was used to visually observe the flow characteristics and path selection in air-water two-phase displacement experiments. The experimental observations agree with the numerical results used to validate the accuracy of phase field model. The effects of channel thickness on the air-water two-phase flow behavior and paths in porous structures were also analyzed. The results indicate that thick channels can induce secondary air flow paths due to the increase in flow resistance; consequently, the flow distribution is different from that in narrow channels. This study provides a new reference for quantitatively analyzing multi-phase flow and predicting the preferential paths of immiscible fluids in porous structures.

42 ENGINEERING↗
DOE OSTI2023

Convection-Induced Compositional Patterning at Grain Boundaries in Irradiated Alloys

We consider the stability of precipitates formed at grain boundaries (GBs) by radiation-induced segregation in dilute alloys subjected to irradiation. Here, the effects of grain size and misorientation of symmetric-tilt GBs are quantified using phase field modeling. A novel regime is identified where, at long times, GBs are decorated by precipitate patterns that resist coarsening. Maps of the chemical Péclet number indicate that arrested coarsening takes place when solute advection dominates over thermal diffusion right up to the precipitate/matrix interface, preventing interfacial local equilibrium and overriding capillary effects. This contrasts with liquid-solid mixtures where convection always accelerates coarsening.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Model Development in MARMOT for High Burn-up Structure Formation

A phase field model for the formation of High Burn-up Structure (HBS) was developed and implemented in MARMOT. The model takes into consideration the chemical energy of gas atoms, interfacial energies of grain boundaries and bubble surfaces, and strain energy associated with dislocations, which renders the model capable of simulating the concurrent formation and growth of both bubbles and sub-grains. Moreover, in contrast to existing phase-field models of recrystallization, the current model makes no ad hoc or a priori assumptions for the nucleation rate or recrystallized grain size and shape. All the model parameters were quantified in terms of thermodynamic and kinetic properties. The model predicts a strong effect of magnitude and distribution of dislocation density, grain boundary energy, and bubble surface energy on the formation of sub-grains. The model results were shown to agree well with theoretical predictions. For polycrystalline UO 2 at 1200K, the recrystallized average grain size was found to be on the order of of 0.6-1.3 microns at a burn-up in the range of 98-105 GWD/tU. These values lie well within the range of reported data in literature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Development of Macroscale Models of UO 2 Fuel Sintering and Densification using Multiscale Modeling and Simulation

Light water reactor fuel pellets are fabricated using sintering to final densities of 95% or greater. During reactor operation, the porosity remaining in the fuel after fabrication decreases further due to irradiation-assisted densification. While empirical models have been developed to describe this densification process, a mechanistic model is needed as part of the ongoing work by the NEAMS program to develop a more predictive fuel performance code. In this work we will develop a phase field model of sintering of UO 2 in the MARMOT code, and validate it by comparing to published sintering data. We will then add the capability to capture irradiation effects into the model, and use it to develop a mechanistic model of densification that will go into the BISON code and add another essential piece to the microstructure-based materials models. The final step will be to add the effects of applied fields, to model field-assisted sintering of UO 2 . The results of the phase field model will be validated by comparing to data from field-assisted sintering. Tasks over three years: 1. Develop a sintering model for UO 2 in MARMOT 2. Expand model to account for irradiation effects 3. Develop a mechanistic macroscale model of densification for BISON

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Effect of patterned inclusions on the fracture behavior of ceramic composites

Developments in additive manufacturing have enabled the reinforcement of ceramic composites across the scales. Understanding the fracture behavior of such composites is critical for structural design and engineering. However, systematic investigations on the relationship between the patterning characteristics and the effective fracture properties have been scarce. In this work, we propose a mixed-mode phase field model to investigate the effects of elastic modulus and pattern characteristics of circular inclusions on the fracture behavior of a ceramic composite beam. The modeling results have shown that the stiff and compliant inclusions can either repel or attract the crack. The extent of repelling and attraction can be controlled by altering the pattern characteristics. In addition, utilizing the relationship between the pattern characteristics and resulting fracture path, fracture growth can be controlled with appropriate inclusion properties. The targeted enhancement in fracture toughness can also be predicted using an effective fracture model linked with the underlying fracture path computed by the proposed mixed-mode phase field model.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Reduced-order model for microstructure evolution prediction in the electrodes of solid oxide fuel cell with dynamic discrepancy reduced modeling

Microstructure evolution in the electrodes of solid oxide fuel cell is an important degradation mechanism which reduces active sites for redox reaction and the electric conductivity. Phase field models for microstructure evolution simulation are usually expensive for large scale simulations. In this paper, a reduced-order coarsening model is developed using dynamic discrepancy reduced modeling, which reduces the model order by inserting Gaussian process stochastic functions into the dynamic equations of Ostwald ripening. The reduced order model has been calibrated on a dataset generated by a phase field model that has been well validated to experiments. A validating dataset has also been generated with which the model prediction show good agreement. This model is further applied to predict long term microstructure evolution in different SOFC electrodes. Finally, this work is the first attempt of building a degradation model of SOFC using data science techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Implementation and testing of physics-based pulverization model in BISON

This report summarizes lower length scale computational research conducted to improve the pulverization criterion for high-burnup UO 2 fuel in the BISON fuel performance code. This research was sponsored by the NEAMS program during FY22. Efforts to improve the model primarily focus on calculating the current pressure of bubbles in the high-burnup structure (HBS) region, as well as the critical bubble pressure at which pulverization occurs in the HBS. The phase-field model for predicting initial bubble pressure in the HBS was improved by implementing a more realistic model for defect production, and by coupling the phase-field model for inter-granular bubble evolution with the spatially-resolved cluster dynamics code Xolotl to simulate intra-granular fission gas evolution. An evolution equation for current bubble pressure was added to BISON, an improvement over the previous model that assumed that bubble pressure was static following HBS formation. To improve calculations of critical bubble pressure for pulverization, 3-D phase-field fracture simulations were performed, and a function for critical bubble pressure was fit—based on the simulation results—to replace the previous function, which had been determined using 2-D simulations. The impacts of these modifications to the existing pulverization assessment cases will be reported in the forthcoming engineering scale milestone on high-burnup UO 2 pulverization. A poromechanics-based approach was used to include the effect of bubble over pressurization on the stress state of the pellet at the engineering scale, and the initial strategy for integration of the pulverization criterion with pellet-scale mechanical degradation using a smeared cracking model was developed. This report also describes the initial research that was conducted to inform a forthcoming transient fission gas release (tFGR) model in BISON. Focusing on the fission gas release caused by pulverization of the outer pellet rim, a model for the amount of fission gas release was developed as a function of fuel porosity, bubble size, and fragment size. Implementation of this model in BISON and its impact of engineering-scale predictions will also be described in the aforementioned engineering-scale milestone on high-burnup UO 2 pulverization.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Three-dimensional imaging of dislocation dynamics during the hydriding phase transformation

Crystallographic imperfections can significantly alter material properties and responses to external stimuli, including solute induced phase transformations and crystal growth and dissolution . Despite recent progress in imaging defects using both electron and x-ray techniques, in situ three-dimensional imaging studies of defect dynamics, necessary to understand and engineer nanoscale processes, remains challenging. Here, we report in situ three-dimensional imaging of defect dynamics during the hydriding phase transformation of individual palladium nanocrystals by Bragg Coherent Diffractive Imaging (BCDI) . During constant pressure experiments, we observed that the phase transformation begins after the nucleation of dislocations in large (300 nm) particles. The 3D dislocation network shows that dislocations are close to the phase boundary. The 3D phase morphology resolved by BCDI suggests that the hydrogen-rich phase is more similar to a spherical cap on the hydrogen-poor phase than the core-shell model commonly assumed. We substantiate this conclusion using 3D phase field modeling and demonstrate how phase morphology affects the critical size for dislocation nucleation. We determine the size dependence of the transformation pressure for large (150-300 nm) palladium nanocrystals using variable pressure experiments. Our results reveal a pathway for solute induced structural phase transformations in nanocrystals and demonstrate BCDI as a novel method for understanding dislocation dynamics in phase transforming systems at the nanoscale.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mesoscale Modeling for Restructuring and Fragmentation in High Burnup UO 2

This report summarizes the mesoscale modeling conducted in fiscal year 2025 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, focusing on the microstructural evolution and restructuring in high burnup UO 2 nuclear fuel and its impact on fuel fragmentation. We developed a pioneering phase-field model to simulate restructuring behavior across different regions of high burnup fuel, including the dark zone and rim region. A grand-potential-based phase-field model is employed to concurrently evaluate subgrain formation and the growth of fission gas bubbles within the fuel. An energy-based subgrain formation criterion was introduced to simulate the restructuring process. The effects of temperature and burnup rate were studied to capture how each of these parameters influences the characteristics of the restructured fuel. Subgrain formation was observed to initiate around existing fission gas bubbles and proceed toward triple junctions, grain boundaries, and grain interiors. Under a given subgrain formation rate, the rate of restructuring increases with rising fuel temperature. The restructuring occurs faster with higher burnup rate. A restructuring bias was observed within the microstructure, due to the variation in defect accumulation when comparing different grains. Microstructures corresponding to the dark zone and rim region can be obtained by parameterizing the model with the appropriate defect production rate, as determined based on the burnup rate and temperature. The predicted microstructures are consistent with experimental observations of the restructured regions. Based on the mesoscale simulations, a mechanistic model for restructuring and grain size evolution was implemented in BISON. Thus, this work provides a first-of-its-kind restructuring model for different regions of high-burnup fuel to BISON. It was shown that the model predicts appropriate grain size evolution along fuel radius as those observed in experiments. This model enables BISON to predict the effect of restructuring on the fission gas release. Finally, the phase-field fracture simulation with dark zone specific microstructures were presented to provide the fragmentation criteria for the dark zone. This work introduces a first-of-its-kind restructuring model for high burnup fuel, significantly enhancing BISON's predictive capabilities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

A phase-field study of stainless-steel oxidation from high-temperature carbon dioxide exposure

An electrochemical phase-field model has been developed to investigate the oxidation mechanisms of the 21-2N valve stainless steel alloy exposed to carbon dioxide at 973 K. Three oxide phases observed in oxidation experiments are included in the model: Mn 3 O 4 , Cr 2 O 3 , and MnCr 2 O 4 . Local charge neutrality and the conserved current condition are assumed to include the impact of the electric potential on oxidation, though it is found to be negligible when the system is electrically isolated. The sensitivity of the oxidation processes to the values of the diffusion mobilities is examined, and the oxidation rate is calibrated against experimental data by modifying the sensitive mobilities. We find that both inward oxygen and outward metal diffusion are important for oxidation. Furthermore, by investigating the impact of the order of the initial oxide layers, we find that Cr 2 O 3 serves as a barrier to outward Mn diffusion. Moreover, controlling the initial oxide layer order can be used to limit oxidation.

36 MATERIALS SCIENCE↗
DOE OSTI2011

Phase-field investigation on the non-equilibrium interface dynamics of rapid alloy solidification

The research program reported here is focused on critical issues that represent conspicuous gaps in current understanding of rapid solidification, limiting our ability to predict and control microstructural evolution (i.e. morphological dynamics and microsegregation) at high undercooling, where conditions depart significantly from local equilibrium. More specifically, through careful application of phase-field modeling, using appropriate thin-interface and anti-trapping corrections and addressing important details such as transient effects and a velocity-dependent (i.e. adaptive) numerics, the current analysis provides a reasonable simulation-based picture of non-equilibrium solute partitioning and the corresponding oscillatory dynamics associated with single-phase rapid solidification and show that this method is a suitable means for a self-consistent simulation of transient behavior and operating point selection under rapid growth conditions. Moving beyond the limitations of conventional theoretical/analytical treatments of non-equilibrium solute partitioning, these results serve to substantiate recent experimental findings and analytical treatments for single-phase rapid solidification. The departure from the equilibrium solid concentration at the solid-liquid interface was often observed during rapid solidification, and the energetic associated non-equilibrium solute partitioning has been treated in detail, providing possible ranges of interface concentrations for a given growth condition. Use of these treatments for analytical description of specific single-phase dendritic and cellular operating point selection, however, requires a model for solute partitioning under a given set of growth conditions. Therefore, analytical solute trapping models which describe the chemical partitioning as a function of steady state interface velocities have been developed and widely utilized in most of the theoretical investigations of rapid solidification. However, these solute trapping models are not rigorously verified due to the difficulty in experimentally measuring under rapid growth conditions. Moreover, since these solute trapping models include kinetic parameters which are difficult to directly measure from experiments, application of the solute trapping models or the associated analytic rapid solidification model is limited. These theoretical models for steady state rapid solidification which incorporate the solute trapping models do not describe the interdependency of solute diffusion, interface kinetics, and alloy thermodynamics. The phase-field approach allows calculating, spontaneously, the non-equilibrium growth effects of alloys and the associated time-dependent growth dynamics, without making the assumptions that solute partitioning is an explicit function of velocity, as is the current convention. In the research described here, by utilizing the phase-field model in the thin-interface limit, incorporating the anti-trapping current term, more quantitatively valid interface kinetics and solute diffusion across the interface are calculated. In order to sufficiently resolve the physical length scales (i.e. interface thickness and diffusion boundary length), grid spacings are continually adjusted in calculations. The full trajectories of transient planar growth dynamics under rapid directional solidification conditions with different pulling velocities are described. As a validation of a model, the predicted steady state conditions are consistent with the analytic approach for rapid growth. It was confirmed that rapid interface dynamics exhibits the abrupt acceleration of the planar front when the effect of the non-equilibrium solute partitioning at the interface becomes signi ficant. This is consistent with the previous linear stability analysis for the non-equilibrium interface dynamics. With an appropriate growth condition, the continuous oscillation dynamics was able to be simulated using continually adjusting grid spacings. This oscillatory dynamics including instantaneous jump of interface velocities are consistent with a previous phenomenological model by and a numerical investigation, which may cause the formation of banded structures. Additionally, the selection of the steady state growth dynamics in the highly undercooled melt is demonstrated. The transition of the growth morphology, interface velocity selection, and solute trapping phenomenon with increasing melt supersaturations was described by the phase-field simulation. The tip selection for the dendritic growth was consistent with Ivantsov's function, and the non-equilibrium chemical partitioning behavior shows good qualitative agreement with the Aziz's solute trapping model even though the model parameter(V D ) remains as an arbitrary constant. This work is able to show the possibility of comprehensive description of rapid alloy growth over the entire time-dependent non-equilibrium phenomenon.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Interpreting Electrochemical and Chemical Sodiation Mechanisms and Kinetics in Tin Antimony Battery Anodes Using in Situ Transmission Electron Microscopy and Computational Methods

Intermetallic compounds such as SnSb are promising anode materials for sodium ion batteries; however, their nanoscale sodiation mechanisms are not well understood. Here, we used a combination of in situ transmission electron microscopy (TEM), first-principles electronic structure calculations, computational thermodynamic modeling, and phase-field simulations to reveal the sodiation mechanisms and to quantify microstructural effects contributing to the underlying reaction kinetics in SnSb electrodes. During in situ sodiation experiments, the nanocrystalline SnSb thin films underwent a rapid amorphous phase transformation upon sodiation, as determined by in situ TEM and electron diffraction experiments. The Na + diffusion coefficients were measured with and without an external electrical bias, and the data showed that an applied potential increased Na + diffusion by an order of magnitude compared to solid-state diffusion. Furthermore, there was a distinct decrease in sodium diffusion upon the formation of the amorphous phase that resulted from a change in the local structure and grain boundaries. To further understand how the Na + transport mechanism correlated with the changes observed in the SnSb thin films, phase-field modeling was used, which considered sodium diffusion within the grain boundaries together with their evolution and stress–strain state. As a result, these findings enhance our understanding of sodiation mechanisms within intermetallic anode materials for sodium ion battery applications.

25 ENERGY STORAGE↗
DOE OSTI2015

An atomistically informed mesoscale model for growth and coarsening during discharge in lithium-oxygen batteries

An atomistically informed mesoscale model is developed for the deposition of a discharge product in a Li-O{sub 2} battery. This mescocale model includes particle growth and coarsening as well as a simplified nucleation model. The model involves LiO{sub 2} formation through reaction of O{sub 2}{sup −} and Li{sup +} in the electrolyte, which deposits on the cathode surface when the LiO{sub 2} concentration reaches supersaturation in the electrolyte. A reaction-diffusion (rate-equation) model is used to describe the processes occurring in the electrolyte and a phase-field model is used to capture microstructural evolution. This model predicts that coarsening, in which large particles grow and small ones disappear, has a substantial effect on the size distribution of the LiO{sub 2} particles during the discharge process. The size evolution during discharge is the result of the interplay between this coarsening process and particle growth. The growth through continued deposition of LiO{sub 2} has the effect of causing large particles to grow ever faster while delaying the dissolution of small particles. The predicted size evolution is consistent with experimental results for a previously reported cathode material based on activated carbon during discharge and when it is at rest, although kinetic factors need to be included. The approach described in this paper synergistically combines models on different length scales with experimental observations and should have applications in studying other related discharge processes, such as Li{sub 2}O{sub 2} deposition, in Li-O{sub 2} batteries and nucleation and growth in Li-S batteries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗