Search NASA⌕ Search

SEARCH · Search NASA

Results for “process simulation 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 901 records · Page 50

Matrix Approach to Accelerate Spin-Up of CLM5

Numerical models have been developed to investigate and understand responses of biogeochemical cycle to global changes. Steady state, when a system is in dynamic equilibrium, is generally required to initialize these model simulations. However, the spin-up process that is used to achieve steady state pose a great burden to computational resources, limiting the efficiency of global modeling analysis on biogeochemical cycles. This study introduces a new Semi-Analytical Spin-Up (SASU) to tackle this grand challenge. We applied SASU to Community Land Model version 5 and examined its computational efficiency and accuracy. At the Brazil site, SASU is computationally 7 times more efficient than (or saved up to 86% computational cost in comparison with) the traditional native dynamics (ND) spin-up to reach the same steady state. Globally, SASU is computationally 8 times more efficient than the accelerated decomposition spin-up and 50 times more efficient than ND. In summary, SASU achieves the highest computational efficiency for spin-up on site and globally in comparison with other spin-up methods. It is generalizable to wide biogeochemical models and thus makes computationally costly studies (e.g., parameter perturbation ensemble analysis and data assimilation) possible for a better understanding of biogeochemical cycle under climate change.

54 ENVIRONMENTAL SCIENCES↗

A Gaussian process autoregressive model capturing microstructure evolution paths in a Ni–Mo–Nb alloy

Additive manufacturing is increasingly being employed to produce components of complex geometries in structural alloys because of the expected energy savings associated with the near-net-shape capability and the ability to build in novel internal features that are not possible with many conventional manufacturing approaches. However, because of the extreme thermal conditions encountered, the non-equilibrium microstructures produced during powder bed-based additive manufacturing processes must be subjected to custom post-heat treatment processes to recover the target mechanical properties. Phase-field models and simulation techniques have matured to a state where the microstructure evolution paths, and the morphologies of the resulting precipitate phases can be predicted reasonably accurately, considering alloy-specific thermodynamic and kinetic aspects of the nucleation and growth processes. However, phase-field simulations are computationally intensive, which precludes the ability to apply the simulations directly to the length scale of the entire component. Therefore, it is highly desirable to develop low-computational-cost surrogate models that effectively capture the physics at the microstructural length scale, while facilitating the design of optimized processing conditions resulting in location-specific targeted microstructures at the component scale. The work presented here demonstrates the application of the materials knowledge system framework to develop a surrogate model that effectively captures the microstructural path during annealing of a Ni–Mo–Nb alloy containing different Mo and Nb compositions known to segregate during solidification under additive manufacturing conditions. Specifically, the surrogate model built in this work is based on a Gaussian process autoregressive model informed by statistical representation of simulated microstructures using two-point correlations and dimensionality reduction through principal component analysis. In conclusion, this surrogate model is shown to capture the bifurcation of the microstructural path during precipitation, which yields a microstructure dominated by the $\gamma^{\prime\prime}$ phase at high Nb concentrations and the $\delta$ phase at low Nb concentrations.

36 MATERIALS SCIENCE↗

Virtual Engineering: Python framework for engineering process design

Virtual Engineering (VE) is a Python software framework designed to accelerate the research and development of engineering processes that are fundamentally defined by multiple unit operations executed in series. VE supports a wide variety of different multi-physics models and integrates them to simulate a complete end-to-end process. To automate the execution of this model sequence, VE provides (i) a robust method to communicate between models, (ii) a high-level, user-friendly interface to set model parameters and enable optimization, and (iii) an overall model-agnostic approach that allows new computational units to be swapped in and out of workflows. Although the VE framework was developed to support the biochemical conversion of biomass to fuel, we have designed each component to easily accommodate new domains and unit models.

09 BIOMASS FUELS↗

How Galaxies Form Stars: The Connection between Local and Global Star Formation in Galaxy Simulations

