Large Eddy Simulation Model Research Articles

Multi-Physics Coupled Acoustic-Mechanics Analysis and Synergetic Optimization for a Twin-Fluid Atomization Nozzle

Fine particulate matter produced during the rapid industrialization over the past decades can cause significant harm to human health. Twin-fluid atomization technology is an effective means of controlling fine particulate matter pollution. In this paper, the influences of the main parameters on the droplet size, effective atomization range and sound pressure level (SPL) of a twin-fluid nozzle (TFN) are investigated, and in order to improve the atomization performance, a multi-objective synergetic optimization algorithm is presented. A multi-physics coupled acoustic-mechanics model based on the discrete phase model (DPM), large eddy simulation (LES) model, and Ffowcs Williams-Hawkings (FW-H) model is established, and the numerical simulation results of the multi-physics coupled acoustic-mechanics method are verified via experimental comparison. Based on the analysis of the multi-physics coupled acoustic-mechanics numerical simulation results, the effects of the water flow on the characteristics of the atomization flow distribution were obtained. A multi-physics coupled acoustic-mechanics numerical simulation result was employed to establish an orthogonal test database, and a multi-objective synergetic optimization algorithm was adopted to optimize the key parameters of the TFN. The optimal parameters are as follows: A gas flow of 0.94 m3/h, water flow of 0.0237 m3/h, orifice diameter of the self-excited vibrating cavity (SVC) of 1.19 mm, SVC orifice depth of 0.53 mm, distance between SVC and the outlet of nozzle of 5.11 mm, and a nozzle outlet diameter of 3.15 mm. The droplet particle size in the atomization flow field was significantly reduced, the spray distance improved by 71.56%, and the SPL data at each corresponding measurement point decreased by an average of 38.96%. The conclusions of this study offer a references for future TFN research.

Open boundary conditions for atmospheric large-eddy simulations and their implementation in DALES4.4

Abstract. Open boundary conditions were developed for atmospheric large-eddy simulation (LES) models and implemented into the Dutch Atmospheric Large-Eddy Simulation model. The implementation was tested in a “Big Brother”-like setup, in which the simulation with open boundary conditions was forced by an identical control simulation with periodic boundary conditions. The results show that the open boundary implementation has minimal influence on the solution. Both the mean state and the turbulent structures are close to the control simulation, and disturbances at the in- and outflow boundaries are negligible. To emulate a setup in which the LES is coupled to a coarser model, the influence of coarse boundary input was tested by smoothing the output of the periodic control simulation both temporally and spatially before feeding it as input to the simulation with open boundary conditions. When smoothing is applied over larger spatial and longer temporal scales, disturbances start to form at the inflow boundary and an area exists where turbulence needs to develop. Adding synthetic turbulence to the smoothed input reduces the size of this area and the magnitude of the disturbances.

Evaluation of flow and flame characteristics of turbulent bluff-body CH4-H2 flame using LES-FPV approach

ABSTRACT This study investigated the flow and flame structure of a bluff-body stabilised CH4-H2 flame using the large eddy simulation (LES) and the flamelet progress variable (FPV) model. The performance of three subgrid-scale (SGS) LES models, i.e. the standard Smagorinsky model (SSM), the wall-adapting local eddy viscosity (WALE) model, and the dynamic Smagorinsky model (DSM) in predicting turbulent and reactive flow characteristics was assessed. Parameters, namely y+ , Pope index (M), vorticity field, and resolved turbulent kinetic energy (kt ), were analysed to determine the predictive capability of the three models. The results demonstrated that DSM and WALE successfully resolved 80% of kt for chosen grid configurations, whereas SSM exhibited discrepancies. The overall trends captured by all three models matched well with the experimental trend with DSM and WALE outperforming SSM in predicting reacting flow fields in the recirculation zone, while SSM produced better predictions in the regions away from the recirculating zone.

Numerical prediction of mean flow and acoustic field of a supersonic impinging jet

