首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 125 毫秒
1.
海洋水温垂直分布数据同化方法研究   总被引:5,自引:1,他引:5  
以一维海洋水温模型为例,利用伴随法进行海洋观测数据同化试验,以便为水温的数值预报提供较准确的初始场.文中利用泛函的Gâteaux微分和Hilbert空间上伴随算子的概念讨论了连续的伴随模型的建立,并通过选择适当的差分格式离散伴随模型,使其保持连续时的伴随关系,同时给出了水温初始场最优化过程及相应的同化试验数值结果.  相似文献   

2.
针对一组近似到二阶完全非线性,四阶色散的Boussinesq方程,在交错网格下建立了数学模型.计算中时间层不交错,模型的求解利用混合四阶Adams-Bashforth-Moulton格式的有限差分法.数值模拟了波浪在潜堤上的演化过程,再现了波浪的浅化、反射以及非线性波能量传递等现象.对数值计算结果采用Friouer变换...  相似文献   

3.
在非交错网格下采用有限差分法首次对一组非线性精确至O(μ2)阶的全非线性Boussinesq方程数学模型进行了二维数值模拟分析.首先通过在方程的非线性项中引入缓坡假定,考察了其对模型数值精度的影响;其次,在模型中对二阶非线性项采用不同精度,考察了其对模型数值结果的影响.数值模拟结果表明,所建立的二阶完全非线性Boussinesq方程二维数值模型具有良好的适用性,模型非线性项中引入缓坡假定以及在二阶非线性项选用不同的精度对数值模拟结果影响不明显.  相似文献   

4.
为模拟潜堤上破碎波浪传播时产生能量的耗散这一特性,在改进的具有四阶色散的Boussinesq水波方程中中入二阶紊动粘性项,建立了考虑波浪破碎的水波数学模型.在非交错网格下建立了有限差分数值模型,并利用三阶Adams-Bash forth格式预报、四阶Adams-Mouton格式校正对数值模型进行求解.通过数值试验,模拟...  相似文献   

5.
二维物体在具有自由表面的定深水域中,作简谐摇荡产生的辐射波二阶解的传播模式,包括一个与一阶正弦进行波对应的二阶Stokes波;和另一个与二阶波数对应的二阶正弦进行波,后者的传播速度比前者慢,因此,距摇荡物体不同距离处测得的波面时间历程是不同的。在物体运动中添加一个适当的二阶量,使二阶运动的一阶响应和一阶运动产生的与二阶波数对应的二阶正弦进行波相抵消,可使二维辐射波的传播模式成为在二阶意义上纯正的Stokes波。二维辐射波二阶解中还有许多非传播模式,它们会对摇荡物体附近的辐射波波形产生比较显著的影响。  相似文献   

6.
本研究基于最优控制理论,采用变分数据同化法,通过建立伴随模型,把观测资料同化到陆架海域潮汐数值模型中去,优化开边界条件,以便提高数值预报的精度.潮汐模型的控制方程为考虑平流项、非线性底摩擦和侧向涡动粘性项的非线性浅水方程组.在第Ⅰ部分建立伴随模型和进行“孪生”数值试验的基础上,给出利用验潮站的水位资料以及TOPEX/Poseidon卫星测高数据在黄海、东海进行变分数据同化试验的数值结果.试验表明利用上述资料对模型进行变分同化校正是可行的.  相似文献   

7.
二维水波辐射问题的二阶时域解   总被引:1,自引:0,他引:1  
用二阶理论在时域范围内研究二维二阶非线性水波的辐射问题,一阶、二阶问题分别满足各自的自由表面条件和物面条件,采用人工阻尼来吸收反射波.计算采用有限元法,对水平半圆柱体在水面的振荡问题进行了计算,计算了波高时间历程和水动力,计算结果与有关文献中的数据相比符合较好.  相似文献   

8.
基于势流理论和时域高阶边界元方法,建立了三维完全非线性数值波浪水槽模型.利用源造波法产生入射波浪,应用五阶斯托克斯波理论给定波浪速度;采用混合欧拉-拉格朗日方法追踪流体瞬时水面,将二阶泰勒级数展开法应用于更新下一时间步的波面和速度势;通过加速势的方法准确计算自由水面速度的法向导数和物面速度势的时间导数.对完全非线性波浪进行了模拟,得到了稳定的波形.当波浪非线性较小时,与四阶Runge- Kutta法(RK4)计算结果和五阶斯托克斯波理论解均吻合良好;随着波浪非线性的增大,计算结果误差逐渐增大.通过数值试验分析,在满足精度要求的基础上,本方法计算时间略大于四阶Runge- Kutta法的四分之一,大大减少了计算量.  相似文献   

