Dynamics and spectra of cumulant update closures for two-dimensional turbulence

Abstract

Non-Markovian closure theories are compared with ensemble averaged direct numerical simulations (DNS) for decaying two-dimensional turbulence at large scale Reynolds numbers ranging from ≈ 50 to ≈ 4000. The closures, as well as DNS, are formulated for discrete wave numbers relevant to flows on the doubly periodic domain and are compared with the results of continuous wave number closures. The direct interaction approximation (DIA), self-consistent field theory (SCFT) and local energy-transfer theory (LET) closures are also compared with cumulant update versions of these closures (CUDIA, CUSCFT, CULET). The cumulant update closures are shown to have comparable performance to the standard closures but are much more efficient allowing long time integrations.

The discrete wave number closures perform considerably better than continuous wave number closures as far as evolved energy and transfer spectra and skewness are concerned. The discrete wave number closures are in reasonable agreement with DNS in the energy containing range of the large scales for Reynolds numbers ranging from ≈ 50 to ≈ 4000. The closures tend to underestimate the enstrophy flux to high wave numbers, increasingly so with increasing Reynolds number, resulting in underestimation of small-scale kinetic energy.     

Modeling high-order statistics in the turbulent Ekman layer

Results from a direct numerical simulation (DNS) of the neutral and unstable turbulent Ekman layer at a Reynolds number of 1000 were used to evaluate turbulence closure models. For the neutrally stratified Ekman layer, the higher-order moments of velocity were examined and the accuracy of a kurtosis model was assessed. For the unstable Ekman layer, the analysis of higher-order moments was extended to temperature-velocity correlations. Model coefficients were optimised using DNS data and it was shown that the optimised models accurately captured the distributions of all fourth-order moments. These low-Reynolds number results can be extrapolated to higher Reynolds numbers to parameterise turbulence in other flow fields with rotational effects such as the atmospheric boundary layer.     

Numerical modeling of environmental flows using DAM: Some preliminary results

Development of closures and parameterizations for subgrid scale effects is a significant and longstanding problem in the numerical simulation of environmental flows. The model described herein uses a rigorous approach for developing double-averaged governing equations — first a traditional Reynolds averaging to derive the Reynolds averaged Navier-Stokes equation (RANS), then a volume average to derive a set of double-averaged equations (DANS). An existing finite element flow model is then modified to accommodate these equations. This process gives rise to several new terms that require closures, as well as a new equation for free surface elevation. This paper is directed toward model development and uses several existing closure schemes as test cases.     

广义散射层析成像反演

本文详细地给出了基于非均匀介质的体散射广义散射层析成像反演的基本理论.广义散射层析成像反演可以描述为波场的反向传播和对成像场进行局部波数域滤波的过程.在数值算例中,利用背景速度沿深度方向均匀变化的v（z）介质中的简单的方块作为速度异常体的模型,通过对该模型产生的低频的Born数据和声波的正演数据的测试,在对采集系统进行有限频率带宽和空间孔径的校正来进行局部成像矩阵谱的恢复中,可以看出模型中各点的谱在恢复后的质量无论从覆盖的面积范围还是幅值的均一性上都有着明显的提高;在对速度模型的重建中,广义散射层析成像反演能够很好地恢复速度模型的低频分量,即便是方块速度异常体相对于背景速度的平均速度扰动是23%也能很好地重建模型中的速度,且对于不同的背景速度模型基本上都能很好地恢复Marmousi速度模型的低频分量.所以该方法将基于Born模型的层析成像反演适应范围进行了一定程度的扩展.     

Storm surge and tidal interaction in the Tjeldsund channel, northern Norway

The aim of this study is to investigate tide–surge interaction in narrow channels with complex and relatively shallow topography. A high-resolution depth-integrated tidal and storm surge model has been implemented for the Tjeldsund channel which is an important sailing lane in northern Norway. A horizontal grid resolution down to 50 m is applied in order to represent the complex bottom topography and the formation of jets and small-scale eddies. Two typically storm surge events in December 2004 have been examined in detail. The tide–surge interaction is found to influence the generation of higher harmonics and the formation of eddies in the current field. In some cases, the magnitude of storm surge currents may reach the same magnitude as the tidal currents enhancing the formation of jets and eddies.     

Turbulence flux measurement above the overstory of a subtropical Pinus plantation over the hilly region in southeastern China

