Conditionnement à la peur est un paradigme réussi pour comprendre les substrats neuronaux de l’apprentissage et de l’émotion. Dans cette thèse, je présente deux modèles informatiques des structures du cerveau qui sous-tendent l'acquisition de la peur conditionnée. Le première modèle est utilisé pour enquêter sur l’effet des changements de l’inhibition tonique sur le traitement des informations reçues. On confirme que la diminution de l’inhibition tonique d’une population augmente la réactivité du réseau. Ensuite, le modèle est analysé d’une perspective fonctionnelle et des prédictions qui découlent de cette proposition sont discutées. En outre, je présenterai un modèle systématique, basé sur un type de modèle de conditionnement récemment introduit utilisant des variables latentes. Je propose que l’interaction entre les neurones dans l’amygdale basale code pour l’interface entre ces variables latentes. Le modèle couvre une large gamme d’effets et l’analyse produit un certain nombre de prédictions vérifiables. / Fear conditioning is a successful paradigm for studying neural substrates of emotional learning. In this thesis, two computational models of the underlying neural circuitry are presented. First, the effects of changes in neuronal membrane conductance on input processing are analyzed in a biologically realistic model. We show that changes in tonic inhibitory conductance increase the responsiveness of the network to inputs. Then, the model is analyzed from a functional perspective and predictions that follow from this proposition are discussed. Next, a systems level model is presented based on a recent high-level approach to conditioning. It is proposed that the interaction between fear and extinction neurons in the basal amygdala is a neural substrate of the switching between latent states, allowing the animal to infer causal structure. Important behavioral and physiological results are reproduced and predictions and questions that follow from the main hypothesis are considered.

