Le travail réalisé durant cette thèse a eu pour objectif de développer un modèle théorique spécifiquement dédié au contexte du développement tumoral tel qu'il peut être observé sous une fenêtre dorsale implantée sur une souris. Le modèle développé est un modèle hybride multi-physique et multi-échelle qui couple deux modules principaux. Le premier module modélise la croissance tumorale par un automate cellulaire qui permet de différencier l'état de chaque cellule en fonction de son histoire (cycle cellulaire), et de son environnement (espace disponible pour proliférer, présence d'oxygène). Le second module modélise le réseau vasculaire et le flux sanguin et rend compte de l'angiogenèse (apparition de nouveaux vaisseaux) et de l'adaptation du diamètre des vaisseaux, en fonction de l'évolution des contraintes hémodynamiques, nettement visible sous la fenêtre dorsale. L'ensemble des processus diffusifs (diffusion de l'oxygène et des facteurs de croissance vasculaire) sont décrits par des équations aux dérivées partielles, couplées à des automates cellulaires qui permettent de localiser à chaque instant pour chaque équation les termes sources (production) et les termes puits (consommation) pour chaque entité diffusive. Les simulations numériques réalisées montrent dans quelle mesure il est possible de rendre compte des observations expérimentales sur le plan qualitatif, qui nécessite la neutralisation des biais numériques ; et sur le plan quantitatif, pour reproduire la cinétique de croissance tumorale et l'évolution de la densité vasculaire. Le modèle numérique de l'évolution tumorale sous fenêtre dorsale est ensuite utilisé pour tester les effets de deux types de molécules : cytotoxiques et anti-vasculaires. Les simulations numériques de ces deux types de traitement explorent différents protocoles, définis par le mode d'action de la molécule, la dose administrée et la fréquence d'administration. Les résultats montrent comment il est alors possible de définir un protocole optimum pour une tumeur donnée en direction d'une individualisation de la thérapie. Ce modèle intégré a permis de poser de façon satisfaisante les bases d'un clone numérique du modèle expérimental d'évolution tumorale sous fenêtre dorsale même si certains aspects nécessitent encore quelques améliorations. La validation des aspects thérapeutiques restera encore à accomplir avant de pouvoir envisager à terme le remplacement (au moins partiel) de l'animal par l'ordinateur.

La "cryopolymérisation" permet d'obtenir des gels de polymère macroporeux ou "cryogels". Cette méthode a été utilisée pour la synthèse d'hydrogels thermosensibles à base de pNIPA. La température critique TC correspondant à la transition de volume a été déterminée par des mesures de taux de gonflement et par DSC. La macroporosité (distribution de la taille des pores et épaisseur des parois) et son évolution en fonction de T ont été étudiées par la microscopie biphotonique donnant des informations à l'échelle du µm à plusieurs dizaines de µm. La diffusion de rayons X (SAXS et WAXS) a été utilisée pour caractériser la structure multi-échelle (de quelques dixièmes à quelques dizaines de nm) du gel constituant les parois des macropores. Les courbes de diffusion ont été décrites analytiquement. L'évolution des dix paramètres contenus dans l'équation a été étudiée en fonction de T et discutée. Enfin, des expériences utilisant les phonons hyperfréquences générés par la technique des réseaux transitoires avec détection hétérodyne (HD-TG) ont été réalisées. Ces mesures ont permis de déterminer la vitesse de propagation de l'onde ultra-sonore (à 340 MHz), son atténuation, et la constante de diffusion thermique à différentes températures.

