Application of a SEEK filter to a 1D biogeochemical model of the Ligurian Sea: Twin experiments and real in-situ data assimilation

The Singular Evolutive Extended Kalman (SEEK) filter has been implemented to assimilate in-situ data in a 1D coupled physical-ecosystem model of the Ligurian Sea. The biogeochemical model describes the partly decoupled nitrogen and carbon cycles of the pelagic food web. The GHER hydrodynamic model (1D version) is used to represent the physical forcings. The data assimilation scheme (SEEK filter) parameterizes the error statistics by means of a set of empirical orthogonal functions (EOFs). Twin experiments are first performed with the aim to choose the suitable experimental protocol (observation and estimation vectors, number of EOFs, frequency of the assimilation,…) and to assess the SEEK filter performances. This protocol is then applied to perform real data assimilation experiments using the DYFAMED data base. By assimilating phytoplankton observations, the method has allowed to improve not only the representation of the phytoplankton community, but also of other variables such as zooplankton and bacteria that evolve with model dynamics and that are not corrected by the data assimilation scheme. The validation of the assimilation method and the improvement of model results are studied by means of suitable error measurements.     

Data assimilation with the Ensemble Kalman Filter and the SEIK filter applied to a finite element model of the North Atlantic

Currently there are different approaches to filter algorithms based on the Kalman filter. One of the most used filter algorithms is the Ensemble Kalman Filter (EnKF). It uses a Monte Carlo approach to the filtering problem. Another approach is given by the Singular Evolutive Extended Kalman (SEEK) and Singular Evolutive Interpolated Kalman (SEIK) filters. These filters operate explicitly on a low-dimensional error space which is represented by an ensemble of model states. The EnKF and the SEIK filter have been implemented within a parallel data assimilation framework in the Finite Element Ocean Model FEOM. In order to compare the filter performances of the algorithms, several data assimilation experiments are performed. The filter algorithms have been applied with a model configuration of FEOM for the North Atlantic to assimilate the sea surface height in twin experiments. The dependence of the filter estimates on the represented error subspace is discussed. In the experiments the SEIK algorithm provides better estimates than the EnKF. Furthermore, the SEIK filter is much cheaper in terms of computing time.     

基础设计中的可靠度评估

阐述基于可靠度理论的基础设计的评估。首先总结具有结构抗力和荷载效应2个基础变量的结构可靠度基本理论。结构抗力的不确定性可以由统计学上的平均值和变异系数（cov或）Ω来描述。变异系数指变量的标准差和平均值的比值。敏感度分析的结果显示结构抗力的变异系数（ΩR）在其应用范围之内时,在分析结构可靠度方面扮演着相当重要的角色。基于这些阐述,在预先指定的风险水平（pf）上ΩR一定有其上限。当结构抗力呈正态分布的时候,这个极限ΩR独立于荷载效应随机性,和安全指数β成反比。安全指数可以定义为在标准正态分布区间极限状态到原点之间的最小距离。在这个极限ΩR之下,结构可以在预定风险水平之下安全工作。中心安全系数（FS）可以由结构抗力和荷载效应的变异系数根据平方关系求得。然而,一些情况下结构抗力为非正态分布的情况并不少见。因此,等效正态分布的概念可以用来得到非正态分布结构抗力的ΩR极限。地质方面的随机变量可能是正态分布,也可能是非正态分布,结构抗力中基本变量之间的关系可能是线性也可能是非线性,或者非常复杂以致于结构抗力只能通过有限元分析才能得到。在此情况下,随机数可以通过蒙特卡洛模拟技术获得。拟合的结构抗力的分布可以在随机试验的基础上通过配合度检验确定。现实中,土壤的力学特性不是各向同性的,同样也不能认为是单一材料的,它们的不确定性是不可以被忽略的。简便的设计方式认为不确定参数是常数,并且通过使用定值的安全系数来确定结构截面,设计原则没有将土壤参数对安全系数的影响考虑进去。参考计算出来的失效概率表明,确定值的安全系数法无法保证足够的安全。因此,在某种情况下,安全系数大于等于3,对于结构容许承受能力,并不能认为太保守。     

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.     

Hinterland strategies of port authorities: A case study of the port of Barcelona

Port authorities generally focus on the development of the local port area and play a minor role in the development of port hinterlands, whereas shippers, forwarders, barge and rail operators have always been involved in the port-hinterland connection. The increasing importance of intermodal hinterland networks for the competitive position of ports has urged port authorities to become active in the hinterland. This new role has already been suggested by different academics. However, limited empirical evidence exists of port authorities taking stakes in inland terminals or developing transport services. Barcelona, as one of the leading port authorities in this respect, is used as a case study in this paper. The case study provides insight in the components and execution of the hinterland strategy of Barcelona. It shows that the strategy of the port authority of Barcelona and the consequent active involvement in the hinterland has had a significant impact on attracting container volumes from distant hinterlands and improving the accessibility of the port.     

