A la frontière entre l'atmosphère neutre et l'espace, le plasma ionosphérique est le siège de réactions d'ionisation complexes. Le champ magnétique terrestre et les champs électriques induits causent des fluctuations spatiales et temporelles de la concentration électronique. Ces irrégularités ionosphériques entrainent des variations rapides de l'amplitude et de la phase des signaux radioélectriques les traversant, notamment aux hautes latitudes. Ce phénomène est appelé scintillation ionosphérique et il est particulièrement craint par la communauté utilisatrice d'applications GNSS qui nécessite une disponibilité et une intégrité optimales des signaux. Le travail présenté dans cette thèse propose une modélisation complète, à 3 axes d'anisotropie, de la scintillation ionosphérique. Ce modèle est basé sur une approche numérique 3D et 2D, de type écrans de phase, et sur la résolution analytique des équations de propagation, en 3D et en 2D. Ces dérivations originales des variances et des spectres de log-amplitude et de phase ont mis en relief les limites de validité d'un modèle numérique 2D. L'étude de sensibilité menée sur les variances et les spectres ouvre également des perspectives d'inversion des données GNSS pour remonter aux caractéristiques du milieu ionosphérique. / The ionospheric plasma is located at the border between neutral atmosphere and outer space and many complex ionization reactions occur inside this turbulent medium. The Earth magnetic field and induced electric fields cause rapid fluctuations of electron density, both spatially and temporarily. When crossing this turbulent layer, RF signals show fast variations of amplitude and phase, especially at high latitudes. This phenomenon is called ionospheric scintillation and it is particularly feared by air navigation using GNSS services, since it degrades the availability and the integrity of signals. This PhD dissertation presents a complete modeling of the effects of ionospheric scintillation, with 3 anisotropy axes. It is based on a numerical approach using the multiple phase screens technique, both in 3D and 2D schemes, and on the analytical resolution of electromagnetic propagation equation, also both in 3D and 2D configurations. The limits of use of a 2D numerical scheme have been outlined by these original formulations of phase and log-amplitude variances and spectra. This complete modeling associated with a sensitivity study on these variances and spectra opens up interesting perspectives on data inversion, in order to better estimate the physical characteristics of the ionospheric medium.

Cette thèse a pour objet l'analyse mathématique de lois de conservation scalaires dont la fonction flux présente une discontinuité par rapport à la variable d'espace. Nous nous intéressons plus particulèrement au problème du raccord le long d'une interface commune des solutions de deux équations quasi linéaires hyperboliques du premier ordre, posées dans deux ouverts disjoints. <br /> En premier lieu nous considérons un problème couplé hyperbolique/hyperbolique. Sous une condition de non dégénérescence du flux, nous avons obtenu un résultat d'existence et d'unicité d'une solution faible entropique d'abord en dimension 1 d'espace puis en dimension quelconque. La preuve de l'unicité est basée sur la méthode de dédoublement des variables due à S.N. Kruzkov puis sur un raisonnement presque partout à l'interface. Dans le cas particulier de la dimension 1 l'existence s'obtient par une régularisation adéquate du coefficient discontinu dans le terme de convection alors que nous utilisons la méthode de viscosité artificielle dans le cas général. <br />En second lieu nous traitons le cas de termes de convection qui apparaissent dans l'ingénierie pétrolière pour lesquels la condition de non dégénérescence de la non linéarité n'est pas vérifiée. Nous ne pouvons donc pas adapter les méthodes précédemment utilisées. Nous nous sommes donc intéressés à un problème couplé perturbé où sur l'un des deux ouverts un terme de diffusion est ajouté. Sous l'hypothèse que les caractéristiques provenant de la zone hyperbolique sont sortantes à l'interface, l'unicité d'une solution faible entropique est établie. La méthode de viscosité artificielle et la notion de processus entropique nous permettent de prouver le résultat d'existence .