L'étude du comportement du béton en température est un problème majeur qui vise in fine à évaluer le niveau de sécurité des structures sous des sollicitations thermiques sévères, lors d'incendies par exemple. Pour cela, de nombreux modèles sont développés dans un cadre de couplage thermo-hydro-mécanique (THM), pour prendre en compte les différents processus physico-chimiques et mécaniques mobilisés par ces sollicitations et conditionnant la stabilité du matériau béton.Cependant, ces modélisations ne prennent souvent pas en compte explicitement la nature hétérogène du matériau. En effet, le béton est un matériau constitué d'inclusions noyées dans une matrice cimentaire pouvant être schématisée comme un milieu poreux ouvert partiellement saturé en eau. Les inclusions sont caractérisées par leurs natures minéralogiques, leurs morphologies et leurs tailles classées en fuseaux granulaires. Cette hétérogénéité introduit une complexité supplémentaire : la nécessité de prendre en compte la microstructure pour quantifier l'effet de l'incompatibilité (thermique, hydrique et mécanique) inclusion-matrice sur le comportement THM du béton. Ce travail constitue une première avancée dans ce sens. A ce titre, une modélisation élément fini multi-échelle tridimensionnelle (3D) est développée permettant d'affecter des comportements spécifiques à la matrice et aux inclusions. Pour la matrice, siège des transports de masse dans son réseau poreux connecté, une approche thermo-hydrique à trois fluides (eau, vapeur et air sec) est adoptée et est couplée à une loi de comportement poro-mécanique endommageable. Les inclusions (granulats) sont considérées hydriquement inertes, une approche thermo-mécanique avec endommagement est alors retenue.Une analyse, par simulations numériques, des effets de la nature minéralogique des granulats (calcaires ou silico-calcaires) de leurs distributions et de leurs morphologies a été menée sur des configurations 2D et 3D. Les effets étudiés ont notamment concerné l'influence de ces paramètres sur les fluctuations locales des champs de température, de pression de gaz et d'endommagement au regard de la dispersion des mesures expérimentales. L'analyse est limitée à l'échelle mésoscopique, celle où les granulats de taille caractéristique supérieure à 5 mm sont considérés, le reste (stable et pâte de ciment) étant une matrice homogène. Enfin, cette analyse a mis en évidence le besoin de mettre en œuvre une approche expérimentale cohérente avec une analyse mutli-échelle, à la fois pour la caractérisation des propriétés (thermiques, hydriques et mécaniques) de chaque constituant et pour l'étude des évolutions des champs lors des changements d'échelles. Un protocole expérimental a été définit et des premiers résultats d'essais sont présentés et analysés au travers de résultats obtenus dans la partie modélisation

Cette thèse se concentre sur le développement d'une méthode précise et efficace pour la simulation des grandes échelles (LES) des écoulements turbulents. Une approche de la LES basée sur la méthode variationnelle multi-échelles (VMS) est considérée. La VMS applique aux équations de la dynamique des fluides une séparation d'échelles a priori sans recours à des hypothèses sur les conditions aux limites ou sur l'uniformité du maillage. Afin d'assurer effectivement une séparation d'échelles dans l'espace des nombres d'onde associé, nous choisissons d'utiliser les ondelettes de deuxième génération (SGW), une base polynomiale qui présente des propriétés de localisation spatiale-fréquence optimales. A partir de la séparation d'échelles ainsi réalisée, l'action du modèle sous-maille est limitée à un intervalle de nombres d'onde proche de la coupure spectrale. Cette approche VMS-LES basée sur les ondelettes est désignée par WAVVMS-LES. Elle est incorporée dans un solveur d'ordre élevé pour la simulation des écoulements incompressibles sur la base d'une méthode de Galerkin discontinue (DG-FEM) stabilisée pour la pression. La méthode est évaluée par réalisation de LES sur des maillages fortement sous-résolus pour le cas test du tourbillon de Taylor-Green 3D à deux nombres de Reynolds différents. / This thesis focuses on the development of an accurate and efficient method for performing Large-Eddy Simulation (LES) of turbulent flows. An LES approach based upon the Variational Multiscale (VMS) method is considered. VMS produces an a priori scale-separation of the governing equations, in a manner which makes no assumptions on the boundary conditions and mesh uniformity. In order to ensure that scale-separation in wavenumber is achieved, we have chosen to make use of the Second Generation Wavelets (SGW), a polynomial basis which exhibits optimal space-frequency localisation properties. Once scale-separation has been achieved, the action of the subgrid model is restricted to the wavenumber band closest to the cutoff. We call this approach wavelet-based VMS-LES (WAV-VMS-LES). This approach has been incorporated within the framework of a high-order incompressible flow solver based upon pressure-stabilised discontinuous Galerkin FEM (DG-FEM). The method has been assessed by performing highly under-resolved LES upon the 3D Taylor-Green Vortex test case at two different Reynolds numbers.

