Galerkin-Bernstein Approximations for the System of Third-Order Nonlinear Boundary Value Problems

This paper is devoted to find the numerical solutions of one dimensional general nonlinear system of third-order boundary value problems (BVPs) for the pair of functions using Galerkin weighted residual method. We derive mathematical formulations in matrix form, in detail, by exploiting Bernstein polynomials as basis functions. A reasonable accuracy is found when the proposed method is used on few examples. At the end of the study, a comparison is made between the approximate and exact solutions, and also with the solutions of the existing methods. Our results converge monotonically to the exact solutions. In addition, we show that the derived formulations may be applicable by reducing higher order complicated BVP into a lower order system of BVPs, and the performance of the numerical solutions is satisfactory. .

采用重组模板的权重优化WENO-Z格式

针对精确模拟含激波等复杂流动结构的流场对高精度格式的低耗散低色散要求,基于5阶有限差分WENO-Z格式,提出一种模板重组技术。在计算WENO非线性权时,引入一个由3点模板重新组合的4点模板,优化原格式中各模板的权重分配,进而提出了两种改进WENO-Z格式。采用近似色散关系分析方法对改进前后格式色散与耗散特性进行了对比与分析。分析表明:两种改进格式耗散有不同程度的降低。数值实验表明:改进格式具有更优越的激波捕捉性能,对小尺度流场结构具有更高的分辨率。

求解双曲守恒律方程的高阶WENO型熵相容格式的对比研究

本文采用五阶WENO重构方法,结合熵相容格式,得到求解双曲守恒律方程的高阶WENO型熵相容格式.数值实验证明了所得格式的鲁棒性.

求解双曲守恒律的修正模板近似的五阶WENO格式

针对经典的五阶加权本质无振荡(WENO)格式在间断附近耗散过大以及临界点不能保精度的问题,本文提出了一种新的修正模板近似方法。改进了经典五阶WENO-JS格式中各候选子模板上数值通量的二阶多项式逼近,通过加入三次修正项使模板逼近达到四阶精度,并且通过引入可调函数φ使得新的格式具有ENO性质,理论分析新的格式具有保精度特性,通过一系列数值算例说明了新格式的高效性。

Existence and Uniqueness Results for a Fully Third-Order Boundary Value Problem

The boundary value problems of the third-order ordinary differential equation have many practical application backgrounds and their some special cases have been studied by many authors. However, few scholars have studied the boundary value problems of the complete third-order differential equations u′′′(t) = f (t,u(t),u′(t),u′′(t)). In this paper, we discuss the existence and uniqueness of solutions and positive solutions of the fully third-order ordinary differential equation on [0,1] with the boundary condition u(0) = u′(1) = u′′(1) = 0. Under some inequality conditions on nonlinearity f some new existence and uniqueness results of solutions and positive solutions are obtained.

基于混合模板改进的五阶WENO格式在求解双曲守恒律方程中的应用

本文在WENO-MZQ3格式的基础上,将其模板推广到五阶WENO,即在经典五阶WENO格式的三个三点模板的基础上引入五点模板,并引入新的改进的全局光滑因子,对其进行加权处理得到新的五阶WENO格式.理论分析表明新的五阶WENO格式在光滑区域比五阶WENO-JS格式具有更高阶的计算精度,通过一系列一维和二维数值算例也验证了新的五阶WENO格式具有更高的精度和分辨率.

一种基于欧拉方程通量分裂的五阶有限差分共权多分辨率WENO格式

本文研究了欧拉方程的高精度数值解法,在多分辨率WENO数值格式中引入了共权思想,获得了一种新的五阶有限差分共权多分辨率WENO格式.数值实验表明,该方法在解的光滑区域能获得一致的高精度,在强间断附近能保持数值解基本无振荡的性质,从而验证算法的有效性.

非均匀结构网格上MUSCL和WENO格式的精度