Dans cette recherche, un modèle de propagation d'ondes de choc sur grandes distances sur un environnement urbain est construit et validé. L'approche consiste à utiliser l'Equation Parabolique Nonlinéaire (NPE) comme base. Ce modèle est ensuite étendu afin de prendre en compte d'autres effets relatifs à la propagation du son en milieu extérieur (surfaces non planes, couches poreuses, etc.). La NPE est résolue en utilisant la méthode des différences finies et donne des résultats en accord avec d'autres méthodes numériques. Ce modèle déterministe est ensuite utilisé comme base pour la construction d'un modèle stochastique de propagation sur environnements urbains. La Théorie de l'Information et le Principe du Maximum d'Entropie permettent la construction d'un modèle probabiliste d'incertitudes intégrant la variabilité du système dans la NPE. Des résultats de référence sont obtenus grâce à une méthode exacte et permettent ainsi de valider les développements théoriques et l'approche utilisée.

Dans cette thèse, nous nous intéressons à des méthodes numériques pour un modèle d'écoulements incompressibles et miscibles ayant des application dans l'hydrogéologie et l'ingénierie pétrolière. Nous étudions et analysons un schéma numérique combinant une méthode d'éléments finis mixtes (EFM) et une méthode des volumes finis (VF) pour approcher le système couplé entre une équation elliptique (pression-vitesse) et une équation de convection-diffusion-réaction (concentration). Le schéma VF considérée est de type "vertex centred" semi-implicite en temps : explicite pour la convection et implicite pour la diffusion. On utilise un schéma de Godunov pour approcher le terme convectif et une approximation élément fini P1 pour le terme de diffusion. Nous montrons des résultats de stabilité L≂ estimations BV et le principe du maximum discret sous une condition CFL appropriée. Ensuite, nous montrons la convergence de la solution approchée obtenue par le schéma combiné EFM-VF vers la solution du problème couplé. La démonstration de la convergence se fait en plusieurs étapes : premièrement, on déduit la convergence forte de la solution approchée de la concentration dans L2(Q), en utilisant la stabilité L≂, les estimations BV et des arguments de compacité. Dans l'étape suivante, on étudie le schéma découplé EFM, en donnant des résultats de convergence pour la pression et la vitesse. Enfin, le processus de convergence de la solution approchée du schéma combiné EFM-VF vers la solution exacte est obtenu par passage à la limite et par unicité de solution pour le problème continu. Des simulations numériques académiques et réalistes pour des problèmes bidimensionnels confirment la stabilité et l'efficacité du schéma combiné. Enfin, nous étudions des estimateurs d'erreur a posteriori de type résiduel pour une équation de convection-diffusion-réaction discrétisée par un schéma VF "vertex centred" semi-implicite en temps. Nous introduisons deux sortes d'indicateurs. Le premier est local en temps et en espace et constitue un outil efficace pour l'adaptation du maillage à chaque pas de temps. Le second est global en espace mais local en temps et peut être utilisé pour l'adaptation en temps. Nous montrons que l'estimateur est une borne supérieure de l'erreur. Des résultats numériques d'adaptations de maillage sont présentés et montrent l'efficacité de la méthode. La partie logiciels de ce travail porte sur deux volets. Le premier a permis de réaliser un code de calcul 2D, MFlow, écrit en C++, pour la résolution du système des écoulements miscibles considérés dans cette thèse. Le second volet concerne la collaboration avec un groupe de chercheurs pour l'élaboration de la plate-forme Homogenizer++ réalisée dans le cadre du GDR MoMaS (http://momas.univ-lyon1.fr/).

