参数激励圆柱形容器中的非线性Faraday波

在柱坐标系下，通过奇异摄动理论的多尺度展开法求解势流方程，研究了垂直强迫激励圆柱形容器中的单一模式水表面驻波模式．假设流体是无粘、不可压且运动是无旋的，在忽略了表面张力的影响下，用两变量时间展开法得到一个具有立方项以及底部驱动项影响的非线性振幅方程．对上述方程进行了数值计算，计算的结果显示了在不同驱动振幅和驱动频率下，会激发不同自由水表面驻波模式，从等高线的图像来看，和以往的实验结果相当吻合．

Application of the characteristic CIP method to a shallow water model on the sphere

Semi-implicit algorithms are popularly used to deal with the gravitational term in numerical models. In this paper, we adopt the method of characteristics to compute the solutions for gravity waves on a sphere directly using a semi-Lagrangian advection scheme instead of the semi-implicit method in a shallow water model, to avoid expensive matrix inversions. Adoption of the semi-Lagrangian scheme renders the numerical model always stable for any Courant number, and which saves CPU time. To illustrate the efficiency of the characteristic constrained interpolation profile (CIP) method, some numerical results are shown for idealized test cases on a sphere in the Yin-Yang grid system.

Influence of the earth rotation on the interfacial wave solutions and wave-induced tangential stress in a two-layer fluid system

Interfacial waves and wave-induced tangential stress are studied for geostrophic small amplitude waves of two-layer fluid with a top free surface and a flat bottom. The solutions were deduced from the general form of linear fluid dynamic equations of two-layer fluid under the f-plane approximation, and wave-induced tangential stress were estimated based on the solutions obtained. As expected; the solutions derived from the present work include as special cases those obtained by Sun et al. (2004. Science in China, Set. D, 47(12): 1147-1154) for geostrophic small amplitude surface wave solutions and wave-induced tangential stress if tire density of the upper layer is much smaller than that of the lower layer. The results show that the interface and the surface will oscillate synchronously, and the influence of the earth's rotation both on the surface wave solutions and the interfacial wave solutions should be considered.

Ocean wave imaging mechanism by imaging radar

Analytical representations of the high frequency spectra of ocean wave and its variation due to the variation of ocean surface current are derived from the wave-number spectrum balance equation. The ocean surface imaging formulation of real aperture radar (RAR) is given using electromagnetic wave backscattering theory of ocean surface and the modulations of ocean surface winds, currents and their variations to RAR are described. A general representation of the phase modulation induced by the ocean surface motion is derived according to standard synthetic aperture radar (SAR) imaging theory. The detectability of ocean current and sea bottom topography by imaging radar is discussed. The results constitute the theoretical basis for detecting ocean wave fields, ocean surface winds, ocean surface current fields, sea bottom topography, internal wave and so on.

Effect of Stokes drift on upper ocean mixing

Stokes drift is the main source of vertical vorticity in the ocean mixed layer. In the ways of Coriolis - Stokes forcing and Langmuir circulations, Stokes drift can substantially affect the whole mixed layer. A modified Mellor-Yamada 2.5 level turbulence closure model is used to parameterize its effect on upper ocean mixing conventionally. Results show that comparing surface heating with wave breaking, Stokes drift plays the most important role in the entire ocean mixed layer, especially in the subsurface layer. As expected, Stokes drift elevates both the dissipation rate and the turbulence energy in the upper ocean mixing. Also, influence of the surface heating, wave breaking and wind speed on Stokes drift is investigated respectively. Research shows that it is significant and important to assessing the Stokes drift into ocean mixed layer studying. The laboratory observations are supporting numerical experiments quantitatively.

A new method for the estimation of oceanic mixed-layer depth using shipboard X-band radar images

Shipboard X-band radar images acquired on 24 June 2009 are used to study nonlinear internal wave characteristics in the northeastern South China Sea. The studied images show three nonlinear internal waves in a packet. A method based on the Radon Transform technique is introduced to calculate internal wave parameters such as the direction of propagation and internal wave velocity from backscatter images. Assuming that the ocean is a two-layer finite depth system, we can derive the mixed-layer depth by applying the internal wave velocity to the mixed-layer depth formula. Results show reasonably good agreement with in-situ thermistor chain and conductivity-temperature-depth data sets.