Ce travail porte sur l’étude de la dégradation des convertisseurs électrochimiques haute température à oxydes solides. Une approche couplant des tests électrochimiques, des caractérisations post-mortem avancées et une modélisation multi-échelle a été mise en place afin d’établir les liens entre les performances, la microstructure des électrodes et leur dégradation. Dans ce but, des essais de durabilité de plus de mille heures ont été menés dans différentes conditions opératoires. La microstructure des électrodes a été reconstruite par nano-holotomographie des rayons X pour la cellule de référence avant et après vieillissement. Une attention particulière a été apportée à la mesure de la résolution spatiale et à la fiabilisation du protocole expérimental. Grâce aux volumes 3D, les propriétés microstructurales de l’électrode H2 en Ni-YSZ ont été quantifiées pour les cellules à l’état initial et vieillies. Un modèle physique d’agglomération des particules de Nickel a ensuite été ajusté sur les analyses tridimensionnelles et intégré dans une structure de modélisation multi-échelle développée au laboratoire. Il a auparavant été nécessaire de compléter l’outil numérique avec un module spécifique dédié aux matériaux composant l’électrode à oxygène fait avec un conducteur mixte ionique-électronique. Une fois le modèle validé sur des courbes de polarisation expérimentales, il a été utilisé pour quantifier la contribution de l’agglomération du Nickel sur les pertes de performances mesurées expérimentalement en mode pile à combustible et électrolyse. / This work aims at a better understanding of the high temperature Solid Oxide Cells degradation. An approach based on electrochemical tests, advanced post-test characterizations and multi-scale models has been used to investigate the links between the performances, the electrodes microstructure and their degradation. In that goal, long-term durability tests have been performed over thousand hours in different operating conditions. Electrode microstructures have been reconstructed by X-ray nano-holotomography for the pristine and the aged cells. It is worth noting that a special attention has been paid to improve both the process reliability for the tomographic experiments as well as the spatial resolution of the 3D reconstructed images. Thanks to the valuable 3D volumes, the Ni-YSZ microstructural properties of the H2 electrode have been quantified for the fresh and the aged samples. Then, a physically-based model for Nickel particle agglomeration has been adjusted on the microstructural parameters obtained by the 3D analysis and implemented in an in-house multi-scale modelling framework. Beforehand, it has been necessary to enrich the available numerical tool with a specific module dedicated to the oxygen electrode made in Mixed Ionic Electronic Conducting materials. Once validated on polarisation curves, the completed model has been used to quantify the contribution of Nickel agglomeration on the experimental degradation rates recorded in fuel cell and electrolysis modes.

Une carte multi-échelle est un ensemble de cartes à différentes échelles, dans lequel l’utilisateur peut naviguer via un géoportail. Chacune de ces cartes est préalablement construite par généralisation cartographique, processus qui adapte la représentation cartographique à une échelle donnée. Les changements de représentations qu’implique la généralisation entre deux cartes à différentes échelles sont susceptibles de perturber l’utilisateur, rendant sa navigation plus difficile. Nous proposons dans cette thèse d’ajouter des représentations intermédiaires dans une carte multi-échelle existante, pour créer une évolution plus fluide du contenu cartographique au fil des échelles. Alors que de solides connaissances théoriques existent pour la conception cartographique traditionnelle, on ne sait pas encore comment concevoir une carte multi-échelle efficace. Pour formaliser des connaissances à ce sujet, nous avons étudié un panel de seize cartes multi-échelles existantes. Nous avons analysé les systèmes de zoom utilisés ainsi que l’évolution des représentations cartographiques au fil des échelles, en particulier les changements de niveaux d’abstraction pour les objets bâtis et routiers. Nous avons aussi évalué la variation de complexité visuelle du contenu cartographique au fil des échelles, en utilisant des mesures de clutter visuel. Nous avons ainsi identifié les tendances générales en termes de représentations multi-échelles (comme l’application du standard WMTS), certains facteurs que nous considérons comme ayant une influence négative sur la navigation de l’utilisateur (comme l’utilisation d’une même carte à différentes échelles), ainsi que des pratiques intéressantes visant à la faciliter (comme les représentations mixtes). A partir de ces constats, nous avons formulé des hypothèses sur l’influence des variables de construction des représentations intermédiaires sur la fluidité de navigation. Nous avons construit un matériel de test à partir d’un extrait de la carte multi-échelle Scan Express de l’IGN, entre les cartes existant au 1 : 25k et au 1 : 100k. Nous avons ainsi produit quatre versions différentes de représentations intermédiaires entre ces deux cartes, implémentant nos différentes hypothèses. Cet exercice nous a permis de mieux cerner les verrous techniques que soulève la production de représentations intermédiaires. Nous avons enfin conduit un test utilisateurs contrôlé, en demandant à 15 participants de réaliser une tâche cartographique sur ces différentes cartes multi-échelles, pour évaluer la pertinence de nos hypothèses / A multi-scale map is a set of maps at different scales, displayed on mapping applications, in which users may navigate by zooming in or out. Each of these maps is produced beforehand by cartographic generalization, which aims to adapt the cartographic representation for a target scale. Due to generalization, the representation changes between maps at different scales may disturb the user during its navigation. We assume that adding intermediate representations in an existing multi-scale map may enable a smooth evolution of cartographic content across scales. While theoretical knowledge exists for traditional cartography, we still do not know how to design efficient multi-scale maps. To formalize knowledge on that subject, we studied sixteen existing multi-scale maps. We focused on the used zooming system (zoom levels and display scales) and on the evolution of cartographic representations across scales, in particular for building and road entities. We also analyzed the variation of visual complexity of the map content across scales, using visual clutter measures. We thus identified general trends in terms of multi-scale representation (i.e. use of WMTS standard), some potential disturbing factors (i.e. use of a same map at different scales), but also good practices which may ease the user navigation (i.e. mixed representations). Based on these findings, we made assumptions on the influence of intermediate representations design on user navigation. We built test material from an extract of the Scan Express multi-scale map of the French IGN, between the existing maps at 1:25k and 1:100k scales. We thus produced four different versions of intermediate representations between these two maps, implementing our different hypotheses. This way, we highlighted the technical issues that we faced when producing intermediate representations. Finally, we conducted a controlled user study, asking 15 participants to perform a cartographic task on these different multi-scale maps, to evaluate our hypotheses

