A 3D wavelet-based incompressible Navier-Stokes solver for fully adaptive computations in time-varying geometries

12/11/2019
by   Thomas Engels, et al.
0

We present a wavelet-based adaptive method for computing 3D flows in complex, time-dependent geometries, implemented on massively parallel computers. The incompressible fluid is modeled with an artificial compressibility approach in order to avoid the solution of elliptical problems. No-slip and in/outflow boundary conditions are imposed using volume penalization. The governing equations are discretized on a locally uniform Cartesian grid with centered finite differences, and integrated in time with a Runge–Kutta scheme, both of 4th order. The domain is partitioned into cubic blocks with equidistant grids and, for each block, biorthogonal interpolating wavelets are used as refinement indicators and prediction operators. Thresholding of wavelet coefficients allows to introduce dynamically evolving grids and the adaption strategy tracks the solution in both space and scale. Blocks are distributed among MPI processes and the global topology of the grid is encoded in a tree-like data structure. Analyzing the different physical and numerical parameters allows balancing their individual error contributions and thus ensures optimal convergence while minimizing computational effort. Different validation tests score accuracy and performance of our new open source code, WABBIT (Wavelet Adaptive Block-Based solver for Interactions with Turbulence), on massively parallel computers using fully adaptive grids. Flow simulations of flapping insects demonstrate its applicability to complex, bio-inspired problems.

READ FULL TEXT

page 14

page 16

page 19

page 20

page 28

research
12/14/2021

MURPHY – A scalable multiresolution framework for scientific computing on 3D block-structured collocated grids

We present the derivation, implementation, and analysis of a multiresolu...
research
02/23/2021

Adaptive two- and three-dimensional multiresolution computations of resistive magnetohydrodynamics

Fully adaptive computations of the resistive magnetohydrodynamic (MHD) e...
research
11/21/2018

Multi-Block Grid deformation Method in 3D

Multi-block grids provide the computational efficiency of structured gri...
research
07/01/2018

Ensuring domain consistency in an adaptive framework with distributed topology for fluid flow simulations

Top-tier parallel computing clusters continue to accumulate more and mor...
research
06/16/2023

Stable nodal projection method on octree grids

We propose a novel collocated projection method for solving the incompre...
research
03/09/2017

ForestClaw: A parallel algorithm for patch-based adaptive mesh refinement on a forest of quadtrees

We describe a parallel, adaptive, multi-block algorithm for explicit int...
research
12/02/2019

An asynchronous incomplete block LU preconditioner for computational fluid dynamics on unstructured grids

We present a study of the effectiveness of asynchronous incomplete LU fa...

Please sign up or login with your details

Forgot password? Click here to reset