Modélisation mathématique de la ventilation des plantes. Problèmes modernes de la science et de l'éducation Modèle mathématique de la ventilation d'alimentation et d'échappement

Le travail considère les processus de modélisation de la ventilation et de la dispersion de ses émissions dans l'atmosphère. La modélisation est basée sur la résolution du système d'équations de Navier-Stokes, les lois de conservation de la masse, de la quantité de mouvement, de la chaleur. Divers aspects de la résolution numérique de ces équations sont considérés. Un système d'équations est proposé pour calculer la valeur du coefficient de turbulence de fond. Pour l'approximation hypersonique, une solution est proposée, ainsi que les équations de la dynamique des fluides présentées dans l'article, de l'équation de position d'un gaz et d'une vapeur réels parfaits. Cette équation est une modification de l'équation de van der Waals et prend en compte plus précisément la taille des molécules de gaz ou de vapeur et leur interaction. Sur la base de la condition de stabilité thermodynamique, une relation est obtenue qui permet d'exclure les racines physiquement irréalisables lors de la résolution de l'équation par rapport au volume. L'analyse des modèles informatiques connus et des packages informatiques de la dynamique des fluides est effectuée.

la modélisation

ventilation

turbulence

équations de transfert de chaleur et de masse

équation d'état

du vrai gaz

dissipation

1. Berlyand M. Ye. Problèmes contemporains diffusion atmosphérique et pollution de l'air. - L. : Gidrometeoizdat, 1975.-- 448 p.

2. Belyaev NN Modélisation du processus de dispersion des gaz toxiques dans les conditions du bâtiment // Bulletin du DIIT. - 2009. - N° 26 - S. 83-85.

3. Byzova NL Études expérimentales de la diffusion atmosphérique et calculs de la diffusion des impuretés / NL Byzova, EK Garger, VN Ivanov. - L. : Gidrometeoizdat, 1985.-- 351 p.

4. Datsyuk TA Modélisation de la dispersion des émissions de ventilation. - SPb : SPBGASU, 2000 .-- 210 p.

5. Sauts AV Application d'algorithmes de graphiques cognitifs et de méthodes d'analyse mathématique pour étudier les propriétés thermodynamiques de l'isobutane R660A sur la raie de saturation : Subvention n° 2C/10 : rapport de recherche (conclusion) / GOUVPO SPBGASU ; mains. Gorokhov V.L., isp. : Sauts A.V.- SPb, 2011.- 30 p. : ill.- Bibliographie : p. 30.- N°GR 01201067977.-Inv. N° 02201158567.

introduction

Lors de la conception de complexes industriels et d'installations uniques, les problèmes liés à la garantie de la qualité de l'environnement de l'air et des paramètres de microclimat normalisés doivent être justifiés de manière exhaustive. Compte tenu du coût élevé de fabrication, d'installation et d'exploitation des systèmes de ventilation et de climatisation, des exigences accrues sont imposées à la qualité des calculs techniques. Choisir le rationnel solutions de conception dans le domaine de la ventilation, il faut pouvoir analyser la situation dans son ensemble, c'est-à-dire pour révéler la relation spatiale des processus dynamiques qui se déroulent à l'intérieur des locaux et dans l'atmosphère. Évaluer l'efficacité de la ventilation, qui dépend non seulement de la quantité d'air fournie à la pièce, mais également du schéma de distribution d'air et de la concentration adoptés produits dangereux dans l'air extérieur aux emplacements des prises d'air.

Objet de l'article- l'utilisation de dépendances analytiques, à l'aide desquelles sont effectués des calculs de quantité d'émissions nocives, pour déterminer les dimensions des canaux, des conduits d'air, des mines et le choix d'une méthode de traitement de l'air, etc. Dans ce cas, il est conseillé d'utiliser le produit logiciel Potok avec le module VSV. Pour préparer les données initiales, il est nécessaire de disposer de schémas des systèmes de ventilation conçus indiquant les longueurs des tronçons et les débits d'air aux extrémités des tronçons. Les données d'entrée pour le calcul sont la description des systèmes de ventilation et leurs exigences. Grâce à la modélisation mathématique, les problèmes suivants sont résolus :

sélection des meilleures options pour l'alimentation et l'évacuation de l'air ;
répartition des paramètres du microclimat par volume de locaux;
évaluation du régime aérodynamique du bâtiment;
sélection d'emplacements pour l'admission et l'évacuation de l'air.

Les champs de vitesse, de pression, de température, de concentrations dans le local et l'atmosphère se forment sous l'influence de nombreux facteurs dont la combinaison est difficile à prendre en compte dans les méthodes de calcul d'ingénierie sans l'utilisation d'un ordinateur.

L'utilisation de la modélisation mathématique dans les problèmes de ventilation et d'aérodynamique est basée sur la résolution des équations de Navier - Stokes.

Pour simuler des écoulements turbulents, il faut résoudre le système d'équations de conservation de la masse et de Reynolds (conservation de la quantité de mouvement) :

(2)

où t- temps, X= X je , j , k- coordonnées spatiales, vous=tu je , j , k - les composantes du vecteur vitesse, R- pression piézométrique, ρ - densité, τ je- composantes du tenseur des contraintes, s m- une source de masse, si je- composants de la source d'impulsions.

Le tenseur des contraintes s'exprime sous la forme :

(3)

où s ij- tenseur des vitesses de déformation ; δ je- tenseur des contraintes supplémentaires dues à la présence de turbulences.

Pour plus d'informations sur les champs de température T et concentration avec substances nocives, le système est complété par les équations suivantes :

équation de conservation de la chaleur

équation passive de conservation des impuretés avec

(5)

où CR- coefficient de capacité calorifique, - coefficient de conductivité thermique, k= k je , j , k est le coefficient de turbulence.

Coefficient de turbulence de base k bases est déterminée à l'aide d'un système d'équations :

(6)

où k F - coefficient de turbulence de fond, k f = 1-15 m 2 / s; = 0,1-04 ;

Les coefficients de turbulence sont déterminés à l'aide des équations :

(7)

Dans une zone ouverte à faible dissipation, la valeur k z est déterminé par l'équation :

k k = k 0 z /z 0 ; (8)