Using a suite of isolated $L_\star$ galaxy simulations, we show that global depletion times and star-forming gas mass fractions in simulated galaxies exhibit systematic and well-defined trends as a function of the local star formation efficiency per freefall time, $\epsilon_{\rm ff}$, strength of stellar feedback, and star formation threshold. We demonstrate that these trends can be reproduced and explained by a simple physical model of global star formation in galaxies. Our model is based on mass conservation and the idea of gas cycling between star-forming and non-star-forming states on certain characteristic time scales under the influence of dynamical and feedback processes. Both the simulation results and our model predictions exhibit two limiting regimes with rather different dependencies of global galactic properties on the local parameters. When $\epsilon_{\rm ff}$ is small and feedback is inefficient, the total star-forming mass fraction, $f_{\rm sf}$, is independent of $\epsilon_{\rm ff}$ and the global depletion time, $\tau_{\rm dep}$, scales inversely with $\epsilon_{\rm ff}$. When $\epsilon_{\rm ff}$ is large or feedback is very efficient, these trends are reversed: $f_{\rm sf} \propto \epsilon_{\rm ff}^{-1}$ and $\tau_{\rm dep}$ is independent of $\epsilon_{\rm ff}$ but scales linearly with the feedback strength. We also compare our results with the observed depletion times and mass fractions of star-forming and molecular gas and show that they provide complementary constraints on $\epsilon_{\rm ff}$ and the feedback strength. In conclusion, we show that useful constraints on $\epsilon_{\rm ff}$ can also be obtained using measurements of the depletion time and its scatter on different spatial scales.

79 ASTRONOMY AND ASTROPHYSICS↗

Surrogate-driven design optimization with uncertainty constraints in Monte Carlo simulations

In multi-objective design tasks, the computational cost increases rapidly when high-fidelity simulations are used to evaluate objective functions. Surrogate models help mitigate this cost by approximating the simulation output, simplifying the design process. However, under high uncertainty, surrogate models trained on noisy data can produce inaccurate predictions, as their performance depends heavily on the quality of training data. This study investigates the impact of data uncertainty on two multi-objective design problems modelled using Monte Carlo transport simulations: a neutron moderator and an ion-to-neutron converter. For each, a grid search was performed using five different tally uncertainty levels to generate training data for neural network surrogate models. These models were then optimized using NSGA-III. The recovered Pareto-fronts were analyzed across uncertainty levels: in the moderator problem, normalized hypervolume dropped from 0.886 at 1.0% uncertainty to 0.748 at 10% uncertainty, while in the converter problem it remained near 0.50 for all cases. Average simulation times were also compared to evaluate the trade-off between accuracy and computational cost. Results show that the influence of simulation uncertainty is strongly problem-dependent. In the neutron moderator case, higher uncertainties led to exaggerated objective sensitivities and distorted Pareto-fronts, reducing normalized hypervolume. In contrast, the ion-to-neutron converter task was less affected—low-fidelity simulations produced results similar to those from high-fidelity data. These findings suggest that a fixed-fidelity approach is not optimal. Surrogate models can recover the Pareto-front under noisy conditions, and multi-fidelity studies help identify suitable uncertainty levels for each problem to balance efficiency and accuracy.

07 ISOTOPE AND RADIATION SOURCES↗

Representation of Nucleation Mode Microphysics in a Global Aerosol Model with Sectional Microphysics

