Comparison and Modification: TVD Schemes for Scalar Transport on An Unstructured Grid

In most TVD schemes, the r-factors were proposed according to the cell-centered (CC) finite volume method (FVM) framework for the numerical approximation to the convective term. However, it is questionable whether those r-factors would be appropriate and effective for the vertex-centered (VC) FVM. In the paper, we collected five kinds of r-factor formulae and found out that only three of those, respectively by Bruner (1996), Darwish and Moukalled (2003) and Cassuli and Zanolli (2005) can be formally extended to a context of the VC FVM. Numerical tests indicate that the TVD schemes and r-factors, after being extended and introduced to a context of the VC FVM, maintained their similar characteristics as in a context of the CC FVM. However, when the gradient-based r-factors and the SUPERBEE scheme were applied simultaneously, non-physical oscillations near the sharp step would appear. In the transient case, the oscillations were weaker in a context of the VC FVM than those in a context of the CC FVM, while the effect was reversed in the steady case. To eliminate disadvantages in the gradient-based r-factor formula, a new modification method by limiting values on the virtual node, namely in the paper, was validated by the tests to effectively dissipate spurious oscillations.     

An Explicit High Resolution Scheme for Nonlinear Shallow Water Equations

The present study develops a numerical model of the two-dimensional fully nonlinear shallow water equations （NSWE） for the wave run-up on a beach. The finite volume method （FVM） is used to solve the equations, and a second-order explicit scheme is developed to improve the computation efficiency. The numerical fluxes are obtained by the two dimensional Roe＇ s flux function to overcome the errors caused by the use of one dimensional fluxes in dimension splitting methods. The high-resolution Godunov-type TVD upwind scheme is employed and a second-order accuracy is achieved based on monotonic upstream schemes for conservation laws （MUSCL） variable extrapolation; a nonlinear limiter is applied to prevent unwanted spurious oscillation. A simple but efficient technique is adopted to deal with the moving shoreline boundary. The verification of the solution technique is carried out by comparing the model output with documented results and it shows that the solution technique is robust.     

A Compact Explicit Difference Scheme of High Accuracy for Extended Boussinesq Equations

Presented here is a compact explicit difference scheme of high accuracy for solving the extended Boussinesq equations.For time discretization,a three-stage explicit Runge-Kutta method with TVD property is used at predicting stage,a cubic spline function is adopted at correcting stage,which made the time discretization accuracy up to fourth order;For spatial discretization,a three-point explicit compact difference scheme with arbitrary order accuracy is employed.The extended Boussinesq equations derived by Beji and Nadaoka are solved by the proposed scheme.The numerical results agree well with the experimental data.At the same time,the comparisons of the two numerical results between the present scheme and low accuracy difference method are made,which further show the necessity of using high accuracy scheme to solve the extended Boussinesq equations.As a valid sample,the wave propagation on the rectangular step is formulated by the present scheme,the modelled results are in better agreement with the experimental data than those of Kittitanasuan.     

Normal Mode as a Cause of Large-Scale Variations in the Troposphere and Stratosphere

Izvestiya, Atmospheric and Oceanic Physics - In a study of zonally averaged meteorological fields in the middle atmosphere, synchronous oscillations (in phase or antiphase) occupying extended...     

不同TVD 格式对内孤立波数值模拟结果影响研究

为了研究不同TVD格式对内孤立波模拟结果的影响,利用改进后的SUNTANS三维非静压海洋模型,通过理想算例和南海东北部海域内孤立波的模拟,比较分析了求解温盐方程TVD格式的4种经典通量限制函数(superbee, minmod, van Leer和MUSCL)对计算结果的影响。综合理想算例和南海东北部海域模拟结果,在所讨论的4种通量限制函数中,利用MUSCL限制函数所得模拟结果最优,建议在南海东北部海域内孤立波模拟中采用MUSCL限制函数。     

Efficient higher-order finite-difference schemes for parabolic models

