Numerical Solution of Lock-Release Gravity Current with Viscous Self-Similar Regime

Lock-release gravity currents with a viscous self-similar regime are simulated by use of the renormalization group(RNG) k - ε model for Reynolds-stress closure. Besides the turbulent regime with initially a slumping phase of a conslant current front speed and later an inviseid self-similar phase of front speed decreasing as t^-1/3(where t is the time measured from release), the viseous self-similar regime is satisfactorily reproduced with front speed decreasing as t^-4/5,consistent with well known experimental observations.     

The impact of the M2 internal tide on data-assimilative model estimates of the surface tide

We use a synthetic data experiment to assess the accuracy of ocean tides estimated from satellite altimetry data, with emphasis on the impact of the phase-locked internal tide, which has a surface expression of several centimeters near its sites of genesis. Previous tidal estimates have regarded this signal as a random measurement error; however, it is deterministic and not scale-separated from the barotropic (surface) tide around complex bathymetric features. The synthetic data experiments show that the internal tide has a negligible impact on the barotropic tidal fields inferred under these circumstances, and the barotropic dissipation (a quadratic functional of the tidal fields) is in good agreement with the energetics of the three-dimensional regional primitive equations model which is the source of the synthetic data.     

“数字海底”系统建模中的关键技术研究

基于“数字海底”系统的理论基础研究和在海底科学应用中的具体实践,提出了“数字海底”的涵义、系统组成和体系结构,对空间数据管理技术、基于XML的WebGIS空间数据表达技术、复杂海底三维建模技术、多源信息的一体化集成显示技术、虚拟现实技术等“数字海底”系统建模中的关键技术进行了探讨,为从事海底资源研究和开发的人员提供各个层次的参考或帮助.     

水下滑翔器的运动建模与分析

介绍了水下滑翔器的工作机理，对其沉浮阶段的滑翔过程进行了动力学分析，推导了滑翔器在垂直剖面上的动力学方程。论文深入分析了水下滑翔器稳态时的运动规律，以水下滑翔器试验模型为例，推导了其稳态运动参数，通过线性化与适当的简化，得到模型在垂直剖面上的运动状态方程，讨论了系统的可控性与可观测性，为水下滑翔器系统的开发设计和控制提供了理论依据，具有重要的指导意义。     

Inverse reconstructions of ecosystem flows in investigating regime shifts: impact of the choice of objective function

Inverse analysis is increasingly used in ecosystem modelling to objectively reconstruct a large number of unknown flows or interactions from a small number of observations. This type of analysis may be useful in relating observed regime shifts in ecosystem structure to underlying processes. Inversions of ecosystem flow networks currently use a constrained least-squares solution which at the same time minimizes the squared norm (the sum of squares) of the reconstructed flows. This minimum norm (MN) inversion is thought to be a parsimonious solution to the ecosystem flow inverse problem, but it may well not reflect how ecosystems are organised. It has been proposed instead that ecosystems evolve to maximize energy/mass flows or that they maximize the information content of the network weighted by ecosystem flows (ascendancy). We used simulated inverse experiments, where inverse analyses are applied to simulations of flow networks, to explore objective functions different than the MN generally used. We could not compute inverse solutions that maximize ascendancy because the objective function is unbounded. We could calculate inversions that maximize flows; however, these generally overestimated the simulated flows, even though the simulations were designed to maximize flows. It appears that the ecosystem flow inverse problem is too under-determined (too few data relative to the number of unknowns) to allow the use of these maximizing goal functions. We introduce a new minimization that simultaneously minimizes the squared flows and the squared differences between flows. This smoothing minimization makes the inverse flows as even as possible and it helps with some technical issues with MN inversions. The simulated inverse experiments indicated that this smoothed norm (SM) is the most robust in comparative analyses of contrasting ecosystem states, such as those that can be associated with regime shifts. Like the MN inversion, the SM inversion has no ecological basis. However, it is a conservative norm that is less likely to produce false differences between the dynamics of regimes.     

Uncertainty in predictions of oil spill trajectories in open sea

A method is introduced to calculate and to account for the uncertainties in the predictions of oil spill trajectories using a classic oil spill model. The method considers the output of the oil spill model as a function of random variables, which are the input parameters, and calculates the standard deviation of the output results which gives a measure of the uncertainty of the model given the uncertainties of the input parameters.Instead of a single trajectory that is calculated by the oil spill model using the mean values of the parameters, a band of trajectories can be defined when various simulations are done taking into account the uncertainties of the input parameters. This band of trajectories defines envelopes of the trajectories that are likely to be followed by the spill given the uncertainties of the input.The method is applied to an oil spill that occurred in open sea near Madeira Islands, in the Atlantic Ocean, in December 1989. The simulations allow the understanding of how a change in the wind direction drove the spill towards the Islands.The envelope of likely trajectories that is obtained with the uncertainty modelling shows a band of trajectories that is in better agreement with the observations than the single trajectory simulated by the oil spill model, based on mean parameters.     

Numerical study of the meridional overturning circulation with “mixing hotspots” in the Pacific Ocean

Using an idealized ocean general circulation model, we examine the effect of “mixing hotspots” (localized regions of intense diapycnal mixing) predicted based on internal wave-wave interaction theory (Hibiya et al., 2006) on the meridional overturning circulation of the Pacific Ocean. Although the assumed diapycnal diffusivity in the mixing hotspots is a little larger than the predicted value, the upwelling in the mixing hotspots is not sufficient to balance the deep-water production; out of 17 Sv of the downwelled water along the southern boundary, only 9.2 Sv is found to upwell in the mixing hotspots. The imbalance as much as 7.8 Sv is compensated by entrainment into the surface mixed layer in the vicinity of the downwelling region. As a result, the northward transport of the deep water crossing the equator is limited to 5.5 Sv, much less than estimated from previous current meter moorings and hydrographic surveys. One plausible explanation for this is that the magnitude of the meridional overturning circulation of the Pacific Ocean has been overestimated by these observations. We raise doubts about the validity of the previous ocean general circulation models where diapycnal diffusivity is assigned ad hoc to attain the current magnitude suggested from current meter moorings and hydrographic surveys.     

The Structure of A Turbulent Jet in A Crossflow-Effect of Jet-Crossflow Velocity

A comprehensive numerical study on the three-dimensional structure of a turbulent jet in crossflow is performed. The jet-to-crossflow velocity ratio (R) varies in the range of 2 - 16; both vertical jets and inclined jets without excess streamwise momentum are considered. The numerical results of the Standard two-equation k-ε model show that the turbulent structure can be broadly categorised according to the jet-to-crossflow velocity ratio. For strong to moderate jet discharges, i.e. R> 4, the jet is characterized by a longitudinal transition through a bent-over phase during which the jet becomes almost parallel with the main freestream, to a sectional vortex-pair flow with double concentration maxima; the computed flow details and scalar mixing characteristics can be described by self-similar relations beyond a dimensionless distance of around 20-60. The similarity coefficients are only weakly dependent on R. The cross-section scalar field is kidney-shaped and bifurcated, vvith distinct double concentr