Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been developed to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

The variability of eddy kinetic energy in the South China Sea deduced from satellite altimeter data

We used fifteen years (1993-2007) of altimetric data, combined from different missions (ERS-1/2, TOPEX/Poseidon, Jason-1, and Envisat), to analyze the variability of the eddy kinetic energy (EKE) in the South China Sea (SCS). We found that the EKE ranged from 64 cm(2)/s(2) to 1 390 cm(2)/s(2) with a mean value of 314 cm(2)/s(2). The highest EKE center was observed to the east of Vietnam (with a mean value of 509 cm(2)/s(2)) and the second highest EKE region was located to the southwest of Taiwan Island (with a mean value of 319 cm(2)/s(2)). We also found that the EKE structure is the consequence of the superposition of different variability components. First, interannual variability is important in the SCS. Spectral analysis of the EKE interannual signal (IA-EKE) shows that the main periodicities of the IA-EKE to the east of Vietnam, to the southwest of Taiwan Island, and in the SCS are 3.75, 1.87, and 3.75 years, respectively. It is to the south of Taiwan Island that the IA-EKE signal has the most obvious impact on EKE variability. In addition, the IA-EKE exhibit different trends in different regions. An obvious positive trend is observed along the east coast of Vietnam, while a negative trend is found to the southwest of Taiwan Island and in the east basin of Vietnam. Correlation analysis shows that the IA-EKE has an obvious negative correlation with the SSTA in Nio3 (5A degrees S-5A degrees N, 90A degrees W-150A degrees W). El Nio-Southern Oscillation (ENSO) affects the IA-EKE variability in the SCS through an atmospheric bridge-wind stress curl over the SCS. Second, the seasonal cycle is the most obvious timescale affecting EKE variability. The locations of the most remarkable EKE seasonal variabilities in the SCS are to the east of Vietnam, to the southwest of Taiwan, and to the west of Philippines. To the east of Vietnam, the seasonal cycle is the dominant mechanism controlling EKE variability, which is attributed primarily to the annual cycle there of wind stress curl. In this area, the maximum EKE is observed in autumn. To the southwest of Taiwan Island, the EKE is enlarged by the stronger SCS circulation, which is caused by the intrusion branch from the Kuroshio in winter. Finally, intra-annual and mesoscale variability, although less important than the former, cannot be neglected. The most obvious intra-annual and mesoscale variability, which may be the result of baroclinic instability of the background flow, are observed to the southwest of Taiwan Island. Sporadic events can have an important effect on EKE variability.

Eddy generation and evolution in the North Pacific Subtropical Countercurrent (NPSC) zone

The seasonal generation and evolution of eddies in the region of the North Pacific Subtropical Countercurrent remain poorly understood due to the scarcity of available data. We used TOPEX/POSEIDON altimetry data from 1992 to 2007 to study the eddy field in this zone. We found that velocity shear between this region and the neighboring North Equatorial Current contributes greatly to the eddy generation. Furthermore, the eddy kinetic energy level (EKE) shows an annual cycle, maximum in April/May and minimum in December/January. Analyses of the temporal and spatial distributions of the eddy field revealed clearly that the velocity shear closely related to baroclinic instability processes. The eddy field seems to be more zonal than meridional, and the energy containing length scale shows a surprising lag of 2-3 months in comparison with the 1-D and 2-D EKE level. A similar phenomenon is observed in individual eddies in this zone. The results show that in this eddy field band, the velocity shear may drive the EKE level change so that the eddy field takes another 2-3 months to grow and interact to reach a relatively stable state. This explains the seasonal evolution of identifiable eddies.

Vertical structure and evolution of the Luzon Warm Eddy

Eddies are frequently observed in the northeastern South China Sea (SCS). However, there have been few studies on vertical structure and temporal-spatial evolution of these eddies. We analyzed the seasonal Luzon Warm Eddy (LWE) based on Argo float data and the merged data products of satellite altimeters of Topex/Poseidon, Jason-1 and European Research Satellites. The analysis shows that the LWE extends vertically to more than 500 m water depth, with a higher temperature anomaly of 5A degrees C and lower salinity anomaly of 0.5 near the thermocline. The current speeds of the LWE are stronger in its uppermost 200 m, with a maximum speed of 0.6 m/s. Sometimes the LWE incorporates mixed waters from the Kuroshio Current and the SCS, and thus has higher thermohaline characteristics than local marine waters. Time series of eddy kinematic parameters show that the radii and shape of the LWE vary during propagation, and its eddy kinetic energy follows a normal distribution. In addition, we used the empirical orthogonal function (EOF) here to analyze seasonal characteristics of the LWE. The results suggest that the LWE generally forms in July, intensifies in August and September, separates from the coast of Luzon in October and propagates westward, and weakens in December and disappears in February. The LWE's westward migration is approximately along 19A degrees N latitude from northwest of Luzon to southeast of Hainan, with a mean speed of 6.6 cm/s.

