Numerical analysis of the Navier-Stokes equations
This paper discusses some conceptional questions of the numerical simulation of viscous incompressible flow which are related to the presence of boundaries.
This paper discusses some conceptional questions of the numerical simulation of viscous incompressible flow which are related to the presence of boundaries.
In this paper the numerical solution of two dimensional fluid-structure interaction problem is addressed. The fluid motion is modelled by the incompressible unsteady Navier-Stokes equations. The spatial discretization by stabilized finite element method is used. The motion of the computational domain is treated with the aid of Arbitrary Lagrangian Eulerian (ALE) method. The time-space problem is solved with the aid of multigrid method. The method is applied onto a problem of interaction of channel...
Numerical schemes are presented for a class of fourth order diffusion problems. These problems arise in lubrication theory for thin films of viscous fluids on surfaces. The equations being in general fourth order degenerate parabolic, additional singular terms of second order may occur to model effects of gravity, molecular interactions or thermocapillarity. Furthermore, we incorporate nonlinear surface tension terms. Finally, in the case of a thin film flow driven by a surface active agent (surfactant),...
Modeling the movement of a rigid particle in viscous fluid is a problem physicists and mathematicians have tried to solve since the beginning of this century. A general model for an ellipsoidal particle was first published by Jeffery in the twenties. We exploit the fact that Jeffery was concerned with formulae which can be used to compute numerically the velocity field in the neighborhood of the particle during his derivation of equations of motion of the particle. This is our principal contribution...
We introduce in this paper some elements for the mathematical and numerical analysis of algebraic turbulence models for oceanic surface mixing layers. In these models the turbulent diffusions are parameterized by means of the gradient Richardson number, that measures the balance between stabilizing buoyancy forces and destabilizing shearing forces. We analyze the existence and linear exponential asymptotic stability of continuous and discrete equilibria states. We also analyze the well-posedness...
This work presents the numerical solution of laminar incompressible viscous flow in a three dimensional branching channel with circular cross section for generalized Newtonian fluids. This model can be generalized by cross model in shear thinning meaning. The governing system of equations is based on the system of balance laws for mass and momentum. Numerical tests are performed on a three dimensional geometry, the branching channel with one entrance and two outlet parts. Numerical solution of the...
An algorithm for approximation of an unsteady fluid-structure interaction problem is proposed. The fluid is governed by the Navier-Stokes equations with boundary conditions on pressure, while for the structure a particular plate model is used. The algorithm is based on the modal decomposition and the Newmark Method for the structure and on the Arbitrary lagrangian Eulerian coordinates and the Finite Element Method for the fluid. In this paper, the continuity of the stresses at the interface was...
An algorithm for approximation of an unsteady fluid-structure interaction problem is proposed. The fluid is governed by the Navier-Stokes equations with boundary conditions on pressure, while for the structure a particular plate model is used. The algorithm is based on the modal decomposition and the Newmark Method for the structure and on the Arbitrary Lagrangian Eulerian coordinates and the Finite Element Method for the fluid. In this paper, the continuity of the stresses at the interface...
We propose a model for a medical device, called a stent, designed for the treatment of cerebral aneurysms. The stent consists of a grid, immersed in the blood flow and located at the inlet of the aneurysm. It aims at promoting a clot within the aneurysm. The blood flow is modelled by the incompressible Navier-Stokes equations and the stent by a dissipative surface term. We propose a stabilized finite element method for this model and we analyse its convergence in the case of the Stokes...
This paper focuses on the mathematical modelling and the numerical approximation of the flow of two immiscible incompressible fluids. The surface tension effects are taken into account and mixed boundary conditions are used. The weak formulation is introduced, discretized in time, and the finite element method is applied. The free surface motion is treated with the aid of the level set method. The numerical results are shown.
This work deals with the numerical solution of generalized Newtonian and Oldroyd-B fluids flow. The governing system of equations is based on the system of balance laws for mass and momentum for incompressible laminar viscous and viscoelastic fluids. Two different definition of the stress tensor are considered. For viscous case Newtonian model is used. For the viscoelastic case Oldroyd-B model is tested. Both presented models can be generalized. In this case the viscosity is defined as a shear rate...
The aim of this work is to present numerical results of non-Newtonian fluid flow in a model of bypass. Different angle of a connection between narrowed channel and the bypass graft is considered. Several rheology viscosity models were used for the non-Newtonian fluid, namely the modified Cross model and the Carreau-Yasuda model. The results of non-Newtonian fluid flow are compared to the results of Newtonian fluid. The fundamental system of equations is the generalized system of Navier-Stokes equations...
We propose here a model and a numerical scheme to compute the motion of rigid particles interacting through the lubrication force. In the case of a particle approaching a plane, we propose an algorithm and prove its convergence towards the solutions to the gluey particle model described in [B. Maury, ESAIM: Proceedings 18 (2007) 133–142]. We propose a multi-particle version of this gluey model which is based on the projection of the velocities onto a set of admissible velocities. Then, we describe...
We propose here a model and a numerical scheme to compute the motion of rigid particles interacting through the lubrication force. In the case of a particle approaching a plane, we propose an algorithm and prove its convergence towards the solutions to the gluey particle model described in [B. Maury, ESAIM: Proceedings18 (2007) 133–142]. We propose a multi-particle version of this gluey model which is based on the projection of the velocities onto a set of admissible velocities. Then, we describe...
Flow of particles suspended in a fluid can be found in numerous industrial processes utilizing sedimentation, fluidization and lubricated transport such as food processing, catalytic processing, slurries, coating, paper manufacturing, particle injection molding and filter operation. The ability to understand rheology effects of particulate flows is elementary for the design, operation and efficiency of the underlying processes. Despite the fact that particle technology is widely used, it is still...
In this work we present new numerical methods to simulate the mechanics of head-tape magnetic storage devices. The elastohydrodynamic problem is formulated in terms of a coupled system which is governed by a nonlinear compressible Reynolds equation for the air pressure over the head, and a rod model for the tape displacement. A fixed point algorithm between the solutions of the elastic and hydrodynamic problems is proposed. For the nonlinear Reynolds equation, a characteristics method and a...
In this contribution, we present the problem of shape optimization of the plunger cooling which comes from the forming process in the glass industry. We look for a shape of the inner surface of the insulation barrier located in the plunger cavity so as to achieve a constant predetermined temperature on the outward surface of the plunger. A rotationally symmetric system, composed of the mould, the glass piece, the plunger, the insulation barrier and the plunger cavity, is considered. The state problem...