The three-dimensional turbulent mean flow and acoustic field of a supersonic jet impinging on a solid plate is studied computationally using the general purpose CFD code Ansys Fluent. A pressure-based coupled solver formulation with the second order weighted central-upwind spatial discretization is applied to compute transient solutions. Cold and hot jet thermal conditions are considered. Mean flow characteristics are investigated by a steady-state modeling approach. Acoustic radiation of impingement tones is simulated using a transient time-domain formulation. The effects of turbulence in steady-state are modeled by the SST k-ω turbulence model. The Wall-Modeled Large-Eddy Simulation (WMLES) model is applied to compute transient solutions. The near-wall mesh on the impingement plate is fine enough to resolve the viscosity-affected near-wall region all the way to the laminar sublayer. Nozzle-to-plate distance is parameterized in the model for automatic re-generation of the mesh and results. Steady-state predictions of hover lift loss and mean jet velocity distributions are compared with experimental data, and favorable agreement is reported. The transient solution reproduces the mechanism of impingement tone generation by the interaction of large scale vortical structures with the impingement plate. The acoustic near-field is directly resolved by Computational Aeroacoustics (CAA) to accurately propagate impingement tone waves to near-field microphone locations. Calculated impingement tone frequencies and sound pressure levels agree with experimental values.

Adaptive detached eddy simulation of turbulent combustion with the subgrid dissipation concept

Detached eddy simulation has become a widely used method in eddy simulations due to its balance between cost and accuracy. The recently developed subgrid dissipation concept (SDC) combustion model [Liu et al., “On the subgrid dissipation concept for large eddy simulation of turbulent combustion,” Combust. Flame 258, 113099 (2023)] is found to be more reasonable and accurate than the conventional eddy dissipation concept model in large eddy simulation (LES). In this paper, the SDC model is adapted to the ℓ2-ω adaptive detached eddy simulation framework, named DES-SDC. The required key quantities, including the fine structure mass fraction and dissipation rate, are appropriately blended across Reynolds-averaged Navier–Stokes and LES regions. The DES-SDC approach is validated using premixed bluff body stabilized flame, non-premixed swirl flame, and premixed swirl flame with complex geometry. It is much more tolerant to coarse mesh resolution than pure LES, yet it preserves the capability of resolving the key unsteady feature critical for the combustion process, as it is designed to be. The DES-SDC approach is relatively insensitive to the grid resolution. The present research provides a promising approach for accurately simulating practical unsteady turbulent combustion problems at an affordable computational cost.

Characteristics of very-large-scale motions in natural rivers based on large eddy simulation

Exploring very-large-scale motions (VLSMs) in open channel flows (OCFs) is crucial for comprehensively understanding material transport and energy exchange. While previous studies focused on OCFs in ideal flumes with simple boundaries, this paper presents large eddy simulation (LES) results on the existence and scale characteristics of VLSMs in complex river morphologies with the Minjiang River in southwestern China as a case study. This research demonstrates that the OpenFOAM-based LES model is capable of accurately reconstructing the time-averaged flow field and providing instantaneous velocity data that capture VLSMs with sufficient resolution. Spectrum analysis of the streamwise fluctuating velocity shows that VLSMs are present in the straight segment of the river, with streamwise wavelengths approximately (16–22) times the water depth but are absent in the bend and confluence segments due to the suppression by secondary currents. This article validates the effectiveness of LES in examining the characteristics of VLSMs in natural rivers, thereby laying a good foundation for further studies on the impacts of such structures on sediment transport and pollutant dispersion.

Assessment on fast simulation of wind-driven pollutant dispersion around a street canyon with regime-switching Markov chain

High-density cities shorten commutes but also create proximity to air pollution sources. Sensors coupled with simulations can achieve real-time monitoring and risk assessment, and the recently emerging data-driven models provide opportunities that allow minimizing the computation burden of physical models. This study investigates the use of a time inhomogeneous Markov chain (CLF-MC) to fast simulation of wind-driven pollutant dispersion around a street canyon. The CLF-MC model is trained and tested using synthetic datasets obtained from large eddy simulations (LES). The LES model is validated using a benchmark wind tunnel test, which yields correlation coefficients above 0.95. The results show that the CLF-MC can faithfully reproduce the cumulative distribution of air pollution concentrations at the receptor site using the optimal model configuration identified in this study. For the hold-out test dataset, the release condition is a time-varying release strength, and the CLF-MC shows improved prediction results compared to the baseline models, i.e. the weighted average deviation is 15 % and 51 % smaller, respectively. In comparison to LES, the CLF-MC model is 4.8 × 106 times faster. These findings provide scientific references for fast and accurate wind-driven pollutant dispersion simulation using the probabilistic modelling approach.

Multi-scale simulation of typhoon wind field at building scale utilizing mesoscale model with nested large eddy simulation

