Application of altimetry data assimilation on mesoscale eddies simulation

Mesoscale eddy plays an important role in the ocean circulation. In order to improve the simulation accuracy of the mesoscale eddies, a three-dimensional variation (3DVAR) data assimilation system called Ocean Variational Analysis System (OVALS) is coupled with a POM model to simulate the mesoscale eddies in the Northwest Pacific Ocean. In this system, the sea surface height anomaly (SSHA) data by satellite altimeters are assimilated and translated into pseudo temperature and salinity (T-S) profile data. Then, these profile data are taken as observation data to be assimilated again and produce the three-dimensional analysis T-S field. According to the characteristics of mesoscale eddy, the most appropriate assimilation parameters are set up and testified in this system. A ten years mesoscale eddies simulation and comparison experiment is made, which includes two schemes: assimilation and non-assimilation. The results of comparison between two schemes and the observation show that the simulation accuracy of the assimilation scheme is much better than that of non-assimilation, which verified that the altimetry data assimilation method can improve the simulation accuracy of the mesoscale dramatically and indicates that it is possible to use this system on the forecast of mesoscale eddies in the future.

Anti-cyclonic eddies northwest of Luzon in summer-fall observed by satellite altimeters

Anti-cyclonic eddies northwest of Luzon of the Philippines in summer-fall are identified in the merged data products of satellite altimeters of Topex/Poseidon, Jason-1 and European Research Satellites. The generation and propagation of the anti-cyclonic eddies, which are confirmed by satellite ocean color data, are found to be a seasonal phenomenon that is phase-locked to the onset of the southwesterly monsoon and the relaxation of the cyclonic wind curl in the northeastern South China Sea. The eddies originate from northwest of Luzon in summer, move across the northeastern South China Sea to reach the China continental slope in fall, and propagate southwestward along the continental slope in fall-winter, inducing shelfbreak current variations in the western South China Sea in fall-winter. The anti-cyclonic eddy discovered by Li et al. (1998) in the northern South China Sea is found to originate from northwest of Luzon and carry primarily the South China Sea waters. It does not appear to be an eddy shed from the Kuroshio in the Luzon Strait area as alluded by Li et al. (1998) and others.

Variability of surface circulation in the South China Sea from satellite altimeter data

11-year satellite altimeter sea surface height (SSH) anomaly data from January 1993 to December 2003 are used to present the dominant spatial patterns and temporal variations of the South China Sea (SCS) surface circulation through Empirical Orthogonal Function (EOF) analysis. The first three EOF modes show the obvious seasonal variations of SSH in the SCS. EOF mode one is generally characterized by a basin-wide circulation. Mode two describes the double-cell basin scale circulation structure. The two cells were located off west of the Luzon Island and southeast of Vietnam, respectively. EOF mode three presents the mesoscale eddy structure in the western SCS, which develops into a strong cyclonic eddy rapidly from July to September. EOF mode one and mode three are also embedded with interannual signals, indicating that the SCS surface circulation variation is influenced by El Nino events prominently. The strong El Nino of 1997/98 obviously changed the SCS circulation structure. This study also shows that there existed a series of mesoscale eddies in the western SCS, and their temporal variation indicates intra-seasonal and interannual signals.

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.

冬季东中国海环流及中尺度涡旋数值模拟

该文在POM模式基础上,建立一个б坐标系下的三维斜压预报模式,考虑了海底地形、外来流、长江径流、海面风应力、海面热通量等多方面因素的影响,模拟结果较好地体现了冬季东中国海环流及涡旋的特征.在此数值模拟的基础上,分别在改变风场、以及黑潮流量的条件下进行数值实验,以考察它们对冬季东中国海环流和中尺度涡旋的影响.数值模拟和数值实验的结果表明,黑潮在流经东海时沿东海陆坡流动,在其两侧出现一些涡旋.冬季台湾暖流的表层水主要来自台湾东北的黑潮水,少部分来自台湾海峡;次表层,台湾东北的黑潮水和台湾海峡的入流对台湾暖流的贡献相当;深层,台湾暖流主要来源于台湾东北的黑潮水.冬季对马暖流的主要来源是黑潮分支和台湾暖流的分支,朝鲜沿岸流也是对马暖流的来源之一.

大洋中尺度涡旋与源区黑潮的相互作用研究