Étant donné l'urgence environnementale, le développement des énergies renouvelables s'est fortement accru ces dernières années. L'implantation de champs d'éoliennes est notamment en pleine expansion dans toute l'Europe. Ces éoliennes, de structure diélectrique et métallique et de grande taille, peuvent avoir un impact significatif sur les systèmes radiofréquences. En particulier, les systèmes de radionavigation et de surveillance opérés par la DGAC (VOR, radar) doivent cohabiter avec de nouveaux champs d'éoliennes. En effet, ces dernières influent sur le champ électromagnétique des systèmes et peuvent dégrader leurs performances (multitrajets, masquages... ). Il est alors important de pouvoir quantifier ces dégradations, en particulier dans le cas du VOR où les multitrajets engendrent une erreur de relèvement. Dans ce travail de thèse, un simulateur électromagnétique appelé VERSO (VOR ERror SimulatOr) est développé. Il permet d'estimer l'impact d'objets diffractants, en particulier d'éoliennes, sur le signal VOR. Dans la littérature, différentes techniques de modélisation sont proposées pour prédire ces phénomènes. Certaines sont trop approximatives, d'autres trop coûteuses en temps. Ainsi, le choix des techniques utilisées dans ce simulateur a été guidé par le compromis entre précision et temps de calcul. L'équation parabolique est utilisée pour modéliser la propagation de la source jusqu'aux éoliennes afin de prendre en compte le relief. Ensuite, une méthode d'optique physique sur matériaux diélectriques est mise en œuvre pour calculer le champ diffracté par ces objets. Le modèle électromagnétique de l'éolienne et les hypothèses inhérentes aux méthodes utilisées par VERSO ont été validés aux fréquences VHF (VOR) par comparaison avec la méthode des moments qui fait office de référence. Une extension de VERSO pour les systèmes radars a été abordée. Par conséquent des validations similaires à ces fréquences ont été réalisées. Plus spécifiquement, un modèle de pale d'éolienne prenant en compte la présence du parafoudre est proposé aux fréquences VOR et radar. Les effets de masquage sont eux aussi quantifiés pour ces 2 domaines de fréquences. Il est notamment démontré que l'effet de masquage pour des éoliennes alignées radialement à un VOR est négligeable. Cette approximation ainsi que le modèle de pale sont ensuite utilisés dans le programme VERSO. Ce dernier est validé à l'aide de contrôles en vol sur un scénario de 9 éoliennes implantées à 5 km du VOR de Boulogne-sur-Mer. Une étude quantitative de l'impact de chaque partie des éoliennes est menée afin de discriminer la source majoritaire d'erreur VOR. On constate notamment qu'à grande distance du VOR dans le scénario d'observation considéré, le mât constitue le principal contributeur en terme de champ diffracté et d'erreur VOR. Enfin, une étude statistique sur l'erreur VOR a permis d'obtenir un simulateur de scénarios qui donne l'erreur maximale avec une confiance fixée, en minimisant le nombre de simulations à effectuer. Pour accélérer cette méthode, une expression analytique approchée de l'erreur VOR maximum a été développée en fonction de la distance d'implantation et de la hauteur du mât de l'éolienne. Cette dernière étude fournit une méthode rapide pour évaluer l'impact de la construction d'un champ d'éoliennes quelconque à proximité de systèmes de l'aviation civile. / Considering the ecological emergency, the renewable energy development has greatly increased for a decade. In particular, the windfarms implantation rapidly expands in Europe. These windturbines are large obstacles composed by dielectric and metallic materials. So their impact on electromagnetic devices is significant. The radionavigation systems for the civil aviation services are particularly concerned. However, they have to work side with new windfarms. Actually, these latter cause scattering effects on the electromagnetic signals and can degrade the performances of these equipments (multipaths, shadowing effects etc.). Thus, quantifying these degradations is crucial, particularly on the VOR devices where multipath effects cause an error on the azimuth. In this thesis work, an electromagnetic simulator called VERSO (VOR ERror SimulatOr) is developed. It can estimate the impact of scattering objects, especially windturbines, on the VOR signal. In literature, several techniques are proposed to model these phenomena: some make coarse approximations and some others are memory intensive. Thus, the choice over the methods used in VERSO is a compromise between precision and memory requirement. The parabolic equation is used to model the propagation from the source to the windturbines so as to take the relief into account. A physical optic based method is used to compute the field scattered by these objects. The electromagnetic model of the windturbine and the hypothesis due to the methods used by VERSO have been validated in the VHF (VOR) frequency by comparison with the method of moments, which is the reference. An extension of VERSO for the radar systems is introduced. Consequently similar validations have been performed at radar frequencies. In particular, a windturbine blade model taking into account the lightning protection is proposed for the VOR and the radar frequencies. The shadowing effects are also quantified in both frequency domains. Especially, a demonstration that the shadowing effects due to radially implemented windturbines can be neglected around a VOR beacon is proposed. This approximation and the blade model are used for the implementation of VERSO. This simulator is validated by comparison with measurements on 9 windturbines built 5~km far from a VORC in Boulogne-sur-Mer (France). A study is performed to quantify the influence of each part of the windturbine. The mast is shown to be the main contributor regarding to the electromagnetic field and the VOR error. Finally, parametric simulations are performed and analytic expressions are proposed to describe the evolution of the maximum VOR error with respect to the mast size and the distance VOR-windturbine. The latter study gives some key parameters that need to be considered for the elaboration of a windfarm building plan close to civil aviation systems for the project to be viable.