A Numerical Study for Turbulent Flow and Thermal Influence Over Inhomogenous Canopy of Roughness Elements

A large-eddy simulation with transitional structure function(TSF) subgrid model we previously proposed was performed to investigate the turbulent flow with thermal influence over an inhomogeneous canopy, which was represented as alternative large and small roughness elements. The aerodynamic and thermodynamic effects of the presence of a layer of large roughness elements were modelled by adding a drag term to the three-dimensional Navier-Stokes equations and a heat source/sink term to the scalar equation, respectively. The layer of small roughness elements was simply treated using the method as described in paper (Moeng 1984, J. Atmos Sci. 41, 2052-2062) for homogeneous rough surface. The horizontally averaged statistics such as mean vertical profiles of wind velocity, air temperature, et al., are in reasonable agreement with Gao et al.(1989, Boundary layer meteorol. 47, 349-377) field observation (homogeneous canopy). Not surprisingly, the calculated instantaneous velocity and temperature fields show that the roughness elements considerably changed the turbulent structure within the canopy. The adjustment of the mean vertical profiles of velocity and temperature was studied, which was found qualitatively comparable with Belcher et al. (2003, J Fluid Mech. 488, 369-398)'s theoretical results. The urban heat island(UHI) was investigated imposing heat source in the region of large roughness elements. An elevated inversion layer, a phenomenon often observed in the urban area (Sang et al., J Wind Eng. Ind. Aesodyn. 87, 243-258)'s was successfully simulated above the canopy. The cool island(CI) was also investigated imposing heat sink to simply model the evaporation of plant canopy. An inversion layer was found very stable and robust within the canopy.

Multiscale Coupling in Complex Mechanical Systems

Multiscale coupling attracts broad interests from mechanics, physics and chemistry to biology. The diversity and coupling of physics at different scales are two essential features of multiscale problems in far-from-equilibrium systems. The two features present fundamental difficulties and are great challenges to multiscale modeling and simulation. The theory of dynamical system and statistical mechanics provide fundamental tools for the multiscale coupling problems. The paper presents some closed multiscale formulations, e.g., the mapping closure approximation, multiscale large-eddy simulation and statistical mesoscopic damage mechanics, for two typical multiscale coupling problems in mechanics, that is, turbulence in fluids and failure in solids. It is pointed that developing a tractable, closed nonequilibrium statistical theory may be an effective approach to deal with the multiscale coupling problems. Some common characteristics of the statistical theory are discussed.

Derivation Of A Nonlinear Reynolds Stress Model Using Renormalization Group Analysis And Two-Scale Expansion Technique

Adopting Yoshizawa's two-scale expansion technique, the fluctuating field is expanded around the isotropic field. The renormalization group method is applied for calculating the covariance of the fluctuating field at the lower order expansion. A nonlinear Reynolds stress model is derived and the turbulent constants inside are evaluated analytically. Compared with the two-scale direct interaction approximation analysis for turbulent shear flows proposed by Yoshizawa, the calculation is much more simple. The analytical model presented here is close to the Speziale model, which is widely applied in the numerical simulations for the complex turbulent flows.

Short-and resonant-range interactions between scales in turbulence and their applications

Interactions between different scales in turbulence were studied starting from the incompressible Navier-Stokes equations. The integral and differential formulae of the short-range viscous stresses, which express the short-range interactions between contiguous scales in turbulence, were given. A concept of the resonant-range interactions between extreme contiguous scales was introduced and the differential formula of the resonant-range viscous stresses was obtained. The short- and resonant-range viscous stresses were applied to deduce the large-eddy simulation (LES) equations as well as the multiscale equations, which are approximately closed and do not contain any empirical constants or relations. The properties and advantages of using the multiscale equations to compute turbulent flows were discussed. The short-range character of the interactions between the scales in turbulence means that the multiscale simulation is a very valuable technique for the calculation of turbulent flows. A few numerical examples were also given.

关于Maunder极小期太阳活动的统计分析