Implicit finite-difference schemes for use in parabolic equation models are developed. Like the familiar Crank-Nicolson scheme, which has hitherto been used almost exclusively for the solution of these equations, these schemes are unconditionally stable and use a computational molecule of only six points on two “time” levels. However, they are accurate to a higher order than the Crank-Nicolson scheme, thus allowing the solution grid to be coarser and the solution time to be (approximately) halved. Examples of computations on constant depth are shown, in which significant reductions in time and grid-point density are achieved, for two different parabolic models. The schemes are then extended to refraction and diffraction, and are shown to have a similar effect in this more general case too. It is recommended that finite-difference schemes based on these higher-order (or Hermitian) methods replace the more commonly used Crank-Nicolson scheme in all physical domain parabolic equation models, but especially in minimax (wide-angle) equation models.     

Radiative transfer modeling and analysis of spatially variant and coherent illumination for undersea object detection

Increasing the optical range of target detection and recognition continues to be an area of great interest in the ocean environment. Light attenuation limits radiative and information transfer for image formation in water. In this paper, the authors briefly review current methods of imaging and then describe a variation of the spatial interferometric technique that relies upon projected spatial gratings with subsequent detection against a coherent return signal for the purpose of noise reduction and image enhancement. A model is developed that simulates the projected structured illumination through turbid water to a target and its return to a detector. The model shows an unstructured backscatter superimposed upon a structured return signal. The model has been extended to predict what a camera would actually see, so that various noise-reduction schemes can be modeled. Finally, some water-tank tests are presented, validating original hypothesis and model predictions. The method is advantageous in not requiring temporal synchronization between reference and signal beams and may use a continuous illumination source. Spatial coherency of the beam allows for the detection of the direct return, while scattered light appears as a noncoherent noise term.     

无结构网格二维河口海岸水动力数值模式的建立及其应用

为完全拟合河口近海复杂岸线和工程结构以及有效局部加密,设计并建立了一个无结构三角形网格二维河口海岸水动力数值模式。空间离散主要基于有限体积法以保证守恒性,时间积分采用预估修正法以提高精度。水位在三角形网格中心通过连续方程求解;水平x方向和y方向的流速U和V均在网格边中点上通过动量方程求解。流速平流项的求解中采用了TVD格式。TVD流速平流通量为一个一阶迎风格式通量和一个二阶格式通量的组合,一阶格式通量和二阶格式通量根据流速的局部分布情况得出配比,最终组合得到TVD通量。TVD格式具有低耗散和无频散的优点,提高了模式的稳定性。应用实测资料验证建立的模式,结果显示水位、流速和流向的计算值与实测值均符合良好。     

Statistical analysis of series of dendrochronological indices

We analyze series of tree-ring indices to understand whether they can be indicators of long-term climate changes. We calculated the synchronous cross-correlation coefficients (CCs) for all possible pairs composed of nine series, of which only one CC (for a pair of the longest series) proved to be sufficiently high (0.72). We revealed a time interval for which this CC is equal to 0.88. We failed to locate sufficiently larger time intervals with a higher CC for other pairs of series. Similar calculations for a series using a low-frequency filter (with allowance for periods of 100 or more years) led to a noticeable decrease in almost all CCs. This confirms the statements circulating in the literature that tree-ring chronologies include climate signals with characteristic times of several decades and that the secular and supersecular oscillations are significantly distorted as a result of standardization and cross-linking. The results presented in this paper show that the complex of series being considered cannot be used as a climate-change indicator without an additional analysis. Here, we describe a possible procedure for such analysis. As a result, we found that only two time series at some time subinterval can provide useful information about climate changes.     

Second-Order Solutions for Random Interfacial Waves Under Steady Uniform Currents