Based on the mesoscale WRF (Weather Research Forecast) and LES (Large Eddy Simulation), an integrated numerical simulation framework aiming to resolve typhoon wind fields around urban building blocks has been developed. The holistic multi-scale simulation spans four scales: macroscale (∼1000 km), mesoscale (∼100 km), microscale (∼1 km), and building scale (0.1–100m). By interactively sharing the meteorological data among different scales, this simulation framework addressed the challenge of downscaling typhoon effects in urban environments. Typhoon weather reanalysis and urban wind field simulation were conducted by the WRF model with nested computational domains incorporating high-resolution topography data. For the downscaled LES model with urban buildings, the wind profiles produced by the WRF model were used at the interface with provided turbulent fluctuations in the flow. K11 building is a 270-m-tall building located in a densely built-up area of Hong Kong. Wind velocity and pressure fields surrounding the K11 building in Hong Kong during Typhoon Kammuri were obtained by the proposed framework, and the simulation results were validated by comparing with the meteorological data and full-scale measurements at the top of the K11 building. The interactions of the wind fields and building clusters resulting in vortex shedding, separation, and channeling effects have been resolved at high resolutions. Probabilistic distributions and higher-order statistics of wind pressure, i.e., skewness and kurtosis, were presented to highlight non-Gaussian characteristics of the simulated local wind pressure on the K11 building. Based on the simulated wind pressure, wind-induced vibration of the K11 building during the typhoon was analyzed and compared to the measurements.

Mixing Time Prediction in a Ladle Furnace

This paper presents a study on the effectiveness of two turbulence models, the large eddy simulation (LES) model and the k-ε turbulence model, in predicting mixing time within a ladle furnace using the computational fluid dynamics (CFD) technique. The CFD model was developed based on a downscaled water ladle from an industrial ladle. Corresponding experiments were conducted to provide insights into the flow field, which were used for the validation of CFD simulations. The correlation between the flow structure and turbulence kinetic energy in relation to mixing time was investigated. Flow field results indicated that both turbulence models aligned well with time-averaged velocity data from the experiments. However, the LES model not only offered a closer match in magnitude but also provided a more detailed representation of turbulence eddies. With respect to predicting mixing time, increased flow rates resulted in extended mixing times in both turbulence models. However, the LES model consistently projected longer mixing times due to its capability to capture a more intricate distribution of turbulence eddies.

Analyzing nitrogen dioxide to nitrogen oxide scaling factors for data-driven satellite-based emission estimation methods: A case study of Matimba/Medupi power stations in South Africa

In this paper, we propose improved nitrogen dioxide (NO2) to nitrogen oxide (NOx) scaling factors for several data-driven methods that are used for the estimation of NOx power plant emissions from satellite observations of NO2. The scaling factors are deduced from high-resolution simulations of power plant plumes with the MicroHH large-eddy simulation model with a simplified chemistry and then applied to Sentinel-5 Precursor (S5P) TROPOspheric Monitoring Instrument (TROPOMI) NO2 satellite observations over the Matimba/Medupi power stations in South Africa. We show that due to the non-linear chemistry the optimal NO2 to NOx scaling factors depend on both the method employed and the specific segments of the plume from which emission estimate is derived. The scaling factors derived from the MicroHH simulations in this study are substantially (more than 50%) higher than the typical values used in the literature with actual NO2 observations. The results highlight the challenge in appropriately accounting for the conversion from NO2 to NOx when estimating point source emissions from satellite NO2 observations.

Comparing the simulated influence of biomass burning plumes on low-level clouds over the southeastern Atlantic under varying smoke conditions

Abstract. Biomass burning plumes are frequently transported over the southeast Atlantic (SEA) stratocumulus deck during the southern African fire season (June–October). The plumes bring large amounts of absorbing aerosols and enhanced moisture, which can trigger a rich set of aerosol–cloud–radiation interactions with climatic consequences that are still poorly understood. We use large-eddy simulation (LES) to explore and disentangle the individual impacts of aerosols and moisture on the underlying stratocumulus clouds, the marine boundary layer (MBL) evolution, and the stratocumulus-to-cumulus transition (SCT) for three different meteorological situations over the southeast Atlantic during August 2017. For all three cases, our LES shows that the SCT is driven by increased sea surface temperatures and cloud-top entrainment as the air is advected towards the Equator. In the LES model, aerosol indirect effects, including impacts on drizzle production, have a small influence on the modeled cloud evolution and SCT, even when aerosol concentrations are lowered to background concentrations. In contrast, local semi-direct effects, i.e., aerosol absorption of solar radiation in the MBL, cause a reduction in cloud cover that can lead to a speed-up of the SCT, in particular during the daytime and during broken cloud conditions, especially in highly polluted situations. The largest impact on the radiative budget comes from aerosol impacts on cloud albedo: the plume with absorbing aerosols produces a total average 3 d of simulations. We find that the moisture accompanying the aerosol plume produces an additional cooling effect that is about as large as the total aerosol radiative effect. Overall, there is still a large uncertainty associated with the radiative and cloud evolution effects of biomass burning aerosols. A comparison between different models in a common framework, combined with constraints from in situ observations, could help to reduce the uncertainty.