L'objectif de cette thèse était d'apporter une approche géophysique non linéaire à l'hydrologie urbaine. Elle a visé l'étude de la mise à l'échelle et de l'intermittence de la non-linéarité des précipitations, la réalisation d'une méthode de prévision stochastique à très court terme ("nowcast"), ainsi que son application aux processus hydrologiques dans les environnements (semi-) urbains. La partie modélisation hydrologique globale concerne la vallée de la Bièvre, zone semi-urbanisée de 110 km2 dans le sud-ouest de la région parisienne. Par conséquent, trois études différentes ont été réalisées dans cette zone à l'aide de deux modèles hydrologiques : le modèle conceptuel semi-distribué InfoWorks CS appliqué sur tout le bassin versant de Bièvre ; et le modèle physique complètement distribué Multi-Hydro, développé à l'École des Ponts ParisTech, appliqué sur deux sous-bassins versants de la Bièvre. Les principaux objectifs étaient de mieux comprendre les impacts de la variabilité spatio-temporelle des données pluviométriques en utilisant deux produits (les données radar bande-C de Météo-France avec une résolution de 1 km x 1 km x 5 min, et les données radar DPSRI band-X de l'ENPC à une résolution de 250 m x 250 m x 3.41 min) comme entrées pour les modèles, et d'identifier les capacités de chaque modèle pour traiter des données à une meilleur résolution, telles que la bande-X. Ensuite, les résultats obtenus démontrent que la fiabilité des simulations hydrologiques dépend intrinsèquement des caractéristiques des données pluviométriques. De plus, les données du radar bande-X pourraient mesurer des pics de précipitations plus élevés et le modèle complètement distribué était plus sensible à une meilleure résolution des données pluviométriques. Par la suite, des données de pluie provenant des radars météorologiques situés à des endroits complètement différents (Brésil, France, Japon) ont été analysées et comparées statistiquement afin d'améliorer la compréhension générale du comportement scalant des précipitations. De plus, le théorème d'intersection a été appliqué pour mettre en évidence les impacts de la variabilité spatiale d'un réseau virtuel de pluviomètres, qui a été généré en considérant l'emplacement des centres de masse de chaque sous-bassin versant de la vallée de la Bièvre. Ainsi, il a été possible d'identifier que la fractalité du réseau virtuel a conduit à une perte d'information importante des champs de pluie, biaisant leurs statistiques. Cela indique que le processus commun (largement retrouvé dans la littérature) de calibration des données radar à l'aide de pluviomètres devrait correctement prendre en compte cette fractalité. Enfin, une nouvelle approche de prévision stochastique immédiate a été proposée, à l'aide du modèle des multifractals universels (UM) en cascades continues. Cette méthode a été appliquée aux données des radars pluviométriques de la région amazonienne brésilienne et de Paris. Bien qu'il soit encore en développement et nécessite quelques améliorations, les premiers résultats obtenus avec ce modèle de prévision présenté ici sont vraiment encourageants et corroborent une fois de plus le besoin de données à haute résolution spatio-temporelle pour faire face aux crues soudaines / The focus of this thesis was to bring a nonlinear geophysical approach to urban hydrology. It aimed the study of rainfall non-linearity scaling and intermittency, achieving a stochastic very short-range forecast (nowcast) method, as well as its application to hydrological processes in (semi-) urban environments. The overall hydrological modelling part concerned the Bièvre Valley, which is a 110 km2 semi-urbanized area in the southwest of Paris region. Therefore, three different studies were performed within this area using two hydrological models: the conceptually-based semi-distributed model InfoWorks CS over the total Bièvre catchment, and the physically-based fully-distributed model developed at École des Ponts ParisTech called Multi-Hydro over two sub-catchments. The main goals were to better understand the impacts of spatio-temporal variability of rainfall data by using two products (the Météo-France C-band radar data with a resolution of 1 km x 1 km x 5 min; and the ENPC DPSRI X-band radar data at a 250 m x 250 m x 3.41 min resolution) as input to the models, and to identify the capacities of each model to deal with better resolution data, such as the X-band one. Then, the obtained results demonstrate that the reliability of the hydrological simulations are intrinsically dependent on rainfall data features. Moreover, the X-band radar data could measure higher peaks of rainfall rates and the fully-distributed model was more sensitive to better resolution rainfall data. Afterwards, different weather rainfall radar data from completely different sites (Brazil, France, Japan) were statistically analysed and compared in order to improve the general comprehension of rainfall scaling behaviour. In addition, the Intersection Theorem was applied to highlight the impacts of spatial variability of a virtual rain gauge network. The latter was generated by considering the location of each Bièvre Valley sub-catchment mass centre. Thus, it was possible to identify that the fractality of the virtual network led to an important information loss of the rainfall fields, biasing their statistics. This indicates that the common process (largely found in literature) of radar data calibration using rain gauges should be properly take into account this fractality. Finally, a new stochastic nowcast approach was proposed, using the continuous in scale cascade Universal Multifractals (UM) model. This method was applied to weather rainfall radar data from the Brazilian Amazon region and Paris. Although it is still under development and needs some improvements, the first results obtained with this forecast model presented here in this thesis are really encouraging and once more corroborate to the need of high spatio-temporal resolution data to cope flash floods

