Divergence-free radial kernel for surface Stokes equations based on the surface Helmholtz decomposition

2020 ◽  
Vol 256 ◽  
pp. 107408
Author(s):  
Jingwei Li ◽  
Zhiming Gao ◽  
Zihuan Dai ◽  
Xinlong Feng
2020 ◽  
Vol 54 (1) ◽  
pp. 105-128 ◽  
Author(s):  
Hwi Lee ◽  
Qiang Du

Nonlocal gradient operators are prototypical nonlocal differential operators that are very important in the studies of nonlocal models. One of the simplest variational settings for such studies is the nonlocal Dirichlet energies wherein the energy densities are quadratic in the nonlocal gradients. There have been earlier studies to illuminate the link between the coercivity of the Dirichlet energies and the interaction strengths of radially symmetric kernels that constitute nonlocal gradient operators in the form of integral operators. In this work we adopt a different perspective and focus on nonlocal gradient operators with a non-spherical interaction neighborhood. We show that the truncation of the spherical interaction neighborhood to a half sphere helps making nonlocal gradient operators well-defined and the associated nonlocal Dirichlet energies coercive. These become possible, unlike the case with full spherical neighborhoods, without any extra assumption on the strengths of the kernels near the origin. We then present some applications of the nonlocal gradient operators with non-spherical interaction neighborhoods. These include nonlocal linear models in mechanics such as nonlocal isotropic linear elasticity and nonlocal Stokes equations, and a nonlocal extension of the Helmholtz decomposition.


2003 ◽  
Vol 70 (1) ◽  
pp. 44-49 ◽  
Author(s):  
V. Sarin ◽  
A. H. Sameh

The paper presents an algebraic scheme to construct hierarchical divergence-free basis for velocity in incompressible fluids. A reduced system of equations is solved in the corresponding subspace by an appropriate iterative method. The basis is constructed from the matrix representing the incompressibility constraints by computing algebraic decompositions of local constraint matrices. A recursive strategy leads to a hierarchical basis with desirable properties such as fast matrix-vector products, a well-conditioned reduced system, and efficient parallelization of the computation. The scheme has been extended to particulate flow problems in which the Navier-Stokes equations for fluid are coupled with equations of motion for rigid particles suspended in the fluid. Experimental results of particulate flow simulations have been reported for the SGI Origin 2000.


2010 ◽  
Vol 51 ◽  
pp. 654 ◽  
Author(s):  
Steven William Armfield ◽  
Nicholas Williamson ◽  
Michael Kirkpatrick ◽  
Robert Street

2018 ◽  
Vol 20 (03) ◽  
pp. 1750046 ◽  
Author(s):  
M. Schonbek ◽  
G. Seregin

In this note, we study the behavior of Lebesgue norms [Formula: see text] of solutions [Formula: see text] to the Cauchy problem for the Stokes system with drift [Formula: see text], which is supposed to be a divergence free smooth vector-valued function satisfying a scale invariant condition.


2018 ◽  
Vol 18 (3) ◽  
pp. 353-372 ◽  
Author(s):  
Naveed Ahmed ◽  
Alexander Linke ◽  
Christian Merdon

AbstractIn this contribution, we review classical mixed methods for the incompressible Navier–Stokes equations that relax the divergence constraint and are discretely inf-sup stable. Though the relaxation of the divergence constraint was claimed to be harmless since the beginning of the 1970s, Poisson locking is just replaced by another more subtle kind of locking phenomenon, which is sometimes called poor mass conservation and led in computational practice to the exclusion of mixed methods with low-order pressure approximations like the Bernardi–Raugel or the Crouzeix–Raviart finite element methods. Indeed, divergence-free mixed methods and classical mixed methods behave qualitatively in a different way: divergence-free mixed methods are pressure-robust, which means that, e.g., their velocity error is independent of the continuous pressure. The lack of pressure robustness in classical mixed methods can be traced back to a consistency error of an appropriately defined discrete Helmholtz projector. Numerical analysis and numerical examples reveal that really locking-free mixed methods must be discretely inf-sup stable and pressure-robust, simultaneously. Further, a recent discovery shows that locking-free, pressure-robust mixed methods do not have to be divergence free. Indeed, relaxing the divergence constraint in the velocity trial functions is harmless, if the relaxation of the divergence constraint in some velocity test functions is repaired, accordingly. Thus, inf-sup stable, pressure-robust mixed methods will potentially allow in future to reduce the approximation order of the discretizations used in computational practice, without compromising the accuracy.


Sign in / Sign up

Export Citation Format

Share Document