1.Introduction Owingtothesimilaritiesbetweenthesurfacewaveandtheinterfacialwave,itisnaturaltoapply themethodsdevelopedforsurfacewavestothestudyofinterfacialwavesasreviewedbyUmeyama (1998;2000).Recently,Song(2004)derivedsecond ordersolutionsforrandominterfacialwavesat aconstantdepthinatwo layerfluidsystemwitharigidlidusinganexpansiontechnique,analogousto thatusedbyLonguet Higgins(1963)andSharmaandDean(1979),tostudyrandomsurfacewaves. Inthispaper,Song’sresultsareextendedtoamoregeneralcaseoft…     

Convective adjustment in box models

A convective adjustment (CA) algorithm is thought to be responsible for grid-scale oceanic-state sustained oscillations seen in oceanic general circulation models (OGCM), an effect that is most evident in simulations with coarse spatio-temporal scales. The CA algorithm is thought to inadvertently create a salt oscillator. Several studies have confirmed that a flip-flop type salt oscillator, which is reminiscent in some respects of simple CA schemes, can develop sustained oscillations. Subsequently, several researchers were able to show how coupled salt oscillators, reacting in a particular temporal sequence, are capable of producing large-scale oscillations not unlike those found in the OGCM simulations. However, the proxy models used to study how CA can create these oscillations in large-scale simulations were never directly related to OGCM results.Here we couple hydrodynamics to the CA and look at zonally-driven flows in the low-frequency, large-scale limit. Adding flow is a step in the direction of developing an analytically tractable model with which to understand the basics of OGCMs. We analytically determine whether, and under what circumstances, the CA scheme is responsible for sustained oscillations.We carry out this program for four basic box-model configurations, each inspired by the general shape of the eigenfunctions and constraints of the large-scale zonally-averaged forced flow over a hemisphere. Furthermore, in order to make our results relevant to the Meridional Ocean Circulation, we also investigate the effect of replacing the usual assumption of a linear relation between thermohaline flow rate and horizontal density gradient with a nonlinear hydraulic relationship.We find that a salt oscillator does not occur in the most common box-model configurations. In one of our models, however, we find wide parameter ranges in which all steady states calculated for the model fail to satisfy the CA scheme, the situation which is expected to result in CA-induced oscillations. The model in question corresponds to a hemispheric shallow thermohaline flow over a deep reservoir. However, we find that oscillations occur in these parameter ranges only if the density threshold for convection is negative, i.e., if the CA scheme turns on convection between vertically adjacent boxes when the density stratification between them is still slightly stable. In this situation, the amplitude and period of the oscillations depend strongly on the size of the density threshold, both vanishing as the threshold is taken toward zero. We also show that the same is true in the Welander flip-flop model of a single salt oscillator. For positive values of the threshold, that is, when the CA scheme is allowed to ignore small unstable stratification changes, oscillations do not occur in the limit of integration time step going to zero, but can still be seen when the time step is finite, even if small. Moreover, the system evolves toward a new steady state, one in which the stratification in one box is exactly the threshold value itself. We show how to calculate these new steady states, and explain why they give way to oscillations when the density threshold is negative.     

The application of a Godunov-type shock capturing scheme for the simulation of waves from deep water up to the swash zone

A two-dimensional vertical (2DV) non-hydrostatic boundary fitted model based on a Godunov-type shock-capturing scheme is introduced and applied to the simulation of waves from deep water up to the swash zone. The effects of shoaling, breaking, surf zone dissipation and swash motions are considered. The application of a Godunov-type shock-capturing algorithm together with an implicit solver on a standard staggered grid is proposed as a new approach in the 2DV simulation of large gradient problems such as wave breaking and hydraulic jumps. The complete form of conservative Reynolds averaged Navier–Stokes (RANS) equations are solved using an implicit finite volume method with a pressure correction technique. The horizontal advection of the horizontal velocity is solved by an explicit predictor–corrector method. Fluxes are predicted by an exact Riemann solver and corrected by a downwind scheme. A simple total variation diminishing (TVD) method with a monotonic upstream-centered scheme for conservation laws (MUSCL) limiter function is employed to eliminate undesirable oscillations across discontinuities. Validation of the model is carried out by comparing the results of the simulations with several experimental test cases of wave breaking and run-up and the analytical solution to linear short waves in deep water. Promising performance of the model has been observed.     

