Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermodynamic consistency”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Experimental and theoretical analysis of carbon driven detonation waves in a heterogeneously premixed Rotating Detonation Engine

Coal dust explosions can be hazardous; however, they can also generate a significant rise in stagnation pressure if adequately harnessed. Rotating detonation combustors seek to take advantage of the stagnation pressure rise phenomenon in a more sustained and controlled manner via confinement to a physical annulus, leading to increased overall thermodynamic efficiency. Here this investigation presents an analysis of detonations fueled by Carbon Black, a solid particulate consisting of virtually pure carbon molecules and lean Hydrogen-Air mixtures. It is realized that with the addition of Carbon Black, an increase of lean mixture detonability and detonation velocities extending the operating limit over that of a pure hydrogen-air mixture is experienced. For all testing conditions, the total equivalence ratio is held at φ = 1, while the fuel mixture's carbon mass fraction is increased from 0 to 0.7 while the hydrogen is decreased. Detonation wave velocities are extracted from high-speed imaging through applying a Discrete Fourier Transform algorithm to determine changes to the wave speed as Carbon Black particles are introduced. As a result, due to the addition of Carbon Black as an auxiliary fuel source, detonations were formed instead of deflagrations in operating conditions where one would expect deflagrations at the same hydrogen-air equivalence ratios without Carbon Black addition. The detonation formation provides evidence that the coal particles are reacting within the detonation wave in a large enough capacity to support a detonation wave within the annulus. Furthermore, the wave speed is shown to increase with the additional of carbon particles. At a constant global equivalence ratio, the detonation wave velocities were found to decrease with hydrogen's incremental replacement with coal particles. Whereby, through a theoretical comparison of the heat of combustion as computed from the experimentally derived detonation wave velocities, a linear relationship of the two was shown to exist. Therefore, the heat of combustion has the potential to describe an operational limit to sustaining a detonation wave.

42 ENGINEERING↗
DOE OSTI2021

Prediction of Solution Behavior via Calorimetric Measurements Allows for Detailed Elucidation of Polyoxometalate Transformation

The solution behavior of a polyoxometalate cluster, LiNa-U 24 Pp 12 (Li 24 Na 24 [(UO 2 O 2 ) 24 (P 2 O 7 ) 12 ]) that consists of 24 uranyl ions, peroxide groups, and 12 pyrophosphate linkers, was successfully predicted based on new thermodynamic results using a calorimetric method recently described for uranyl peroxide nanoclusters (UPCs), molybdenum blues, and molybdenum browns. The breakdown of LiNa-U 24 Pp 12 and formation of U 24 (Li 24 [UO 2 O 2 OH] 24 ) was monitored in situ via Raman spectroscopy using a custom heating apparatus. A combination of analytical techniques confirmed the simultaneous existence of U 24 Pp 12 and U 24 midway through the conversion process and U 24 as the single end product. The application of a molecular weight filter resulted in a complete and successful separation of UPCs from solution and, in conjunction with DOSY results, confirmed the presence of large intermediate cluster building blocks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Molecular Origins of the Barriers to Proton Transport in Acidic Aqueous Solutions

The self-consistent iterative multistate empirical valence bond (SCI-MS-EVB) method is used to analyze the structure, thermodynamics, and dynamics of hydrochloric acid solutions. The reorientation timescales of irreversible proton transport are elucidated by simulating 0.43 M, 0.85 M, 1.68 M, and 3.26 M HCl solutions at 270 K, 285 K, 300 K, 315 K, and 330 K. The results indicate increased counter-ion pairing with increasing concentration, which manifests itself via a reduced hydronium oxygen–chloride (O*-Cl) structuring in the radial distribution functions. Increasing ionic concentration also reduces the diffusion of the hydrated excess protons, principally by reducing the contribution of the Grotthuss proton hopping (shuttling) mechanism to the overall diffusion process. In agreement with prior experimental findings, a decrease in the activation energy of reorientation timescales was also observed, which is explicitly explained using activated rate theory and an energy-entropy decomposition of the state-averaged radial distribution functions. Furthermore, these results provide atomistic verification of suggestions from recent two-dimensional infrared spectroscopy experiments that chloride anions (as opposed to hydrated excess protons) create entropic barriers to proton transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Robust increase in South Asian monsoon rainfall under warming driven by extratropical clouds and ocean