où k 0 - valeur k k en haut z 0 (k 0 = 0,1 m 2 / s à z 0 = 2 mètres).

En zone dégagée, le profil de vitesse du vent n'est pas déformé ;

Avec une stratification atmosphérique inconnue dans une zone ouverte, le profil de vitesse du vent peut être déterminé :

; (9)

où z 0 est une hauteur donnée (hauteur d'une girouette) ; vous 0 - vitesse du vent en altitude z 0 ; B = 0,15.

Sous réserve de la condition (10), le critère de Richardson local Ri défini comme:

(11)

Dérivons l'équation (9) des équations (7) et (8), à partir de là on exprime k socles

(12)

Assimilons l'équation (12) aux équations du système (6). On substitue (11) et (9) dans l'égalité obtenue, dans la forme finale on obtient le système d'équations :

(13)

Le terme palpitant, suivant les idées de Boussinesq, est représenté par :

(14)

où t- la viscosité turbulente et les termes supplémentaires dans les équations de transfert d'énergie et les composants d'impuretés sont modélisés comme suit :

(15)

(16)

Le système d'équations est fermé en utilisant l'un des modèles de turbulence décrits ci-dessous.

Pour les écoulements turbulents étudiés en pratique ventilatoire, il convient d'utiliser soit l'hypothèse de Boussinesq sur la petitesse des changements de densité, soit l'approximation dite « hypersonique ». Les contraintes de Reynolds sont supposées proportionnelles aux vitesses de déformation moyennes dans le temps. Le coefficient de viscosité turbulente est introduit, ce concept s'exprime ainsi :

. (17)

Le coefficient de viscosité effectif est calculé comme la somme des coefficients moléculaire et turbulent :

(18)

L'approximation "hypersonique" suppose la solution, avec les équations ci-dessus, de l'équation de stationnaire d'un gaz parfait :

ρ = p/(RT) (19)

où p - pression dans environnement; R- constante de gaz.

Pour des calculs plus précis, la densité des impuretés peut être déterminée à l'aide de l'équation de van der Waals modifiée pour les gaz et vapeurs réels

(20)

où les constantes N et M- prendre en compte l'association/dissociation des molécules de gaz ou de vapeur ; mais- prend en compte les autres interactions ; b" - prise en compte de la taille des molécules de gaz ; = 1 / .

En séparant de l'équation (12) la pression R et en le différenciant en volume (en tenant compte de la stabilité thermodynamique), on obtiendra la relation suivante :

. (21)

Cette approche permet de réduire significativement le temps de calcul par rapport au cas d'utilisation des équations complètes pour un gaz compressible sans diminuer la précision des résultats obtenus. Il n'y a pas de solution analytique aux équations ci-dessus. À cet égard, des méthodes numériques sont utilisées.

Pour résoudre les problèmes de ventilation associés au transfert de substances scalaires par un écoulement turbulent, lors de la résolution d'équations différentielles, un schéma de division des processus physiques est utilisé. Selon les principes de dédoublement, l'intégration aux différences finies des équations de l'hydrodynamique et le transfert convectif-diffus d'une substance scalaire à chaque pas de temps Δ t s'effectue en deux étapes. Dans un premier temps, les paramètres hydrodynamiques sont calculés. A la deuxième étape, les équations de diffusion sont résolues sur la base des champs hydrodynamiques calculés.

L'influence des transferts de chaleur sur la formation du champ de vitesse de l'air est prise en compte à l'aide de l'approximation de Boussinesq : un terme supplémentaire est introduit dans l'équation du mouvement pour la composante de vitesse verticale, qui prend en compte les forces de flottabilité.

Il existe quatre approches pour résoudre les problèmes de mouvement de fluide turbulent :

modélisation directe "DNS" (solution des équations de Navier-Stokes non stationnaires) ;
solution des équations de Reynolds moyennées "RANS", dont le système, cependant, n'est pas fermé et nécessite des relations de fermeture supplémentaires ;
méthode des grands tourbillons « LES » , qui est basé sur la résolution d'équations de Navier - Stokes non stationnaires avec paramétrisation des tourbillons à l'échelle du sous-réseau ;
méthode "DES" , qui est une combinaison de deux méthodes : dans la zone de flux séparés - "LES", et dans la zone de flux "lisse" - "RANS".

La plus intéressante du point de vue de la précision des résultats obtenus est sans aucun doute la méthode de simulation numérique directe. Cependant, à l'heure actuelle, les capacités de la technologie informatique ne permettent pas encore de résoudre des problèmes avec la géométrie et les nombres réels. Ré, et avec la résolution de tourbillons de toutes tailles. Par conséquent, lors de la résolution d'un large éventail de problèmes d'ingénierie, des solutions numériques des équations de Reynolds sont utilisées.

Actuellement, des packages certifiés tels que "STAR-CD", "FLUENT" ou "ANSYS / FLOTRAN" sont utilisés avec succès pour simuler des tâches de ventilation. Avec un problème correctement formulé et un algorithme de résolution rationnelle, la quantité d'informations obtenues permet au stade de la conception de choisir la meilleure option, mais effectuer des calculs à l'aide de ces programmes nécessite une formation appropriée, et leur utilisation incorrecte peut conduire à des résultats erronés.

Comme "cas de base", on peut considérer les résultats des méthodes de calcul d'équilibre généralement acceptées, qui permettent de comparer les valeurs intégrales caractéristiques du problème considéré.

Un des les points importants lorsqu'on utilise des logiciels universels pour résoudre les problèmes de ventilation, c'est le choix d'un modèle de turbulence. A présent, il est connu un grand nombre de divers modèles de turbulence qui sont utilisés pour fermer les équations de Reynolds. Les modèles de turbulence sont classés en fonction du nombre de paramètres pour les caractéristiques de la turbulence, respectivement, un paramètre, deux et trois paramètres.

La plupart des modèles semi-empiriques de turbulence, d'une manière ou d'une autre, utilisent "l'hypothèse de la localité du mécanisme de transfert turbulent", selon laquelle le mécanisme de transfert de quantité de mouvement turbulent est complètement déterminé en spécifiant les dérivées locales des vitesses moyennes et propriétés physiques liquides. Cette hypothèse ne prend pas en compte l'influence des processus se produisant loin du point considéré.

Les plus simples sont les modèles à un paramètre utilisant le concept de viscosité turbulente « n t», Et la turbulence est supposée isotrope. Une version modifiée du "n t-92" est recommandé pour la modélisation des jets et des flux séparés. Le modèle à un paramètre "S-A" (Spalart - Almaras), qui contient l'équation de transfert pour la quantité, donne également un bon accord avec les résultats expérimentaux.

Le manque de modèles avec une seule équation de transport est associé au fait qu'ils manquent d'informations sur la distribution de l'échelle de turbulence L... Par le montant L les processus de transfert, les méthodes de formation de la turbulence, la dissipation de l'énergie turbulente sont influencés. Dépendance universelle à déterminer L n'existe pas. Équation pour l'échelle de turbulence L s'avère souvent être exactement l'équation qui détermine la précision du modèle et, par conséquent, le domaine de son applicabilité. Fondamentalement, la portée de ces modèles est limitée à des écoulements de cisaillement relativement simples.

Dans les modèles à deux paramètres, sauf pour l'échelle de turbulence L, le taux de dissipation de l'énergie turbulente est utilisé comme deuxième paramètre . De tels modèles sont le plus souvent utilisés dans la pratique informatique moderne et contiennent les équations de transfert d'énergie de turbulence et de dissipation d'énergie.

Un modèle bien connu comprend des équations pour le transfert d'énergie de turbulence k et le taux de dissipation de l'énergie turbulente . Des modèles comme " k- e " peut être utilisé à la fois pour des écoulements proches des parois et pour des écoulements séparés plus complexes.

Des modèles à deux paramètres sont utilisés dans la version bas et haut Reynolds. Dans le premier, le mécanisme d'interaction du transport moléculaire et turbulent à proximité d'une surface solide est directement pris en compte. Dans la version à haut Reynolds, le mécanisme de transfert turbulent près d'une frontière solide est décrit par des fonctions spéciales près de la paroi qui relient les paramètres d'écoulement à la distance à la paroi.

À l'heure actuelle, les modèles les plus prometteurs sont les modèles SSG et Gibson-Launder, qui utilisent une relation non linéaire entre le tenseur des contraintes turbulentes de Reynolds et le tenseur des vitesses de déformation moyennes. Ils ont été conçus pour améliorer la prévision des courants de séparation. Étant donné que toutes les composantes tensorielles y sont calculées, elles nécessitent des ressources informatiques importantes par rapport aux modèles à deux paramètres.

Pour les écoulements séparés complexes, certains avantages ont été révélés par l'utilisation de modèles à un paramètre « n t-92 "," S-A " dans la précision de la prédiction des paramètres d'écoulement et dans le taux de comptage par rapport aux modèles à deux paramètres.

Par exemple, le programme "STAR-CD" prévoit l'utilisation de modèles tels que " k- e ”, Spalart - Almaras,“ SSG ”,“ Gibson-Launder ”, ainsi que la méthode des grands vortex“ LES ”, et la méthode“ DES ”. Les deux dernières méthodes sont mieux adaptées pour calculer le mouvement de l'air dans des conditions géométriques complexes, où de nombreuses régions de vortex séparées apparaîtront, mais elles nécessitent de grandes ressources de calcul.

Les résultats du calcul dépendent fortement du choix de la grille de calcul. Actuellement, des programmes de maillage spéciaux sont utilisés. Les cellules maillées peuvent être de différentes formes et tailles pour s'adapter au mieux à votre application spécifique. Le type de grille le plus simple est lorsque les cellules sont les mêmes et ont une forme cubique ou rectangulaire. Les programmes de calcul universels actuellement utilisés dans la pratique de l'ingénierie permettent de travailler sur des grilles non structurées arbitraires.

Pour effectuer des calculs pour la simulation numérique des problèmes de ventilation, il est nécessaire de définir les conditions aux limites et initiales, c'est-à-dire valeurs des variables dépendantes ou leurs gradients normaux aux limites du domaine de calcul.

Spécification avec un degré de précision suffisant des caractéristiques géométriques de l'objet à l'étude. À ces fins, il est possible de recommander pour la construction de modèles tridimensionnels des packages tels que "SolidWorks", "Pro / Engeneer", "NX Nastran". Lors de la construction d'une grille de calcul, le nombre de cellules est choisi de manière à obtenir une solution fiable avec un temps de calcul minimum. L'un des modèles de turbulence semi-empiriques les plus efficaces pour l'écoulement considéré doit être choisi.

DANS conclusion nous ajoutons qu'une bonne compréhension de l'aspect qualitatif des processus en cours est nécessaire pour formuler correctement les conditions aux limites du problème et évaluer la fiabilité des résultats. La modélisation des émissions de ventilation au stade de la conception des installations peut être considérée comme l'un des aspects de la modélisation de l'information visant à assurer la sécurité environnementale de l'installation.

Réviseurs :

Volikov Anatoly Nikolaevich, docteur en sciences techniques, professeur au département de l'approvisionnement en chaleur et en gaz et de la protection des bassins d'air, FGBOU VPOI "SPBGASU", Saint-Pétersbourg.
Polushkin Vitaly Ivanovich, docteur en sciences techniques, professeur, professeur du département de chauffage, ventilation et climatisation, FGBOU VPO "SPbGASU", Saint-Pétersbourg.

Référence bibliographique

Datsyuk T.A., Sauts A.V., Yurmanov B.N., Taurit V.R. MODÉLISATION DES PROCESSUS DE VENTILATION // Problèmes modernes de la science et de l'éducation. - 2012. - N°5.;
URL : http://science-education.ru/ru/article/view?id=6744 (date de consultation : 17/10/2019). Nous portons à votre connaissance les revues publiées par l'Académie des Sciences Naturelles

Envoyez votre bon travail dans la base de connaissances est simple. Utilisez le formulaire ci-dessous

Les étudiants, les étudiants diplômés, les jeunes scientifiques qui utilisent la base de connaissances dans leurs études et leur travail vous seront très reconnaissants.

Documents similaires

Les bases du fonctionnement du système contrôle automatique afflux ventilation d'échappement, sa construction et sa description mathématique. Équipement processus technologique... Sélection et calcul du régulateur. Etude de la stabilité de l'ATS, indicateurs de sa qualité.

dissertation, ajouté le 16/02/2011

caractéristiques générales et but, champ d'application pratique du système de contrôle automatique ventilation d'alimentation et d'extraction... Automatisation du processus de régulation, de ses principes et étapes de mise en œuvre. Le choix des fonds et leur justification économique.

thèse, ajoutée le 04/10/2011

Analyse des schémas types existants d'automatisation de la ventilation ateliers de fabrication. Modèle mathématique le processus de ventilation des locaux industriels, la sélection et la description des équipements d'automatisation et des contrôles. Calcul du coût d'un projet d'automatisation.

thèse, ajoutée le 11/06/2012

Analyse comparative caractéristiques techniques dessins typiques tours de refroidissement. Éléments des systèmes d'approvisionnement en eau et leur classification. Modèle mathématique du processus d'alimentation en eau en circulation, sélection et description des équipements d'automatisation et des éléments de contrôle.

thèse, ajoutée le 04/09/2013

Caractéristiques générales de l'oléoduc. Caractéristiques climatiques et géologiques du site. Disposition générale de la station de pompage. Stations de pompage de tronc et parc de réservoirs PS-3 "Almetyevsk". Calcul du système de ventilation d'alimentation et d'extraction de l'atelier de pompage.

thèse, ajoutée le 17/04/2013

Analyse de l'élaboration d'un projet de design pour une canne décorative. L'héraldique en tant que discipline spéciale traitant de l'étude des armoiries. Méthodes de fabrication d'outillage pour modèles en cire. Étapes de calcul de la ventilation d'alimentation et d'échappement pour le compartiment de fusion.

thèse, ajoutée le 26/01/2013

Description de l'installation en tant qu'objet d'automatisation, options pour améliorer le processus technologique. Calcul et sélection des éléments d'un ensemble de moyens techniques. Calcul du système de contrôle automatique. Développement de logiciels applicatifs.

thèse, ajoutée 24/11/2014

Prévision conditions thermiques dans les zones desservies est une tâche multifactorielle. On sait que le régime thermique est créé à l'aide de systèmes de chauffage, de ventilation et de climatisation. Cependant, lors de la conception des systèmes de chauffage, l'effet des courants d'air créés par d'autres systèmes n'est pas pris en compte. Ceci est en partie dû au fait que l'effet des flux d'air sur le régime thermique peut être insignifiant avec la mobilité de l'air standard dans les zones desservies.

Application des systèmes chauffage radiant nécessite de nouvelles approches. Cela inclut la nécessité de se conformer aux normes d'exposition humaine sur les lieux de travail et de prendre en compte la répartition de la chaleur rayonnante sur les surfaces intérieures des structures enveloppantes. En effet, avec le chauffage radiant, ces surfaces sont majoritairement chauffées, qui, à leur tour, dégagent de la chaleur dans la pièce par convection et rayonnement. C'est grâce à cela que la température de l'air interne requise est maintenue.

En règle générale, pour la plupart des types de locaux, ainsi que les systèmes de chauffage, des systèmes de ventilation sont nécessaires. Ainsi, lors de l'utilisation de systèmes de chauffage au gaz radiant, la pièce doit être équipée de systèmes de ventilation. L'échange d'air minimum dans les locaux avec dégagement de gaz et de vapeurs nocifs est stipulé par SP 60.13330.12. Le chauffage, la ventilation et la climatisation sont au moins ponctuels et à une hauteur de plus de 6 m - au moins 6 m 3 pour 1 m 2 de surface au sol. Par ailleurs, les performances des systèmes de ventilation sont également déterminées par la destination des locaux et sont calculées à partir des conditions d'assimilation des émissions de chaleur ou de gaz ou de compensation de l'aspiration locale. Bien entendu, il faut également contrôler le taux d'échange d'air pour l'état d'assimilation des produits de combustion. La compensation des volumes d'air extrait est effectuée par des systèmes aération d'alimentation... Dans ce cas, un rôle important dans la formation du régime thermique dans les zones desservies appartient aux jets d'alimentation et à la chaleur apportée par ceux-ci.

Méthode de recherche et résultats

Ainsi, il devient nécessaire de développer un modèle mathématique approximatif des processus complexes de transfert de chaleur et de masse se produisant dans une pièce avec chauffage et ventilation radiants. Le modèle mathématique est un système d'équations pour les bilans air-chaleur pour les volumes et surfaces caractéristiques d'une pièce.

La solution du système vous permet de déterminer les paramètres de l'air dans les zones desservies lorsque différentes options placement d'appareils de chauffage radiant, en tenant compte de l'influence des systèmes de ventilation.

Considérons la construction d'un modèle mathématique à partir de l'exemple d'une salle de production équipée d'un système de chauffage radiant et n'ayant pas d'autres sources de dégagement de chaleur. Les flux de chaleur provenant des radiateurs sont répartis comme suit. Les courants convectifs montent vers la zone supérieure sous le plafond et dégagent de la chaleur vers la surface intérieure. La composante radiante du flux de chaleur de l'émetteur est perçue par les surfaces intérieures des structures enveloppantes extérieures de la pièce. À leur tour, ces surfaces dégagent de la chaleur par convection vers l'air interne et par rayonnement vers d'autres surfaces internes. Une partie de la chaleur est transférée à l'air extérieur à travers les structures d'enceinte extérieures. Le schéma de transfert de chaleur calculé est illustré à la Fig. 1a.

Considérons la construction d'un modèle mathématique à partir de l'exemple d'une salle de production équipée d'un système de chauffage radiant et n'ayant pas d'autres sources de dégagement de chaleur. Les courants convectifs montent vers la zone supérieure sous le plafond et dégagent de la chaleur vers la surface intérieure. La composante radiante du flux de chaleur de l'émetteur est perçue par les surfaces intérieures des structures enveloppantes extérieures de la pièce

Ensuite, nous examinerons la construction du schéma de circulation du flux d'air (Fig. 1b). Prenons le schéma d'organisation de l'échange d'air complémentaire. L'air est fourni en quantité M pr en direction de la zone desservie et est retiré de la zone supérieure avec un débit M dans = M etc. Au niveau du haut de la zone desservie, le débit d'air dans le jet est M p.L'augmentation du débit d'air dans le jet d'alimentation se produit en raison de l'air en circulation, qui est déconnecté du jet.

Introduisons les frontières conditionnelles des écoulements - surfaces sur lesquelles les vitesses n'ont que des composantes qui leur sont normales. En figue. 1b, les limites d'écoulement sont représentées par une ligne pointillée. Ensuite, nous sélectionnerons les volumes calculés : la zone desservie (espace avec le séjour constant des personnes) ; volumes du jet d'alimentation et des écoulements convectifs proches de la paroi. La direction des flux convectifs à proximité des parois dépend du rapport des températures de la surface interne des structures d'enceinte externes et de l'air ambiant. En figue. 1b montre un diagramme avec un flux convectif descendant près de la paroi.

Ainsi, la température de l'air dans la zone desservie t wz est formé en raison du mélange d'air des jets d'alimentation, des flux convectifs près de la paroi et de l'apport de chaleur convective de surfaces internes sol et murs.

En tenant compte des schémas développés d'échange thermique et de circulation des flux d'air (Fig. 1), nous composerons les équations des bilans chaleur-air pour les volumes choisis :

Ici avec- capacité calorifique de l'air, J / (kg · ° С); Q de est la puissance du système de chauffage radiant au gaz, W; Q avec et Q* c - transfert de chaleur par convection sur les surfaces intérieures du mur à l'intérieur de la zone desservie et le mur au-dessus de la zone desservie, W ; t pp, t c et t wz est la température de l'air dans le jet d'alimentation à l'entrée de la zone de travail, dans le flux convectif proche de la paroi et dans la zone de travail, ° C ; Q TP est la perte de chaleur de la pièce, W, égale à la somme des pertes de chaleur à travers les structures d'enceinte externes :

Le débit d'air dans le flux d'alimentation à l'entrée de la zone desservie est calculé à l'aide des dépendances obtenues par M.I.Grimitlin.

Par exemple, pour les distributeurs d'air qui créent des jets compacts, le débit dans le jet est :

où m- coefficient d'amortissement de la vitesse ; F 0 - section transversale du tuyau d'admission du distributeur d'air, m 2; X- distance du répartiteur d'air au point d'entrée dans la zone desservie, m ; À n - coefficient de non-isothermie.

La consommation d'air dans le flux convectif proche de la paroi est déterminée par :

où tс - température de la surface intérieure des parois extérieures, ° C.

Équations bilan thermique pour les surfaces limites sont :

Ici Q c, Q*c, Q pl et Q Fri - transfert de chaleur par convection sur les surfaces intérieures du mur dans la zone desservie - le mur au-dessus de la zone desservie, le sol et le revêtement, respectivement ; Q tp.s, Q* tp.s, Q TP.pl, Q tp.pt - perte de chaleur à travers les structures correspondantes; W avec, W*c, W PL, W Fri - flux de chaleur rayonnante de l'émetteur entrant dans ces surfaces. Le transfert de chaleur par convection est déterminé par la relation connue :

où m J - coefficient déterminé en tenant compte de la position de la surface et de la direction du flux de chaleur; F J - superficie, m 2; Δ t J est la différence de température entre la surface et l'air ambiant, °C ; J- indice du type de surface.

Perte de chaleur QтJ peut être exprimé comme

où t n - température de l'air extérieur, ° C; t J - températures des surfaces internes des structures d'enceinte externes, ° C; R et R n - résistance thermique et transfert de chaleur de la clôture extérieure, m 2 · ° C / W.

Un modèle mathématique des processus de transfert de chaleur et de masse sous l'action combinée du chauffage radiant et de la ventilation a été obtenu. Les résultats de la solution permettent d'obtenir les principales caractéristiques du régime thermique dans la conception des systèmes de chauffage radiant pour les bâtiments à usages divers, équipés de systèmes de ventilation.

Flux de chaleur radiante des radiateurs des systèmes de chauffage radiant Wj sont calculés à travers les zones de rayonnement mutuel selon la méthode pour une orientation arbitraire des émetteurs et des surfaces environnantes :

où avec 0 - émissivité d'un corps absolument noir, W / (m 2 · K 4); ε IJ - émissivité réduite des surfaces impliquées dans le transfert de chaleur je et J; H IJ - zone mutuelle de rayonnement des surfaces je et J, m2; T JE - température moyenne surface émettrice, déterminée à partir du bilan thermique de l'émetteur, K ; T J est la température de la surface absorbant la chaleur, K.

En substituant des expressions pour les flux de chaleur et les débits d'air dans les jets, nous obtenons un système d'équations qui sont un modèle mathématique approximatif des processus de transfert de chaleur et de masse pendant le chauffage radiant. Des programmes informatiques standard peuvent être utilisés pour résoudre le système.

Un modèle mathématique des processus de transfert de chaleur et de masse sous l'action combinée du chauffage radiant et de la ventilation est obtenu. Les résultats de la solution permettent d'obtenir les principales caractéristiques du régime thermique dans la conception des systèmes de chauffage radiant pour les bâtiments à usages divers, équipés de systèmes de ventilation.

Daria Denisikhina, Maria Lukanina, Mikhaïl Samoletov

DANS monde moderne il n'est plus possible de se passer de la modélisation mathématique du flux d'air dans la conception des systèmes de ventilation.

Dans le monde moderne, il est déjà impossible de se passer de la modélisation mathématique des flux d'air lors de la conception des systèmes de ventilation. Les techniques d'ingénierie conventionnelles sont bien adaptées aux pièces typiques et aux solutions de distribution d'air standard. Lorsqu'un designer est confronté à des objets non standards, des méthodes de modélisation mathématique doivent lui venir en aide. L'article est consacré à l'étude de la répartition de l'air pendant la saison froide dans l'atelier de fabrication de tuyaux. Cet atelier fait partie d'un complexe industriel situé dans un climat fortement continental.

Au XIXe siècle, équations différentielles pour décrire l'écoulement des liquides et des gaz. Ils ont été formulés par le physicien français Louis Navier et le mathématicien britannique George Stokes. Les équations de Navier - Stokes sont parmi les plus importantes en hydrodynamique et sont utilisées dans modélisation mathématique nombreux phénomènes naturels et problèmes techniques.

Par dernières années accumulé une grande variété d'objets géométriquement et thermodynamiquement complexes dans la construction. L'utilisation de méthodes numériques de dynamique des fluides augmente considérablement les possibilités de conception de systèmes de ventilation, permettant de prédire avec un haut degré de précision les distributions de vitesse, pression, température, concentration des composants en tout point d'un bâtiment ou de l'un de ses locaux .

L'utilisation intensive des méthodes numériques de dynamique des fluides a commencé en 2000, avec l'apparition des coques logicielles universelles (packages CFD), qui ont permis de trouver des solutions numériques au système d'équations de Navier - Stokes par rapport à un objet d'intérêt. Depuis cette époque, "BURO TEKHNIKI" s'est engagé dans la modélisation mathématique en relation avec les problèmes de ventilation et de climatisation.

Description de la tâche

Dans cette étude, des simulations numériques ont été réalisées à l'aide de STAR-CCM+, un progiciel CFD développé par CD-Adapco. La performance de cet ensemble dans la résolution des problèmes de ventilation a été
Il a été testé à plusieurs reprises sur des objets de complexité variable, des bureaux aux salles de théâtre et aux stades.

Le problème est d'un grand intérêt tant du point de vue de la conception que de la modélisation mathématique.

Température de l'air extérieur -31°C. Il existe des objets à apport calorifique important dans la pièce : four de trempe, four de revenu, etc. Ainsi, il existe de grandes différences de température entre les structures extérieures d'enceinte et les objets intérieurs générateurs de chaleur. Par conséquent, la contribution du transfert thermique radiatif ne peut pas être négligée dans la simulation. Une difficulté supplémentaire dans la formulation mathématique du problème réside dans le fait qu'un train lourd avec une température de -31°C est amené dans le bâtiment plusieurs fois par quart de travail. Il se réchauffe progressivement, refroidissant l'air qui l'entoure.

Pour maintenir la température de l'air requise dans le volume de l'atelier (en saison froide au moins 15 ° C), le projet prévoit des systèmes de ventilation et de climatisation. Au stade de la conception, le débit et la température de l'air fourni nécessaires pour maintenir les paramètres requis ont été calculés. La question restait - comment alimenter en air le volume de l'atelier afin d'assurer la répartition la plus uniforme de la température dans tout le volume. La modélisation a permis pendant un temps relativement court (deux à trois semaines) de voir le modèle de flux d'air pour plusieurs options d'alimentation en air, puis de les comparer.

ÉTAPES DE LA MODÉLISATION MATHÉMATIQUE

Construire une géométrie solide.
Diviser l'espace de travail en cellules de la grille de calcul. Des dispositions doivent être prises à l'avance pour les zones dans lesquelles un raffinement supplémentaire des cellules est nécessaire. Lors de la construction d'une grille, il est très important de trouver un terrain d'entente où la taille des cellules est suffisamment petite pour obtenir des résultats corrects, tandis que le nombre total de cellules ne sera pas si important que de prolonger le temps de calcul à un moment inacceptable. Par conséquent, la construction de grilles est tout un art qui vient avec l'expérience.
Définir les conditions aux limites et initiales conformément à l'énoncé du problème. Une compréhension des spécificités des tâches de ventilation est requise. Un rôle important dans la préparation du calcul est joué par bon choix modèles de turbulence.
Sélection d'un modèle physique et d'un modèle de turbulence appropriés.

Résultats de la simulation

Pour résoudre le problème considéré dans cet article, toutes les étapes de la modélisation mathématique ont été franchies.

Pour comparer l'efficacité de la ventilation, trois options d'alimentation en air ont été sélectionnées : aux angles par rapport à la verticale 45°, 60° et 90°. L'air était fourni par des grilles de distribution d'air standard.

Champs de température et de vitesse calculés à différents angles d'avance air soufflé sont représentés sur la Fig. un.

Après analyse des résultats, l'angle de soufflage d'air égal à 90° a été choisi comme la plus réussie des options envisagées pour la ventilation de l'atelier. Avec cette méthode d'alimentation, des vitesses accrues ne sont pas créées dans la zone de travail et il est possible d'obtenir une image assez uniforme de la température et de la vitesse dans tout le volume de l'atelier.

Décision finale

Les champs de température et de vitesse dans trois sections transversales passant à travers les grilles de soufflage sont illustrés à la Fig. 2 et 3. La répartition de la température dans toute la pièce est uniforme. Ce n'est que dans la zone où les fours sont concentrés que les températures sont plus élevées sous le plafond. Il y a une zone plus froide dans le coin droit de la pièce la plus éloignée des fours. C'est ici qu'entrent les voitures froides de la rue.

Figure. 3 on voit clairement comment les jets horizontaux de l'air fourni se propagent. Avec ce mode d'alimentation, le jet d'alimentation a une portée suffisamment longue. Ainsi, à une distance de 30 m du réseau, la vitesse actuelle est de 0,5 m/s (à la sortie du réseau, la vitesse est de 5,5 m/s). Dans le reste de la pièce, la mobilité de l'air est faible, au niveau de 0,3 m/s.

L'air chauffé du four de trempe dévie le flux d'air d'alimentation vers le haut (Fig. 4 et 5). Le poêle chauffe beaucoup l'air qui l'entoure. La température au sol est plus élevée ici qu'au milieu de la pièce.

Le champ de température et les lignes de courant dans deux sections de l'atelier chaud sont illustrés à la Fig. 6.

conclusions

Les calculs ont permis d'analyser l'efficacité différentes façons alimentation en air de l'atelier de tuyauterie. Il a été constaté que lorsqu'il est alimenté par un flux horizontal, l'air d'alimentation se propage plus loin dans la pièce, contribuant à son chauffage plus uniforme. Cela ne crée pas de zones avec une trop grande mobilité de l'air dans la zone de travail, comme cela se produit lorsque l'air soufflé est fourni à un angle vers le bas.

L'utilisation de méthodes de modélisation mathématique dans les problèmes de ventilation et de climatisation est une direction très prometteuse, permettant au stade du projet de corriger la solution, d'éviter la nécessité de corriger les solutions de conception infructueuses après la mise en service des objets. ●

Daria Denisikhina - Chef du Département "Modélisation Mathématique" ;
Maria Lukanina - Ingénieur en chef du département « Modélisation mathématique » ;
Mikhaïl Samoletov - Directeur Exécutif de la SARL "MM-Technologies"

Glebov R.S., étudiant de troisième cycle Tumanov M.P., candidat en sciences techniques, professeur agrégé

Antyushin S.S., étudiant de troisième cycle (Moscou institut d'état Electronique et Mathématiques (Université Technique)

ASPECTS PRATIQUES DE L'IDENTIFICATION DU MODÈLE MATHÉMATIQUE

UNITÉ DE VENTILATION

Dans le cadre de l'émergence de nouvelles exigences pour les systèmes de ventilation, les méthodes expérimentales de réglage des boucles de régulation fermées ne peuvent pas résoudre pleinement les problèmes d'automatisation du processus technologique. Les méthodes de réglage expérimentales ont des critères d'optimisation inhérents (critères de qualité de contrôle), ce qui limite le champ de leur application. Synthèse paramétrique d'un système de contrôle qui prend en compte toutes les exigences Termes de référence, nécessite un modèle mathématique de l'objet. L'article propose une analyse des structures des modèles mathématiques unité de ventilation, la méthode d'identification de l'unité de ventilation est considérée, la possibilité d'utiliser les modèles obtenus pour une application pratique est évaluée.

Mots clés : identification, modèle mathématique, unité de ventilation, étude expérimentale modèle mathématique, critères de qualité du modèle mathématique.

ASPECTS PRATIQUES DE L'IDENTIFICATION DU MODÈLE MATHÉMATIQUE

D'INSTALLATION DE VENTILATION

En relation avec l'apparition de nouvelles exigences aux systèmes de ventilation, les méthodes expérimentales d'ajustement des contours fermés de gestion ne peuvent "pas résoudre complètement un problème d'automatisation du processus technologique. Les méthodes expérimentales d'ajustement ont les critères d'optimisation mis de gestion) qui limite le domaine de leur application. La synthèse paramétrique du système de contrôle, le projet technique tenant compte de toutes les exigences, exige un modèle mathématique de l'objet. Dans l'article à résulter l'analyse des structures des modèles mathématiques de l'installation de ventilation, la méthode de l'identification de l'installation de ventilation est considérée, la possibilité de l'application des modèles reçus pour l'application dans la pratique est estimée.

Mots clés : identification, modèle mathématique, installation de ventilation, recherche expérimentale du modèle mathématique, critères de qualité du modèle mathématique.

introduction

Le contrôle des systèmes de ventilation est l'une des tâches principales de l'automatisation systèmes d'ingénierie imeuble. Les exigences relatives aux systèmes de contrôle des unités de ventilation sont formulées sous forme de critères de qualité dans le domaine temporel.

Principaux critères de qualité :

1. Temps transitoire (tnn) - temps nécessaire à la centrale de traitement d'air pour atteindre le mode de fonctionnement.

2. Erreur de régime permanent (eust) - l'écart maximal admissible de la température de l'air soufflé par rapport à celle définie.

Critères de qualité indirects :

3. Dépassement (Ah) - dépassement de puissance lors du contrôle de la centrale de traitement d'air.

4. Degré d'oscillation (y) - usure excessive des équipements de ventilation.

5. Degré d'atténuation (y) - caractérise la qualité et la vitesse d'établissement du régime de température requis.

La tâche principale de l'automatisation du système de ventilation est la synthèse paramétrique du contrôleur. La synthèse paramétrique consiste à déterminer les coefficients du régulateur pour assurer les critères de qualité du système de ventilation.

Pour la synthèse du régulateur de l'unité de ventilation, on choisit des méthodes d'ingénierie adaptées à l'application pratique, qui ne nécessitent pas l'étude du modèle mathématique de l'objet : la méthode de Ncho18-21gier (W), la méthode Chien -Hrope8-Re8, wsk (SNK). À systèmes modernes automatisation de la ventilation, des exigences élevées sont imposées aux indicateurs de qualité, les conditions limites admissibles des indicateurs sont rétrécies, des problèmes de contrôle multicritères apparaissent. Les méthodes d'ingénierie de réglage du régulateur ne permettent pas de modifier les critères de qualité du contrôle qui y sont énoncés. Par exemple, en utilisant la méthode N2 pour régler le régulateur, le critère de qualité est un décrément d'amortissement égal à quatre, et en utilisant la méthode SAE, le critère de qualité est la vitesse de balayage maximale en l'absence de dépassement. L'utilisation de ces méthodes pour résoudre des problèmes de contrôle multicritères nécessite un ajustement manuel supplémentaire des coefficients. Le temps et la qualité du réglage des boucles de contrôle, dans ce cas, dépendent de l'expérience de l'ingénieur de service.

Application moyens modernes la modélisation mathématique pour la synthèse du système de contrôle de l'unité de ventilation améliore significativement la qualité des processus de contrôle, réduit le temps de mise en place du système, et permet également de synthétiser des moyens algorithmiques de détection et de prévention des accidents. Pour simuler le système de contrôle, il est nécessaire de créer un modèle mathématique adéquat de l'unité de ventilation (objet de contrôle).

L'utilisation pratique de modèles mathématiques sans en évaluer l'adéquation soulève un certain nombre de problèmes :

1. Les réglages du régulateur obtenus au cours de la modélisation mathématique ne garantissent pas la conformité des indicateurs de qualité dans la pratique.

2. L'application en pratique de régulateurs avec un modèle mathématique intégré (contrôle forcé, extrapolateur de Smith, etc.) peut entraîner une détérioration des indicateurs de qualité. Si la constante de temps ne correspond pas ou que le gain est trop faible, le temps nécessaire à la centrale de traitement d'air pour atteindre le mode de fonctionnement augmente, lorsque le gain est surestimé, une usure excessive des équipements de ventilation se produit, etc.

3. L'application pratique de contrôleurs adaptatifs avec une estimation selon le modèle de référence entraînera également une détérioration des indicateurs de qualité, similaire à l'exemple ci-dessus.

4. Les réglages du régulateur obtenus par les méthodes de contrôle optimal ne garantissent pas la conformité des indicateurs de qualité en pratique.

Le but de cette étude est de déterminer la structure du modèle mathématique de l'unité de ventilation (le long de la boucle de régulation régime de température) et l'évaluation de son adéquation aux processus physiques réels de chauffage de l'air dans les systèmes de ventilation.

L'expérience de la conception de systèmes de contrôle montre qu'il est impossible d'obtenir un modèle mathématique adéquat à un système réel uniquement sur la base d'études théoriques des processus physiques du système. Par conséquent, lors du processus de synthèse du modèle de l'unité de ventilation, parallèlement aux études théoriques, des expériences ont été menées pour déterminer et affiner le modèle mathématique du système - son identification.

Le processus technologique du système de ventilation, l'organisation de l'expérience

et identification structurelle

L'objet de contrôle du système de ventilation est le climatiseur central, dans lequel le flux d'air est traité et fourni aux locaux ventilés. La tâche du système de contrôle local de la ventilation est de maintenir automatiquement la température de l'air soufflé dans le conduit. La valeur actuelle de la température de l'air est évaluée par un capteur installé dans la gaine de soufflage ou dans le local habité. La température de l'air soufflé est contrôlée par un chauffe-eau électrique ou à eau. Lors de l'utilisation d'un chauffe-eau, l'actionneur est une vanne à trois voies, lors de l'utilisation d'un radiateur électrique - un régulateur de puissance à largeur d'impulsion ou à thyristor.

L'algorithme de contrôle standard pour la température de l'air soufflé est un système de contrôle automatique en boucle fermée (ACS), avec un contrôleur PID comme dispositif de contrôle. La structure du système automatisé de contrôle de la température de l'air soufflé par ventilation est illustrée (Fig. 1).

Riz. 1. Schéma fonctionnel du système de contrôle automatisé de l'unité de ventilation (canal de contrôle de la température de l'air soufflé). Wreg - PF du régulateur, Zhio - PF de l'organe exécutif, Wcal - PF du réchauffeur, Wvv - fonction de transfert du conduit. u1 - consigne de température, XI - température dans la gaine, XI - relevés du capteur, E1 - erreur de régulation, U1 - action de contrôle du régulateur, U2 - traitement du signal du régulateur par l'actionneur, U3 - chaleur transférée par le réchauffeur au canal.

La synthèse d'un modèle mathématique du système de ventilation suppose que la structure de chaque fonction de transfert comprise dans sa composition est connue. L'utilisation d'un modèle mathématique contenant les fonctions de transfert des éléments individuels du système est une tâche difficile et ne garantit pas en pratique la superposition des éléments individuels avec le système d'origine. Pour identifier le modèle mathématique, il convient de diviser la structure du système de contrôle de la ventilation en deux parties : a priori connue (contrôleur) et inconnue (objet). La fonction de transfert de l'objet ^ environ) comprend : la fonction de transfert de l'organe exécutif ^ uo), la fonction de transfert du réchauffeur ^ cal), la fonction de transfert du conduit d'air ^ vv), la fonction de transfert du capteur ^ Rendez-vous). La tâche d'identifier l'unité de ventilation lors du contrôle de la température du flux d'air se réduit à déterminer la relation fonctionnelle entre le signal de commande envoyé à l'actionneur du réchauffeur U1 et la température du flux d'air XI.

