## Section: Overall Objectives

### Overall Objectives

The overall objectives of the NACHOS project-team are the formulation, analysis and evaluation of numerical methods and high performance algorithms for the solution of first order linear systems of partial differential equations (PDEs) with variable coefficients pertaining to electrodynamics and elastodynamics. In both domains, the applications targeted by the team involve the interaction of the underlying physical fields with media exhibiting space and time heterogeneities. Moreover, in most of the situations of practical relevance, the computational domain is irregularly shaped or/and it includes geometrical singularities. Both the heterogeneity and the complex geometrical features of the underlying media motivate the use of numerical methods working on non-uniform discretizations of the computational domain. In this context, the research efforts of the team aim at the development of unstructured (or hybrid structured/unstructured) mesh based methods with activities ranging from the mathematical analysis of numerical methods for the solution of the systems of PDEs of electrodynamics and elastodynamics, to the development of prototype 3d simulation software that efficiently exploits the capabilities of modern high performance computing platforms.

In the case of electrodynamics, the mathematical model of interest is the full system of unsteady Maxwell equations [50] which is a first-order hyperbolic linear system of PDEs (if the underlying propagation media is assumed to be linear). This system can be numerically solved using so-called time-domain methods among which the Finite Difference Time-Domain (FDTD) method introduced by K.S. Yee [55] in 1996 is the most popular and which often serves as a reference method for the works of the team. In the vast majority of existing time-domain methods, time advancing relies on an explicit time scheme. For certain types of problems, a time-harmonic evolution can be assumed leading to the formulation of the frequency-domain Maxwell equations whose numerical resolution requires the solution of a linear system of equations (i.e in that case, the numerical method is naturally implicit). Heterogeneity of the propagation media is taken into account in the Maxwell equations through the electrical permittivity, the magnetic permeability and the electric conductivity coefficients. In the general case, the electrical permittivity and the magnetic permeability are tensors whose entries depend on space (i.e heterogeneity in space) and frequency. In the latter case, the time-domain numerical modeling of such materials requires specific techniques in order to switch from the frequency evolution of the electromagnetic coefficients to a time dependency. Moreover, there exist several mathematical models for the frequency evolution of these coefficients (Debye model, Drude model, Drude-Lorentz model, etc.).

In the case of elastodynamics, the mathematical model of interest is the system of elastodynamic equations [45] for which several formulations can be considered such as the velocity-stress system. For this system, as with Yee's scheme for time-domain electromagnetics, one of the most popular numerical method is the finite difference method proposed by J. Virieux [54] in 1986. Heterogeneity of the propagation media is taken into account in the elastodynamic equations through the Lamé and mass density coefficients. A frequency dependence of the Lamé coefficients allows to take into account physical attenuation of the wave fields and characterizes a viscoelastic material. Again, several mathematical models are available for expressing the frequency evolution of the Lamé coefficients.