首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 127 毫秒
1.
The Von-Karman vortex street flow behind a circular cylinder issuccessfully solved at R_e=100 for incompressible N-S equations in theform of stream function and vorticity.The method adopted here involvesADI-BGE schemes,a recent treatment of body boundary conditions anda new acceptable perturbation model.The solutions of the flow fieldobtained here are better than previous numerical ones and agree with theexperimental data closer.Therefore,it can be used to study sound fieldnumerically.The Curle equation is directly and numerically solved to ana-lyze the distribution of sound sources and the characteristics of their radi-ated sound field at R_e=100 and to compare with the solution fromthe theory of vortex sound.It is shown that this method is successful.  相似文献   

2.
李兴国  戴保东  王灵卉 《中国物理 B》2010,19(12):120202-120202
In this paper,a meshfree boundary integral equation(BIE) method,called the moving Kriging interpolationbased boundary node method(MKIBNM),is developed for solving two-dimensional potential problems.This study combines the BIE method with the moving Kriging interpolation to present a boundary-type meshfree method,and the corresponding formulae of the MKIBNM are derived.In the present method,the moving Kriging interpolation is applied instead of the traditional moving least-square approximation to overcome Kronecker’s delta property,then the boundary conditions can be imposed directly and easily.To verify the accuracy and stability of the present formulation,three selected numerical examples are presented to demonstrate the efficiency of MKIBNM numerically.  相似文献   

3.
柴振华  施保昌  郑林 《中国物理》2006,15(8):1855-1863
By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number (Re) flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50 000 in the literature; however, this paper can readily extend the maximum to 1000 000 with the above combination.  相似文献   

4.
The central problem of the lattice Boltzmann method (LBM) is to construct a discrete equilibrium. In this paper, a multi-speed 1D cell-model of Boltzmann equation is proposed, in which the cell-population equilibrium, a direct non- negative approximation to the continuous Maxwellian distribution, plays an important part. By applying the explicit one-order Chapman-Enskog distribution, the model reduces the transportation and collision, two basic evolution steps in LBM, to the transportation of the non-equilibrium distribution. Furthermore, 1D dam-break problem is performed and the numerical results agree well with the analytic solutions.  相似文献   

5.
陈丽  程玉民 《中国物理 B》2010,19(9):90204-090204
On the basis of the reproducing kernel particle method (RKPM), a new meshless method, which is called the complex variable reproducing kernel particle method (CVRKPM), for two-dimensional elastodynamics is presented in this paper. The advantages of the CVRKPM are that the correction function of a two-dimensional problem is formed with one-dimensional basis function when the shape function is obtained. The Galerkin weak form is employed to obtain the discretised system equations, and implicit time integration method, which is the Newmark method, is used for time history analysis. And the penalty method is employed to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional elastodynamics are obtained. Three numerical examples of two-dimensional elastodynamics are presented, and the CVRKPM results are compared with the ones of the RKPM and analytical solutions. It is evident that the numerical results of the CVRKPM are in excellent agreement with the analytical solution, and that the CVRKPM has greater precision than the RKPM.  相似文献   

6.
A lattice Boltzmann model (LBM) has been developed for simulating magnetohydrodynamics (MHD) along the line of Dellar [J. Comput. Phys. 179 (2002)95]. In this model the magnetic field is presented by a vector valued magnetic distribution function which obeys a vector Boltzmann equation. The truncated error of the equilibrium distribution in the present model is up to order O(u^4) in velocity u rather than the usual 0(u^3). For verification, the model is applied to solve the shock tube problem and the main features of the flow predicted by the model are found to compare well with the corresponding results obtained with high-order semi-discrete schemes [J. Comput. Phys. 201 (2004) 261]. The numerical experiments have also shown that the present LBM model with the equilibrium distribution truncated at O(u^4) performs much better in terms of numerical stability than those truncated at O(u^3).  相似文献   

7.
A generalized Fisher equation (GFE) relates the time derivative of the average of the intrinsic rate of growth to its variance. The exact mathematical result of the GFE has been widely used in population dynamics and genetics, where it originated. Many researchers have studied the numerical solutions of the GFE, up to now. In this paper, we introduce an element-free Galerkin (EFG) method based on the moving least-square approximation to approximate positive solutions of the GFE from population dynamics. Compared with other numerical methods, the EFG method for the GFE needs only scattered nodes instead of meshing the domain of the problem. The Galerkin weak form is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. In comparison with the traditional method, numerical solutions show that the new method has higher accuracy and better convergence. Several numerical examples are presented to demonstrate the effectiveness of the method.  相似文献   