L'observation des systèmes complexes artificiels à grande échelle est difficile en raison de leur dimension et de leur ouverture. Cette thèse propose le modèle SMA-R (Système Multi-Agent Récursif) basé sur la récursivité pour l'observation multi-niveau des systèmes complexes artificiels. À partir d'un Système Multi-Agent donné, ce modèle est capable de construire une représentation abstraite multi-niveau. Le modèle d'agent récursif proposé possède un module contenant les connaissances, les capacités et le contexte de récursivité; un module pour observer les changements; deux mécanismes pour construire ou détruire les niveaux abstraits; un module d'interaction récursive prenant en charge des comportements collectifs et individuels. Pour appliquer ce modèle à des SMA physiquement décentralisés, nous avons proposé une architecture générique décentralisée d'agents récursifs en adoptant les concepts du modèle OSI qui offre en avant les capacités que nous cherchons: multi niveau, encapsulation, échange de messages virtuels et physiques... Un framework décentralisé générique a été développé permettant aux agents réels de construire les niveaux d'observation. Ce framework est appliqué à l'observation d'un réseau de capteurs sans fil. / Observation of large scale artificials systems is difficult because of their dimension and their openness. This thesis proposes a model SMA-R (Recursive Multi-Agent Systems) based on recursion for multi-level observation of artificial complex systems. From a given SMA, this model is able to build multi-level of abstractions. The model's agent has a module containning knowledge, skills and context of recursion, an eye to observe changes; two mechanisms to build or destroy an abstract level, a module of recursive interaction to support collective and individual behaviors. For applying this model to SMA physically decentralized, we proposed a generic decentralized architecture for recursive agent by adopting the concepts of the OSI model which offerts forward capabilities that we look for: multi-level encapsulation, exchange of virtual and physical messages ... A generic decentralized framework was developed allowing applicatives agents to build multi-levels observation. This framework is applied to the observation of a wireless sensor network.