黑潮是西北太平洋一支高温、高盐、高流速的西边界流，它起源于菲律宾以东海域，沿吕宋岛、台湾东岸、东海大陆架外缘流动，并穿越吐噶喇海峡返回太平洋。其中吕宋岛至台湾岛这段西边界流命名为源区黑潮。源区黑潮在流经吕宋海峡以东时，由于失去岸界支撑会发生形变，并通过吕宋海峡与南海环流体系进行质量、动量和能量的交换，并影响着南海的环流特征。因此研究黑潮通过吕宋海峡入侵南海的关键动力过程和机制，对我们国家的军事、生产、环境等具有十分重大的意义，一直是国内外海洋学家关注的焦点。 本文利用高度计资料和涡旋精度的高分辨率数值模式（OFES）输出结果分析了南海东北部和吕宋海峡东侧海域的上层环流特征，并独立建立1.5层浅水约化重力模式下的理想模型，分析了黑潮从吕宋海峡入侵南海的非线性特征和西传的大洋中尺度涡旋与黑潮在吕宋海峡处的相互作用。 对高度计资料和模式结果资料分析的结果表明季风是影响南海北部上层环流的主要因素，黑潮以各种方式入侵南海也是影响南海北部环流特征的重要组成部分。通过数值模拟，我们得到在不同的参数下，黑潮入侵南海会存在稳定态、分叉、周期解等非线性特征，对应着黑潮南海流套、分离流环、无入侵等多种流态。通过讨论大洋中尺度涡旋与黑潮在吕宋海峡处的相互作用，发现涡旋在西传遇到黑潮后基本局地耗散或随黑潮继续向北移动，但当黑潮较弱时，涡旋可能将大洋中的能量通过吕宋海峡传入南海。

THE VARIATION AND MOVEMENT OF A KIND OF EDDIES IN SHELF SLOPE WATER

Considering the characteristics of the time and space scales of the eddies we established a quasi-static and quasi-geostrophic model to describe their variation and movement in shelf slope water. The analytical solution revealed the main properties of the variation: slow expansion and fast stagnation processes and the law of the eddy motion affected under the background field. All theoretical results are proved by satellite image measurements.

Close coupling between phytoplankton growth and microzooplankton grazing in the western South China Sea

We conducted 28 dilution experiments during August-September 2007 to investigate the coupling of growth and microzooplankton grazing rates among ultraphytoplankton populations and the phytoplankton community and their responses to habitat variability (open-ocean oligotrophy, eddy-induced upwelling, and the Mekong River plume) in the western South China Sea. At the community level, standing stocks, growth, and grazing rates were strongly and positively correlated, and were related to the higher abundance of larger phytoplankton cells (diatoms) at stations with elevated chlorophyll concentration. Phytoplankton growth rates were highest (> 2 d(-1)) within an eastward offshore jet at 13 degrees N and at a station influenced by the river plume. Among ultraphytoplankton populations, Prochlorococcus dominated the more oceanic and oligotrophic stations characterized by generally lower biomass and phytoplankton community growth, whereas Synechococcus became more important in mesotrophic areas (eddies, offshore jet, and river plume). The shift to Synechococcus dominance reflected, in part, its higher growth rates (0.87 +/- 0.45 d(-1)) compared to Prochlorococcus (0.65 +/- 0.29 d(-1)) or picophytoeukaryotes (0.54 +/- 0.50 d(-1)). However, close coupling of microbial mortality rates via common predators is seen to play a major role in driving the dominance transition as a replacement of Prochlorococcus, rather than an overprinting of its steady-state standing stock.

Mesofracture mechanics: a necessary link