基于一维均匀网格条件下构造的差分格式,在实际应用中须推广到非均匀或者曲线网格上,坐标变换过程引入几何诱导误差。目前常用收敛解误差随着网格细化变化的精度测试方法评估差分格式的精度。在二维柱坐标均匀网格上,采用1阶迎风、2阶MUSCL和5阶WENO计算流场参数为常数的自由流问题,按照精度测试方法比较收敛曲线斜率,发现1阶迎风的网格收敛精度是2阶的,5阶WENO的网格收敛精度不到1阶。理论分析表明,这种精度测试方法与差分格式精度定义不等价,而且所采用的数据无法反映差分格式的固有缺陷,因此,不能用来作为差分格式精度评价指标。很多研究WENO的文献经常模拟双Mach反射问题、二维Riemann问题等经典算例,把接触间断是否演变成不稳定涡结构作为特征,理论上可以证明涡结构是非物理现象,因此用是否出现涡结构作为算法高精度的论据并不合适。

Hypersonic Shock Wave/Boundary Layer Interactions by a Third-Order Optimized Symmetric WENO Scheme

A novel third-order optimized symmetric weighted essentially non-oscillatory(WENO-OS3)scheme is used to simulate the hypersonic shock wave/boundary layer interactions.Firstly,the scheme is presented with the achievement of low dissipation in smooth region and robust shock-capturing capabilities in discontinuities.The Maxwell slip boundary conditions are employed to consider the rarefied effect near the surface.Secondly,several validating tests are given to show the good resolution of the WENO-OS3 scheme and the feasibility of the Maxwell slip boundary conditions.Finally,hypersonic flows around the hollow cylinder truncated flare(HCTF)and the25°/55°sharp double cone are studied.Discussions are made on the characteristics of the hypersonic shock wave/boundary layer interactions with and without the consideration of the slip effect.The results indicate that the scheme has a good capability in predicting heat transfer with a high resolution for describing fluid structures.With the slip boundary conditions,the separation region at the corner is smaller and the prediction is more accurate than that with no-slip boundary conditions.

Efficient Finite Difference WENO Scheme for Hyperbolic Systems withNon-conservativeProducts

Higher order finite difference weighted essentially non-oscillatory(WENO)schemes have been constructed for conservation laws.For multidimensional problems,they offer a high order accuracy at a fraction of the cost of a finite volume WENO or DG scheme of the comparable accuracy.This makes them quite attractive for several science and engineering applications.But,to the best of our knowledge,such schemes have not been extended to non-linear hyperbolic systems with non-conservative products.In this paper,we perform such an extension which improves the domain of the applicability of such schemes.The extension is carried out by writing the scheme in fluctuation form.We use the HLLI Riemann solver of Dumbser and Balsara(J.Comput.Phys.304:275-319,2016)as a building block for carrying out this extension.Because of the use of an HLL building block,the resulting scheme has a proper supersonic limit.The use of anti-diffusive fluxes ensures that stationary discontinuities can be preserved by the scheme,thus expanding its domain of the applicability.Our new finite difference WENO formulation uses the same WENO reconstruction that was used in classical versions,making it very easy for users to transition over to the present formulation.For conservation laws,the new finite difference WENO is shown to perform as well as the classical version of finite difference WENO,with two major advantages:(i)It can capture jumps in stationary linearly degenerate wave families exactly.(i)It only requires the reconstruction to be applied once.Several examples from hyperbolic PDE systems with non-conservative products are shown which indicate that the scheme works and achieves its design order of the accuracy for smooth multidimensional flows.Stringent Riemann problems and several novel multidimensional problems that are drawn from compressible Baer-Nunziato multiphase flow,multiphase debris flow and twolayer shallow water equations are also shown to document the robustness of the method.For some test problems that require well-balancing we have even been able to apply the scheme without any modification and obtain good results.Many useful PDEs may have stiff relaxation source terms for which the finite difference formulation of WENO is shown to provide some genuine advantages.

