Effective block diagonal preconditioners for Biot's consolidation equations in piled‐raft foundations     

K. B. Chaudhary  K. K. Phoon  K. C. Toh 《国际地质力学数值与分析法杂志》2013,37(8):871-892

The finite element (FE) simulation of large‐scale soil–structure interaction problems (e.g. piled‐raft, tunnelling, and excavation) typically involves structural and geomaterials with significant differences in stiffness and permeability. The symmetric quasi‐minimal residual solver coupled with recently developed generalized Jacobi, modified symmetric successive over‐relaxation (MSSOR), or standard incomplete LU factorization (ILU) preconditioners can be ineffective for this class of problems. Inexact block diagonal preconditioners that are inexpensive approximations of the theoretical form are systematically evaluated for mitigating the coupled adverse effects because of such heterogeneous material properties (stiffness and permeability) and because of the percentage of the structural component in the system in piled‐raft foundations. Such mitigation led the proposed preconditioners to offer a significant saving in runtime (up to more than 10 times faster) in comparison with generalized Jacobi, modified symmetric successive over‐relaxation, and ILU preconditioners in simulating piled‐raft foundations. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

非线性层结与对流运动的稳定性     

陈忠明 《成都信息工程学院学报》1991,(Z1)

本文讨论了温度层结为普遍化三次分布情况下对流运动的分岔和突变特征,得到了一些有意义的结果,加深了我们对非线性对流运动特征的一些认识。  相似文献   

Simulation of nonlinear wave run-up with a high-order Boussinesq model   总被引：2，自引：0，他引：2  

David R. Fuhrman  Per A. Madsen 《Coastal Engineering》2008

This paper considers the numerical simulation of nonlinear wave run-up within a highly accurate Boussinesq-type model. Moving wet–dry boundary algorithms based on so-called extrapolating boundary techniques are utilized, and a new variant of this approach is proposed in two horizontal dimensions. As validation, computed results involving the nonlinear run-up of periodic as well as transient waves on a sloping beach are considered in a single horizontal dimension, demonstrating excellent agreement with analytical solutions for both the free surface and horizontal velocity. In two horizontal dimensions cases involving long wave resonance in a parabolic basin, solitary wave evolution in a triangular channel, and solitary wave run-up on a circular conical island are considered. In each case the computed results compare well against available analytical solutions or experimental measurements. The ability to accurately simulate a moving wet–dry boundary is of considerable practical importance within coastal engineering, and the extension described in this work significantly improves the nearshore versatility of the present high-order Boussinesq approach.  相似文献   

The Wiener–Levinson algorithm and ill-conditioned normal equations     

R. J. O'Dowd 《Geophysical Journal International》1991,106(2):399-406

  相似文献   

附有限制条件的间接平差秩亏时解法初探     

秦永宽  黄声享  张书毕 《海洋测绘》2009,29(1):7-9

在附有限制条件的间接平差的应用中，若所列观测方程中实际用到的未知参数的个数少于所选参数的总数u时，就会出现法方程系数阵秩亏的问题。为了解决此问题，在理论公式推导的基础上，给出了一种新的解决方法，并结合算例进行了验证。  相似文献   

A two-dimensional finite volume hydrodynamic model for coastal areas: Model development and validation     

A.H. Daoud  K.A. Rakha  A.G. Abul-azm   《Ocean Engineering》2008,35(1):150-164

A two dimensional implicit finite volume scheme for solving the shallow-water equations is developed. The effects of the Coriolis force, surface wind stress, and waves are included. A non-uniform rectilinear forward staggered grid is used with Cartesian coordinates. The time integration is performed using the Euler implicit technique. The convective flux is treated using the deferred correction method. The viscous terms are discretized using a second order central difference approximation. The SIMPLE (Semi-Implicit Method for Pressure-Linked Equations) algorithm is used for coupling the velocity components and the water elevation gradient for the water level correction. The system of equations is solved sequentially using the Strongly Implicit Procedure (SIP). To simulate wave driven current, a phase averaged wave model is used first to simulate wave transformation and calculate radiation stresses. The performance of the developed model is validated for different sources of external forces and different combinations of boundary conditions. The validation cases include tidal circulation in a harbor and wave induced currents behind a breakwater parallel to the coastline. The model is finally applied to simulate the flow pattern in a closed artificial lagoon and along the coastline near Damietta Port located along the Northern coast of Egypt. Results of the developed model agree well with the published results for the considered cases.  相似文献   

