Search NASA⌕ Search

SEARCH · Search NASA

Results for “eddying model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Three-Dimensional, Analytical Wind Turbine Wake Model: Flow Acceleration, Empirical Correlations, and Continuity

A new, three-dimensional, analytical, steady-state wake model is presented that includes local flow acceleration near the rotor, improving the wake description compared to existing models. Wake structures such as the momentum deficit and regions of accelerated flow are concisely described with compound and normal Gaussian functions. Large-eddy simulations (LES) are used as training data to develop the model using two in-line turbines under various inflow conditions parameterized by hub-height wind speed and turbulence intensity. Mass conservation is considered by fixing two components of the wake velocity model and optimizing the third to best satisfy continuity; after which, the model performs comparably if not better than existing work with regards to both relative error and mass consistency. The final model demonstrates a high degree of flexibility, making use of empirical correlations to scale across different inflow conditions. The inclusion of these effects is capable of revealing unused opportunities for enhanced power generation by aligning wake trajectories with these regions of accelerated flow.

modeling↗
DOE OSTI2020

Global evidence for the acclimation of ecosystem photosynthesis to light

Photosynthesis responds quickly to changes in light, increasing with incoming photosynthetic photon flux density (PPFD) until the leaves become light saturated. This instantaneous response to PPFD, which is widely studied and incorporated into models of photosynthesis, is overlaid on non-instantaneous photosynthetic changes resulting from the acclimation of plants to average PPFD over intermediate timescales of a week to months ($\overline{PPFD}$). Such photosynthetic light acclimation is not typically incorporated into models, due to the lack of observational constraints. In this work, we use eddy covariance observations from globally distributed and automated sensor networks, along with photosynthesis estimates from nine terrestrial biosphere models (TBMs), to quantify and assess photosynthetic acclimation to light in natural environments. We also use recent theoretical developments to incorporate light acclimation in a TBM. Our results show widespread light acclimation of ecosystem photosynthesis. On average, a 1 μmol m -2 s -1 increase in $\overline{PPFD}_{10}$ (ten-day average PPFD) leads to a 0.031 ± 0.013 μmol C m -2 s -1 increase in the maximum photosynthetic assimilation rate (A max ), with croplands having stronger acclimation rates than grasslands and forests. Our analysis shows that the TBMs examined either neglect or substantially underestimate light acclimation. By updating a TBM to include photosynthetic acclimation, successfully reproducing the $\overline{PPFD}_{10}$-A max relationship, we provide a robust method for the incorporation of photosynthetic light acclimation in future models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Feasibility of full-core pin resolved CFD simulations of small modular reactor with momentum sources

Complex flow structure interactions and heat transfer processes take place in nuclear reactor cores. Given the extreme pressure/temperature and radioactive conditions inside the core, numerical simulations offer an attractive and sometimes more feasible approach to study the related flow and heat transfer phenomena in addition to the experiments. Under the Exascale Computing Project, the full-core simulation of a small modular reactor (SMR) has been pursued coupling Computational Fluid Dynamics (CFD) and neutronics. A key aspect of the modeling of SMR fuel assemblies is the presence of spacer grids and the mixing promoted by mixing vanes or the equivalent. A reduced order methodology is adopted based on momentum sources to mimic the mixing of the vanes. The momentum sources have been carefully calibrated with detailed Large Eddy Simulations (LES) of spacer grids performed with Nek5000. Modeling the spacer grid and mixing vanes (SGMV) effect without body-fitted computational grid avoids the excessive costs in resolving the local geometric details, and thus supports the simulation to be scaled up to the full core. Besides the progress on momentum source modeling, this paper also features the first full-core pin resolved CFD simulation ever performed to the authors' knowledge. This represents a significant advancement in capability for the CFD of nuclear reactors, which will hopefully serve as an inspiration for further integrating high-fidelity numerical simulations in actual engineering designs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Effects of Wind Shear and Aerosol Conditions on the Organization of Precipitating Marine Stratocumulus Clouds