Pour déterminer la structure du modèle mathématique de l'unité de ventilation, il est nécessaire de réaliser une expérience d'identification. L'obtention des caractéristiques souhaitées est possible grâce à une expérimentation passive et active. La méthode d'expérimentation passive est basée sur l'enregistrement de paramètres de processus contrôlés dans le fonctionnement normal de l'objet sans y introduire de perturbations délibérées. Pendant la phase d'installation, le système de ventilation n'est pas en fonctionnement normal, la méthode d'expérimentation passive n'est donc pas adaptée à nos besoins. La méthode de l'expérimentation active repose sur l'utilisation de certaines perturbations artificielles introduites dans l'objet selon un programme préétabli.

Il existe trois méthodes fondamentales d'identification d'objet actif : la méthode des caractéristiques transitoires (la réaction de l'objet au "pas"), la méthode de perturbation de l'objet avec des signaux de forme périodique (la réaction de l'objet aux perturbations harmoniques avec fréquences différentes) et la méthode de réaction de l'objet à l'impulsion delta. En raison de la forte inertie des systèmes de ventilation (TOB est de quelques dizaines de secondes à plusieurs minutes), l'identification par des signaux de péri

Pour une lecture plus approfondie de l'article, vous devez acheter le texte intégral. Les articles sont envoyés au format PDF au courrier indiqué lors du paiement. Le délai de livraison est moins de 10 minutes... Coût d'un article - 150 roubles.

Ouvrages scientifiques similaires sur le thème "Problèmes généraux et complexes des sciences naturelles et exactes"

COMMANDE D'UNITÉ D'AIR ADAPTATIVE AVEC FLUX D'AIR D'APPROVISIONNEMENT DYNAMIQUE
R.S. GLEBOVM.P. TUMANOV - 2012
La problématique de la gestion et de la modélisation des situations d'urgence dans les mines de pétrole
M. Yu. Liskova et I. S. Naumov - 2013
SUR L'APPLICATION DE LA THÉORIE DU CONTRLE PARAMÉTRIQUE POUR LES MODÈLES CALCULABLES D'ÉQUILIBRE GÉNÉRAL
ADILOV ZHEKSENBEK MAKEEVITCH, ASHIMOV ABDYKAPPAR ASHIMOVITCH, ASHIMOV ASKAR ABDYKAPPAROVITCH, BOROVSKY NIKOLAY YURIEVICH, BOROVSKY YURI VYACHESLAVOVICH, SULTANOV BAKHYTVI TURCHLY - 2010
MODÉLISATION D'UNE TOITURE BIOCLIMATIQUE PAR VENTILATION NATURELLE
OUEDRAOGO A., OUEDRAOGO I., PALM K., ZEGHMATI B. - 2008