Continuous turbulence flux measurement using the eddy covariance (EC) technique was made from January 1 to December 31 in 2003 at two and three canopy heights of a subtropical Pinus plantation on the red earth hilly region in southeastern China. To be able to make sure that the measured turbulence flux will equal the net ecosystem/atmosphere exchange, the quality of the data has to be assessed. Three criteria were investigated here, including the power spectra and cospectra analyses, flux variance similarity (integral turbulence test) and energy balance closure. The spectral analyses suggested that above-canopy power spectral slopes for all velocity components and scalars such as CO2, H2O and air temperature followed the expected -2/3 power law in the inertial subrange, and their cospectral slopes were close to -4/3 power law in the inertial subrange. The important contribution of large-scale motions to energy and mass transfer above the canopy at higher measurement level was also confirmed by the spectral analyses. The eddy covariance systems have the ability to resolve fluctuations associated with small-scale eddies and did not induce an obvious underestimation of the measured turbulence flux. The Monin-Obukhov similarity functions for the normalized standard deviation of vertical wind speed and air temperature were well-defined functions of atmospheric stability at two heights above the forest canopy, which indicated that turbulence flux measurements made at two heights were within the surface layer. Nocturnal flux underestimation and departures of this normalized standard deviation of vertical wind speed similarity function from that expected from Monin-Obukhov theory were a function of friction velocity. Thus, an optimal criterion of friction velocity was determined to be greater than 0.2-0.3 m s-1 for nocturnal fluxes so that the eddy covariance flux measurement was under high turbulent mixing conditions. Energy balance closure reached about 72%-81% at the studied site, which was comparable to the 10%-30% of energy imbalance reported in the literature. However, the energy balance closure could be only used as a useful reference criterion.     

Global stability of the ballooning mode in a cylindrical model

Ballooning disturbances in a finite-pressure plasma in a curvilinear magnetic field are described by the system of coupled equations for the Alfvén and slow magnetosonic modes. In contrast to most previous works that locally analyzed the stability of small-scale disturbances using the dispersion relationship, a global analysis outside a WKB approximation but within a simple cylindrical geometry, when magnetic field lines are circles with constant curvature, is performed in the present work. This model is relatively simple; nevertheless, it has the singularities necessary for the formation of the ballooning mode: field curvature and non-uniform thermal plasma pressure. If the disturbance finite radial extent is taken into account, the instability threshold increases as compared to a WKB approximation. The simplified model used in this work made it possible to consider the pattern of unstable disturbances at arbitrary values of the azimuthal wavenumber (k _y ). Azimuthally large-scale disturbances can also be unstable, although the increment increases with decreasing azimuthal scale and reaches saturation when the scales are of the order of the pressure nonuniformity dimension.     

Atmospheric boundary layer over steep surface waves

Turbulent air-sea interactions coupled with the surface wave dynamics remain a challenging problem. The needs to include this kind of interaction into the coupled environmental, weather and climate models motivate the development of a simplified approximation of the complex and strongly nonlinear interaction processes. This study proposes a quasi-linear model of wind-wave coupling. It formulates the approach and derives the model equations. The model is verified through a set of laboratory (direct measurements of an airflow by the particle image velocimetry (PIV) technique) and numerical (a direct numerical simulation (DNS) technique) experiments. The experiments support the central model assumption that the flow velocity field averaged over an ensemble of turbulent fluctuations is smooth and does not demonstrate flow separation from the crests of the waves. The proposed quasi-linear model correctly recovers the measured characteristics of the turbulent boundary layer over the waved water surface.     

Seasonal closures as a measure of trawling effort control in two Mediterranean trawling grounds: effects on epibenthic communities

Within the framework of ecosystem-based management, we focused on the use of seasonal closures as effective measures to minimise the degradation of benthic communities by trawling. These closures imply the complete cessation of trawling fleet activity and are commonly used in the Mediterranean to reduce the annual fishing effort, with the ultimate goal of effective resource management. In this study, we aimed to investigate how epibenthic communities respond to seasonal closures. The potential benefits of short-term annual closures in two Mediterranean fishing grounds were evaluated by analysing changes in community structure and composition that were linked to the closure. A decrease of faunal abundance was observed with the resumption of fishing activity after the closure at both fishing grounds. Remarkably, results indicated that some large and mobile fauna were able to respond to these closures. We concluded that the currently planned closures are too short to benefit benthic communities.     

On the merger of subsurface isolated vortices

Vortex merger is a phenomenon characterizing the whole class of geophysical vortices, from atmospheric storms and large oceanic eddies up to small scale turbulence. Here we focus on the merger of subsurface oceanic anticyclones in an idealized primitive equations model. This study has been motivated by past and recent observations of colliding lens-like anticyclones off of Gibraltar Strait. The critical conditions for merger (critical merger distance and time needed for merger) are determined. We will show that the predictions of classical two-dimensional merger are not verified for subsurface isolated vortices. For instance, critical merger distances will be reduced because of the vortex potential vorticity (PV) structure. The post-merger characteristics of the vortex (radius, extension and PV), are also determined. Merger-related effects, like production of peripheral filaments and small-scale eddies are also investigated and suggest the contribution of merger in both direct and inverse energy cascades.     