A mass transfer-based LES modelling methodology for analyzing the movement of submarine sediment flows with extensive shear behavior

This study employs the large-eddy simulation (LES) method to investigate the interaction of submarine sediment flows with extensive shear behavior and ambient water. This method is validated with good accuracy by simulating the head velocity and evolution geometries of a submarine mud flow with non-Newtonian fluid characteristics and a submarine turbidity current with Newtonian fluid characteristics in inclined and regular lock-exchange flume experiments. This study finds that the violent mass transfer at the interface of the submarine sediment flow and seawater is caused by numerous eddies, which significantly alter the submarine sediment flow's geometry, resulting in varying degrees of undulation, including the deposition pattern of the submarine sediment flow tail. The acceleration zone, where the velocity of the submarine sediment flow increases significantly at this undulation, propels the sediment flow forward, supporting its long-distance transport. Furthermore, the turbulence and mass transport characteristics of submarine turbidity currents with low dynamic viscosity Newtonian fluid characteristics are stronger than those of submarine mud flows with high dynamic viscosity non-Newtonian fluid characteristics. Therefore, when submarine landslides develop into later stages, such as submarine turbidity currents with high velocity, large volume, and long run-out distance characteristics, more attention must be given to the mass transport process.

Challenges of constructing and selecting the “perfect” boundary conditions for the large-eddy simulation model PALM

Abstract. We present the process of and difficulties in acquiring the proper boundary conditions (BCs) for the state-of-the-art large-eddy simulation (LES)-based PALM model system. We use the mesoscale Weather Research and Forecasting (WRF) model as a source of inputs for the PALM preprocessor and investigate the influence of the mesoscale model on the performance of the PALM model. A total of 16 different WRF configurations were used as a proxy for a multi-model ensemble. We developed a technique for selecting suitable sets of BCs, performed PALM model simulations driven by these BCs, and investigated the consequences of selecting a sub-optimal WRF configuration. The procedure was tested for four episodes in different seasons of the year 2019, during which WRF and PALM outputs were evaluated against the atmospheric radiosounding observations. We show that the PALM model outputs are heavily dependent on the imposed BCs and have different responses at different times of the day and in different seasons. We demonstrate that the main driver of errors is the mesoscale model and that the PALM model is capable of attenuating but not fully correcting them. The PALM model attenuates the impact of errors in BCs in wind speed, while for the air temperature, PALM shows variable behavior with respect to driving conditions. This study stresses the importance of high-quality driving BCs and the complexity of the process of their construction and selection.

A-priori evaluation of data-driven models for large-eddy simulations in Rayleigh–Bénard convection

Natural convection is a commonly occurring heat-transfer problem in many industrial flows and its prediction with conventional large eddy simulations (LES) at higher Rayleigh numbers using progressively coarser grids leads to increasingly inaccurate estimates of important performance indicators, such as Nusselt number (Nu). Thus, to improve the heat transfer predictions, we utilize Gene Expression Programming (GEP) to develop sub-grid scale (SGS) stress and SGS heat-flux models simultaneously for LES. The models’ development is performed using reference direct numerical simulation data of a typical natural convection case, i.e. the Rayleigh–Bénard convection (RBC). The training frameworks are built by using alignment as the evaluation metrics of the basis functions when comparing to the Gaussian-filtered SGS stress and SGS heat-flux. The trained models in isotropic form, by utilizing the L2 norm of the grid cell ΔL2 as the length scale, demonstrate good performance in the bulk region, but less improved performance in the near wall region. It is shown, that for LES of wall-bounded flow, the GEP models in anisotropic form, i.e. using different grid length scale for the different spatial directions, are required to obtain generalized models suitable for different regions. Consequently, the a-priori results demonstrate a significant improvement in the prediction of both instantaneous and mean quantities for a wide range of filter widths.

On the Sensitivity of Large-Eddy Simulations of the Atmospheric Boundary Layer Coupled with Realistic Large-Scale Dynamics