9.
船舶在随机波浪中操纵运动预报   总被引:1,自引:0,他引:1  
将船舶操纵运动模型和船舶耐波性理论结合起来,建立了适用于求解船舶在随机波浪中运动的六自由度数学模型.该模型采用水平随船坐标系,与操纵运动有关的水动力采用经验公式估算;与耐波性相关的水动力除二阶力外采用从频域到时域的方法,其中一阶波浪力通过脉冲响应函数的卷积求取,记忆效应力由时延函数的卷积表示;对于二阶波浪力则只考虑其定常部分.最后,应用该模型对一艘集装箱船在静水、规则波、随机波下的回转运动进行了预报,并就波浪对回转运动的影响做了讨论.认为一阶力对船舶操纵运动影响不明显;在考察波浪对船舶操纵运动的影响,必须考虑二阶力的作用.  相似文献   

10.
散射波的二阶幅射边界条件及其应用   总被引:1,自引:0,他引:1       下载免费PDF全文
黄河宁  李鉴初 《海洋学报》1988,10(2):233-239
本文借助于一个线性微分算子序列,对线性波浪绕射问题,推导出了散射波的一阶和二阶幅射边界条件,作为数值算例,在内域采用有限元,在有限距离幅射边界上分别采用一阶和二阶幅射条件及Sommerfeld条件,计算了作用在大尺度圆柱上的波浪荷载,结果表明,当采用二阶幅射条件时,精度得到很大的提高,因此,可以将内域划得很小,从而避免了以往采用Sommerfeld条件时,为获得足够的精度必须增加单元数目而导致计算量很大和占用很多计算机内存的弊病.  相似文献   

11.
The adjoint approach, one of the variational data assimilation (VDA) methods, is now widely used for fitting numerical models of meteorology and oceanography tO the observaions. The fundamental idea is to minimize a cost function, which is sum of Squares in the differences betweenthe data and their model counterparts, by adjllsting the independent model variables such as initialvalue, boundary value and Parameters. The numerical model is an operator maPPing the independent model variables i…  相似文献   

12.
We present the derivation of the discrete Euler–Lagrange equations for an inverse spectral element ocean model based on the shallow water equations. We show that the discrete Euler–Lagrange equations can be obtained from the continuous Euler–Lagrange equations by using a correct combination of the weak and the strong forms of derivatives in the Galerkin integrals, and by changing the order with which elemental assembly and mass averaging are applied in the forward and in the adjoint systems. Our derivation can be extended to obtain an adjoint for any Galerkin finite element and spectral element system.We begin the derivations using a linear wave equation in one dimension. We then apply our technique to a two-dimensional shallow water ocean model and test it on a classic double-gyre problem. The spectral element forward and adjoint ocean models can be used in a variety of inverse applications, ranging from traditional data assimilation and parameter estimation, to the less traditional model sensitivity and stability analyses, and ensemble prediction. Here the Euler–Lagrange equations are solved by an indirect representer algorithm.  相似文献   

13.
浅水潮波模式变分同化共轭码技术研究   总被引:4,自引:0,他引:4  
以浅水潮波模式为例,详细讨论了共轭码技术的使用方法以及代码检验,并建立了海洋浅水模式的共轭模式。利用浅水潮波模式及其共轭模式进行了流速和水位的初始场优化试验。试验结果表明,初始场优化对于潮波系统数值模拟具有重要的作用,同时也说明利用共轭码技术可以有效地设计共轭模式,进行各种同化试验研究,显示了共轭码技术的诸多优点。  相似文献   

14.
特征线计算格式下共轭方程两种导出途径的比较   总被引:1,自引:0,他引:1  
共轭方程的导出是建立资料同化模型的关键,其导出方式有两种途径:AFD形式与FDA形式。在特征线计算格式基础上针对一类较广泛海洋动力控制方程分析了其两种共轭方程(AFD形式与FDA形式)之间的关系,并将理论结果应用于波谱共轭方程的讨论。  相似文献   

