Showing papers on "Convection published in 2014"

A scale and aerosol aware stochastic convective parameterization for weather and air quality modeling

Abstract: . A convective parameterization is described and evaluated that may be used in high resolution non-hydrostatic mesoscale models as well as in modeling system with unstructured varying grid resolutions and for convection aware simulations. This scheme is based on a stochastic approach originally implemented by Grell and Devenyi (2002). Two approaches are tested on resolutions ranging from 20 km to 5 km. One approach is based on spreading subsidence to neighboring grid points, the other one on a recently introduced method by Arakawa et al. (2011). Results from model intercomparisons, as well as verification with observations indicate that both the spreading of the subsidence and Arakawa's approach work well for the highest resolution runs. Because of its simplicity and its capability for an automatic smooth transition as the resolution is increased, Arakawa's approach may be preferred. Additionally, interactions with aerosols have been implemented through a cloud condensation nuclei (CCN) dependent autoconversion of cloud water to rain as well as an aerosol dependent evaporation of cloud drops. Initial tests with this newly implemented aerosol approach show plausible results with a decrease in predicted precipitation in some areas, caused by the changed autoconversion mechanism. This change also causes a significant increase of cloud water and ice detrainment near the cloud tops. Some areas also experience an increase of precipitation, most likely caused by strengthened downdrafts.

Meteorological conditions for the persistent severe fog and haze event over eastern China in January 2013

Abstract: In January 2013, a severe fog and haze event (FHE) of strong intensity, long duration, and extensive coverage occurred in eastern China. The present study investigates meteorological conditions for this FHE by diagnosing both its atmospheric background fields and daily evolution in January 2013. The results show that a weak East Asian winter monsoon existed in January 2013. Over eastern China, the anomalous southerly winds in the middle and lower troposphere are favorable for more water vapor transported to eastern China. An anomalous high at 500 hPa suppresses convection. The weakened surface winds are favorable for the fog and haze concentrating in eastern China. The reduction of the vertical shear of horizontal winds weakens the synoptic disturbances and vertical mixing of atmosphere. The anomalous inversion in near-surface increases the stability of surface air. All these meteorological background fields in January 2013 were conducive to the maintenance and development of fog and haze over eastern China. The diagnosis of the daily evolution of the FHE shows that the surface wind velocity and the vertical shear of horizontal winds in the middle and lower troposphere can exert dynamic effects on fog and haze. The larger (smaller) they are, the weaker (stronger) the fog and haze are. The thermodynamic effects include stratification instability in middle and lower troposphere and the inversion and dew-point deficit in near-surface. The larger (smaller) the stratification instability and the inversion are, the stronger (weaker) the fog and haze are. Meanwhile, the smaller (larger) the dewpoint deficit is, the stronger (weaker) the fog and haze are. Based on the meteorological factors, a multi-variate linear regression model is set up. The model results show that the dynamic and thermodynamic effects on the variance of the fog and haze evolution are almost the same. The contribution of the meteorological factors to the variance of the daily fog and haze evolution reaches 0.68, which explains more than 2/3 of the variance.

Evaluation of the convection‐resolving regional climate modeling approach in decade‐long simulations

Abstract: The uncertainties in current global and regional climate model integrations are partly related to the representation of clouds, moist convection, and complex topography, thus motivating the use of convection-resolving models. On climate time scales, convection-resolving methods have been used for process studies, but application to long-term scenario simulations has been very limited. Here we present a convection-resolving simulation for a 10 yearlong period (1998–2007) integrated with the Consortium for Small-Scale Modeling in Climate Mode model. Two one-way nested grids are used with horizontal resolutions of 2.2 km for a convection-resolving model (CRM2) on an extended Alpine domain (1100 km × 1100 km) and 12 km for a convection-parametrizing model (CPM12) covering Europe. CPM12 is driven by lateral boundary conditions from the ERA-Interim reanalysis. Validation is conducted against high-resolution surface data. The CRM2 model strongly improves the simulation of the diurnal cycles of precipitation and temperature, despite an enhanced warm bias and a tendency for the overestimation of precipitation over the Alps. The CPM12 model has a poor diurnal cycle associated with the use of parameterized convection. The assessment of extreme precipitation events reveals that both models adequately represent the frequency-intensity distributions for daylong events in summer, but large differences occur for hourly precipitation. The CPM12 model underestimates the frequency of heavy hourly events, while CRM2 shows good agreement with observations in the summer season. We also present results on the scaling of precipitation extremes with local daily mean temperatures. In accordance with observations, CRM2 exhibits adiabatic scaling for intermediate hourly events (90th percentile) and superadiabatic scaling for extreme hourly events (99th and 99.9th percentiles) during the summer season. The CPM12 model partly reproduces this scaling as well. The excellent performance of CRM2 in representing hourly precipitation events in terms of intensity and scaling is highly encouraging, as this addresses a previously untested (and untuned) model capability.