Cette thèse est dédiée a l''etude de certaines propriétés des équations d' évolutions non-autonomes $u'(t)+A(t)u(t)=f(t), u(0)=x.$ Il s'agit précisément de la propriété de la régularité maximale $L^p$: étant donnée $fin L^{p}(0,tau;H)$, montrer l'existence et unicité de la solution $u in W^{1,p}(0,tau;H)$. Ce problème a 'et'e intensivement étudie dans le cas autonome, i.e., $A(t)=A$ pour tout $t$. Dans le cas non-autonome, le problème a été considéré par J.L.Lions en 1960. Nous montrons divers résultats qui étendent tout ce qui est connu sur ce problème. On suppose ici que la famille des opérateurs $(mathcal{A}(t))_{tin [0,tau]}$ est associée à des formes quasi-coercives, non autonomes $(fra(t))_{t in [0,tau]}.$ Nous considérons également le problème de régularité maximale pour les d'ordre 2 (équations des ondes). Plusieurs exemples et applications sont considérés. / This Thesis is devoted to certain properties of non-autonomous evolution equations $u'(t)+A(t)u(t)=f(t), u(0)=x.$ More precisely, we are interested in the maximal $L^p$-regularity: given $fin L^{p}(0,tau;H),$ prove existence and uniqueness of the solution $u in W^{1,p}(0,tau;H)$. This problem was intensively studied in the autonomous cas, i.e., $A(t)=A$ for all $t.$ In the non-autonomous cas, the problem was considered by J.L.Lions in 1960. We prove serval results which extend all previously known ones on this problem. Here we assume that the familly of the operators $(mathcal{A}(t))_{tin [0,tau]}$ is associated with quasi-coercive, non-autonomous forms $(fra(t))_{t in [0,tau]}.$ We also consider the problem of maximal regularity for second order equations (the wave equation). Serval examples and applications are given in this Thesis.