The responses of South Asian Monsoon (SAM) circulation under global warming are known to be highly uncertain, leading to the wide spread of SAM rainfall projections among models. Here, we show that the uncertain SAM circulation in Coupled Model Intercomparison Project Phase 6 models consists of two robust components that partly offset each other: a weakening component linked to a global thermodynamic constraint and a northward shift component understood through a regional 2D energetic perspective. We further attribute the robust northward shift of SAM circulation to positive cloud feedback over the Eurasia Continent and heat uptake in the Southern Ocean. A set of climate model simulations supports the finding that SAM rainfall increase is primarily due to the northward shift of circulation driven by extratropical processes. This energetic perspective opens new avenues for predicting monsoon rainfall by connecting circulation changes to radiative forcing, feedbacks, and ocean heat uptake.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Modeling the electrical double layer at solid-state electrochemical interfaces

Models of the electrical double layer (EDL) at electrode/liquid-electrolyte interfaces no longer hold for all-solid-state electrochemistry. Here we show a more general model for the EDL at a solid-state electrochemical interface based on the Poisson–Fermi–Dirac equation. By combining this model with density functional theory predictions, the interconnected electronic and ionic degrees of freedom in all-solid-state batteries, including the electronic band bending and defect concentration variation in the space-charge layer, are captured self-consistently. Along with a general mathematical solution, the EDL structure is presented in various materials that are thermodynamically stable in contact with a lithium metal anode: the solid electrolyte Li 7 La 3 Zr 2 O 12 (LLZO) and the solid interlayer materials LiF, Li 2 O and Li 2 CO 3 . Here, the model further allows design of the optimum interlayer thicknesses to minimize the electrostatic barrier for lithium ion transport at relevant solid-state battery interfaces.

25 ENERGY STORAGE↗
DOE OSTI2025

Observation of Unprecedented Fractional Magnetization Plateaus in a New Shastry-Sutherland Ising Compound

Geometrically frustrated magnetic systems, such as those based on the Shastry-Sutherland lattice (SSL), offer a rich playground for exploring unconventional magnetic states. The delicate balance between competing interactions in these systems leads to the emergence of novel phases. We present the characterization of Er 2 ⁢Be 2⁢ GeO 7 , an SSL compound with Er 3+ ions forming orthogonal dimers separated by nonmagnetic layers whose structure is invariant under the 𝑃⁢$\bar{4}21$⁢𝑚 space group. Neutron scattering reveals an antiferromagnetic dimer structure at zero field, typical of Ising spins on that lattice and consistent with the anisotropic magnetization observed. However, magnetization measurements exhibit fractional plateaus at 1/4 and 1/2 of saturation, in contrast to the expected 1/3 plateau of the SSL Ising model. By comparing the energy of candidate states with ground-state lower bounds we show that this behavior requires spatially anisotropic interactions, leading to an anisotropic Shastry-Sutherland Ising model symmetric under the 𝐶⁢𝑚⁢𝑚⁢2 space group. This anisotropy is consistent with the small orthorhombic distortion observed with single-crystal neutron diffraction. The other properties, including thermodynamics, which have been investigated theoretically using tensor networks, point to small residual interactions, potentially due to further couplings and quantum fluctuations. This study highlights Er 2 ⁢Be 2 ⁢GeO 7 as a promising platform for investigating exotic magnetic phenomena.