Representing Equilibrium and Nonequilibrium Convection in Large-Scale Models

Abstract: A new diagnostic convective closure, which is dependent on convective available potential energy (CAPE), is derived under the quasi-equilibrium assumption for the free troposphere subject to boundary layer forcing. The closure involves a convective adjustment time scale for the free troposphere and a coupling coefficient between the free troposphere and the boundary layer based on different time scales over land and ocean. Earlier studies with the ECMWF Integrated Forecasting System (IFS) have already demonstrated the model's ability to realistically represent tropical convectively coupled waves and synoptic variability with use of the "standard" CAPE closure, given realistic entrainment rates. A comparison of low-resolution seasonal integrations and high-resolution short-range forecasts against complementary satellite and radar data shows that with the extended CAPE closure it is also possible, independent of model resolution and time step, to realistically represent nonequilibrium convection such as the diurnal cycle of convection and the convection tied to advective boundary layers, although representing the late night convection over land remains a challenge. A more in-depth regional analysis of the diurnal cycle and the closure is provided for the continental United States and particularly Africa, including comparison with data from satellites and a cloud-resolving model (CRM). Consequences for global numerical weather prediction (NWP) are not only a better phase representation of convection, but also better forecasts of its spatial distribution and local intensity.

Cessation of deep convection in the open Southern Ocean under anthropogenic climate change

Abstract: In 1974, newly available satellite observations unveiled the presence of a giant ice-free area, or polynya, within the Antarctic ice pack of the Weddell Sea, which persisted during the two following winters1. Subsequent research showed that deep convective overturning had opened a conduit between the surface and the abyssal ocean, and had maintained the polynya through the massive release of heat from the deep sea2, 3. Although the polynya has aroused continued interest1, 2, 3, 4, 5, 6, 7, 8, 9, the presence of a fresh surface layer has prevented the recurrence of deep convection there since 19768, and it is now largely viewed as a naturally rare event10. Here, we present a new analysis of historical observations and model simulations that suggest deep convection in the Weddell Sea was more active in the past, and has been weakened by anthropogenic forcing. The observations show that surface freshening of the southern polar ocean since the 1950s has considerably enhanced the salinity stratification. Meanwhile, among the present generation of global climate models, deep convection is common in the Southern Ocean under pre-industrial conditions, but weakens and ceases under a climate change scenario owing to surface freshening. A decline of open-ocean convection would reduce the production rate of Antarctic Bottom Waters, with important implications for ocean heat and carbon storage, and may have played a role in recent Antarctic climate change.

Moist Static Energy Budget of the MJO during DYNAMO