8.
We consider the density dependent diffusion Nagumo equation, where the diffusion coefficient is a simple power function. This equation is used in modelling electrical pulse propagation in nerve axons and in population genetics (amongst other areas). In the present paper, the δ-expansion method is applied to a travelling wave reduction of the problem, so that we may obtain globally valid perturbation solutions (in the sense that the perturbation solutions are valid over the entire infinite domain, not just locally; hence the results are a generalization of the local solutions considered recently in the literature). The resulting boundary value problem is solved on the real line subject to conditions at z →±∞. Whenever a perturbative method is applied, it is important to discuss the accuracy and convergence properties of the resulting perturbation expansions. We compare our results with those of two different numerical methods (designed for initial and boundary value problems, respectively) and deduce that the perturbation expansions agree with the numerical results after a reasonable number of iterations. Finally, we are able to discuss the influence of the wave speed c and the asymptotic concentration value α on the obtained solutions. Upon recasting the density dependent diffusion Nagumo equation as a two-dimensional dynamical system, we are also able to discuss the influence of the nonlinear density dependence (which is governed by a power-law parameter m) on oscillations of the travelling wave solutions.  相似文献   

9.
王聚丰  孙凤欣  程荣军 《中国物理 B》2010,19(6):60201-060201
The present paper deals with the numerical solution of the third-order nonlinear KdV equation using the elementfree Galerkin (EFG) method which is based on the moving least-squares approximation.A variational method is used to obtain discrete equations,and the essential boundary conditions are enforced by the penalty method.Compared with numerical methods based on mesh,the EFG method for KdV equations needs only scattered nodes instead of meshing the domain of the problem.It does not require any element connectivity and does not suffer much degradation in accuracy when nodal arrangements are very irregular.The effectiveness of the EFG method for the KdV equation is investigated by two numerical examples in this paper.  相似文献   

10.
The present paper deals with the numerical solution of the coupled Schrodinger-KdV equations using the elementfree Galerkin (EFG) method which is based on the moving least-square approximation. Instead of traditional mesh oriented methods such as the finite difference method (FDM) and the finite element method (FEM), this method needs only scattered nodes in the domain. For this scheme, a variational method is used to obtain discrete equations and the essential boundary conditions are enforced by the penalty method. In numerical experiments, the results are presented and compared with the findings of the finite element method, the radial basis functions method, and an analytical solution to confirm the good accuracy of the presented scheme.  相似文献   

11.
周晓阳 《物理学报》2008,57(1):238-248
The central problem of the lattice Boltzmann method (LBM) is to construct a discrete equilibrium. In this paper, a multi-speed 1D cell-model of Boltzmann equation is proposed, in which the cell-population equilibrium, a direct non-negative approximation to the continuous Maxwellian distribution, plays an important part. By applying the explicit one-order Chapman--Enskog distribution, the model reduces the transportation and collision, two basic evolution steps in LBM, to the transportation of the non-equilibrium distribution. Furthermore, 1D dam-break problem is performed and the numerical results agree well with the analytic solutions.  相似文献   

12.
The formulation of a consistent thermohydrodynamics with a discrete model of the Boltzmann equation requires the representation of the velocity moments up to the fourth order. Space-filling discrete sets of velocities with increasing accuracy were obtained using a systematic approach in accordance with a quadrature method based on prescribed abscissas (Philippi et al., Phys. Rev. E, 73 (5), n. 056702, 2006). These sets of velocities are suitable for collision-propagation schemes, where the discrete velocity and physical spaces are coupled and the Courant number is unitary. The space-filling requirement leads to sets of discrete velocities which can be large in thermal models. In this work, although the discrete sets of velocities are also obtained with a quadrature method based on prescribed abscissas, the lattices are not required to be space-filling. This leads to a reduced number of discrete velocities for the same approximation order but requires the use of an alternative numerical scheme. The use of finite difference schemes for the advection term in the continuous Boltzmann equation has shown to have some advantages with respect to the collision-propagation LBM method by freeing the Courant number from its unitary value and reducing the discretization error. In this work, a second order Runge-Kutta method was used for the simulation of the Sod's shock tube problem, the Couette flow and the Lid-driven cavity flow. Boundary conditions without velocity slip and temperature jumps were written for these discrete Boltzmann equation by splitting the velocity distribution function into an equilibrium and a non-equilibrium part. The equilibrium part was set using the local velocity and temperature at the wall and the non-equilibrium part by extrapolating the non-equilibrium moments to the wall sites.  相似文献   

13.
In this article, a discrete effect in the thermal Lattice BGK two-speed model is studied. These effects are due to the non-equilibrium state in the particle distribution function, and the non-equilibrium occurs near walls. The mechanism of the LBM counterpart of the thermal creep flow, which appears due to the temperature gradient of the boundary in rarefied gases, is clarified analytically and numerical calculations are performed for some cases. A technique for eliminating this effect is also shown.  相似文献   