Yadav, Lalit [Duke University, Durham, NC (United ↗
DOE OSTI2020

Electronic and Vibrational Contributions to the Bulk Stabilities of Trivalent 3d Transition Metal Oxyhydroxides from Electronic Structure Calculations

Any understanding of the thermodynamic driving forces active in the global cycles of earth abundant metals requires attendant knowledge of the bulk mineral free energies. This can be particularly challenging for metal (oxyhydr)oxide phases which at low temperature aqueous conditions can be produced at a variety of hydration states and particle sizes. Here we present high quality density functional theory calculations of the relative bulk stabilities of a geologically important family of polymorphic trivalent metal oxyhydroxide crystal structures MOOH (M=Cr, Al, Fe). To improve understanding of their relative stabilities we permute the metal cations into the various possible canonical structure types including goethite (isostructural to diaspore and bracewillite), lepidocrocite (isosctructural to boehmite), guyaniite, and grimaldiite. Electronic, vibrational, and magnetic contributions are considered on a consistent theoretical footing using hybrid functionals and converged plane-wave bases. We find the considered bulk structures to be thermodynamically degenerate within the precision of the available theory, despite the fact that the vibrational and magnetic contributions vary substantially among the cations. A significant conclusion is that Cr3+ supports strong hydrogen bonds, low O-H stretching frequencies and facile H movement relative to the other cations.

Prange, Micah P.↗
DOE OSTI2022

Coupling Warm Rain With an Eddy Diffusivity/Mass Flux Parameterization: 1. Model Description and Validation

A new version of the stochastic multiplume Jet Propulsion Laboratory Eddy-Diffusivity/Mass-Flux (JPL-EDMF) parameterization which consistently couples the simplified Khairoutdinov and Kogan (2000), https://doi.org/10.1175/1520-0493(2000)128<0229:ANCPPI>2.0.CO;2, warm phase cloud microphysical parameterization with the parameterization of cloud macrophysical and subgrid scale dynamical processes is described. The new parameterization combines the EDMF approach with an assumed shape of a joint probability density function of thermodynamic and kinematic variables which provide the basis for the computation of all parameterized processes. As far as we are aware this is the first attempt to consistently couple all of these parameterized processes in the EDMF framework. This paper is part one of a two paper series. Here, the JPL-EDMF parameterization is described and benchmark simulations of precipitating stratocumulus and cumulus convection are performed in a single-column-model framework. The parameterization results compare favorably to the reference large-eddy-simulation results. In the second part (Smalley et al., 2022, https://doi.org/10.1029/2021MS002729) the JPL-EDMF parameterization is validated for a wide range of observation-based scenarios covering the continuous transition from subtropical stratocumulus to cumulus convection derived from global reanalysis, and parameterization uncertainties are studied in detail.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Modeling regional precipitation over the Indus River basin of Pakistan using statistical downscaling

Complex processes govern spatiotemporal distribution of precipitation within the high-mountainous headwater regions (commonly known as the upper Indus basin (UIB)), of the Indus River basin of Pakistan. Reliable precipitation simulations particularly over the UIB present a major scientific challenge due to regional complexity and inadequate observational coverage. Here, we present a statistical downscaling approach to model observed precipitation of the entire Indus basin, with a focus on UIB within available data constraints. Taking advantage of recent high altitude (HA) observatories, we perform precipitation regionalization using K-means cluster analysis to demonstrate effectiveness of low-altitude stations to provide useful precipitation inferences over more uncertain and hydrologically important HA of the UIB. We further employ generalized linear models (GLM) with gamma and Tweedie distributions to identify major dynamic and thermodynamic drivers from a reanalysis dataset within a robust cross-validation framework that explain observed spatiotemporal precipitation patterns across the Indus basin. Final statistical models demonstrate higher predictability to resolve precipitation variability over wetter southern Himalayans and different lower Indus regions, by mainly using different dynamic predictors. The modeling framework also shows an adequate performance over more complex and uncertain trans-Himalayans and the northwestern regions of the UIB, particularly during the seasons dominated by the westerly circulations. However, the cryosphere-dominated trans-Himalayan regions, which largely govern the basin hydrology, require relatively complex models that contain dynamic and thermodynamic circulations. Furthermore, we also analyzed relevant atmospheric circulations during precipitation anomalies over the UIB, to evaluate physical consistency of the statistical models, as an additional measure of reliability. Overall, our results suggest that such circulation-based statistical downscaling has the potential to improve our understanding towards distinct features of the regional-scale precipitation across the upper and lower Indus basin. Additionally, such understanding should help to assess the response of this complex, data-scarce, and climate-sensitive river basin amid future climatic changes, to serve communal and scientific interests.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Formic Acid: A Hydrogen-Bonding Cocatalyst for Formate Decomposition

Hydrogen bonding accelerates many catalytic reactions by orienting intermediates, stabilizing transition states, and even opening reaction pathways. However, most mechanistic studies regarding the decomposition of formic acid (FA), a promising hydrogen storage material, neglect hydrogen-bonding interactions even though FA is a strong hydrogen-bond donor and acceptor. Here, we probe the formation of bimolecular hydrogen-bonded complexes between FA and formate (FA–HCOO complexes) adsorbed on metal surfaces and how these complexes affect HCOO* decomposition. Using first-principles density functional theory (DFT) calculations on 12 close-packed (111)/(0001) and 8 open (100) surfaces of 12 transition metals—Ag, Au, Co, Cu, Ir, Ni, Os, Re, Pd, Pt, Rh, and Ru, we—show that FA–HCOO complexes are generally thermodynamically stable, even at elevated temperatures and pressures. We then illustrate that these complexes produce infrared spectroscopic signatures consistent with as yet unassigned experimental peaks. We last demonstrate that by stabilizing the dangling bond of monodentate HCOO*, these complexes significantly lower the barriers for rotation of HCOO* from a bidentate to a monodentate configuration, the rate-limiting step for HCOO* decomposition on many surfaces. FA thus acts as a cocatalyst for HCOO* decomposition. Our results may guide the community toward improved catalysts for reactions involving HCOO* such as FA decomposition, methanol steam reforming, and the water gas shift reaction. More broadly, our work highlights the ability of hydrogen bonding to modify the adsorbed structures of intermediates and lower the barriers for their reaction on heterogeneous catalysts. Lastly, this phenomenon can be relevant for other reactions involving ammonia, alcohols, and carboxylic acids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Design Criteria of Dilute Ether Electrolytes toward Reversible and Fast Intercalation Chemistry of Graphite Anode in Li-Ion Batteries

To date, dilute ether electrolytes have been believed to be incompatible with graphite in Li-ion batteries due to the detrimental solvent cointercalation and graphite exfoliation. Here, in this paper, we provide design criteria of dilute ether electrolytes for a reversible graphite anode based on tailoring the solvation structures and thermodynamic properties. We clarify that ether solvents can support graphite reversibly by modulating the anion. Our redesigned electrolyte consisting of a single-solvent 1,3-dioxolane (DOL) and 1 M single-salt lithium bis(fluorosulfonyl)imide (LiFSI) shows weakened Li-solvent interaction and results in an inorganic-rich solid-electrolyte interphase. Consequently, we achieved ~99.9% Coulombic efficiency with >96% capacity retention (~350 mAh/g) after 300 cycles at C/5 using natural graphite. The weakly solvated electrolyte maintains desirable transport properties, enabling better rate capability than carbonate electrolytes with an areal capacity of 2–4 mAh/cm 2 . We have demonstrated the potential of dilute ether electrolytes for facile desolvation-based intercalation chemistry in graphite, creating a viable path toward fast-charge Li batteries.

25 ENERGY STORAGE↗
DOE OSTI2023

Interface atomic structures in a cadmium arsenide/III–V semiconductor heterostructure

The interface atomic structure between an epitaxial thin film of the prototype topological semimetal cadmium arsenide (Cd 3 As 2 ) and a III–V semiconductor layer is investigated using high-angle annular dark-field imaging in an aberration-corrected scanning transmission electron microscope. We find that the interface unit cell adopts a defined stoichiometry that is CdSb-like, which is achieved through the insertion of periodically arranged Cd vacancies in the terminating Cd-plane of Cd 3 As 2 . This interface stoichiometry is consistent with the Sb-termination of the III–V layer and the fact that CdSb is the thermodynamically stable phase in the Cd–Sb binary system. We find at least two distinct alignments of the film with respect to the buffer layer, which are characterized by a $\frac{1}{4}$ <100> Cd 3 As 2 shift parallel to the interface. We show that steps of half unit cell height in the III–V layer can produce these distinct interface structures.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Dirac spectral density in N f = 2 + 1 QCD at T = 230 MeV

We compute the renormalized Dirac spectral density in N f = 2 + 1 QCD at physical quark masses, temperature T = 230 MeV , and system size L s = 3.4 fm . To that end, we perform a pointwise continuum limit of the staggered density in lattice QCD with staggered quarks. We find, for the first time, that a clear infrared structure (IR peak) emerges in the density of the Dirac operator describing dynamical quarks. We also provide numerical evidence that a component of this peak, which becomes dominant in the thermodynamic limit, is due to a nontrivial accumulation of near-zero modes. Features of this structure are consistent with those previously attributed to the recently proposed IR phase of thermal QCD. Our results (i) provide the only complete first-principles evidence that these IR features exist and are physical, (ii) improve the upper bound for IR phase transition temperature T IR so that the new window is 200 < T IR < 230 MeV , and (iii) are consistent with the nonrestoration of anomalous U A ( 1 ) symmetry (chiral limit) below T = 230 MeV . Published by the American Physical Society 2024

Alexandru, Andrei (ORCID:0000000345471554)↗
DOE OSTI2024

Salt Hydrate Eutectic Thermal Energy Storage for Building Thermal Regulation (Final Technical Report)

Thermal energy storage is anticipated to play an important role in developing the power grid of the future - a power grid that meets increasing demands of users, is resistant to disruptions, but also allows for greater penetration of renewable resources. Specifically, thermal energy storage materials can be integrated into HVAC systems and building envelopes, where they can be used to shift power demands for building climate control from periods of peak demand to periods of low demand. Phase change materials (PCMs) are compelling as low-cost, high energy density thermal energy storage materials for building thermal management. However, there is a lack of high performance low-cost PCMs within the specific temperature ranges which would most effectively allow for power load shifting. Inorganic salt hydrates represent a promising class of PCMs, but their inherent limitations cause them to be currently unavailable for reliable building applications. The overarching goals of this research effort are to: 1) Discover low-cost, high volumetric density salt hydrate eutectic PCMs to store low-quality heat (10 to 40 °C); 2) Introduce a high thermal conductivity matrix to reduce the time constant for energy storage to ~0.1 to 1 hr, incorporate nucleation catalysts to decrease undercooling, and utilize microencapsulation and shape stabilization approaches, minimizing moisture loss/gain, mitigating phase separation, and maintaining stable melting behavior over the lifetime of the compounds; 3) Evaluate the impact these systems have on peak load shifting, and the potential for overall energy savings under different climatic scenarios and building configurations. These goals will be achieved by an integrated research program consisting of six cohesive research subtasks: 1) Materials discovery of eutectic salt hydrate PCMs by using computationally predicted thermodynamic equilibria, coupled with high-throughput experimental validation, 2) Rapid experimental screening of nucleation catalysts identified through robust computational databases, 3) Embedding salt hydrate PCM into a low cost and scalable high conductivity matrix, 4) Microencapsulation of salt hydrate microspheres using hybrid inorganic-polymer microencapsulation approach, 5) Shape stabilization by thermoreversible salt hydrate salogels, and 6) Analysis of end-use using thermal simulations, and characterization of mock-up energy storage finished components.