Sparse-Grid Implementation of Fixed-Point Fast Sweeping WENO Schemes for Eikonal Equations

Fixed-point fast sweeping methods are a class of explicit iterative methods developed in the literature to efficiently solve steady-state solutions of hyperbolic partial differential equations(PDEs).As other types of fast sweeping schemes,fixed-point fast sweeping methods use the Gauss-Seidel iterations and alternating sweeping strategy to cover characteristics of hyperbolic PDEs in a certain direction simultaneously in each sweeping order.The resulting iterative schemes have a fast convergence rate to steady-state solutions.Moreover,an advantage of fixed-point fast sweeping methods over other types of fast sweeping methods is that they are explicit and do not involve the inverse operation of any nonlinear local system.Hence,they are robust and flexible,and have been combined with high-order accurate weighted essentially non-oscillatory(WENO)schemes to solve various hyperbolic PDEs in the literature.For multidimensional nonlinear problems,high-order fixed-point fast sweeping WENO methods still require quite a large amount of computational costs.In this technical note,we apply sparse-grid techniques,an effective approximation tool for multidimensional problems,to fixed-point fast sweeping WENO methods for reducing their computational costs.Here,we focus on fixed-point fast sweeping WENO schemes with third-order accuracy(Zhang et al.2006[41]),for solving Eikonal equations,an important class of static Hamilton-Jacobi(H-J)equations.Numerical experiments on solving multidimensional Eikonal equations and a more general static H-J equation are performed to show that the sparse-grid computations of the fixed-point fast sweeping WENO schemes achieve large savings of CPU times on refined meshes,and at the same time maintain comparable accuracy and resolution with those on corresponding regular single grids.

7阶WENO-S格式的计算效率研究

WENO-S格式是一类适合于含间断问题数值模拟的加权本质无振荡格式.这类格式的光滑因子满足对单频波为常数,这使得其近似色散关系与线性基底格式一致,并且具有良好的小尺度波动模拟能力.计算效率是数值方法性能指标的一个重要方面.由于WENO-S格式的光滑因子在各子模板上的计算公式除下标不同外形式一致,在计算线性对流方程相邻数值通量时,部分光滑因子完全相同.为此提出一种消除WENO-S格式冗余光滑因子计算的方法.该方法要求一条网格线上用于重构或插值的量可以表示为一个序列.基于此要求分析其对于几种不同物理问题的可行性和使用方法.以7阶WENO-S格式为例介绍了格式性质和去除冗余光滑因子计算的方法.该方法中预先计算和存储一条网格线上的所有光滑因子,在网格点较多的情况下,光滑因子计算次数约为原7阶WENO-S格式的1/4.对一维对流问题、球面波传播问题、二维旋转问题、二维小扰动传播问题及一维和二维无黏流动问题进行了数值模拟.结果表明该格式对多种流动结构具有良好的捕捉能力,并且同时具有良好的计算效率,去除冗余计算后又降低了约20%的计算时间.

MULTIPLE POSITIVE SOLUTIONS OF SINGULAR THIRD-ORDER PERIODIC BOUNDARY VALUE PROBLEM

This paper deals with the singular nonlinear third-order periodic boundary value problem u'' + p(3)u = f (t, u), 0 less than or equal to t less than or equal to 2pi, with u((i)) (0) = u((i)) (2pi), i = 0, 1, 2, where p is an element of (Graphics) and f is singular at t = 0, t = 1 and u = 0. Under suitable weaker conditions than those of [1], it is proved by constructing a special cone in C[0, 2pi] and employing the fixed point index theory that the problem has at least one or at least two positive solutions.

求解理想磁流体方程的四阶WENO型熵稳定格式

构造了一种用于求解理想磁流体方程的四阶熵稳定半离散有限体积格式.该格式空间方向上将高阶熵守恒通量与采用WENO重构的耗散项结合,得到高阶熵稳定通量.通过在耗散项中添加开关函数,使得数值通量具有更低的耗散并且高阶WENO重构满足符号性质.对用来控制磁场散度的源项采用中心格式离散,最终得到与熵守恒通量一致的高阶精度.几个一维、二维算例表明该格式无振荡,鲁棒性强,可以精确捕捉间断.

