Approximation of multi-scale elliptic problems using patches of finite elements

In this paper we present a method to solve numerically elliptic problems with multi-scale data using multiple levels of not necessarily nested grids. The method consists in calculating successive corrections to the solution in patches whose discretizations are not necessarily conforming. It resembles the FAC method (see Math. Comp. 46 (174) (1986) 439–456) and its convergence is obtained by a domain decomposition technique (see Math. Comp. 57 (195) (1991) 1–21). However it is of much more flexible use in comparison to the latter. To cite this article: R. Glowinski et al., C. R. Acad. Sci. Paris, Ser. I 337 (2003).     

Operator splitting for the numerical solution of free surface flow at low capillary numbers

Free surface flows are pervasive in engineering and biomedical applications. In many interesting cases—particularly when small length scales are involved—surface forces (capillarity) dominate the flow dynamics. In these cases, computing the flow together with the shape of the surfaces, requires specialized solution techniques. This article investigates the capabilities of an operator splitting/finite elements method at handling accurately incompressible viscous flow with free surfaces at low capillary numbers. The test case of flow in the downstream section of a slot coater is used for three reasons: (1) it is an established benchmark; (2) it represents an idealized, yet industrially relevant flow; (3) high-fidelity results obtained with monolithic algorithms are available in literature. The flow and free surface shape attained with the new operator splitting scheme agree very satisfactorily with the results obtained with monolithic solvers. Because of its inherent computational simplicity, the new operator splitting scheme is attractive for large-scale simulations, three-dimensional flows, and flows of complex fluids.     

Direct simulation of the motion of neutrally buoyant balls in a three-dimensional Poiseuille flow

In a previous article the authors introduced a Lagrange multiplier based fictitious domain method. Their goal in the present article is to apply a generalization of the above method to: (i) the numerical simulation of the motion of neutrally buoyant particles in a three-dimensional Poiseuille flow; (ii) study – via direct numerical simulations – the migration of neutrally buoyant balls in the tube Poiseuille flow of an incompressible Newtonian viscous fluid. Simulations made with one and several particles show that, as expected, the Segré–Silberberg effect takes place. To cite this article: T.-W. Pan, R. Glowinski, C. R. Mecanique 333 (2005).     

A kinematically coupled time-splitting scheme for fluid–structure interaction in blood flow

We present a new time-splitting scheme for the numerical simulation of fluid–structure interaction between blood flow and vascular walls. This scheme deals in a successful way with the problem of the added mass effect. The scheme is modular and it embodies the stability properties of implicit schemes at the low computational cost of loosely coupled ones.     

A penalty/Newton/conjugate gradient method for the solution of obstacle problems

Motivated by the search for non-negative solutions of a system of Eikonal equations with Dirichlet boundary conditions, we discuss in this Note a method for the numerical solution of parabolic variational inequality problems for convex sets such as $\text{K={v∣v∈H}{}_0{}^1\text{(Ω)}$, v?ψ a.e. on $\text{Ω}.}$ The numerical methodology combines penalty and Newton's method, the linearized problems being solved by a conjugate gradient algorithm requiring at each iteration the solution of a linear problem for a discrete analogue of the elliptic operator I?μΔ. Numerical experiments show that the resulting method has good convergence properties, even for small values of the penalty parameter. To cite this article: R. Glowinski et al., C. R. Acad. Sci. Paris, Ser. I 336 (2003).     

Models for capacitated lot-sizing problem with backlogging,setup carryover and crossover

Setup operations are significant in some production environments. It is mandatory that their production plans consider some features, as setup state conservation across periods through setup carryover and crossover. The modelling of setup crossover allows more flexible decisions and is essential for problems with long setup times. This paper proposes two models for the capacitated lot-sizing problem with backlogging and setup carryover and crossover. The first is in line with other models from the literature, whereas the second considers a disaggregated setup variable, which tracks the starting and completion times of the setup operation. This innovative approach permits a more compact formulation. Computational results show that the proposed models have outperformed other state-of-the-art formulation.     

A parallel finite element sliding mesh technique for the simulation of viscous flows in agitated tanks

A parallel sliding mesh algorithm for the finite element simulation of viscous fluid flows in agitated tanks is presented. Lagrange multipliers are used at the sliding interfaces to enforce the continuity between the fixed and moving subdomains. The novelty of the method consists of the coupled solution of the resulting velocity–pressure‐Lagrange multipliers system of equations by an ILU(0)‐QMR solver. A penalty parameter is introduced for both the interface and the incompressibility constraints to avoid pivoting problems in the ILU(0) algorithm. To handle the convective term, both the Newton–Raphson scheme and the semi‐implicit linearization are tested. A penalty parameter is introduced for both the interface and the incompressibility constraints to avoid the failure of the ILU(0) algorithm due to the lack of pivoting. Furthermore, this approach is versatile enough so that it allows partitioning of sliding and fixed subdomains if parallelization is required. Although the sliding mesh technique is fairly common in CFD, the main advantage of the proposed approach is its low computational cost due to the inexpensive and parallelizable calculations that involve preconditioned sparse iterative solvers. The method is validated for Couette and coaxial stirred tanks. Copyright © 2011 John Wiley & Sons, Ltd.