Yuri V. Vassilevski

Alexey Y. Chernyshenko, Maxim A. Olshanskii, Yuri V. Vassilevski

The paper develops a hybrid method for solving a system of advection--diffusion equations in a bulk domain coupled to advection--diffusion equations on an embedded surface. A monotone nonlinear finite volume method for equations posed in the bulk is combined with a trace finite element method for equations posed on the surface. In our approach, the surface is not fitted by the mesh and is allowed to cut through the background mesh in an arbitrary way. Moreover, a triangulation of the surface into regular shaped elements is not required. The background mesh is an octree grid with cubic cells. As an example of an application, we consider the modeling of contaminant transport in fractured porous media. One standard model leads to a coupled system of advection--diffusion equations in a bulk (matrix) and along a surface (fracture). A series of numerical experiments with both steady and unsteady problems and different embedded geometries illustrate the numerical properties of the hybrid approach. The method demonstrates great flexibility in handling curvilinear or branching lower dimensional embedded structures.

Alexander Lozovskiy, Maxim A. Olshanskii, Yuri V. Vassilevski

The paper develops a finite element method for the Navier-Stokes equations of incompressible viscous fluid in a time-dependent domain. The method builds on a quasi-Lagrangian formulation of the problem. The paper provides stability and convergence analysis of the fully discrete (finite-difference in time and finite-element in space) method. The analysis does not assume any CFL time-step restriction, it rather needs mild conditions of the form $Δt\le C$, where $C$ depends only on problem data, and $h^{2m_u+2}\le c\,Δt$, $m_u$ is polynomial degree of velocity finite element space. Both conditions result from a numerical treatment of practically important non-homogeneous boundary conditions. The theoretically predicted convergence rate is confirmed by a set of numerical experiments. Further we apply the method to simulate a flow in a simplified model of the left ventricle of a human heart, where the ventricle wall dynamics is reconstructed from a sequence of contrast enhanced Computed Tomography images.

Kirill D. Nikitin, Maxim A. Olshanskii, Kirill M. Terekhov et al.

The paper develops a method for the numerical simulation of a free-surface flow of incompressible viscous fluid around a streamlined body. The body is a rigid stationary construction partially submerged in the fluid. The application we are interested in the paper is a flow around a surface mounted offshore oil platform. The numerical method builds on a hybrid finite volume / finite difference discretization using adaptive octree cubic meshes. The mesh is dynamically refined towards the free surface and the construction. Special care is taken to devise a discretization for the case of curvilinear boundaries and interfaces immersed in the octree Cartesian background computational mesh. To demonstrate the accuracy of the method, we show the results for two benchmark problems: the sloshing 3D container and the channel laminar flow passing the 3D cylinder of circular cross-section. Further, we simulate numerically a flow with surface waves around an offshore oil platform for the realistic set of geophysical data.