champ d'application de cette bibliothèque est donc aussi vaste que le domaine de validité de l'équation précédente. La première partie, essentiellement technique, est une description de chacun des éléments de cette bibliothèque. Ecrits en langage FORTRAN IV, ses programmes traitent les aspects suivants : - Introduction et modification des données - Calcul des caractéristiques des écoulements (charges, vitesses, débits) - Dessin automatique des maillages utilisés et des courbes d'isovaleurs des variables définies aux noeuds de ces maillages. La seconde partie présente deux exemples d'utilisation de la bibliothèque ELFI. Le premier, en milieu poreux, se rapporte à 1 'influence de 1 * exploitation de gravières sur l'écoulement des eaux souterraines dans le cas d'une nappe alluviale. Le second, en milieu faillé, concerne le problème des venues d'eau dans une exploitation minière. Outre la représentation de systèmes hydrogéologiques réels, cette bibliothèque permet la constitution de modèles correspondant à des configurations types (série d'abaques) fréquemment rencontrées par 1 'hydrogéologue et faisant intervenir, en particulier, des écoulements tridimensionnels. Un document complémentaire de ce mémoire réunit la description détaillée des fichiers et des ordres d'appel nécessaires à l'utilisation de cette bibliothèque. Ce document doit bien entendu être adapté à chaque configuration informatique.

La compréhension clinique et le traitement des maladies cardiovasculaires est extrêmement complexe. Pour chaque patient, les cardiologues sont confrontés à des difficultés pour déterminer la pathologie, choisir une thérapie ou encore sélectionner les patients susceptibles de bien répondre à un traitement donné. Afin de fournir une aide additionnelle aux cardiologues, de nombreuses équipes de recherche étudient la possibilité de planifier de telles thérapies grâce à des modèles biophysiques du cœur. Ils formulent l'hypothèse que l'on peut combiner les données fonctionnelles et anatomiques afin de créer des modèles cardiaques personnalisés à chaque patient qui auraient le potentiel de prédire les bénéfices des différentes thérapies. Les simulations électromécaniques du cœur sont basées sur des modèles informatiques qui peuvent représenter la géométrie, le mouvement et la propagation d'ondes électriques pendant un cycle cardiaque avec suffisamment de précision. L'intégration d'information anatomique, mécanique et électrophysiologique pour un patient donné est essentielle pour construire ce type de modèle.Dans cette thèse, nous présentons tout d'abord les méthodes de personnalisations géométriques, cinématiques et électrophysiologiques nécessaires à toutes modélisations mécaniques. Nous utilisons ensuite le modèle électromécanique de Bestel-Clément-Sorine qui a déjà prouvé avoir un bon réalisme sans être trop complexe au vu des données disponibles. Nous commençons par détailler la nouvelle implémentation de ce modèle dans une plateforme efficace de simulation médicale ayant l'avantage d'être libre et interactive, puis nous analysons les résultats de la simulation par une étude de sensibilité complète.Dans un deuxième temps, nous étudions la possibilité de personnaliser les paramètres mécaniques de ce modèle à partir d'images médicales (IRM). Pour cela, nous proposons en premier lieu une méthode automatique de calibration qui estime les paramètres mécaniques globaux à partir de courbes de pressions et volumes. Cette technique testée sur 6 volontaires et 2 cas pathologiques nous a permis de faire une étude de spécificité qui consiste à déterminer les paramètres pertinents capables de différencier les cas pathologiques des cas sains.Une fois initialisés à ces valeurs calibrées, les paramètres sont personnalisés localement avec un algorithme d'optimisation plus complexe. Le " Reduced Order Unscented Kalman Filtering " est utilisé pour estimer les contractilités de toutes les zones AHA du ventricule gauche à partir des volumes régionaux extraits des séquences d'images IRM. Cette stratégie de personnalisation a été validée et testée sur plusieurs cas pathologiques et volontaires. Ces différentes contributions ont montré des résultats prometteurs tout au long de cette thèse et certains sont déjà utilisés pour quelques études de recherche.

Des modifications du gène PMP22 (Peripheral Myelin Protein 22) sont responsables de neuropathies du système nerveux périphérique : l’Hypersensibilité à la Pression (HNPP : Hereditary Neuropathy with liability to Pressure Palsy) lorsqu’il est délété, CMT1A (Charcot-Marie-Tooth type 1A) lorsqu’il est dupliqué et CMT1E ou HNPP lors de mutations ponctuelles. Cependant, le rôle de la protéine PMP22 dans ces neuropathies demeure obscur. Dans un premier temps, nous nous sommes intéressés aux partenaires d’interactions potentiels de PMP22 : l’intégrine α6β4 (un récepteur des laminines), pourrait être impliqué dans le CMT1A. Dans cette étude, nous avons utilisé un modèle de rat transgénique portant des copies supplémentaires de PMP22 de souris. Chez ce modèle, nous avons mis en évidence des variations d’expression génique des intégrines, ainsi qu’une mauvaise localisation cellulaire de celles-ci. Ces variations des niveaux d’expression des intégrines sont les témoins d’un retard de la maturation des cellules de Schwann myélinisantes, expliquant la diminution de l’épaisseur des gaines de myéline, observée chez les rats CMT1A. Dans un second temps, nous avons étudié le cas particulier d’un patient sans expression de PMP22, en raison de deux mutations composites sur les deux allèles de PMP22. Nous avons observé que l’absence de PMP22 chez l’homme entraine une absence complète de myéline due à un blocage de l’initialisation de la myélinisation lors de la formation du mésaxone. Dans un troisième temps, nous avons effectué une étude comparative de modèles animaux et de patients atteints de CMT1A / 1E et d’HNPP permettant de valider l’utilisation de tels modèles dans l’étude de ses neuropathies. Enfin, nous nous sommes intéressés d’un point de vue informatique à la structure tridimensionnelle de différentes protéines dont PMP22 grâce à la dynamique moléculaire. Ce modèle tridimensionnel de PMP22 est le point de départ de l’étude des mutations ponctuelles ainsi que des interactions de PMP22 avec son environnement. Grâce aux animaux modèles et à l’étude de patients, nous avons montré le rôle indispensable de PMP22 dans l’initialisation de la myélinisation ainsi que son effet sur les intégrines dans le CMT1A. L’utilisation de modèles informatiques tridimensionnels créés de PMP22 permettra de comprendre les effets des mutations ponctuelles sur sa structure, et ses interactions. / Interactions study of PMP22 protein with integrins in Charcot-Marie-Tooth disease type 1A pathogenesisChanges in the PMP22 gene (Peripheral Myelin Protein 22) are responsible for peripheral nervous system neuropathies: Hereditary Neuropathy with liability to Pressure Palsy (HNPP) when PMP22 is deleted, Charcot Marie Tooth disease subtype 1A (CMT1A) when PMP22 is duplicated and CMT1E or HNPP when point mutations are present on the PMP22 gene. However, the PMP22 role in these neuropathies remains unclear. Firstly, we studied one of the PMP22 interaction partner: the α6β4 integrin (a laminin receptor), which could be involved in the CMT1A. During this study, we used a transgenic rat model carrying supplementary copies of PMP22 gene from mouse. With this model, we showed variations of expression of integrins genes and a mislocalization of integrins proteins. These variations of integrins witnesses a delay in myelinating Schwann cells maturation, explaining the reduction of myelin sheath thickness observed on CMT1A rats. In a second time, we studied the case of a patient without PMP22 expression, carrying compound mutations one the two alleles of PMP22. We showed that the lack of PMP22 on human leads to a complete lack of myelin due to a blocking in mesaxon formation. In a third time, we conducted a comparative study of animal models and patients with CMT1A / 1E and HNPP. This study allows us to validate the use of such models in the study of neuropathies. Finally, thanks to computational tools we studied the three-dimensional structure of different protein, including PMP22 by using molecular dynamics. This three-dimensional model is the starting point of point mutations study as well as PMP22 interactions with its environments. With animal models and through the study of patients, we have demonstrated the indispensable role of PMP22 in myelin initiation as well as, its effect on integrins expression in CMT1A. The use of PMP22 three-dimensional computational model will help us to understand effects of point mutation on PMP22 structure and interactions.

La compréhension clinique et le traitement des maladies cardiovasculaires est extrêmement complexe. Pour chaque patient, les cardiologues sont confrontés à des difficultés pour déterminer la pathologie, choisir une thérapie ou encore sélectionner les patients susceptibles de bien répondre à un traitement donné. Afin de fournir une aide additionnelle aux cardiologues, de nombreuses équipes de recherche étudient la possibilité de planifier de telles thérapies grâce à des modèles biophysiques du cœur. Ils formulent l'hypothèse que l'on peut combiner les données fonctionnelles et anatomiques afin de créer des modèles cardiaques personnalisés à chaque patient qui auraient le potentiel de prédire les bénéfices des différentes thérapies. Les simulations électromécaniques du cœur sont basées sur des modèles informatiques qui peuvent représenter la géométrie, le mouvement et la propagation d'ondes électriques pendant un cycle cardiaque avec suffisamment de précision. L'intégration d'information anatomique, mécanique et électrophysiologique pour un patient donné est essentielle pour construire ce type de modèle.Dans cette thèse, nous présentons tout d'abord les méthodes de personnalisations géométriques, cinématiques et électrophysiologiques nécessaires à toutes modélisations mécaniques. Nous utilisons ensuite le modèle électromécanique de Bestel-Clément-Sorine qui a déjà prouvé avoir un bon réalisme sans être trop complexe au vu des données disponibles. Nous commençons par détailler la nouvelle implémentation de ce modèle dans une plateforme efficace de simulation médicale ayant l'avantage d'être libre et interactive, puis nous analysons les résultats de la simulation par une étude de sensibilité complète.Dans un deuxième temps, nous étudions la possibilité de personnaliser les paramètres mécaniques de ce modèle à partir d'images médicales (IRM). Pour cela, nous proposons en premier lieu une méthode automatique de calibration qui estime les paramètres mécaniques globaux à partir de courbes de pressions et volumes. Cette technique testée sur 6 volontaires et 2 cas pathologiques nous a permis de faire une étude de spécificité qui consiste à déterminer les paramètres pertinents capables de différencier les cas pathologiques des cas sains.Une fois initialisés à ces valeurs calibrées, les paramètres sont personnalisés localement avec un algorithme d'optimisation plus complexe. Le « Reduced Order Unscented Kalman Filtering » est utilisé pour estimer les contractilités de toutes les zones AHA du ventricule gauche à partir des volumes régionaux extraits des séquences d'images IRM. Cette stratégie de personnalisation a été validée et testée sur plusieurs cas pathologiques et volontaires. Ces différentes contributions ont montré des résultats prometteurs tout au long de cette thèse et certains sont déjà utilisés pour quelques études de recherche. / The clinical understanding and treatment of cardiovascular diseases is highly complex. For each patient, cardiologists face issues in determining the pathology, choosing a therapy or selecting suitable patients for the therapy. In order to provide additional guidance to cardiologists, many research groups are investigating the possibility to plan such therapies based on biophysical models of the heart. The hypothesis is that one may combine anatomical and functional data to build patient-specific cardiac models that could have the potential to predict the benefits of different therapies. Cardiac electromechanical simulations are based on computational models that can represent the heart geometry, motion and electrophysiology patterns during a cardiac cycle with sufficient accuracy. Integration of anatomical, mechanical and electrophysiological information for a given subject is essential to build such models.In this thesis, we first introduce the geometry, kinematics and electrophysiology personalizations that are necessary inputs to mechanical modeling. We propose to use the Bestel-Cl'ement-Sorine electromechanical model of the heart, which is sufficiently accurate without being over-parametrized for the available data. We start by presenting a new implementation of this model in an efficient opensource framework for interactive medical simulation and we analyze the resulting simulations through a complete sensitivity analysis.In a second step, the goal is to personalize the mechanical parameters from medical images (MRI data). To this end, we first propose an automatic calibration algorithm that estimates global mechanical parameters from volume and pressure curves. This technique was tested on 7 volunteers and 2 heart failure cases and allowed to perform a preliminary specificity study that intends to determine the relevant parameters able to differentiate the pathological cases from the control cases.Once initialized with the calibrated values, the parameters are then locally personalized with a more complex optimization algorithm. Reduced Order Unscented Kalman Filtering is used to estimate the contractilities on all of the AHA zones of the Left Ventricle, matching the regional volumes extracted from cine MRI data. This personalization strategy was validated and tested on several pathological and healthy cases. These contributions have led to promising results through this thesis and some are already used for various research studies.

La thèse est structurée en trois chapitres. Dans le premier, on rappelle brièvement le contexte historique et législatif de l'évolution de l'assainissement en Europe. Le deuxième chapitre s'intéresse à l'étude diagnostic du réseau d'assainissement de la ville de Saint-Etienne, tant sur les plans quantitatif que qualitatif. Pour cela, nous développons un logiciel de simulation dont l'originalité réside dans sa structure informatique. Nous effectuons une étude de la pluviométrie locale, nous localisons les points de dysfonctionnement par la simulation et nous analysons l'impact des rejets stéphanois en phosphore sur la qualité des eaux de la Loire. Le troisième chapitre essaie par des développements méthodologiques de répondre aux questions soulevées précédemment, tout d'abord par une étude visant à caractériser les structures spatiales des pluies en utilisant des techniques issues de l'analyse de données, ensuite par une étude de sensibilité du modèle d'écoulement à la distribution spatiale de la pluie, au coefficient de ruissellement et au temps de réponse de bassins versants élémentaires, enfin par une dernière étude visant à déterminer la localisation optimale du réseau de pluviographes, en utilisant la technique du krigeage, en y intégrant la structure spatiale de la pluie et les coefficients de ruissellement des bassins versants

La fibrillation auriculaire, l'arythmie la plus fréquente en clinique, affecte 2.3 millions de patients en Amérique du Nord. Pour en étudier les mécanismes et les thérapies potentielles, des modèles animaux de fibrillation auriculaire ont été développés. La cartographie électrique épicardique à haute densité est une technique expérimentale bien établie pour suivre in vivo l'activité des oreillettes en réponse à une stimulation électrique, à du remodelage, à des arythmies ou à une modulation du système nerveux autonome. Dans les régions qui ne sont pas accessibles par cartographie épicardique, la cartographie endocardique sans contact réalisée à l'aide d'un cathéter en forme de ballon pourrait apporter une description plus complète de l'activité auriculaire. Dans cette étude, une expérience chez le chien a été conçue et analysée. Une reconstruction électro-anatomique, une cartographie épicardique (103 électrodes), une cartographie endocardique sans contact (2048 électrodes virtuelles calculées à partir un cathéter en forme de ballon avec 64 canaux) et des enregistrements endocardiques avec contact direct ont été réalisés simultanément. Les systèmes d'enregistrement ont été également simulés dans un modèle mathématique d'une oreillette droite de chien. Dans les simulations et les expériences (après la suppression du nœud atrio-ventriculaire), des cartes d'activation ont été calculées pendant le rythme sinusal. La repolarisation a été évaluée en mesurant l'aire sous l'onde T auriculaire (ATa) qui est un marqueur de gradient de repolarisation. Les résultats montrent un coefficient de corrélation épicardique-endocardique de 0.8 (expérience) and 0.96 (simulation) entre les cartes d'activation, et un coefficient de corrélation de 0.57 (expérience) and 0.92 (simulation) entre les valeurs de ATa. La cartographie endocardique sans contact apparait comme un instrument expérimental utile pour extraire de l'information en dehors des régions couvertes par les plaques d'enregistrement épicardique. / Atrial fibrillation is the most common clinical arrhythmia currently affecting 2.3 million patients in North America. To study its mechanisms and potential therapies, animal models of atrial fibrillation have been developed. Epicardial high-density electrical mapping is a well-established experimental instrument to monitor in vivo the activity of the atria in response to pacing, remodeling, arrhythmias and modulation of the autonomic nervous system. In regions that are not accessible by epicardial mapping, noncontact endocardial mapping performed through a balloon catheter may provide a more comprehensive description of atrial activity. In this study, a dog experiment was designed and analyzed in which electroanatomical reconstruction, epicardial mapping (103 electrodes), noncontact endocardial mapping (2048 virtual electrodes computed from a 64-channel balloon catheter), and direct-contact endocardial catheter recordings were simultaneously performed. The recording system was also simulated in a computer model of the canine right atrium. For simulations and experiments (after atrio-ventricular node suppression), activation maps were computed during sinus rhythm. Repolarization was assessed by measuring the area under the atrial T wave (ATa), a marker of repolarization gradients. Results showed an epicardial endocardial correlation coefficient of 0.8 (experiment) and 0.96 (simulation) between activation times, and a correlation coefficient of 0.57 (experiment) and 0.92 (simulation) between ATa values. Noncontact mapping appears to be a valuable experimental device to retrieve information outside the regions covered by epicardial recording plaques.

La fibrillation auriculaire, l'arythmie la plus fréquente en clinique, affecte 2.3 millions de patients en Amérique du Nord. Pour en étudier les mécanismes et les thérapies potentielles, des modèles animaux de fibrillation auriculaire ont été développés. La cartographie électrique épicardique à haute densité est une technique expérimentale bien établie pour suivre in vivo l'activité des oreillettes en réponse à une stimulation électrique, à du remodelage, à des arythmies ou à une modulation du système nerveux autonome. Dans les régions qui ne sont pas accessibles par cartographie épicardique, la cartographie endocardique sans contact réalisée à l'aide d'un cathéter en forme de ballon pourrait apporter une description plus complète de l'activité auriculaire. Dans cette étude, une expérience chez le chien a été conçue et analysée. Une reconstruction électro-anatomique, une cartographie épicardique (103 électrodes), une cartographie endocardique sans contact (2048 électrodes virtuelles calculées à partir un cathéter en forme de ballon avec 64 canaux) et des enregistrements endocardiques avec contact direct ont été réalisés simultanément. Les systèmes d'enregistrement ont été également simulés dans un modèle mathématique d'une oreillette droite de chien. Dans les simulations et les expériences (après la suppression du nœud atrio-ventriculaire), des cartes d'activation ont été calculées pendant le rythme sinusal. La repolarisation a été évaluée en mesurant l'aire sous l'onde T auriculaire (ATa) qui est un marqueur de gradient de repolarisation. Les résultats montrent un coefficient de corrélation épicardique-endocardique de 0.8 (expérience) and 0.96 (simulation) entre les cartes d'activation, et un coefficient de corrélation de 0.57 (expérience) and 0.92 (simulation) entre les valeurs de ATa. La cartographie endocardique sans contact apparait comme un instrument expérimental utile pour extraire de l'information en dehors des régions couvertes par les plaques d'enregistrement épicardique. / Atrial fibrillation is the most common clinical arrhythmia currently affecting 2.3 million patients in North America. To study its mechanisms and potential therapies, animal models of atrial fibrillation have been developed. Epicardial high-density electrical mapping is a well-established experimental instrument to monitor in vivo the activity of the atria in response to pacing, remodeling, arrhythmias and modulation of the autonomic nervous system. In regions that are not accessible by epicardial mapping, noncontact endocardial mapping performed through a balloon catheter may provide a more comprehensive description of atrial activity. In this study, a dog experiment was designed and analyzed in which electroanatomical reconstruction, epicardial mapping (103 electrodes), noncontact endocardial mapping (2048 virtual electrodes computed from a 64-channel balloon catheter), and direct-contact endocardial catheter recordings were simultaneously performed. The recording system was also simulated in a computer model of the canine right atrium. For simulations and experiments (after atrio-ventricular node suppression), activation maps were computed during sinus rhythm. Repolarization was assessed by measuring the area under the atrial T wave (ATa), a marker of repolarization gradients. Results showed an epicardial endocardial correlation coefficient of 0.8 (experiment) and 0.96 (simulation) between activation times, and a correlation coefficient of 0.57 (experiment) and 0.92 (simulation) between ATa values. Noncontact mapping appears to be a valuable experimental device to retrieve information outside the regions covered by epicardial recording plaques.

Cotutelle internationale avec l'université de Caen / La complexité des patients de réanimation justifie le recours à des systèmes d’aide à la décision thérapeutique. Ces systèmes rassemblent des protocoles automatisés de prise en charge permettant le respect des recommandations et des simulateurs physiologiques ou patients virtuels, utilisables pour personnaliser de façon sécuritaire les prises en charge. Ces dispositifs fonctionnant à partir d’algorithmes et d’équations mathématiques ne peuvent être développés qu’à partir d’un grand nombre de données de patients. Le principal objectif de cette thèse était la mise en place d’une base de données haute résolution automatiquement collectée de patients de réanimation pédiatrique dont le but sera de servir au développement et à la validation d’un simulateur physiologique : SimulResp©. Ce travail présente l’ensemble du processus de mise en place de la base de données, du concept jusqu’à son utilisation. / The complexity of the patients in the intensive care unit requires the use of clinical decision support systems. These systems bring together automated management protocols that enable adherence to guidelines and virtual physiological or patient simulators that can be used to safely customize management. These devices operating from algorithms and mathematical equations can only be developed from a large number of patients’ data. The main objective of the work was the elaboration of a high resolution database automatically collected from critically ill children. This database will be used to develop and validate a physiological simulator called SimulResp© . This manuscript presents the whole process of setting up the database from concept to use.

L'ablation par radiofréquence (ARF) de tumeurs abdominales est rendue difficile par l’influence des vaisseaux sanguins et les variations de la conductivité thermique, compliquant la planification spécifique à un patient donné. En fournissant des outils prédictifs, les modèles biophysiques pourraient aider les cliniciens à planifier et guider efficacement la procédure. Nous introduisons un modèle mathématique détaillé des mécanismes impliqués dans l’ARF des tumeurs du foie comme la diffusion de la chaleur et la nécrose cellulaire. Il simule l’étendue de l’ablation à partir d’images médicales, d’après lesquelles des modèles personnalisés du foie, des vaisseaux visibles et des tumeurs sont segmentés. Dans cette thèse, une nouvelle approche pour résoudre ces équations basée sur la méthode de Lattice Boltzmann est introduite. Le modèle est d’abord évalué sur des données de patients qui ont subi une ARF de tumeurs du foie. Ensuite, un protocole expérimental combinant des images multi-modales, anatomiques et fonctionnelles pré- et post-opératoires, ainsi que le suivi de la température et de la puissance délivrée pendant l'intervention est présenté. Il permet une validation totale du modèle qui considère des données les plus complètes possibles. Enfin, nous estimons automatiquement des paramètres personnalisés pour mieux prédire l'étendu de l’ablation. Cette stratégie a été validée sur 7 ablations dans 3 cas cliniques. A partir de l'étude préclinique, la personnalisation est améliorée en comparant les simulations avec les mesures faites durant la procédure. Ces contributions ont abouti à des résultats prometteurs, et ouvrent de nouvelles perspectives pour planifier et guider l’ARF. / The outcome of radiofrequency ablation (RFA) of abdominal tumors is challenged by the presence of blood vessels and time-varying thermal conductivity, which make patient-specific planning extremely difficult. By providing predictive tools, biophysical models may help clinicians to plan and guide the procedure for an effective treatment. We introduce a detailed computational model of the biophysical mechanisms involved in RFA of hepatic tumors such as heat diffusion and cellular necrosis. It simulates the extent of ablated tissue based on medical images, from which patient-specific models of the liver, visible vessels and tumors are segmented. In this thesis, a new approach for solving these partial differential equations based on the Lattice Boltzmann Method is introduced. The model is first evaluated against clinical data of patients who underwent RFA of liver tumors. Then, a comprehensive pre-clinical experiment that combines multi-modal, pre- and post-operative anatomical and functional images, as well as the interventional monitoring of the temperature and delivered power is presented. This enables an end-to-end validation framework that considers the most comprehensive data set for model validation. Then, we automatically estimate patient-specific parameters to better predict the ablated tissue. This personalization strategy has been validated on 7 ablations from 3 clinical cases. From the pre-clinical study, we can go further in the personalization by comparing the simulated temperature and delivered power with the actual measurements during the procedure. These contributions have led to promising results, and open new perspectives in RFA guidance and planning.