Spelling suggestions: "subject:"grille""
101 |
Couplage fluide-structure d'ordre (très) élevé pour des schémas volumes finis 2D Lagrange-projection / High-order fluid-structure coupling with conservative Lagrange-remap finite volume schemes on Cartesian gridsDakin, Gautier 09 November 2017 (has links)
Ce travail est consacré à l’étude numérique de l’interaction entre un fluide compressible et une structure indéformable, en adaptant une famille récente de schémas d’ordre très élevé à la prise en compte de conditions aux bords particulières entre le fluide et la structure. Plus précisément,on évalue l’apport de schémas d’ordre strictement supérieur à 3 par rapport à des stratégies plus classiques dans la littérature restreintes aux ordres 1 et 2. Un résultat important est qu’il est possible de réaliser le couplage à tout ordre et qu’il existe des configurations pour lesquelles on observe un gain important pour les ordres élevés. Une revue bibliographique est faite rappelant les résultats théoriques concernant les systèmes hyperboliques et décrivant les méthodes utilisées dans la littérature pour la simulation de la dynamique des gaz et la prise en compte des conditions aux bords. Un schéma sur grilles cartésiennes décalées et d’ordre très élevé est proposé pour la résolution des équations d’Euler en 1D/2D. Ce schéma est basé sur le formalisme Lagrange-projection et bien que formulé en énergie interne assure conservation et consistance faible grâce à un correctif en énergie interne. Parallèlement, l’étude pour les systèmes hyperboliques linéaires de discrétisation à l’ordre très élevé des conditions aux bords est faite. Elle met en évidence la nécessité pour l’ordre élevé de s’intéresser à la stabilité des schémas ainsi obtenus. À partir de ces travaux, la prise en compte de conditions aux bords en vitesse normale imposée est réalisée pour les équations d’Euler en 1D et 2D. Enfin, une procédure de couplage entre fluide compressible et structure indéformable est proposée. / This work is devoted to the construction of stable and high-order numerical methods in order to simulate fluid - rigid body interactions. In this manuscript, a bibliographic overview is done, which highlights theoretical results about hyperbolic system of conservation laws, as well as the methods available in the literature for the hydrodynamics simulation and the numericalboundary treatment. A high-order accurate scheme is proposed on staggered Cartesian grids to approximate the solution of Euler equations in 1D and 2D. The scheme relies on Lagrange-remap formalism, and although formulated in internal energy, ensures both conservation and weak consistency thanks to an internal energy corrector. In the same time, the study of high-order numerical boundary treatment for linear hyperbolic system is done. It highlights the necessity to focus especially on the linear stability of the effective scheme. Starting from the linear results, the numerical boundary treatment with imposed normal velocity is done for Euler equations in 1D and 2D. Last, the coupling between a compressible fluid and a rigid body is realized, using the designed procedure for numerical boudary treatment.
|
102 |
Contribution à la modélisation numérique des flammes turbulentes : comparaison DNS-EEM-Expériences / Contribution to numerical modelling of turbulent flames : DNS-EEM-comparisonsAlbin, Eric 27 April 2010 (has links)
La dynamique des flammes de prémélange est étudiée par deux approches numériques différentes. La première résout les équations compressibles de Navier-Stokes avec une chimie simplifiée (DNS). Afin de réduire les coûts de calcul, nous analysons et développons un schéma numérique à grille décalée. Le traitement des ondes acoustiques aux sorties est connu pour rendre les flammes cylindriques légèrement carrées. Ces déformations non-physiques sont expliquées en mettant en évidence la modélisation insuffisamment précise de l'accélération du fluide lorsque l'écoulement est oblique à la sortie. Une étude paramétrique et statistique de flammes turbulentes est menée en 2D et une simulation parallèle 3D est réalisée dans un domaine de (3cm)3. En considérant la flamme infiniment mince, l'approche EEM diminue considérablement les coûts de calcul. Les mêmes simulations sont réalisées et comparées aux résultats de DNS pour tester la capacité du modèle EEM à fournir des résultats quantitatifs. / We study premixed flame dynamics using DNS and EEM approaches. DNS solves compressible Navier-Stokes equations with simplified chemistry. To reduce computational costs and increase efifciency, we analyse and develop a modified staggered scheme. Treatment of acoustic waves at boundaries is known to slightly square cylindrical flames. We try to explain these unphysical distortions by highlighting the poor modeling of fluid acceleration when mainstrean is transverse to outow. A parametric and statistical study of expanding flames is carried out in 2D and also for an expanding (3cm)3 flame. The EEM approach models the flame as an infinitely thin interface. This perturbative strategy dramatically decreases cpu costs. Simulations are carried out and compared to DNS results to check the ability of EEM modeling to give quantitative results.
|
103 |
Influência do gradiente médio de velocidade na zona de contato na eficiência de clarificação de água para abastecimento por flotação / Influence of the velocity gradient in the contact zone in the efficiency of clarification of potable water by flotationSouza, Irene Monteiro da Franca 13 September 2012 (has links)
O presente trabalho teve como objetivo estudar a influência do gradiente médio de velocidade (Gzc), associado ao tempo de detenção hidráulico (tzc) e à taxa de aplicação superficial (TASzc), na zona de contato, na eficiência de clarificação de uma unidade retangular de flotação por ar dissolvido, em escala piloto, tratando água para abastecimento. Foram estudadas duas configurações da Unidade Piloto de Flotação por Ar Dissolvido (UPFAD) (A e B), as quais diferem no comprimento da zona de contato A (90 mm) e B (140 mm). A TASzc - A (133 ± 3 m/h) e B (86 ± 1 m/h) - e o tzc - A (54 ± 1s) e B (84 ± 1 s) - são diferentes. Foram avaliados três valores do Gzc A (1,8; 6,5;10,2) e B (0,9; 4,7; 6,3), através da inserção de telas com características diferentes no interior da zona de contato. E, Para cada valor de Gzc, foram investigados três valores de A/V (3,73; 4,90; 6,41). Em todos os ensaios realizados, foram fixados os seguintes parâmetros: vazão da água bruta preparada (5,4 m³/h); temperatura da água bruta preparada (26 ± 1ºC); pH de coagulação (6,5 ± 0,1); tempo de floculação (17,2 min); gradiente médio de velocidade de floculação (110 s-1); vazão de descarte após a floculação (0,8 m³/h); vazão de entrada para a flotação (Qe) (4,6 m³/h); velocidade de entrada para a flotação (ve) (180 m/h); velocidade cross-flow (45 m/h); TASzs (14,6 ± 0,2 m/h); pressão de saturação (4,5 ± 0,1 Kgf/cm²). A unidade piloto FAD foi alimentada com água preparada, que apresentou turbidez em torno de 7,4 ± 0,6 UNT, e cor aparente, em torno de 39,4 ± 4,3 UH. A fração removida dos parâmetros analisados para cada Gzc é menos sensível a variações de A/V na configuração A que na configuração B. Na configuração A, a tela que apresentou Gzc em torno de 6,5 s-1 gerou melhores resultados para cada A/V analisado. Na configuração B, a tela que gerou Gzc em torno de 4,7 s-1 proporcionou melhor fração removida de turbidez e cor para cada A/V estudado. Além disso, na configuração B, a fração removida de turbidez e cor foi menos sensível a variações de A/V, para Gzc de 4,7 s-1 que para os demais valores de Gzc. A tela que gerou resultados mais eficientes na configuração A foi a mesma para a configuração B. Portanto o uso da referida tela favoreceu o desempenho das configurações A e B. Comparando todas as situações estudadas nas duas configurações, o ensaio com a configuração A (com tzc de 54 s) com Gzc de 6,3 s-1 e A/V de 4,9 g/m³, apresentou as maiores frações de remoção de cor e turbidez. / The present work aimed to study the influence of velocity gradient (Gzc) associated with the hydraulic detention time (tzc) and the rate of surface application (TASzc) in the contact zone in the efficiency of clarification of a unit rectangular dissolved air flotation pilot scale, treating water supplies. Two different kinds of configurations of the reactor dissolved air flotation in pilot scale (UPFA) (A and B) were studied, which differ in the length of the contact zone - A (90 mm) and B (140 mm). The TASzc - A (133 ± 3 m/h) and B (86 ± 1 m/h) - and tzc - A (54 ± 1 s) and B (84 ± 1 s) were different. Three different values of Gzc were analysed - A (1.8, 6.5, 10.2 s-1) and B (0.9, 4.7, 6.3 s-1) by insertion of grilles with different characteristics within the zone contact. And for each value of Gzc, three values of A/V (3.73, 4.90, 6.41 g/m³) were investigated. In all trials, the following parameters were set: flow of synthetic raw water (5.4 m³/h), synthetic raw water temperature (26 ± 1ºC), coagulation pH (6.5 ± 0.1); flocculation time (17.2 minutes), the speed gradient flocculation (110 s-1); flow disposal after flocculation (0.8 m³/h), input flow to the flotation (Qe) (4, 6 m³/h), input speed for flotation (ve) (180 m/h), cross-flow velocity (45 m/h), TASzs (14.6 ± 0.2 m/h), saturation pressure (4.5 ± 0.1 kgf/cm ²). The pilot plant was fed with synthetic water, which showed turbidity of approximately 7.4 ± 0.6 NTU, and apparent color, about 39.4 ± 4.3 HU. The removal efficiency of the analyzed parameters for each Gzc is less sensitive to variations in the A/V in the configuration A than in the configuration B. In configuration A, the grille that showed Gzc about 6.5 s-1 gave better results for each A/V analysis. In configuration B, the grille that has generated Gzc around 4.7 s-1 provided the best removal efficiency of turbidity and color for each A/V studied. Furthermore, in the configuration B, the removal efficiency of color and turbidity was less sensitive to variations in A/V for Gzc 4.7 s-1 than for other values of Gzc. The grille that generated the best results in the configuration A was the same for the configuration B. Therefore the use of that grille favored the performance of configurations A and B. Comparing all the situations studied in two configurations, the test with the configuration A (with tzc the 54 s) with Gzc of 6.3 s-1 and A/V of 4.9 g/m³, had the highest removal efficiencies of color and turbidity.
|
104 |
Étude et exploitation de bolomètres de nouvelle génération à électrodes concentriques pour la recherche de matière noire froide non-baryonique dans l'expérience Edelweiss IIDomange, Jocelyn 30 September 2011 (has links) (PDF)
EDELWEISS est une expérience de détection directe de matière noire froide non-baryonique sous forme de particules massives et faiblement interagissantes (connues sous l'acronyme de WIMPs), qui constituent actuellement les candidats les plus populaires pour rendre compte de la masse manquante de l'Univers. Dans ce but, EDELWEISS utilise des bolomètres de germanium opérés à température cryogénique (20 mK environ) dans le Laboratoire Souterrain de Modane (LSM) à la frontière franco-italienne. En particulier, depuis 2008, un nouveau type de détecteur est en fonctionnement, équipé d'électrodes concentriques pour optimiser le rejet des évènements de surface (détecteurs à grilles coplanaires). Cette thèse se décompose en plusieurs axes de recherche. Tout d'abord, nous avons réalisé des mesures concernant la collecte des charges dans les cristaux. Les lois de vitesse des porteurs (électrons et trous) ont été déterminées dans le germanium à 20 mK dans la direction <100>, et une étude complète de la répartition des charges a été menée, avec une évaluation de l'anisotropie du transport et de la diffusion transverse des porteurs. Ces résultats permettent d'avoir une meilleure compréhension du fonctionnement interne des détecteurs d'Edelweiss. Ensuite, des études portant sur l'amélioration des performances ont été effectuées. Nous avons en particulier permis d'optimiser la procédure de régénération des cristaux et améliorer le rejet passif des évènements de surface (β). Le volume utile de détection des détecteurs a été évalué en utilisant les raies de deux radio-isotopes activés cosmiquement, le 68Ge et le 65Zn. Enfin, une étude exhaustive portant sur l'étude des spectres à basse énergie a été menée, ce qui permet de mettre au point une méthode d'analyse systématique pour la recherche de WIMPs de basse masse dans EDELWEISS.
|
105 |
Optimisation de la gestion des données pour les applications MapReduce sur des infrastructures distribuées à grande échelleMoise, Diana 16 December 2011 (has links) (PDF)
Les applications data-intensive sont largement utilisées au sein de domaines diverses dans le but d'extraire et de traiter des informations, de concevoir des systèmes complexes, d'effectuer des simulations de modèles réels, etc. Ces applications posent des défis complexes tant en termes de stockage que de calcul. Dans le contexte des applications data-intensive, nous nous concentrons sur le paradigme MapReduce et ses mises en oeuvre. Introduite par Google, l'abstraction MapReduce a révolutionné la communauté data-intensive et s'est rapidement étendue à diverses domaines de recherche et de production. Une implémentation domaine publique de l'abstraction mise en avant par Google a été fournie par Yahoo à travers du project Hadoop. Le framework Hadoop est considéré l'implémentation de référence de MapReduce et est actuellement largement utilisé à des fins diverses et sur plusieurs infrastructures. Nous proposons un système de fichiers distribué, optimisé pour des accès hautement concurrents, qui puisse servir comme couche de stockage pour des applications MapReduce. Nous avons conçu le BlobSeer File System (BSFS), basé sur BlobSeer, un service de stockage distribué, hautement efficace, facilitant le partage de données à grande échelle. Nous étudions également plusieurs aspects liés à la gestion des données intermédiaires dans des environnements MapReduce. Nous explorons les contraintes des données intermédiaires MapReduce à deux niveaux: dans le même job MapReduce et pendant l'exécution des pipelines d'applications MapReduce. Enfin, nous proposons des extensions de Hadoop, un environnement MapReduce populaire et open-source, comme par example le support de l'opération append. Ce travail inclut également l'évaluation et les résultats obtenus sur des infrastructures à grande échelle: grilles informatiques et clouds.
|
106 |
Schémas à deux-grilles pour la résolution du problème de Navier-Stokes instationnaire incompressibleAbboud, Hyam 03 July 2006 (has links) (PDF)
Dans ce travail nous nous intéressons à la résolution du problème d'évolution de Navier-Stokes incompressible totalement discrétisé en temps et en espace, en dimension deux par une méthode à deux grilles. Dans un premier temps, nous étendons la méthode à deux grilles, appliquée par V. Girault et J.-L. Lions au problème de Navier-Stokes instationnaire semi-discrétisé au problème de Navier-Stokes totalement discrétisé en temps (par un schéma d'ordre un) et en espace (par une méthode d'éléments finis d'ordre un). Dans la première étape, le problème non-linéaire est discrétisé en espace et en temps sur une grille grossière de pas d'espace H avec un pas de temps Delta t. Puis dans la deuxième étape, le problème, linéarisé autour de la vitesse u_H calculée à l'étape précédente, est discrétisé en espace sur une grille fine de pas d'espace h et le même pas de temps. L'idée de la méthode à deux grilles est que, sous des hypothèses adéquates, la contribution de u_H à l'erreur dans le terme non-linéaire en espace, est mesurée en norme L^2 en espace et en temps et a un ordre plus élevé que si elle était mesurée en norme H^1. Dans un deuxième temps, vu que le but est de gagner en ordre de convergence de l'erreur totale du schéma ainsi qu'en complexité nous étudions un schéma à deux grilles d'ordre deux en temps du problème totalement discrétisé en temps et en espace de Navier-Stokes. Nous présentons les résultats suivants: dans le cas de la résolution du schéma d'ordre un en temps, si h = H^2 = Delta t, alors l'erreur globale de l'algorithme à deux grilles est de l'ordre de h. Dans le cas du schéma d'ordre deux en temps, si h^2 = (Delta t)^2 = H^3, alors l'erreur globale de l'algorithme à deux grilles est de l'ordre de h^2: résultats identiques à ceux de la résolution directe du problème non-linéaire sur une grille fine.
|
107 |
Méthodes de décomposition de domaine espace-temps pour la formulation mixte de problèmes d'écoulement et de transport en milieu poreuxHoang, Thi Thao Phuong 11 December 2013 (has links) (PDF)
Cette thèse présente une contribution aux développements de méthodes numériques pour la simulation d'écoulements en milieu poreux, en particulier par des méthodes de décomposition de domaine espace--temps qui permettent l'utilisation de pas de temps différents dans les différents sous--domaines. Nous étudions deux types de méthodes: la première est basée sur une généralisation de l'opérateur de Steklov-Poincaré au cas de problèmes dépendants du temps, et la seconde est basée sur la méthode de Relaxation d'Onde Optimisée de Schwarz (OSWR) dans laquelle des conditions de transmission plus générales (Robin ou Ventcell) sont utilisées pour accélérer la convergence de l'algorithme. Ces deux méthodes sont étudiées sur une formulation mixte qui est bien adaptée à la modélisation de l'écoulement et du transport en milieu poreux. Nous considérons tout d'abord un problème de diffusion et formulons, pour chaque méthode, un problème sur l'interface espace -temps entre les sous-domaines. Le caractère bien posé de ces problèmes, avec des conditions aux limites de Dirichlet ou de Robin, est démontré. Les preuves de convergence de l'algorithme OSWR et de sa version semi-discrète sous forme mixte sont également données. Des expériences numériques sont menées en 2D pour comparer les performances des deux méthodes sur des problèmes fortement hétérogènes, et un préconditionneur de Neumann--Neumann dépendant du temps permet d'accélérer la première méthode. Les deux méthodes sont ensuite étendues au cas d'une équation d'advection-diffusion, l'advection et la diffusion étant traitées séparément grâce une technique de séparation d'opérateurs, ce qui permet d'utiliser des pas de temps différents pour les deux phénomènes dans chaque sous-domaine. Des conditions de transmission sont proposées séparément pour l'advection et pour la diffusion. La convergence des méthodes est étudiée sur des exemples numériques, pour des problèmes en régime d'advection dominante ou de diffusion dominante, et leur précision en temps est étudiée dans le cas de grilles non-conformes en temps. Deux exemples inspirés de la simulation du stockage de déchets nucléaires sont étudiés, et la simulation sur des temps longs est réalisée par l'intermédiaire de fenêtres en temps. Nous considérons également la méthode OSWR avec des conditions de transmission de Ventcell, étendues à la formulation mixte. Nous démontrons que les problèmes de sous--domaine avec des conditions aux limites de Ventcell sont bien posés. Nous comparons les performances des paramètres optimisés pour Ventcell et Robin dans le cas de problèmes hétérogènes pour une décomposition en deux sous-domaines. Enfin, nous étudions l'extension des deux méthodes au cas où l'interface représente une fracture pour un modèle réduit d'écoulement dans un milieu poreux fracturé.
|
108 |
Exécution d'applications parallèles en environnements hétérogènes et volatils : déploiement et virtualisationMiquée, Sébastien 25 January 2012 (has links) (PDF)
La technologie actuelle permet aux scientifiques de divers domaines d'obtenir des données de plus en plus précises et volumineuses, Afin de résoudre ces problèmes associés à l'obtention de ces données, les architectures de calcul évoluent, en fournissant toujours plus de ressources, notamment grâce à des machines plus puissantes et à leur mutualisation. Dans cette thèse, nous proposons d'étudier dans un premier temps le placement des tâches d'applications itératives asynchrones dans des environnements hétérogènes et volatils. Notre solution nous permet également de s'affranchir de l(hétérogénéité des machines hôtes tout en offrent une implantation facilitée de politiques de tolérance aux pannes, les expérimentations que nous avons menées sont encourageantes et montrent qu'il existe un réel potentiel quand à l'utilisation d'une telle plateforme pour l'exécution d'applications scientifiques.
|
109 |
Etude et modélisation compacte du transistor FinFET ultime / Study and compact modeling of ultimate FinFET transistorChevillon, Nicolas 13 July 2012 (has links)
Une des principales solutions technologiques liées à la réduction d’échelle de la technologie CMOS est aujourd’hui clairement orientée vers les transistors MOSFET faiblement dopés à multiples grilles. Ceux-ci proposent une meilleure immunité contre les effets canaux courts comparés aux transistors MOSFET bulk planaires (cf. ITRS 2011). Parmi les MOSFETs à multiples grilles, le transistor FinFET SOI est un candidat intéressant de par la similarité de son processus de fabrication avec la technologie des transistors planaires. En parallèle, il existe une réelle attente de la part des concepteurs et des fonderies à disposer de modèles compacts efficaces numériquement, précis et proches de la physique, insérés dans les « design tools » permettant alors d’étudier et d’élaborer des circuits ambitieux en technologie FinFET. Cette thèse porte sur l’élaboration d’un modèle compact orienté conception du transistor FinFET valide aux dimensions nanométriques. Ce modèle prend en compte les effets canaux courts, la modulation de longueur de canal, la dégradation de la mobilité, leseffets de mécanique quantique et les transcapacités. Une validation de ce modèle est réalisée par des comparaisons avec des simulations TCAD 3D. Le modèle compact est implémenté en langage Verilog-A afin de simuler des circuits innovants à base de transistors FinFET. Une modélisation niveau-porte est développée pour la simulation de circuits numériques complexes. Cette thèse présente également un modèle compact générique de transistors MOSFET SOI canaux long faiblement dopés à multiple grilles. La dépendance à la température est prise en compte. Selon un concept de transformation géométrique, notre modèle compact du transistor MOSFET double grille planaire est étendu pour s’appliquer à tout autre type de transistor MOSFET à multiple grille (MuGFET). Une validation expérimentale du modèle MuGFET sur un transistor triple grille est proposée. Cette thèse apporte enfin des solutions pour la modélisation des transistors MOSFET double grille sans jonction. / One of the main technological solutions related to downscaling of CMOS technology is now clearly oriented to lightly doped multigate MOSFETs. They offer better immunity against short channel effects compared to planar bulk MOSFETs (see ITRS 2011). Among the multigate MOSFETs, the SOI FinFET transistor is an interesting candidate because of the similarity of its manufacturing process with the planar transistor technology. In parallel, there is a real expectation on the part of designers and foundries to have compact models numerically efficient, accurate and close to the physics, and then inserted in to the design tools in order to study and develop ambitious circuits in FinFET technology. This thesis focuses on the development of a design-oriented compact model of FinFET transistor valid to nanoscale dimensions. This model takes into account the short channel effects, the channel length modulation, the mobility degradation, the quantum mechanic effects and the transcapacitances. A validation of this model is carried out by comparisons with 3DTCAD simulations. The compact model is implemented in Verilog-A to simulate innovative FinFET-based circuits. A gate-level modeling is developed for the simulation of complex digital circuits. This thesis also presents a generic compact modeling of multigate SOI MOSFETs with lightly doped channels and temperature dependent. According to a concept of geometric transformation, our compact model of the planar double-gate MOSFET is extended to be applied to any other type of multigate MOSFETs (MuGFET). An experimental validation of the MuGFET compact model with a triple gate transistor is proposed. This thesis finally brings solutions for the modeling of junction less double-gate MOSFET.
|
110 |
Discrétisation et commande frontière de systèmes vibro-acoustiques, une approche hamiltonienne à ports / Discretization and boundary control of vibroacoustic systems, a port-Hamiltonian approachTrenchant, Vincent 27 November 2017 (has links)
Cette thèse répond à une problématique de commande frontière d’une conduite acoustique dont l’actionnement est assuré par un réseau d’actionneurs/capteurs co-localisés constituant une peau active. Pour faire face au caractère intrinsèquement multiphysique de ce problème vibro-acoustique, nous avons choisi dans cette thèse d’employer une approche hamiltonienne à ports, approche structurée basée sur la représentation des échanges entre différents domaines énergétiques au sein d’un système et entre différents systèmes. Nous avons proposé une modélisation hamiltonienne à ports de l’équation d’onde interconnectée à la frontière au système d’actionnement distribué, correspondant à une formulation 2D du problème physique. Nous avons développé une méthode de discrétisation spatiale basée sur l’utilisation de différences finies sur plusieurs grilles en quinconce qui préserve la structure hamiltonienne à ports de l’équation d’onde. Cette méthode permet en outre d’interconnecter facilement le système discrétisé avec d’autres sous-systèmes, dans le but de mettre en place un actionnement par exemple. Son principal avantage sur d’autres méthodes préservatives de structure réside dans sa simplicité de mise en œuvre qui découle de l’utilisation de différences finies. Concernant la commande du système vibro-acoustique, nous avons proposé une méthode de synthèse de loi de commande distribuée pour les systèmes régis par deux lois de conservation en 1D. L’originalité de cette méthode réside en le fait qu’elle repose sur le calcul d’invariants structuraux (fonctions de Casimir) exploités afin de modifier la structure du système en boucle fermée. Les conditions d’application sur un système 2D sont étudiées et des résultats numériques valident les lois de commande synthétisées. / This thesis deals with the boundary control of an acoustic by a network of co-localised sensors/actuators which constitutes a smart skin. In order to cope with this multiphysical problem, we chose to place our study in the framework of port-Hamiltonian systems, a structured approach based on the representation of energy exchanges between different energy domains between different systems of subsystems. We proposed a port-Hamiltonian model of the wave equation interconnected through its boundary to the distributed actuation system, which corresponds to a 2D formulation of the physical problem. We developed a spatial discretization method based on the use of finite differences on several staggered grids that preserve the port-Hamiltonian structure of the wave equation. This method also permits to easily interconnect the discretized system with other subsystems, which is convenient for instance for control purposes. Its main advantage over other structure preserving methods is its simplicity of implementation which stems from the use of finite differences. In order to control the vibro-acoustic system, we proposed a control law synthesis method for systems governed by two conservation laws in 1D. The originality of this method lies in the fact that it relies on the computation of structural invariants (Casimir functions) exploited in order to modify the structure of the system in closed loop. The conditions of application of these laws on a 2D system are studied and numerical results validate the synthesized control laws.
|
Page generated in 0.0749 seconds