共查询到20条相似文献，搜索用时 0 毫秒
1.
为捕捉双曲守恒律方程间断,提高算法求解精度,本文应用区域压缩PINN(physics-informed neural networks)算法对双曲守恒律方程近似求解。首先,对物理方程添加速度梯度监测函数,以此识别和压缩大梯度区域;随后,针对不同初始条件的双曲守恒律方程,设定相应的大梯度区域压缩控制系数,降低其在损失函数中占的比重;最后,将带有速度梯度权重项的损失函数放入神经网络中训练,通过最小化损失函数学习方程在整个区域上的解。利用区域压缩PINN算法求解各种经典双曲守恒律问题,通过对满足不同初始条件的一维和二维双曲守恒律方程进行数值模拟,并与经典PINN算法结果进行比较,验证了区域压缩PINN算法的良好性能。 相似文献
2.
R. Bürger K. H. Karlsen E. M. Tory W. L. Wendland 《ZAMM - Journal of Applied Mathematics and Mechanics / Zeitschrift für Angewandte Mathematik und Mechanik》2002,82(10):699-722
The one‐dimensional kinematical sedimentation theory for suspensions of small spheres of equal size and density is generalized to polydisperse suspensions and several space dimensions. The resulting mathematical model, obtained by introducing constitutive assumptions and performing a dimensional analysis, is a system of first‐order conservation laws for the concentrations of the solids species coupled to a variant of the Stokes system for incompressible flow describing the mixture. Various flux density vectors for the first‐order system have been proposed in the literature. Some of them cause the first‐order system of conservation laws to be non‐hyperbolic, or to be of mixed hyperbolic‐elliptic type in the bidisperse case. The criterion for ellipticity is equivalent to a well‐known instability criterion predicting phenomena like blobs and viscous fingering in bidisperse sedimentation. We show that loss of hyperbolicity, that is the occurrence of complex eigenvalues of the Jacobian of the first‐order system, can be viewed as an instability criterion for arbitrary polydisperse suspensions, and that for tridisperse mixtures this criterion can be evaluated by a convenient calculation of a discriminant. We determine instability regions (or alternatively prove stability) for three different choices of the flux vector of the first‐order system of conservation laws. Consequently, mixed or non‐hyperbolic, rather than hyperbolic, systems of conservation laws are the appropriate general mathematical framework for polydisperse sedimentation. The stability analysis examines a first‐order system of conservation laws, but its predictions are applicable to the full multidimensional system of model equations. The findings are consistent with experimental evidence and are appropriately embedded into the current state of knowledge of non‐hyperbolic systems of conservation laws. 相似文献
3.
4.
The mathematical model describing transients in natural gas pipelines constitutes a non‐homogeneous system of non‐linear hyperbolic conservation laws. The time splitting approach is adopted to solve this non‐homogeneous hyperbolic model. At each time step, the non‐homogeneous hyperbolic model is split into a homogeneous hyperbolic model and an ODE operator. An explicit 5‐point, second‐order‐accurate total variation diminishing (TVD) scheme is formulated to solve the homogeneous system of non‐linear hyperbolic conservation laws. Special attention is given to the treatment of boundary conditions at the inlet and the outlet of the pipeline. Hybrid methods involving the Godunov scheme (TVD/Godunov scheme) or the Roe scheme (TVD/Roe scheme) or the Lax–Wendroff scheme (TVD/LW scheme) are used to achieve appropriate boundary handling strategy. A severe condition involving instantaneous closure of a downstream valve is used to test the efficacy of the new schemes. The results produced by the TVD/Roe and TVD/Godunov schemes are excellent and comparable with each other, while the TVD/LW scheme performs reasonably well. The TVD/Roe scheme is applied to simulate the transport of a fast transient in a short pipe and the propagation of a slow transient in a long transmission pipeline. For the first example, the scheme produces excellent results, which capture and maintain the integrity of the wave fronts even after a long time. For the second example, comparisons of computational results are made using different discretizing parameters. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
5.
We propose a fully conservative high‐order upwind multi‐moment method for the conservation equation. The proposed method is based on a third‐order polynomial interpolation function and semi‐Lagrangian formulation and is a variant of the constrained interpolation profile conservative semi‐Lagrangian scheme with third‐order polynomial function method. The third‐order interpolation function is constructed based on three constraints in the upwind cell (two boundary values and a cell average) and a constraint in the downwind cell (a cell center value). The proposed method shows fourth‐order accuracy in a benchmark problem (sine wave propagation). We also propose a less oscillatory formulation of the proposed method. The less oscillatory formulation can minimize numerical oscillations. These methods were validated through scalar transport problems, and compressible flow problems (shock tube and 2D explosion problems). Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
6.
A traveling wave solution to the Aw-Rascle traffic flow model that includes the relaxation and diffusion terms is investigated. The model can be approximated by the well-known Kortweg-de Vries (KdV) equation. A numerical simulation is conducted by the first-order accurate Lax-Friedrichs scheme, which is known for its ability to capture the entropy solution to hyperbolic conservation laws. Periodic boundary conditions are applied to simulate a lengthy propagation, where the profile of the derived KdV solution is taken as the initial condition to observe the change of the profile. The simulation shows good agreement between the approximated KdV solution and the numerical solution. 相似文献
7.
将Jin's的界面方法应用到求解双曲守恒型方程的半离散中心迎风方法中,给出了一种新的求解浅水波方程的半离散中心迎风差分方法。对于源项,不是采用传统的单元均值而是采用单元界面处的值来近似,使所得格式对稳定态的求解是均衡的。且已证明所给的二阶精度的求解格式保持水深的非负性,这一特性使其能够较好的处理干河床问题。使用该方法产生的数值粘性(与O(Δ2r-1)同阶)要比交错的中心格式小(与O(Δx2r/Δt)同阶),而且由于数值粘性与时间步长无关,从而时间步长可根据稳定性需要尽可能的小,因此适用于稳定态的求解。 相似文献
8.
Yu. E. Ivanova V. E. Ragozina 《Journal of Applied Mechanics and Technical Physics》2006,47(6):892-898
The method of matched asymptotic expansions was employed to obtain approximate solutions to the one-dimensional boundary-value
problems of nonlinear dynamic elasticity theory of impact loading on the surface of a cylindrical cavity of an incompressible
medium that causes antiplane motion or torsion of the medium. The expansion of the solution in the near-front region is based
on solutions of evolution equations different from the equations for quasi-simple waves.
__________
Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 6, pp. 144–151, November–December, 2006. 相似文献
9.
Summary A numerical scheme is presented which employs the characteristic surfaces in space-time for solving Navier-Stokes equations for compressible fluid flow. We consider the general case of a three-dimensional flow, a simplification of which yields the equations of the two-dimensional case. Emphasis is put on the method itself. We apply it to simulate a laminar hypersonic flow around a circular cylinder of a five-components gas mixture of nitrogen and oxygen with thermally perfect constituents and at chemical nonequilibrium. First, the partial differential equations are transformed into a standard form with directional derivatives, enabling to attain the compatibility conditions, including the viscosity terms. These conditions are discretized by approximating their integrals along the corresponding characteristic surfaces. The result is an explicit time-marching numerical scheme. Using a grid fitted between the shock and the cylinder, and starting from roughly estimated initial conditions, a steady solution is searched. A comparison is made with the solution obtained under the assumption of a perfect gas. Received 6 April 1999; accepted for publication 13 May 1999 相似文献
10.
L. J. Crane A. G. McVeigh 《ZAMM - Journal of Applied Mathematics and Mechanics / Zeitschrift für Angewandte Mathematik und Mechanik》2015,95(12):1439-1447
This work gives a mathematical account of the propagation of disturbances produced by a sudden draught of air impinging on a cylindrical fibre produced by the melt‐spinning process. Accounting for varying tension and aerodynamic drag, the non‐dimensional wave‐type equation of motion is derived in dimensionless form; the solution of which is sought using the Riemann method of characteristics. Using this procedure, the solutions along the leading forward and reflected characteristics are obtained in closed‐form and enable numerical solutions to be obtained via a finite difference routine along the entire computational domain. The analysis shows how such disturbances may be prevented from penetrating beyond the material crystallisation point (where it is extremely susceptible to disturbances) and discusses the application of these findings by optimising the location of a protective shroud. 相似文献
11.
12.
JérÔme Jaffré 《Transport in Porous Media》1995,21(3):195-207
When simulating two-phase flow in porous media, one has to consider the case where there is a discontinuity in the medium. There relative permeabilities and capillary pressure functions may change and we address the problem of calculating the convective part of the numerical flux at the interface between the two rock types. Several solutions are compared. 相似文献
13.
Sonic boom focusing phenomenon can be predicted using the solution to the nonlinear Tricomi equation which is a hybrid (hyperbolic‐elliptic) second‐order partial differential equation. In this paper, the hyperbolic conservation law form is derived, which is valid in the entire domain. In this manner, the presence of two regions where the equation behaves differently (hyperbolic in the upper and elliptic in the lower half‐plane) is avoided. On the upper boundary, a new mixed boundary condition for the acoustic pressure is employed. The discretization is carried out using a discontinuous Galerkin (DG) method combined with a Runge–Kutta total‐variation diminishing scheme. The results show the accuracy of DG methods to solve problems involving sharp gradients and discontinuities. Comparisons with analytical results for the linear case, and other numerical results using classical explicit and compact finite difference schemes and weighted essentially non‐oscillatory schemes are included. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
14.
基于波前动量守恒理论和位移不连续方法所提出的时域分析新方法,引入岩石非线性法向本构关系,对弹性纵波在岩石非线性节理中的传播特性进行了理论分析。采用节理变形的双曲线模型(BB模型),获得纵波P波斜入射非线性节理的传播波动方程,并通过参数研究分析了在岩石节理中节理非线性系数、节理初始刚度、应力波入射角和入射波幅值等因素对纵波传播规律的影响。结果表明:所推导的应力波传播方程在考虑多种非线性问题时,通过迭代计算即可方便求出透射波和反射波的数值解,避免了复杂的数学运算;当波斜入射节理面时,产生了波型转换,节理变形的非线性对透射波和反射波有较大影响,透射系数和反射系数并非随着非线性参数的变化而单调变化。时域内所推导的波传播方程更有益于波斜入射时非线性参数的广泛研究,为开展该方面的理论研究工作提供了借鉴。 相似文献
15.
Xian Wen;Xuhong Yu ;Zhongqing Wang 《Discrete and continuous dynamical systems》2023,28(6):3349-3372
An efficient and accurate Legendre dual-Petrov-Galerkin spectral element method for solving the Kawahara-type equations is proposed. The Sobolev bi-orthogonal basis functions are constructed, which reduce the number of nonzero elements in the matrix and improve the computational efficiency. The generalized stability of the fully discrete spectral element scheme is analyzed. For applications, we consider the Kawahara, KdV-Kawahara and modified Kawahara equations, and discuss numerically the motion of single solitary wave, conservation laws, the phenomena of wave generation and interaction of two and three solitons. Numerical results illustrate the effectiveness of the suggested approach. 相似文献
16.
A new third‐order WENO scheme is proposed to achieve the desired order of convergence at the critical points for scalar hyperbolic equations. A new reference smoothness indicator is introduced, which satisfies the sufficient condition on the weights for the third‐order convergence. Following the truncation error analysis, we have shown that the proposed scheme achieves the desired order accurate for smooth solutions with arbitrary number of vanishing derivatives if the parameter ε satisfies certain conditions. We have made a comparative study of the proposed scheme with the existing schemes such as WENO‐JS, WENO‐Z, and WENO‐N3 through different numerical examples. The result shows that the proposed scheme (WENO‐MN3) achieves better performance than these schemes. 相似文献
17.
In this article, we have devised a new reference smoothness indicator for third‐order weighted essentially non‐oscillatory (WENO) scheme to achieve desired order of convergence at critical points. In the context of the weighted essentially non‐oscillatory scheme, reference smoothness indicator is constructed in such a way that it satisfies the sufficient condition on the weights for the third‐order convergence. The goal is to construct a reference smoothness indicator such that the resulted scheme have to achieve the required order of accuracy even if the first two derivatives vanish but not the third derivative. The construction of such reference smoothness indicator is not possible through a linear combination of local smoothness indicators only. We have proposed a reference smoothness indicator to be of the fourth order of accuracy on three‐point stencil that contains the linear combination of the first derivative information of the local and global stencils. The performance enhancement of the WENO scheme through this reference smoothness indicator is verified through the standard numerical experiments. Numerical results indicate that the new scheme provides better results in comparison with the earlier third‐order WENO schemes like WENO‐JS and WENO‐Z. Copyright © 2017 John Wiley & Sons, Ltd. 相似文献
18.
In this article, we present two improved third‐order weighted essentially nonoscillatory (WENO) schemes for recovering their design‐order near first‐order critical points. The schemes are constructed in the framework of third‐order WENO‐Z scheme. Two new global smoothness indicators, τL3 and τL4, are devised by a nonlinear combination of local smoothness indicators (ISk) and reference values (ISG) based on Lagrangian interpolation polynomial. The performances of the proposed schemes are evaluated on several numerical tests governed by one‐dimensional linear advection equation or one‐ and two‐dimensional Euler equations. Numerical results indicate that the presented schemes provide less dissipation and higher resolution than the original WENO3‐JS and subsequent WENO3‐N scheme. 相似文献
19.
In a previous paper (Transp. Porous Media,55(1): 47–70), algorithms are given for computing the analytical solution to the three-phase Riemann problem. Application of those algorithms requires that the wave configuration is known. The purpose of this note is to provide a procedure to determine the wave structure for any initial and injected saturation states. 相似文献
20.
近几年,在计算流体力学中,高精度、高分辨率的加权本质无振荡(weighted essentially non-oscillatory , WENO)格式得到很大的发展.WENO格式的主要思想是通过低阶的数值流通量的凸组合重构得到高阶的逼近,并且在间断附近具有本质无振荡的性质.本文综合介绍了双曲守恒律方程的有限差分和有限体积迎风型WENO,中心WENO,紧致中心WENO以及优化的WENO格式等,讨论了负权的处理和多维问题的解决方法.最后,通过一些算例证明WENO格式的高精度,本质无振荡的性质.图6参40 相似文献