共查询到19条相似文献,搜索用时 78 毫秒
1.
《中国海洋大学学报(自然科学版)》2017,(Z1)
近年来,复非对称代数Riccati方程越来越受到专家和学者的关注,本文主要考虑复非对称代数Riccati方程的求解问题。本文首先讨论了复非对称代数Riccati方程极值解的存在唯一性,证得了在Q的比较矩阵Q是非奇异M矩阵的假设条件下,复非对称代数Riccati方程存在唯一的极值解,并且该解可由牛顿法和不动点迭代法求得,最后通过数值实验表明这两种算法是有效可行的。 相似文献
2.
分步分层拟三维水流数学模型及其在廉州湾潮流计算中的应用 总被引:5,自引:0,他引:5
以Navier-Stokes方程为基础,经代数坐标变换导出模式方程,然后在该坐标系下将三维流场进行垂直分层,并在每一层内积分模式方程。在数值模拟求解方面,采用分步、特征差分近似与Galerkin集中质量有限元相结合的方法进行求解。在实际应用方面,对广西廉州湾潮流进行了三维模拟计算。 相似文献
3.
在曲线坐标系下,建立了缓变水深水域波浪传播的数值模拟模型.模型适宜于复杂变化的边界形状,克服了各种代数坐标变换的局限性.在建立模型时,将原始的椭圆型缓坡方程的近似型式——依赖时间变化的抛物型方程,作为控制方程,既克服了一般抛物近似方法的缺点,又便利了方程的求解;从开边界条件、不同反射特性的固壁边界条件相统一的表达式出发,对边界条件进行处理;用ADI法数值求解控制方程.对模型的验证表明,数值解与物模实验值吻合良好,模型对于具有复杂边界的工程实际有较强的适应性. 相似文献
4.
5.
6.
N-S方程的数值解法及其在水波动力学中应用的综述 总被引:7,自引:0,他引:7
对于描述流体运动的基本方程N-S方程,现已开发出了多种不同的数值求解方法。其中以MAC法和VOF法为代表的网格数值方法较为成熟,并被逐渐用于实际工程的水动力研究中:以SPH和MPS为代表的粒子方法则刚刚起步.具有很大的发展空间。针对N-S方程中的能量耗散问题,采用雷诺时均化的方程形式仍是目前解决水动力问题的主要途径,但需要引入相应湍流模型,以封闭方程。本文首先对N-S方程求解方法的发展过程进行简要回顾,对一些重要的计算方法进行评述,简要介绍其在水波动力学中的应用;然后对时均化的N-S方程的主要封闭模式进行总结介绍;最后对数值求解N-S方程的发展趋势进行展望。 相似文献
7.
非线性波浪时域计算的三维耦合模型 总被引:3,自引:1,他引:2
将计算区域Ω划分为内域Ω1和外域Ω2(Ω2=Ω-Ω1),外域控制方程采用改进线性频散特性的二维Boussinesq方程,用预报一校正法数值求解;结构物附近的内域控制方程为三维Navier-Stokes方程,由VOF方法数值求解。通过在外域和内域相匹配的交界面上设置合适的速度和波面边界条件,建立了三维非线性波浪时域计算的耦合模型。模拟试验表明:(1)耦合模型数值波浪水池可以产生稳定的、重复性较好的波动过程;(2)用耦合模型数值波浪水池求解较大浅水区域上的非线性波浪数值计算问题可以取得较高的计算效率,同时又能得出结构物附近的复杂流场。 相似文献
8.
9.
研究并行算法解决应用并行计算机完成规模尽可能大的偏微分方程的数值求解问题。利用Hopf-Cole变换,将一维非线性Burgers方程转化为线性扩散方程,基于第二类Saul’yev型非对称格式和Crank-Nicolson格式对扩散方程进行差分离散,建立解Burgers方程的交替分段并行差分格式,并讨论该方法的稳定性,给出了数值算例。此算法把剖分节点分成若干组,在每组上构造能够独立求解的差分方程,因此具有并行本性,适合在高性能多处理器的并行计算机上使用。数值试验的结果表明此方法是有效的,且有较高的精度。 相似文献
10.
研究在持续外界扰动作用下,具有控制时滞线性系统的动态输出反馈扰动抑制问题。首先利用模型转换将控制时滞系统转化为形式上无时滞的系统,通过求解Riccati方程和Sylvester方程,推导出前馈—反馈最优扰动抑制控制律。然后构造能同时预估状态和扰动的降维观测器,不仅解决前馈控制和状态反馈的物理不可实现问题,而且得到了近似于最优扰动抑制控制律的动态输出反馈扰动抑制控制器。仿真实例证明此控制律的有效性。 相似文献
11.
An exact analytical method is described to solve the diffraction problem of a group of truncated vertical cylinders. In order to account for the interaction between the cylinders, Kagemoto and Yue's exact algebraic method is utilised. The isolated cylinder diffraction potential due to incident waves is obtained using Garret's solution and evanescent mode solutions are derived in a similar manner.Numerical results are presented for arrays of two and four cylinders. Comparisons between the results obtained from the method presented here and those obtained from numerical methods show excellent agreement. 相似文献
12.
A finite differences formulation for the linear and nonlinear dynamics of 2D catenary risers 总被引:1,自引:0,他引:1
A finite differences (FD) solution method is proposed for the numerical treatment of the dynamic equilibrium problem of 2D catenary risers. The method is based on the so-called Box approximation, which in the scope of the present contribution is applied to the complete nonlinear model as well as to the reduced linearized formulation. The application of the Box method transforms the original governing systems into convenient sets of algebraic equations, which in turn are solved efficiently by the relaxation method. Extensive numerical calculations are presented that describe the dynamic behaviour of the structure and evaluate the amplification in loading due to the dynamic components. The effect of the geometric nonlinearities is assessed through comparative calculations that concern both mathematical formulations examined in the present, i.e. the complete nonlinear, and the reduced linearized model. Special attention is paid to the heave excitations as they amplify significantly the magnitudes of the loading components. 相似文献
13.
We study the problem of determination of the sound field of a point harmonic source in the coastal zone and the influence of a cylindrical body floating above the source on the sound field formed in the marine medium. A numerical-analytic method is proposed for the determination of the velocity potential. According to this method, the unknown coefficients in the general solution of the problem are determined from the corresponding infinite system of linear algebraic equations by the method of reduction. We present results of numerical calculations for a special case of a waveguide whose parameters are typical of the coastal part of the sea and perform the comparative analysis of the data obtained as a result of variation of the indicated parameters. 相似文献
14.
By using the hydroacoustic field generated by a point source in a nonuniform marine medium with a cylindrical body floating
on the surface and determined as a solution of the corresponding problem by the method of partial domains, we study the energy
characteristics of the near and far fields. The asymptotics capable of improving the accuracy of calculations is proposed
for the analysis of the posed problem for an infinite system of linear algebraic equations. We also present some results of
numerical calculations for some special cases and their comparative analysis in the case of variation of the parameters of
a waveguide within the ranges typical of the coastal zone of the sea.
__________
Translated from Morskoi Gidrofizicheskii Zhurnal, No. 3, pp. 52–59, May–June, 2006. 相似文献
15.
To examine the circulation in shallow water with tidal flat, a finite element model for the numerical solution of the shallow water equations was developed by means of standard Galerkin's method. The domain computed was covered with triangular finite elements, and water elevation and velocity were approximated by linear interpolation functions, and the lumped coefficients were used to substitute for solving the high order algebraic equation system. The time-dependent land-water boundary changes are treated mathematically by interrelating the location of the land-water boundary with the instantaneous tidal level. The implicit scheme was adopted for the terms of the bottom friction and the Coriolis effect in the motion equation so that the numerical stability of the model has been improved.The model was applied to the tidal current on shoaling water with large tidal flat off Pikou, and a comparison between observed and calculated values showed good agreement, the flow pattern being reproduced. The result 相似文献
16.
P.D. Scarlatos 《Coastal Engineering》1982,6(1):27-45
The Saint-Venant system of partial differential equations is solved by a pure finite-element method, in which integrations in both space and time are performed by utilizing Galerkin's procedure. With a special treatment of the non-linear terms, the problem is finally reduced to a linear system of algebraic equations that is solved by the conjugate gradient algorithm. This implicit scheme is proved, by numerical experiments, to be unconditionally stable. The reliability of the method is investigated by comparison of the numerical results with experimental data. Also the accuracy of the model is tested against analytical solutions for simplified cases of the unsteady free surface flow equations. 相似文献
17.
State-dependent Riccati equation-based robust dive plane control of AUV with control constraints 总被引:1,自引:0,他引:1
The paper treats the question of suboptimal dive plane control of autonomous underwater vehicles (AUVs) using the state-dependent Riccati equation (SDRE) technique. The SDRE method provides an effective mean of designing nonlinear control systems for minimum as well as nonminimum phase AUV models. It is assumed that the hydrodynamic parameters of the nonlinear vehicle model are imprecisely known, and in order to obtain a practical design, a hard constraint on control fin deflection is imposed. The problem of depth control is treated as a robust nonlinear output (depth) regulation problem with constant disturbance and reference exogenous signals. As such an internal model of first-order fed by the tracking error is constructed. A quadratic performance index is chosen for optimization and the algebraic Riccati equation is solved to obtain a suboptimal control law for the model with unconstrained input. For the design of model with fin angle constraints, a slack variable is introduced to transform the constrained control input problem into an unconstrained problem, and a suboptimal control law is designed for the augmented system using a modified performance index. Using the center manifold theorem, it is shown that in the closed-loop system, the system trajectories are regulated to a manifold (called output zeroing manifold) on which the depth tracking error is zero and the equilibrium state is asymptotically stable. Simulation results are presented which show that effective depth control is accomplished in spite of the uncertainties in the system parameters and control fin deflection constraints. 相似文献
18.
The problem of propagation of obliquely incident surface water waves over a pair of asymmetrical rectangular trenches in a channel of finite depth is examined for its solution. The mathematical problem is handled for its numerical solution with the aid of a system of singular integral equations of first kind. The resulting integral equations are solved numerically by using suitably designed polynomial approximations of the unknown functions. System of linear algebraic equations is obtained by utilizing the zeros of Chebyshev polynomial of first kind as collocation points and hence the occurrence of ill-conditioned matrix is avoided. The effectiveness of the pair of trenches is studied by analyzing the physical quantities namely the reflection and transmission coefficients. As a special limiting case, the results for a single trench problem are derived and are found to be in excellent agreement with the results available in the literature. The effect of various parameters are analyzed through different graphs for a pair of asymmetrical trenches. The energy balance relation is derived and used to check the accuracy of numerical results. 相似文献
19.
A. Hayir B. Seseogullari İ. Kilinc A. Erturk H.K. Cigizoglu M.S. Kabdasli O. Yagci K. Day 《Coastal Engineering》2008
In this paper the tsunamis resulting from a submarine mass failure such as slides and slumps triggered by earthquakes or other environmental effects, which is settled at the bottom of the north eastern Sea of Marmara are examined in one sample region. As the solution method, one hybrid method is developed. The main objective of this method is to combine an analytical solution presenting near-field tsunami amplitudes above the submarine mass failure with a numerical solution indicating the tsunami amplitudes in the coastal regions. For this purpose, one common linear boundary between analytical and numerical solution domains is defined. Movements of Submarine Mass Failures (SMF) are modeled using one simple kinematics source model and the amplitudes of the tsunamis at the region that are closer to the landslide are computed by using the analytical method. SMF is modeled approximately from the bottom geometry, and an average depth is used. Scenarios of SMF are established depending on the velocities and thicknesses of the failure, and near-field tsunami amplitudes are obtained in the open sea during the source time. After the source times, the solutions are found in the numerical region using TELEMAC-2D software system with the mentioned boundary above. In this boundary, the output of the analytical solutions is taken as the boundary conditions or the disturbances for the numerical method. With these disturbances, the numerical method is performed and the amplitudes are calculated in the coastal area. The generation, propagation and coastal amplifications of the tsunamis are illustrated at some certain points and regions both in the open sea and near the coast line. The results have been visualized and discussed. 相似文献