In models, nucleation mode (1 nm<Dp <10 nm) particle microphysics can be represented explicitly with aerosol microphysical processes or can be parameterized to obtain the growth and survival of nuclei to the model's lower size boundary. This study investigates how the representation of nucleation mode microphysics impacts aerosol number predictions in the TwO-Moment Aerosol Sectional (TOMAS) aerosol microphysics model running with the GISS GCM II-prime by varying its lowest diameter boundary: 1 nm, 3 nm, and 10 nm. The model with the 1 nm boundary simulates the nucleation mode particles with fully resolved microphysical processes, while the model with the 10 nm and 3 nm boundaries uses a nucleation mode dynamics parameterization to account for the growth of nucleated particles to 10 nm and 3 nm, respectively.We also investigate the impact of the time step for aerosol microphysical processes (a 10 min versus a 1 h time step) to aerosol number predictions in the TOMAS models with explicit dynamics for the nucleation mode particles (i.e., 3 nm and 1 nm boundary). The model with the explicit microphysics (i.e., 1 nm boundary) with the 10 min time step is used as a numerical benchmark simulation to estimate biases caused by varying the lower size cutoff and the time step. Different representations of the nucleation mode have a significant effect on the formation rate of particles larger than 10 nm from nucleated particles (J10) and the burdens and lifetimes of ultrafinemode (10 nm≤Dp ≤70 nm) particles but have less impact on the burdens and lifetimes of CCN-sized particles. The models using parameterized microphysics (i.e., 10 nm and 3 nm boundaries) result in higher J10 and shorter coagulation lifetimes of ultrafine-mode particles than the model with explicit dynamics (i.e., 1 nm boundary). The spatial distributions of CN10 (Dp ≥10 nm) and CCN(0.2 %) (i.e., CCN concentrations at 0.2%supersaturation) are moderately affected, especially CN10 predictions above ∼700 hPa where nucleation contributes most strongly to CN10 concentrations. The lowermost-layer CN10 is substantially improved with the 3 nm boundary (compared to 10 nm) in most areas. The overprediction in CN10 with the 3 nm and 10 nm boundaries can be explained by the overprediction of J10 or J3 with the parameterized microphysics, possibly due to the instantaneous growth rate assumption in the survival and growth parameterization. The errors in CN10 predictions are sensitive to the choice of the lower size boundary but not to the choice of the time step applied to the microphysical processes. The spatial distribution of CCN(0.2 %) with the 3 nm boundary is almost identical to that with the 1 nm boundary, but that with the 10 nm boundary can differ more than 10-40% in some areas. We found that the deviation in the 10 nm simulations is partly due to the longer time step (i.e., 1 h time step used in the 10 nm simulations compared to 10 min time step used in the benchmark simulations), but, even with the same time step, the 10 nm cutoff showed noticeably higher errors than the 3 nm cutoff. In conclusion, we generally recommend using a lower diameter boundary of 3 nm for studies focused on aerosol indirect effects but down to 1 nm boundary for studies focused on CN10 predictions or nucleation.

microphysics↗

Two-stage membrane-based process utilizing highly CO 2 -selective membranes for cost and energy efficient carbon capture from coal flue gas: A process simulation study

Membrane technology for CO 2 capture has become an attractive strategy due to its cost and energy efficiency and low materials costs. In the past decade, membrane-based process designs for post-combustion power plant CO 2 capture have been developed, utilizing existing highly CO 2 -permeable membranes with relatively low CO 2 /N 2 selectivity (<50), and have obtained reasonably economic carbon capture. However, few membrane-based process designs were proposed for moderate to highly CO 2 -selective membranes (CO 2 /N 2 selectivity of 50–300, and >300, respectively), which have vastly emerged in recent years, such as various facilitated transport membranes (FTMs). Herein, we proposed a two-stage membrane-base process design targeting economic carbon capture from coal-fired flue gas. This process design features the utilization of highly CO 2 -selective membranes for one-stage CO 2 enrichment to 95% dry-base purity in the first stage and recycle of the remaining CO 2 by a highly CO 2 -permeable membrane in the second stage, in order to achieve economic CO 2 capture with 90% capture rate and >95% CO 2 product purity. Through an integration-iteration membrane model and the Aspen Plus process simulation, a sensitivity study of operating pressures (feed and permeate pressures) and membrane properties (CO 2 permeance and CO 2 /N 2 selectivity) was conducted. Critical CO 2 /N 2 selectivity of 300–400 was found for the highly CO 2 -selctive membranes to meet the demand for cost and energy efficient results. The lowest possible membrane area of 4.8 × 10 5 m 2 and fractional energy of 19.3% were obtained, which is comparable to or even more attractive than reported membrane-based process designs. Here, this work provides a new membrane process design option for highly CO 2 -selective membranes and gives insights on the influence of membrane performance and operation condition.

42 ENGINEERING↗

Thermodynamic properties of the heliospheric plasma

A global-scale study is conducted of the radial evolution of the solar wind's temperature (T), number density (n), and entropy (s) thermodynamic properties, on the basis of plasma data from such sources as IMP and ISEE 3, in conjunction with an MHD simulation model. The radial evolution of T and s indicates that the heliospheric plasma is subjected to a heating process; this heating is interpreted to be largely due to the cumulative effect of the shock process with increasing distance from the sun. This simulation model is also used to extrapolate the thermodynamic properties of the solar wind from Voyager 2 to the region of the outer heliosphere that is bounded by the termination shock.

Whang, Y. C.↗

Impact of Asian Aerosols on Precipitation Over California: An Observational and Model Based Approach