Abstract: Theauthors analyzethe column-integratedmoist static energybudgetover theregionof the tropical Indian Ocean covered by the sounding array during the Cooperative Indian Ocean Experiment on Intraseasonal Variability in the Year 2011 (CINDY2011)/Dynamics of the Madden‐Julian Oscillation (DYNAMO) field experiment in late 2011. The analysis is performed using data from the sounding array complemented by additional observational datasets for surface turbulent fluxes and atmospheric radiative heating. The entire analysis is repeated using the ECMWF Interim Re-Analysis (ERA-Interim). The roles of surface turbulent fluxes, radiative heating, and advection are quantified for the two MJO events that occurred in October and November using the sounding data; a third event in December is also studied in the ERA-Interim data. These results are consistent with the view that the MJO’s moist static energy anomalies grow and are sustained to a significant extent by the radiative feedbacks associated with MJO water vapor and cloud anomalies and that propagation of the MJO is associated with advection of moist static energy. Both horizontal and vertical advection appear to play significant roles in the events studied here. Horizontal advection strongly moistens the atmosphere during the buildup to the active phase of the October event when the lowlevel winds switch from westerly to easterly. Horizontal advection strongly dries the atmosphere in the wake of the active phases of the November and December events as the westerlies associated with off-equatorial cyclonic gyres bring subtropical dry air into the convective region from the west and north. Vertical advection provides relative moistening ahead of the active phase and drying behind it, associated with an increase of the normalized gross moist stability.

Propagating versus Nonpropagating Madden–Julian Oscillation Events

Abstract: Basinwide convective anomalies over the Indian Ocean (IO) associated with the Madden–Julian oscillation (MJO) sometimes propagate eastward and reach the west Pacific (WP), but sometimes do not. Long-term observations and reanalysis products are used to investigate the difference between the propagating and nonpropagating MJO events. IO convection onset events associated with the MJO are grouped into three categories based on the strengths of the simultaneous dry anomalies over the eastern Maritime Continent and WP. The IO convection anomaly preferentially makes eastward propagation and reaches the WP when the dry anomaly is stronger.Analysis of the column-integrated moist static energy (MSE) budget shows that horizontal advection moistens the atmosphere to the east of the positive MSE anomaly associated with the active convection over the IO and is of sufficient magnitude to explain the eastward propagation of the positive MSE anomaly. Interpretation is complicated, however, by lack of closure in ...

Improvements to stellar structure models, based on a grid of 3D convection simulations. II. Calibrating the mixing-length formulation

Abstract: Weperformacalibrationofthemixinglengthofconvectioninstellarstructuremodelsagainstrealistic3Dradiation-coupledhydrodynamics(RHD)simulationsofconvectioninstellarsurfacelayers,determiningtheadiabatdeepinconvectivestellarenvelopes.The mixing-length parameter α is calibrated by matching averages of the3Dsimulationsto1Dstellarenvelopemodels,ensuringidenticalatomicphysicsinthetwocases.Thisisdoneforapreviouslypublishedgridofsolar-metallicityconvectionsimulations,coveringfrom4200Kto6900Konthemainsequence,and4300–5000Kforgiantswithlogg =2.2.Ourcalibrationresultsinanα varyingfrom1.6forthewarmestdwarf,whichisjustcoolenoughtoadmitaconvectiveenvelope,andupto2.05forthecoolestdwarfsinourgrid.Inbetweentheseisatriangularplateauofα ∼1.76.TheSunislocatedonthisplateauandhasseenlittlechangeduringitsevolutionsofar.Whenstarsascendthegiantbranch,theylargelydosoalongtracksofconstantα,withα decreasingwithincreasingmass.Key words: Stars:atmospheres–stars:evolution–convection 1 INTRODUCTIONDue to the lack of a better theory of convection in stars,the mixing-length theory (B¨ohm-Vitense 1958, MLT) hasbeen in use for more thanhalf a century.By far the largestpart of the solar convection zone is very close to adia-batic, and the stratification in the bulk of the convectionzone is therefore determined by the adiabatic gradient,∇

A simulation of convective dynamo in the solar convective envelope: maintenance of the solar-like differential rotation and emerging flux

Abstract: We report the results of a magnetohydrodynamic (MHD) simulation of a convective dynamo in a model solar convective envelope driven by the solar radiative diffusive heat flux. The convective dynamo produces a large-scale mean magnetic field that exhibits irregular cyclic behavior with oscillation time scales ranging from about 5 to 15 yr and undergoes irregular polarity reversals. The mean axisymmetric toroidal magnetic field is of opposite signs in the two hemispheres and is concentrated at the bottom of the convection zone. The presence of the magnetic fields is found to play an important role in the self-consistent maintenance of a solar-like differential rotation in the convective dynamo model. Without the magnetic fields, the convective flows drive a differential rotation with a faster rotating polar region. In the midst of magneto-convection, we found the emergence of strong super-equipartition flux bundles at the surface, exhibiting properties that are similar to emerging solar active regions.

