一种河道紊流三维数值模拟离散格式的应用探讨

通过引用气体动力学中的Osher格式计算控制体交界面通量,建立了一种基于有限体积法的河道紊流三维数学模型。该格式是对准确的非线性黎曼问题的近似解法,具有很好的守恒性、逆风性、单调性和高分辨率捕捉间断的能力,解决了河床变化剧烈处水流计算易发散的难题,在紊流模型的选取上,考虑到紊动能量和紊动尺度在水平方向和垂直方向上的差异,平面上采用大涡模型,垂向上采用k-ε模型。通过桥墩和丁坝局部水流实测资料校核模型模拟局部复杂水流结构的能力,模型计算得到的马蹄涡、回流结构及流速分布与实测基本一致。     

A numerical model of the long term flow along the Malin-Hebrides shelf

A three dimensional hydrodynamic model of the Malin-Hebrides shelf region is used to investigate the spatial variability of the wind and tidally induced residual flow in the region and the influence of flow from the Irish Sea and along the shelf edge. By this means it is possible to understand the spatial variability in the long term observed flow fields in the region and the range of driving forces producing this flow. The model uses a sigma coordinate grid in the vertical with a finer grid in the near surface and near bed shear layers. The vertical diffusion of momentum in the model is parameterised using an eddy viscosity coefficient which is derived from turbulence energy closure models. Two different turbulence models are used to compute the eddy viscosity, namely a two-equation (itq²−q²ℓ) model which has prognostic equations for both turbulence energy and mixing length and a simpler model in which the mixing length is a specified algebraic function of the water depth.The wind induced response to spatially and temporally constant orthogonal wind stresses, namely westerly and southerly winds of 1 N m⁻², are derived from the model. By using orthogonal winds and assuming linearity, then to first order the response to any wind direction can be derived. Computed flows show a uniform wind driven surface layer of magnitude about 3% of the wind speed and direction 15 ° to the right of the wind, in deep water. Currents at depth particularly in the shelf edge and near coastal region show significant spatial variability which is related to variations in bottom topography and the coastline.Calculations show that tidal residual flows are only significant in the near coastal regions where the tidal current is strong and exhibits spatial variability. Flow into the region from the Irish Sea through the North Channel although having its greatest influence in the near coastal region, does affect currents near the shelf edge region. Again the spatial variability of the flow is influenced by topographic effects.A detailed examination of wind induced current profiles together with turbulence, mixing length and viscosity, at a number of locations in the model from deep ocean to shallow near coastal, shows that both turbulence models yield comparable results, with the mixing length in the two equation model showing a similar dependence to that specified in the simpler turbulence model.Calculations clearly show that flow along the shelf edge area to the west of Ireland and from the Irish Sea entering the region, together with local wind forcing can have a major effect upon currents along the Malin-Hebrides shelf. The flow fields show significant spatial variability in the region, comparable to those deduced from long term tracer measurements. The spatial variability found in the calculations suggests that a very intense measurement programme together with inflow measurements into the area is required to understand the circulation in the region, and provide data sets suitable for a rigorous model validation.     

Spatial distribution of cadmium over a cyclonic eddy in the southern East China Sea

Eddy systems are a unique ecosystem, usually having high biological masses and primary production in the sea. In this study, both particulate and dissolved phases of cadmium in the water column of 15 stations over a cyclonic eddy in the southern East China Sea were determined to obtain their spatial distributions. This allows us to confirm that cyclonic eddy systems play the role of a Cd pump in the sea.Results showed that particulate Cd (PCd) and dissolved Cd (DCd) concentrations in water varied greatly, ranging over two orders and one order of magnitude, respectively. Large spatial variability was found not only for PCd but also for DCd in the upper water, apparently due to the effects of the cyclonic eddy system on the Cd distributions over the study area. DCd accounts for about 99% of the total Cd. For the surface water, DCd concentration at the eddy center was about five times the average of the water surrounding the eddy center. The depth distributions of DCd exhibited a typical surface depletion and a subsequent increase with depth; however, the PCd distribution showed the opposite, i.e. a surface maximum and a subsequent decline with depth. In general, the DCd maxima were found at depths of 600–1000 m, agreeing well with the literature. It reflected the internal biogeochemical cycling of Cd in the water column, which was driven by the utilization of Cd by plankton in the euphotic zone and by the regeneration of Cd at depth. In addition, a remarkably high DCd concentration existed in near-bottom water around the cyclonic eddy center.The horizontal distributions of both PCd and DCd in the upper water shared a common feature with elevated concentrations centering around the eddy center and a decline in concentration with distance from the eddy center. This shows that the cyclonic eddy could bring up the Cd-rich deep water to the surface water around the eddy center and could then expand toward eddy's vicinity via advection–diffusion. It is thus likely that it produces ample biological masses over the eddy system. Hence, this work can demonstrate that the ascending nutrient-rich water driven by the cyclonic eddies can serve as an important source not only for many nutrients but also for Cd in the sea.     