Dust and pollution emissions from Asia are often transported across the Pacific Ocean to over the western United States. Therefore, it is essential to fully understand the impact of these aerosols on clouds and precipitation forming over the eastern Pacific and western United States, especially during atmospheric river events that account for up to half of California's annual precipitation and can lead to widespread flooding. In order for numerical modeling simulations to accurately represent the present and future regional climate of the western United States, we must account for the aerosol-cloud-precipitation interactions associated with Asian dust and pollution aerosols. Therefore, we have constructed a detailed study utilizing multi-sensor satellite observations, NOAA-led field campaign measurements, and targeted numerical modeling studies where Asian aerosols interacted with cloud and precipitation processes over the western United States. In particular, we utilize aerosol optical depth retrievals from the NASA Moderate Resolution Imaging Spectroradiometer (MODIS), NOAA Geostationary Operational Environmental Satellite (GOES-11), and Japan Meteorological Agency (JMA) Multi-functional Transport Satellite (MTSAT) to effectively detect and monitor the trans-Pacific transport of Asian dust and pollution. The aerosol optical depth (AOD) retrievals are used in assimilating the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) in order to provide the model with an accurate representation of the aerosol spatial distribution across the Pacific. We conduct WRF-Chem model simulations of several cold-season atmospheric river events that interacted with Asian aerosols and brought significant precipitation over California during February-March 2011 when the NOAA CalWater field campaign was ongoing. The CalWater field campaign consisted of aircraft and surface measurements of aerosol and precipitation processes that help extensively validate our WRF-Chem model simulations. After validating the capability of the WRF-Chem in realistically simulating the aerosol-cloud precipitation interactions, we conduct sensitivity studies where the AOD is doubled to diagnose whether an increasing concentration of Asian aerosols over the western United States will lead to further impacts on the cloud and precipitation processes over California. We also perform sensitivity studies where the aerosols will be partitioned into dust-only and pollution-only in order to separate the impacts of the differing Asian aerosol species. The results of our WRF-Chem model simulations aim to show that the trans-Pacific transport of Asian aerosols influence the precipitation associated with atmospheric river events that can ultimately impact the regional climate of the western United States. 1 University

Naeger, Aaron R.↗

Meson exchange current (MEC) models in neutrino interaction generators

Understanding of the so-called 2 particle-2 hole (2p-2h) effect is an urgent program in neutrino interaction physics for current and future oscillation experiments. Such processes are believed to be responsible for the event excesses observed by recent neutrino experiments. The 2p-2h effect is dominated by the meson exchange current (MEC), and is accompanied by a 2-nucleon emission from the primary vertex, instead of a single nucleon emission from the charged-current quasi-elastic (CCQE) interaction. Current and future high resolution experiments can potentially nail down this effect. For this reason, there are world wide efforts to model and implement this process in neutrino interaction simulations. In these proceedings, I would like to describe how this channel is modeled in neutrino interaction generators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Evaluation of top-down and bottom-up global terrestrial respiration estimates and their mismatch with model simulations

Terrestrial respiration is one of the most poorly understood processes in the global carbon cycle, making respiration predictions uncertain. However, expanding observations and machine learning approaches have led to a proliferation of estimates. We compiled total ecosystem and heterotrophic respiration estimates derived from top-down atmospheric inversions and bottom-up upscaling of ecosystem observations and compared them with dynamic vegetation models (DGVM) simulations over the 1980-2020 period. Our analysis revealed a convergence in mean annual global total ecosystem respiration estimates between top-down 97.1 (± SD 6.8) PgC yr-1 and bottom-up 98.5 (+/-13.4) PgC yr-1, which were both significantly lower than the ensemble mean from DGVMs estimates 133.7 (±4.7) PgC yr-1. We also found similar temporal trends between top-down estimates with a mean of 0.075 (±0.05) PgC yr-2, and bottom-up estimates of 0.05 (±0.05) PgC yr-2, which were 5 to 7 times smaller than the ensemble mean trend of 0.34 PgC yr-2 simulated by DGVMs. Global heterotrophic respiration showed much less agreement, ranging from top-down estimates of 42.7 (±4.0) PgC yr-1 to bottom-up estimates of 51.5 (±4.0) PgC yr-1 and a significantly larger ensemble model mean estimate of 60.8 PgC yr-1 (±1.9). The temporal trends in observation-based bottom-up estimates of heterotrophic respiration of 0.03 PgC yr-2 were five times lower than the model ensemble mean 0.15 PgC yr-2. Large regional disagreements in heterotrophic respiration estimates and simulations were evident in tropical and boreal latitudes. Therefore, improved regional and heterotrophic respiration estimates are necessary to reduce uncertainties regarding the future vulnerability of soil carbon.