L'environnement urbain est le siège de fortes nuisances sonores notamment générées par les moyens de transport. Afin de lutter contre ces nuisances, la réglementation européenne impose la réalisation de cartographies de bruit. Dans ce contexte, des travaux fondamentaux sont menés autour de la propagation d'ondes acoustiques basses fréquences en milieu urbain. Différents travaux de recherche récents portent sur la mise en œuvre de méthodes ondulatoires pour la propagation d'ondes acoustiques dans de tels milieux. Le coût numérique de ces méthodes limite cependant leur utilisation dans un contexte d'ingénierie. L'objectif de ces travaux de thèse porte sur l'approximation unidirectionnelle de la propagation des ondes, appliquée à l'acoustique urbaine. Cette approximation permet d'apporter des simplifications à l'équation d'onde afin de limiter le temps de calcul lors de sa résolution. La particularité de ce travail de thèse réside dans la prise en compte des variations, continues ou discontinues, de la largeur des rues. Deux formalismes sont utilisés : l'équation parabolique et une approche multimodale. L'approche multimodale sert de support à une étude théorique sur les mécanismes de couplages de modes dans des guides d'ondes irréguliers bidimensionnels. Pour cela, le champ de pression est décomposé en fonction du sens de propagation des ondes à la manière d'une série de Bremmer. La contribution particulière de l'approximation unidirectionnelle est étudiée en fonction des paramètres géométriques du guide d'ondes, ce qui permet de mieux cerner les limites de validité de cette approximation. L'utilisation de l'équation parabolique a pour but une application à l'acoustique urbaine. Une transformation de coordonnées est associée à l'équation parabolique grand angle afin de prendre en compte l'effet de la variation de la section du guide d'ondes. Une méthode de résolution est alors spécifiquement développée et permet une évaluation précise du champ de pression. D'autre part, une méthode de résolution de l'équation parabolique grand angle tridimensionnelle est adaptée à la modélisation de la propagation acoustique en milieu urbain. Cette méthode permet de tenir compte des variations brusques ou continues de la largeur de la rue. Une comparaison avec des mesures sur maquette de rue à échelle réduite permet de mettre en avant les possibilités de la méthode. / The urban environment is the seat of loud noise generated by means of transportation. To fight against these nuisances, European legislation requires the achievement of noise maps. In this context, fundamental work is carried around the propagation of acoustic low-frequency waves in urban areas. Several recent research focuses on the implementation of wave methods for acoustic wave propagation in such environments. The computational cost of these methods, however, limits their use in the context of engineering. The objective of this thesis focuses on the one-way approximation of wave propagation, applied to urban acoustics. This approximation allows to make simplifications on the wave equation in order to limit the computation time. The particularity of this thesis lies in the consideration of variations, continuous or discontinuous, of the width of streets. Two formalisms are used: parabolic equation and a multimodal approach. The multimodal approach provides support for a theoretical study on the mode-coupling mechanisms in two-dimensional irregular waveguides. For this, the pressure field is decomposed according to the direction of wave propagation in the manner of a Bremmer series. The specific contribution of the one-way approximation is studied as a function of the geometric parameters of the waveguide, which helps identify the limits of validity of this approximation. Use of the parabolic equation is intended for application to urban acoustic. A coordinate transformation is associated with the wide-angle parabolic-equation in order to take into account the variation effect of the waveguide section. A resolution method is developed specifically and allows an accurate assessment of the pressure field. On the other hand, a solving method of the three-dimensional parabolic-equation is suitable for the modeling of acoustic propagation in urban areas. This method takes into account sudden or continuous variations of the street width. A comparison with measurements on scaled model of street allows to highlight the possibilities of the method.

