中子输运方程数值解与Burgers方程格子Boltzmann方法研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
本文分两部分。第一部分研究中子输运方程的数值方法,主要讨论
    输运方程离散纵标方法,积分方程方法和球谐近似方法中当前大家
    关注的研究课题;第二部分针对Burgers方程研究了格子Boltzmann方
    法的收敛性和稳定性问题。
     全文共七章。第一章为前言,简要介绍了中子输运方程和格子
    Boltzmann理论的研究概况以及本文所讨论的基本内容。第二章研究
    离散纵标方法。采用将角变量与空间变量的离散组合起来的方法,
    给出了常见的输运方程数值格式的误差阶。结果表明,对标量通量
    而言,这些数值方法的精度均达不到二阶。在第三章,对平板几何输
    运方程,我们将解的奇异部分分离出来,应用算子逼近理论构造了
    二阶精度的数值格式,并证明了解的存在唯一性。第四章研究球几
    何输运方程的数值解问题,也分为独立的两部分。前一部分利用调
    和分析方法精细地研究了解在球心附近的正则性和在边界附近的奇
    异性,将解的奇异部分分离出来,并且构造了二阶精度数值格式。后
    一部分考虑离散纵标方法中的内迭代问题,给出了步函数格式内迭
    代的收敛性证明。第五章应用奇异摄动理论和Case的奇异本征值理
    论,给出了P_2方程的边界条件,数值结果表明我们的条件较Marshak
    的边界条件更精确。第六,七章分别研究了格子Boltzmann(LB)方法对
    一维和二维Burgers方程的模拟,构造了带有修正项的BGK型的LB
    方程。第六章,针对一维问题的特殊性,应用离散泛函分析方法,证
    明了LB方程的解在变换后收敛到Burgers方程的解。第七章,我们证
    明了二维LB方程满足极值原理和L~1压缩性,从而证明了解在L~1意
    义下的稳定性。章末的数值例子表明LB方法的数值精度与经典的二
    阶守恒方法符合得非常好。
This paper consists of two parts. The first part is devoted to the study of nu-
     merical solutions of neutral transport equations, and the other to the investigation
     of Lattice Boltzmann (LB) methods for modeling Burgers equation.
    
     The developments of the transport equation and their numerical methods,
     and Lattice Boltzmann methods have been introduced in Chapter 1, in which
     main work of this paper is also described. In Chapter 2 the error estimates for
     some common used schemes of the transport equations with combined spatial and
     angular approximations are considered. The conclusions show that error order of
     scalar flux in all of these schemes can not get second order accurate. In Chapter
     3, we consider the properties of solution of transport equation and decompose
     the singularities of solution near boundary and interface. Furthermore, we con-
     struct the second order accurate scheme and prove the existence and uniqueness
     of numerical solutions.
    
     Numerical solution problems of transport equation in sphere geometry are
     discussed in Chapter 4 which is divided into two parts too. In the first part,
     we research the regularities of solution near the spherical center and the singu-
     larities near the spherical surface. We get the coefficients of weak singular term
     and construct a second order accurate numerical scheme. In the second part we
     discuss the inner iteration schemes of discrete ordinate equations and prove the
     convergence of iterative sequence.
    
     In Chapter 5 we study the boundary conditions of P2 equations using singu-
     lar perturbation method and Case singular eigenvalue technologies. The numer-
     ical examples show that asymptotic boundary condition is more accurate than
     Marshak boundary condition.
    
     In Chapter 6 and 7, we research the mathematics theories of Lattice Boltz
    
    
    mann methods. LB equations fOr modeling 1D Burgers equations are constructed
    using modified BGK model in Chapter 6. By discrete functionaI analysis method.
    we prove the convergence of solutions for LB equation. In Chapter 7, we prove the
    solutions of 2D LB equations satisfy maxmum value principle and Ll constraction,
    and so we get the stability of numerical solutions.
