Ce travail doctoral porte sur la synthèse expérimentale de champs de pression acoustique aléatoires sur des surfaces planes. L'objectif principal de cette reproduction est de pouvoir tester le comportement vibroacoustique de panneaux plans sous l'action de deux champs de pression aléatoires : le champ acoustique diffus et la couche limite turbulente. Une étude théorique et numérique est tout d'abord présentée, qui aborde la formulation de ce problème de reproduction sur la base de la théorie de l'holographie acoustique plane de champ proche, et qui permet de vérifier que celle-ci permet une résolution de ce problème. L'influence de nombreux paramètres est également étudiée et des points critiques dans la mise en oeuvre expérimentale de l'approche sont identifiés, parmi eux la densité élevée de sources nécessaire pour la reproduction d'une couche limite turbulente. Une première mise en oeuvre expérimentale est détaillée. Trois approches, celle basée sur l'holographie, une seconde sur la synthèse de champs sonores (ou Wave Field Synthesis ) ainsi qu'une troisième basée sur les moindres carrés sont utilisées pour le calcul des amplitudes complexes des sources de reproduction. Le concept d'antenne synthétique, pour lequel un petit élément d'un réseau est déplacé pour créer un réseau de grandes dimensions en post-traitement, sera utilisé pour résoudre une majorité des points critiques dans la mise en oeuvre. Des indicateurs vibroacoustiques fréquentiels tels que le Transmission Loss (TL) sont obtenus pour un panneau d'aluminium simplement supporté, et la comparaison avec des résultats de calculs numériques montrent que ces approches permettent l'estimation de ce TL sous les deux champs aléatoires d'intérêt. La possibilité d'estimer le TL dans le domaine des nombres d'onde est également illustrée. Une seconde mise en oeuvre expérimentale est réalisée sous la forme d'une comparaison entre la méthode décrite ci-dessus et la méthode des chambres couplées, couramment utilisée pour l'estimation du TL. Pour un panneau composite aéronautique représentatif et pour un champ acoustique diffus, les écarts entre les résultats obtenus par les deux méthodes sont jugés acceptables, et permettent une seconde validation. Parallèlement, une autre application a concerné la reproduction d'un champ acoustique diffus sur la surface d'un matériau poreux, afin d'en estimer le coefficient d'absorption sous ce champ excitateur. L'approche théorique de cette estimation est décrite et validée expérimentalement sur une mousse de mélamine, prouvant qu'elle est réalisable sous des champs acoustiques synthétisés et montrant un potentiel pour des mesures in situ. Suite à la conclusion de ce mémoire qui rappelle les résultats importants obtenus dans le cadre de ce travail doctoral, une possibilité de mesure des fluctuations de pression pariétale liées à la couche limite turbulente est finalement présentée en annexe.