Ballantyne, Ashley↗

Contrasting the Biophysical and Radiative Effects of Rising CO 2 Concentrations on Ozone Dry Deposition Fluxes

Abstract The dry deposition of ozone from the atmosphere to ecosystems is an important coupling mechanism between atmospheric chemistry and terrestrial biogeochemical processes. In most Earth system models, dry deposition is simulated using a resistor‐in‐series approach that aims to parameterize the governing biological, chemical, and physical processes through a series of functional approximations. Here, we evaluate the influence of carbon cycle‐climate responses on this parameterization using the results of the Energy Exascale Earth System Model v1.1 Biogeochemistry simulation campaign. This simulation campaign was designed in part to explore the biophysical and radiative effects of rising historical CO 2 concentrations on the Earth system. We find that while the global annual ozone dry deposition is relatively insensitive to these effects, regionally the influence can be up to 10%. The strongest regional sensitivities in ozone dry deposition are predominantly in higher latitudes over land in the northern hemisphere and are dominated by the radiative effect of CO 2 , with little net influence of biophysical responses. Of all the impacts of the radiative effect of CO 2 , we point to the potential importance of accurately representing ozone deposition to snow in Earth System Models and provide recommendations for future simulation campaigns.

54 ENVIRONMENTAL SCIENCES↗

COCOMO2: A Coarse-Grained Model for Interacting Folded and Disordered Proteins

Biomolecular interactions are essential in many biological processes, including complex formation and phase separation processes. Coarse-grained computational models are especially valuable for studying such processes via simulation. Here, we present COCOMO2, an updated residue-based coarse-grained model that extends its applicability from intrinsically disordered peptides to folded proteins. This is accomplished with the introduction of a surface exposure scaling factor, which adjusts interaction strengths based on solvent accessibility, to enable the more realistic modeling of interactions involving folded domains without additional computational costs. COCOMO2 was parametrized directly with solubility and phase separation data to improve its performance on predicting concentration-dependent phase separation for a broader range of biomolecular systems compared to the original version. COCOMO2 enables new applications including the study of condensates that involve IDPs together with folded domains and the study of complex assembly processes. COCOMO2 also provides an expanded foundation for the development of multiscale approaches for modeling biomolecular interactions that span from residue-level to atomistic resolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Correction to “COCOMO2: A Coarse-Grained Model for Interacting Folded and Disordered Proteins”

Biomolecular interactions are essential in many biological processes, including complex formation and phase separation processes. Coarse-grained computational models are especially valuable for studying such processes via simulation. Here, we present COCOMO2, an updated residue-based coarse-grained model that extends its applicability from intrinsically disordered peptides to folded proteins. This is accomplished with the introduction of a surface exposure scaling factor, which adjusts interaction strengths based on solvent accessibility, to enable the more realistic modeling of interactions involving folded domains without additional computational costs. COCOMO2 was parametrized directly with solubility and phase separation data to improve its performance on predicting concentration-dependent phase separation for a broader range of biomolecular systems compared to the original version. COCOMO2 enables new applications including the study of condensates that involve IDPs together with folded domains and the study of complex assembly processes. COCOMO2 also provides an expanded foundation for the development of multiscale approaches for modeling biomolecular interactions that span from residue-level to atomistic resolution.

Molecular interactions↗

Physics-informed machine-learning model of temperature evolution under solid phase processes

We model temperature dynamics during Shear Assisted Proccess Extrusion (ShAPE), a solid phase process that plasticizes feedstock with a rotating tool and subsequently extrudes it into a consolidated tube, rod, or wire. Control of temperature is critical during ShAPE processing to avoid liquefaction, ensure smooth extrusion, and develop desired material properties in the extruded products. Accurate modeling of the complicated thermo-mechanical feedbacks between process inputs, material temperature, and heat generation presents a significant barrier to predictive modeling and process design. In particular, connecting micro-structural scale mechanisms of heat generation to macro-scale predictions of temperature can become computationally intractable. In this work we use a neural network (NN) model of heat generation to bridge this gap, by combining it with a simplified model of the temperature dynamics due to conduction and convection to capture the macro scale evolution of temperature. We inform the construction of the NN heat generation model using crystal plasticity simulations at the micro-structural scale to model the effects of process inputs on generation of heat. We achieved close fits of the temperature dynamics model to a diverse experimental data-set. Further, the relationships learned by the NN model between process inputs and heat generation showed qualitative agreement with those predicted by crystal plasticity simulations.