船舶涂装车间VOCs废气排除效果数值模拟

针对大型船舶涂装车间，建立用于两种不同通风方案和挥发性有机物（Volatile Organic Compounds，VOCs）废气扩散的数学模型。运用数值模拟方法，研究不同通风方案的气流速度场和气流压力场的分布特点及其对VOCs废气排除效果的影响规律。结果表明：侧进侧排（Side Intake and Side Discharge，SISD）式的原有通风方案导致车间VOCs废气不断向车间顶部累积，通风效率较低；而侧进底排顶部补风式的优化通风方案可有效提升车间湍流区域面积和风速，对VOCs废气排除效果作用显著。研究结果可为船舶涂装车间通风系统的优化设计提供参考依据。     

Variability of the shelf water off the northeast Spanish coast

During the seasonal stratification of 1983 and 1984, intensive sampling of the shelf region of the northeast Spanish coast was carried out. This unique dataset is used to define the most salient features of the shelf water dynamics. We show that region has a strong spatial and temporal variability associated with the southward spreading of continental waters from the Rhone river, in the northern gulf of Lyons. Southward spreading begins in early spring and was particularly obvious during June 1983 when strong horizontal density gradients were observed. Spreading almost completely disappears in August when observed upper layer salinities are uniform. We also show that significant shelf/slope water exchange takes place associated with the abrupt topographic canyons found in the region. In particular we observed in June 1983 a small positive geostrophic vorticity region at the Palamós canyon while a larger anticyclonic eddy was observed over the shelf south of this canyon. The biological implications of these two features have been demonstrated. We conclude indicating that in this region where no dominant winds exist, the shelf water dynamics is dominated by the southward spreading of low salinity continental waters and by the intrusion of high salinity open ocean water through submarine canyons.     

Modelling physical processes of haline stratification in ROFIs with application to the Rhine outflow region

A physical and numerical study is made of the processes governing the stratification and circulation in ROFIs (Regions of Freshwater Influence) where there is an important impact of wind and tides. Observations in the Rhine ROFI showed that the salinity field consists of a mean and a tidally oscillating part. The physical processes are first analysed using the analytical solutions from a one-dimensional two-layer model. A justification is given for the neglection of non-linear advective terms in the equations of momentum and salinity. The dimensionless forms of the solutions can be expressed in terms of a series of dimensionless numbers. It is shown in particular that stratification and cross-shore circulation largely depend on the balance between rotation and turbulent diffusion, which depends in turn on parameters such as the Ekman number, the bottom friction coefficient, the eddy viscosity ratio and the depth of the layer interface. Surface winds either enhance or destroy stratification depending on the wind angle. The response to wind forcing is discussed using classical Ekman theory. To verify the analytical theory numerical tests are performed with a point model including an advanced turbulence closure scheme. Differences arise due to the non-linear interaction between turbulence on the one hand and current shear and stratification on the other hand. It is shown in particular that the amplitude of the tidal forcing and the off-shore horizontal salinity gradient strongly affect the semi-diurnal and semi-monthly variation of stratification. The effect of the wind is found to be in good agreement with the analysis of the two-layer model. Finally, the numerical model is compared with existing observational data in the Rhine ROFI for October 1990.     

7500t起重船压载舱的优化研究

针对7500吨起重船的作业进行了简单的受力分析,并且对原起重船压载舱的特点进行了详细的总结。结合实际经验和相关理论,对原起重船的压载舱提出了四套改进方案,选取四种典型的工况,采用传统的多方案比较方法,从总压载量和压载水的调拨量两个主要方面进行综合评价,得出最优的压载方案。最后,总结了起重船压载舱优化的一些规律。     

