The Richardson's Law in Large-Eddy Simulations of Boundary Layer flows

Relative dispersion in a neutrally stratified planetary boundary layer (PBL) is investigated by means of Large-Eddy Simulations (LES). Despite the small extension of the inertial range of scales in the simulated PBL, our Lagrangian statistics turns out to be compatible with the Richardson $t^3$ law for the average of square particle separation. This emerges from the application of nonstandard methods of analysis through which a precise measure of the Richardson constant was also possible. Its values is estimated as $C_2\sim 0.5$ in close agreement with recent experiments and three-dimensional direct numerical simulations.Comment: 15 LaTex pages, 4 PS figure

Comparison of two- and three-dimensional simulations of miscible Rayleigh-Taylor instability

A comparison of two-dimensional and three-dimensional high-resolution numerical large-eddy simulations of planar, miscible Rayleigh-Taylor instability flows are presented. The resolution of the three-dimensional simulation is sufficient to attain a fully turbulent state. A number of different statistics from the mixing region (e.g., growth rates, PDFs, mixedness measures, and spectra) are used to demonstrate that two-dimensional flow simulations differ substantially from the three-dimensional one. It is found that the two-dimensional flow grows more quickly than its three-dimensional counterpart at late times, develops larger structures, and is much less well mixed. These findings are consistent with the concept of inverse cascade in two-dimensional flow, as well as the influence of a reduced effective Atwood number on miscible flow

Interactions between downslope flows and a developing cold-air pool

A numerical model has been used to characterize the development of a region of enhanced cooling in an alpine valley with a width of order (Formula presented.) km, under decoupled stable conditions. The region of enhanced cooling develops largely as a region of relatively dry air which partitions the valley atmosphere dynamics into two volumes, with airflow partially trapped within the valley by a developing elevated inversion. Complex interactions between the region of enhanced cooling and the downslope flows are quantified. The cooling within the region of enhanced cooling and the elevated inversion is almost equally partitioned between radiative and dynamic effects. By the end of the simulation, the different valley atmospheric regions approach a state of thermal equilibrium with one another, though this cannot be said of the valley atmosphere and its external environment.Peer reviewe

Evaluation of large-Eddy simulations via observations of nocturnal marine stratocumulus

Data from the first research flight (RF01) of the second Dynamics and Chemistry of Marine Stratocumulus (DYCOMS-II) field study are used to evaluate the fidelity with which large-eddy simulations (LESs) can represent the turbulent structure of stratocumulus-topped boundary layers. The initial data and forcings for this case placed it in an interesting part of parameter space, near the boundary where cloud-top mixing is thought to render the cloud layer unstable on the one hand, or tending toward a decoupled structure on the other hand. The basis of this evaluation consists of sixteen 4-h simulations from 10 modeling centers over grids whose vertical spacing was 5 m at the cloud-top interface and whose horizontal spacing was 35 m. Extensive sensitivity studies of both the configuration of the case and the numerical setup also enhanced the analysis. Overall it was found that (i) if efforts are made to reduce spurious mixing at cloud top, either by refining the vertical grid or limiting the effects of the subgrid model in this region, then the observed turbulent and thermodynamic structure of the layer can be reproduced with some fidelity; (ii) the base, or native configuration of most simulations greatly overestimated mixing at cloud top, tending toward a decoupled layer in which cloud liquid water path and turbulent intensities were grossly underestimated; (iii) the sensitivity of the simulations to the representation of mixing at cloud top is, to a certain extent, amplified by particulars of this case. Overall the results suggest that the use of LESs to map out the behavior of the stratocumulus-topped boundary layer in this interesting region of parameter space requires a more compelling representation of processes at cloud top. In the absence of significant leaps in the understanding of subgrid-scale (SGS) physics, such a representation can only be achieved by a significant refinement in resolution—a refinement that, while conceivable given existing resources, is probably still beyond the reach of most centers

Revival of the magnetar PSR J1622-4950: observations with MeerKAT, Parkes, XMM-Newton, Swift, Chandra, and NuSTAR

New radio (MeerKAT and Parkes) and X-ray (XMM-Newton, Swift, Chandra, and NuSTAR) observations of PSR J1622-4950 indicate that the magnetar, in a quiescent state since at least early 2015, reactivated between 2017 March 19 and April 5. The radio flux density, while variable, is approximately 100x larger than during its dormant state. The X-ray flux one month after reactivation was at least 800x larger than during quiescence, and has been decaying exponentially on a 111+/-19 day timescale. This high-flux state, together with a radio-derived rotational ephemeris, enabled for the first time the detection of X-ray pulsations for this magnetar. At 5%, the 0.3-6 keV pulsed fraction is comparable to the smallest observed for magnetars. The overall pulsar geometry inferred from polarized radio emission appears to be broadly consistent with that determined 6-8 years earlier. However, rotating vector model fits suggest that we are now seeing radio emission from a different location in the magnetosphere than previously. This indicates a novel way in which radio emission from magnetars can differ from that of ordinary pulsars. The torque on the neutron star is varying rapidly and unsteadily, as is common for magnetars following outburst, having changed by a factor of 7 within six months of reactivation.Comment: Published in ApJ (2018 April 5); 13 pages, 4 figure

Large-Eddy Simulation of Flow and Pollutant Transport in Urban Street Canyons with Ground Heating

Our study employed large-eddy simulation (LES) based on a one-equation subgrid-scale model to investigate the flow field and pollutant dispersion characteristics inside urban street canyons. Unstable thermal stratification was produced by heating the ground of the street canyon. Using the Boussinesq approximation, thermal buoyancy forces were taken into account in both the Navier–Stokes equations and the transport equation for subgrid-scale turbulent kinetic energy (TKE). The LESs were validated against experimental data obtained in wind-tunnel studies before the model was applied to study the detailed turbulence, temperature, and pollutant dispersion characteristics in the street canyon of aspect ratio 1. The effects of different Richardson numbers (Ri) were investigated. The ground heating significantly enhanced mean flow, turbulence, and pollutant flux inside the street canyon, but weakened the shear at the roof level. The mean flow was observed to be no longer isolated from the free stream and fresh air could be entrained into the street canyon at the roof-level leeward corner. Weighed against higher temperature, the ground heating facilitated pollutant removal from the street canyon.Singapore-MIT Alliance for Research and Technology. Center for Environmental Sensing and Monitorin