L'objectif principal de cette thèse est dirigé vers l'établissement et l'explication des relations entre les conditions de synthèse des couches minces de CZTSSe, ses propriétés physiques et les performances des dispositifs photovoltaïques. Pour faire face à cette tâche la première approche était de comprendre le mécanisme de formation de la matière par rapport aux conditions de croissance du matériau. Le CZTSSe est synthétisé par un processus de sélénisation en deux étapes, où une première étape de dépôt par PVD de précurseurs est nécessaire, suivie d'une seconde étape de recuit sous atmosphère de sélénium. Différents ordres d'empilement de précurseurs ont été étudiés afin de comprendre la séquence de réactions qui, à partir de leur dépôt, conduise à la couche finale de CZTSSe. Cette étude, fait en plusieurs étapes, a nécessité de un effort important sur la caractérisation du matériau à chaque étape de la synthèse. Le résultat a montré que dans le cas du procédé en deux étapes, le matériau final est indépendant du dépôt de précurseurs. Les possibles implications bénéfiques en raison de l'incorporation de sodium dans le CZTSSe sont également décrites. Cette étude est réalisée en synthétisant la couche de CZTSSe sur différents substrats contenant diffèrent taux de sodium: de cette manière, pendant la synthèse, le sodium migre de substrats vers l'absorbeur. Après quantification du Na dans le CZTSSe juste après la croissance, le matériau est caractérise afin d'évaluer sa qualité. Ensuite il est employé dans une cellule solaire complète pour vérifier ses propriétés photovoltaïques. Les résultats ont montré que, comme pour la technologie CIGS, le sodium est bénéfique pour le CZTSSe, permettant l'augmentation de la tension à circuit ouvert et le rendement de cellule. Le molybdène est le contact arrière le plus utilisé pour les cellules solaires à base CZTSSe. Cependant, il a été suggéré récemment que le Mo n'est pas stable à l'interface avec le CZTSSe. En outre, à ma connaissance, aucune étude expérimentale n'a été effectuée à ce jour pour tester si les cellules solaires construites sur un autre contact arrière pourraient présenter de meilleures propriétés photovoltaïques. A cet effet, divers métaux (Au, W, Pd, Pt et Ni) sont déposées sur le dessus de Mo et testés comme contacts arrières dans les cellules solaire à base de CZTSSe. Il est démontré qu'il est possible synthétiser de films minces de CZTSSe de qualité quand le tungstène, l'or et le platine sont employé comme contacts arrière. Il est démontré que les contacts en W et Au permettent d'augmenter le courant photogénéré, mais aussi que le Mo reste le meilleur contact arrière en termes d'efficacité de conversion. Les effets de la variation du rapport [S]/([S]+[Se]) sur les performances des cellules solaires à base CZTSSe ont été étudiés. Cette étude a été faite par simulations des cellules solaires à base de CZTSSe, où le taux de chalcogènes dans l'absorbeur est varié, avec l'objective de trouver la composition optimale de l'absorbeur. Deux types d'approche différente ont été étudiés: la variation linéaire du rapport des chalcogènes, et une variation parabolique. Les simulations conduisent à un rendement de 16,5% (avec une tension en circuit ouvert de 0,56 V, courant de court-circuit de 37,0 mA/cm2 et un facteur de forme de 79,0%) lorsque la teneur en soufre est diminué linéairement à partir du contact arrière en direction de la couche tampon. Sur la base de ces résultats, nous proposons que l'ingénierie de bande interdite sur la base de la variation du taux [S]/([S]+[Se]) dans l'absorbeur est un outil puissant qui permet d'augmenter les performances des cellules solaires à base CZTSSe sans changer la qualité de l'absorbeur en lui-même.

Le développement de couches minces dopées terres rares a suscité un regain d'intérêt au cours des dernières années. Dans ce mémoire nous présentons l'élaboration des couches minces d'yttria (Y2O3), d'alumine (Al2O3) et les couches mixtes Y2O3-Al2O3 dopées erbium. La technique utilisée est le procédé de dépôt chimique en phase vapeur à partir de précurseurs organométallique (MOCVD) assisté par aérosol. Un dispositif d'irradiation UV est appliqué afin d'assister le processus de réaction avec une modification de l'hygrométrie de l'air vecteur. Les meilleures propriétés sont obtenues pour les couches déposées sous une forte humidité de l'air vecteur et avec l'assistance UV. Dans ces conditions les couches d'yttria présentent une faible vitesse de croissance, une faible contamination organique et une bonne cristallinité dans la phase cubique de l'yttria. Plusieurs phénomènes d'Up-conversion ont été mis en évidence dans les spectres de fluorescence visible de l'erbium dans l'yttria. Une durée de vie du niveau 4I13/2 de l'erbium de 3.07 ms a été mesurée pour ce matériau après recuit à 800°C. Cette valeur est supérieure à celle obtenue pour l'échantillon déposé sous une faible humidité de l'air et sans l'assistance UV après recuit à 1000°C. Les couches d'alumine déposées dans les conditions optimales présentent des vitesses de croissance élevées et se caractérisent par une grande stabilité thermique, permettant l'élimination complète des impuretés tout en restant amorphe. Enfin, l'étude du système Y2O3-Al2O3 montre que les conditions de dépôt jouent un rôle important sur la composition et les propriétés physico-chimiques des dépôts.

Le Nafion dans la couche catalytique des PEMFCs peut former un revêtement de film ultra-mince enrobant la surface du catalyseur et ses supports. La morphologie du Nafion se révèle être très sensibles à la nature du matériau sur lequel le film est déposé, et en particulier le caractère hydrophobe/phile de ces matériaux. Notre travail consiste à effectuer une enquête complète sur les effets hydrophiles du substrat sur les propriétés physiques du film ultra-mince de Nafion à des niveaux différents d'hydratation. Par conséquent, nous étudions selon un cadre unique une variété d'environnements spécifiques de la couche du catalyseur de la PEMFC, pouvant aller d'un substrat hydrophobe (carbone) à hydrophile (platine). La méthode numérique choisie pour ce travail est une simulation par Dynamique Moléculaire classique. Les configurations de films ultra-minces correctement thermalisés ont été décrites en détail en fonction de leurs propriétés structurales et dynamiques.

