Comparison of Explicit Method of Solution for CFD Euler Problems using MATLAB® and FORTRAN 77
| Autor: | Hassan Abbas Khawaja, Anders Samuelsen Nordli |
|---|---|
| Jazyk: | angličtina |
| Rok vydání: | 2019 |
| Předmět: |
VDP::Teknologi: 500::Informasjons- og kommunikasjonsteknologi: 550::Datateknologi: 551
VDP::Technology: 500::Information and communication technology: 550::Computer technology: 551 Computer science Fortran Computational Mechanics 02 engineering and technology Computational fluid dynamics computer.software_genre 01 natural sciences 010305 fluids & plasmas Physics::Fluid Dynamics symbols.namesake 020401 chemical engineering ComputingMethodologies_SYMBOLICANDALGEBRAICMANIPULATION 0103 physical sciences Code (cryptography) 0204 chemical engineering MATLAB ComputingMilieux_MISCELLANEOUS ComputingMethodologies_COMPUTERGRAPHICS computer.programming_language Fluid Flow and Transfer Processes Numerical Analysis Programming language business.industry Explicit method lcsh:QC1-999 Mechanics of Materials Modeling and Simulation Euler's formula symbols Computer Science::Mathematical Software business computer lcsh:Physics |
| Zdroj: | International Journal of Multiphysics, Vol 13, Iss 2 (2019) Scopus-Elsevier |
| ISSN: | 1750-9548 |
| Popis: | Published version, available at: http://dx.doi.org/10.21152/1750-9548.13.2.203 This work presents a comparison of an explicit method of solution for an inviscid compressible fluid mechanics problem using Euler equations for two-dimensional internal flows. The same algorithm was implemented in both FORTRAN 77 and MATLAB®. The algorithm includes Runge‒Kutta time marching scheme with smoothing. Both solvers were initialized in the same manner. In addition, it was ensured that both solvers have the exact same values for time step, convergence criteria, boundary conditions, and the grid. The only difference between the two solvers was the precision of variables. The problem solved was a two-dimensional dual bump with an accelerating flow through a duct. The same algorithm solving the Euler equations of fluid flow is implemented in both FORTRAN 77 and MATLAB®, and applied to identical input. While the solutions look qualitativly the same, a 20% difference in the stationary solution is observed. No claim is made of the relevance of the computations to actual fluid flow, rather the key takeaway being that two finite and deterministic computations of the same algorithm on the same input in FORTRAN 77 and MATLAB® produce different output. |
| Databáze: | OpenAIRE |
| Externí odkaz: |
|