Classical fracture mechanics is based on the premise that small scale features could be averaged to give a larger scale property such that the assumption of material homogeneity would hold. Involvement of the material microstructure, however, necessitates different characteristic lengths for describing different geometric features. Macroscopic parameters could not be freely exchanged with those at the microscopic scale level. Such a practice could cause misinterpretation of test data. Ambiguities arising from the lack of a more precise range of limitations for the definitions of physical parameters are discussed in connection with material length scales. Physical events overlooked between the macroscopic and microscopic scale could be the link that is needed to bridge the gap. The classical models for the creation of free surface for a liquid and solid are oversimplified. They consider only the translational motion of individual atoms. Movements of groups or clusters of molecules deserve attention. Multiscale cracking behavior also requires the distinction of material damage involving at least two different scales in a single simulation. In this connection, special attention should be given to the use of asymptotic solution in contrast to the full field solution when applying fracture criteria. The former may leave out detail features that would have otherwise been included by the latter. Illustrations are provided for predicting the crack initiation sites of piezoceramics. No definite conclusions can be drawn from the atomistic simulation models such as those used in molecular dynamics until the non-equilibrium boundary conditions can be better understood. The specification of strain rates and temperatures should be synchronized as the specimen size is reduced to microns. Many of the results obtained at the atomic scale should be first identified with those at the mesoscale before they are assumed to be connected with macroscopic observations. Hopefully, "mesofracture mechanics" could serve as the link to bring macrofracture mechanics closer to microfracture mechanics.

Numerical Study Of Heat Transfer Mechanism In Turbulent Supercritical Co2 Channel Flow

Direct numerical simulation (DNS) of supercritical CO2 turbulent channel flow has been performed to investigate the heat transfer mechanism of supercritical fluid. In the present DNS, full compressible Navier-Stokes equations and Peng-Robison state equation are solved. Due to effects of the mean density variation in the wall normal direction, mean velocity in the cooling region becomes high compared with that in the heating region. The mean width between high-and low-speed streaks near the wall decreases in the cooling region, which means that turbulence in the cooling region is enhanced and lots of fine scale eddies are created due to the local high Reynolds number effects. From the turbulent kinetic energy budget, it is found that compressibility effects related with pressure fluctuation and dilatation of velocity fluctuation can be ignored even for supercritical condition. However, the effect of density fluctuation on turbulent kinetic energy cannot be ignored. In the cooling region, low kinematic viscosity and high thermal conductivity in the low speed streaks modify fine scale structure and turbulent transport of temperature, which results in high Nusselt number in the cooling condition of the supercritical CO2.

Microstructure, deformation and failure of polymer bonded explosives

Polymer bonded explosives (PBXs) are highly particle filled composite materials comprised of explosive crystals and a polymeric binder (ca. 5-10% by weight). The microstructure and mechanical properties of two pressed PBXs with different binder systems were studied in this paper. The initial microstructure of the pressed PBXs and its evolution under different mechanical aggressions were studied, including quasi-static tension and compression, ultrasonic wave stressing and long-pulse low-velocity impact. Real-time microscopic observation of the PBXs under tension was conducted by using a scanning electron microscope equipped with a loading stage. The mechanical properties under tensile creep, quasi-static tension and compression were studied. The Brazilian test, or diametrical compression, was used to study the tensile properties. The influences of pressing pressures and temperatures, and strain rates on the mechanical properties of PBXs were analyzed. The mesoscale damage modes in initial pressed samples and the samples insulted by different mechanical aggressions, and the corresponding failure mechanisms of the PBXs under different loading conditions were analyzed.

Multiscale Method for Turbulent-Flow Computations

A main method of predicting turbulent flows is to solve LES equations, which was called traditional LES method. The traditional LES method solves the motions of large eddies of size larger than filtering scale An while modeling unresolved scales less than Delta_n. Hughes et al argued that many shortcomings of the traditional LES approaches were associated with their inabilities to successfully differentiate between large and small scales. One may guess that a priori scale-separation would be better, because it can predict scale-interaction well compared with posteriori scale-separation. To this end, a multi-scale method was suggested to perform scale-separation computation. The primary contents of the multiscale method are l) A space average is used to differentiate scale. 2) The basic equations include the large scale equations and fluctuation equations. 3) The large-scale equations and fluctuation equations are coupled through turbulent stress terms. We use the multiscale equations of n=2, i.e., the large and small scale (LSS) equations, to simulate 3-D evolutions of a channel flow and a planar mixing layer flow Some interesting results are given.

Numerical Study of Heat Transfer Mechanism in Turbulent Supercritical CO2 Channel Flow

