A grid-free abstraction of the Navier-Stokes equations in Fortran 95/2003

Author:

Rouson Damian W. I.1,Rosenberg Robert1,Xu Xiaofeng2,Moulitsas Irene3,Kassinos Stavros C.3

Affiliation:

1. U.S. Naval Research Laboratory, Washington, DC

2. University of Maryland, College Park, Nova, MI

3. University of Cyprus, Nicosia, Cyprus

Abstract

Computational complexity theory inspires a grid-free abstraction of the Navier-Stokes equations in Fortran 95/2003. A novel complexity analysis estimates that structured programming time grows at least quadratically with the number of program lines. Further analysis demonstrates how an object-oriented strategy focused on mathematical objects renders the quadratic estimate scale-invariant, so the time required for the limiting factor in program development (debugging) no longer grows as the code grows. Compared to the coordinate-free C++ programming of Grant et al. [2000], grid-free Fortran programming eliminates a layer of procedure calls, eliminates a related need for the C++ template construct, and offers a shorter migration path for Fortran programmers. The grid-free strategy is demonstrated by constructing a physical-space driver for a Fourier-space Navier-Stokes solver. Separating the expression of the continuous mathematical model from the discrete numerics clarifies issues that are otherwise easily conflated. A run-time profile suggests that grid-free design substantially reduces the fraction of the procedures that significantly impact runtime, freeing more code to be structured in ways that reduce development time. Applying Amdahl's law to the total solution time (development time plus run time) leads to a strategy that negligibly impacts development time but achieves 58% of the maximum possible speedup.

Funder

Office of Naval Research

Sixth Framework Programme

Publisher

Association for Computing Machinery (ACM)

Subject

Applied Mathematics,Software

Reference37 articles.

1. Object-Oriented Programming via Fortran 90/95

2. Barton J. J. and Nackman L. R. 1994. Scientific and Engineering C++: An Introduction with Advanced Examples. Addison-Wesley. Boston MA. Barton J. J. and Nackman L. R. 1994. Scientific and Engineering C++: An Introduction with Advanced Examples. Addison-Wesley. Boston MA.

3. A Component Architecture for High-Performance Scientific Computing

Cited by 9 articles. 订阅此论文施引文献 订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献

1. Proving a core code for FDM correct by 2 + dw tests;Proceedings of the 5th ACM SIGPLAN International Workshop on Libraries, Languages, and Compilers for Array Programming;2018-06-19

2. High-Performance Design Patterns for Modern Fortran;Scientific Programming;2015

3. Dispersed-phase structure in sheared MHD turbulence;Journal of Turbulence;2012-01

4. Object construction and destruction design patterns in Fortran 2003;Procedia Computer Science;2010-05

5. Design patterns for multiphysics modeling in Fortran 2003 and C++;ACM Transactions on Mathematical Software;2010-01

同舟云学术

1.学者识别学者识别

2.学术分析学术分析

3.人才评估人才评估

"同舟云学术"是以全球学者为主线,采集、加工和组织学术论文而形成的新型学术文献查询和分析系统,可以对全球学者进行文献检索和人才价值评估。用户可以通过关注某些学科领域的顶尖人物而持续追踪该领域的学科进展和研究前沿。经过近期的数据扩容,当前同舟云学术共收录了国内外主流学术期刊6万余种,收集的期刊论文及会议论文总量共计约1.5亿篇,并以每天添加12000余篇中外论文的速度递增。我们也可以为用户提供个性化、定制化的学者数据。欢迎来电咨询!咨询电话:010-8811{复制后删除}0370

www.globalauthorid.com

TOP

Copyright © 2019-2024 北京同舟云网络信息技术有限公司
京公网安备11010802033243号  京ICP备18003416号-3