关于K-S方程的非线性伽辽金方法(英文)

该文讨论了关于 K- S方程的伽辽金方法和非线性伽辽金方法的收敛性和 L2 误差估计 ,并得出误差阶一致的结论     

双参数弹性地基上圆底扁球壳的非线性弯曲

本文从Krmn型非线性基本微分方程出发,提出了将修正迭代法和伽辽金法联合应用,分析了Pasternak弹性地基上周边固定凹圆底扁球壳在均匀压力作用下的非线性弯曲问题,给出了荷载与挠度间的数学表达式,其所得结果与已有文献结果吻合较好,且简明、计算量小.     

大挠度圆柱壳在温度场中的热弹耦合振动分析

对温度场与与应力场耦合时的圆柱壳的非线性热弹耦合的振动问题,推导得到了基本的振动方程,热传导方程和协调方程,对短圆柱壳运用伽辽金(Galerkin)法求解,得出振幅随时间变化的数值解,得到一些有价值的结论．即随着温度幅值和耦合系数的增大,振动衰减的速度变缓,热弹耦合效应减弱．随着长径比、长厚比的增大,振幅衰减的速度变快,同时热振动频率也随之增大,即热弹耦合效应增强．耦合系数越大,轴向应力、轴向力以及轴向弯矩越小．     

小波伽辽金方法应用于变系数波动方程

我们考虑问题K（x）uxx=ua.0＜X〈1，t≥0，其中K（x）≥a≥0，u（0，t）=g，ix（0，t）=0．这是一个不适当的方程，因为当解存在时在边界g上一个小的扰动将对它的解造成很大的改变．我们考虑存在解u（x，·）∈L^2（R）用小波伽辽金方法和Meyer多分辨分析去滤掉高频部分，从而在尺度空间Vj上得到适定的近似解．我们也可以得到问题的准确解与它在Vj上的正交投影之间的误差估计．     

REFLECTION/TRANSMISSION CHARACTERISTICS OF A DISCONTINUOUS GALERKIN METHOD FOR MAXWELL＇S EQUATIONS IN DISPERSIVE INHOMOGENEOUS MEDIA

In this paper, we analyze the transmission and reflection properties of a high order discontinuous Galerkin method for dispersive Maxwell＇s equations, originally proposed by Lu et al. [J. Comput. Phys. 200 （2004）, pp. 549-580]. We study the reflection and transmission properties of the numerical method for up to second-order polynomial elements for one- and two-dimensional Maxwell＇s equations with rectangular meshes. High order accuracy has been shown for reflection and transmission coefficients near material interfaces.     

无单元Galerkin法在热传导中的应用

对热传导问题的微分方程采用无单元Galerkin法进行数值求解.首先,将微分方程用Galerkin加权残量法转化为等效的积分形式.然后,先将时间变量看作参数,对空间变量进行离散化,得到方程的半离散形式,接着,对时间采用向后Euler-Galerkin格式进行离散,得到方程的全离散形式最后,编制MATLAB程序,上机计算...     

求解m阶非线性Volterra-Fredholm型积分微分方程的一种算法

针对m阶非线性Volterra-Fredholm型积分微分方程,利用勒让德-伽辽金方法进行求解.勒让德多项式被选作基函数,通过基函数与残差正交得到有限维方程组,求解有限维方程组得到待定系数,便能求出方程的近似解.一些数值算例的给出证明了方法的可行性和有效性.     

动脉血管流动计算的伽辽金有限元法研究

得到大动脉三维模型的过二重分叉的二维截定常流的NS方程有限元解,采用了物理坐标系统换到曲线边界贴休坐标系的数学技巧,以支流至主动脉流率为参数,计算了雷诺数为1000的壁面切应力,所得结果与前人的工作（包括实验数据）进行了比较,发现与他们的结果非常接近,改进了Sharma和Kapoor(1995)的工作,相比之下,所用的数值方法上更经济,适用的雷诺数更大。     

基于单位分解积分的伽辽金无网格方法研究

数值积分是伽辽金无网格方法实施的一个重要环节,提出了一种适合于伽辽金无网格方法的单位分解积分技术．该积分技术建立在有限覆盖和单位分解基础之上,不需要对积分区域进行分解,具有较高的积分精度．并以无单元伽辽金方法为例,详细说明了基于单位分解积分的伽辽金无网格方法的实现过程．这样,在近似函数建立和数值积分过程中都不需要进行网格划分,从而形成一种“真正的”无网格方法．     

一种新的计算Timoshenko梁截面剪切系数的方法

Timoshenko梁理论中考虑了截面剪切变形的影响,推导了一种新的计算剪切系数的方法,首先采用悬臂梁纯弯曲变形条件下截面剪应力分布的精确解,基于能量原理得到了各种梁截面剪切系数新的表达式,然后推导了弯扭耦合变形条件下截面剪应力分布的精确解,进一步获得了该条件下截面的剪切系数.结果表明,悬臂梁端面作用力偏离截面的弯曲中心将使剪切系数变小,通过与Cowper计算结果的对比发现结果偏小,其原因是Cowper没有考虑与外力垂直的剪应力的影响,因此新的计算结果更优越.     

弹性波方程有限元模拟的曲边地表处理

本文对比研究了关于弹性波模拟中的曲边地表形状处理的两种方法,一种是用给定的实际介质数值划定的地表形状,另一种是用样条插值逼近地表形状.本文采用有限元方法进行弹性波数值模拟,给出了基于这两种方法计算的数值例子,并对结果进行了分析比较.结果表明使用后一种方法对地表进行处理时,地表人工离散产生的干扰明显减少,优于前一种方法.     

