Ocean current estimation using a Multi-Model Ensemble Kalman Filter during the Grand Lagrangian Deployment experiment (GLAD)

In the summer and fall of 2012, during the GLAD experiment in the Gulf of Mexico, the Consortium for Advanced Research on Transport of Hydrocarbon in the Environment (CARTHE) used several ocean models to assist the deployment of more than 300 surface drifters. The Navy Coastal Ocean Model (NCOM) at 1 km and 3 km resolutions, the US Navy operational NCOM at 3 km resolution (AMSEAS), and two versions of the Hybrid Coordinates Ocean Model (HYCOM) set at 4 km were running daily and delivering 72-h range forecasts. They all assimilated remote sensing and local profile data but they were not assimilating the drifter’s observations. This work presents a non-intrusive methodology named Multi-Model Ensemble Kalman Filter that allows assimilating the local drifter data into such a set of models, to produce improved ocean currents forecasts. The filter is to be used when several modeling systems or ensembles are available and/or observations are not entirely handled by the operational data assimilation process. It allows using generic in situ measurements over short time windows to improve the predictability of local ocean dynamics and associated high-resolution parameters of interest for which a forward model exists (e.g. oil spill plumes). Results can be used for operational applications or to derive enhanced background fields for other data assimilation systems, thus providing an expedite method to non-intrusively assimilate local observations of variables with complex operators. Results for the GLAD experiment show the method can improve water velocity predictions along the observed drifter trajectories, hence enhancing the skills of the models to predict individual trajectories.     

Investigation of drag crisis phenomenon using CFD methods

When fluid flow passes a cylinder, the drag crisis phenomenon occurs between the sub-critical and the super-critical Reynolds numbers. The focus of the present studies was on the numerical prediction of the drag crisis based on CFD methods. In this work, block structured meshes with refined grids near the cylinder surface and in the downstream were employed. Both 2D and 3D simulations were performed using various turbulence models, including the SST k − ω model, the k − ϵ model, the SST with LCTM, the DES model, and the LES model. In the convergence studies, the effects of the grid size, the time step, the first grid size and the aspect ratio (for 3D simulations) on the solutions were examined. The errors due to spatial and time discretizations were quantified according to a V&V procedure. Validation studies were carried out for various Reynolds numbers between Re = 6.31 × 10⁴ and 7.57 × 10⁵. The averaged drag force, the RMS of lift force and the Strouhal number were compared with experimental data. The studies indicated that standard 2D and 3D RANS methods were inadequate to capture the drag crisis phenomenon. The LES method however has the potential to address the problem.     

海洋沉积地质过程模拟:性质与问题及前景

探讨了海洋沉积地质数值模拟的功能问题.地球科学研究以正演方法(以数值模拟为代表)和反演方法(以地层记录分析为代表)的结合为特征.数值模拟作为一种重要的研究工具,具有获得过程和机理分析结果、形成工作假说、指导现场观测和采样的作用.与此同时,在正反演交融中模型本身也得到了提高和完善.在海洋沉积地质领域,数值模拟在沉积物输运...     

Numerical modeling of hyperpycnal flows in an idealized river mouth

Numerical experiments in an idealized river mouth are conducted using a three-dimensional hydrodynamics model (EFDC model) to examine the impacts of suspended sediment concentration (SSC), settling velocity of sediment and tidal mixing on the formation and maintenance of estuarine hyperpycnal flows. The standard experiment presents an illustrative view of hyperpycnal flows that carry high-concentrated sediment and low-salinity water in the bottom layer (>1.0 m in thickness) along the subaqueous slope. The structure and intra-tidal variation of the simulated hyperpycnal flows are quite similar to those previously observed off the Huanghe (Yellow River) mouth. Results from the three control experiments show that SSC of river effluents is the most important parameter to the formation of hyperpycnal flows. High SSC will increase the bulk density of river effluents and thus offset the density difference between freshwater and seawater. Low SSC of river effluents will produce a surface river plume, as commonly observed in most large estuaries. Both the settling velocity of sediment particles and the tidal mixing play an important role in maintaining the hyperpycnal flows. Increasing settling velocity enhances the deposition of sediment from the hyperpycnal layer and thus accelerates the attenuation of hyperpycnal flows, whereas increasing tidal mixing destroys the stratification of water column and therefore makes the hyperpycnal flows less evident. Our results from numerical experiments are of importance to understand the initiation and maintenance of hyperpycnal flows in estuaries and provide a reference to the rapidly decaying hyperpycnal flows off the Huanghe river mouth due to climatic and anthropogenic forcing over the past several decades.     

渤海悬浮颗粒物的三维输运模式 I.模式

为了长时间,大范围研究渤海中悬浮物的输运规律,研究应用并改造了备国汉堡大学的粒子追踪悬浮物输运的三维模式,模式中考虑了风、气压等气象要素以外海传入的潮波作用下的三维正压海流对悬浮颗粒物的输运．同时考虑了悬浮颗粒物的学降及再悬浮机制和底质中的细颗粒物的运动,将风浪的作用引入了悬浮物输运的计算,数值模式应用的粒五追踪方法,能较好地反映悬浮物浓度的迅速变化。     

