Les problèmes de valeurs propres pour les modèles d'installations pilotes rectangulaires sont formulés en séparant l'opérateur de transport spatial du temps. Pour des propriétés constantes et en l'absence de génération interne, la température ou la concentration transitoire satisfait une équation de diffusion telle que (\partial \theta/\partial t = D\nabla^2\theta). Dans un domaine rectangulaire, la séparation des variables réduit chaque direction spatiale à un problème de valeurs propres de la forme (\mathcal{L}_x\Phi=-\lambda^2\Phi), avec des fonctions de base sinus-cosinus choisies pour satisfaire les conditions aux limites de l'équipement.

La géométrie détermine l'opérateur spatial, mais les conditions aux limites déterminent les valeurs propres et les fonctions propres admissibles. Une fois ces fonctions choisies de manière à satisfaire les conditions aux limites homogènes, le champ transitoire est représenté comme une somme de modes spatiaux à décroissance exponentielle.

Commencer par le modèle de diffusion directeur

La conduction thermique et la diffusion de matière ont la même structure

Pour un solide rectangulaire, une plaque, une membrane ou une région remplie présentant des propriétés constantes, l'équation transitoire de conduction thermique est

[ \frac{\partial \theta}{\partial t}

D_H \left( \frac{\partial^2\theta}{\partial x^2} + \frac{\partial^2\theta}{\partial y^2} + \frac{\partial^2\theta}{\partial z^2} \right), ]

où (\theta) représente la température et

[ D_H=\frac{k}{\rho C_p} ]

est la diffusivité thermique.

La diffusion de matière prend une forme analogue,

[ \frac{\partial c}{\partial t}

D \left( \frac{\partial^2c}{\partial x^2} + \frac{\partial^2c}{\partial y^2} + \frac{\partial^2c}{\partial z^2} \right), ]

où (c) représente la concentration et (D) la diffusivité massique.

Le domaine rectangulaire fournit la structure des coordonnées

Pour une région rectangulaire de dimensions (L_x), (L_y) et (L_z),

[ 0<x<L_x,\qquad 0<y<L_y,\qquad 0<z<L_z. ]

Le modèle peut représenter une plaque de séchage, un lit catalytique, un bloc de stockage thermique, une section de membrane ou un canal rectangulaire d'installation pilote, à condition que les propriétés et la géométrie permettent l'approximation à coefficients constants.

Séparer le temps et l'espace

La solution séparée définit les fonctions de base

Un mode séparé peut s'écrire

[ \theta(x,y,z,t)

X(x)Y(y)Z(z)T(t). ]

La substitution dans l'équation de diffusion tridimensionnelle donne

[ \frac{1}{D_H}\frac{1}{T}\frac{dT}{dt}

\frac{1}{X}\frac{d^2X}{dx^2} + \frac{1}{Y}\frac{d^2Y}{dy^2} + \frac{1}{Z}\frac{d^2Z}{dz^2}. ]

Comme chaque terme dépend d'une variable indépendante différente, chaque terme spatial reçoit une constante de séparation.

Chaque coordonnée devient un problème de valeurs propres unidimensionnel

Pour la direction (x), l'équation spatiale est

[ \frac{d^2X}{dx^2}+\lambda_x^2X=0, ]

ou, en utilisant une notation opératorielle,

[ \mathcal{L}_xX=-\lambda_x^2X, \qquad \mathcal{L}_x=\frac{d^2}{dx^2}. ]

La solution générale est

[ X(x)=A_x\sin(\lambda_x x)+B_x\cos(\lambda_x x). ]

La même formulation s'applique aux directions (y) et (z), avec les valeurs propres (\lambda_y) et (\lambda_z).

Le facteur temporel décroît selon le nombre d'onde spatial total

Pour un mode possédant trois valeurs propres spatiales, l'équation temporelle est

[ \frac{dT}{dt} + D_H \left( \lambda_x^2+\lambda_y^2+\lambda_z^2 \right)T =0. ]

Par conséquent,

[ T(t)

\exp\left[ -D_H \left( \lambda_x^2+\lambda_y^2+\lambda_z^2 \right)t \right]. ]

Les modes spatiaux d'ordre élevé décroissent plus rapidement. Ainsi, la réponse à long terme est dominée par les plus petites valeurs propres admissibles.