超高加速度抛射体内电子芯片的两维欧拉-拉格朗日模拟

文章建立两维欧拉—拉格朗日模型以研究抛射体发射过程对其内部电子芯片的冲击及其机理。该模型计入了火药燃烧的两相流动、抛射体结构运动及流固相互作用,其中火药的固体相和气体相分别采用了线性和指数状态方程描述。通过数值结果分析了以位移、弹口速度、加速度表征的抛射体过程中的冲击响应,并讨论了火药比能对抛射体的加速度冲击及对其内部电子芯片的冲击的影响规律。该模拟方法也可供船用智能火炮系统抗冲击分析及船上电子设备生存能力设计时参考。     

Effect of WO<Subscript>3</Subscript> nanoparticle loading on the microstructural,mechanical and corrosion resistance of Zn matrix/TiO<Subscript>2</Subscript>-WO<Subscript>3</Subscript> nanocomposite coatings for marine application

In this study, for marine application purposes, we evaluated the effect of process parameter and particle loading on the microstructure, mechanical reinforcement and corrosion resistance properties of a Zn-TiO₂-WO₃ nanocomposite produced via electrodeposition. We characterized the morphological properties of the composite coatings with a Scanning Electron Microscope (SEM) equipped with an Energy Dispersive Spectrometer (EDS). We carried out mechanical examination using a Dura Scan hardness tester and a CERT UMT-2 multi-functional tribological tester. We evaluated the corrosion properties by linear polarization in 3.5% NaCl. The results show that the coatings exhibited good stability and the quantitative particle loading greatly enhanced the structural and morphological properties, hardness behavior and corrosion resistance of the coatings. We observed the precipitation of this alloy on steel is greatly influenced by the composite characteristics.     

Effect of injection angle on artificial cavitation using the design of experiment method

Using the supercavitation phenomenon is necessary to reach high velocities underwater. Supercavitation can be achieved in two ways: natural and artificial. In this article, the simulation of flows around a torpedo was studied naturally and artificially. The validity of simulation using theoretical and practical data in the natural and artificial phases was evaluated. Results showed that the simulations were consistent with the laboratory results. The results in different injection coefficient rates, injection angles, and cavitation numbers were studied. The obtained results showed the importance of cavitation number, injection rate coefficient, and injection angle in cavity shape. At the final level, determining the performance conditions using the Design of Experiment (DOE) method was emphasized, and the performance of cavitation number, injection rate coefficient, and injection angle in drag and lift coefficient was studied. The increase in injection angle in the low injection rate coefficient resulted in a diminished drag coefficient and that in the high injection rate coefficient resulted in an enhanced drag coefficient.     

Scattering of oblique surface water waves by thin vertical barrier over undulating bed topography

The present study deals with the scattering of oblique surface water waves by small undulation on the bottom in the presence of a thin vertical barrier. Here, three different configurations of vertical barriers are investigated. Perturbation analysis is employed to determine the physical quantities, namely, the reflection and transmission coefficients. In this analysis, many different Boundary Value Problems (BVPs) are obtained out of which the first two bvps are considered. The zeroth order bvp is solved with the aid of eigenfunction expansion method. The first order reflection and transmission coefficients are derived in terms of the integrals by the method of the Green’s integral theorem. The variation of these coefficients is plotted and analyzed for different physical parameters. Furthermore, the energy balance relation, an important relation in the study of water wave scattering, is derived and checked for assuring the correctness of the numerical results for the present problem.     

The use of digital filter initialization to diagnose the mesoscale circulation and vertical motion in the California coastal transition zone

A dynamical method of initializing the primitive equations is tested and used to diagnose the three-dimensional circulation associated with jets and eddies as found in the California coastal transition zone (CTZ). The initialization method, referred to as digital filter initialization (DFI), was recently developed by [Monthly Weather Review 120 (1992) 1019] for use in an intermittent data assimilation system in the atmosphere. The ability of DFI to recover the mesoscale ageostrophic circulation associated with finite amplitude jets and eddies in the ocean is first demonstrated using control data produced by simulations with a primitive equation model. The DFI method is then applied to synoptic hydrographic data collected during several California CTZ surveys in the summer of 1988. The diagnostic results indicate the existence of jets, eddies, and filaments in the CTZ domain with maximum horizontal currents of the order of 0.6 m/s at the surface. Currents associated with such jets and filaments are coherent to a depth of over 500 m. The surface currents associated with a prominent cool filament are generally confluent, and weakly convergent on average, along the 270 km offshore extent of the filament. Meanders in the jet display convergence and downwelling upstream of pressure troughs and divergence and upwelling downstream of the troughs. Maximum vertical velocities at 100 m are of the order of 10 m/day. This result is consistent with independent estimates of subduction rates made from biological studies in this and similar coastal filaments in the CTZ program.