25 ENERGY STORAGE↗
DOE OSTI2024

Thermodynamics and kinetics of interconversion between platinum nanoparticles and cations in zeolites

In metal-containing zeolites, sintering and redispersion processes exercise control over the identity of metal site structures, but the thermodynamic and kinetic factors that influence these molecular-level processes are not completely understood. Furthermore, we assess the ability of first principles informed free energy models (for supported and unsupported nanoparticles) and kinetic models integrating Ostwald ripening with atom trapping to describe the interconversion between Pt cations and nanoparticles encapsulated in chabazite (CHA) zeolites. Density functional theory-derived thermodynamic phase diagrams show that the interconversion between cations, favored in oxidizing environments, and particles, favored in reducing environments, is fully reversible within a wide range of their respective conditions (temperatures and pressures) for CHA and several other zeolite topologies. Kinetic Monte Carlo simulations of Pt redispersion are consistent with experimentally observed redispersion kinetics of encapsulated Pt nanoparticles in CHA zeolites, and model results suggest the zeolite host imparts additional stability for Pt nanoparticles. We envision our thermodynamic and kinetic models for Pt-CHA are also capable of describing nanoparticle and cation interconversions for other zeolite frameworks under similar conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Single-Electron Transfer Stabilizes Metastable Alane in a Bipyridine-Functionalized MOF Nanopore