A 3D, finite element model for baroclinic circulation on the Vancouver Island continental shelf

This paper describes the development and application of a 3-dimensional model of the barotropic and baroclinic circulation on the continental shelf west of Vancouver Island, Canada. A previous study with a 2D barotropic model and field data revealed that several tidal constituents have a significant baroclinic component (the K1 in particular). Thus we embarked on another study with a 3D model to study the baroclinic effects on the residual and several selected tidal constituents.The 3D model uses a harmonic expansion in time and a finite element discretization in space. All nonlinear terms are retained, including quadratic bottom stress, advection and wave transport (continuity nonlinearity). The equations are solved as a global and a local problem, where the global problem is the solution of the wave equation formulation of the shallow water equations, and the local problem is the solution of the momentum equation for the vertical velocity profile. These equations are coupled to the advection-diffusion equation for density so that density gradient forcing is included in the momentum equations. However, the study presented here describes diagnostic calculations for the baroclinic residual circulation only.The model is sufficiently efficient that it encourages sensitivity testing with a large number of model runs. In this sense, the model is akin to an extension of analytical solutions to the domain of irregular geometry and bottom topography where this parameter space can be explored in some detail.In particular, the consequences of the sigma coordinate system used by the model are explored. Test cases using an idealized representation of the continental shelf, shelf break and shelf slope, lead to an estimation of the velocity errors caused by interpolation errors inherent in the sigma coordinate system. On the basis of these estimates, the computational grid used in the 2D model is found to have inadequate resolution. Thus a new grid is generated with increased accuracy in the region of the shelf break. However, even with increased resolution, spurious baroclinic circulation seaward of the shelf break and in the vicinity of Juan de Fuca canyon remained a significant problem when the pressure gradient terms were evaluated using the σ coordinate system and using a realistic density profile.With the new grid, diagnostic calculations of the barotropic and baroclinic residual circulation are performed using forcing from the observed σ_t (density) field and from the gradient of this field.     

通州湾港区围填海工程起步方案比较

围填海工程在促进区域经济发展的同时会对周边海域滩槽形态、生态环境产生一定的影响。以通州湾港区起步工程为例，通过建立潮流数学模型，模拟南黄海辐射沙脊群海域的潮流形态。在此基础上，比较了两套不同围填海开发方案对工程区域滩面的冲淤影响。从整体上看，两套方案对工程周边海域的海底冲淤产生明显影响，但随着时间的推移趋于平衡。基于目前的模拟结果，先期实施离岸式的二港池更有利于维持小庙洪水道的畅通以及滩涂资源的保护。     

Evolution of the open-sea eddy ALGERS'98 in the Algerian Basin with Lagrangian trajectories and remote sensing observations

The westward evolution of an open-sea anticyclonic eddy along the western Algerian Basin is shown, for the first time, by means of 15 buoy trajectories and remote sensing observations. For 3 months, the buoy trajectories described several anticyclonic loops in periods of 4–21 days. The eddy's movement, translation, and rotation were separated with a kinematic model, resulting in a mean translation speed of 2 km/day, which fits the self-propulsion speed predicted on theoretical models for isolated eddies on a beta plane. Fluctuations in translation speed were associated with advection of the mean flow and topographic interactions. Both mechanisms changed the eddy's horizontal shape from circular to elliptical, inducing fluctuations in its swirl velocity and solid-body rotation. The initial stage of the eddy is an isolated asymmetric dipole, comprised by a small cyclone and a large anticyclone, the latter generated from a frontal instability, which under the Coriolis term acquires anticyclonic relative vorticity. During its first days of life, the anticyclonic eddy was shallow Ro=0.9 and small (diameter less than 50 km). Later on, it reached a diameter of 150 km and a vertical structure of 3 km (Ro=0.1). A retrospective analysis with infrared images shows that the eddy's generation took place at about 3–4°E. Then, the eddy completed a counterclockwise circuit never before reported in other studies and ended up at the entrance of the Algerian Basin, where the interaction with the topography and the coastal instability induced its decay. The eddy's life span was 10 months. Computations of the heating rate following clusters of buoy trajectories show fluctuations throughout the eddy's journey, induced by advection and a seasonal warming.     