引文
[1]J.C. Maxwell, Phil. Trans. R. Soc. 157(1866); reprinted in: The Scientific Papers of J. C. Maxwell, Dover Publications, New York, 1965.
    [2]L. Boltzmann, Lectures on Gas Theory, English translation from the German, University of California Press, Berkeley, 1964.
    [3]D. Hilbert, Math Ann. 72(1912),562-572. (见周毓麟,输运方程讲义,未发表)
    [4]S. Chandrasekhar, Radiative Transfer, Oxford, Clarendon Press, 1950.
    [5]K.M. Case and P. F. Zweifel, Linear Transport Theory, Addison-Wesley, Reading, Mass., 1967.
    [6]R. Dautray and J. L. Lions, Mathematical Analysis and Numerical Methods for Science and Technology, Vol 6, Evolution Problems II, Springer-Verlag,Berlin Heidelbeng, 1993.
    [7]J.J. Duderstadt and W. R. Martin, Transport Theory, Wiley-Interscience,New York, 1979.
    [8]V. I. Agoshkov, Boundary Value Problems for Transport Equations, Birkhanser, 1998.
    [9]M. Borysiewicz and G. Spiga, Smoothness of the Solution to 2-D Neutron Transport Equation. Preprint. Warsaw: Inst. of Nuclear Research Computing Center, 1974.
    [10]V. P. Shutyaev, Some Regularity Properties in Energy Variable of a Solution of Transport Equation, In: Numerical methods and Statistical modeling in transport theory, Comp. Center of Siberian Branch of USSR Academy of Sci.,Novosibirsk, 1980, pp. 81-82.
    [11]V.P. Shutyaev, Regularity Properties of Solutions of nonsatationary Transport Problem in a Plane Layer, In: Problems of Theory and Numerical Solution of Transport Problems, Dept. of Numerical Mathematics, USSR Academy of Sci., Moscow, 1983, pp. 148-159.
    
    
    [12]A. Douglis, The Solutions of multidimensional generalized transport equations and Their Calculation by difference Methods. Proc. Symp. Numberical Solutions of Partial Diff. Eq. Maryland, Ed. 1. II. Brambi. L., 1966.
    [13]R. Sanchez and N. J. McCormick, A Review of Neutron Transport Approximations, Nucl. Sci. Eng., 80(1982), pp.481-535.
    [14]B. G. Carlson, Solution of the Transport Equation by Approximations, Los Alamos Scientic Laboratory Report LA-1599, 1953.
    [15]M. A. Golberg, Numerical Solution of Integral Equations, Plenum Press, New York, 1990.
    [16]C. V. Pao, On Nonlinear Neutron Transport Equations, J. Math. Anal. Appl.42(1973), pp. 578-593.
    [17]郭柏灵,沈隆钧,非线性中子输运方程的Galerkin有限元解,应用数学与计算数学学报,3(1989),pp.8-20.
    [18]G. W. Yuan, Convergence of P_N Approximation for the Neutron Transport Equation with Reflective Boundary Condition, J. Math. Phys., 41(2),2000,pp. 867-874.
    [19]E. W. Larsen and W. F. Miller, Jr. , Convergence rates of spatial difference equations for the discrete-ordinates neutron transport equations in slab geomety, Nucl. Sci. Eng., 73(1980), pp.481-535.
    [20]S. M. Lee and R. Vaidyanathan, Comparison of the Order of Approximation in Several Spatial Difference Schemes for the Discrete-ordinates Transport Equation in One-Dimensional Plane Geometry, Nucl. Sci. Eng., 76(1980), pp.1-9.
    [21]B. Neta and H. D. Victory, Jr., The Converengece Analysis for Sixth-Order Methods for Solving Discrete-Ordinates Slab Transport Equations, Numer.Funct. Anal. and Optimiz., 5(1)(1982), pp. 85-126.
    [23]J. Pitk(?)ranta, Asmptotic Behavior of the Solution of the Integral Transport Equation in the Vicinity of a Curved Material Interface, SIAMJ. Appl. Math.36(1979), pp. 200-218.
    
    
    [24] H. B. Keller and B. Wendroff, On the Formulation and Analysis of Numerical Methods for Time Dependent Transport Equations, Comm. Pure Appl.Math., 10(1957), pp.567-582.
    [25] K. Kobayashi and H. Nishihara, The Solution of Integral Transport Equations in Cylindrical Geometry Using the Gaussian Quadrature Formula, J. Nucl.Energy, 18(1964), pp.513-522.
    [26] G.I. Marchuk and V. I. Lebedev, Numerical Methods in the Theory of Neutron Transport, Revised Second Edition, London, harwood Academic, 1986.
    [27] G. I. Marchuk, Methods of Numerical Mathematics, Berlin, Springer, 1982.
    [28] H. B. Keller, On the pointwise convergence of the discrete-ordinates solution method, J. Soc. Indust. Appl. Math., 8(1960), pp.43-73.
    [29] B. Wendroff, On the convergence of the discrete-ordinates method, J. Soc.Indust. Appl. Math., 8(1960), pp.508-513.
    [30] P. M. Anselone, Convergence of Chandrasekhar's method for inhomogeneous transfer problems, J. Math. Mech., 10(1961), pp. 537-546.
    [31] P. Nelson, Jr. and H. D. Victory, Jr., Theoretical properties of one-dimensional discrete-ordinates, SIAM J.Numer. Anal., 16(1979), pp. 270-283.
    [32] P. Nelson, Jr. and H. D. Victory, Jr., Convergence of the discrete-ordinates method for anisotropically scattering multiplying particles in a subcritical slab, SIAM J.Numer. Anal., 10(1973), pp. 175-181.
    [33] A. Marsak, On the solution of the radiation transfer eqaution by the method of the discrete ordinates in non-homogeneous medium, Trans. Estonian Acad. Sci., 30(1981), pp. 191-201.
    [34] W.Larsen and P. Nelson, Jr., Finite-difference approximations and supercon vergence for the discrete-ordinates in slab geometry, SIAM J. Numer. Anal. 19(1982), pp. 334-348.
    
    
    [36] M.Asadzadeh, Analysis of a Fully Discrete Scheme for Neutron Transport in Two-Dimensional Geometry, SIAM J.Numer. Anal.23(1986),pp.543-561.
    [37] H.D.Victory, Jr and K.Ganguly, On finite-difference methods for solving discrete-okinates transport equations, SIAM J.Numer. Anal.23(1986),pp.78-108.
    [38] 匡继昌,常用不等式,湖南教育出版社,1993.
    [39] R.Kress, Linear Integral Equations, Spinger-Verlag, World Publishing Corp,1991.
    [40] G.Vainikko, Multidimensional wdakly singular integral equation, SpringerVerlag:Berlin Heidelberg,1993.
    [41] G.J.Mitsis, Transport Solutions to the Monoenergetic Critical Problems, ANS-6787(TiD-4500)1963.
    [42] D.E.Cullen, Numerical Solutions of the Linear Integral Boltzmann Equation, Nucl. Sci. Eng. 53(1974),pp.93-106.
    [43] S.K.Loyalka, A Numerical method for solving Integral Equations of Neutron Transport, Nucl. Sci.Eng.56(1975),pp.317-319.
    [44] H.G.Kaper,Integral Operators Generated by Exponential integrals, J.Math. Anal. Appl.(1977),pp,216-226.
    [45] H.G.Kaper and R.B.Kellogg, Asmptotic Behavior of the Solution of the Integral Transport Equation in Slab Geometry, SIAM J. Appl. Math. 32(1977),pp.191-200.
    [46] H.G.Kaper, Asymptotic Calculation of Two Families of Integrals, J.Math.Anal.Appl.59(1977),pp.461-471.
    [47] A.V.Voronkov and E.P.Sychugova, CDSN-Method for Solving the Transport Equation, Trans. Theory Stat. Phys.,22(23)(1993),pp.221-245.
    [48] I.Ohnishi, Finith Element Solution Technique for Neutron Transprot Equation(F_N Approximation), Proc.Sem. Numerical Reactor Calculations, 1972, pp.629-642.
    
    
    [49]夏道行,吴卓人,严绍宗和舒五昌,实变函数与泛函分析(下册),高等教育出版社,1979.
    [50]P. G. Ciarlet, The finite Element Mothod for Elliptic Problems, North- Holland. Amsterdam, 1978.
    [51]周民强,实变函数(第二版),北京大学出版社,1995.
    [52]周民强,调和分析讲义(实变方法),北京大学出版社,1999.
    [53]H. Greenspan, C. N. Kelber and D. Okrent, Computing Methods in Reactor Physics, 1968, Chp. 3.
    [54]P. Nelson, Convergence of inner iterations for finite-difference approximations to the linear transport equation, SIAM J. Number. Anal., 23(1986), pp. 1017-1022.
    [55]徐树芳,矩阵计算的理论与方法(第二版),北京大学出版社,1999.
    [56]杜书华,张树华,冯庭桂,刑静茹编著,输运问题的计算机模拟,湖南科学技术出版社,1989.
    [57]G.I.贝尔和S.格拉斯登,核反应堆理论.原子能出版社,1979.
    [58]R. S. Varga, Matrix Iterative Analysis, Printice-Hall, Engle-Wood Cliffs, NJ,1962.
    [59]S. V. G. Menon and D. C. Sahni, Convergence of Discrete Ordinates Iteration Scheme, Transp. Theory Stat. Phys., 14(1985), pp. 353-372.
    [60]Guang-wei Yuan, Zhi-jun Shen and Yu-lin Zhou, Convergence of Inner Iterations Scheme of the Discrete Ordinate Method for Neutron Transport Equations, preprint.
    [61]袁光伟,沈智军,沈隆钧,输运理论中扩散近拟的边界条件,计算机物理实验室年报,1998,pp.303-309.
    
    
    [62]袁光伟,沈智军,沈隆钧,周毓麟,中子输运方程SP_2方法稳定性分析,第七届反应堆数值计算与粒子输运学术交流会论文集,1998.pp.49-54.(第三章)
    [63]沈智军,袁光伟,沈隆钧,周毓麟,中子输运方程SP_2方法的并行计算,第七届反应堆数值计算与粒子输运学术交流会论文集,1998.pp.55-58.(第三章).
    [64]G. J. Habetler and B. J. Matkowsky, Uniform Asympototic Expansions in Transport Theory with Small Mean Free Paths, and the Diffusion Approximation, J. Math. Phys., 16(1975), pp. 846-854.
    [65]U. Shin and W. F. Miller, Jr., The Time-Dependent Simplified P_2 Equations:Asymptoic Analyses and Numerical Ecpriments, Nucl. Sci. Eng., 128(1998),pp.27-46.
    [66]E. W. Larsen, J.E. Morel and W.F.Miller,Jr., Asymptotic Solution of Numerical Transport Problems in Optically Thick, Diffusive Regimes, J. Comp.Phys., 69(1987), pp.283-324.
    [67]G. C. Pomraning, Variational Boundary Conditions for the Spherical Harmonics Approximation to the Neutron Transport Equation, Ann. Phys., 27(1964),pp.193-215.
    [68]G. C. Pomraning, An Improved Free-Surface Boundary Condition for the P-3 Approximation, Nucl. Sci. Eng., 18(1964), pp.528-530.
    [69]E. M. Gelbard, Application of Spherical Harmonics Method to Reactor Problems, WAPD-BT-20(Sep. 1960).
    [70]E. W. Larsen, J. M. McGhee and J. E. Morel, The Simplified P_N Equations as an Asmptotic Limit of the Transport Equation, Trans. Am. Nucl. Soc., 66(1992), pp.232-234.
    [71]E. W. Larsen and G. C. Pomraning, The P_N Theory as an Asymptotic Limit of Transport Theory in Planar Geometry-I: Analysis, Nucl. Sci. Eng., 109(1991), pp. 49-75.
    [72]R. P. Rulko and E. W. Larsen, Variational Derivation and Numerical Analysis of P_2 Theory in Planar Geometry, Nucl. Sci. Eng., 114(1993), pp.271-285.
    
    
    [73]R. P. Rulko and E. W. Larsen, Variational Derivation of Multigroup P_1 Equations and Boundary Conditions, Proc. Topl. Mtg. Advances in Mathematics,Computations, and Reactor Physics, Pittsburgh, 4(1991), pp. 2-10.
    [74]李家春,周显初,数学物理中的渐进方法,科学出版社,1998.
    [75]G. C. Pomraning, Near-Infinite-Medium solutions of the Equation of transfer,J. Quant. Spectrose. Radiat. Transfer, 44(3)(1990), pp.317-338.
    [76]F. Malvagi and G.C. Pomraning, Initial and Boundary Conditions for Diffusive Linear Transport Problems, J. Math. Phys., 32(3)(1991), pp.805-820.
    [77]G. C. Pomraning, Diffusion Theory via Asymptotics, Trans. Theory. Stats Phys, 18(1990), pp. 383-438.
    [78]T. Lefvert, A Multiple Scatting and Collision Probability Method for Neutron Transport Problems, Nucl. Sci. Eng., 42(1970), pp.267-271.
    [79]R. Daniel, V. C. Badham and G. C. Pomraning, A Review of Transport Theory Methods in Nonhomogeneous Media, Trans. Theory Stats. Phys.,14(2)(1985), pp.125-166.
    [80]B. R. Wienke and R. E. Hiromoto, Parallel S_n Transport Algorithms, Transp.Theory & Stat. Phys., 15(1&2), 1986, pp. 49-59.
    [81]杜明笙,冯铁凯,傅连祥,曹常恕,刘玉兰,解二维输运问题的三角网间断有限元方法,计算物理,1(1984),pp.84-107.
    [82]W. F. Miller, jr., E. E. Lewis and E. C. Rossow, Nucl. Sci. Eng., 51(1973),148-156.
    [84]Yuan Guangwei, Convergence of DP_N Approximation to Neutron Transport Equations, preprint.
    [85]S. Chapman and T. G. Cowling, The Mathematical Theory of Non-Uniform Gases, Cambridge University Press. Cambridge, 1952.
    [86]C. Cercignani, Mathematical Methods in Kinetic Theory, Plenum Press, 1990.(Second Edition).
    
    
    [87]G. R. McNamara and G. Zanetti, Use of Boltzmann Equation to Sirnulate Lattice Gas Automata, Phys Rev. Lett., 61(1988), pp. 2332-2335.
    [88]P. Bhatnagar, E. P. Gross and M. K. Krook, A Model for Collision Processes in Gases. I. Small Amplitude Processes in Charged and Neutral One-Component Systems, Phys. Rev. 94(1954), pp.511-525.
    [89]H. D. Chen, S. Y. Chen and W. Matthaeus, Recovery of the Navier-Stokes Euations Using a Lattice-Gas Boltzmann Method, Phys. Rev. A 45(1992),pp. R5339-R5342.
    [90]Y. H. Qian, D. d'Humieres and P. Lallemand, Lattice BGK Models for NavierStokes Equation, Europhys. Lett., 17(1992), pp. 479-484.
    [91]U. Frisch, B. Hasslacher and Y. Pomeau, Lattice Gas Automata for the NovierStokes Equation, Phys. Rev. Lett., 56(1986), pp. 1505-1508.
    [92]H.Bateman, Some recent researches on the motion of fluids, Monthly Weather Review, 43(1915), pp. 163-170.
    [93]J. Burgers, A Mathematical model ilustrating the theory of turbulence, Advances in Applied Mechanics, Academic Press, 1948.
    [94]P. Lagerstrom, J. D. Cole and L. Trilling, Problems in the Theory of Viscous Compressible Fluids, Monograph, California Institute of Technology, 1949.
    [95]J. D. Cole, On a Quasi-linear Parabolic Equation Occurring in Aerodynamics,Q. Appl. Math. 9(1951), pp. 225-236.
    [96]E. Hopf, The Partial Differential Equation u_t+uu_x=μu_xx. Comm. Pure. Appl. Math., 3(1950), pp. 201-230.
    [97]O. Oleinik, Discontinuous solutions nonlinear Differential Equations, English transl, in Amer. Math. Soc. Transl. Ser. 2, 26(1957), 95-172.1991.
    [98]应隆安,滕振寰,双曲型守恒律方程及其差分方法,科学出版社,1991.
    [99]M. N. Kogan, Rarefied Gas Dynamics, Plenum, New York,1969.
    
    
    [100]Zhou Yulin, Applications of the Discrete Functional Analysis to the Finite Difference Method, International Academic Publishers. 1990.
    [101]Zhou Ynlin,Theory of Difference Schemes with Nonuniform Meshes, Annual Report (1993), Lab. Comp. Phys., IAPCM, pp. 216-238.
    [102]Shen Zhijun, Yuan Guangwei and Shen Long, jun, Lattice Boltzmann Method for Burgers Equation, Chinese J. Comp. Phys., 3(2000), pp. 56-64.
    [103]G. A. Sod, Numerical Methods in Fluid Dynamics: Initial and Initial Boundary-Value Problems, Cambridge University Press, Cambridge, 1988.
    [104]Y. Chen, H.Ohashi and M. Akiyama, Thermal Lattice Bhatnagar-GrossKrook Model Without Nonlinear Deviations in Macrodynamic Equations, Phys. Rev. E., 50(1994), pp. 2776-2783.
    [105]J. Sterling and S. Chen, Stability Analysis of Lattice Boltzmann Methods, preprint.
    [106]M. G. Ancona, Fully-Lagrangian and Lattice-Boltzmann Methods for Solving Systems of Conservation Equations, J. Comp. Phys., 115(1994), pp. 107-120.
    [107]S. P. Dawson, S. Chen and G. D. Doolen, Lattice Boltzmann Computations for Reaction-Diffusion Equation, J. Chem. Phys. 98(1993), pp. 1514-1523.
    [108]D. Grunau, S. Chen and K. Eggert, A Lattice Boltzmarnn Model for Multiphase Fluid Flows, Phys. Fluids., 5(10)(1993), pp. 2557-2562.
    [109]B. H. Elton, C. D. Levermore and G. H. Rodrigue, Convergence of ConvectiveDiffusion Lattice Boltzmann Methods, SIAM J. Numer. Anal. 32(1995), pp.1327-1354.
    [110]B. H. Elton, Comparisons of Lattice Boltzmannand Finite Difference Methods for a Two -Dimensional Viscous Burgers Equation, SIAM J. Sci. Comput.,17(1996), pp. 783-813.
    [111]J. L. Steger and R. F. Warming, Flux Vector Splitting of the Inviscid Gasdynamic Equations with Application to Finite Difference Methods, J. Comput.Phys., 40(1981), pp. 263-293.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700