MATHICSE Technical Report : Multiscale methods for wave problems in heterogeneous media

In this paper we give a survey on various multiscale methods for the numerical solution of second order hyperbolic equations in highly heterogenous media. We concentrate on the wave equation and distinguish between two classes of applications. First we discuss numerical methods for the wave equation in heterogeneous media without scale separation. Such a setting is for instance encountered in the geosciences, where natural structures often exhibit a continuum of different scales, that all need to be resolved numerically to get meaningful approximations. Approaches tailored for these settings typically involve the construction of generalized finite element spaces, where the basis functions incorporate information about the data variations. In the second part of the paper, we discuss numerical methods for the case of structured media with scale separation. This setting is for instance encountered in engineering sciences, where materials are often artificially designed. If this is the case, the structure and the scale separation can be explicitly exploited to compute appropriate homogenized/upscaled wave models that only exhibit a single coarse scaleand that can be hence solved at significantly reduced computational costs.

Multiscale methods for wave problems in heterogenous media

Assyr Abdulle, Patrick Henning

In this paper we give a survey on various multiscale methods for the numerical solution of second order hyperbolic equations in highly heterogenous media. We concentrate on the wave equation and distinguish between two classes of applications. First we discuss numerical methods for the wave equation in heterogeneous media without scale separation. Such a setting is for instance encountered in the geosciences, where natural structures often exhibit a continuum of different scales, that all need to be resolved numerically to get meaningful approximations. Approaches tailored for these settings typically involve the construction of generalized finite element spaces, where the basis functions incorporate information about the data variations. In the second part of the paper, we discuss numerical methods for the case of structured media with scale separation. This setting is for instance encountered in engineering sciences, where materials are often artificially designed. If this is the case, the structure and the scale separation can be explicitly exploited to compute appropriate homogenized/upscaled wave models that only exhibit a single coarse scale and that can be hence solved at significantly reduced computational costs.

2017

Analysis and numerical simulation of viscoelastic flows

Andrea Bonito

Mathematical and numerical aspects of viscoelastic flows are investigated here. Two simplified mathematical models are considered. They are motivated by a splitting algorithm for solving viscoelastic flows with free surfaces. The first model is a simplified Oldroyd-B model. Existence on a fixed time interval is proved in several Banach spaces provided the data are small enough. Short time existence is also proved for arbitrarily large data in Hölder spaces for the time variable. These results are based on the maximal regularity property of the Stokes operator and on the analycity behavior of the corresponding semi-group. A finite element discretization in space is then proposed. Existence of the numerical solution is proved for small data, as well as a priori error estimates, using an implicit function theorem framework. Then, the extension of these results to a stochastic simplified Hookean dumbbells model is discussed. Because of the presence of the Brownian motion, existence in a fixed time interval, provided the data are small enough, is proved only in some of the Banach spaces considered previously. The dumbbells' elongation is split in two parts, one satisfying a standart stochastic differential equation, the other satisfying a partial differential equation with a stochastic source term. A finite element discretization in space is also proposed. Existence of the numerical solution is proved for small data, as well as a priori error estimates. A numerical algorithm for solving viscoelastic flows with free surfaces is also described. This algorithm is based on a splitting method in time and two different meshes are used for the space discretization. Convergence of the numerical model is checked for the pure extensional flow and the filling of a pipe. Then, numerical results are reported for the stretching of a filament and for jet buckling.

EPFL2006

Modélisation mathématique et simulation numérique des phénomènes dynamiques et thermiques apparaissant dans un glacier

Pierre Kirner

In this work, two problems linked to glacier modeling are investigated. We propose an optimisation method for studying the flow of the ice and we present a numerical study about glacier thermal phenomena. In the first chapter of this thesis, we expose the models of these two problems. On one hand, we note that the boundary conditions on the bedrock are misunderstood, which explains it is difficult to obtain an accurate simulation of the motion of the ice. Also we establish a mathematical model where the bedrock boundary conditions depend on a control parameter. The aim of this study is to minimize a cost functional describing the difference between the computed velocity at the surface and the measure done. We study the cost functional with respect to the control parameter and we detail an optimisation method to solve the optimal control problem. On the other hand, we introduce two thermodynamical model governing the temperature and the water content field. The models correspond to a Stefan problem for the temperature and a convection-diffusion equation for the water content. The second chapter deals with the numerical resolution of the optimisation problem. First, a Finite Element Method (FEM) is described to solve the partial differential equations. Then, the algorithms used for the optimal control problem are detailed. Finally, this techniques are applied on two glaciers : Griesgletcher for 2D and Storglaciaren for 3D. The third chapter deals with the numerical resolution of the temperature and the water content models. A FEM is used for each problem. Concerning the temperature problem, the Stefan problem is numerically solved and the results allow to detect a free surface between the temperated ice and the cold ice. The water content field is also simulated. Numerical results are discussed on the Storglaciaren.

EPFL2007