15.
Higher order Boussinesq equations   总被引:2,自引:0,他引:2  
A new form of Boussinesq-type equations accurate to the third order are derived in this paper to improve the linear dispersion and nonlinearity characteristics in deeper water. Fourth spatial derivatives in the third order terms of the equations are transformed into second derivatives and present no difficulty in numerical computations. With the increase in accuracy of the equations, the nonlinear and dispersion characteristics of the equations are of one order of magnitude higher accuracy than those of the classical Boussinesq equations. The equations can serve as a fully nonlinear model for shallow water waves. The shoaling property of the equations is also of high accuracy through shallow water to deep water by introducing an extra source term into the second order continuity equation. An approach to increase the accuracy of the nonlinear characteristics of the new equations is introduced. The expression for the vertical distribution of the horizontal velocities is a fourth order polynomial.  相似文献   

16.
根据湍流封闭理论,建立一种适用于正压浅海湍流运动的雷诺应力封闭模型(RSM),以代替目前三维浅海动力学模型中普遍采用的湍粘性系数的传统假设。通过直接建立并模化f—平面上正压海洋的雷诺应力传输方程,分别得到的微分形式和代数形式的RSM方程组。并讨论了进行数值计算所需要的边界条件。利用该模型可以进一步研究浅海潮流、风暴潮流及风海流等浅海流动的三维结构和湍流特性。  相似文献   

17.
We adopt a parameterized internal tide dissipation term to the two-dimensional (2-D) shallow water equations, and develop the corresponding adjoint model to investigate tidal dynamics in the South China Sea (SCS). The harmonic constants derived from 63 tidal gauge stations and 24 TOPEX/Poseidon (T/P) satellite altimeter crossover points are assimilated into the adjoint model to minimize the deviations of the simulated results and observations by optimizing the bottom friction coefficient and the internal tide dissipation coefficient. Tidal constituents M2, S2, K1 and O1 are simulated simultaneously. The numerical results (assimilating only tidal gauge data) agree well with T/P data showing that the model results are reliable. The co-tidal charts of M2, S2, K1 and O1 are obtained, which reflect the characteristics of tides in the SCS. The tidal energy flux is analyzed based on numerical results. The strongest tidal energy flux appears in the Luzon Strait (LS) for both semi-diurnal and diurnal tidal constituents. The analysis of tidal energy dissipation indicates that the bottom friction dissipation occurs mainly in shallow water area, meanwhile the internal tide dissipation is mainly concentrated in the LS and the deep basin of the SCS. The tidal energetics in the LS is examined showing that the tidal energy input closely balances the tidal energy dissipation.  相似文献   

18.
一种基于Roe格式的有限体积法在二维溃坝问题中的应用   总被引:3,自引:0,他引:3  
溃坝属于典型的非恒定含间断的浅水问题。应用有限体积法离散二维浅水控制方程的守恒型方程组,将基于近似黎曼解的Roe格式用于数值计算溃坝问题,并利用MUSCL方法构造二阶空间积分格式和预测-校正二步法构造二阶时间格式,从而使数值解的整体达到二阶,提高了精度。文中算法在一维溃坝的Stoker问题的数值结果与解析解进行对比,结果证明了此方法的可行性。应用此方法在二维溃坝问题上的结果,说明了此算法可有效模拟溃坝水流的演进过程。  相似文献   

19.
A new method of assimilating sea surface height (SSH) data into ocean models is introduced and tested. Many features observable by satellite altimetry are approximated by the first baroclinic mode over much of the ocean, especially in the lower (but non-equatorial) and mid latitude regions. Based on this dynamical trait, a reduced-dynamics adjoint technique is developed and implemented with a three-dimensional model using vertical normal mode decomposition. To reduce the complexity of the variational data assimilation problem, the adjoint equations are based on a one-active-layer reduced-gravity model, which approximates the first baroclinic mode, as opposed to the full three-dimensional model equations. The reduced dimensionality of the adjoint model leads to lower computational cost than a traditional variational data assimilation algorithm. The technique is applicable to regions of the ocean where the SSH variability is dominated by the first baroclinic mode. The adjustment of the first baroclinic mode model fields dynamically transfers the SSH information to the deep ocean layers. The technique is developed in a modular fashion that can be readily implemented with many three-dimensional ocean models. For this study, the method is tested with the Navy Coastal Ocean Model (NCOM) configured to simulate the Gulf of Mexico.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号