Le paradigme d'Internet des objets (IdO) envisage d'augmenter l'Internet actuel avec un grand nombre de dispositifs intelligents communicants. Les réseaux de capteurs sans fil (RCsF) déploient les appareils fonctionnant sur les approvisionnements énergétiques maigres et mesurent les phénomènes environnementaux (comme la température, la radioactivité, ou CO2). Les applications populaires de réseaux de capteurs sans fil comprennent la surveillance, la télémétrie et la prévention des catastrophes naturelles. Les majeurs défis de réseaux de capteurs consistent à: comment obtenir l'efficacité énergétique? Comment surmonter les déficiences de support sans fil? et comment les faire fonctionner d'une manière auto-organisée? L'intégration de réseaux de capteurs dans IoD s'appuiera sur un ensemble de standards ouverts qui s'efforcent d'offrir le passage à l'échelle, la fiabilité dans une variété des conditions, et les scénarios d'exploitation. Néanmoins, les standards dans leur état actuel ont des problèmes d'interopérabilité et peuvent bénéficier de nouvelles améliorations. Les contributions de la thèse sont les suivantes: * Nous avons mené une analyse expérimentale et identifié les caractéristiques d'un environnement de réseau de capteurs. Notre analyse a porté sur la modélisation des liens et la dynamique du réseau, ainsi que la corrélation des liens avec les paramètres environnementaux. L'étude analytique nous a permis d'identifier les principales faiblesses de l'environnement du RCsF et a permis aussi d'obtenir une meilleure compréhension de la dynamique de lien radio et le voisinage des nœuds. * Nous avons confrontés la question de l'interopérabilité des deux principaux standards - IEEE 802.15.4 dans la couche de contrôle d'accès au medium, et RPL dans la couche de transport. Nous proposons d'adapter la structure originale de cluster tree et de construire un framework élégant, qui permettra le fonctionnement conjoint de l'IEEE 802.15.4 et RPL sans l'apparition des collisions. En outre, nous évaluons à travers de nombreuses simulations les deux algorithmes distribués d'attribution de slots. Ces algorithmes permettront d'atteindre un fonctionnement presque sans collisions. * Nous proposons une solution inter-couche distribuée dans le but de construire une topologie afin que IEEE 802.15.4 et RPL puissent fonctionner ensemble. Notre algorithme auto-organisé obtient une structure topologique obéissante aux recommandations globales, cela, par la seule utilisation des métriques mesurées localement. De nombreuses simulations démontrent les avantages de la structure résultante en termes de convergence de temps, de stabilité, et d'efficacité énergétique dans le long terme. Les simulations démontrent aussi un impact positif sur les performances de routage. * Nous proposons un ensemble de nouveaux mécanismes qui permettront d'améliorer les performances de RPL et d'utiliser la Qualité de Service pour gérer les problèmes de trafic sensible. Notre extension de routage opportuniste et multi-chemin contribue à améliorer la livraison des paquets avant une date limite, tout en minimisant le surcout et la consommation d'énergie par rapport à la version de base de RPL.

Face à la volonté d'intégrer une quantité toujours plus importante de nouveaux services au sein des terminaux mobiles de nouvelle génération et afin de répondre à leurs contraintes d'encombrement, des nouveaux concepts d'antennes intelligentes font l'objet de nombreuses recherches. Parmi les solutions proposées dans la littérature, la technique consistant à charger l'antenne par un matériau aux propriétés commandables apparaît particulièrement intéressante puisque elle cumule les effets de miniaturisation et d'agilité. Le travail de cette thèse concerne l'intégration des films minces La-Ti-O-N et BST dans des antennes miniatures et agiles en fréquence. Pour cela, une étude systématique des propriétés diélectriques (constante, tangente de pertes et accordabilité) des films La-Ti-O-N a été réalisée en basses et hautes fréquences. Ces propriétés sont fonction des caractéristiques structurales des films, elles-mêmes issues de la nature du substrat utilisé et des conditions de dépôt. En parallèle, une étude sur l'intégration de ces matériaux dans des structures rayonnantes pour atteindre l'agilité souhaitée a été menée. De premiers démonstrateurs d'antennes miniatures et reconfigurables à base des films minces La-Ti-O-N et BST ont été réalisés et caractérisés.