应用数值模拟方法探讨河口最大浑浊带若干机理

利用平面二维潮流方程结合悬沙输运方程模拟了河口最大浑浊带现象。并且对河口边界进行进一步概化 ,比较了恒定流与非恒定流、稳定源与非稳定源、矩形河口与线性河口等不同条件下河口悬沙浓度的平面分布特点。结果表明河口地形边界和非恒定潮流作用对河口最大浑浊带的悬沙富集有重要贡献。     

The Monagas Fold-Thrust Belt of Eastern Venezuela. Part I: Structural and thermal modeling

The Eastern Venezuelan Basin (EVB) contains one of the largest hydrocarbon accumulations in the world. Main petroleum targets are buried structures of the Monagas Fold-Thrust Belt (MFTB) which forms the northeastern edge of the EVB. The objective of this study is to integrate the seismic and well data that has been acquired over the last 10 years across the MFTB and EVB, to create an updated structural model. Three regional cross sections 60-75 km long are presented across an area of 4000 km².Five structural domains are described: Amarilis, Furrial, Jusepín, Cotoperí and Pirital. They are characterized by thrusts and high-angle reverse faults. Structural style changes along strike are related to variations in depth of detachment levels and to the strike-slip component of the deformation. We have estimated a shortening between 43 and 59 km that increases eastward over a distance of 40 km.The evolution of the MFTB is divided in four episodes based on stratigraphic, structural and thermal maturity evidences: (1) Oligocene-early Miocene initial movement of Pirital thrust. (2) Early Miocene simultaneous movement on Pirital, Furrial and Cotoperí thrusts. (3) Middle Miocene increases in velocity and change in geometry of Pirital thrust, during an out of sequence period of thrusting. (4) Late Miocene to Holocene minor thrust activity. This evolution is consistent with the oblique convergence between the Caribbean and South American plates and the convergence between North and South America that affected Eastern Venezuela during the Cenozoic.By analyzing the along-strike variations in structural style, new exploratory opportunities have been identified. Under the Orocual and Santa Bárbara fields two untested duplex structures are proposed; they were developed during the middle Miocene. Other prospective hydrocarbon traps are associated to oblique transpressive faults that create anticline structures.     

Numerical simulation of groundwater flowing to horizontal seepage wells under a river

A horizontal seepage well, consisting of an interconnected vertical well, galleries, chambers and small-diameter radiating bores, is used to acquire relatively clean water that has been filtered through natural alluvial deposits in a riverbed. It has wide application, especially in arid and semi-arid areas. The lack of calculation formulae or models for horizontal seepage wells, up until now, has resulted in several false applications. Based on the analysis of groundwater flow characteristics, it has been concluded that several flow regimes coexist and hydraulic head loss exists in the horizontal seepage well. To avoid the difficulty of confirming the flux or head distribution in such a complex system, the model boundary of the whole horizontal seepage well has been moved to that of just the vertical well, and the well-aquifer system was treated as a heterogeneous medium, where the horizontal seepage well itself is a highly permeability medium. A mathematical model has been developed, based on the coupled seepage-pipe flow, by the introduction of equivalent hydraulic conductivity according to different flow regimes. Then a three-dimensional finite difference numerical model, based on the mathematical model, was developed and applied to a horizontal seepage well in China. The numerical model verified the groundwater flow characteristics of the horizontal seepage well. An erratum to this article can be found at     

Numerical algorithm for the dynamic analysis of base-isolated structures with dry friction

This paper deals with the techniques for solving ordinary differential equations with essential nonlinearity arising from the representation of frictional force by the sign function of the relative velocity. This problem is connected with dynamic behavior studies of base-isolated structures with dry friction devices.Mathematical evidence is presented concerning the inaccuracies of the solution which are connected with the representation of the frictional force by the sign function. An approximation of this sign model is made on the basis of the appropriate linearization of a small part of the frictional force function. A numerical algorithm for the dynamic analysis of base-isolated structures with dry friction is proposed. Numerical tests are carried out for establishing the exactness of the proposed techniques. The obtained results prove the accuracy of the techniques.     

Generalized plasticity state parameter‐based model for saturated and unsaturated soils. Part 1: Saturated state

The behavior of granular materials is known to depend on its loose or dense nature, which in turns depends both on density and confining pressure. Many models developed in the past require the use of different sets of constitutive parameters for the same material under different confining pressures. The purpose of this paper is to extend a basic generalized plasticity model for sands proposed by Pastor, Zienkiewicz and Chan by modifying the main ingredients of the model flow—rule, loading–unloading discriminating direction and plastic modulus—to include a dependency on the state parameter. The proposed model is tested against the available experimental data on three different sands, using for each of them a single set of material parameters, finding a reasonably good agreement between experiments and predictions. Copyright © 2010 John Wiley & Sons, Ltd.