本文通过自相关分析探讨在“Maunder极小期”期间太阳活动是否仍然存在11年活动周期.利用Eddy根据欧洲望远镜在“Maunder极小期”时观测的太阳黑子记录所作的相对黑子数年平均值进行了自相关分析,也与中国地方志在“Maunder极小期”时的有关黑子记录一起进行了同样的自相关分析.在作自相关分析时采用几种不同的方法进行比较,结果发现,11年活动周期性的存在与否与所选用的统计方法有关,因此不能绝对地说在“Maunder极小期”时太阳活动完全停止了,或者说在“Maunder极小期”时仍然存在着11年周期的活动性质.

An improvement on the theoretical model for the large-scale coherent structures in the outer region of a turbulent boundary layer

According to the experimental results, there exist large-scale coherent structures in the outer region of a turbulent boundary layer, which have been studied by many authors.As experimental results, Antonia (1990) showed the phase- aver aged streamlines and isovorticity lines of the large-scale coherent structures in a turbulent boundary layer for different Reynolds numbers. Based on the hydrodynamic stability theory, the 2-D theoretical model for the large-scale structures was proposed by Luo and Zhou, in which the eddy viscosity was defined as a complex function of the position in the normal direction. The theoretical results showed in ref. were in agreement with those in ref. However, there were two problems in the results. One is that in the experimental results, there were divergent focuses between two saddle points in the streamlines, but in the theoretical results, there were centers. The other is that the stretched parts of the isovorticity lines appear at the location of centers in the theoretical results, while in the experimental results they located somewhere between the focuses and saddle points. The reason is that the computations were based on a 2-D model.

Assessment of bursting detection methods with LES database of turbulent channel flows

Validated by comparison with DNS, numerical database of turbulent channel flows is yielded by Large Eddy Simulation (LES). Three conventional techniques: uv quadrant 2, VITA and mu-level techniques for detecting turbulent bursts are applied to the identification of turbulent bursts. With a grouping parameter introduced by Bogard & Tiedemann (1986) or Luchik & Tiederman (1987), multiple ejections detected by these techniques which originate from a single burst can be grouped into a single-burst event. The results are compared with experimental results, showing that all techniques yield reasonable average burst period. However, uv quadrant 2 and mu-level are found to be superior to VITA in having large threshold-independent range.

Statistics of Turbulence in Benthal Boundary Layers

This paper deals with turbulence behavior inbenthalboundarylayers by means of large eddy simulation (LES). The flow is modeled by moving an infinite plate in an otherwise quiescent water with an oscillatory and a steady velocity components. The oscillatory one aims to simulate wave effect on the flow. A number of large-scale turbulence databases have been established, based on which we have obtained turbulencestatisticsof the boundarylayers, such as Reynolds stress, turbulence intensity, skewness and flatness ofturbulence, and temporal and spatial scales of turbulent bursts, etc. Particular attention is paid to the dependences of those statistics on two nondimensional parameters, namely the Reynolds number and the current-wave velocity ratio defined as the steady current velocity over the oscillatory velocity amplitude. It is found that the Reynolds stress and turbulence intensity profile differently from phase to phase, and exhibit two types of distributions in an oscillatory cycle. One is monotonic occurring during the time when current and wave-induced components are in the same direction, and the other inflectional occurring during the time when current and wave-induced components are in opposite directions. Current component makes an asymmetrical time series of Reynolds stress, as well as turbulence intensity, although the mean velocity series is symmetrical as a sine/cosine function. The skewness and flatness variations suggest that the turbulence distribution is not a normal function but approaches to a normal one with the increasing of Reynolds number and the current-wave velocity ratio as well. As for turbulent bursting, the dimensionless period and the mean area of all bursts per unit bed area tend to increase with Reynolds number and current-wave velocity ratio, rather than being constant as in steady channel flows.

DNS and LES of turbulent channel flow with hydrophobic surface

Hydrophobic surface benefits for drag reduction. Min and Kim[1] do the first Direct Numerical Simulation on drag reduction in turbulent channel flow. And Fukagata and Kasagi[2] make some theoretical analysis based on Dean[3]'s formula and some observations in the DNS results. Using their theory, they conclude that drag reduction is possible in large Reynolds number. Both Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) are performed in our research. How the LES behaving in the turbulent channel flow with hydrophobic surface is examined. Original Smagorinsky model and its Dynamical model are used in LES. The slip velocities predicted by LES using Dynamical model are in good agreement with DNS as shown in the Figure. Although the percentage of drag reduction predicted by LES shows some discrepancies, it is in the error limit for industrial flow. First order and second order moments of LES are also examined and compared with DNS's results. The first-order moments is calculated well by LES. But there are some discrepancies of second-order moments between LES and DNS. [GRAPHICS]