L'étude effectuée au cours de cette thèse s'inscrit dans le contexte du contrôle des écoulements en turbomachines. Suite à l'intérêt avéré de l'aspiration de la couche limite en compresseurs axiaux pour en améliorer les performances, elle vise à évaluer si cette méthodologie de contrôle est pertinente en géométrie radiale. Elle s'appuie pour cela sur un compresseur centrifuge transsonique Turbomeca. Le fonctionnement du compresseur est d'abord analysé sur la base de simulations numériques stationnaires et instationnaires chorochroniques, réalisées grâce au code de calcul elsA développé par l'Onera et en utilisant le modèle de turbulence k-l de Smith. La validité des résultats numériques est vérifiée par comparaison avec les données expérimentales disponibles : performances mono-dimensionnelles et signaux de pressions instationnaires. Ces données sont issues à la fois de la caractérisation expérimentale du compresseur effectuée par Turbomeca et de campagnes de mesures menées a l'ISAE, ou le compresseur est monté sur un banc de recherche dédié à l'étude des interactions entre le rouet et le diffuseur radial. L'analyse détaillée des résultats numériques révèle le développement d'un décollement de coin dans le diffuseur, entre la paroi moyeu et la face en dépression des aubages, lorsque le point de fonctionnement se rapproche du pompage. L'étude topologique du spectre du frottement pariétal permet d'en préciser le lieu et la structure, et de localiser les points singuliers ainsi que la ligne de séparation principale. En stationnaire, la croissance du décollement aboutit au décrochage du diffuseur, et compromet les performances de l'étage de compression et finalement la stabilité numérique. L'exploitation des résultats instationnaires permet de préciser le comportement temporel de l'écoulement : l'étendue du décollement fluctue sous l'effet du défilement des ondes de pression générées par l'interaction rouet-diffuseur, mais les trajectoires des particules fluides impliquées dans le décollement sont en accord avec les lignes de courant du champ de l'écoulement moyenné temporellement. Le décollement instationnaire est fixe, et sa topologie correspond à celle prévue par les simulations stationnaires. En particulier, le lieu du décollement est prévu de façon similaire par les modèles stationnaires et instationnaires. Une stratégie d'aspiration est ensuite développée sur la base de ces observations : une fente de prélèvement est positionnée au voisinage du col de séparation principal, dont le lieu correspond à celui du maximum du gradient de pression adverse stationnaire. Cette stratégie est implémentée au sein des modèles numériques stationnaires et instationnaires. En stationnaire, un prélèvement de 1 % de la valeur totale du débit traversant le rouet permet un contrôle total du décollement dans le diffuseur, et conduit à une augmentation significative de la plage de stabilité numérique. La réduction du débit de prélèvement à 0,3 % du débit total permet un contrôle partiel du décollement, et conduit également à une augmentation de la plage de stabilité numérique. En instationnaire, le décollement de coin initial est contrôlé. Mais les résultats mettent en évidence le rôle majeur joue par le défilement des ondes de pression le long des aubages du diffuseur. Celles-ci se renforcent au passage du col du diffuseur, et engendrent l'existence d'un important gradient de pression adverse instantané. Ce maximum du gradient de pression instantané conduit au développement d'un nouveau décollement de coin, en aval de la fente de prélèvement, et la plage de stabilité des calculs instationnaires n'est pas augmentée. Ces résultats mettent en évidence la possibilité d'agir sur les décollements se produisant dans les diffuseurs radiaux à l'aide de la technique d'aspiration de la couche limite. [...]

Cette étude est consacrée à l'analyse des écoulements réactifs supersoniques cisailléset, plus particulièrement, des couches de mélange compressibles pouvant se développerdans les moteurs ramjet et scramjet. Des méthodes numériques appropriées ont été implémentéeset vérifiées pour aboutir au développement d'un code de calcul numériquemassivement parallèle, appelé CREAMS (compressible reactive multi-species solver). Cedernier a été spécialement conçu pour conduire des simulations numériques haute précision(simulations numériques directes ou DNS) de ce type d'écoulements. Une attentionparticulière a été portée à la description des termes de transport moléculaire et des termessources chimiques de façon à considérer la description physique la plus fidèle possible desmélanges des gaz réactifs à haute vitesse, au sein desquelles les temps caractéristiqueschimiques et de mélange aux petites échelles sont susceptibles d'être du même ordre degrandeur. Les simulations des couches de mélange bidimensionnelles et tridimensionnelles,inertes et réactives, confirment l'importance des effets associés à la compressibilité et autaux de dégagement de chaleur. Les résultats ainsi obtenus diffèrent en certains points deceux issus d'autres simulations qui introduisaient certaines hypothèses simplificatrices :développement temporel, emploi d'une chimie globale ou encore lois de transport simplifiées.En revanche, ils reproduisent certains tendances déjà observées dans un certainnombre d'études expérimentales conduites dans des conditions similaires.