Probability distribution of random wave-current forces

Based on the second-order solutions obtained for the three-dimensional weakly nonlinear random waves propagating over a steady uniform current in finite water depth, the joint statistical distribution of the velocity and acceleration of the fluid particle in the current direction is derived using the characteristic function expansion method. From the joint distribution and the Morison equation, the theoretical distributions of drag forces, inertia forces and total random forces caused by waves propagating over a steady uniform current are determined. The distribution of inertia forces is Gaussian as that derived using the linear wave model, whereas the distributions of drag forces and total random forces deviate slightly from those derived utilizing the linear wave model. The distributions presented can be determined by the wave number spectrum of ocean waves, current speed and the second order wave-wave and wave-current interactions. As an illustrative example, for fully developed deep ocean waves, the parameters appeared in the distributions near still water level are calculated for various wind speeds and current speeds by using Donelan-Pierson-Banner spectrum and the effects of the current and the nonlinearity of ocean waves on the distribution are studied. (c) 2006 Elsevier Ltd. All rights reserved.

澳大利亚阿德莱德地区大气重力波的气辉观测研究

气辉是研究中高层大气及电离层高度的重力波的传播特性的一种重要手段。本文利用澳大利亚阿德莱德地区一个三镜头气辉仪7年的OH和OI气辉观测数据，结合同一地区中频雷达所测日平均风场数据，采用最大熵谱分析方法，分析研究了该地区上空中层顶区域重力波的统计传播特性，并研究了重力波在中层顶区域的垂直传播特征。结果表明: （1）在中层顶气辉观测区，通过OH和OI分别探测到49个和28个周期大于40分钟的重力波事件。统计结果表明，在OH和OI两个高度上观测到的较大尺度重力波主要集中在周期为68分钟，水平相速度为36米/秒。同时这些重力波接近水平方向传播，具有东南向的优势传播方向。OH仰角平均为10度，OI为8度。表明传播到中层顶气辉层区域的这些重力波可能产生于相距观测点很远的地方，并且受到背景风场的方向滤波和拉伸的作用。 （2）在1999年1月12日晚上的气辉观测中， 我们通过OH和OI气辉在两个不同高度上观测到了同一次重力波事件，根据两个高度上观测到的扰动时序相位差， 我们判断该重力波为上行波。然而，我们所得到的大部分事例中，符合在两个高度上观测到的重力波事例之间具有明显相关性的极少。分析表明，由于该气辉仪及其观测方法本身存在固有的误差，同时也因为重力波在中层顶高度受到强烈的风剪切作用而导致传播相位的改变，使得我们很少能在两个气辉探测高度上观测到同一个重力波事件。

电离层行进式扰动的GPS台网监测研究

