共查询到19条相似文献,搜索用时 78 毫秒
1.
2.
热传导方程的一类无网格方法 总被引:1,自引:0,他引:1
构造求解热传导方程的一类无网格方法,只要选择好每个节点的适当的邻点集合,便可利用节点信息顺利进行计算.作为特殊情形,也可在各种结构或非结构网格上进行计算.在矩形或均匀平行四边形网格上进行计算时具有二阶精度,当在任意的不规则四边形或三角形网格上计算时仍然是守恒的和相容的,且至少具有一阶精度.作为数值试验,将该方法用于在不规则四边形网格上及四边形与三角形混合网格上求解二维非线性抛物型方程,并在不规则四边形网格上求解二维三温辐射热传导方程,均获得了较为理想的数值结果. 相似文献
3.
4.
对半圆域上的自然对流采用控制体有限元方法进行数据模拟,采用非结构三角形网格与质量加权迎风格式对计算域及控制方法进行离散,对压力离散方程采用波前法求解,保证了对所有Ra≤10^7时迭代过程的收敛性。 相似文献
5.
讨论了用分层求解的思想来快速构造二维数值网格。把以往文献中针对于外场区域的解抛物型方程的方法推广到封闭的二维区域,并把分层求解的思想应用于更为复杂的用椭圆型方程构造网格的数值求解中,给出了一些闭域的网格构造图。 相似文献
6.
7.
本文将非均匀网格直接离散的高阶紧致格式从二维推广到三维,结合附加修正多重网格方法提高了传统迭代方法的收敛效率,并且验证了该格式在不同边界条件的数值表现。结果表明:该方法可以有效的求解NS方程中的压力泊松方程. 相似文献
8.
9.
研究了将无覆盖区域分裂与精确控制法相结合、用于气动外形表面雷达吸波涂层问题电磁散射场的数值求解方法.无覆盖区域分裂是按空气和涂层两种不同介质物理分区要求进行的,子域交界处离散网格彼此之间可以是非匹配的,以便处理不同介质的不同网格尺度要求.电磁波通过子域交界面的约束条件是通过引入Lagrange乘子而弱满足.整个求解过程包含基于精确控制的外迭代和基于Lagrange乘子控制的区域分裂内迭代.最后计算给出了二维翼型气动外形涂层问题的电磁散射场,并与对应的无涂层外形的计算结果进行了比较,反映出表面涂层对电磁散射场的影响. 相似文献
10.
11.
提出一种数值模拟凝聚炸药爆轰问题的单元中心型Lagrange方法.利用有限体积离散爆轰反应流动方程组,基于双曲型偏微分方程组的特征理论获得离散网格节点的速度与压力,获得的网格节点速度与压力用于更新网格节点位置以及计算网格单元边的数值通量.以这种方式获得的网格节点解是一种"真正多维"的理论解,是一维Godunov格式在二维Riemann问题的推广.有限体积离散得到的爆轰反应流动的半离散系统使用一种显-隐Runge-Kutta格式来离散求解:显式格式处理对流项,隐式格式处理化学反应刚性源项.算例表明,提出的单元中心型Lagrange方法能够较好地模拟凝聚炸药的爆轰反应流动. 相似文献
12.
等离子体中Fokker-Planck方程的有限元模拟 总被引:1,自引:0,他引:1
利用有限元方法设计了一套相对简单明了的求解Fokker-Planck方程的方案.这个方案不必严格限制计算格点的步长和时间步长,就可以确保分布函数的非负性和粒子数的守恒.通过一维程序模拟,进一步证实了这个方案的可靠性.对于多维问题的分析和一维问题完全一样,所以非常容易将其推广到多维问题. 相似文献
13.
大时间步长叠波格式最初思想为LeVeque提出的大时间步长Godunov格式,通过叠加间断分解发出的强波来构造数值格式.原方法只给出了间断强波的穿越叠加方法,文章对其进行了完善,并推广到多维.针对膨胀波提出了一种网格单元分解法可以自动满足熵条件,避免出现非物理解.给出了格式的具体计算公式,并用单个守恒律方程、一维/多维Euler方程组进行了数值计算.计算结果表明,新格式除了可以采用大时间步长的优点外,在一定范围内随CFL数增加其耗散反而更低,因而对激波接触间断膨胀波的分辨率更高. 相似文献
14.
K.-M. Shyue 《Journal of computational physics》2010,229(23):8780-8801
We describe a simple mapped-grid approach for the efficient numerical simulation of compressible multiphase flow in general multi-dimensional geometries. The algorithm uses a curvilinear coordinate formulation of the equations that is derived for the Euler equations with the stiffened gas equation of state to ensure the correct fluid mixing when approximating the equations numerically with material interfaces. A γ-based and a α-based model have been described that is an easy extension of the Cartesian coordinates counterpart devised previously by the author [30]. A standard high-resolution mapped grid method in wave-propagation form is employed to solve the proposed multiphase models, giving the natural generalization of the previous one from single-phase to multiphase flow problems. We validate our algorithm by performing numerical tests in two and three dimensions that show second order accurate results for smooth flow problems and also free of spurious oscillations in the pressure for problems with interfaces. This includes also some tests where our quadrilateral-grid results in two dimensions are in direct comparisons with those obtained using a wave-propagation based Cartesian grid embedded boundary method. 相似文献
15.
K.R. Arun M. Kraft M. Luká?ová-Medvid’ová Phoolan Prasad 《Journal of computational physics》2009,228(2):565-590
We present a generalization of the finite volume evolution Galerkin scheme [M. Luká?ová-Medvid’ová, J. Saibertov’a, G. Warnecke, Finite volume evolution Galerkin methods for nonlinear hyperbolic systems, J. Comp. Phys. (2002) 183 533– 562; M. Luká?ová-Medvid’ová, K.W. Morton, G. Warnecke, Finite volume evolution Galerkin (FVEG) methods for hyperbolic problems, SIAM J. Sci. Comput. (2004) 26 1–30] for hyperbolic systems with spatially varying flux functions. Our goal is to develop a genuinely multi-dimensional numerical scheme for wave propagation problems in a heterogeneous media. We illustrate our methodology for acoustic waves in a heterogeneous medium but the results can be generalized to more complex systems. The finite volume evolution Galerkin (FVEG) method is a predictor–corrector method combining the finite volume corrector step with the evolutionary predictor step. In order to evolve fluxes along the cell interfaces we use multi-dimensional approximate evolution operator. The latter is constructed using the theory of bicharacteristics under the assumption of spatially dependent wave speeds. To approximate heterogeneous medium a staggered grid approach is used. Several numerical experiments for wave propagation with continuous as well as discontinuous wave speeds confirm the robustness and reliability of the new FVEG scheme. 相似文献
16.
Singular terms in differential equations pose severe challenges for numerical approximations on regular grids. Regularization of the singularities is a very useful technique for their representation on the grid. We analyze such techniques for the practically preferred case of narrow support of the regularizations, extending our earlier results for wider support. The analysis also generalizes existing theory for one dimensional problems to multi-dimensions. New high order multi-dimensional techniques for differential equations and numerical quadrature are introduced based on the analysis and numerical results are presented. We also show that the common use of distance functions in level-set methods to extend one dimensional regularization to higher dimensions may produce O(1) errors. 相似文献
17.
针对双曲型守恒律方程问题,发展一种有效的自适应多分辨分析方法.通过对嵌套网格上的数值解构造离散多分辨分析,建立小波系数与多层嵌套网格点之间的对应关系.对于小波系数较大的网格点采用高精度WENO格式计算,其余区域则直接采用多项式插值.数值试验表明,该方法在保持原规则网格方法的精度和分辨率的同时,显著地减少计算的CPU时间. 相似文献
18.
Jürgen Steinacker Remi HackertAdriane Steinacker Aurore Bacmann 《Journal of Quantitative Spectroscopy & Radiative Transfer》2002,73(6):557-569
We investigate the explicit numerical solution strategies of multi-dimensional radiative transfer equations which are commonly used, e.g., to determine the radiation emerging from astrophysical objects surrounded by absorbing and scattering matter. For explicit grid solvers, we identify numerical diffusion as a severe source of error in first-order discretization schemes, underestimated in former work about radiative transfer. Using the simple example of a beam propagating through vacuum, we illustrate the influence of the diffusion on the solution and discuss various techniques to reduce it. In view of the large required storage for implicit solvers, we propose to use second-order explicit grid techniques to solve 3D radiative transfer problems. 相似文献
19.
While the approximate solutions of one-dimensional nonlinear Volterra-Fredholm integral equations with smooth kermels are now well understood,no systematic studies of the numerical solutions of their multi-dimensional counterparts exist.In this paper,we provide an efficient numerical approach for the multi-dimensional nonlinear Volterra-Fredholm integral equations based on the multi-variate Legendre-collocation approach.Spectral collocation methods for multi-dimensional nonlinear integral equations are known to cause major difficulties from a convergence analysis point of view.Consequently,rigorous error estimates are provided in the weighted Sobolev space showing the exponential decay of the numerical errors.The existence and uniqueness of the numerical solution are established.Numerical experiments are provided to support the theoretical convergence analysis.The results indicate that our spectral collocation method is more flexible with better accuracy than the existing ones. 相似文献