Imposer les conditions physiques aux limites

Les limites de Dirichlet produisent des contraintes de type sinus

Une température ou une concentration imposée devient homogène après soustraction d'une référence ou d'une solution en régime permanent appropriée. La condition résultante est généralement

[ X(0)=0 \quad\text{or}\quad X(L_x)=0. ]

Par exemple, l'imposition de

[ X(0)=0 ]

élimine le terme en cosinus :

[ X(x)=A_x\sin(\lambda_x x). ]

Si l'autre limite satisfait également

[ X(L_x)=0, ]

alors

[ \sin(\lambda_x L_x)=0, ]

de sorte que les valeurs propres sont

[ \lambda_{x,n}=\frac{n\pi}{L_x}, \qquad n=1,2,3,\ldots ]

Les fonctions de base correspondantes sont des fonctions sinus.

Les limites de Neumann produisent des contraintes de type cosinus

Une paroi isolée ou une limite à flux nul est représentée par

[ \frac{dX}{dx}=0. ]

En (x=0),

[ X'(0)=0 ]

élimine le terme en sinus et donne

[ X(x)=B_x\cos(\lambda_x x). ]

Si l'autre paroi est également isolée,

[ X'(L_x)=0, ]

les valeurs propres satisfont

[ \sin(\lambda_xL_x)=0. ]

Cela donne

[ \lambda_{x,n}=\frac{n\pi}{L_x}, \qquad n=0,1,2,\ldots ]

Le mode (n=0) est le mode spatial constant. Il est physiquement important pour les systèmes isolés, car la température ou la concentration moyenne totale peut rester non nulle.

Les limites de Robin représentent une résistance au transfert de surface

Une limite de transfert thermique ou massique convectif s'écrit généralement

[ -k\frac{\partial \theta}{\partial x}

h(\theta-\theta_\infty), ]

ou, après nondimensionnalisation et homogénéisation,

[ X'(L_x)+A X(L_x)=0. ]

Ici, (A) contient le coefficient de transfert pertinent ainsi que l'échelle de conductivité ou de diffusivité. La même forme mathématique décrit le transfert thermique convectif et le transfert de matière interfacial.

L'application de cette condition à

[ X(x)=A_x\sin(\lambda_xx)+B_x\cos(\lambda_xx) ]

produit une équation caractéristique faisant intervenir (\lambda_x). Contrairement aux cas purs de Dirichlet ou de Neumann, les valeurs propres sont généralement les racines d'une équation transcendante et doivent être calculées numériquement.

Les conditions aux limites doivent être homogénéisées avant de choisir la base

Une température ou une concentration de surface imposée non nulle n'est pas représentée directement par une fonction de base qui s'annule à la limite. La procédure habituelle consiste à écrire

[ \theta=\theta_{\mathrm{ss}}+\theta', ]

où (\theta_{\mathrm{ss}}) satisfait les conditions aux limites non homogènes et (\theta') satisfait les conditions homogènes.

Les fonctions propres sont construites pour (\theta'), car cela permet aux termes de bord apparaissant lors de l'intégration par parties de s'annuler ou de se réduire à la forme de Robin appropriée.

Construire la base rectangulaire

Les fonctions propres multidimensionnelles sont des produits

Une fois les fonctions propres unidimensionnelles déterminées, une fonction propre rectangulaire tridimensionnelle est formée comme suit :

[ \Phi_{mnp}(x,y,z)

X_m(x)Y_n(y)Z_p(z). ]

Sa valeur propre totale est

[ \Lambda_{mnp}^2

\lambda_{x,m}^2 + \lambda_{y,n}^2 + \lambda_{z,p}^2. ]

La solution transitoire complète peut alors s'écrire

[ \theta'(x,y,z,t)

\sum_{m,n,p} C_{mnp} \Phi_{mnp}(x,y,z) \exp\left(-D_H\Lambda_{mnp}^2t\right). ]

Pour la diffusion de matière, remplacer (D_H) par la diffusivité massique (D).

Les modèles bidimensionnels utilisent la même construction

Si le composant de l'installation pilote est modélisé comme une plaque ou un canal rectangulaire présentant une variation négligeable dans une direction, la solution se réduit à

[ \theta'(x,y,t)

\sum_{m,n} C_{mn} X_m(x)Y_n(y) \exp\left[ -D_H \left( \lambda_{x,m}^2+\lambda_{y,n}^2 \right)t \right]. ]

Un modèle de plaque unidimensionnel est la version correspondante à une seule coordonnée.

Les coefficients proviennent des produits scalaires

Les coefficients sont déterminés en projetant le champ initial de température ou de concentration sur les fonctions propres. En une dimension,

[ C_n

\frac{\displaystyle\int_0^{L_x} f(x)X_n(x),dx} {\displaystyle\int_0^{L_x} X_n^2(x),dx}, ]

où (f(x)) représente la perturbation initiale après retrait de toute composante en régime permanent.

Pour un domaine rectangulaire tridimensionnel, la projection utilise le produit scalaire multidimensionnel

[ C_{mnp}

\frac{\iiint f(x,y,z) \Phi_{mnp}(x,y,z),dV} {\iiint \Phi_{mnp}^2(x,y,z),dV}. ]

Cette orthogonalité rend l'expansion spectrale exploitable sur le plan numérique.

Comprendre le rôle des valeurs propres

Les valeurs propres déterminent le temps de réponse

Chaque mode possède un temps de décroissance caractéristique approximativement égal à

[ \tau_{mnp}

\frac{1}{D_H\Lambda_{mnp}^2}. ]

La plus petite valeur propre non nulle contrôle la réponse dominante à long terme et influence donc fortement les estimations du temps de chauffage, de refroidissement, de séchage ou de séjour associé à la diffusion.

Les groupes adimensionnels relient le modèle au fonctionnement de l'équipement

Pour les limites convectives, les équations de valeurs propres dépendent de groupes adimensionnels tels que le nombre de Biot,

[ \mathrm{Bi}=\frac{hL}{k}. ]

Le nombre de Biot compare la résistance à la conduction interne à la résistance convective externe. Il modifie donc les valeurs propres et détermine si le composant rectangulaire se comporte approximativement comme un corps spatialement uniforme ou développe d'importants gradients internes.

Les fonctions de base peuvent soutenir les démarches analytiques et numériques

Le développement en fonctions propres fournit des solutions analytiques de référence pour les modèles de procédés et la vérification expérimentale. Il peut également être utilisé pour contrôler les prédictions obtenues par différences finies, éléments finis ou CFD pour des équipements pilotes rectangulaires soumis à des conditions aux limites contrôlées.

Comprendre les compromis

Les solutions analytiques exigent des hypothèses restrictives

La formulation rectangulaire séparée est plus simple lorsque la géométrie est rectangulaire, que les propriétés sont constantes et que le coefficient de diffusion ne varie pas avec la température ou la concentration. Une forte variation des propriétés, le rayonnement non linéaire, les limites mobiles, les changements de phase ou la génération interne peuvent invalider la formulation simple en sinus-cosinus.

Les conditions de Robin nécessitent généralement une recherche numérique des racines

Pour des limites mixtes, les valeurs propres ne sont généralement pas disponibles sous la forme de simples multiples de (\pi/L). L'équation caractéristique doit être résolue numériquement et il faut veiller à identifier toutes les racines pertinentes sans manquer des modes rapprochés ou répétés.

Une géométrie idéale peut ne pas représenter l'ensemble de l'installation pilote

Une installation pilote réelle peut comporter des angles, des supports, des chicanes, des résistances de contact, un écoulement non uniforme ou des coefficients de transfert variables dans l'espace. Un modèle à fonctions propres rectangulaires reste utile, mais ses prédictions doivent être comparées à des mesures ou à un modèle numérique de plus haute fidélité lorsque ces effets sont importants.

Des fonctions de base incorrectes produisent un comportement incorrect aux limites

Utiliser des fonctions sinus pour une paroi isolée ou des fonctions cosinus pour une paroi à valeur imposée revient à imposer une physique erronée. Dans les méthodes de transformation, cette incohérence laisse également les termes de bord non résolus lors de l'intégration par parties. Les équations transformées obtenues ne représentent donc pas correctement l'équipement pilote.

D'autres géométries nécessitent d'autres fonctions propres

La base sinus-cosinus est spécifique à la géométrie cartésienne rectangulaire. Les systèmes cylindriques nécessitent généralement des fonctions de Bessel, tandis que les systèmes sphériques utilisent des fonctions de Bessel sphériques et, lorsqu'une variation angulaire existe, des polynômes de Legendre.

Comment appliquer cela à votre projet

La formulation correcte dépend de la géométrie, du champ initial mesuré et des conditions réelles aux limites du transfert thermique ou massique.

Si votre objectif principal est une plaque ou un bloc rectangulaire : Formulez des problèmes de valeurs propres unidimensionnels indépendants dans chaque direction cartésienne et combinez leurs fonctions propres sinus-cosinus sous forme de produits.
Si votre objectif principal est une température ou une concentration de surface fixée : Homogénéisez la condition aux limites et utilisez des fonctions propres de type sinus s'annulant aux limites concernées.
Si votre objectif principal concerne des surfaces isolées : Utilisez des problèmes de valeurs propres de Neumann et incluez le mode constant de valeur propre nulle lorsque la physique le permet.
Si votre objectif principal est le transfert thermique ou massique convectif : Appliquez des conditions de Robin, dérivez l'équation caractéristique transcendante et calculez ses racines numériquement.
Si votre objectif principal est de prévoir le temps de réponse transitoire : Identifiez la plus petite valeur propre combinée non nulle, car son mode de décroissance exponentielle régit le comportement à long terme.
Si votre objectif principal est de valider une simulation d'installation pilote : Utilisez la série de fonctions propres comme référence, tout en vérifiant si la variation des propriétés, la génération interne ou les détails géométriques nécessitent un modèle numérique.

Un modèle de transport rectangulaire fiable est construit en faisant correspondre l'opérateur, les fonctions propres, les valeurs propres et les conditions aux limites à la physique réelle de l'installation pilote.

Tableau récapitulatif :

Aspect clé	Formulation ou enseignement
Modèle directeur	Équation de diffusion de la température ou de la concentration avec des propriétés constantes et sans génération interne
Domaine spatial	Coordonnées rectangulaires de dimensions (L_x), (L_y) et (L_z)
Séparation des variables	Exprimer le champ sous la forme (\theta=X(x)Y(y)Z(z)T(t))
Problème de valeurs propres	Chaque coordonnée satisfait (\mathcal{L}X=-\lambda^2X), généralement (X''+\lambda^2X=0)
Limite de Dirichlet	Une température ou une concentration fixée produit des fonctions propres de type sinus
Limite de Neumann	Les surfaces isolées ou à flux nul produisent des fonctions propres de type cosinus et peuvent inclure un mode constant
Limite de Robin	Le transfert thermique ou massique convectif produit des équations de valeurs propres transcendantes nécessitant des racines numériques
Base multidimensionnelle	Fonctions propres produits : (\Phi_{mnp}=X_mY_nZ_p)
Décroissance transitoire	Chaque mode décroît selon (\exp[-D\Lambda_{mnp}^2t])
Réponse dominante	La plus petite valeur propre non nulle contrôle le chauffage, le refroidissement, le séchage ou le comportement diffusif à long terme
Coefficients	Déterminés en projetant le champ initial sur les fonctions propres à l'aide de produits scalaires
Limites du modèle	La variation des propriétés, les changements de phase, la génération interne, la géométrie complexe et les transferts non linéaires peuvent nécessiter des méthodes numériques

Améliorer votre modèle de transport en installation pilote

Vous avez besoin d'un équipement éducatif ou industriel fiable pour étudier la conduction thermique transitoire, la diffusion de matière ou les processus de transport couplés ? Contactez LABPARK pour discuter de votre projet.

LABPARK fournit des installations pilotes éducatives et professionnelles d'opérations unitaires aux universités, aux instituts de recherche et aux entreprises, notamment des systèmes destinés au génie chimique, aux bioprocédés et aux biotechnologies, ainsi qu'au traitement de l'environnement et de l'eau. Nos installations pilotes vous aident à valider les modèles en fonctions propres, à étudier les effets des conditions aux limites, à mesurer les temps de réponse et à relier les équations théoriques de transport aux données pratiques d'exploitation.

Indiquez-nous votre application, vos conditions opératoires, votre géométrie ainsi que vos objectifs de formation ou de recherche. LABPARK peut vous aider à identifier une configuration d'installation pilote adaptée à l'enseignement, à la recherche en laboratoire, au développement de procédés et à la démonstration industrielle.