Automatic indexing of underwater survey video: algorithm and benchmarking method

It is often the case that only a few sparse sequences of long videos from scientific underwater surveys actually contain important information for the expert. Locating such sequences is time consuming and tedious. A system that automatically detects those critical parts, online or during post-mission tape analysis, would alleviate the expert workload and improve data exploitation. In this paper, a methodology for evaluating the performance of such a system on real data is presented. Interesting sequences are started by changes of visual context. An algorithm to detect significant context changes in benthic videos in real time has been presented by Lebart et al. in 2000. It is used as an illustration for this methodology - its performance is studied and benchmarked on real underwater data, ground truthed by an expert biologist. Various issues relating to the complexity of the problems of automatically analyzing underwater video are also discussed.     

Cause and countermeasure of long-period oscillations of moored ships and the quantification of surge and heave amplitudes

Long-period oscillations of moored ships whose periods are about 1 or 2 min cause many troubles in many ports and harbours. It is necessary to investigate these phenomena and verify their causes and countermeasures in each case because they are strongly dependent on the environment of each port and harbour. From this point of view, long-period oscillations of moored ships in the Port of Shibushi in Japan were investigated by means of wave observations, the image processing of moored ship motions using the video camera and motion-capture software and numerical simulations. From observation results, the relationship between offshore long-period waves and long-period oscillations of moored ships was recognized and surge and heave amplitudes were quantified by using wave data in order to forecast moored ship motions. Furthermore, from observation and numerical results, it was revealed that long-period waves with the peak period of 120 s from the offshore typhoon kept or exaggerated the local harbour oscillation of 60–70 s and it caused long-period oscillations of moored ships. Numerical results in case of reducing the reflection coefficient of the target berth implied that it ceased the local harbour oscillation and it would give an effective countermeasure to reduce long-period oscillations of moored ships in the Port of Shibushi.     

Numerical Simulation of Self-Excited and Forced Vibration of Circular Cylinders in Current

Numerical simulations of a low-mass-damping circular cylinder which can oscillate freely at transverse and streamwise directions are presented in this work.The Navier-Stokes equations are solved with finite volume method,and large eddy simulation of vortex is also performed in the calculation.In order to implement dynamic mesh,overlapping grids are generated to lessen the computation for mesh field itself.Self-excited vibrations are firstly calculated to obtain the average amplitudes and frequencies of the target circular cylinder in the current flow situation,and then forced oscillations are implemented with parameters obtained in vortex-induced vibrations previously.With slight amplitude modulation,time series of displacements in vortex-induced vibrations are essentially harmonic.Regarding the fluid force,which are larger in forced oscillations than those in corresponding self-excited cases because the fluid subtracts energy from the forced cylinders.The phase angles between forces and displacements are 0° and 180° for self-excited case and forced case respectively.In vortex-induced vibrations,the interactions between fluid and structure produce some weakly energetic vortices which induce the modulations of amplitude and frequency.     

Numerical Simulation of Self-Excited and Forced Vibrating Circular Cylinders in Current

Numerical simulations of a low-mass-damping circular cylinder which can oscillate freely at transverse and streamwise directions are presented in this work. The Navier-Stokes equations are solved with finite volume method, and large eddy simulation of vortex is also performed in the calculation. In order to implement dynamic mesh, overlapping grids are generated to lessen the computation for mesh field itself. Self-excited vibrations are firstly calculated to obtain the average amplitudes and frequencies of the target circular cylinder in the current flow situation, and then forced oscillations are implemented with parameters obtained in vortex-induced vibrations previously. With slight amplitude modulation, time series of displacements in vortex-induced vibrations are essentially harmonic. Regarding the fluid force, which are larger in forced oscillations than those in corresponding self-excited cases because the fluid subtracts energy from the forced cylinders. The phase angles between forces and displacements are 0o and 180o for self-excited case and forced case respectively. In vortex-induced vibrations, the interactions between fluid and structure produce some weakly energetic vortices which induce the modulations of amplitude and frequency.     

