Large-eddy simulation of plant canopy flows using plant-scale representation

Turbulent flow in a corn canopy is simulated using large-eddy simulation (LES) with a Lagrangian dynamic Smagorinsky model. A new numerical representation of plant canopies is presented that resolves approximately the local structure of plants and takes into account their spatial arrangement. As a validation, computational results are compared with experimental data from recent field particle image velocimetry (PIV) measurements and two previous experimental campaigns. Numerical simulation using the traditional modelling method to represent the canopy (field-scale approach) is also conducted as a comparison to the plant-scale approach. The combination of temporal PIV data, LES and spatial PIV data allows us to couple a wide range of relevant turbulence scales. There is good agreement between experimental data and numerical predictions using the plant-scale approach in terms of various turbulence statistics. Within the canopy, the plant-scale approach also allows the capture of more details than the field-scale approach, including instantaneous gusts that penetrate deep inside the canopy.     

Experimental observation on the characteristics of the near-surface turbulence over the Antarctic ice sheets during the polar day period

Based on the ultrasonic anemometer/thermometer data in the East Antarctic coastal area ice sheets ob-tained first by Chinese scientists, turbulent intensity, kinetic energy and sensible heat of turbulence, surface roughness height drag coefficient and normalized variation were calculated and analysed using the eddy-correlation method. The results show that the values of roughness height and drag coefficient are 4.3 X 10 m and 1.8x 10 -³, respectively. These turbulent parameters have apparent diurnal variations. Project supported by the National Natural Science Foundation of China (Grant No. 49675252).     

Driving Galactic Turbulence by Supernova Explosions

We investigate the general properties of supernova driven interstellar turbulence using local three-dimensional MHD simulations under Galactic conditions. Our model includes the effects of large-scale shear due to Galactic differential rotation, density stratification, compressibility, magnetic fields, heating via supernova explosions and parameterized radiative cooling of the interstellar medium. In addition to investigating isolated supernova explosions we allow for multiple supernovae distributed randomly in the Galactic disc and exponentially in the vertical direction. Single supernova explosions drive a strong shock, the lifetime of which is approximately 2 Myr in our model. This stage is found to be characterized by a kinetic energy spectrum in the diffuse gas with spectral index consistent with k = –2. Large-scale shear and Coriolis force act on the supernova remnant producing some vorticity inside it, but this process was found to be very weak. In the case of multiple supernova explosions, older remnants have an important role causing density fluctuations in the interstellar medium. In this clumpy medium, the propagation velocity of the shock fronts changes due to the changing density, and vorticity is generated. In the absence of these supernova interactions the kinetic energy spectrum shows a relatively wide shock spectrum with spectral index k = –2, but when the supernova interactions become dominant the classical k = –5/3 spectrum is observed.     

Galactic Spiral Arms and Dynamo Control Parameters

We discuss the effects of galactic spiral arms on the -coefficient, turbulent diffusivity and turbulent energy density of the interstellar turbulence. We argue that the -coefficient and the dynamo number are larger in the interarm regions, whereas the kinetic energy density of turbulence is larger in the arms; the turbulent magnetic diffusivity can be only weakly affected by the spiral pattern.     

Lattice dynamical system modelling of molecular clouds

Because a comprehensive microscopic treatment of interstellar molecular clouds is out of reach, an alternative approach is proposed in which most of the crucial ingredients of the problem are considered, but at some 'minimal' level of modelling. This leads to the elaboration of a lattice dynamical system , i.e. a time-dependent, spatially extended, deterministic system of macroscopic cells coupled through radiative transfer. Each cell is characterized by a small set of variables and supports a caricatural chemistry possessing the essential dynamical features of more realistic reaction schemes.   This approach naturally precludes quantitative results, but allows heretofore unavailable insights into some of the basic mechanisms at play. We focus on the response of the transfer process and the chemistry to a frozen 'turbulent' velocity field. It is shown that the system settles generically into a state where the effective coupling between cells is neither local nor global, and for which no single length-scale exists.   The spectral lines reconstructed from the spatiotemporal evolution of our model may, depending on the velocity field, exhibit profiles ranging from Gaussian to bimodal with strong realization effects. In the bimodal case, the model intrinsically displays an energy cascade transport mechanism to the cells that cool most efficiently: the feedback of chemistry on radiative transfer cannot be neglected.   Finally, extensions of this work are discussed and future developments are outlined.     

潮汐河口二维动床紊流模型

给出了正交曲线坐标系下平面二维非恒定动床紊流模型基本方程组,即水流运动方程及潮汐水流不平衡输沙方程。对计算中的有关问题提出了处理方法,如基本资料及代表潮过程选取、糙率的确定、悬沙粒径及沉速、挟沙能力恢复饱和系数及河床冲淤判别条件等。以甬江口为例,计算了镇海弯道至大游山口门约5km河段潮汐水流潮位及流速过程、丁坝作用下的回流域及疏浚与整治相结合的河床变形,计算结果与野外资料吻合良好。在此基础上,计算了宁波大桥桥墩作用下的水流结构变化及对镇海港航道的冲淤影响。     

The linkage between velocity patterns and sediment entrainment in a forced‐pool and riffle unit

A field‐based project was initiated in order to characterize velocities and sediment entrainment in a forced‐pool and riffle sequence. Three‐dimensional velocities and turbulence intensities were measured with an acoustic Doppler velocimeter at 222 different points at three similar flows that averaged approximately 4·35 m³ s⁻¹ within a large pool–riffle unit on North Saint Vrain Creek, Colorado. Sediment‐sorting patterns were observed with the introduction of 500 tracer particles painted according to initial seeding location. Tracer particles moved sporadically during a 113 day period in response to the annual snowmelt peak flow, which reached a maximum level of 14·8 m³ s⁻¹. Velocity data indicate high instantaneous velocities and turbulence levels in the centre of pools. Patterns of sediment deposition support the notion that stream competence is higher in the pool than the downstream riffle. Flow convergence around a large channel constriction appears to play a major role in multiple processes that include helical flow development and sediment routing, and backwater development with low velocities and turbulence levels above the constriction that may locally limit sediment supply. Jet flow, flow separation, vortex scour and turbulence generation enhance scour in the centre of pools. Ultimately, multiple processes appear to play some role in maintenance of this forced pool and the associated riffle. Copyright © 2008 John Wiley & Sons, Ltd.     

Application of Two Turbulence Closure Schemes in the Modelling of Tidal Currents and Salinity in the Hooghly Estuary

This paper has two purposes. The first is to study the circulation and salinity in Hooghly Estuary, along the east coast of India and the second is to compare the performance of two turbulence closure schemes by modelling it. A breadth averaged numerical model using a sigma co-ordinate system in the vertical is briefly described. Vertical diffusion of momentum and salt are parameterized by a simple first-order turbulent closure or by a one equation model for turbulent kinetic energy (TKE) which uses a specified mixing length. The results are compared with the available neap and spring tide observations along the estuary for both low and high discharge periods.The computed elevations and currents are in reasonable agreement with the observations showing no major differences in vertical current profiles by both the turbulent schemes. However, there is a slight under-prediction of bottom currents. The salinity profiles predicted by TKE model show better matching with observations. Statistical tests are also conducted to study the comparative performance of the turbulent closure schemes. The maintenance of two layer structure in residual currents and salt variability are also studied by using the model.