Here this study examines how wind shear affects precipitating marine stratocumulus clouds under different cloud droplet number concentrations (N d ). We performed a series of large eddy simulations (LES) of nocturnal marine stratocumulus clouds using Cloud Model 1 (CM1). The simulations show that N d is the dominant factor for cloud cellular organization in this cloud regime rather than wind shear. Low N d characterizes the open cellular structure with a high in-cloud liquid water path (LWP). When wind shear is increased, the cloud fraction tends to decrease along with LWP, suggesting the cloud top region is significantly influenced by the entrainment and mixing of dry air from the free troposphere. We also examine cold pools in open and closed cellular clouds. Open-cell clouds produce larger and deeper cold pools compared to closed-cell clouds. Interestingly, cold pools can exist without surface precipitation and are produced by evaporation of light precipitation (drizzle) below the cloud base with weak downdrafts. The evaporation of raindrops and drizzle play an important role in initiating new convection, particularly where colliding outflows occur downstream of the cloud. This secondary convection contributes to the development and maintenance of the cloud cellular organization and formation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Two-Phase Turbulent Kinetic Energy Budget Computation in Co-Current Taylor Bubble Flow

Simulation of two-phase flows is relevant for reactor design and safety at normal operation or during accident scenarios. Often, the two-phase flow is in a regime in which slugs are formed or where the flow stratifies. Modeling such situations using standard single-phase Reynolds-averaged Navier-Stokes (RANS) turbulence models fails due to an overestimation of the eddy viscosity at the resolved two-phase interface. To solve this, an ad hoc turbulence damping term has been proposed in the literature that reduces the turbulence production locally at a two-phase interface, analogously to turbulence wall functions. However, this approach must be tailored to the specific setting and does not consider physical contributions such as surface tension or flow topology. Therefore, the problem of two-phase interfacial turbulence must be studied more in-depth. In this work, we consider co-current turbulent Taylor bubble flow using high-fidelity numerical simulation. The Basilisk code is used to simulate a Taylor bubble rising in a vertical pipe. By simulating the bubble in a moving frame of reference, we may study the turbulent kinetic energy (TKE) budgets ahead of the bubble, in its wake, and across the interface. The implementation of the TKE budget computation and the underlying averaging techniques are validated for the single-phase region ahead of the Taylor bubble using reference direct numerical simulation data. The analysis of the TKE budgets in the setting of Taylor bubble flow allows for the study of how turbulence behaves due to the presence of a two-phase interface and, in turn, supports the improvement of two-phase RANS models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A Macroalgal Cultivation Modeling System (MACMODS): Evaluating the Role of Physical-Biological Coupling on Nutrients and Farm Yield

Offshore aquaculture has the potential to expand the macroalgal industry. However, moving into deeper waters requires suspended structures that will present novel farm-environment interactions. Here, we present a computational modeling framework, the Macroalgal Cultivation Modeling System (MACMODS), to explore within-farm modifications to light, seawater flow, and nutrient fields across time and space scales relevant to macroalgae. A regional ocean model informs the site-specific setting, the Santa Barbara Channel in the Southern California Bight. A fine-scale hydrodynamic model predicts modified flows and turbulent mixing within the farm. A spatially resolved macroalgal growth model, parameterized for giant kelp, Macrocystis pyrifera , predicts kelp biomass. Key findings from model integration are that regional ocean conditions set overall farm performance, while fine-scale within-farm circulation and nutrient delivery are important to resolve variation in within-farm macroalgal performance. Therefore, we conclude that models resolving within-farm dynamics can provide benefit to farmers with insight on how farm design and regional ocean conditions interact to influence overall yield. Here, the presence of repeating longlines aligned with the mean current generate flow diversions around the farm as well as attached Langmuir circulations and increased turbulence intensity. These flow-induced phenomena lead to less biomass in the interior portion of the farm relative to the edges. We also find that there is an effluent “footprint” that extends as much as 20 km beyond the farm. In this regard, MACMODS can be used to not only evaluate farm design and cultivation practices that maximize yield but also explore interactions between the farm and ecosystem in order to minimize impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Turbulent Kinetic Energy Budgets over Gentle Topography Covered by Forests

Large-eddy simulations of flow over a “horizontally” uniform model forest are used to investigate the effects of gentle topography on the turbulent kinetic energy (TKE) budget within the canopy roughness sublayer. Despite significant differences between simulations using idealized sinusoidal topography and real topography of the Amazon forest, results indicate that the effects of topography are located predominantly in the upper canopy and above, and are mostly caused by mean advection of TKE. The “horizontally” averaged TKE budget from idealized and real gentle topographies are almost identical to that for flat terrain, including a clear inertial layer above the roughness sublayer in which shear production is balanced by local dissipation. At topography crests, where observational towers are usually located, mean vertical advection of TKE can be as important as horizontal advection. Here, we propose the use on an approximate TKE balance equation to estimate mean advection from single tower measurements, and introduce a new advection index that can be used as a proxy to quantify the importance of the topography on the TKE budget.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Subsetted model output for surface-related variables from a LASSO-CACTI WRF simulation

