Large-eddy simulation (LES) is used to model turbulent winds in a nominally neutral atmospheric boundary layer at varying mesh resolutions. The boundary layer is driven by wind shear with zero surface heat flux and is capped by a stable inversion. Because of entrainment the boundary layer is in a weakly stably stratified regime. The simulations use meshes varying from 1282× 64 to 10242× 512 grid points in a fixed computational domain of size (2560, 2560, 896) m. The subgrid-scale (SGS) parameterizations used in the LES vary with the mesh spacing. Low-order statistics, spectra, and structure functions are compared on the different meshes and are used to assess grid convergence in the simulations. As expected, grid convergence is primarily achieved in the middle of the boundary layer where there is scale separation between the energy-containing and dissipative eddies. Near the surface second-order statistics do not converge on the meshes studied. The analysis also highlights differences between one-dimensional and two-dimensional velocity spectra; differences are attributed to sampling errors associated with aligning the horizontal coordinates with the vertically veering mean wind direction. Higher-order structure functions reveal non-Gaussian statistics on all scales, but are highly dependent on the mesh resolution. A generalized logarithmic law and a k−1spectral scaling regime are identified with mesh-dependent parameters in agreement with previously published results.
- Award ID(s):
- 2103874
- NSF-PAR ID:
- 10331091
- Date Published:
- Journal Name:
- Physical review fluids
- ISSN:
- 2469-9918
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
more » « less
Abstract We present observations and modeling of spatial eigen‐functions of resonating waves within fault zone waveguide, using data recorded on a dense seismic array across the San Jacinto Fault Zone (SJFZ) in southern California. The array consists of 5‐Hz geophones that cross the SJFZ with ~10–30 m spacing at the Blackburn Saddle near the Hemet Stepover. Wavefield snapshots after the
S wave arrival are consistent for more than 50 near‐fault events, suggesting that this pattern is controlled by the fault zone structure rather than source properties. Data from example event with high signal to noise ratio show three main frequency peaks at ~1.3, ~2.0, and ~2.8 Hz in the amplitude spectra of resonance waves averaged over stations near the fault. The data are modeled with analytical expressions for eigen‐functions of resonance waves in a low‐velocity layer (fault zone) between two quarter‐spaces. Using a grid search‐based method, we investigate the possible width of the waveguide, location within the array, and shear wave velocities of the media that fit well the resonance signal at ~1.3 Hz. The results indicate a ~300 m wide damaged fault zone layer with ~65%S wave velocity reduction compared to the host rock. The SW edge of the low‐velocity zone is near the mapped fault surface trace, indicating that the damage zone is asymmetrically located at the regionally faster NE crustal block. The imaging resolution of the fault zone structure can be improved by modeling fault zone resonance modes and trapped waves together. -
more » « less
Context. Molecular outflows are believed to be a key ingredient in the process of star formation. The molecular outflow associated with DR21 Main in Cygnus-X is one of the most extreme molecular outflows in the Milky Way in terms of mass and size. The outflow is suggested to belong to a rare class of explosive outflows formed by the disintegration of protostellar systems.Aims. We aim to explore the morphology, kinematics, and energetics of the DR21 Main outflow, and to compare those properties to confirmed explosive outflows in order to unravel the underlying driving mechanism behind DR21.Methods. We studied line and continuum emission at a wavelength of 3.6 mm with IRAM 30 m and NOEMA telescopes as part of the Cygnus Allscale Survey of Chemistry and Dynamical Environments (CASCADE) program. The spectra include (J = 1−0) transitions of HCO+, HCN, HNC, N2H+, H2CO, and CCH, which trace different temperature and density regimes of the outflowing gas at high velocity resolution (~0.8 km s−1). The map encompasses the entire DR21 Main outflow and covers all spatial scales down to a resolution of 3″ (~0.02 pc).Results. Integrated intensity maps of the HCO+emission reveal a strongly collimated bipolar outflow with significant overlap of the blueshifted and redshifted emission. The opening angles of both outflow lobes decrease with velocity, from ~80 to 20° for the velocity range from 5 to 45 km s−1relative to the source velocity. No evidence is found for the presence of elongated, “filament-like” structures expected in explosive outflows. N2H+emission near the western outflow lobe reveals the presence of a dense molecular structure, which appears to be interacting with the DR21 Main outflow.Conclusions. The overall morphology as well as the detailed kinematics of the DR21 Main outflow are more consistent with a typical bipolar outflow than with an explosive counterpart. -
A realistic approach for gathering high-resolution observations of the rainfall rate, R, in the vertical plane is to use data from vertically pointing Doppler radars. After accounting for the vertical air velocity and attenuation, it is possible to determine the fine, spatially resolved drop size spectra and to calculate R for further statistical analyses. The first such results in a vertical plane are reported here. Specifically, we present results using MRR-Pro Doppler radar observations at resolutions of ten meters in height over the lowest 1.28 km, as well as ten seconds in time, over four sets of observations using two different radars at different locations. Both the correlation functions and power spectra are useful for translating observations and numerical model outputs of R from one scale down to other scales that may be more appropriate for particular applications, such as flood warnings and soil erosion, for example. However, it was found in all cases that, while locally applicable radial power spectra could be calculated, because of statistical heterogeneity most of the power spectra lost all generality, and proper correlation functions could not be computed in general except for one 17-min interval. Nevertheless, these results are still useful since they can be combined to develop catalogs of power spectra over different meteorological conditions and in different climatological settings and locations. Furthermore, even with the limitations of these data, this approach is being used to gain a deeper understanding of rainfall to be reported in a forthcoming paper.more » « less
-
more » « less
ABSTRACT In this study, we explore the statistics of pressure fluctuations in kinetic collisionless turbulence. A 2.5D kinetic particle-in-cell simulation of decaying turbulence is used to investigate pressure balance via the evolution of thermal and magnetic pressure in a plasma with β of order unity. We also discuss the behaviour of thermal, magnetic, and total pressure structure functions and their corresponding wavenumber spectra. The total pressure spectrum exhibits a slope of −7/3 extending for about a decade in the ion-inertial range. In contrast, shallower −5/3 spectra are characteristic of the magnetic pressure and thermal pressure. The steeper total pressure spectrum is a consequence of cancellation caused by density-magnetic field magnitude anti-correlation. Further, we evaluate higher order total pressure structure functions in an effort to discuss intermittency and compare the power exponents with higher order structure functions of velocity and magnetic fluctuations. Finally, applications to astrophysical systems are also discussed.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1103/PhysRevFluids.7.064601