In this dissertation, we investigated two types of traveling ionospheric disturbances (TIDs)/gravity waves (GWs) triggered separately by auroral energy input during super geomagnetic storms and solar terminator (ST) under quiet geomagnetic conditions (kp<3+) using TEC measurements from the global network of GPS receivers. Research into the generation and propagation of TIDs/GWs during storms greatly enhance our understandings on the evolution processes of energy transportation from the high-latitude’s magnetosphere to the low-latitude ionosphere and the conjugated effect of TIDs propagation between the northern and southern hemispheres. Our results revealed that the conjugacy of propagation direction between the northern and southern hemispheres was subject to the influence of Coriolis force. We also figure out the evolution processes of ionospheric disturbances at the global scale. These are important topics that had not been well addressed previously. In addition, we also obtained thee wave structures of medium scale TIDs excited by the solar terminator (ST) moving over the northern America and physical mechanisms involved. Our observations confirm that the ST is a stable and repetitive source of ionospheric wave disturbances and the evidence of solar terminator generated disturbances has been demonstrated experimentally via the GPS TEC measurement. The main researches and results of this dissertation are as follows. First, the global traveling ionospheric disturbances (TIDs) during the drastic magnetic storms of October 29–31, 2003 were analyzed using the Global Position System (GPS) total electron content (TEC) data observed in the Asian-Australian, European and North American sectors. We collected the most comprehensive set of the TEC data from more than 900 GPS stations on the International GNSS Services (IGS) website and introduce here a strategy that combines polynomial fitting and multi-channel maximum entropy spectral analysis to obtain TID parameters. Moreover, in collaboration with my thesis advisor, I have developed an imaging technique of 2-dimensional map of TIDs structures to obtain spatial and temporal maps of large scale traveling ionospheric disturbances (LSTIDs). The clear structures of TEC perturbations map during the passage of TIDs were displayed. The results of our study are summarized as follows: (1) Large-scale TIDs (LSTIDs) and medium-scale TIDs (MSTIDs) were detected in all three sectors after the sudden commencement (SC) of the magnetic storm, and their features showed longitudinal and latitudinal dependences. The duration of TIDs was longer at higher latitudes than at middle latitudes, with a maximum of about 16 h. The TEC variation amplitude of LSTIDs was larger in the North American sector than in the two other sectors. At the lower latitudes, the ionospheric perturbations were more complicated, and their duration and amplitude were relatively longer and larger. (2) The periods and phase speeds of TIDs were different in these three sectors. In Europe, the TIDs propagated southward; in North America and Asia, the TIDs propagated southwestward; in the near-equator region, the disturbances propagated with the azimuth (the angle of the propagation direction of the LSTIDs measured clockwise from due north with 0°) of 210° showing the influence of Coriolis force; in the Southern Hemisphere, the LSTIDs propagated conjugatedly northwestward. Both the southwestward and northeastward propagating LSTIDs are found in the equatorial region. These results mean that the Coriolis effect cannot be ignored for the wave propagation of LSTIDs and that the propagation direction is correlated with the polar magnetic activity. (3) The day (day of year: 301) before the SC (sudden commencement) of magnetic storm, we observed a sudden TEC skip disturbances (±10 TECU). It should be a response for the high flux of proton during the solar flare event, but not the magnetic storms. Next, the most comprehensive and dense GPS network’s data from North-America region were used in this paper to analyze the medium scale traveling ionospheric disturbances (MSTIDs) which were generated by the moving solar terminator during the quiet days in 2005. We applied the multi-channel maximum entropy spectral analysis to calculated TID parameters, and found that the occurrence of ST-MSTIDs depends on the seasonal variations. The results of our study are summarized as follows: (1) MSTIDs stimulated by the moving ST (ST-MSTIDs) are detected at mid-latitudes after the passage of the solar terminator with the life time of 2～3 hours and the variation amplitude of 0.2～0.8 TECU. Spectral analysis indicated that the horizontal wavelength, average period, horizontal phase velocity of the MSTIDs are around 300±150 km，150±80 m/s and 25±15 min, respectively. In addition, ST-MSTIDs have wave fronts elongating the moving ST direction and almost parallel to ST. (2) The statistical results demonstrate that the dusk MSTIDs stimulated by ST is more obvious than the dawn MSTIDs in summer. On the contrary, the more-pronounced dawn MSTIDs occurs in winter. (3) Further analysis indicates that the seasonal variations of ST-MSTIDs occurrence frequency are most probably related to the seasonal differences of the variations of EUV flux in the ionosphere region and recombination process during sunrise and sunset period at mid-latitudes. Statistical study of occurrence characteristics of TIDs using the GPS network in North-American and European during solar maximum, In conclusion, statistical studies of the propagation characteristics of TIDs, which excited by the two common origins including geomagnetic storms and moving solar terminator, were involved with global GPS TEC databasein this thesis. We employed the multichannel maximum entropy spectral analysis method to diagnose the characteristics of propagation and evolvement of ionospheric disturbances, also, the characteristics of their regional distribution and climatological variations were revealed by the statistic analysis. The results of these studies can improve our knowledge about the energy transfer in the solar-terrestrial system and the coupling process between upper and lower atmosphere (thermosphere-ionosphere-mesosphere). On the other hand, our results of the investigation on TIDs generated by particular linear origin such as ST are important for developing ionospheric irregularity physics and modeling the transionosphere radio wave propagation. Besides, the GPS TEC representation of the ST-generated ionospheric structure suggests a better possibility for investigating this phenomenon. Subsequently, there are scientific meaning of the result of this dissertation to deeply discuss the energy transfer and coupling in the ionosphere, as well as realistic value to space weather forecast in the ionosphere region.

Transition waves and nonlinear interactions in the near wake of a circular cylinder

