Le travail de cette thèse porte sur le développement de modèles neutroniques innovants pour le couplage avec la thermohydraulique, associant précision et temps de calcul raisonnable. Un des cas d'application principaux étant le réacteur à sel fondu, à spectre neutronique rapide et en cycle thorium MSFR (Molten Salt Fast Reactor), réacteur de 4ème génération à combustible liquide circulant, la prise en compte du mouvement des précurseurs de neutrons retardés et des phénomènes associés est nécessaire. Les études de conception de ce type de réacteur ont été le point de départ de ces développements, via le besoin d'une représentation multiphysique adaptée pour l'obtention d'une image globale et la réalisation d'études de transitoire.Dans un premier temps un couplage stationnaire a été développé, associant un modèle neutronique basé sur une approche stochastique, et un code de CFD (Computational Fluid Dynamics) résolvant les équations de Navier Stokes des écoulements turbulents ainsi que le transport des précurseurs de neutrons retardés. Ce modèle neutronique intègre l'effet lié au transport de ces précurseurs par une reconstruction de la gerbe prompte qu'ils génèrent. Cette approche dite par gerbe considère le réacteur critique comme un système sous-critique prompt amplifiant la source de neutrons retardés.Dans un second temps, un modèle neutronique basé sur une version temporelle des matrices de fission (Transient Fission Matrix ou TFM) a été développé afin de réaliser des études de transitoires. Le modèle TFM permet, en un premier calcul des matrices avec un code stochastique (MCNP, SERPENT), de réaliser une caractérisation de l'ensemble de la réponse neutronique spatiale et temporelle du réacteur avec une précision proche de celle du calcul Monte Carlo. Dans un second temps cette information est utilisée pour les calculs de transitoires tout en gardant un temps de calcul réduit. Le modèle TFM, utilisable pour différents types de systèmes, permet également le calcul de paramètres cinétiques effectifs tels que la fraction effective de neutrons retardés ou le temps de génération effectif. Différents cas d'application ont été utilisés afin de vérifier et d'illustrer cette approche sur des calculs temporels ou de paramètres cinétiques.Enfin le modèle TFM a été implémenté dans le code de thermohydraulique OpenFOAM. Ce couplage a été testé sur un benchmark numérique à géométrie simplifiée, puis des calculs sur le MSFR ont été réalisés, pour des transitoires normaux (suivis de charge) ou accidentels (insertions de réactivité, sur-refroidissements). / In this PhD thesis, we describe the development of innovative neutronic models for their coupling with thermalhydraulics such that they combine precision and reasonable computational times. One of the main cases where this method is applied is the Molten Salt Fast Reactor (MSFR) whose combines a fast neutron spectrum with a thorium cycle. In this fourth generation reactor, the motion of the delayed neutron precursors and the associated phenomena have to be taken into account due to the liquid fuel circulation. The starting point for these developments was the preliminary design of this type of system where a dedicated multi-physical representation was needed to study the reactor performance in steady and transient conditions.As a first step, a stationary coupling was developed. A neutronic model based on a stochastic approach was associated to a CFD (Computational Fluid Dynamics) code to solve the Navier Stokes equations for turbulent flows and the transport of the delayed neutron precursors. The impact of this precursor motion is taken into account by reconstructing the prompt shower that they generate. This approach, called by shower, views the critical reactor as a prompt subcritical reactor that amplifies a source of delayed neutrons.A second step consisted in developing a neutronic model based on a time dependent version of the fission matrices (Transient Fission Matrix or TFM) so as to enable reactor transient studies. With the TFM model, an initial computation of the matrices with a stochastic code (MCNP, SERPENT) allows the characterization of the global spatial and time dependent neutronic response of the reactor with a precision close to that of a Monte Carlo calculation. The information thus obtained is then used to calculate transients, while retaining the advantage of reduced computational time. The TFM model, which can be used for various system concepts, also allows the evaluation of effective kinetic parameters such as the effective fraction of delayed neutrons or the effective generation time. The method was applied to various cases in order to verify it and demonstrate the approach for time dependent or kinetic parameter calculations.Finally, the TFM model was integrated in the OpenFOAM thermalhydraulic code. The coupling was first tested on a simple geometry numerical benchmark. Subsequently, it was applied to the MSFR to calculate normal (load-following) and accidental (reactivity insertion, over-cooling) transients.