Abstract We present a new ensemble of 36 numerical experiments aimed at comprehensively gauging the sensitivity of nested large-eddy simulations (LES) driven by large-scale dynamics. Specifically, we explore 36 multiscale configurations of the Weather Research and Forecasting (WRF) Model to simulate the boundary layer flow over the complex topography at the Perdigão field site, with five nested domains discretized at horizontal resolutions ranging from 11.25 km to 30 m. Each ensemble member has a unique combination of the following input factors: (i) large-scale initial and boundary conditions, (ii) subgrid turbulence modeling in the gray zone of turbulence, (iii) subgrid-scale (SGS) models in LES, and (iv) topography and land-cover datasets. We probe their relative importance for LES calculations of velocity, temperature, and moisture fields. Variance decomposition analysis unravels large sensitivities to topography and land-use datasets and very weak sensitivity to the LES SGS model. Discrepancies within ensemble members can be as large as 2.5 m s−1 for the time-averaged near-surface wind speed on the ridge and as large as 10 m s−1 without time averaging. At specific time points, a large fraction of this sensitivity can be explained by the different turbulence models in the gray zone domains. We implement a horizontal momentum and moisture budget routine in WRF to further elucidate the mechanisms behind the observed sensitivity, paving the way for an increased understanding of the tangible effects of the gray zone of turbulence problem. Significance Statement Several science and engineering applications, including wind turbine siting and operations, weather prediction, and downscaling of climate projections, call for high-resolution numerical simulations of the lowest part of the atmosphere. Recent studies have highlighted that such high-resolution simulations, coupled with large-scale models, are challenging and require several important assumptions. With a new set of numerical experiments, we evaluate and compare the significance of different assumptions and outstanding challenges in multiscale modeling (i.e., coupling large-scale models and high-resolution atmospheric simulations). The ultimate goal of this analysis is to put each individual assumption into the wider perspective of a realistic problem and quantify its relative importance compared to other important modeling choices.

Comparative analysis of different CFD turbulence models for a diesel pool combustion event in an underground mine: a case study

ABSTRACT Underground mine fires are dangerous, consume resources, and emit pollutants. Therefore, acquiring in-mine data for fire research is hazardous and expensive due to unsteady temperature changes and toxic gas emissions. Computational fluid dynamics (CFD) models are often developed to mimic airflows and combustion-associated heat transfer. However, research on the performance evaluation of the CFD turbulence models to mimic fire events is limited. This paper summarises the modelling results for a diesel pool fire event in an experimental underground mine that uses an exhaust ventilation system. A fire dynamics simulator (FDS) model was developed to show the spatio-temporal temperature evolution. Reynolds-Averaged Navier-Stokes and large eddy simulation (LES) models were developed using ANSYS Fluent to model species transportation in the transient-state combustion event. LES models showed the closest agreement with the experiments but required the most computing resources. Similar exercises will assist the mine operators in emergency preparedness and in training their safety teams.

MUST dome design based on dome seeing quantitative evaluation

ABSTRACT The Multiplexed Survey Telescope (MUST) project is led by Tsinghua University, which entrusted the European Industrial Engineering (EIE) GROUP with the design, manufacture, and assembly of the dome. Located on Saishiteng Mountain in Qinghai Province, MUST benefits from exceptional atmospheric seeing conditions etc. Local dome seeing may be comparable to atmospheric seeing and requires careful consideration. This research, based on numerical simulations, focuses on refining the dome structure and temperature regulation strategies to achieve optimal dome seeing. The existing simulations only consider nighttime dome seeing and overlook the impact of daytime dome heating on nighttime conditions. The Computational Fluid Dynamics (CFD) part mostly relies on a steady-state k − ε model, which cannot simulate transient processes or capture optical turbulence. In this study, a comprehensive multiphysics field coupling simulation was conducted, encompassing radiation heat transfer, fluid heat transfer, CFD, and ray tracing. Simulations include both daytime and nighttime scenarios, taking into account the daytime heating of the dome due to solar irradiation, as well as dome seeing under natural ventilation at night. The CFD utilizes the large eddy simulation model, enabling three-dimensional transient simulation and the simulation of optical turbulence. Ultimately, the broadening of the point spread function was statistically analysed after a certain integration time, facilitating a quantitative evaluation of dome seeing. This numerical simulation approach is closer to real world conditions, improving simulation accuracy and addressing the shortcomings of existing simulations. Some qualitative conclusions are consistent with practical engineering experience. In the end, the dome seeing was successfully regulated to 0.21 arcsec, meeting the observational requirements.