孔穴流激噪声的计算与验证研究

流激噪声的预报问题是流体动力声学领域的重点和难点.在气动声学领域,人们对此问题进行了很多探索,而在水声学领域,该问题尚缺乏深入的研究,文章的目的在于初步建立适用于水中流激噪声预报的数值方法.通过大涡模拟(LES)方法和FW-H声学类比方法,计算了两类孔穴的流激噪声问题.分析了声压谱级,壁面压力功率谱、涡量分布等,并利用试验结果对数值预报进行了验证.研究表明,计算结果量级可靠,符合声学一般规律.     

Using a random displacement model to simulate turbulent particle motion in a baroclinic frontal zone: A new implementation scheme and model performance tests

Numerical artifacts can limit accurate simulation of turbulent particle motion when Lagrangian particle-tracking models are implemented in hydrodynamic models with stratified conditions like fronts. Yet, modeling of individual particle motion in frontal regions is critical for understanding sediment dynamics as well as the transport and retention of planktonic organisms. The objective of this research was to develop a numerical technique to accurately simulate turbulent particle motions in a particle-tracking model embedded within a hydrodynamic model of a frontal zone. A new interpolation scheme, the ‘water column profile’ scheme, was developed and used to implement a random displacement model for turbulent particle motions. A new interpolation scheme was necessary because linear interpolation schemes caused artificial aggregation of particles where abrupt changes in vertical diffusivity occurred. The new ‘water column profile’ scheme was used to fit a continuous function (a tension spline) to a smoothed profile of vertical diffusivities at the x–y particle location. The new implementation scheme was checked for artifacts and compared with a standard random walk model using (1) Well Mixed Condition tests, and (2) dye-release experiments. The Well Mixed Condition tests confirmed that the use of the ‘water column profile’ interpolation scheme for implementing the random displacement model significantly reduced numerical artifacts. In dye-release experiments, high concentrations of Eulerian tracer and Lagrangian particles were released at the same location up-estuary of the salt front and tracked for 4 days. After small differences in initial dispersal rates, tracer and particle distributions remained highly correlated (r = 0.84 to 0.99) when a random displacement model was implemented in the particle-tracking model. In contrast, correlation coefficients were substantially lower (r = 0.07 to 0.58) when a random walk model was implemented. In general, model performance tests indicated that the ‘water column interpolation’ scheme was an effective technique for implementing a random displacement model within a hydrodynamic model, and both could be used to accurately simulate diffusion in a highly baroclinic frontal region. The new implementation scheme has the potential to be a useful tool for investigating the influence of hydrodynamic variability on the transport of sediment particles and planktonic organisms in frontal zones.     

Phytoplankton size classes competitions at sub-mesoscale in a frontal oceanic region

The effects of mesoscale and sub-mesoscale dynamics on the competition between two different phytoplankton size classes are investigated with a 3D primitive equations model. The model reproduces realistic simulations of mesoscale turbulence generated by a westward current in the southern hemisphere at statistical equilibrium in a summer situation. Effects of two different grazing pressures on phytoplankton competitions are compared and the role of eddy variability is quantified comparing high and low resolution simulations.High resolution simulations reveal a filamentary distribution of biomass and nutrients induced by the combination of vertical advection and horizontal stirring. This fine scale variability is observed not only on the horizontal but also on the vertical into the subsurface chlorophyll maximum.One of the key results is that such a dynamics induces a spatial segregation of the phytoplankton in the southern part of the frontal region that is mainly filamentary. This spatial segregation consists in biomass maxima for large phytoplankton in rich nutrients filaments and maxima for small phytoplankton outside these filaments. This anti-correlation is particularly strong when grazing pressure is low and is confirmed by statistical analysis. In the central frontal region, dominated by mesoscale dynamics, the two phytoplankton classes are strongly correlated together and biomass maxima are located close to downwelling regions that are poor in nutrients.It is shown that the effect of grazing is significantly amplified by the fine scale dynamics and that the combination of these two mechanisms is responsible of a switch of the ecosystem dominance in the surface layers.In addition, the effect of frontal dynamics on the detritus export is very sensitive to grazing pressure: increasing grazing induces a significant decrease of the export in the presence of frontal dynamics whereas it induces an increase of the export without small-scale variability.     