Dans cette thèse, nous étudions l’équation parabolique non linéaire ∂ t u = ∆u + a |x|⎺⥾ |u|ᵅ u, t > 0, x ∈ Rᴺ \ {0}, N ≥ 1, ⍺ ∈ R, α > 0, 0 < Ƴ < min(2,N) et avec une donnée initiale u(0) = φ. On établit l’existence et l’unicité locale dans Lq(Rᴺ) et dans Cₒ(Rᴺ). En particulier, la valeur q = N ⍺/(2 − γ) joue un rôle critique. Pour ⍺ > (2 − γ)/N, on montre l’existence de solutions auto-similaires globales avec données initiales φ(x) = ω(x) |x|−(2−γ)/⍺, où ω ∈ L∞(Rᴺ) homogène de degré 0 et ||ω||∞ est suffisamment petite. Nous montrons ainsi que si φ(x)∼ω(x) |x| ⎺(²⎺⥾)/⍺ pour |x| grande, alors la solution est globale et asymptotique dans L∞(Rᴺ) à une solution auto-similaire de l’équation non linéaire. Tandis que si φ(x)∼ω(x) |x| (x)|x|−σ pour des |x| grandes avec (2 − γ)/⍺ < σ < N, alors la solution est globale, mais elle est asymptotique dans L∞(Rᴺ) à eᵗ∆(ω(x) |x|−σ). L’équation avec un potentiel plus général, ∂ t u = ∆u + V(x) |u|ᵅ u, V(x) |x |⥾ ∈ L∞(Rᴺ), est également étudiée. En particulier, pour des données initiales φ(x)∼ω(x) |x| ⎺(²⎺⥾)/⍺, |x| grande, nous montrons que le comportement à grand temps est linéaire si V est à support compact au voisinage de l’origine, alors qu’il est non linéaire si V est à support compact au voisinage de l’infini. Nous étudions également le système non linéaire ∂ t u = ∆u + a |x|⎺⥾ |v|ᴾ⎺¹v, ∂ t v = ∆v + b |x|⎺ ᴾ |u|q⎺¹ u, t > 0, x ∈ Rᴺ \ {0}, N ≥ 1, a,b ∈ R, 0 < y < min(2,N)? 0 < p < min(2,N), p,q > 1. Sous des conditions sur les paramètres p, q, γ et ρ nous montrons l’existence et l’unicité de solutions globales avec données initiales petites par rapport à certaines normes. En particulier, on montre l’existence de solutions auto-similaires avec donnée initiale Φ = (φ₁, φ₂), où φ₁, φ₂ sont des données initiales homogènes. Nous montrons également que certaines solutions globales sont asymptotiquement auto-similaires. Comme deuxième objectif, nous considérons l’équation de la chaleur non linéaire ut = ∆u + |u|ᴾ⎺¹u - |u| q⎺¹u, avec t ≥ 0 et x ∈ Ω, la boule unité de Rᴺ, N ≥ 3, avec des conditions aux limites de Dirichlet. Soit h une solution stationnaire à symétrie radiale avec changement de signe de (E). On montre que la solution de (E) avec donnée initiale λh explose en temps fini si |λ − 1| > 0 est suffisamment petit et si 1 < q < p < Ps = N+2/N−2 et p suffisamment proche de Ps. Ceci prouve que l’ensemble des données initiales pour lesquelles la solution est globale n’est pas étoilé au voisinage de 0. / In this thesis, we study the nonlinear parabolic equation ∂ t u = ∆u + a |x|⎺⥾ |u|ᵅ u, t > 0, x ∈ Rᴺ \ {0}, N ≥ 1, ⍺ ∈ R, α > 0, 0 < Ƴ < min(2,N) and with initial value u(0) = φ. We establish local well-posedness in Lq(Rᴺ) and in Cₒ(Rᴺ). In particular, the value q = N ⍺/(2 − γ) plays a critical role.For ⍺ > (2 − γ)/N, we show the existence of global self-similar solutions with initial values φ(x) = ω(x) |x|−(2−γ)/⍺, where ω ∈ L∞(Rᴺ) is homogeneous of degree 0 and ||ω||∞ is sufficiently small. We then prove that if φ(x)∼ω(x) |x| ⎺(²⎺⥾)/⍺ for |x| large, then the solution is global and is asymptotic in the L∞-norm to a self-similar solution of the nonlinear equation. While if φ(x)∼ω(x) |x| (x)|x|−σ for |x| large with (2 − γ)/α < σ < N, then the solution is global but is asymptotic in the L∞-norm toe t(ω(x) |x|−σ). The equation with more general potential, ∂ t u = ∆u + V(x) |u|ᵅ u, V(x) |x |⥾ ∈ L∞(Rᴺ), is also studied. In particular, for initial data φ(x)∼ω(x) |x| ⎺(²⎺⥾)/⍺, |x| large , we show that the large time behavior is linear if V is compactly supported near the origin, while it is nonlinear if V is compactly supported near infinity. we study also the nonlinear parabolic system ∂ t u = ∆u + a |x|⎺⥾ |v|ᴾ⎺¹v, ∂ t v = ∆v + b |x|⎺ ᴾ |u|q⎺¹ u, t > 0, x ∈ Rᴺ \ {0}, N ≥ 1, a,b ∈ R, 0 < y < min(2,N)? 0 < p < min(2,N), p,q > 1. Under conditions on the parameters p, q, γ and ρ we show the existence and uniqueness of global solutions for initial values small with respect of some norms. In particular, we show the existence of self-similar solutions with initial value Φ = (φ₁, φ₂), where φ₁, φ₂ are homogeneous initial data. We also prove that some global solutions are asymptotic for large time to self-similar solutions. As a second objective we consider the nonlinear heat equation ut = ∆u + |u|ᴾ⎺¹u - |u| q⎺¹u, where t ≥ 0 and x ∈ Ω, the unit ball of Rᴺ, N ≥ 3, with Dirichlet boundary conditions. Let h be a radially symmetric, sign-changing stationary solution of (E). We prove that the solution of (E) with initial value λ h blows up in finite time if |λ − 1| > 0 is sufficiently small and if 1 < q < p < Ps = N+2/N−2 and p sufficiently close to Ps. This proves that the set of initial data for which the solution is global is not star-shaped around 0.