Direct numerical simulation (DNS) of supercritical CO2 turbulent channel flow has been performed to investigate the heat transfer mechanism of supercritical fluid. In the present DNS, full compressible Navier-Stokes equations and Peng-Robison state equation are solved. Due to effects of the mean density variation in the wall normal direction, mean velocity in the cooling region becomes high compared with that in the heating region. The mean width between high-and low-speed streaks near the wall decreases in the cooling region, which means that turbulence in the cooling region is enhanced and lots of fine scale eddies are created due to the local high Reynolds number effects. From the turbulent kinetic energy budget, it is found that compressibility effects related with pressure fluctuation and dilatation of velocity fluctuation can be ignored even for supercritical condition. However, the effect of density fluctuation on turbulent kinetic energy cannot be ignored. In the cooling region, low kinematic viscosity and high thermal conductivity in the low speed streaks modify fine scale structure and turbulent transport of temperature, which results in high Nusselt number in the cooling condition of the supercritical CO2.

Computation of turbulence-generated noise by large-eddy simulation

The hybrid method of large eddy simulation (LES) and the Lighthill analogy is being developed to compute the sound radiated from turbulent flows. The results obtained from the hybrid method are often contaminated by the absence of small scales in LES, since the energy level of sound is much smaller than that of turbulent flows. Previous researches investigate the effects of subgrid sacle (SGS) eddies on the frequency spectra of sound radiated by isotropic turbulence and suggest a SGS noise model to represent the SGS contributions to the frequency spectra. Their investigations are conducted in physical space and are unavoidably influenced by boundary conditions. In this paper, we propose to perform such calculations in Fourier space so that the effects of boundary conditions can be correctly treated. Posteriori tests are carried out to investigate the SGS contribution to the sound. The results obtained recover the -7/2 law within certain wave-number ranges, but under-estimate the amplitudes of the frequency spectra. The reason for the underestimation is also discussed.

Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence

Large-eddy simulation (LES) has emerged as a promising tool for simulating turbulent flows in general and, in recent years,has also been applied to the particle-laden turbulence with some success (Kassinos et al., 2007). The motion of inertial particles is much more complicated than fluid elements, and therefore, LES of turbulent flow laden with inertial particles encounters new challenges. In the conventional LES, only large-scale eddies are explicitly resolved and the effects of unresolved, small or subgrid scale (SGS) eddies on the large-scale eddies are modeled. The SGS turbulent flow field is not available. The effects of SGS turbulent velocity field on particle motion have been studied by Wang and Squires (1996), Armenio et al. (1999), Yamamoto et al. (2001), Shotorban and Mashayek (2006a,b), Fede and Simonin (2006), Berrouk et al. (2007), Bini and Jones (2008), and Pozorski and Apte (2009), amongst others. One contemporary method to include the effects of SGS eddies on inertial particle motions is to introduce a stochastic differential equation (SDE), that is, a Langevin stochastic equation to model the SGS fluid velocity seen by inertial particles (Fede et al., 2006; Shotorban and Mashayek, 2006a; Shotorban and Mashayek, 2006b; Berrouk et al., 2007; Bini and Jones, 2008; Pozorski and Apte, 2009).However, the accuracy of such a Langevin equation model depends primarily on the prescription of the SGS fluid velocity autocorrelation time seen by an inertial particle or the inertial particle–SGS eddy interaction timescale (denoted by $\delt T_{Lp}$ and a second model constant in the diffusion term which controls the intensity of the random force received by an inertial particle (denoted by C_0, see Eq. (7)). From the theoretical point of view, dTLp differs significantly from the Lagrangian fluid velocity correlation time (Reeks, 1977; Wang and Stock, 1993), and this carries the essential nonlinearity in the statistical modeling of particle motion. dTLp and C0 may depend on the filter width and particle Stokes number even for a given turbulent flow. In previous studies, dTLp is modeled either by the fluid SGS Lagrangian timescale (Fede et al., 2006; Shotorban and Mashayek, 2006b; Pozorski and Apte, 2009; Bini and Jones, 2008) or by a simple extension of the timescale obtained from the full flow field (Berrouk et al., 2007). In this work, we shall study the subtle and on-monotonic dependence of $\delt T_{Lp}$ on the filter width and particle Stokes number using a flow field obtained from Direct Numerical Simulation (DNS). We then propose an empirical closure model for $\delta T_{Lp}$. Finally, the model is validated against LES of particle-laden turbulence in predicting single-particle statistics such as particle kinetic energy. As a first step, we consider the particle motion under the one-way coupling assumption in isotropic turbulent flow and neglect the gravitational settling effect. The one-way coupling assumption is only valid for low particle mass loading.