The influence of different turbulence schemes on modelling primary production in a 1D coupled physical–biological model

A one-dimensional (1D) coupled physical–microbiological model has been applied to a site in the central North Sea. The impact of the choice of the turbulence closure scheme on the modelling the primary production has been investigated.The model was run with four different parameterisations of vertical mixing of heat, momentum and dissolved and suspended matters, using M2 tidal forcing and the hourly mean meteorological forcing of 1989 to reproduce the annual thermal structure and primary production. The four mixing parameterisations are: Level 2 turbulence closure scheme [Mellor, G.L., Yamada, T., 1974. A hierarchy of turbulence closure models for planetary boundary layers. J. Atmos. Sci. 31, 1791–1806; Mellor, G.L., Yamada, T., 1982. Development of a turbulence closure model for geophysical Fluid problems. Rev. Geophys. Space Phys. 20 (4) 851–875] using an explicit numerical scheme [Sharples, J., Tett, P., 1994. Modelling the effect of physical variability on the midwater chlorophyll maximum. J. Mar. Res. 52, 219–238]; a version of the Level 2.5 turbulence closure scheme [Galperin, B., Kantha, L.H., Hassid, S., Rosati, A., 1988. A quasi-equilibrium turbulent energy model for geophysical flows. J. Atmos. Sci. 45, 55–62; Ruddick, K.G., Deleersnijder, E., Luyten, P.J., Ozer, J., 1995. Haline stratification in the rhine/meuse freshwater plume: a 3D model sensitivity analysis. Cont. Shelf Res. 15 (13) 1597–1630] simplified to use an algebraic mixing length by Sharples and Simpson [Sharples, J., Simpson, J.H., 1995. Semidiurnal and longer period stability cycles in the Liverpool Bay region of freshwater influence. Cont. Shelf Res. 15, 295–313], also solved explicitly; the same simplified L2.5 scheme with an implicit numerical solution and modified vertical discretisation scheme [Annan, J.D., 1999. Numerical methods for the solution of the turbulence energy equations in the shelf seas. Int. J. Numer. Methods Fluids 29, 193–206]; and another version of the same scheme (but using a different algebraic mixing length) as described by Xing and Davies [Xing, J., Davies, A.M., 1996a. Application of turbulence energy models to the computation of tidal currents and mixing intensities in the shelf edge regions. J. Phys. Oceanogr. 26, 417–447; Xing, J., Davies, A.M., 1996b. Application of a range of turbulence models to the computation of tidal currents and mixing intensities in shelf edge regions. Cont. Shelf. Res. 16, 517–547; Xing, J., Davies, A.M., 1998. Application of a range of turbulence energy models to the computation of the internal tide. Int. J. Numer. Methods Fluids 26, 1055–1084]. Various model outputs at the sea surface and in depth profiles have been compared with data collected in 1989 as part of the North Sea Project [Huthnance, J.M., 1990. Progress on North Sea Project. NERC News, vol. 12, pp. 25–29, UK]. It is shown that the biological results are extremely sensitive to the small changes in the physical conditions, which arise due to the different turbulence schemes tested. The timing of the spring bloom and the maintenance of the midwater chlorophyll maximum all differ greatly between model runs, and the gross primary production varies by a factor of two from the highest to lowest results. The simplified Level 2.5 scheme, implemented using the numerical methods of Annan [Annan, J.D., 1999. Numerical methods for the solution of the turbulence energy equations in the shelf seas. Int. J. Numer. Methods Fluids 29, 193–206], produces results, which give the best agreement with the available data.     

增大泥泵最大球形通道的方案

增大泥泵最大球形通道可以提升疏浚工程连续性,减少堵泵的风险,提高疏浚施工效率。泥泵最大球形通道和泥泵效率是泥泵设计中的矛盾点。针对上述问题通过调研提出了两种增大泥泵最大球形通道方案,采用数值模拟方法计算了3种方案的泥泵内部流场。通过对数值模拟结果的统计分析,对比3种方案泥泵的性能。确定了增大泥泵最大球形通道叶轮设计方案,为后续泥泵改造提供基础。     