异育银鲫准回交世代经济性状优势的遗传基础研究

用RAPD技术分析了方正银鲫、异育银鲫、异育银鲫准回交世代和兴国红鲤4个群体的遗传相似性。方正银鲫、异育银鲫和异育银鲫准回交世代多态位点数接近。且显著低于兴国红鲤;异育银鲫准回交世代存在特异带或特有带;方正银鲫与异育银鲫间的遗传相似性高于它与异育银鲫准回交世代间的遗传相似性,而兴国红鲤与异育银鲫间的相似性低于其与异育银鲫准回交世代间的遗传相似性。结果表明：三倍体群体遗传多样性低于正常二倍体;异育银鲫准回交世代存在遗传重组;随着兴国红鲤精子刺激次数增加。后代遗传背景与方正银鲫差异越来越大,而与兴国红鲤越来越相似。精子刺激三倍体卵发育,也能将一些遗传物质整合进卵细胞核,使后代经济性状呈现一定的优势。     

A prediction technique for the transverse vortex-induced oscillations of tensioned risers

In this paper, a model for the prediction of the vortex-induced transverse oscillations of tensioned risers is proposed.The model, based on the peak amplitudes of the oscillations recorded for cylindrical bodies in a uniform current, has been extended to multimodal oscillations in a sheared flow condition through a modal priority technique. The difference between this model and other similar ones lies in the following features:

the hydrodynamic damping, which has been introduced in a more general form without ‘a priori’ approximations;

the modal priority, which is based on the local amplitudes of the active vibration modes, rather than on their frequencies.

Some comparisons with experimental and numerical results available in the literature are presented and discussed.     

Modelling fish trade liberalisation: Does fish trade liberalisation result in welfare gains or losses?

Recent research has on the basis of general equilibrium models warned that trade liberalisation in non-optimally managed renewable resources might cause over-exploitation and reduced steady-state welfare. Welfare effects of trade liberalisation in open access and optimal management are found case specific and dependent on factors such as the country's status as importer or exporter, the state of the fish stocks and the size of the countries on the world market. The present paper develops an alternative partial equilibrium framework capable of identifying welfare effects of fish trade liberalisation also in the presence of complex but realistic management schemes, such as regulated open access and regulated restricted access. The case dependency known from the general equilibrium analysis is confirmed, but the introduction of realistic fisheries management schemes in the partial equilibrium model extended this case dependency. The welfare effect of trade liberalisation in an exporter country is negative under open access, but was found positive under regulated restricted access in the present paper.     

A better,more discriminating test problem for ocean tracer transport

Many test problems have been posed in which the behavior of schemes may be examined within the simplified context of passive tracer transport. Most often these tests involve transport of a tracer containing sharp or discontinuous gradients within a smooth advecting field. Hecht et al. [J. Geophys. Res. 100 (1995) 20763–20778] presented a test problem, using the simple analytical gyre of Stommel [Trans. Am. Geophys. Union 29 (1948) 202–206], in which an initially smooth tracer field is transported through a highly sheared circulation, representing another scenario of relevance to ocean modeling. We present a modification of that test in which the underlying grid is rotated relative to the gyre, such that the western boundary falls at 45° to the principal grid axes. The advection schemes used by most ocean modelers to date are examined within this test problem. We find, in two cases, that methods which performed acceptably under the original formulation of the problem, and which have seen extensive application to geophysical fluid dynamics, fail to perform acceptably within the reformulated problem. The important qualities of this more discriminating test include the skewing of the fast boundary current flow relative to the underlying grid axes and the high shear in that boundary current.