14.
A new implementation of boundary condition based on the half-covolume and bounce-back rule for the non-equilibrium distribution function for the finite volume LBM is proposed here. The numerical simulation results for the expansion channel flow and driven cavity problem indicate that this method is workable for arbitrary meshes. In addition, the fourth order Runge–Kutta scheme is found to be a practical way in the LBM to accelerate the calculation speed.  相似文献   

15.
马致遥  单锋  章东 《声学学报》2018,43(2):217-223
高强度聚焦超声(HIFU)是一种新型的无创治疗肿瘤新技术,其中换能器声场数值计算能够为HIFU治疗提供重要的依据。传统非线性KZK和SBE模型广泛应用于换能器声场数值计算,但依然存在某些不足。我们采用一种介观尺度的新型流体力学方法,即格子Boltzmann方法(LBM),基于2维9离散速度(D2Q9)格子构建了轴对称多弛豫参数LBM模型,并通过调节弛豫参数分析其对模型的影响;利用该模型对两个具有不同张角的球面聚焦换能器的声场进行数值模拟,并与KZK和SBE模型的计算结果进行比较。结果表明LBM模型能够很好地描述超声波的激发和传播机制,从流体力学的角度描述聚焦声场的分布,具有清晰的物理意义,且计算过程不受换能器张角的限制,在换能器声场的理论分析和模拟计算及其在HIFU治疗中的应用有着积极的意义。   相似文献   

16.
翟旭军  赵凯 《计算物理》2012,29(3):347-353
针对非平衡态外推边界处理方法在某些条件下发生质量泄漏的问题,结合质量守恒定律对其进行修正,建立一种新的质量守恒型的边界处理格式.其基本思想是在边界处定义一个虚拟密度来修正平衡态分布函数,使得系统质量守恒.以重力作用下的密度场及磁流体模拟问题为例,利用该格式进行分析,讨论质量泄露问题的本质及各种因素的影响.通过对三种典型的流动:二维Poiseuille流、Couette流、Womersley流的模拟,验证该边界格式的计算精度和稳定性.  相似文献   

17.
We propose a boundary scheme for addressing multi-mechanism flow in a porous medium in slip and early transition flow regimes, which is frequently encountered in shale gas reservoirs. Micro-gaseous flow in organic-rich shale involves a complex flow mechanism. A self-developed boundary scheme that combines the non-equilibrium extrapolation scheme and the combined diffusive reflection and bounce-back scheme(half-way DBB) to embed the Langmuir slip boundary into the single-relaxation-time lattice Boltzmann method(SRT-LBM) enables us to describe this process, namely, the coupling effect of micro-gaseous flow and surface diffusion in organic-rich nanoscale pores. The present LBM model comes with the careful consideration of the local Knudsen number, local pressure gradient, viscosity correction model, and regularization procedure to account for the rarefied gas flows in irregular pores. Its validity and accuracy are verified by several benchmarking cases, and the calculated results by this boundary scheme accord well with our analytical solutions.This boundary scheme shows a higher accuracy than the existing studies. Additionally, a subiteration strategy is presented to tackle the coupled micro-gaseous flow and surface diffusion, which necessitates the iteration process matching of these two mechanisms. The multi-mechanism flow in the self-developed irregular pores is also numerically investigated and analyzed over a wide range of parameters. The results indicate that the present model can effectively capture the coupling effect of micro-gaseous flow and surface diffusion in a tree-like porous medium.  相似文献   

18.
王佐  刘雁  张家忠 《物理学报》2016,65(1):14703-014703
为提高采用二维九速离散速度模型的格子Boltzmann方法 (LBM)模拟微尺度流动中非线性现象的精度和效率,引入Dongari等提出的有效平均分子自由程对黏性进行修正(Dongari N,Zhang Y H,Reese J M2011 J.Fluids Eng.133 071101);并针对以往研究微尺度流动时采用边界处理格式含有离散误差的问题,采用多松弛系数格子Boltzmann方法结合二阶滑移边界条件,对微尺度Couette流动和周期性Poiseuille流动进行模拟,并将速度分布以及质量流量等模拟结果与直接模拟蒙特卡罗方法模拟数据、线性Boltzmann方程的数值解以及现有的LBM模型模拟结果进行对比.结果表明,相对于现有的LBM模型,引入新的修正函数所建立的有效黏性多松弛系数LBM模型有效提高了LBM模拟过渡区的微尺度流动中的非线性现象的能力.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号