Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary effects”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Selective laser sintering and spark plasma sintering of (Zr,Nb,Ta,Ti,W)C compositionally complex carbide ceramics

Abstract Two advanced manufacturing processes, spark plasma sintering (SPS) and selective laser sintering (SLS), have been developed for synthesis of (Zr,Nb,Ta,Ti,W)C compositionally complex carbide (CCC) via reactive sintering of a powder mixture of constitute monocarbides. X‐ray diffraction analysis confirmed that the single‐phase CCC can be formed by both SPS and SLS. While a homogenous microstructure with uniform metal element distributions was developed during SPS, three‐layer microstructures with a thin TiC‐rich layer and two TaC‐rich layers along with a TiO 2 ‐rich surface layer containing W nanoparticles were formed during SLS. In addition, cellular structures with W, Zr, and Ti element segregation and dislocations on cell boundaries were observed in the SLS‐CCC sample, indicating the effect of nonequilibrium conditions on microstructure formation during laser melting followed by rapid cooling and solidification process. Compared to the SPS‐CCC sample, the SLS‐CCC showed enhanced hardness and reduced thermal conductivity, which may be related to their unique cellular structures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

NMACFoam (Numerical-model for Membrane-spacer Assemblies and Configurations in OpenFOAM) [SWR-23-57]

NMACFoam is a computational fluid dynamics tool for simulating flow and solids transport in membrane-spacer channels pertaining to reverse osmosis or nanofiltration applications. It is developed on the OpenFOAM framework and allows for simulating spacers with complex geometries. It includes porous wall boundary conditions for membranes and can accurately capture concentration polarization effects. One can perform large calculations at module-scale domain sizes using high performance computing resources with NMACFoam, which scales to 1000s of processors for problems with 1-10 million cells.

Sitaraman, Hariswaran↗
DOE OSTI2021

Efficient Hydrogen Delivery for Microbial Electrosynthesis via 3D-Printed Cathodes

The efficient delivery of electrochemically in situ produced H 2 can be a key advantage of microbial electrosynthesis over traditional gas fermentation. However, the technical details of how to supply large amounts of electric current per volume in a biocompatible manner remain unresolved. Here, we explored for the first time the flexibility of complex 3D-printed custom electrodes to fine tune H 2 delivery during microbial electrosynthesis. Using a model system for H 2 -mediated electromethanogenesis comprised of 3D fabricated carbon aerogel cathodes plated with nickel-molybdenum and Methanococcus maripaludis, we showed that novel 3D-printed cathodes facilitated sustained and efficient electromethanogenesis from electricity and CO 2 at an unprecedented volumetric production rate of 2.2 L CH4 /L catholyte /day and at a coulombic efficiency of 99%. Importantly, our experiments revealed that the efficiency of this process strongly depends on the current density. At identical total current supplied, larger surface area cathodes enabled higher methane production and minimized escape of H 2 . Specifically, low current density (<1 mA/cm 2 ) enabled by high surface area cathodes was found to be critical for fast start-up times of the microbial culture, stable steady state performance, and high coulombic efficiencies. Our data demonstrate that 3D-printing of electrodes presents a promising design tool to mitigate effects of bubble formation and local pH gradients within the boundary layer and, thus, resolve key critical limitations for in situ electron delivery in microbial electrosynthesis.

03 NATURAL GAS↗
DOE OSTI2020

Air quality impact of the Northern California Camp Fire of November 2018

The Northern California Camp Fire that took place in November 2018 was one of the most damaging environmental events in California history. Here, we analyze ground-based station observations of airborne particulate matter that has a diameter <2.5 µm (PM 2.5 ) across Northern California and conduct numerical simulations of the Camp Fire using the Weather Research and Forecasting model online coupled with chemistry (WRF-Chem). Simulations are evaluated against ground-based observations of PM 2.5 , black carbon, and meteorology, as well as satellite measurements, such as Tropospheric Monitoring Instrument (TROPOMI) aerosol layer height and aerosol index. The Camp Fire led to an increase in Bay Area PM 2.5 to over 50 µg m –3 for nearly 2 weeks, with localized peaks exceeding 300 µg m –3 . Using the Visible Infrared Imaging Radiometer Suite (VIIRS) high-resolution fire detection products, the simulations reproduce the magnitude and evolution of surface PM 2.5 concentrations, especially downwind of the wildfire. The overall spatial patterns of simulated aerosol plumes and their heights are comparable with the latest satellite products from TROPOMI. WRF-Chem sensitivity simulations are carried out to analyze uncertainties that arise from fire emissions, meteorological conditions, feedback of aerosol radiative effects on meteorology, and various physical parameterizations, including the planetary boundary layer model and the plume rise model. Downwind PM 2.5 concentrations are sensitive to both flaming and smoldering emissions over the fire, so the uncertainty in the satellite-derived fire emission products can directly affect the air pollution simulations downwind. Our analysis also shows the importance of land surface and boundary layer parameterization in the fire simulation, which can result in large variations in magnitude and trend of surface PM 2.5 . Inclusion of aerosol radiative feedback moderately improves PM 2.5 simulations, especially over the most polluted days. Results of this study can assist in the development of data assimilation systems as well as air quality forecasting of health exposures and economic impact studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Decentralized Voltage optimization Method for Distribution Networks with Distributed PVs