Stochastic Models of Quasigeostrophic Turbulence

Atmospheric and oceanic eddies are believed to be manifestations of quasigeostrophic turbulence — turbulence that occurs in rapidly rotating, vertically stratified fluid systems. The heat, momentum, and water transport by these eddies constitute a significant component of the climate balance, without which climate change cannot be understood. A major, unsolved problem is whether the turbulent eddy fluxes can be parameterized in terms of the large-scale, background flow. In the past, stochastic models have been used quite extensively to investigate quasigeostrophic turbulence in the case in which the eddy statistics are isotropic and homogeneous. Unfortunately, these models ignore the background shear which is absolutely essential to maintaining the eddies in the presence of dissipation. Recent attempts to extend stochastic models to shear flows have shown significant skill in predicting the structure of the eddy fluxes in arbitrary, three-dimensionally varying flows. This paper provides an accessible introduction to these models. The topics reviewed include quasigeostrophic turbulence and two-dimensional turbulence, non-modal andoptimal perturbations, mathematical theory of stochastic models, stochastic model simulations with realistic background states, and recent closure theories. A list of unsolved problems concludes this review.     

基于入射角的两项流体阻抗反演方法

本文在孔隙弹性介质理论的指导下,基于入射角AVO近似方程推导了包含Russell流体项的两项AVO近似方程和相应的弹性阻抗方程,通过分析可知其精度符合反演要求.在贝叶斯理论框架下,建立了包含正则化约束的弹性阻抗反演方法,在此基础上直接提取Russell流体项.该方法可在缺少大角度叠前地震资料的情况下进行叠前直接反演得到流体因子,减少传统方法带来的累积误差.模型试算表明,该方法具有较好的准确度和稳定性.实际工区应用取得了良好的效果,表明该方法有实用性.     

On cumulus-scale transport of horizontal momentum in monsoon depression over India

A diagnostie method of cumulus parameterization is suggested in which vertical transport of horizontal momentum by cumulus-scale is derived by making use of large-scale vorticity as well as divergence budget equations. Data for composite monsoon depression over India available from our earlier studies used to test the method. As a first approximation, the results are obtained using only the vorticity budget equation.The results show that in the southwest sector of the monsoon depression, which is characterized by maximum cloudiness and precipitation, there is an excess of cyclonic vorticity in the lower troposphere and anticyclonic vorticity in the upper troposphere associated with the large-scale motion. The distribution of eddy vertical transport of horizontal momentum is such that anticyclonic vorticity is generated in the lower troposphere and cyclonic vorticity aloft. Cumulus-scale eddies thus work against the large-scale system and tend to off-set the large-scale imbalance in vorticity.     

Waves and tides responsible for the intermittent closure of the entrance of a small, sheltered tidal wetland at San Francisco, CA