Surface-related model output from a Weather Research and Forecasting (WRF) simulation from the Large-Eddy Simulation (LES) Atmospheric Radiation Measurement (ARM) Symbiotic Simulation and Observation (LASSO) deep-convection scenario for the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign. The LASSO-CACTI simulations span grid spacings from 7.5 km to 100 m through the use of four nested domains, labeled d1 through d4. The simulations are of convection near the Sierras de Córdoba mountain range, roughly centered on the ARM Mobile Facility. More information can be found at https://www.arm.gov/capabilities/modeling/lasso. This version of the output is a collection of variables subsetted from the raw WRF output to reduce the file overhead for users not needing the full raw dataset. The subset files are in netCDF format with the variables in these files at the surface and sub-surface levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Evaluation of Spray and Combustion Models for Simulating Dilute Combustion in a Direct-Injection Spark-Ignition Engine

Dilute combustion in spark-ignition engines has the potential to improve thermal efficiency by mitigating knock and by reducing throttling and wall heat losses. However, ignition and combustion processes can become unstable for dilute operation due to a lowered laminar flame speed, resulting in excessive cycle-to-cycle variability (CCV) of the combustion process. To compensate for the slower combustion in less reactive mixtures, a modified intake port geometry can be employed to generate a strong tumble flow in the cylinder and elevate turbulence levels around the spark plug, thereby promoting a faster transition to turbulent deflagration. Consequently, optimizing combustion chamber geometry and operating strategy is crucial to maximizing the benefits of using dilute combustion with enhanced in-cylinder turbulence across a wide range of operating conditions. Computational fluid dynamics (CFD) simulations can be utilized for virtual engine optimization tasks, but this would require the models to be truly predictive regarding the impact of changes to the engine design and operational parameters.In this study, multicycle large-eddy simulations (LES) are performed for a direct-injection spark-ignition engine to investigate the model performance in predicting engine combustion characteristics with respect to changes in the intake configuration. A tumble plate that blocks the lower part of the intake port inlet is used to vary the tumble. A set of CFD models that have been recently developed are employed, which takes into account the drag of nonspherical droplets, flash-boiling behavior of liquid sprays, spray-wall interaction, surrogate formulation of a research-grade E10 gasoline, and fast chemical kinetic solvers. Simulation results are compared to experimental engine data in terms of cylinder pressure, apparent heat release rate, mass fraction burned timing, and flame images. It is found that LES employing the state-of-the-art CFD models are capable of properly predicting the spray processes and reproducing the measured mean cylinder pressure for the case with the tumble plate. On the other hand, the LES over-predicts the combustion rate during the early combustion stage and under-estimates the CCV, and these discrepancies become larger when the tumble plate is removed.

computational fluid dynamics simulation↗
DOE OSTI2022

Impact of Advection Schemes on Tracer Interrelationships in Large-Eddy Simulations of Deep Convection

This study investigates the preservation of tracer interrelationships during advection in large-eddy simulations of an idealized deep convective cloud, which is particularly relevant to chemistry, aerosol, and cloud microphysics models. Employing the Cloud Model 1, advection is represented using third-, fifth-, and seventh-order weighted essentially non-oscillatory schemes. As a simplified analogy for cloud hydrometeors and aerosols, several inert passive tracers following linear and nonlinear relationships are initialized after the cloud reaches ~6-km depth. Numerical mixing in the simulated turbulent convective clouds leads to significant deviations from the initial nonlinear relationships between tracers. In these simulations, a considerable fraction of the grid points where the tracers’ nonlinear relationships are altered from advection are classified as unrealistic (e.g., ~13% for the environmental tracers on average), including errors from range-preserving unmixing and overshooting. Errors in the sum of three tracers are also relatively large, ranging between ~1% and 16% for 5% of the grid points in and near the cloud. The magnitude of unrealistic mixing and errors in the sum of three tracers generally increase with the order of accuracy of the advection scheme. These results are consistent across model grid spacings ranging from 50 to 200 m, and across three different flow realizations for each combination of grid spacing and advection scheme tested. Tests employing a previously proposed scalar normalization procedure show substantially reduced errors in the sum of three tracers with a relatively small negative impact on other tracer relationships. This analysis, therefore, suggests efficacy of the normalization procedure when applied to turbulent three-dimensional cloud simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