桂林城区大气能见度与颗粒物浓度和气象因子关系研究     

龙凤翔  张瑀琳 《气象与环境学报》2020,36(1):21-27

利用2015年1月至2017年6月桂林国家基本气象站能见度、相对湿度、气温、气压、降水等气象要素和PM₁₀、PM_2.5、PM_1.0颗粒物质量浓度资料,分析桂林城区大气能见度与颗粒物浓度和气象因子之间关系。结果表明：桂林城区大气能见度和PM₁₀、PM_2.5、PM_1.0呈对数关系,相关系数分别为-0.341、-0.461、-0.509,颗粒物对大气能见度影响在相对湿度为60%—70%时最为显著。在各气象因子中,大气能见度与风速的相关性最好,其次为相对湿度,与风速呈二次函数关系,与相对湿度呈幂指数关系,与温度相关性较小,与气压在秋冬季节呈正相关,相关系数冬季可达0.301,但在春、夏季节相关性不显著;利用颗粒物浓度和气象要素建立8种大气能见度非线性统计回归模型,比较后发现利用PM_1.0、风速、相对湿度、气温等因子建立的不同季节大气能见度拟合公式在实际检验中效果最优,能较好地模拟桂林地区大气能见度的变化。  相似文献   

Parameter estimation for basin-scale ecosystem-linked population models of large pelagic predators: Application to skipjack tuna   总被引：2，自引：1，他引：1  

Inna Senina  John Sibert  Patrick Lehodey 《Progress in Oceanography》2008,78(4):319-335

A Spatial Ecosystem and Population Dynamic Model (SEAPODYM) is used in a data assimilation study aiming to estimate model parameters that describe dynamics of Pacific skipjack tuna population on ocean-based scale. The model based on advection–diffusion–reaction equations explicitly predicts spatial dynamics of large pelagic predators, while taking into account data on several mid-trophic level components, oceanic primary productivity and physical environment. In order to improve its quantitative ability, the model was parameterized through assimilation with commercial fisheries data, and optimization was carried out using maximum likelihood estimation approach. To address the optimization task we implemented an adjoint technique to obtain an exact, analytical evaluation of the likelihood gradient. We conducted a series of computer experiments in order to (i) determine model sensitivity with respect to variable parameters and, hence, investigate their observability; (ii) estimate observable parameters and their errors; and (iii) justify the reliability of the computed solution. Parameters describing recruitment, movement, habitat preferences, natural and fishing mortality of skipjack population were analysed and estimated. Results of the study suggest that SEAPODYM with achieved parameterization scheme can help to investigate the impact of fishing under various management scenarios, and also conduct forecasts of a given species stock and spatial dynamics in a context of environmental and climate changes.  相似文献   

Lambert problem solution in the hill model of motion   总被引：1，自引：0，他引：1  

Alexander Sukhanov  Antonio F. Bertachini A. Prado 《Celestial Mechanics and Dynamical Astronomy》2004,90(3-4):331-354

The goal of this paper is obtaining a solution of the Lambert problem in the restricted three-body problem described by the Hill equations. This solution is based on the use of pre determinate reference orbits of different types giving the first guess and defining the sought-for transfer type. A mathematical procedure giving the Lambert problem solution is described. This procedure provides step-by-step transformation of the reference orbit to the sought-for transfer orbit. Numerical examples of the procedure application to the transfers in the Sun–Earth system are considered. These examples include transfer between two specified positions in a given time, a periodic orbit design, a halo orbit design, halo-to-halo transfers, LEO-to-halo transfer, analysis of a family of the halo-to-halo transfer orbits. The proposed method of the Lambert problem solution can be used for the two-point boundary value problem solution in any model of motion if a set of typical reference orbits can be found.  相似文献   

Behaviour of the SMF method for the numerical integration of satellite orbits     

Pablo Martín  José M. Ferrándiz 《Celestial Mechanics and Dynamical Astronomy》1995,63(1):29-40

The SMF algorithms were recently developed by the authors as a multistep generalization of the ScheifeleG-functions one-step method. Like the last, the proposed codes integrate harmonic oscillations without truncation error and the perturbing parameter appears as a factor of that error when integrating perturbed oscillations. Therefore they seemed to be convenient for the accurate integration of orbital problems after the application of linearizing transformations, such as KS or BF. In this paper we present several numerical experiments concerning the propagation of Earth satellite orbits, that illustrate the performance of the the SMF method. In general, it provides greater accuracy than the usual standard algorithms for similar computational cost.  相似文献