WENO格式的一种新型光滑因子及其应用

光滑因子是构造WENO格式的关键,决定了WENO格式能否达到最优收敛精度以及在间断附近能否保持本质无振荡特性。针对广泛应用的五阶WENO格式,采用算子函数近似导数关系的方法,设计了一种新型光滑因子。在间断区域,与传统的光滑因子相比,新型光滑因子对间断的识别更准确。在光滑区域,新型光滑因子使得新型WENO格式的非线性权更接近线性权。理论分析和数值验证表明新型WENO格式即使在一阶临界点也能保持一致五阶精度。一维典型激波问题的数值结果表明,与现有WENO格式相比,新型WENO格式提高了短波的分辨率,降低了格式的耗散,同时能够准确识别间断。对于典型包含激波和剪切层的流动问题,新型WENO格式不仅能够准确地捕捉强激波,而且能够精细模拟剪切层和声波等流场结构。

Accuracy of the half-power bandwidth method with a third-order correction for estimating damping in multi-DOF systems

A third-order correction was recently suggested to improve the accuracy of the half-power bandwidth method in estimating the damping of single DOF systems.This paper analyzes the accuracy of the half-power bandwidth method with the third-order correction in damping estimation for multi-DOF linear systems.Damping ratios in a two-DOF linear system are estimated using its displacement and acceleration frequency response curves,respectively.A wide range of important parameters that characterize the shape of these response curves are taken into account.Results show that the third-order correction may greatly improve the accuracy of the half-power bandwidth method in estimating damping in a two-DOF system.In spite of this,the half-power bandwidth method may significantly overestimate the damping ratios of two-DOF systems in some cases.

Conformal invariance and conserved quantity of third-order Lagrange equations for non-conserved mechanical systems

This paper studies conformal invariance and conserved quantity of third-order Lagrange equations for non- conserved mechanical systems. Third-order Lagrange equations, the definition and a determining equation of conformal invariance of the system are presented. The conformal factor expression is deduced from conformal invariance and Lie symmetry. The necessary and sufficient condition that conformal invaxiance of the system would have Lie symmetry under single-parameter infinitesimal transformations is obtained. The corresponding conserved quantity of conformal invariance is derived with the aid of a structure equation. Lastly, an example is given to illustrate the application of the results.

Oscillation Criteria for Third-order NonlinearNeutral Dynamic Equations on Time Scales

Many practical problems, such as those from electronic engineering, mechanicalengineering, ecological engineering, aerospace engineering and so on, need to bedescribed by dynamic equations on time scales, so it is important in theory andpractical significance to study these equations. In this paper, the oscillation andasymptotic behavior of third-order nonlinear neutral delay dynamic equations ontime scales are studied by using generalized Riccati transformation technique, integralaveraging methods and comparison theorems. The main purpose of this paperis to establish some new oscillation criteria for such dynamic equations. The newKamenev criteria and Philos criteria are given, and an example is considered toillustrate our main results.

Symmetry Reduction of Initial-Value Problems for a Class of Third-order Evolution Equations

Symmetry reduction of a class of third-order evolution equations that admit certain generalized conditionalsymmetries (GCSs) is implemented.The reducibility of the initial-value problem for an evolution equation to a Cauchyproblem for a system of ordinary differential equations (ODEs) is characterized via the GCS and its Lie symmetry.Complete classification theorems are obtained and some examples are taken to show the main reduction procedure.

On the Existence of Positive Solution for a Nonlinear Third-order Three-point Boundary Value Problem

An existence theorem of positive solution is established for a nonlinear third-order three-point boundary value problem. Here, we concentrated on the case that the nonlinear term is neither superlinear nor sublinear, and is not asymptotic at zero and infinity.