36 MATERIALS SCIENCE↗

Formation of Arctic Stratocumuli Through Atmospheric Radiative Cooling

Stratocumulus clouds are important to the Arctic climate because they are prevalent and exert a strong radiative forcing on the surface. However, relatively little is known about how stratocumulus clouds form in the Arctic. In this study, radiative transfer calculations are used to show that the timescale over which stably stratified Arctic temperature and water vapor profiles cool to saturation is less than typical residence times for individual air parcels in the Arctic. This result is consistent with previous studies in suggesting that elevated stratocumulus can form naturally through clear-sky radiative cooling during all seasons, without assistance from frontal lifting or other atmospheric forcing. Single column model simulations of the cloud formation process, after radiative cooling has resulted in saturation in a stably stratified profile, suggest that stratocumulus cloud properties are sensitive to the characteristics of the environment in which the formation process takes place. For example, sensitivity tests suggest that clouds may attain liquid water paths of over 50 g/m 2 if they form in moist environments but may become locked in a low-liquid water path quasi steady state or dissipate within hours if they form in dry environments. A potential consequence of these sensitivities is that when an Arctic stratocumulus layer forms by radiative cooling, it is more likely to become optically thick, optically thin, or dissipate than it is to obtain an intermediate optical thickness. This could help explain why the cloudy and radiatively clear atmospheric states are so prevalent across the Arctic.

54 ENVIRONMENTAL SCIENCES↗

Data and Scripts Associated with "Modeling Ecohydrological Responses of Vegetation to Urban Microclimates Using the E3SM Land Model"

This dataset supports the study of vegetation ecohydrological responses to urban microclimates using the land component of the Energy Exascale Earth System Model (ELM) at four urban sites in Knoxville, Tennessee, USA. It includes the model inputs, simulation outputs, and associated scripts for running ELM simulations and analyzing the resulting data. The Model_Inputs folder includes static surface data, satellite-derived phenology (i.e., leaf area index), and atmospheric forcing data used to drive ELM simulations. Detailed descriptions of these datasets are provided in Section 2.3.2 of the associated manuscript. The Model_Outputs folder contains simulation results for the baseline, treatment, and ensemble experiments. Outputs from the baseline and treatment simulations are provided as raw ELM NetCDF files. Because the raw outputs from the 4,000-member ensemble are prohibitively large, the ensemble results are provided as summarized CSV files, which also serve as the source data for Figure 5 of the associated manuscript. The Scripts folder contains three components: E3SM, the core codebase of the Energy Exascale Earth System Model (E3SM); elm-olmt, the Offline Land Model Testbed (OLMT) used to perform the simulations; and knoxville_elm, which contains the analysis scripts used to process model outputs and generate the figures and results presented in the associated manuscript. Additional information is provided in Scripts_readme.txt within the Scripts directory.

Lu, Xiaoman [ORNL] (ORCID:0000000306698780)↗

Sensitivity of Precipitation Processes to Microphysics and Resolution in a Cloud-resolving Model

The Goddard Cumulus Ensemble (GCE) model is used to examine the impact of various microphysical schemes, and vertical and horizontal resolution on the development, intensity and rainfall associated with mesoscale convective systems, idealized hurricanes and an ensemble of clouds. The model variables include horizontal and vertical velocities, potential temperature, perturbation pressure, turbulent kinetic energy, and mixing ratios of all water phases (vapor, liquid, and ice). The major characteristics of the GCE model are the explicit representation of warm rain and ice microphysical processes, and their complex interactions with solar and infrared radiative transfer processes. For idealized hurricane, an axisymmetric version of the GCE model was developed and used successfully to simulate the tropical cyclogenesis process using both a Rankin vortex and saturated air within a specified radius as initial conditions. For mesoscale convective systems, the 3-D version of the GCE model was used to simulate squall lines that developed in the western Pacific, South China Sea, eastern Atlantic, South America and central US. For the cloud ensemble, the GCE model was integrated for several days in order to have a good sampling of cloud statistics.

Tao, Wei-Kuo↗