An electric-circuit model on the inter-tape contact resistance and current sharing for REBCO cable and magnet applications

REBCO coated conductor has demonstrated high current capacity that can enable high-field magnets for high energy physics and fusion applications. However, quench protection is still one of the main challenges to be addressed for these applications. In addition, Ic and n value variations along the length of REBCO tapes exist in commercial production. The inter-tape contact resistance plays a key role to develop the self protection capability in cables and magnets by enabling current sharing and suppressing excessive eddy currents. Here we propose an electric-circuit model to describe the inter-tape contact resistance and its impact on the current sharing between REBCO tapes. We report the experiments on a 2-stacked tape REBCO cable with local Ic drop to validate the model. With the developed model, we study the upper limit of the contact resistance which allows current sharing between tapes. We also study the impact of variation in Ic and n values in tapes on the cable performance. Finally, our model is expected to provide useful insight into the current sharing and target values for inter-tape contact resistance in REBCO cables and magnets for various applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Quantifying carbon budget, crop yields and their responses to environmental variability using the ecosys model for U.S. Midwestern agroecosystems

As one of the major agricultural production areas in the world, the United States (U.S.) Midwest plays a vital role in the global food supply and agricultural ecosystem services. Although significant efforts have been made in modeling the carbon cycle dynamics over this area, large uncertainty still exists in the previous simulations in terms of reproducing individual components of the carbon cycle and their responses to environmental variability. Here we evaluated the performance of an advanced agroecosystem model, ecosys, in simulating carbon budgets over the U.S. Midwest, considering both the magnitude of carbon flux/yield and its response to environmental (climate and soil) variability. We conducted model simulations and evaluations at 7 cropland eddy-covariance sites as well as over 293 counties of Illinois, Indiana, and Iowa in the U.S. Midwest. The site-level simulations showed that ecosys captured both the magnitude and seasonal patterns of carbon fluxes (i.e., net ecosystem carbon exchange (NEE), ecosystem gross primary production (GPP), and ecosystem respiration (Reco)), leaf area index (LAI), and dynamic plant carbon allocation processes, with R2 equal to 0.92, 0.87, 0.87, and 0.78 for GPP, NEE, Reco, and LAI, respectively across all the sites compared with the observations. For regional scale simulations, ecosys reproduced the spatial distribution and interannual variability of corn and soybean yields with the constraints of observed yields and a new remotely sensed GPP product, with R2 of multi-year averaged simulated and observed yield equal 0.83 and 0.80 for corn and soybean, respectively. The simulated responses of carbon cycle dynamics to environmental variability were consistent with that from the empirical observations at both site and regional scales. Furthermore, our results demonstrated the applicability of ecosys in simulating the carbon cycle and soil carbon dynamics of the U.S. Midwestern agroecosystems under different climate and soil conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Interactions Between the Amazonian Rainforest and Cumuli Clouds: A Large‐Eddy Simulation, High‐Resolution ECMWF, and Observational Intercomparison Study

Abstract The explicit coupling at meter and second scales of vegetation's responses to the atmospheric‐boundary layer dynamics drives a dynamic heterogeneity that influences canopy‐top fluxes and cloud formation. Focusing on a representative day during the Amazonian dry season, we investigate the diurnal cycle of energy, moisture and carbon dioxide at the canopy top, and the transition from clear to cloudy conditions. To this end, we compare results from a large‐eddy simulation technique, a high‐resolution global weather model, and a complete observational data set collected during the GoAmazon14/15 campaign. The overall model‐observation comparisons of radiation and canopy‐top fluxes, turbulence, and cloud dynamics are very satisfactory, with all the modeled variables lying within the standard deviation of the monthly aggregated observations. Our analysis indicates that the timing of the change in the daylight carbon exchange, from a sink to a source, remains uncertain and is probably related to the stomata closure caused by the increase in vapor pressure deficit during the afternoon. We demonstrate quantitatively that heat and moisture transport from the subcloud layer into the cloud layer are misrepresented by the global model, yielding low values of specific humidity and thermal instability above the cloud base. Finally, the numerical simulations and observational data are adequate settings for benchmarking more comprehensive studies of plant responses, microphysics, and radiation.

