Superconvergence of discontinuous Galerkin finite element method for the stationary Navier‐Stokes equations

This article focuses on discontinuous Galerkin method for the two‐ or three‐dimensional stationary incompressible Navier‐Stokes equations. The velocity field is approximated by discontinuous locally solenoidal finite element, and the pressure is approximated by the standard conforming finite element. Then, superconvergence of nonconforming finite element approximations is applied by using least‐squares surface fitting for the stationary Navier‐Stokes equations. The method ameliorates the two noticeable disadvantages about the given finite element pair. Finally, the superconvergence result is provided under some regular assumptions. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 421–436, 2007     

On weak–strong uniqueness for the compressible Navier–Stokes system with non-monotone pressure law

We show the weak–strong uniqueness property for the compressible Navier–Stokes system with general non-monotone pressure law. A weak solution coincides with the strong solution emanating from the same initial data as long as the latter solution exists.     

The Spin-Coating Process: Analysis of the Free Boundary Value Problem

In this paper, an accurate model of the spin-coating process is presented and investigated from the analytical point of view. More precisely, the spin-coating process is being described as a one-phase free boundary value problem for Newtonian fluids subject to surface tension and rotational effects. It is proved that for T > 0 there exists a unique, strong solution to this problem in (0, T) belonging to a certain regularity class provided the data and the speed of rotation are small enough in suitable norms. The strategy of the proof is based on a transformation of the free boundary value problem to a quasilinear evolution equation on a fixed domain. The keypoint for solving the latter equation is the so-called maximal regularity approach. In order to pursue in this direction one needs to determine the precise regularity classes for the associated inhomogeneous linearized equations. This is being achieved by applying the Newton polygon method to the boundary symbol.     

Application of the finite point method to high‐Reynolds number compressible flow problems

In this work, the finite point method is applied to the solution of high‐Reynolds compressible viscous flows. The aim is to explore this important field of applications focusing on two main aspects: the easiness and automation of the meshless discretization of viscous layers and the construction of a robust numerical approximation in the highly stretched clouds of points resulting in such domain areas. The flow solution scheme adopts an upwind‐biased scheme to solve the averaged Navier–Stokes equations in conjunction with an algebraic turbulence model. The numerical applications presented involve different attached boundary layer flows and are intended to show the performance of the numerical technique. The results obtained are satisfactory and indicative of the possibilities to extend the present meshless technique to more complex flow problems. Copyright © 2014 John Wiley & Sons, Ltd.     

Stability of algebraic multigrid for Stokes problems

An investigation is made of the performance of algebraic multigrid (AMG) solvers for the discrete Stokes problem. The saddle‐point formulations are based on the direct enforcement of the fundamental conservation laws in discrete spaces and subsequently stabilised with the aid of a regular splitting of the diffusion operator. AMG solvers based on an independent coarsening of the fields (the unknown approach) and also on a common coarsening (the point approach) are investigated. Both mixed‐order and equal‐order interpolations are considered. The dependence of convergence on the ‘degree of coarsening’ is investigated by studying the ‘convergence versus coarsening’ characteristics and their variation with mesh resolution. They show a consistency in shape, which reveals two distinct performance zones, one convergent the other divergent. The transition from the convergent to the divergent zones is discontinuous and occurs at a critical coarsening factor that is largely mesh independent. It signals a breakdown in the stability of the smoothing at the coarser levels of coarse grid approximation. It is shown that the previously observed, mesh‐dependent, scaling of convergence factors, which had suggested inconsistencies in the coarse grid approximation, is not a reliable marker of inconsistency. It is an indirect consequence of the breakdown in the stability of smoothing. For stable smoothing, reduction factors are shown to be largely mesh independent. The ability of mixed‐order interpolation to permit stable smoothing and therefore to deliver mesh‐independent convergence is explained. Two expedient options are suggested for obtaining mesh‐independent convergence for those AMG codes that are based on an equal‐order interpolation. Copyright © 2012 John Wiley & Sons, Ltd.     

An operator splitting strategy for fluid–structure interaction problems with thin elastic structures in an incompressible Newtonian flow

We present a computational framework based on the use of the Newton and level set methods to model fluid–structure interaction problems involving elastic membranes freely suspended in an incompressible Newtonian flow. The Mooney–Rivlin constitutive model is used to model the structure. We consider an extension to a more general case of the method described in Laadhari (2017) to model the elasticity of the membrane. We develop a predictor–corrector finite element method where an operator splitting scheme separates different physical phenomena. The method features an affordable computational burden with respect to the fully implicit methods. An exact Newton method is described to solve the problem, and the quadratic convergence is numerically achieved. Sample numerical examples are reported and illustrate the accuracy and robustness of the method.     

Application of Reynolds stress transport turbulence closure models to flows affected by Lorentz and buoyancy forces

The effect of magnetic field strength and orientation on two types of electromagnetically influenced turbulent flows was studied numerically under the Reynolds averaged Navier–Stokes (RANS) framework. Previous work (Wilson et al., 2014) used an electromagnetically extended linear eddy-viscosity model, whilst the current paper focuses on the performance of a more advanced Reynolds stress transport type model both with and without electromagnetic modifications proposed by Kenjereš et al. (2004). First, a fully-developed 2D channel flow is considered with a magnetic field imposed in either the wall-normal or streamwise direction. Both forms of the RSM gave good agreement with the DNS data for the wall-normal magnetic field across the range of Hartmann numbers with the additional electromagnetic terms providing a small, but noticeable, difference. For the streamwise magnetic field, where electromagnetic influence is only through the turbulence, the electromagnetically extended RSM performed well at moderate Hartmann numbers but returned laminar flow at the highest Hartmann number considered, contrary to the DNS. The RSM results were, however, significantly better than the previous eddy-viscosity model predictions. The second case is that of unsteady 3D Rayleigh–Bénard convection with a magnetic field imposed in either a horizontal or vertical direction. Results revealed that a significant reorganization of the flow structures is predicted to occur. For a vertically oriented magnetic field, the plume structures increase in number and become thinner and elongated along the magnetic field lines, leading to an increase in thermal mixing within the core in agreement with Hanjalić and Kenjereš (2000). With a horizontal magnetic field, the structures become two-dimensional and a striking realignment of the roll cells’ axes with the magnetic field lines occurs. The results demonstrate the capability of the Reynolds stress transport approach in modelling MHD flows that are relevant to industry and offer potential for those wishing to control levels of turbulence, heat transfer or concentration without recourse to mechanical means.     

Existence of global weak solutions for Navier‐Stokes‐Poisson equations with quantum effect and convergence to incompressible Navier‐Stokes equations

In this paper, we consider a three dimensional quantum Navier‐Stokes‐Poisson equations. Existence of global weak solutions is obtained, and convergence toward the classical solution of the incompressible Navier‐Stokes equation is rigorously proven for well prepared initial data. Furthermore, the associated convergence rates are also obtained. Copyright © 2014 John Wiley & Sons, Ltd.