Les plasmas inductifs radiofréquence basse pression sont utilisés pour la gravure de motifs nanométriques. Ces plasmas contiennent des ions négatifs et peuvent être soumis à deux types d'instabilités. La première est décrite par des oscillations de relaxation entre les modes de couplage de l'énergie capacitif (E) et inductif (H). Les mesures temporelles au cours de ces oscillations sont conformes aux résultats publiés. L'autre instabilité, liée au transport des espèces chargées, n'existe que lorsque le plasma peut diffuser. Des mesures spatio-temporelles montrent que les fluctuations des paramètres plasma sont liées à la formation et à la propagation périodique d'une double couche d'amplitude moyenne. Cette double couche est une gaine interne séparant un plasma haute densité, forte temp! érature électroniques et sans ions négatifs dans la source d'un plasma faible densité, faible température électronique et de forte fraction d'ions négatifs dans la chambre de diffusion. Cette double couche propagative apparaît lorsque la fraction d'ions négatifs dépasse un seuil. Sous ce seuil, la double couche est statique et disparaît pour des fractions d'ions négatifs très faibles. Les transitions entre ces trois états sont soumises à une forte activité acoustique ionique. L'ajout d'un champ magnétique modifie le chauffage des électrons par absorption de l'onde hélicon (mode W). Des oscillations de relaxation H/W ont été observées dans une fenêtre restreinte de paramètres. D'autre part, le confinement des espèces chargées est modifié: au centre, le plasma ne contient pas d'ions négatifs; au bord, le plasma ne contient plus d'électrons, un plasma ions - ions a été créé.