Vilà‐Guerau de Arellano, J.↗
DOE OSTI2021

Time Evolution and Diurnal Variability of the Parametric Sensitivity of Turbine-Height Winds in the MYNN-EDMF Parameterization

The Mellor-Yamada-Nakanishi-Niino (MYNN) parameterization applied in the Weather Research and Forecasting (WRF) model has been augmented to include the Eddy-Diffusion Mass-Flux (EDMF) approach to better represent transport by boundary-layer eddies. This change includes the addition of new parameters associated with convective updrafts and boundary-layer clouds. These additions lead to new parametric sensitivities in the hub-height wind speed compared to simulations using the standard MYNN parameterization that were explored previously. This work builds on efforts focused on WRF’s MYNN parameterization by examining the sensitivity of hub-height wind speed to parameters in the MYNN-EDMF parameterization as a function of simulation lead time. Summer and winter periods were selected from the second Wind Forecast Improvement Project (WFIP2). Five sets of simulations were completed for each season, with lead times ranging from 2 to 6 days. The results show that the sensitivity to the new parameters associated with the EDMF scheme is generally small compared to other parameters, with the largest impact in convective boundary layers. The spread in the perturbed parameter ensembles was found to grow quickly over the first 8 to 19 hours in the summer simulations and 17 to 24 hours in the winter simulations with little change after that, regardless of the simulation length. A strong diurnal cycle in the parameter sensitivity was also found associated with changes in the atmospheric stability, as well as an increase in the sensitivity to the entrainment parameter used in the EDMF parameterization that is associated with increasing fractional areas covered by plumes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Multifidelity validation of digital surrogates using variable-density turbulent mixing models

In this study, ensembles of experimental data are presented and utilized to compare and validate two models used in the simulation of variable-density (Atwood = 0.22), compressible turbulent mixing. Though models of this kind (Reynolds averaged NavierStokes and large-eddy simulations) have been validated extensively with more canonical flows in previous studies, here the present approach offers novelty in the complexity of the geometry, the ensemble-based validation, and the uniformity of the computational framework on which the models are tested. Moreover, all experimental and computational tasks were completed by the authors which has led to a tightly coupled experimental configuration with its “digital twin.” The experimental divergent-shock-tube facility and its data acquisition methods are described and replicated in simulation space. A 2D Euler model which neglects the turbulent mixing at the interface is optimized to experimental data using a Gaussian process. This model then serves as the basis for both the 2D RANS and 3D LES studies that make comparisons to the mixing-layer data from the experiment. A relatively simple RANS model is shown to produce good agreement with experimental data only at late flow development times. The LES ensembles generally show good agreement with experimental data but display sensitivity to the characterization of initial conditions. Resolution-dependent behavior is also observed for certain higher-order statistics of interest. Overall, the LES model successfully captures the effects of divergent geometry, compressibility, and combined nonlinear instabilities inherent to the problem. The successful prediction of mixing width and its growth rate highlight the existence of three distinct regimes in the development of the instability, each with similarities to previously studied instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Combined Normal and Disruption, Electromagnetic Transient, Thermal, and Structural Analysis of COMPASS Upgrade