Model-based evaluation of cloud geometry and droplet size retrievals from two-dimensional polarized measurements of specMACS

Abstract. Cloud radiative properties play a significant role in radiation and energy budgets and are influenced by both the cloud top height and the particle size distribution. Both cloud top heights and particle size distributions can be derived from 2-D intensity and polarization measurements by the airborne spectrometer of the Munich Aerosol Cloud Scanner (specMACS). The cloud top heights are determined using a stereographic method (Kölling et al., 2019), and the particle size distributions are derived in terms of the cloud effective radius and the effective variance from multidirectional polarized measurements of the cloudbow (Pörtge et al., 2023). In this study, the accuracy of the two methods is evaluated using realistic 3-D radiative transfer simulations of specMACS measurements of a synthetic field of shallow cumulus clouds, and possible error sources are determined. The simulations are performed with the 3-D Monte Carlo radiative transport model MYSTIC (Mayer, 2009) using cloud data from highly resolved large-eddy simulations (LESs). Both retrieval methods are applied to the simulated data and compared to the respective properties of the underlying cloud field from the LESs. Moreover, the influence of the cloud development on both methods is evaluated by applying the algorithms to idealized simulated data where the clouds did not change during the simulated overflight of 1 min over the cloud field. For the cloud top height retrieval, an absolute mean difference of less than 70 m with a standard deviation of about 130 m compared to the expected heights from the model is found. The elimination of the cloud development as a possible error source results in mean differences of (46±140) m. For the effective radius, an absolute average difference of about (-0.2±1.30) µm from the expected effective radius from the LES model input is derived for the realistic simulation and (-0.03±1.28) µm for the simulation without cloud development. The difference between the effective variance derived from the cloudbow retrieval and the expected effective variance is (0.02±0.05) for both simulations. Additional studies concerning the correlations between larger errors in the effective radius or variance and the optical thickness of the observed clouds have revealed that low values in the optical thickness do not have an impact on the accuracy of the retrieval.

Analysis of the air stream flow parameters generated by the positive pressure ventilator—full scale experiment and CFD simulation

Positioning the positive pressure ventilator in front of the door opening affects the effectiveness of the rescue operation carried out during a fire. An important factor determining the effectiveness of the positive pressure ventilator is also the layout of the rooms within the gas exchange path and the obstacles present there. The purpose of this article is to assess the feasibility of using analyses such as large eddy simulation (LES) to verify the efficiency of mobile fans under simulation conditions, without the need for time-consuming experimentation (also for complex room volumes of buildings). The article presents a comparative analysis to assess the degree of convergence of flow parameters obtained during an experiment (in a multi-story building) and computational fluid dynamics (CFD) simulations. For volumetric flow rate, convergence was achieved at levels ranging from 0.4% (for 5 m) to 11.5% (1 m), and for pressure values, the differences achieved ranged from 0.6% (5 m) to 30.1% (4 m). This paper demonstrates that the LES model can be used to perform CFD simulations in the area of assessing the performance of a positive pressure ventilator. The article also describes a test methodology for determining the flow parameters of an air stream, which can be used to perform numerical simulations.

Above-cloud concentrations of cloud condensation nuclei help to sustain some Arctic low-level clouds

Abstract. Previous studies have found that low-level Arctic clouds often persist for long periods even in the face of very low surface cloud condensation nuclei (CCN) concentrations. Here, we investigate whether these conditions could occur due to continuous entrainment of aerosol particles from the free troposphere (FT). We use an idealized large eddy simulation (LES) modeling framework, where aerosol concentrations are low in the boundary layer (BL) but increased up to 50× in the free troposphere. We find that the tests with higher tropospheric aerosol concentrations simulated clouds, which persisted for longer and maintained higher liquid water paths (LWPs). This is due to direct entrainment of the tropospheric aerosol into the cloud layer, which results in a precipitation suppression from the increase in cloud droplet number and in stronger cloud-top radiative cooling, which causes stronger circulations maintaining the cloud in the absence of surface forcing. Together, these two responses result in a more well-mixed boundary layer with a top that remains in contact with the tropospheric aerosol reservoir and can maintain entrainment of those aerosol particles. The surface aerosol concentrations, however, remained low in all simulations. The free-tropospheric aerosol concentration necessary to maintain the clouds is consistent with concentrations that are frequently seen in observations.