Transition waves and interactions between two kinds of instability-vortex shedding and transition wave in the near wake of a circular cylinder in the Reynolds number range 3 000-10 000 are studied by a domain decomposition hybrid numerical method. Based on high resolution power spectral analyses for velocity new results on the Reynolds-number dependence of the transition wave frequency, i.e. f(t)/f(s) similar to Re-0.87 are obtained. The new predictions are in good agreement with the experimental results of Wei and Smith but different from Braza's prediction and some early experimental results f(t)/f(s) similar to Re-0.5 given by Bloor et nl. The multi-interactions between two kinds of vortex are clearly visualized numerically. The strong nonlinear interactions between the two independent frequencies (f(t), f(s)) leading to spectra broadening to form the coupling mf(s) +/- nf(t) are predicted and analyzed numerically, and the characteristics of the transition are described. Longitudinal variations of the transition wave and its coupling are reported. Detailed mechanism of the flow transition in the near wake before occurrence of the three-dimensional evolution is provided.

Short-Term Nonlinear Response of Tension Leg Platform in Random Sea Waves

Most of the existing mathematical models for analyzing the dynamic response of TLP are based on explicit or implicit assumptions that motions (translations and rotations) are small magnitude. However, when TLP works in severe adverse conditions, the a priori assumption on small displacements may be inadequate. In such situation, the motions should be regarded as finite magnitude. This paper will study stochastic nonlinear dynamic responses of TLP with finite displacements in random waves. The nonlinearities considered are: large amplitude motions, coupling the six degrees-of-freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force. The nonlinear dynamic responses are calculated by using numerical integration procedure in the time domain. After the time histories of the dynamic responses are obtained, we carry out cycle counting of the stress histories of the tethers with rain-flow counting method to get the stress range distribution.

Highly nonlinear internal solitary waves over the continental shelf of the northwestern South China Sea

Large amplitude internal solitary waves (ISWs) often exhibit highly nonlinear effects and may contribute significantly to mixing and energy transporting in the ocean. We observed highly nonlinear ISWs over the continental shelf of the northwestern South China Sea (19A degrees 35'N, 112A degrees E) in May 2005 during the Wenchang Internal Wave Experiment using in-situ time series data from an array of temperature and salinity sensors, and an acoustic Doppler current profiler (ADCP). We summarized the characteristics of the ISWs and compared them with those of existing internal wave theories. Particular attention has been paid to characterizing solitons in terms of the relationship between shape and amplitude-width. Comparison between theoretical prediction and observation results shows that the high nonlinearity of these waves is better represented by the second-order extended Korteweg-de Vries (KdV) theory than the first-order KdV model. These results indicate that the northwestern South China Sea (SCS) is rich in highly nonlinear ISWs that are an indispensable part of the energy budget of the internal waves in the northern South China Sea.

Effect of surface tension on the mode selection of vertically excited surface waves in a circular cylindrical vessel

Singular perturbation theory of two-time-scale expansions was developed in inviscid fluids to investigate patternforming, structure of the single surface standing wave, and its evolution with time in a circular cylindrical vessel subject to a vertical oscillation. A nonlinear slowly varying complex amplitude equation, which involves a cubic nonlinear term, an external excitation and the influence of surface tension, was derived from the potential flow equation. Surface tension was introduced by the boundary condition of the free surface in an ideal and incompressible fluid. The results show that when forced frequency is low, the effect of surface tension on the mode selection of surface waves is not important. However, when the forced frequency is high, the surface tension cannot be neglected. This manifests that the function of surface tension is to cause the free surface to return to its equilibrium configuration. In addition, the effect of surface tension seems to make the theoretical results much closer to experimental results.

Surface waves in a vertically excited circular cylindrical container

The nonlinear free surface amplitude equation, which has been derived from the inviscid fluid by solving the potential equation of water waves with a singular perturbation theory in a vertically oscillating rigid circular cylinder, is investigated successively in the fourth-order Runge-Kutta approach with an equivalent time-step. Computational results include the evolution of the amplitude with time, the characteristics of phase plane determined by the real and imaginary parts of the amplitude, the single-mode selection rules of the surface waves in different forced frequencies, contours of free surface displacement and corresponding three-dimensional evolution of surface waves, etc. In addition, the comparison of the surface wave modes is made between theoretical calculations and experimental measurements, and the results are reasonable although there are some differences in the forced frequency.