一种全速域的计算方法及其应用

针对原可压缩流动求解器不能用于低速不可压缩流动预测的缺点,采用预处理技术对控制方程特征系统、隐式求解方法进行修正,并采用预处理修正的AUSM+-up格式离散对流项.采用修正后的求解器对无粘鼓包流动、顶盖驱动粘性方腔流动以及Laval(拉瓦尔)喷管流动等算例进行数值仿真,并将数值仿真结果与基准解进行对比.结果表明将预处理技术应用于全速域流动问题的求解是可行的,经预处理修正后的求解器能够用于低速、亚音速、跨音速以及超音速流动问题的求解.     

定常Kuramoto-Sivashinsky方程关于非线性伽辽金方法的一个注记

本文讨论了定常Ｋ－Ｓ方程关于伽辽金方法和非线性伽辽金方法的收敛性和最大模估计;对相同模数而言,两者的误差阶完全一致,数值结果表明非线性伽辽金方法同样成功地计算出了Ｋ－Ｓ方程的分歧解,并且在计算时间方面非线性伽辽金方法比伽辽金方法要少得多。     

Secondary Buckling Analysis of Thin Rectangular Plates Based on the Wavelet Galerkin Method北大核心CSCD

Application of the wavelet Galerkin method (WGM) to numerical solution of nonlinear buckling problems was studied with classical elastic thin rectangular plates. First, the discretized scheme of the von Kármán equation were introduced, then a simple calculation approach to the Jacobian and Hessian matrices based on the WGM was proposed, and the wavelet discretized scheme-based eigenvalue equation method, the extended equation method and the pseudo arc-length method for nonlinear buckling analysis were discussed. Second, the secondary post-buckling equilibrium paths of elastic thin rectangular plates and the effects of aspect ratios, boundary conditions and bi-directional compression on the mode jumping behaviors, were discussed in detail. Numerical results show that, the WGM possesses good convergence for solving buckling loads on rectangular plates, and the obtained equilibrium paths are in good agreement with those from the stability experiments, the 2-step perturbation method and the nonlinear finite element method. Given the feasibility of combination with different bifurcation computation methods, the WGM makes an efficient spatial discretization method for complex nonlinear stability problems of typical plates and shells. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

On Fully Developed Flows of Fluids with a Pressure Dependent Viscosity in a Pipe

Stokes recognized that the viscosity of a fluid can depend on the normal stress and that in certain flows such as flows in a pipe or in channels under normal conditions, this dependence can be neglected. However, there are many other flows, which have technological significance, where the dependence of the viscosity on the pressure cannot be neglected. Numerous experimental studies have unequivocally shown that the viscosity depends on the pressure, and that this dependence can be quite strong, depending on the flow conditions. However, there have been few analytical studies that address the flows of such fluids despite their relevance to technological applications such as elastohydrodynamics. Here, we study the flow of such fluids in a pipe under sufficiently high pressures wherein the viscosity depends on the pressure, and establish an explicit exact solution for the problem. Unlike the classical Navier-Stokes solution, we find the solutions can exhibit a structure that varies all the way from a plug-like flow to a sharp profile that is essentially two intersecting lines (like a rotated V). We also show that unlike in the case of a Navier-Stokes fluid, the pressure depends both on the radial and the axial coordinates of the pipe, logarithmically in the radial coordinate and exponentially in the axial coordinate. Exact solutions such as those established in this paper serve a dual purpose, not only do they offer solutions that are transparent and provide the solution to a specific but simple boundary value problems, but they can be used also to test complex numerical schemes used to study technologically significant problems.     

A Two-Level Method for Nonsymmetric Eigenvalue Problems

A two-level discretization method for eigenvalue problems is studied. Compared to the standard Galerkin finite element discretization technique performed on a fine grid this method discretizes the eigenvalue problem on a coarse grid and obtains an improved eigenvector (eigenvalue) approximation by solving only a linear problem on the fine grid (or two linear problems for the case of eigenvalue approximation of nonsymmetric problems). The improved solution has the asymptotic accuracy of the Galerkin discretization solution. The link between the method and the iterated Galerkin method is established. Error estimates for the general nonsymmetric case are derived.     

A Sequential Least Squares Method for Elliptic Equations in Non-Divergence Form

We develop a new least squares method for solving the second-order elliptic equations in non-divergence form. Two least-squares-type functionals are proposed for solving the equation in two sequential steps. We first obtain a numerical approximation to the gradient in a piecewise irrotational polynomial space. Then together with the numerical gradient, we seek a numerical solution of the primitive variable in the continuous Lagrange finite element space. The variational setting naturally provides an a posteriori error which can be used in an adaptive refinement algorithm. The error estimates under the $L^2$ norm and the energy norm for both two unknowns are derived. By a series of numerical experiments, we verify the convergence rates and show the efficiency of the adaptive algorithm.     

A Galerkin Method for a Gaseous Ignition Model

Our aim is to investigate the existence of solutions to some second order difference equations of monotone type. Theorems 1.1 and 1.2 below are the discrete versions of some existence results due to V. Barbu [1] corresponding to the continuous case.     

A posteriori estimates for the Bubble Stabilized Discontinuous Galerkin Method

In this paper, two reliable and efficient a posteriori error estimators for the Bubble Stabilized Discontinuous Galerkin (BSDG) method for diffusion-reaction problems in two and three dimensions are derived. The theory is followed by some numerical illustrations.