In order to solve the voltage violations caused by high penetration distributed photovoltaics (PVs), this paper proposes a decentralized voltage optimization method for distribution networks. Through the optimal dispatch of distributed PVs and reactive power resources, it realizes the voltage security and economic operation of distribution systems and facilitates the integration of massive distributed PVs. The proposed method is inspired by the idea of locational marginal price (LMP) and feasible cutting plane. It separates the optimization of boundary bus voltage and power flow in adjacent regions and uses the interaction of boundary variables’ marginal prices to obtain the optimal solution. The effectiveness and accuracy of the proposed method are demonstrated via simulation test on a practical l0kV feeder in China.

Xue, Yaosuo↗
DOE OSTI2022

The encasulation effect on thermal performance of micro-encapsulated phase change materials during energy absorption

Understanding the performance of the thermal storage materials is required to successfully meet the needs of next generation concentrating solar power projects. This study focuses on the thermal characteristics of the micro-encapsulated phase change material (MEPCM) with metallic core under high temperature. A numerical model was developed and verified to study the transient heat absorption process of MEPCMs made of tin, aluminum and copper core. Convective heat transfer rate, total energy absorption and volumetric energy density were presented. The effect of encapsulation on theses variables was analyzed under various boundary and initial conditions. Its impact on characteristic time variables was discussed. The method presented can be adopted to analyzing both manufacturing and application processes for MEPCMs of various compositions, during which melting and solidification happen. Given certain operating conditions, the results from the model can be referenced to obtain optimal thermal performance of these materials.

Benner, Jingru↗
DOE OSTI2026

Coupling Microstructural Evolution Simulations to Material Property Degradation Predictions for Plasma-Facing Materials

Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500 C to 1500 C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 D/m^2-s. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

Connect microstructure evolution to property degradation with validated simulation

Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500$^o$C to 1500$^o$C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 $\frac{D}{m^2s}$. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.

36 - MATERIALS SCIENCE↗
DOE OSTI2023

The Observed Effects of Cold Pools on Convection Triggering and Organization during DYNAMO/AMIE

Cold pools that form from existing convection can generate new clouds nearby, a process that has been suggested to contribute to cloud organization. However, some disagreement exists on the role of cold pools in cloud organization over tropical oceans, and much remains to be understood on how this role depends on the state of large-scale disturbances such as the Madden-Julian Oscillation (MJO). Here, this study addresses this question by examining the intraseasonal variability in the properties of cold pools and their effects on convection triggering using observations. The unique set of surface meteorology data and ground-based radars deployed during the DYNAMO/AMIE field campaign are used to identify cold pools and the associated spatiotemporal evolution of rainfall. The results suggest that cold pools enhance rainfall at their expected speed of 3–8 m s –1 , which includes the Doppler-shifting effect by the background wind. The cold pools also enhance rainfall mainly within a 20 km radius, especially in the direction of boundary-layer vertical wind shear. However, gravity waves appear to contribute predominantly to the propagation of convection over a wider range of distances due to their faster propagation speeds. The interface of convection triggering and propagation at different speeds of cold pools and gravity waves seems to contribute to cloud organization. The effectiveness of cold pools and gravity waves on cloud organization also depends highly on the MJO due to the combined effects of changes in environmental humidity and boundary to lower-tropospheric wind.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Effect of anisotropic thermal transport on tearing mode stability in negative versus positive triangularity plasmas

The combined effects of anisotropic thermal transport and the plasma shaping, including negative triangularity, on the n = 1 (n is the toroidal mode number) tearing mode (TM) stability are numerically investigated utilizing the MARS-F code [Liu et al., Phys. Plasmas 7, 3681–3690 (2000)]. While varying the plasma boundary triangularity, the TM stability is found to be dictated by the competing effects of the Shafranov shift induced stabilization and the bad-curvature induced destabilization. The negative triangularity shape increases the Shafranov shift (stabilizing) in the plasma core but also enlarges bad-curvature regions (destabilizing) near the plasma edge, with the net effect being largely destabilizing for the TM as compared to the positive triangularity counter-part. Large negative triangularity however can also lead to more stabilization for the plasma core-localized TM. Anisotropic thermal transport reduces the stabilizing effect on the TM associated with the favorable averaged curvature, resulting in more unstable core-localized TMs in both negative and positive triangularity plasmas. But the opposite effect can also take place for the edge-localized TM in finite-pressure plasmas with negative triangularity.