Validation of a hybrid optimal interpolation and Kalman filter scheme for sea surface temperature assimilation

A hybrid data assimilation scheme designed for operational assimilation of satellite sea surface temperatures (SST) into an ocean model has been developed and validated against in-situ observations. The scheme consists of an optimal interpolation (OI) part and a greatly simplified Kalman filter (KF) part.The OI is performed only in the longitudinal and latitudinal directions. A climatological field is used as a background field for the interpolation. It is constructed by fitting daily averages of satellite SST to the annual mean, annual, and semiannual harmonics in a 20 km by 20 km grid. The background error covariance is approximated by a spatially varying two-dimensional exponential covariance model. The parameters of the covariance model are fitted to the deviations of the satellite data from the background field using data from a full year.The simplified KF uses ocean model forecasts as a background field. It is based on the assumption that it is possible to neglect horizontal SST covariances in the filter and that the typical time scale for vertical mixing in the mixed layer is much shorter than the average time between observations. We therefore assume that the error variance in a column of water is evenly spread out throughout the mixed layer. The result of these simplifications is a computationally very efficient KF.A one year validation of the scheme is performed for year 2001 using an operational eddy resolving ocean model covering the North Sea and the Baltic Sea. It is found that assimilation of sea surface temperature data reduces the model root mean square error from 1.13 °C to 0.70 °C. The hybrid scheme is found to reduce the root mean square error slightly more than the simplified KF without OI to 0.66 °C. The inclusion of spatially varying satellite error variances does not improve the performance of the scheme significantly.     

液化天然气动力船供气系统动态模拟

本文以液化天然气(LNG)动力船供气系统的工作流程为原型，讨论发动机在工况变化时如何设置控制方案，以达到更好的温度控制效果，满足双燃料发动机的工作需求。使用Aspen HYSYS流程模拟软件对该工作流程进行模拟，研究动态模拟时系统的运行情况。针对供气时的工况变化，改变液化天然气的流量，同时供气温度需要稳定在一定范围内，提出两种维持供气温度相对稳定的调节方案：改变热源的供热量或调节加热介质（水乙二醇）的流量，在HYSYS中分别进行动态模拟。改变热水流量时，输气温度可以控制在20℃-32℃；改变水乙二醇流量时，输气温度控制在22℃-28℃。两种调控方式结合的综合方案中，天然气温度也较好稳定在22℃-28℃之间，并且不需要引入其他外加冷源，响应较快，最适合应用于实际系统中。     

不同时空格式在求解Boussinesq水波模型中的应用

为了解不同时间空间差分格式在常用的二阶或四阶Boussinesq模型中的应用，针对4组近似到0（旷）阶完全非线性的二阶或四阶色散性的高阶Boussinesq水波方程，在非交错网格下，利用Crank-Nicolson格式、蛙跳格式、混合四阶Adams-Bashforth-Moulton格式，建立不同的数学模型。利用这些数值模型模拟波浪在潜堤上的传播变形，通过数值结果与试验结果的比较．考察时间格式及空间格式对模型的影响。结果表明：对同一方程，混合四阶Adams-Bashforth-Mouhon格式和Crank-Nicolson格式均能取得较好模拟效果，蛙跳格式的模拟效果最差；二阶Boussinesq模型采用追赶法求解已能满足要求；对四阶Boussineq模型，二阶空间导数色散项亦采用四阶精度，其数值效果会更好。     

淮河干流刘台子弯道水动力及河床变形数值模拟研究

为分析和评价刘台子弯道河段在不同整治方案下的工程效果,采用SIMHYD-2D平台构建二维泥沙数学模型对淮河干流临淮岗闸下—鲁台子河段2种整治方案进行水动力及河床变形数值模拟研究。结果表明:二维水沙数学模型能够很好地模拟河道的水力因子,也能合理地模拟自然条件下的冲淤分布,可以用于开展工程河段的冲淤演变预测分析;经过冲淤调整10 a后,在2种整治方案下,临鲁段河道主槽仍以冲刷为主,冲刷速度减缓,滩地仍以淤积为主,整体河道表现为冲刷;通过对比分析2种整治方案裁弯对上下游河道水沙运动及冲淤的影响、通航水流条件及新老河段未来冲淤情况,推荐方案2作为整治方案。