Crissy Field Marsh (CFM; http://www.nps.gov/prsf/planyourvisit/crissy-field-marsh-and-beach.htm) is a small, restored tidal wetland located in the entrance to San Francisco Bay just east of the Golden Gate. The marsh is small but otherwise fairly typical of many such restored wetlands worldwide. The marsh is hydraulically connected to the bay and the adjacent Pacific Ocean by a narrow sandy channel. The channel often migrates and sometimes closes completely, which effectively blocks the tidal connection to the ocean and disrupts the hydraulics and ecology of the marsh. Field measurements of waves and tides have been examined in order to evaluate the conditions responsible for the intermittent closure of the marsh entrance. The most important factor found to bring about the entrance channel closure is the occurrence of large ocean waves. However, there were also a few closure events during times with relatively small offshore waves. Examination of the deep-water directional wave spectra during these times indicates the presence of a small secondary peak corresponding to long period swell from the southern hemisphere, indicating that CFM and San Francisco Bay in general may be more susceptible to long period ocean swell emanating from the south or southwest than the more common ocean waves coming from the northwest. The tidal records during closure events show no strong relationship between closures and tides, other than that closures tend to occur during multi-day periods with successively increasing high tides. It can be inferred from these findings that the most important process to the intermittent closure of the entrance to CFM is littoral sediment transport driven by the influence of ocean swell waves breaking along the CFM shoreline at oblique angles. During periods of large, oblique waves the littoral transport of sand likely overwhelms the scour potential of the tidal flow in the entrance channel.     

基于随机介质模型的储层非均质性分析

本文利用随机介质模型对复杂岩性储层进行了非均质性描述.利用模型特征量即非均质纵横比、纵向谱指数、横向谱指数以及扰动标准差等来模拟不同的随机介质.在前人工作的基础上,由某油气田的两口井资料估计储层非均质性能谱,从能谱曲线上提取储层纵向大小尺度非均质谱指数.通过将二维随机介质模拟的合成井记录互相关系数与实际井记录互相关系数进行分析比较,分别得到大小尺度非均质情况下最佳拟合时的横向谱指数和非均质性纵横比.以上求得的各种特征量从不同角度定量分析了储层非均质性的纵横向变化,为储层横向预测提供了依据.     

正则化等效层重力向下延拓方法

近年来,利用时移微重力技术进行储层开发监测受到国内外学者广泛关注.时移微重力观测数据存在信噪比低,信号弱的问题,难以实现储层内物质运移的定量解释.为压制数据噪声,增强有效弱信号,本文研究了利用Tikhonov正则化方法反演等效层（源）,并由等效源实现重力场向下延拓的方法;在此基础上,本文推导了波数域正则化等效源向下延拓算子.针对向下延拓场幅值衰减问题,提出了正则化等效源迭代补偿算法.通过模拟数据实验研究了不同深度正则化等效源滤波算子及向下延拓算子的波数响应;与波数域Tikhonov正则化向下延拓方法相比,正则化等效源向下延拓方法的延拓精度更高、更稳定.最后,将基于迭代补偿的正则化等效源向下延拓技术应用于实测时移微重力数据证实了该方法能够有效增强局部异常,实现时移微重力数据大深度稳定向下延拓.     

Certain spectral properties of cyclonic turbulence in the Earth’s liquid core

Rapid diurnal rotation of planets results in the appearance of cyclonic thermal turbulence in liquid cores, which is the source of generation of the observed magnetic fields. The model that makes it possible to reproduce characteristic features of small-scale geostrophic flows in the physical and wave spaces is considered in the work. The flows of energy and hydrodynamic helicity as a function of the wavenumber have been estimated. Joint existence of direct and inverse cascades has been indicated. The consequences for the Earth’s core and geodynamo problems have been considered.     

Stochastic modelling of unresolved eddy fluxes

A subgrid-scale parameterization scheme motivated by statistical closure theory, but employing statistics obtained from high-resolution direct numerical simulations, is applied to large eddy simulations of two-level quasigeostrophic turbulence on the sphere. It is shown that these parameterizations are consistent with the phenomenology of quasigeostrophic turbulence. The parameterizations consist of 2 × 2 dissipation and stochastic forcing covariance matrices at each wavenumber, with the off-diagonal elements of the matrices representing vertical mixing. Two flow regimes, characterized by their deformation scales, are considered, namely atmospheric and oceanic. In the former, the deformation scale is fully resolved, and the truncation scale is within the enstrophy cascading interial range. In the latter, the deformation scale is not fully resolved, and the truncation scale is within the energy cascading inertial range. It is demonstrated through numerical experiments that both stochastic and deterministic variants of the scheme give comparable results for the energy spectra in the atmospheric regime. In the oceanic regime, the stochastic variant again gives excellent results, but the deterministic variant is found to be numerically unstable.     

Idealised flow past an island in a dynamically adaptive finite element model

The problem of flow separation around islands is investigated using a dynamically adaptive finite element model to allow for resolution of the shear layers that form in the advent of separation. The changes in secondary circulation and vertical motion that occur in both attached and separated flows are documented, as is the degree of closure of the wake eddies. In the numerical experiments presented, the strongest motion always takes place at the sides of the idealised island, where flow curvature and shear act together to induce ascent. In contrast, it is the slower motion within the wake eddies that allow streamlines to extend from the bottom to the surface. We find no evidence for closure of the wake eddies. Rather, all of our separated experiments show that streamlines that pass through the eddies originate outside of the shear layers and frictional boundary layers on the upstream side of the idealised island. The numerical experiments demonstrate the potential for dynamically adaptive, unstructured meshes to resolve the separated shear layers that occur downstream of the idealised island, as well as the narrow boundary layers that form on the island itself.     

Anisotropy of the horizontal velocity fluctuation field in the large wavenumber region

Mesoscale atmospheric wind, temperature, and density fluctuations are described in terms of their power spectra. In the early 1980s, VanZandt[1] showed that the observed spectrum of the atmospheric fluctuations can result, as in the ocean, from an interna…