Approaching the asymptotic regime of rapidly rotating convection: boundary layers versus interior dynamics.

Abstract: Rapidly rotating Rayleigh-Benard convection is studied by combining results from direct numerical simulations (DNS), laboratory experiments, and asymptotic modeling. The asymptotic theory is shown to provide a good description of the bulk dynamics at low, but finite Rossby number. However, large deviations from the asymptotically predicted heat transfer scaling are found, with laboratory experiments and DNS consistently yielding much larger Nusselt numbers than expected. These deviations are traced down to dynamically active Ekman boundary layers, which are shown to play an integral part in controlling heat transfer even for Ekman numbers as small as 10 −7 . By adding an analytical parametrization of the Ekman transport to simulations using stress-free boundary conditions, we demonstrate that the heat transfer jumps from values broadly compatible with the asymptotic theory to states of strongly increased heat transfer, in good quantitative agreement with no-slip DNS and compatible with the experimental data. Finally, similarly to nonrotating convection, we find no single scaling behavior, but instead that multiple well-defined dynamical regimes exist in rapidly rotating convection systems.

Numerical Simulations of Active Region Scale Flux Emergence: From Spot Formation to Decay

Abstract: We present numerical simulations of active region scale flux emergence covering a time span of up to 6 days. Flux emergence is driven by a bottom boundary condition that advects a semi-torus of magnetic field with 1.7 × 1022 Mx flux into the computational domain. The simulations show that, even in the absence of twist, the magnetic flux is able the rise through the upper 15.5 Mm of the convection zone and emerge into the photosphere to form spots. We find that spot formation is sensitive to the persistence of upflows at the bottom boundary footpoints, i.e., a continuing upflow would prevent spot formation. In addition, the presence of a torus-aligned flow (such flow into the retrograde direction is expected from angular momentum conservation during the rise of flux ropes through the convection zone) leads to a significant asymmetry between the pair of spots, with the spot corresponding to the leading spot on the Sun being more axisymmetric and coherent, but also forming with a delay relative to the following spot. The spot formation phase transitions directly into a decay phase. Subsurface flows fragment the magnetic field and lead to intrusions of almost field free plasma underneath the photosphere. When such intrusions reach photospheric layers, the spot fragments. The timescale for spot decay is comparable to the longest convective timescales present in the simulation domain. We find that the dispersal of flux from a simulated spot in the first two days of the decay phase is consistent with self-similar decay by turbulent diffusion.

Soil moisture-temperature feedbacks at meso-scale during summer heat waves over Western Europe

Abstract: This paper investigates the impact of soil moisture-temperature feedback during heatwaves occurring over France between 1989 and 2008. Two simulations of the weather research and forecasting regional model have been analysed, with two different land-surface models. One resolves the hydrology and is able to simulate summer dryness, while the other prescribes constant and high soil moisture and hence no soil moisture deficit. The sensitivity analysis conducted for all heatwave episodes highlights different soil moisture-temperature responses (1) over low-elevation plains, (2) over mountains and (3) over coastal regions. In the plains, soil moisture deficit induces less evapotranspiration and higher sensible heat flux. This has the effect of heating the planetary boundary layer and at the same time of creating a general condition of higher convective instability and a slight increase of shallow cloud cover. A positive feedback is created which increases the temperature anomaly during the heatwaves. In mountainous regions, enhanced heat fluxes over dry soil reinforce upslope winds producing strong vertical motion over the mountain slope, first triggered by thermal convection. This, jointly to the instability conditions, favors convection triggering and produces clouds and precipitation over the mountains, reducing the temperature anomaly. In coastal regions, dry soil enhances land/sea thermal contrast, strengthening sea-breeze circulation and moist cold marine air advection. This damps the magnitude of the heatwave temperature anomaly in coastal areas, expecially near the Mediterranean coast. Hence, along with heating in the plains, soil dryness can also have a significant cooling effect over mountains and coastal regions due to meso-scale circulations.