We report that COMPASS Upgrade is a replacement for the COMPASS device at the Institute of Plasma Physics (IPP) of the Czech Academy of Sciences. It is a cryogenic copper machine with a major radius of 0.89 m, with 5 T at the plasma centerline and 2-MA plasma current. It is intended to develop ITER relevant plasma densities and high power fluxes in the divertor region. The entire vacuum vessel is planned to reach 500 °C. The COMPASS-U TF coil is a bitter plate-like design with radial conductor plates that extend the width of the TF inner leg. Currents redistribute radially based on resistive and inductive effects. This necessitates an electromagnetic (EM) transient analysis of the coil current distribution. For copper coils like COMPASS Upgrade, the solution must be coupled with a thermal solution to properly include the temperature-dependent resistive effects. The EM transient solution includes the TF magnetic loads, and a stress pass on the results produces in-plane loads. The EM solution chosen requires inclusion of all the conducting and nonconducting regions surrounding the TF coil. With the TF current distribution solved, the addition of the poloidal coils and resulting background fields allows the determination of the out-of-plane loads. Integration of the Joule heating produced from the EM solution produces temperature distributions throughout the shot, which can be read in to a structural pass to include thermal stresses in the evaluation. This is true of both TF and PF temperature prediction with appropriate packing fractions. Currents in the PF coils, with computed background fields, produce Lorentz loads. Time transients of the PF currents produce startup eddy currents in the structures included in the model; inclusion of the vessel and other passive structures allows these to be included in a stress pass. The EM transient solution includes all the elements needed for a disruption analysis with prescribed motions and quench of the plasma. From a model originally intended to simulate current diffusion in the TF coil, the analysis can be refined to be used to predict the performance of most of the components of the tokamak or provide boundary conditions for more detailed submodels. In this article, the use of the current diffusion model for bounding assessments of the COMPASS Upgrade TF, PF, vessel, external structure, and disruption analyses will be presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Large Eddy Simulation of Externally Induced Ingress about an Axial Seal by Stator Vanes

Turbine inlet temperatures in advanced gas turbines could be as high as 2000 °C. To prevent ingress of this hot gas into the wheelspace between the stator and rotor disks, whose metals can only handle temperatures up to 850 °C, rim seals and sealing flows are used. This study examines the abilities of large eddy simulation (LES) based on the WALE subgrid model and Reynolds-averaged Navier–Stokes (RANS) based on the SST model in predicting ingress in a rotor–stator configuration with vanes but no blades, a configuration with experimental data for validation. Results were obtained for an operating condition, where the ratio of the external Reynolds number to the rotational Reynolds number is 0.538. At this operating condition, both LES and RANS were found to correctly predict the coefficient of pressure, C p , located downstream of the vanes and upstream of the seal, but only LES was able to correctly predict the sealing effectiveness. This shows C p by itself is inadequate in quantifying externally induced ingress. RANS was unable to predict the sealing effectiveness because it significantly under predicted the pressure drop in the hot gas path along the axial direction, especially about the seal region. This affected the pressure difference across the seal in the radial direction, which ultimately drives ingress.

42 ENGINEERING↗
DOE OSTI2022

Combined Effect of the Wegener–Bergeron–Findeisen Mechanism and Large Eddies on Microphysics of Mixed-Phase Stratiform Clouds

Abstract The process of glaciation in mixed-phase stratiform clouds was investigated by a novel Lagrangian–Eulerian model (LEM) in which thousands of adjoining Lagrangian parcels moved within a turbulent-like velocity field with statistical parameters typical of the Arctic boundary layer. We used detailed bin microphysics to describe the condensation/evaporation processes in each parcel, in which droplets, aerosols, and ice particles were described using size distributions of 500 mass bins. The model also calculated aerosol mass inside droplets and ice particles. Gravitational sedimentation of droplets and ice particles was also accounted for. Assuming that droplet freezing is the primary source of ice particles, the Arctic clouds observed in Indirect and Semi-Direct Aerosol Campaign (ISDAC) were successfully simulated. The model showed that at a low ice particle concentration typical of ISDAC, large vortices (eddies) led to a quasi-stationary regime, in which mixed-phase St existed for a long time. The large eddies controlled the water partitioning in the mixed-phase clouds. Droplets formed and grew in updrafts, typically reaching the cloud top, and evaporated in downdrafts. Ice particles grew in updrafts and downdrafts. The Wegener–Bergeron–Findeisen (WBF) mechanism was efficient in downdrafts and some parts of updrafts, depending on ice concentration and vertical velocity. At low ice concentrations, the effect of ice on the phase partitioning was negligible. In this regime, liquid droplets were found near the cloud top, whereas ice particles precipitated through the cloud base. When ice concentration exceeded about 10 L −1 , the WBF mechanism led to glaciation of almost the entire cloud, with the exception of narrow cloud regions associated with strong updrafts. At ice particle concentrations of a few tens per liter, the oscillatory regime took place due to the ice–liquid interaction. The microphysical structure of mixed-phase St forms as a combined effect of cloud dynamics (large eddies) and the WBF mechanism.

54 ENVIRONMENTAL SCIENCES↗