A vectorizable simulation method for the Boltzmann equation
In this work, we address the numerical solution of fluid-structure interaction problems. This issue is particularly difficulty to tackle when the fluid and the solid densities are of the same order, for instance as it happens in hemodynamic applications, since fully implicit coupling schemes are required to ensure stability of the resulting method. Thus, at each time step, we have to solve a highly non-linear coupled system, since the fluid domain depends on the unknown displacement of the structure....
In this work, we address the numerical solution of fluid-structure interaction problems. This issue is particularly difficulty to tackle when the fluid and the solid densities are of the same order, for instance as it happens in hemodynamic applications, since fully implicit coupling schemes are required to ensure stability of the resulting method. Thus, at each time step, we have to solve a highly non-linear coupled system, since the fluid domain depends on the unknown displacement of...
This paper provides an accelerated two-grid stabilized mixed finite element scheme for the Stokes eigenvalue problem based on the pressure projection. With the scheme, the solution of the Stokes eigenvalue problem on a fine grid is reduced to the solution of the Stokes eigenvalue problem on a much coarser grid and the solution of a linear algebraic system on the fine grid. By solving a slightly different linear problem on the fine grid, the new algorithm significantly improves the theoretical error...
The reduced basis method is a model reduction technique yielding substantial savings of computational time when a solution to a parametrized equation has to be computed for many values of the parameter. Certification of the approximation is possible by means of an a posteriori error bound. Under appropriate assumptions, this error bound is computed with an algorithm of complexity independent of the size of the full problem. In practice, the evaluation of the error bound can become very sensitive...
We present an alternative framework for designing efficient numerical schemes for non-conservative hyperbolic systems. This approach is based on the design of entropy conservative discretizations and suitable numerical diffusion operators that mimic the effect of underlying viscous mechanisms. This approach is illustrated by considering two model non-conservative systems: Lagrangian gas dynamics in non-conservative form and a form of isothermal Euler equations. Numerical experiments demonstrating...
We present an alternative framework for designing efficient numerical schemes for non-conservative hyperbolic systems. This approach is based on the design of entropy conservative discretizations and suitable numerical diffusion operators that mimic the effect of underlying viscous mechanisms. This approach is illustrated by considering two model non-conservative systems: Lagrangian gas dynamics in non-conservative form and a form of isothermal Euler equations. Numerical experiments demonstrating...
This work represents a first step towards the simulation of the motion of water in a complex hydrodynamic configuration, such as a channel network or a river delta, by means of a suitable “combination” of different mathematical models. In this framework a wide spectrum of space and time scales is involved due to the presence of physical phenomena of different nature. Ideally, moving from a hierarchy of hydrodynamic models, one should solve throughout the whole domain the most complex model (with...
Two-scale convergence is a special weak convergence used in homogenization theory. Besides the original definition by Nguetseng and Allaire two alternative definitions are introduced and compared. They enable us to weaken requirements on the admissibility of test functions . Properties and examples are added.
This paper presents an extension to stabilized methods of the standard technique for the numerical analysis of mixed methods. We prove that the stability of stabilized methods follows from an underlying discrete inf-sup condition, plus a uniform separation property between bubble and velocity finite element spaces. We apply the technique introduced to prove the stability of stabilized spectral element methods so as stabilized solution of the primitive equations of the ocean.
This paper presents an extension to stabilized methods of the standard technique for the numerical analysis of mixed methods. We prove that the stability of stabilized methods follows from an underlying discrete inf-sup condition, plus a uniform separation property between bubble and velocity finite element spaces. We apply the technique introduced to prove the sta bi li ty of stabilized spectral element methods so as stabilized solution of the primitive equations of the ocean.