줄리아 언어를 이용한 고성능 해양모델의 개발 = Development of a High Performance Ocean Model using Julia Language

국문 초록 (Abstract)

고성능 해양모델을 개발하기 위하여 적시 컴파일(Just-In-Time) 언어인 줄리아 언어를 사용하였고, 운동량 방정식의 해를 구하기 위해 연속완화법으로 푸아송 방정식을 푸는 코드를 작성하였다. 다음으로, 줄리아 계산 코드를 시험하기 위하여 두 가지의 모델을 구축하였다. 첫 번째로, 일정한 유량의 생성/소멸(source/sink) 조건을 시험하기 위해 단순한 수로 형태를 모델링하였고, 두 번째로, 조석 외력(tidal forcing) 및 전향력(Coriolis force), 난류확산계수로 인한 효과 등을 시험하기 위해 황해(Yellow Sea)를 단순화하여 모델링하였다. 모델은 두 가지 시나리오 안에서 총 8개의 실험안으로 테스트되었다. 테스트 결과, 시나리오 1에서 3가지 실험안의 수심 평균된 유속은 이론 값에 완벽하게 수렴하였고, 해저마찰로 인한 수직적 유속 구배를 잘 보여주었다. 또한 시나리오 2에서는 황해의 무조점과 우리나라 서해 중부와 남부 연안의 조석 특성을 잘 재현하였고, 전향력과 수직 난류확산계수에 따른 결과의 차이를 잘 보여주었다. 따라서, 줄리아 언어를 이용한 해양모델을 개발하는 데에 성공하였다고 판단되며, 이는 해양모델이 고전적인 컴파일 언어에서 적시 컴파일 언어로 성공적으로 넘어가는 단계에 오게 됐다는 것을 시사한다.

다국어 초록 (Multilingual Abstract)

In order to develop a high performance ocean model, we used Julia, a Just-In-Time compile language, and to obtain the solution of the momentum equation, we made the code to solve the Poisson equation by the Successive Over-Relaxation method. And then we made two models to test Julia calculation codes. First, a simple channel form is modeled to test constant source/sink conditions. Second, the simplified Yellow Sea was modeled to test tidal forcing, Coriolis forces, and the effect of vertical eddy diffusivity coefficients. The model has been tested with a total of eight cases in the two scenarios. As a result of the test, the depth-averaged current speed of the three cases in Scenario 1 converged perfectly to the theoretical value, and that showed well a vertical flow velocity gradient due to the bottom friction. Also, the result of Scenario 2 represented well the amphidromic points of Yellow Sea and the tidal characteristics of mid-western and southwestern coast of Korea. Therefore, it is considered that the ocean model using Julia language has developed successfully, this suggests that the ocean model has come to the stage of successful transition from a classical compile language to a Just-In-Time compile language.

참고문헌 (Reference)

1 Väisälä, V., "Über die Wirkung der Windschwankungen aufdie Pilotbeobachtungen" 2 (2): 1-46, 1925

2 Kochergin, V. P., "Three- Dimensional Coastal Ocean Models, Coastal Estuarine Science Series, vol. 4" Union, American Geophysical 201-208, 1987

3 Brunt, D., "The period of simple vertical oscillations in the atmosphere" 53 (53): 30-32, 1927

4 JuliaLang, "The Julia Language"

5 OSU Tidal Data Inversion, "TPX09-ATLAS"

6 Poisson, S.-D., "Remarques sur une équation qui se présente dans a théorie des attractions des sphéroides" 3 : 388-392, 1813

7 Fringer, O. B., "Reducing numerical diffusion in interfacial gravity wave simulations" 49 : 301-329, 2005

8 McCarthy, J., "Recursive Functions of Symbolic Expressions and Their Computation by Machine, Part I" 3 (3): 184-195, 1960

9 The R Foundation, "R"

10 "Python Software Foundation"

11 Kämpf, J., "Ocean Modelling for Beginners (using open-source software)" Springer 59-, 2009

12 MathWorks, "MATLAB"

13 Lattner, C., "LLVM: An Infrastructure for Mulit-Stage Optimization" University of Illinois at Urbana-Champaign 2002

14 JuliaLang, "Julia 1.0 Documentation"

15 Aycock, "JIT Compilation Techniques, 2.1 Genesis"

16 GitHub, Inc, "GitHub"

17 Smagorinsky, J., "General circulation experiments with the primitive equations - I. The basic experiment" 91 : 99-165, 1963

18 Arakawa, A., "Computational Design of the Basic Dynamical Processes of the UCLA General Circulation Model" 17 : 173-265, 1977

19 Kämpf, J., "Advanced Ocean Modelling (using open-source software)" Springer 2010

20 Hamrick, J. M., "A three-dimensional environmental fluid dynamics computer code : Theoretical and computational aspects" The College of William and Mary, Virginia Institute of Marine Science 1-63, 1992