Physics↗
DOE OSTI2025

Pore Resolved Simulations of Joule Heating in Fibrous Media using an Embedded Boundary Method

Joule heating has been regarded as an energy-efficient and sustainable method for heating materials and gases at large scales. The modeling of local temperature effects at pore-resolved scales for such systems, however, has been difficult to achieve due to challenges in coupling thermo-chemical processes in complex porous media and in large representative volume elements (RVEs). To this end, we developed an electro-thermal model at the pore scale to study Joule heating effects in large heterogeneous systems with different microstructures. This was achieved using the level set method to implicitly delineate distinct regions within the domain, and an embedded boundary method to facilitate heat exchange across the fluid-solid interface. Moreover, we applied this method to investigate unsteady non-linear electro-thermal effects in non-woven fibrous graphite conductors for RVEs with characteristic lengths of 2 mm, with different fiber orientations, porosity (80% – 90%) and fiber diameters (10 – 20µm). The coupled equations were solved numerically and they produced peak temperatures greater than 2000 K resulting in heating rates as high as 80,000 K/s. Moreover, the results depended strongly on the microstructure of the fiber skeleton and current density. Geometries with large fibers (∼ 20µm) had the highest average and peak temperatures with the mean temperature increasing by 3.9 % while the peak temperature increased by 9.9 %. Anisotropic domains on the other hand had the lowest mean and peak temperatures with peak and mean temperatures of 2293 K and 1437.7K respectively representing a corresponding 12.1% and 5.1% drop in the temperatures. An increase in porosity from 80% to 90%, however, led to an increase in the peak temperature by 5.1%.

Joule heating↗
DOE OSTI2025

Experimental characterization and atomistic simulation of grain boundary segregation in Mg-Y alloys

As a rare earth solute element in Mg alloys, Y has the beneficial effects of increasing both the strength and the ductility as well as weakening the crystallographic texture. To achieve a more fundamental understanding on how Y addition affects the microstructural evolution and mechanical properties, the Y segregation behavior at grain boundaries was investigated in Mg-1wt.%Y and Mg-7wt.%Y alloys at different conditions. The segregation intensity and its dependence on the grain boundary misorientation angle were experimentally characterized and computationally predicted. Strong segregation at grain boundaries was observed in both low and high Y-containing alloys. Y segregation was found to remain in alloy Mg-7Y after high-temperature annealing heat treatment at 540 °C. No direct correlation between the Y segregation intensity and the grain boundary misorientation angle could be established based on either the experimental characterization or the atomistic simulation with a spectral model. We thus conclude that grain boundary segregation of Y is independent of grain boundary misorientation angle.

Grain boundary↗
DOE OSTI2021

Ribbon Wrinkling Origin and Its Effect on Tape-Wound Core Shape

Improved understanding of advanced manufacturing techniques is needed to fabricate next- generation soft magnetic tape-wound cores. These utilize metal amorphous nano-composites (MANCs) in the form of long thin ribbons. MANCs exhibit attractive soft magnetic properties for use in high frequency and high power density magnetic components. Of particular interest here, the permeability and hysteresis of MANC ribbons can be tuned in a strain annealing process. An unexpected consequence is that wider ribbons develop longitudinal wrinkle patterns. This can substantially reduce the tape-wound core stacking factor. However, when winding a wrinkled ribbon into a toroidal core, the wrinkles can completely flatten, or crinkle. By this phenomenon, the stacking factor can be recovered. To understand the crinkling effect we have developed analytical and finite element analysis (FEA) models. Experiments with known boundary conditions demonstrate the transition from wrinkled to crinkled versus winding curvature. The models utilize an elastic modulus measured via ribbon tensile testing and a ribbon thickness value measurement made by profilometry. The analytical model provides insight into the crinkling effect but cannot explain the length of the transition region from wrinkling to crinkling. With FEA, we successfully predict the transition region and depending on wrinkling amplitude can determine when a ribbon will wind into a useable core. Such mechanics-based insight provides critical feedback for allowable processing parameters in the strain annealing process.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2022

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

On the Interfacial Assembly of Anisotropic Amphiphilic Janus Particles

It has been shown both theoretically and experimentally that amphiphilic Janus particles are the most effective solid surfactants to stabilize interfaces. In most cases, the Janus particles investigated have uniform morphologies with Janus boundaries dividing the particle into halves. However, there are many examples of Janus particles where the hydrophilic and hydrophobic domains are not equally distributed. The effects of this uneven domain distribution on the mechanism and kinetics of Janus particle assembly, and final equilibrium state are not well-understood. Dynamic pendant drop tensiometry offers a means to probe both the equilibrium assembly and the kinetics and mechanism of assembly. Here, in this work, the interfacial kinetics and assembly of spherical anisotropic Janus particles are investigated using dynamic pendant drop tensiometry. Systematic studies quantifying the time-dependent interfacial behavior as a function of Janus particle morphology, chemical composition, particle concentration, and NaOH and HCl concentration are performed. These studies shed light on the assembly mechanism of more complex Janus particle morphologies and highlight their effectiveness as interface stabilizers.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Predicting interface structure using the minima hopping method