Nanoconfinement of metastable hydrides such as alane (AlH 3 ) in functionalized porous frameworks represents a promising yet largely untapped strategy for high-capacity energy storage. Despite its potential, the underlying mechanisms responsible for the thermodynamic stabilization of metastable hydrides are poorly understood. Here, in this study, concepts from solution Lewis acid–base chemistry were applied to a crystalline metal–organic framework (MOF) to stabilize AlH 3 . The long-range order and synthetically versatile pore chemistry enabled us to reveal the intimate details of the hydride-host interactions. Powder X-ray diffraction (PXRD) of AlH 3 -infiltrated UiO-67bpy (Zr 6 O 4 (OH) 4 (bpydc) 6 ; bpydc 2– = 2,2′-bipyridine-5,5′-dicarboxylate) confirms that the framework maintains its crystallinity up to 250 °C and is stable under 70 MPa H 2 pressure. We find that thermodynamic stabilization arises from coupled single-electron and hydrogen-atom transfer from AlH 3 to bipyridine-containing linkers. Electron paramagnetic resonance detects a bipyridyl radical with an anisotropic g-tensor (g values of 2.00132, 2.00215, and 2.00275), consistent with Bader charge analysis predicting 0.82 e – transferred from the hydride to the MOF. These findings establish critical structure–property relationships underpinning AlH 3 stabilization and suggest general strategies for tuning the thermodynamics and kinetics of metastable materials.

