共查询到20条相似文献，搜索用时 296 毫秒
1.
在悬沙输运的数值模拟中,初始场的准确给定至关重要。目前诸多确定初始场的方案均存在一定的缺陷,初始场的准确性有待进一步提高。本文基于一个三维悬沙输运伴随同化模型,通过孪生实验和实际实验,对模型初始场进行了伴随法反演研究。在孪生实验中,首先验证了初始场的相对重要性;其次,探讨了初始场的反演结果对优化算法、初始猜测值、卫星遥感数据数量、同化时间窗口宽度和背景流场误差的敏感性;最后,比较了伴随法和插值法重构初始场的能力。孪生实验结果表明:最速下降法对初始场的优化反演效果要优于三种共轭梯度法和有限记忆BFGS法;初始场的反演效果对初始猜测值、卫星遥感数据数量和背景流场误差不敏感,而对同化窗口宽度较为敏感;与插值法相比,伴随法是重构模型初始场更有效的手段。实际实验中,在杭州湾海域同化典型的小潮时期和大潮时期的GOCI卫星遥感资料所得表层悬沙浓度数据,优化反演了初始场。实际实验结果表明:数据同化后,得到了更符合实际的最优初始场,表明伴随法是实现初始场优化反演的有效手段。该研究对进一步改进悬沙输运模型的初始化方案具有一定的参考价值,也对其他数值模型的初始化方案具有一定的借鉴价值。 相似文献
2.
3.
The behavior of a highly deformable membrane to ocean waves was studied by coupling a nonlinear boundary element model of the fluid domain to a nonlinear finite element model of the membrane. The hydrodynamic loadings induced by water waves are computed assuming large body hydrodynamics and ideal fluid flow and then solving the transient diffraction/radiation problem. Either linear waves or finite amplitude waves can be assumed in the model and thus the nonlinear kinematic and dynamic free surface boundary conditions are solved iteratively. The nonlinear nature of the boundary condition requires a time domain solution. To implicitly include time in the governing field equation, Volterra's method was used. The approach is the same as the typical boundary element method for a fluid domain where the governing field equation is the starting point. The difference is that in Volterra's method the time derivative of the governing field equation becomes the starting point.The boundary element model was then coupled through an iterative process to a finite element model of membrane structures. The coupled model predicts the nonlinear interaction of nonlinear water waves with highly deformable bodies. To verify the coupled model a large scale test was conducted in the OH Hinsdale wave Research Laboratory at Oregon State University on a 3-ft-diameter fabric cylinder submerged in the wave tank. The model data verified the numerical prediction of the structure displacements and of the changes in the wave field.The boundary element model is an ideal modeling technique for modeling the fluid domain when the governing field equations is the Laplace equation. In this case the nonlinear boundary element model was coupled with a finite element model of membrane structures, but the model could have been coupled with other finite element models of more rigid structures, such as a pontoon floating breakwater. 相似文献
4.
近岸局部范围内的波浪是由远处洋面上生成并传播过来的，它的生成发展源地(敏感性海域)对于所关注范围的波浪起决定性影响。共轭方程中的模式变量反映了物理参量的梯度变化，利用它可以达到精确确定敏感海域具体位置的目的。本文在LAGFD-WAM海浪模式基础上建立了海浪谱共轭方程，提出了通过在关注点邻域加一小扰动作为逆向积分海浪谱共轭方程的驱动项，来计算出梯度极值分布的思想。将这一思想应用于东中国海区域，通过一系列试验来确定舟山及长江口附近海域的敏感性区域位置。试验结果符合物理意义，梯度极值分布结果为下一步的多源卫星观测资料复合分析及最优实测方案设计奠定了基础。 相似文献
5.
6.
7.
8.
The construction of stationary expressions for quantities of physical interest such as radiated power and target strength is discussed broadly for acoustic problems involving radiation or scattering from finite objects of arbitrary shape. The Kirchhoff-Helmholtz integral corollaries of the wave equation, which express acoustic pressure at either interior or exterior points in terms of pressure and its normal derivative over any closed surface, yield for both interior and exterior problems two mathematically dissimilar but related functional relations between surface field quantities. One of these is the better known surface Helmholtz integral equation; the other is a differential-integral relation which involves the tangential derivatives of pressure on the surface. The four linear operators involved in these functional relations are studied and it is found that two are self-adjoint, while the other two are an adjoint pair. A general technique for constructing variational expressions recently developed by Gerjuoy et al. [28] is adapted to acoustic radiation and scattering problems with the functional relations taken as the primary governing relations. Included examples are stationary expressions for the radiated power when either the normal velocity or the pressure are specified on the surface (the other quantity being unknown) and the target strength for scattering from a rigid object. The adjoint relations allow simple physical interpretations for the Lagrange multipliers that arise in the theory, such that the guesses for good trial functions can take advantage of existing physical insight. It is demonstrated with a specific example (transversely vibrating disk) that the resulting estimate for radiated power is substantially more accurate than that of the trial function for surface pressure which was inserted into the stationary expression. 相似文献
9.
特征线计算格式下共轭方程两种导出途径的比较 总被引：1，自引：0，他引：1
共轭方程的导出是建立资料同化模型的关键，其导出方式有两种途径：AFD形式与FDA形式。在特征线计算格式基础上针对一类较广泛海洋动力控制方程分析了其两种共轭方程（AFD形式与FDA形式）之间的关系，并将理论结果应用于波谱共轭方程的讨论。 相似文献
10.
三维斜压台风模式 Ⅱ.预报试验 总被引：1，自引：0，他引：1
一种斜压多重移动套网格台网模式在国家海洋环境预报中心已被应用于海洋环境预报。本文第一部分已描述了模式方程组和数值方法。本文继续概述模式网格、变分调整初始化和预报试验结果。最外粗网格域固定,内部细网格域随台风中心轨迹移动。模式中,粗细网格变量采用双向耦合。平衡方程和方程,理想台风场和变分调整方案被用于台风模式初始化。一种简单而有效的资料同化方法,即用第6h台风报和弱约束变分原理调整初始场,被应用于提高预报结果。最后本文给出预报试验结果。预报误差统计显示本模式对台风路径预报具有相当能力,而且可以提供海面风和气压场较好的预报。本模式已经与海浪模式联结,得到满意的波高预报结果。 相似文献
11.
Estimation of Land Subsidence Based on Groundwater Flow Model 总被引：8，自引：0，他引：8
This article presents an approach for estimating land subsidence due to withdrawal of groundwater. The proposed method calculates the groundwater seepage in 3-D-condition and calculates the land subsidence one-dimensionally. The governing equation on groundwater seepage is based on the three-dimensional mass conservation law and the principle of effective stress. The land subsidence calculation method is derived based on the following assumptions: (1) displacements occur only in the vertical direction, and (2) in vertical direction the total stresses do not change. The governing equation is solved by numerical method, i.e., finite element method (FEM) in spatial discretization and finite difference method (FDM) in time series discretization. In FEM Galerkin method is adopted and in FDM, lumped matrix method is employed. The proposed method is calibrated via analyzing 1-D consolidation problem and the results are compared with those from Terzaghi's one-dimensional consolidation theory and oedometer test. The proposed method is employed to analyze the consolidation of a soft layer due to withdrawal of groundwater from an aquifer under it. Moreover, this method is also applied to a field case of land subsidence due to groundwater pumping in a gas production field in Japan. The analytical results are compared with the field observed data. The results show that this approach simulates the field case well. 相似文献
12.
This article presents an approach for estimating land subsidence due to withdrawal of groundwater. The proposed method calculates the groundwater seepage in 3-D-condition and calculates the land subsidence one-dimensionally. The governing equation on groundwater seepage is based on the three-dimensional mass conservation law and the principle of effective stress. The land subsidence calculation method is derived based on the following assumptions: (1) displacements occur only in the vertical direction, and (2) in vertical direction the total stresses do not change. The governing equation is solved by numerical method, i.e., finite element method (FEM) in spatial discretization and finite difference method (FDM) in time series discretization. In FEM Galerkin method is adopted and in FDM, lumped matrix method is employed. The proposed method is calibrated via analyzing 1-D consolidation problem and the results are compared with those from Terzaghi's one-dimensional consolidation theory and oedometer test. The proposed method is employed to analyze the consolidation of a soft layer due to withdrawal of groundwater from an aquifer under it. Moreover, this method is also applied to a field case of land subsidence due to groundwater pumping in a gas production field in Japan. The analytical results are compared with the field observed data. The results show that this approach simulates the field case well. 相似文献
13.
A variational formulation for three-dimensional analysis of extensible marine riser transporting fluid 总被引：1，自引：0，他引：1
Chainarong Athisakul 《Ocean Engineering》2011,38(4):609-620
This paper presents a model formulation for static and dynamic analysis of three-dimensional extensible marine riser transporting fluid. A variational model formulation is developed based on the principle of virtual work-energy and the extensible elastica theory. The virtual work-energy functional is composed of the virtual strain energy due to axial stretching, bending, and torsion and the virtual work done by the external and internal fluid. The governing dynamic equilibrium equations are derived in the Cartesian coordinate. The finite element method is used to obtain the numerical solutions. The numerical examples are provided to demonstrate interesting effects of fluid transportation and axial deformation on large displacements and dynamic properties of the three-dimensional extensible marine riser. 相似文献
14.
15.
The elliptic mild slope equation is used to simulate linear wave propagation over variable seabed topography with mild slopes.The governing equation is discretized by the finite difference method.Ba-sed on the BI-CGSTAB technique.an attractive variant of BI-Conjugate Gradients(BI-CG)method,theobtained linear algebraic system of equations is solved.Numerical experiments show that the BI-CGSTABmethod is efficient for solving the elliptic mild slope equation.The results obtained by the BI-CGSTAB-Ba-sed method are much the same as those obtained by other authors with different solution methods,but theconvergence rate is much faster than that of other methods. 相似文献
16.
传统的四维变分数据同化方法在同化观测资料的同时可以对数值模式参数进行优化,然而传统的四维变分方法需要针对不同的数值模式编写特有的伴随模式,因此算法的可移植性差,同时计算时耗费大量资源。本文提出了一种新的基于解析四维集合变分的参数优化方法,该方法以迭代搜索得到的模式参数为基准展开扰动并构建样本集合,由此显式地计算协方差矩阵,并得到代价函数极小值的解析解,从而避免了伴随模式的使用。基于Lorenz-63模型对该方法进行单参数和多参数数值试验和优化效果检验,并在不同的同化时间窗口长度和观测采样间隔情况下,采用传统四维变分方法与之进行对比,结果显示,新方法表现出与传统四维变分相同的优化性能,都能有效收敛到真值,而新方法不需要计算伴随模式,可移植性好。本文还测试了不同的集合成员个数和模式参数真值的情况下新方法的同化效果,结果表明,新方法对集合样本个数及模型参数真值不敏感,采用较少的集合样本即可完成数据同化。 相似文献
17.
本文基于时间分布参数设置,利用伴随同化方法,反演了Ekman模型中随时间变化的风应力拖曳系数,并在孪生实验和实际实验中对该方法进行了验证。在孪生实验中,研究了参数反演结果对不同影响因素的响应,包括:风速分布、风应力系数分布、风应力系数初始猜测值、风应力系数独立变量个数、观测数据误差和观测的深度。孪生实验结果验证了伴随同化方法反演Ekman模型中时变风应力系数的有效性,具体包括如下五个方面结论:1)不同风速分布下均能成功反演出不同风应力拖曳系数分布; 2)反演结果对初始猜测值较为敏感,风应力系数初始猜测值越接近给定值,反演结果越好;3)风应力系数独立点个数的选取会显著影响反演结果,合理的选择有利于提高反演效率及减小观测数据误差;4)观测误差能够影响反演结果,观测数据误差在20%以下时能取得合理的反演结果; 5)反演结果对观测数据的表层和次表层流速更为敏感,这是由Ekman流的物理性质决定的。实际实验,利用百慕大锚系试验平台的风速和流速数据,去除周期性潮流和地转流成分后得到Ekman流成分,并作为观测输入到该同化模型,反演出了适用于该区域和该时段的随时间变化的风应力系数。通过比较模拟流速和观测流速,证明利用伴随同化方法能从实测数据中反演出合理的时变风应力系数,对于海洋模型风应力系数的确定是一项有益的尝试。 相似文献
18.
The study investigates perspectives of the parameter estimation problem with the adjoint method in eddy-resolving models. Sensitivity to initial conditions resulting from the chaotic nature of this type of model limits the direct application of the adjoint method by predictability. Prolonging the period of assimilation is accompanied by the appearance of an increasing number of secondary minima of the cost function that prevents the convergence of this method. In the framework of the Lorenz model it is shown that averaged quantities are suitable for describing invariant properties, and that secondary minima are for this type of data transformed into stochastic deviations. An adjoint method suitable for the assimilation of statistical characteristics of data and applicable on time scales beyond the predictability limit is presented. The approach assumes a greater predictability for averaged quantities. The adjoint to a prognostic model for statistical moments is employed for calculating cost function gradients that ignore the fine structure resulting from secondary minima. Coarse resolution versions of eddy-resolving models are used for this purpose. Identical twin experiments are performed with a quasigeostrophic model to evaluate the performance and limitations of this approach in improving models by estimating parameters. The wind stress curl is estimated from a simulated mean stream function. A very simple parameterization scheme for the assimilation of second-order moments is shown to permit the estimation of gradients that perform efficiently in minimizing cost functions. 相似文献
19.
提出了利用共轭方程研究海洋对于局域大气平均气温的作用的方法。在现今海洋基本动力方程基础上导出了球坐标系下的海水温度共辊方程。分析了海水共轭温度的意义，其量值(或垂直梯度量值)表征了不同区域海洋对于局域气温的相应贡献，它在全球海洋中的分布对于局域气候形成与变化研究具有重要意义。 相似文献
20.
WANG Lixia SUN Yinglan ZHENG Lianyuan Associate Professor Insititute of Environmental Oceanography Ocean University of Qingdao Qingdao China.Professor Insititute of Environmental Oceanography Ocean University of Qingdao Qingdao China. Lecturer Insititute of Environmental Oceanography Ocean University of Qingdao Qingdao China. 《中国海洋工程》1998,(3)
A three-dimensional,first order turbulence closure,thermal diffusion model is described inthis paper.The governing equations consist of an equation of continuity,three components of momentum,conservation equations for salt,temperature and subgridscale energy,and an equation of state.In the mod-el,according to the hypothesis of Kolmogorov and Prandtl,the viscosity coefficient of turbulent flow ofhomogeneous fluid is related to the local turbulent energy,and the horizontal and vertical exchangecoefficients of mass,heat and momentum are computed with the introduction of subgridscale turbulenceenergy.The governing equations are solved by finite difference techniques.This model is applied to theJiaozhou bay to predict thermal pollution by the Huangdao power plant.An instantaneous tidal currentfield is computed,then the distribution of temperature increment is predicted,and finally the effect of windstress on thermal discharge is discussed. 相似文献