Journal of computational fluids engineering (한국전산유체공학회지)

Korean Society of Computational Fluids Engineering (한국전산유체공학회)
권호 표지

Numerical Simulation for the Advection Equation on the Sphere by Sphere-Lagrangian Method

Semi-Lagrangian법을 이용한 구 좌표계에서의 이류 방정식 해석

  • 윤성영 (영남대학교 기계공학과)
  • Published : 2004.09.01
Download PDF
⟨ Previous Next ⟩

Abstract

A Semi-Lagrangian method based on CIP(Cubic Interpolated Pseudoparticle)method is proposed and it is applied to solve the two dimensional advection equation. Especially the attentions are given to settle the pole problem and to enhance the accuracy in solving the advection equation on the spherical coordinate system. Tn this algorithm, the CU method is employed as the Semi-Lagrangian method and extended to the spherical coordinate system. To enhance the accuracy of the solution, the spatial discretization is made by CIP method. The mathematical formulation and numerical results are also described. To verify the efficiency, accuracy and capability of proposed algorithm, two dimensional rotating cosine bell problem and the frontogenesis problem are simulated by the present scheme. As results, it is confirmed that the present scheme gives an accurate solution and settles the pole problem in the advection equation on the sphere.

Keywords

References

  1. D.L.Williamson and P.J.Rasch, 'Two- Dimensional Semi-Lagrangian Transport with Shape-Preserving Interpolation,', Monthly Weather Review, 117, pp.102-129, 1989
  2. P.K.Smolarkiewicz and P.J.Rasch, 'Monotone Advection on the Sphere: An Euler Versus Semi-Lagrangian Approach,' Journal of the Atmospheric Sciences, 48, pp.793-810, 1991
  3. P.Rasch, 'Conservative shape-preseving two-dimensional transport on a spherical reduced grid,' Monthly Weather Review, 122, pp.1337-1350,1994
  4. Seong Y.Yoon, Y.Tanaka, F.Xiao, and T.Yabe 'Numerical simulation of shallow water equations using CIP Method with high order spatial resolution,' 2001 Workshop on the Solution of Partial Differential Equations on the Sphere, Montreal, Canada, p.19, May15-18, 2001
  5. T.Yabe and T.Aoki, 'A universal solver for hyperbolic equations by cubic-polynomial interpolation,' Computer Physics Communications, 66, pp.219-232,1991
  6. Seong Y.Yoon and T.Yabe 'The unified simulation for incompressible and compressible flow by the predictor- corrector scheme based on the CIP method', Computer Physics Communications, 119, pp.149-158, 1999
  7. T.Yabe and F.Xiao, 'The Constrained Interpolation Profile Method for Multiphase Analysis', Journal of Computational Physics, 169, pp.556-593, 2001
  8. H.Akima,'A new method of interpolation and smooth curve fitting based on local procedures,' J. Assoc. Comput. Mach., 17, pp.589-602,1970
  9. J.M.Hyman, 'Accurate monotonicity preserving cubic interpolations,' SIAM J.Sci.Stat.Comput., 5,pp.300-304,1984
  10. D.L.Williamson,J.B.Drake,J.J.Hack,R.Jakob and P.Swarztrauber, 'A standard test set for numerical approximations to the shallow water equations in spherical geometry,' Journal of computational physics,102, pp.211-224.1992
  11. C.A.Doswell III, 'A kinematic analysis of frontogenesis associated with a nondivergent flow,' Journal of the atmospheric sciences, 41, pp.1242-1248, 1984
  12. R.D.Nair, B.Machenhauer, 'The Mass Convervative Cell-Integrated Semi- Lagrangian Advection Scheme on the Sphere,' Monthly Weather Review, 130, pp.649-667, 2002
  13. R.D.Nair, J.S.Scroggs, and F.H.M.Semazzi, 'A forward-trajectory global semi- Lagrangian transport scheme,' Journal of Computational Physics, 190, pp.275-294, 2003
  14. R.D.Nair, J.S.Scroggs, and F.H.M.Semazzi, 'A forward-trajectory global semi- Lagrangian transport scheme,' Journal of Computational Physics, 54, pp.174-201, 1984