La couche-limite atmosphérique turbulente stable, particulièrement en zone de relief, n'est pas totalement comprise. Elle est, donc, mal représentée par les modèles atmosphériques. En présence de pente et d'un refroidissement du sol, l'augmentation locale de masse volumique génère un écoulement catabatique. En région de montagne, le maximum de vent est généralement enregistré à une hauteur (z_j) de 1-10 m. Le jet de paroi engendre un changement de signe du flux de qdm ainsi qu'une variabilité du flux de chaleur sensible proche du sol. Ces variabilités de flux turbulents contredisent l'applicabilité de la théorie des similitudes de Monin-Obukhov (TSMO), pourtant utilisée de manière universelle dans les modèles atmosphériques. Si la TSMO est discutable pour les cas très stables, c'est en présence de pente qu'elle n'est naturellement plus valide puisqu'elle néglige le couplage entre les équations de vitesse et de température. Il est donc impossible de représenter correctement un écoulement catabatique (z_j O(1m)) par un modèle utilisant la TSMO et avec une résolution verticale de l'ordre de la hauteur du maximum de vent. L'objectif du travail de thèse est d'apporter une contribution dans la compréhension et la modélisation de ces écoulements.Afin de compléter les observations peu nombreuses, une campagne de mesure a été réalisée sur une pente raide (20-40 deg) : la pente ouest du Grand Colon (chaîne de Belledonne, Alpes). Les analyses spectrales témoignent de la sensibilité de l'écoulement local aux perturbations externes, même faibles. Les caractéristiques turbulentes classiques sont observées à haute-fréquence alors que des comportements moins standards sont observés aux fréquences intermédiaires ou basses et expliqués par la présence de perturbations turbulentes d'énergie du même ordre de grandeur que l'injection locale. Les cospectres montrent un comportement propre aux écoulements catabatiques: recouvrement progressif selon z des corrélations croisées <0 et >0. La TSMO est mise en défaut pour l'écoulement observé et une solution alternative est utilisée pour estimer les flux en surface, permettant une bonne description de la vitesse de frottement.Le modèle 1D de surface de ISBA (Météo-France) est modifié pour répondre à la modélisation des écoulements catabatiques. Dans un premier temps, le modèle est validé sur un cas standard: en comparaison avec un modèle de Prandtl adapté. Dans un second temps, les données in-situ sont modélisées, d'abord en fournissant des profils de diffusivités effectives puis en utilisant un modèle modifié de turbulence d'ordre 1.5. Les modélisations 1D représentent correctement les champs moyens de vitesse et température mais montrent cependant des comportements trop diffusifs. Le modèle de longueur de mélange est principalement remis en cause, y compris en utilisant des paramétrisations adaptées.Des simulations LES 3D réalistes (Meso-NH, Météo-France) sont effectuées à haute résolution pour représenter le cas d'étude. Ces modélisations représentent finement les variabilités spatiales de l'écoulement catabatique. Cependant, des biais sont engendrés principalement par l'utilisation de la TSMO en condition aux limites de surface. Malgré la forte résolution spatiale, l'utilisation de la TSMO repousse à seulement z=2 m la perception des termes sources de l'écoulement catabatique par le modèle, alors que la source de l'écoulement atteint son maximum précisément en surface. Les modèles analytiques d'écoulement catabatique (de type Prandtl, qui pourraient aisément être intégrés en conditions aux limites) nécessitent de connaître "a priori" les profils de diffusivité. Ceci implique l'utilisation d'un modèle de turbulence. Le couplage du modèle 1D de surface (précédemment modifié et validé "off-line") est donc proposé pour répondre au manque de description de la physique par les CaL classiques de surface. Le travail préliminaire du couplage est présenté et des solutions sont proposées en perspective. / The stable atmospheric boundary layer, particularly in complex terrain, is not yet fully understood and it is thus still inadequately modelled. A surface cooling of a sloping terrain generates katabatic wind due to local density increase. This flow behaves as a wall-bounded turbulent jet, often simply modelled by a local balance between the buoyancy force and the turbulent friction. In mountainous regions, the wind maximum is typically observed at a height (z_j) of 1-10 m above the ground. The wall-bounded jet is responsible for a momentum-flux sign change and a heat-flux variability close to the ground. Those turbulent-flux variabilities are fully conflicting with the aplicability of the Monin-Obukhov similarity theory (MOST), which is nevertheless universally used in the atmospheric models to provide the surface boundary condition. If the MOST is already questionable for the very stable cases, it is obviously not valide over sloping surfaces because it neglectes the coupling of the wind and temperature equations, which constitutes the katabatic source. Hence, it is not possible to adequatly model a katabatic flow (z_j O(1m)) using the MOST, especially with a vertical resolution of the order of magnitude of z_j. The aim of the this PhD work is thus to improve the current understanding and modelling capacity of the katabatic winds.Since data sets of turbulent-katabatic-flow measurements are still scarce, a new field campain was carried out on a steep slope (20-40 deg): the west face of the Grand-Colon mountain (Belledonne ridge, French Alps). The experimental setup was mainly composed of a 6m mast with four sonic-anemometer levels (1, 2, 4 and 6m) to measure the turbulence on both sides of the katabatic jet. The spectral analysis shows the hight sensitivity of the local flow to external perturbations, even when these are weak. The hight-frequency subrange shows a classical behaviour (energy-injection frequency, inertial subrange), but the spectra of the intermediate and low-frequency subranges are less typical: turbulent perturbations with an energy of the order of magnitude of the local injection are present. A specific cospectra behaviour of the katabatic flows is shown: negative and positive cross-correlations overlap gradually, increasing z. The MOST fails in representing the observed flow and a surface-flux alternative estimation is succesfully used to describe the friction vellocity.The 1D surface model of ISBA (Météo-France) is modified to model katabatic flows. The model is firstly validated with a standard calibrated Prandtl model (with variable eddy difusivity). Secondly, the field data are modelled both with a prescribed effective diffusivity (from data) and using the 1.5-order turbulence scheme. The mean velocity and temperature fields are well reproduced, but it appears that the model is over-diffusive (which generates excessive fluxes), even when an adapted mixing-length is used.Realistic 3D LES simulations (Meso-NH, Météo-France) are computed with high resolution to model the field data. Spatial flow variabilities over sloping terrain are finely represented, but are biased, mainly due to the using of MOST for the surface boundary counditions. The using of MOST shifts the start of the katabatic source detection by the atmospheric model to a height of 2 m, while the katabatic source reaches its maximum at the surface. Analytical katabatic models (of the Prandtl type, which could be easily used to feed surface boundary counditions) need an "apriori" definition of the eddy and heat diffusivities. Currently, the general definition of these diffusivities is only possible by the use of turbulent models that include closures. The coupling of the previously-presented 1D surface model (validated off-line) is suggested to overcome the lack of physics description in the classic surface boundary counditions. Preliminary work on this coupling is developed and perspective solutions are proposed.