Radiative-convective instability

Abstract: [1] Radiative-moist-convective equilibrium (RCE) is a simple paradigm for the statistical equilibrium the earth's climate would exhibit in the absence of lateral energy transport. It has generally been assumed that for a given solar forcing and long-lived greenhouse gas concentration, such a state would be unique, but recent work suggests that more than one stable equilibrium may be possible. Here we show that above a critical specified sea surface temperature, the ordinary RCE state becomes linearly unstable to large-scale overturning circulations. The instability migrates the RCE state toward one of the two stable equilibria first found by Raymond and Zeng (2000). It occurs when the clear-sky infrared opacity of the lower troposphere becomes so large, owing to high water vapor concentration, that variations of the radiative cooling of the lower troposphere are governed principally by variations in upper tropospheric water vapor. We show that the instability represents a subcritical bifurcation of the ordinary RCE state, leading to either a dry state with large-scale descent, or to a moist state with mean ascent; these states may be accessed by finite amplitude perturbations to ordinary RCE in the subcritical state, or spontaneously in the supercritical state. As first suggested by Raymond (2000) and Sobel et al. (2007), the latter corresponds to the phenomenon of self-aggregation of moist convection, taking the form of cloud clusters or tropical cyclones. We argue that the nonrobustness of self-aggregation in cloud system resolving models may be an artifact of running such models close to the critical temperature for instability.

Solutal-convection regimes in a two-dimensional porous medium

Abstract: We numerically characterize the temporal regimes for solutal convection from almost first contact to high dissolved solute concentration in a two-dimensional ideal porous layer for Rayleigh numbers -dependent shut-down regime. For the constant-flux and shut-down regimes, we rigourously derive upscaled equations connecting the horizontally averaged concentration, vertical advective flux and plume widths. These are partially complete; a universal expression for the plume width remains elusive. We complement these governing equations with phenomenological boundary conditions based on a marginally stable diffusive boundary layer at the top and zero advective flux at the bottom. Making appropriate approximations in each regime, we find good agreement between predictions from this model and simulated results for both solutal and thermal convection. In the partially permeable upper boundary case, fluid from the convecting layer can penetrate an overlying separate-phase-solute bearing layer where it immediately saturates. The regime diagram remains almost the same as for the impermeable case, but the dissolution flux is significantly augmented. Our work is motivated by dissolution of carbon dioxide relevant to geological storage, and we conclude with a simple flux parameterization for inclusion in gravity current models and suggest that the upscaled equations could lay the foundation for accurate inclusion of dissolution in reservoir simulators.

Investigation of refrigeration efficiency for fully wet circular porous fins with variable sections by combined heat and mass transfer analysis

Abstract: Temperature distribution equation and refrigeration efficiency for fully wet circular porous fins with variable sections are introduced in this study by a new modified wet fin parameter presented by Sharqawy and Zubair. This parameter can be calculated without knowing the fin tip condition by considering the temperature and humidity ratio differences for the driving forces of heat and mass transfer, respectively. It's assumed that heat and mass convective coefficients vary with fin temperature and heat transfer through porous media is simulated using passage velocity from the Darcy's model. After presenting the governing equation, Least Square Method (LSM) and fourth order Runge-Kutta method (NUM) are applied for predicting the temperature distribution in the sample aluminum porous fins. After that, effects of porosity, Darcy number, Rayleigh number, Lewis number and etc. on fin efficiency are examined. As a main outcome, for reaching to high values of fin efficiency, rectangular fin should be used instead of convex and triangular sections.

Upscale energy transfer in three-dimensional rapidly rotating turbulent convection.