Ce mémoire de thèse est consacrée à l'étude des problèmes d'évolution où la dynamique est régi par l'opérateur de diffusion de sous-gradient. Nous nous intéressons à deux types de problèmes d'évolution. Le premier problème est régi par un opérateur local de type Leray-Lions avec un domaine borné. Dans ce problème, l'opérateur est maximal monotone et ne satisfait pas la condition standard de contrôle de la croissance polynomiale. Des exemples typiques apparaît dans l'étude de fluide non-Neutonian et aussi dans la description de la dynamique du flux de sous-gradient. Pour étudier le problème nous traitons l'équation dans le contexte de l'EDP non linéaire avec le flux singulier. Nous utilisons la théorie de gradient tangentiel pour caractériser l'équation d'état qui donne la relation entre le flux et le gradient de la solution. Dans le problème stationnaire, nous avons l'existence de la solution, nous avons également l'équivalence entre le problème minimisation initial, le problème dual et l'EDP. Dans l'équation de l'évolution, nous proposons l'existence, l'unicité de la solution. Le deuxième problème est régi par un opérateur discret. Nous étudions l'équation d'évolution discrète qui décrivent le processus d'effondrement du tas de sable. Ceci est un exemple typique de phénomènes auto-organisés critiques exposées par une slope critique. Nous considérons l'équation d'évolution discrète où la dynamique est régie par sous-gradient de la fonction d'indicateur de la boule unité. Nous commençons par établir le modèle, nous prouvons existence et l'unicité de la solution. Ensuite, en utilisant arguments de dualité nous étudions le calcul numérique de la solution et nous présentons quelques simulations numériques. / This thesis is devoted to the study of evolution problems where the dynamic is governed by sub-gradient diffusion operator. We are interest in two kind of evolution problems. The first problem is governed by local operator of Leray-Lions type with a bounded domain. In this problem, the operator is maximal monotone and does not satisfied the standard polynomial growth control condition. Typical examples appears in the study of non-Neutonian fluid and also in the description of sub-gradient flows dynamics. To study the problem we handle the equation in the context of nonlinear PDE with singular flux. We use the theory of tangential gradient to characterize the state equation that gives the connection between the flux and the gradient of the solution. In the stationary problem, we have the existence of solution, we also get the equivalence between the initial minimization problem, the dual problem and the PDE. In the evolution one, we provide the existence, uniqueness of solution and the contractions. The second problem is governed by a discrete operator. We study the discrete evolution equation which describe the process of collapsing sandpile. This is a typical example of Self-organized critical phenomena exhibited by a critical slop. We consider the discrete evolution equation where the dynamic is governed by sub-gradient of indicator function of the unit ball. We begin by establish the model, we prove existence and uniqueness of the solution. Then by using dual arguments we study the numerical computation of the solution and we present some numerical simulations.

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.