Shivanna, Mohana [Sandia National Laboratories (SN↗
DOE OSTI2023

Self-consistent and detailed opacities from a non-equilibrium average-atom model

Modern density functional theory (DFT) is a powerful tool for accurately predicting self-consistent material properties such as equations of state, transport coefficients and opacities in high energy density plasmas, but it is generally restricted to conditions of local thermodynamic equilibrium (LTE) and produces only averaged electronic states instead of detailed configurations. We propose a simple modification to the bound-state occupation factor of a DFT-based average-atom model that captures essential non-LTE effects in plasmas—including autoionization and dielectronic recombination—thus extending DFT-based models to new regimes. Finally, we then expand the self-consistent electronic orbitals of the non-LTE DFT-AA model to generate multi-configuration electronic structure and detailed opacity spectra. This article is part of the theme issue ‘Dynamic and transient processes in warm dense matter’.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Model Intercomparison Study to Investigate Mixing Characteristics in Non‐Precipitating Stratocumulus Clouds

Recent aircraft observations of marine stratocumulus clouds consistently showed that cloud microphysical relationships vary with altitude, indicating inhomogeneous mixing characteristics near cloud top and homogeneous mixing characteristics in mid-levels of clouds. Here, we conduct model intercomparison of an idealized, non-precipitating stratocumulus cloud to evaluate model consistency and examine whether simulations can reproduce the observed mixing characteristics. The results show that eleven large-eddy simulations with various dynamics and microphysics schemes show good agreement on the thermodynamical, microphysical, and dynamical properties of the stratocumulus-topped boundary layer in a steady state. The inter-model spread in steady-state liquid water path is significantly reduced compared to previous model intercomparison studies. This improvement might be due to better models and more consistent initial conditions than those used decades ago. In addition, most simulations, including a low-dimensional simulation, capture inhomogeneous mixing characteristics near the cloud top and homogeneous mixing characteristics inside the cloud. Moreover, simulations using Lagrangian microphysics schemes agree better with the observed mixing characteristics compared with those using the bin microphysics schemes. Since most simulations do not fully resolve the entrainment process, the apparent mixing characteristics arise from the variations in the resolved cloud properties. Our results support the vertical circulation mixing hypothesis, which suggests that homogeneous mixing characteristics in mid-levels of clouds are due to the vertical circulation of entrainment-affected and diluted parcels from the cloud top moved to lower levels.

54 ENVIRONMENTAL SCIENCES↗