Abstract: Rotating Rayleigh-Benard convection exhibits, in the limit of rapid rotation, a turbulent state known as geostrophic turbulence. This state is present for sufficiently large Rayleigh numbers representing the thermal forcing of the system, and is characterized by a leading order balance between the Coriolis force and pressure gradient. This turbulent state is itself unstable to the generation of depth-independent or barotropic vortex structures of ever larger scale through a process known as spectral condensation. This process involves an inverse cascade mechanism with a positive feedback loop whereby large-scale barotropic vortices organize small scale convective eddies. In turn, these eddies provide a dynamically evolving energy source for the large-scale barotropic component. Kinetic energy spectra for the barotropic dynamics are consistent with a k-3 downscale enstrophy cascade and an upscale cascade that steepens to k-3 as the box-scale condensate forms. At the same time the flow maintains a baroclinic convective component with an inertial range consistent with a k-5/3 spectrum. The condensation process resembles a similar process in two dimensions but is fully three-dimensional.

A Unified Convection Scheme (UNICON). Part I: Formulation

Abstract: The author develops a unified convection scheme (UNICON) that parameterizes relative (i.e., with respect to the grid-mean vertical flow) subgrid vertical transport by nonlocal asymmetric turbulent eddies. UNICON is a process-based model of subgrid convective plumes and mesoscale organized flow without relying on any quasi-equilibrium assumptions such as convective available potential energy (CAPE) or convective inhibition (CIN) closures. In combination with a relative subgrid vertical transport scheme by local symmetric turbulent eddies and a grid-scale advection scheme, UNICON simulates vertical transport of water species and conservative scalars without double counting at any horizontal resolution.UNICON simulates all dry–moist, forced–free, and shallow–deep convection within a single framework in a seamless, consistent, and unified way. It diagnoses the vertical profiles of the macrophysics (fractional area, plume radius, and number density) as well as the microphysics (production and evaporati...

Large-scale vortices in rapidly rotating Rayleigh-Bénard convection

Abstract: Using numerical simulations of rapidly rotating Boussinesq convection in a Cartesian box, we study the formation of long-lived, large-scale, depth-invariant coherent structures. These structures, which consist of concentrated cyclones, grow to the horizontal scale of the box, with velocities significantly larger than the convective motions. We vary the rotation rate, the thermal driving and the aspect ratio in order to determine the domain of existence of these large-scale vortices (LSV). We find that two conditions are required for their formation. First, the Rayleigh number, a measure of the thermal driving, must be several times its value at the linear onset of convection; this corresponds to Reynolds numbers, based on the convective velocity and the box depth, . Second, the rotational constraint on the convective structures must be strong. This requires that the local Rossby number, based on the convective velocity and the horizontal convective scale, . Simulations in which certain wavenumbers are artificially suppressed in spectral space suggest that the LSV are produced by the interactions of small-scale, depth-dependent convective motions. The presence of LSV significantly reduces the efficiency of the convective heat transport.

Heat and mass transfer in unsteady rotating fluid flow with binary chemical reaction and activation energy.

Abstract: In this study, the Spectral Relaxation Method (SRM) is used to solve the coupled highly nonlinear system of partial differential equations due to an unsteady flow over a stretching surface in an incompressible rotating viscous fluid in presence of binary chemical reaction and Arrhenius activation energy. The velocity, temperature and concentration distributions as well as the skin-friction, heat and mass transfer coefficients have been obtained and discussed for various physical parametric values. The numerical results obtained by (SRM) are then presented graphically and discussed to highlight the physical implications of the simulations.

A multi-model comparison of Atlantic multidecadal variability

Abstract: A multi-model analysis of Atlantic multidecadal variability is performed with the following aims: to investigate the similarities to observations; to assess the strength and relative importance of the different elements of the mechanism proposed by Delworth et al. (J Clim 6:1993–2011, 1993) (hereafter D93) among coupled general circulation models (CGCMs); and to relate model differences to mean systematic error. The analysis is performed with long control simulations from ten CGCMs, with lengths ranging between 500 and 3600 years. In most models the variations of sea surface temperature (SST) averaged over North Atlantic show considerable power on multidecadal time scales, but with different periodicity. The SST variations are largest in the mid-latitude region, consistent with the short instrumental record. Despite large differences in model configurations, we find quite some consistency among the models in terms of processes. In eight of the ten models the mid-latitude SST variations are significantly correlated with fluctuations in the Atlantic meridional overturning circulation (AMOC), suggesting a link to northward heat transport changes. Consistent with this link, the three models with the weakest AMOC have the largest cold SST bias in the North Atlantic. There is no linear relationship on decadal timescales between AMOC and North Atlantic Oscillation in the models. Analysis of the key elements of the D93 mechanisms revealed the following: Most models present strong evidence that high-latitude winter mixing precede AMOC changes. However, the regions of wintertime convection differ among models. In most models salinity-induced density anomalies in the convective region tend to lead AMOC, while temperature-induced density anomalies lead AMOC only in one model. However, analysis shows that salinity may play an overly important role in most models, because of cold temperature biases in their relevant convective regions. In most models subpolar gyre variations tend to lead AMOC changes, and this relation is strong in more than half of the models.

Inverse cascade and symmetry breaking in rapidly rotating Boussinesq convection

Abstract: In this paper, we present numerical simulations of rapidly rotating Rayleigh-Benard convection in the Boussinesq approximation with stress-free boundary conditions. At moderately low Rossby number and large Rayleigh number, we show that a large-scale depth-invariant flow is formed, reminiscent of the condensate state observed in two-dimensional flows. We show that the large-scale circulation shares many similarities with the so-called vortex, or slow-mode, of forced rotating turbulence. Our investigations show that at a fixed rotation rate the large-scale vortex is only observed for a finite range of Rayleigh numbers, as the quasi-two-dimensional nature of the flow disappears at very high Rayleigh numbers. We observe slow vortex merging events and find a non-local inverse cascade of energy in addition to the regular direct cascade associated with fast small-scale turbulent motions. Finally, we show that cyclonic structures are dominant in the small-scale turbulent flow and this symmetry breaking persists in the large-scale vortex motion.

Large-scale vortices in rapidly rotating Rayleigh-B\'enard convection

Abstract: Using numerical simulations of rapidly rotating Boussinesq convection in a Cartesian box, we study the formation of long-lived, large-scale, depth-invariant coherent structures. These structures, which consist of concentrated cyclones, grow to the horizontal size of the box, with velocities significantly larger than the convective motions. We vary the rotation rate, the thermal driving and the aspect ratio in order to determine the domain of existence of these large-scale vortices (LSV). We find that two conditions are required for their formation. First, the Rayleigh number, a meaure of the thermal driving, must be several times its value at the linear onset of convection; this corresponds to Reynolds numbers, based on the convective velocity and the box depth, $\gtrsim 100$. Second, the rotational constraint on the convective structures must be strong. This requires that the local Rossby number, based on the convective velocity and the horizontal convective scale, $\lesssim 0.15$. Simulations in which certain wavenumbers are artificially suppressed in spectral space suggest that the LSV are produced by the interactions of small-scale, depth-dependent convective motions. The presence of LSV significantly reduces the efficiency of the convective heat transport.

Understanding and controlling plasmon-induced convection.

Abstract: The heat generation and fluid convection induced by plasmonic nanostructures is attractive for optofluidic applications. However, previously published theoretical studies predict only nanometre per second fluid velocities that are inadequate for microscale mass transport. Here we show both theoretically and experimentally that an array of plasmonic nanoantennas coupled to an optically absorptive indium-tin-oxide (ITO) substrate can generate >micrometre per second fluid convection. Crucially, the ITO distributes thermal energy created by the nanoantennas generating an order of magnitude increase in convection velocities compared with nanoantennas on a SiO2 base layer. In addition, the plasmonic array alters absorption in the ITO, causing a deviation from Beer-Lambert absorption that results in an optimum ITO thickness for a given system. This work elucidates the role of convection in plasmonic optical trapping and particle assembly, and opens up new avenues for controlling fluid and mass transport on the micro- and nanoscale.