With the rapidly growing computational power of modern computers, physically based rendering has found its way into real world applications. Real-time simulations and renderings of fire and smoke had become one major research interest in modern video game industry, and will continue being one important research direction in computer graphics. To visually recreate realistic dynamic fire and smoke is a complicated problem. Furthermore, to solve the problem requires knowledge from various areas, ranged from computer graphics and image processing to computational physics and chemistry. Even though most of the areas are well-studied separately, when combined, new challenges will emerge. This thesis focuses on three aspects of the problem, dynamic, real-time and realism, to propose a solution in form of a GPGPU pipeline, along with its implementation. Three main areas with application in the problem are discussed in detail: fluid simulation, volumetric radiance estimation and volumetric rendering. The weights are laid upon the first two areas. The results are evaluated around the three aspects, with graphical demonstrations and performance measurements. Uniform grids are used with Finite Difference (FD) discretization scheme to simplify the computation. FD schemes are easy to implement in parallel, especially with ComputeShader, which is well supported in Unity engine. The whole implementation can easily be integrated into any real-world applications in Unity or other game engines that support DirectX 11 or higher.

Dans cette thèse, tout d’abord, nous faisons l’Analyse Mathématique du modèle exact du chauffage radiatif d’un corps semi-transparent $\Omega$ par une source radiative noire qui l’entoure. Il s’agit donc d’étudier le couplage d’un système d’Equations de Transfert Radiatif avec condition au bord de réflectivité indépendantes avec une équation de conduction de la chaleur non linéaire avec condition limite non linéaire de type Robin. Nous prouvons l’existence et l’unicité de la solution et nous démontrons des bornes uniformes sur la solution et les intensités radiatives dans chaque bande de longueurs d’ondes pour laquelle le corps est semi-transparent, en fonction de bornes sur les données, Deuxièmement, nous considérons le problème du contrôle optimal de la température absolue à l’intérieur du corps semi-transparent $\Omega$ en agissant sur la température absolue de la source radiative noire qui l’entoure. À cet égard, nous introduisons la fonctionnelle coût appropriée et l’ensemble des contrôles admissibles $T_{S}$, pour lesquels nous prouvons l’existence de contrôles optimaux. En introduisant l’espace des états et l’équation d’état, une condition nécessaire de premier ordre pour qu’un contrôle $T_{S}$ : t ! $T_{S}$ (t) soit optimal, est alors dérivée sous la forme d’une inéquation variationnelle en utilisant le théorème des fonctions implicites et le problème adjoint. Ensuite, nous considérons le problème de l’existence et de l’unicité d’une solution faible des équations de la thermoviscoélasticité dans une formulation mixte de type Hellinger- Reissner, la nouveauté par rapport au travail de M.E. Rognes et R. Winther (M3AS, 2010) étant ici l’apparition de la viscosité dans certains coefficients de l’équation constitutive, viscosité qui dépend dans ce contexte de la température absolue T(x, t) et donc en particulier du temps t. Enfin, nous considérons dans ce cadre le problème du contrôle optimal de la déformation du corps semi-transparent $\Omega$, en agissant sur la température absolue de la source radiative noire qui l’entoure. Nous prouvons l’existence d’un contrôle optimal et nous calculons la dérivée Fréchet de la fonctionnelle coût réduite. / This thesis begins with a rigorous mathematical analysis of the radiative heating of a semi-transparent body made of glass, by a black radiative source surrounding it. This requires the study of the coupling between quasi-steady radiative transfer boundary value problems with nonhomogeneous reflectivity boundary conditions (one for each wavelength band in the semi-transparent electromagnetic spectrum of the glass) and a nonlinear heat conduction evolution equation with a nonlinear Robin boundary condition which takes into account those wavelengths for which the glass behaves like an opaque body. We prove existence and uniqueness of the solution, and give also uniform bounds on the solution i.e. on the absolute temperature distribution inside the body and on the radiative intensities. Now, we consider the temperature $T_{S}$ of the black radiative source S surrounding the semi-transparent body $\Omega$ as the control variable. We adjust the absolute temperature distribution (x, t) 7! T(x, t) inside the semi-transparent body near a desired temperature distribution Td(·, ·) during the time interval of radiative heating ]0, tf [ by acting on $T_{S}$. In this respect, we introduce the appropriate cost functional and the set of admissible controls $T_{S}$, for which we prove the existence of optimal controls. Introducing the State Space and the State Equation, a first order necessary condition for a control $T_{S}$ : t 7! $T_{S}$ (t) to be optimal is then derived in the form of a Variational Inequality by using the Implicit Function Theorem and the adjoint problem. We come now to the goal problem which is the deformation of the semi-transparent body $\Omega$ by heating it with a black radiative source surrounding it. We introduce a weak mixed formulation of this thermoviscoelasticity problem and study the existence and uniqueness of its solution, the novelty here with respect to the work of M.E. Rognes et R. Winther (M3AS, 2010) being the apparition of the viscosity in some of the coefficients of the constitutive equation, viscosity which depends on the absolute temperature T(x, t) and thus in particular on the time t. Finally, we state in this setting the related optimal control problem of the deformation of the semi-transparent body $\Omega$, by acting on the absolute temperature of the black radiative source surrounding it. We prove the existence of an optimal control and we compute the Fréchet derivative of the associated reduced cost functional.