共查询到19条相似文献,搜索用时 46 毫秒
1.
2.
粒子输运方程的线性间断有限元方法 总被引:1,自引:0,他引:1
将空间线性间断有限元方法应用于动态粒子输运方程的求解.数值算例表明,空间线性间断有限元方法在网格边界的数值精度方面明显高于指数格式和菱形格式,并且通量在时间上的微分曲线相对光滑,避免了指数格式、菱形格式数值解的非物理振荡现象. 相似文献
3.
4.
讨论了谱有限元方法在中子测井数值模拟中的应用,用球谐函数谱展开和有限元耦合方法求解Boltzmann中子输运方程,得到了这种耦合方法的收敛性。研制了三维有限元程序,实现了中子测井问题数值模拟的正演计算,实际算例表明此方法是有效的。 相似文献
5.
6.
7.
将龙格库塔间断有限元方法(RDDG)与自适应方法相结合,求解三维欧拉方程.区域剖分采用非结构四面体网格,依据数值解的变化采用自适应技术对网格进行局部加密或粗化,减少总体网格数目,提高计算效率.给出四种自适应策略并分析不同自适应策略的优缺点.数值算例表明方法的有效性. 相似文献
8.
研究自适应Runge-Kutta间断Galerkin (RKDG)方法求解双曲守恒律方程组,并提出两种生成相容三角形网格的自适应算法.第一种算法适用于规则网格,实现简单、计算速度快.第二种算法基于非结构网格,设计一类基于间断界面的自适应网格加密策略,方法灵活高效.两种方法都具有令人满意的计算效果,而且降低了RKDG的计算量. 相似文献
9.
中子输运方程的计算量非常大。在现有的计算机条件下,进行精密物理的数值模拟所需要的中子计算仍是非常的费时间和费内存的,不采用并行计算是难以承受的。并且,由于中子输运隐式离散纵标方法引起的数据强相关,以及计算过程必须严格沿中子运动方向进行(否则会出现计算不稳定),因此会出现相当严重的算法同步的问题,使得隐式格式在大型并行计算机上实施时所能得到的并行度十分有限,严格限制了其实现具有高并行度的迭代计算的可能性。因此,对中子输运方程隐式差分格式进行并行改造是十分必要的。 相似文献
10.
高分辨率间断有限元方法 总被引:2,自引:0,他引:2
间断有限元方法是集高分辨率有限差分方法和有限体积方法的优点发展起来的一种数值方法,在计算流体动力学问题上显示了优良的效能.利用守恒问题给出间断有限元方法的基本概念和过程,利用简单算例给出该方法的精度分析和限制器对精度的影响,并给出浅水波问题、交通流问题和波传播问题的数值模拟结果,进一步,综合评介该方法在椭圆、抛物、对流扩散、Hamilton-Jacobi方程、Navier-Stokes方程等的实际应用进展. 相似文献
11.
We construct a new nonlinear finite volume scheme for diffusion equation on polygonal meshes and prove that the scheme satisfies the discrete extremum principle. Our scheme is locally conservative and has only cell-centered unknowns. Numerical results are presented to show how our scheme works for preserving discrete extremum principle and positivity on various distorted meshes. 相似文献
12.
13.
P. Surya Mohan Tanja Tarvainen Martin Schweiger Aki Pulkkinen Simon R. Arridge 《Journal of computational physics》2011,230(19):7364-7383
We propose the PN approximation based on a finite element framework for solving the radiative transport equation with optical tomography as the primary application area. The key idea is to employ a variable order spherical harmonic expansion for angular discretization based on the proximity to the source and the local scattering coefficient. The proposed scheme is shown to be computationally efficient compared to employing homogeneously high orders of expansion everywhere in the domain. In addition the numerical method is shown to accurately describe the void regions encountered in the forward modeling of real-life specimens such as infant brains. The accuracy of the method is demonstrated over three model problems where the PN approximation is compared against Monte Carlo simulations and other state-of-the-art methods. 相似文献
14.
15.
In this paper, we use finite element methods to simulate the hydrodynamical systems governing the motions of nematic liquid crystals in a bounded domain Ω. We reformulate the original model in the weak form which is consistent with the continuous dissipative energy law for the flow and director fields in W1,2+σ(Ω) (σ > 0 is an arbitrarily small number). This enables us to use convenient conformal C0 finite elements in solving the problem. Moreover, a discrete energy law is derived for a modified midpoint time discretization scheme. A fixed iterative method is used to solve the resulted nonlinear system so that a matrix free time evolution may be achieved and velocity and director variables may be solved separately. A number of hydrodynamical liquid crystal examples are computed to demonstrate the effects of the parameters and the performance of the method. 相似文献
16.
Bou Matar O Guerder PY Li Y Vandewoestyne B Van Den Abeele K 《The Journal of the Acoustical Society of America》2012,131(5):3650-3663
A nodal discontinuous Galerkin finite element method (DG-FEM) to solve the linear and nonlinear elastic wave equation in heterogeneous media with arbitrary high order accuracy in space on unstructured triangular or quadrilateral meshes is presented. This DG-FEM method combines the geometrical flexibility of the finite element method, and the high parallelization potentiality and strongly nonlinear wave phenomena simulation capability of the finite volume method, required for nonlinear elastodynamics simulations. In order to facilitate the implementation based on a numerical scheme developed for electromagnetic applications, the equations of nonlinear elastodynamics have been written in a conservative form. The adopted formalism allows the introduction of different kinds of elastic nonlinearities, such as the classical quadratic and cubic nonlinearities, or the quadratic hysteretic nonlinearities. Absorbing layers perfectly matched to the calculation domain of the nearly perfectly matched layers type have been introduced to simulate, when needed, semi-infinite or infinite media. The developed DG-FEM scheme has been verified by means of a comparison with analytical solutions and numerical results already published in the literature for simple geometrical configurations: Lamb's problem and plane wave nonlinear propagation. 相似文献
17.
We suggest here a least-change correction to available finite element (FE) solution. This postprocessing procedure is aimed at recovering the monotonicity and some other important properties that may not be exhibited by the FE solution. Although our approach is presented for FEs, it admits natural extension to other numerical schemes, such as finite differences and finite volumes. For the postprocessing, a priori information about the monotonicity is assumed to be available, either for the whole domain or for a subdomain where the lost monotonicity is to be recovered. The obvious requirement is that such information is to be obtained without involving the exact solution, e.g. from expected symmetries of this solution.The postprocessing is based on solving a monotonic regression problem with some extra constraints. One of them is a linear equality-type constraint that models the conservativity requirement. The other ones are box-type constraints, and they originate from the discrete maximum principle. The resulting postprocessing problem is a large scale quadratic optimization problem. It is proved that the postprocessed FE solution preserves the accuracy of the discrete FE approximation.We introduce an algorithm for solving the postprocessing problem. It can be viewed as a dual ascent method based on the Lagrangian relaxation of the equality constraint. We justify theoretically its correctness. Its efficiency is demonstrated by the presented results of numerical experiments. 相似文献
18.
In this paper, a compact finite difference scheme for the fractional sub-diffusion equations is derived. After a transformation of the original problem, the L1 discretization is applied for the time-fractional part and fourth-order accuracy compact approximation for the second-order space derivative. The unique solvability of the difference solution is discussed. The stability and convergence of the finite difference scheme in maximum norm are proved using the energy method, where a new inner product is introduced for the theoretical analysis. The technique is quite novel and different from previous analytical methods. Finally, a numerical example is provided to show the effectiveness and accuracy of the method. 相似文献
19.
The theory of quantum fluid dynamics (QFD) helps nanotechnology engineers to understand the physical effect of quantum forces. Although the governing equations of quantum fluid dynamics and classical fluid mechanics have the same form, there are two numerical simulation problems must be solved in QFD. The first is that the quantum potential term becomes singular and causes a divergence in the numerical simulation when the probability density is very small and close to zero. The second is that the unitarity in the time evolution of the quantum wave packet is significant. Accurate numerical evaluations are critical to the simulations of the flow fields that are generated by various quantum fluid systems.A finite volume scheme is developed herein to solve the quantum hydrodynamic equations of motion, which significantly improve the accuracy and stability of this method. The QFD equation is numerically implemented within the Eulerian method. A third-order modified Osher–Chakravarthy (MOC) upwind-centered finite volume scheme was constructed for conservation law to evaluate the convective terms, and a second-order central finite volume scheme was used to map the quantum potential field. An explicit Runge–Kutta method is used to perform the time integration to achieve fast convergence of the proposed scheme.In order to meet the numerical result can conform to the physical phenomenon and avoid numerical divergence happening due to extremely low probability density, the minimum value setting of probability density must exceed zero and smaller than certain value. The optimal value was found in the proposed numerical approach to maintain a converging numerical simulation when the minimum probability density is 10?5 to 10?12. The normalization of the wave packet remains close to unity through a long numerical simulation and the deviations from 1.0 is about 10?4.To check the QFD finite difference numerical computations, one- and two-dimensional particle motions were solved for an Eckart barrier and a downhill ramp barrier, respectively. The results were compared to the solution of the Schrödinger equation, using the same potentials, which was obtained using by a finite difference method. Finally, the new approach was applied to simulate a quantum nanojet system and offer more intact theory in quantum computational fluid dynamics. 相似文献