A Stable Unstructured Finite Volume Method with Multigrid for Parallel Large-Scale Incompressible Viscous Fluid Flow Computations

Author(s):  
Mehmet Sahin
2020 ◽  
Vol 50 (3) ◽  
pp. 287-302
Author(s):  
Róbert ČUNDERLÍK ◽  
Matej MEDĽA ◽  
Karol MIKULA

The paper presents local quasigeoid modelling in Slovakia using the finite volume method (FVM). FVM is used to solve numerically the fixed gravimetric boundary value problem (FGBVP) on a 3D unstructured mesh created above the real Earth's surface. Terrestrial gravimetric measurements as input data represent the oblique derivative boundary conditions on the Earth's topography. To handle such oblique derivative problem, its tangential components are considered as surface advection terms regularized by a surface diffusion. The FVM numerical solution is fixed to the GOCE-based satellite-only geopotential model on the upper boundary at the altitude of 230 km. The horizontal resolution of the 3D computational domain is 0.002 × 0.002 deg and its discretization in the radial direction is changing with altitude. The created unstructured 3D mesh of finite volumes consists of 454,577,577 unknowns. The FVM numerical solution of FGBVP on such a detailed mesh leads to large-scale parallel computations requiring 245 GB of internal memory. It results in the disturbing potential obtained in the whole 3D computational domain. Its values on the discretized Earth's surface are transformed into the local quasigeoid model that is tested at 404 GNSS/levelling benchmarks. The standard deviation of residuals is 2.8 cm and decreases to 2.6 cm after removing 9 identified outliers. It indicates high accuracy of the obtained FVM-based local quasigeoid model in Slovakia.


REAKTOR ◽  
2017 ◽  
Vol 5 (1) ◽  
pp. 1
Author(s):  
L. Buchori ◽  
Y. Bindar ◽  
D. Sasongko ◽  
IGBN Makertihartha

Generally, the momentum equation of fluid flow in porous media was solved by neglecting the terms of diffusion and convection such as Ergun, Darcy, Brinkman and Forchheimer models. Their model primarily applied for laminar flow. It is true that these model are limited to condition whether the models can be applied. Analytical solution for the model type above is available only for simple one-dimensional cases. For two or three-dimentional problem, numerical solution is the only solution. This work advances the flow model in porous media and provide two-dimentional flow field solution in porous media, which includes the diffusion and convection terms. The momentum lost due to flow and porous material interaction is modeled using the available  Brinkman-Forchheimer equation. The numerical method to be used is finite volume method. This method is suitable for the characteristic of fluid  flow in porous media which is averaged by a volume base. The effect of the solid and fluid interaction in porous  media is the basic principle of the flow model in morous media. The Brinkman-Forchheimer consider the momentum lost term to be determined by a quadratic function of the velocity component. The momentum and the continuity equation are solved for two-dimentional cylindrical coordinat . the result were validated with the experimental data. The velocity of the porous media was treated to be radially oscillated. The result of velocity profile inside packing show a good agreement in their trend with the Stephenson and Steward experimental data. The local superficial  velocity attains its global maximum and minimum at distances near 0.201 and 0.57 particle diameter, dp. velocity profile below packing was simulated. The result were validated with Schwartz and Smith experimental data. The result also show an excellent agreement with those experimental data.Keywords : finite volume method, porous media, flow distribution, velocity profile


Sign in / Sign up

Export Citation Format

Share Document