Here, we adapt the minima hopping method (MHM) to the problem of interfacial structure prediction and apply it to study a canonical problem, the tilt grain boundaries in SrTiO 3 . Our method employs a hybrid approach by first exploring the potential energy surface (PES) of different grain boundary samplings with an empirical force field, among which the fifteen candidates with lower energies are then refined using ab initio density functional theory (DFT) calculations. During the exploratory stage, we bias the search using a local order parameter to primarily sample various reconstructions in the vicinity of the interface, while preserving the crystallinity of the bulk regions. We further enhance the search by incorporating initial structures with rigid body displacements to account for translational variations between bulk phases, enabling the MHM to effectively generate both stoichiometric and nonstoichiometric SrTiO 3 Σ⁢3(111)[110] and Σ⁢3(112)[110] grain boundaries. From an algorithmic standpoint, MHM outperforms earlier studies based on genetic algorithms (GA) by identifying more stable interfacial structures of several SrTiO 3 grain boundaries. The performance of the present implementation of the MHM approach is primarily limited by exploring an approximate description of the PES with a rather simple Buckingham potential. This limitation leads to variations in performance when compared to approaches utilizing more advanced surrogate PES models, such as direct DFT-PES sampling or GA with the embedded atom method (EAM). Despite the present limitations, the MHM approach is able to yield interfacial structures with comparable or lower interfacial energies in specific cases, such as Σ⁢3(111)[110] Γ=1, ±0.5 and Σ⁢3(112)[110] Γ= ±1, −2, underscoring the robustness of the MHM approach even with a simple approximation of the DFT PES. The MHM interfacial structure prediction method thus offers an efficient approach to understanding the grain boundaries and heterointerfaces at the atomic scale, providing an important prerequisite for effective materials design.

density functional theory↗
DOE OSTI2026

The quantum evolutions of the diffractive transverse-momentum dependent gluon distribution

Using the Colour Glass Condensate description of electron-nucleus collisions at high energy, we study the diffractive production of a pair of jets with transverse momenta much larger than the nuclear saturation momentum Q s . At leading order in the QCD coupling, the di-jet cross-section exhibits transverse-momentum dependent (TMD) factorisation, with a gluon diffractive TMD distribution (DTMD) which is controlled by gluon saturation and describes the transverse-momentum imbalance between the produced jets. The next-to-leading corrections generate the various quantum evolutions of the diffractive gluon distribution. We focus on the Collins-Soper-Sterman (CSS) evolution which describes the change in the gluon DTMD when increasing the “hard scale” (the typical transverse momentum of the di-jets). We consider two different representations for this equation, one in transverse-momentum space, the other one in transverse-coordinate space. They are not fully equivalent with each other (despite being related by a Fourier transform) because of the respective boundary conditions. These conditions encode the essential physics of gluon saturation together with the effects of two other types of quantum evolution: the BK/JIMWLK evolution over the rapidity gap (“inside the Pomeron”) and the DGLAP evolution outside the rapidity gap (“within the diffractive system”). We demonstrate that, due to gluon saturation, one can compute both the boundary conditions and the CSS solutions mostly from first principles, without the need for a non-perturbative Sudakov. We numerically find a good agreement between the CSS solutions in the two aforementioned representations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Statistically Steady State Large-Eddy Simulations Forced by an Idealized GCM: 1. Forcing Framework and Simulation Characteristics

Using large-eddy simulations (LES) systematically has the potential to inform parameterizations of subgrid-scale (SGS) processes in general circulation models (GCMs), such as turbulence, convection, and clouds. Here we show how LES can be run to emulate grid columns of GCMs to generate a library of LES across a cross-section of dynamical regimes. The LES setup replicates the thermodynamic and water budgets in GCM grid columns. Resolved horizontal and vertical transports of heat and water and large-scale pressure gradients from the GCM are prescribed as forcing in the LES. The LES satisfies the same (slab-ocean) surface boundary conditions as the GCM, leaving the LES temperatures free to adjust. Radiative transfer is treated in a unied but highly idealized manner (a semi-gray atmosphere without cloud radiative effects) in both the GCM and LES. We show that the LES with these forcing and boundary conditions reaches statistically steady states without nudging to reference profiles. These steady states provide a training dataset for developing GCM parameterizations. The same LES setup also provides a good basis for studying the cloud response to global warming.

54 ENVIRONMENTAL SCIENCES↗