La modélisation de la rupture des géo-matériaux constitue un important défi pour les applications telles que la séquestration du CO2 , le stockage de déchets nucléaires, la production des hydrocarbures ainsi que les projets de génie civil concernant les tunnels ou les excavations. L'objectif de cette thèse est de développer des lois d'évolution macroscopiques d'endommagement à partir des descriptions explicites de la rupture à l'échelle microscopique en vue de la modélisation du comportement d'endommagement à long terme des sites de stockage géologique. L'approche adoptée est basée sur l'homogénéisation par développements asymptotiques et la description énergétique de la propagation des micro-fissures, qui permettent l'obtention des lois d'endommagement et conduisent à une quantification explicite de l'énergie de l'émission acoustique associée à la rupture. Les modèles obtenus sont capables de prédire la dégradation des modules d'élasticité en raison de l'évolution des micro-fissures. Cette représentation permet de modéliser la propagation des ondes dans un milieu à endommagement évolutif. Deux types de modèles d'endommagement seront proposés: indépendants de temps et dépendants de temps. Les modèles dépendants de temps décrivent l'évolution progressive quasi-fragile de la micro-fissuration. Dans les modèles dépendants de temps, l'évolution des micro-fissures est décrite à travers un critère sous-critique et la propagation mixte, par branchement. En utilisant le modèle dépendant de temps, des simulations seront faites à trois niveaux: du laboratoire, du tunnel et du réservoir. / Accurate modeling of failure of geomaterials is the key to the success of a diverse range of engineering challenges including the topic of CO2 sequestration, nuclear waste disposal and hydrocarbon production plus civil engineering projects for tunnels or excavations. The aim of this thesis is to develop macroscopic damage evolution laws based on explicit descriptions of fracture at the micro-scale level which can be successfully employed to describe long term damage behavior of geologic storage sites. The approach taken is based on homogenization through asymptotic developments combined with micro-crack propagation energy analysis which leads to an explicit quantification of the acoustic emission (AE) energy associated with damage. Proposed damage models are capable of modeling the degradation of elastic moduli due to the micro-crack evolution. This representation allows the modeling of wave propagation in a medium with evolving damage. Two types of damage models will be considered: time independent and time dependent. Time independent damage models capable of describing progressive micro-cracking propagation (i.e. quasi-brittle type damage law) are considered. In the case of time-dependent damage models, the evolution of the micro-crack length during propagation is described through a sub-critical criterion and mixed mode propagation by branching. Using the time dependent damage model including rotational micro-cracks, simulations will be made at three levels: laboratory, tunnel and reservoir scales.

Dans un objectif de réduction de l'utilisation des pesticides, des méthodes alternatives doivent assurer la protection des cultures. Les dégâts causés par les insectes ravageurs sont les manifestations de dynamiques de population dont certains déterminants paysagers sont des éléments structuraux potentiellement relocalisables. Nous explorons ici, par la simulation numérique, la modification du paysage comme stratégie de protection des cultures.Les éléments paysagers peuvent impacter les populations d'insectes par leur influence sur les comportements individuels. Nous avons développé un modèle de population tenant compte de l'approvisionnement individuel. Ce processus résulte (1) de la perception des ressources distantes par l'individu, (2) de la localisation de la population sur une dimension additionnelle quantifiant les réserves énergétiques, et (3) d'une procédure d'optimisation qui définit des mouvements dirigés adaptés.Une heuristique évolutionnaire de modification de paysages a été développée. Elle produit des réarrangements sous contraintes de composition et de structure. Les modèles de paysages et de dynamique de population ont été appliqués à un problème théorique d'aménagement du paysage. La description mécaniste de l'approvisionnement optimal offre des leviers d'action efficaces, basés sur la perturbation des comportements d'approvisionnement. Nos travaux montrent la réactivité de la population aux modifications du paysage, ainsi que la capacité des algorithmes évolutionnaires à proposer des paysages résistants, sous contraintes agronomiques. / Environmental and health issues call for a switch in crop protection towards less chemically driven strategies. Pest damages on crops result of population dynamics that are influenced by landscape features. Those features may be relocated or dismissed to improve landscapes resistance to pest damages. Here we explore numerically the landscape modification as a crop protection strategy.Landscape features can influence population dynamics through their impact on individuals. Therefore, we developed a multiscale population dynamics model accounting for individual foraging. The foraging process results of (1) the perception of distant features used as resources by the individuals, (2) the localisation of the population along an additional dimension quantifying individual energy supply, (3) an optimisation procedure defining adapted directed motionsWe developed an evolutionary heuristic for landscape modification. It is able to rearrange landscapes with respect to compositional and structural constraints. Population and landscape models were applied to a theoretical landscape planning problem. The mechanistic description of the optimal foraging process enables new and efficient levers for crop protection, building on the disruption of the foraging behaviours. This application enlightens the responsiveness of the simulated population to landscape modifications, as well as the ability of evolutionary algorithms to produce resilient landscapes under agronomic productive constraints.