Sélectionnez les fonctions de base à partir du problème aux valeurs propres défini par la géométrie de l'équipement et ses conditions aux limites. Pour les domaines d'usine pilote rectangulaires, utilisez des fonctions propres sinus et cosinus, généralement sous forme de produits à travers les coordonnées cartésiennes. Pour les domaines cylindriques, utilisez des fonctions de Bessel dans la direction radiale, des fonctions sinus ou cosinus de Fourier dans la direction angulaire, et des fonctions trigonométriques appropriées le long de l'axe. Les fonctions correctes sont celles qui satisfont la forme homogène des conditions aux limites physiques, car elles forment alors une base orthogonale pour le calcul des coefficients de développement.

La géométrie détermine la famille de fonctions propres, tandis que les conditions aux limites déterminent les valeurs propres admissibles et les combinaisons spécifiques de sinus, cosinus ou de Bessel. Les fonctions de base qui correspondent aux deux produisent des modes indépendants et orthogonaux, rendant les modèles de diffusion ou de conduction beaucoup plus faciles à résoudre.

Commencer par le problème spatial aux valeurs propres

Séparer l'équation directrice par système de coordonnées

Les modèles de conduction et de diffusion transitoires se réduisent couramment, après séparation des variables, à un problème spatial aux valeurs propres tel que

[ \mathcal{L}_x\Phi=-\lambda^2\Phi, ]

où (\mathcal{L}_x) est la partie spatiale de l'opérateur de transfert de chaleur ou de masse.

Les fonctions de base appropriées sont les fonctions propres de cet opérateur sous les conditions aux limites de l'usine pilote. Elles ne sont pas sélectionnées uniquement parce qu'elles sont mathématiquement pratiques.

Faire correspondre la base aux conditions aux limites homogènes

Une condition aux limites peut prescrire la valeur de la température ou de la concentration, prescrire son flux, ou relier le flux à la valeur de surface.

Celles-ci correspondent globalement à :

Dirichlet : (\Phi=0), représentant une valeur fixe après que toute composante stationnaire a été supprimée.
Neumann : (\Phi'=0), représentant une limite isolée ou à flux nul.
Robin : (\Phi'+A\Phi=0), représentant un transfert convectif ou interfacial.

Pour les conditions non homogènes, le champ est souvent décomposé en une solution stationnaire ou de levée plus un reste transitoire. Les fonctions propres sont alors choisies pour satisfaire les conditions aux limites homogènes de ce reste.

Sélection des fonctions pour les géométries rectangulaires

Utiliser des fonctions propres sinus et cosinus dans les directions cartésiennes

Pour une région rectangulaire unidimensionnelle, l'équation spatiale a la forme

[ \frac{d^2\Phi}{dx^2}=-\lambda^2\Phi, ]

avec une solution générale

[ \Phi(x)=a\sin(\lambda x)+b\cos(\lambda x). ]

Les conditions aux limites déterminent quelle combinaison est autorisée et quelles valeurs discrètes de (\lambda) satisfont les deux limites.

Utiliser des produits pour les équipements rectangulaires multidimensionnels

Pour une dalle, un canal, une section de catalyseur ou un composant à plaque plane rectangulaire, les modes séparés ont généralement la forme

[ \Phi_{mn}(x,y)=X_m(x)Y_n(y), ]

ou, en trois dimensions,

[ \Phi_{mnp}(x,y,z)=X_m(x)Y_n(y)Z_p(z). ]

Chaque facteur est sélectionné indépendamment à partir des conditions aux limites sur sa paire de surfaces correspondante. Le mode complet est orthogonal car les facteurs unidimensionnels sont orthogonaux.

Appliquer la condition aux limites pour choisir les modes sinus ou cosinus

Une limite à valeur nulle en (x=0) élimine le terme cosinus, produisant des modes sinus. Une limite à flux nul en (x=0) élimine le terme sinus, produisant des modes cosinus.

Par exemple, sur un intervalle (0\leq x\leq L) :

(\Phi(0)=\Phi(L)=0) conduit à des modes sinus tels que (\sin(n\pi x/L)).
(\Phi'(0)=\Phi'(L)=0) conduit à des modes cosinus tels que (\cos(n\pi x/L)), incluant un mode à valeur propre nulle constante.
Une condition fixe à une extrémité et une condition à flux nul à l'autre produit une séquence de valeurs propres décalée, impliquant couramment des multiples demi-entiers de (\pi/L).
Les limites convectives produisent des valeurs propres de Robin déterminées par une équation caractéristique plutôt que par les multiples entiers les plus simples.

Les valeurs propres exactes dépendent des dimensions et des paramètres de limite de l'équipement.

Sélection des fonctions pour les géométries cylindriques

Utiliser les fonctions de Bessel pour la variation radiale

En coordonnées cylindriques, la partie radiale de l'opérateur de diffusion ou de conduction contient le facteur (1/r). L'équation aux valeurs propres radiale diffère donc de l'équation cartésienne et conduit aux fonctions de Bessel.

Pour les problèmes axisymétriques, la solution radiale est généralement basée sur

[ J_0(\lambda r), ]

où (J_0) est la fonction de Bessel de première espèce d'ordre zéro. Selon le domaine et les conditions aux limites, d'autres combinaisons ou ordres peuvent être nécessaires.

Un cylindre solide qui inclut l'axe exclut normalement la solution de Bessel singulière (Y_\nu(\lambda r)), car elle n'est pas finie en (r=0). Une région annulaire, en revanche, peut nécessiter une combinaison linéaire de (J_\nu) et (Y_\nu), puisque le domaine n'inclut pas l'axe.

Ajouter des fonctions angulaires de Fourier lorsque le processus n'est pas axisymétrique

Si la température ou la concentration varie avec la coordonnée angulaire, les fonctions propres angulaires sont

[ \cos(m\theta) \quad\text{et}\quad \sin(m\theta), ]

avec un numéro de mode azimutal entier (m).

Pour un modèle de réacteur tubulaire, de coque ou de colonne entièrement axisymétrique, seul le mode (m=0) est conservé, de sorte que la dépendance angulaire disparaît. Le chauffage non axisymétrique, les conditions de paroi ou la distribution de l'alimentation nécessitent les modes angulaires supérieurs.

Utiliser des fonctions trigonométriques le long de l'axe du cylindre

Pour un cylindre fini ou une unité pilote tubulaire, la coordonnée axiale produit généralement des fonctions propres sinus et cosinus, sélectionnées en fonction des conditions d'entrée, de sortie, de paroi d'extrémité ou de symétrie.

Un mode cylindrique séparé peut donc avoir la structure

[ \Phi_{mnp}(r,\theta,z)

R_{mn}(r), \Theta_m(\theta), Z_p(z), ]

où (R_{mn}) est une fonction radiale basée sur Bessel, (\Theta_m) est une fonction angulaire sinus ou cosinus, et (Z_p) est une fonction trigonométrique axiale.

Laisser les conditions aux limites déterminer les valeurs propres

Les conditions de Dirichlet sélectionnent les zéros de la fonction propre

Si une paroi cylindrique a une température ou une concentration prescrite, la fonction propre radiale doit s'annuler à cette paroi après homogénéisation.

Pour un cylindre solide avec (\Phi(R)=0), les valeurs propres radiales sont liées aux zéros de la fonction de Bessel pertinente. C'est la contrepartie cylindrique des valeurs propres sinus obtenues à partir d'un domaine rectangulaire avec deux limites à valeur fixe.

Les conditions de Neumann sélectionnent les zéros de la dérivée

Une paroi cylindrique isolée nécessite une dérivée radiale nulle :

[ \left.\frac{d\Phi}{dr}\right|_{r=R}=0. ]

Les valeurs propres admissibles sont alors déterminées par les zéros de la dérivée de la fonction de Bessel, plutôt que par les zéros de la fonction elle-même.

Le même principe s'applique aux géométries rectangulaires, où les conditions à flux nul sélectionnent des modes de type cosinus ou des équations caractéristiques basées sur la dérivée.

Les conditions de Robin créent des équations caractéristiques transcendantales

Une condition de paroi convective a la forme

[ -k\left.\frac{d\Phi}{dr}\right|_{r=R}

h\Phi(R), ]

ou une forme adimensionnelle équivalente. Dans un cylindre, la substitution de la solution radiale de Bessel produit une équation caractéristique impliquant à la fois une fonction de Bessel et sa dérivée.

Ces valeurs propres doivent généralement être calculées numériquement. C'est attendu et cela n'invalide pas le développement analytique.

Préserver correctement l'orthogonalité

Utiliser le produit scalaire approprié

Les fonctions propres rectangulaires sont généralement orthogonales selon le produit scalaire cartésien standard,

[ \langle f,g\rangle

\int f(x)g(x),dx. ]

Pour les domaines rectangulaires multidimensionnels, le produit scalaire s'étend sur tout le volume.

Les fonctions radiales cylindriques sont orthogonales avec le poids radial (r) :

[ \langle f,g\rangle_r

\int_0^R r,f(r)g(r),dr. ]

Ce poids est essentiel car il provient de l'élément de volume cylindrique (r,dr,d\theta,dz).

Normaliser les modes avant de calculer les coefficients

Après avoir identifié les fonctions propres, normalisez-les de sorte que

[ \langle \Phi_i,\Phi_j\rangle=\delta_{ij}, ]

où (\delta_{ij}) est le delta de Kronecker.

Le coefficient de développement pour une distribution initiale de température, de concentration ou de forçage (u) est alors obtenu à partir d'une projection scalaire,

[ c_i=\langle u,\Phi_i\rangle. ]

Pour les modèles cylindriques, cette projection doit inclure le poids radial correct et, le cas échéant, l'intégration sur l'angle et la position axiale.

Inclure le mode zéro lorsque la physique le permet

Les problèmes de Neumann purs peuvent contenir une fonction propre constante avec (\lambda=0). Ce mode représente la conservation de la moyenne spatiale lorsqu'il n'y a pas de flux net de chaleur ou de masse à travers les limites.

L'omettre peut produire une solution à long terme incorrecte, en particulier pour les solides isolés ou les systèmes de diffusion fermés.

Comprendre les compromis

Les bases rectangulaires sont plus simples à mettre en œuvre

Les bases sinus et cosinus ont des valeurs propres et des formules de normalisation simples pour les conditions aux limites rectangulaires standard. Elles sont souvent le choix le plus efficace pour les plaques planes, les dalles, les conduits rectangulaires et les domaines de calcul cartésiens.

Leur simplicité diminue lorsque différentes limites utilisent différents paramètres de Robin ou lorsque les propriétés des matériaux varient spatialement.

Les bases cylindriques sont physiquement précises mais plus complexes

Les fonctions de Bessel représentent correctement l'étalement radial et le changement de section transversale de l'équipement cylindrique. Les remplacer par des fonctions sinus cartésiennes peut fausser la solution radiale, en particulier près de l'axe ou lorsque les gradients radiaux dominent.

La complexité supplémentaire réside dans la nécessité de résoudre les équations caractéristiques de Bessel, d'appliquer une pondération radiale et de gérer la régularité en (r=0).

Les fonctions non correspondantes créent des erreurs de limite et de coefficient

Une base qui ne satisfait pas les conditions aux limites homogènes requises laisse des termes de limite indésirables après intégration par parties. L'équation de diffusion ou de conduction transformée contient alors des termes non résolus ou des contributions de forçage incorrectes.

Ce n'est pas seulement une perte d'élégance cosmétique. Cela peut changer le flux de surface prédit, les taux de décroissance transitoire et l'équilibre à long terme.

Les méthodes de fonctions propres idéales ont des limites de modélisation

Les constructions standard sinus, cosinus et de Bessel supposent un opérateur linéaire avec des propriétés appropriées, incluant souvent des propriétés de matériau constantes ou suffisamment régulières. Des propriétés fortement variables, une géométrie irrégulière, des coefficients de transfert non linéaires ou des sources internes compliquées peuvent nécessiter des méthodes numériques telles que les volumes finis, les éléments finis ou des techniques spectrales numériques.

La sélection de base basée sur la géométrie reste utile dans ces cas, mais elle peut servir de référence, de modèle réduit ou d'approximation localement appropriée plutôt que de solution exacte complète.

Comment appliquer cela à votre usine pilote

Commencez par homogénéiser les conditions aux limites et écrivez l'opérateur spatial dans le système de coordonnées qui correspond à l'équipement.

Si votre objectif principal est une dalle, un canal ou une unité à plaque plane rectangulaire : Utilisez des fonctions propres sinus et cosinus, en prenant des produits à travers les coordonnées et en sélectionnant chaque facteur pour satisfaire les limites de Dirichlet, Neumann ou Robin correspondantes.
Si votre objectif principal est une unité cylindrique solide ou tubulaire : Utilisez des fonctions de Bessel pour les modes radiaux, éliminez les solutions radiales singulières à l'axe et combinez-les avec des modes sinus ou cosinus axiaux.
Si votre objectif principal est une région annulaire ou côté coque : Autorisez à la fois les branches de fonction de Bessel régulières et singulières là où le domaine exclut l'axe, puis imposez les deux conditions de paroi radiale pour déterminer leur combinaison et leurs valeurs propres.
Si votre objectif principal est un modèle axisymétrique : Ne conservez que le mode angulaire zéro et utilisez le poids radial (r) lors de l'imposition de l'orthogonalité et du calcul des coefficients.
Si votre objectif principal est le transfert de paroi convectif : Formez l'équation caractéristique à partir de la condition de Robin et calculez ses valeurs propres numériquement plutôt que de forcer une séquence de Fourier standard.
Si votre objectif principal est la conservation dans un équipement isolé : Vérifiez explicitement la valeur propre zéro et incluez le mode constant lorsque les conditions aux limites le permettent.

Choisir une base qui correspond à la fois à la géométrie de l'usine pilote et à ses conditions aux limites transforme le modèle de transport en une représentation modale valide et efficace.

Tableau récapitulatif :

Géométrie de l'usine pilote	Fonctions de base recommandées	Considérations clés
Dalle, canal ou plaque plane rectangulaire	Fonctions sinus et cosinus, combinées sous forme de produits de coordonnées	Sélectionnez chaque facteur selon les limites de Dirichlet, Neumann ou Robin.
Cylindre solide ou unité tubulaire	Fonctions de Bessel radialement avec des fonctions sinus ou cosinus axiales	Excluez les solutions de Bessel singulières à l'axe et appliquez les valeurs propres radiales correctes.
Région annulaire ou côté coque	Combinaisons de fonctions de Bessel de première et deuxième espèce	Le domaine exclut l'axe, donc les deux branches radiales peuvent être nécessaires.
Équipement axisymétrique	Mode angulaire de Fourier zéro	Utilisez le poids radial `r` pour l'orthogonalité et les calculs de coefficients.
Parois convectives	Fonctions propres spécifiques à la géométrie avec valeurs propres de Robin calculées numériquement	Appliquez la relation flux-valeur pour former l'équation caractéristique.
Équipement isolé	Modes basés sur le cosinus ou la dérivée, incluant le mode constant le cas échéant	Conservez la valeur propre zéro pour préserver la conservation de la moyenne spatiale.

Vous planifiez un système de transport ou de séparation à l'échelle pilote ? Contactez LABPARK pour un soutien expert dans la sélection et l'application des bonnes usines pilotes d'opérations unitaires éducatives et professionnelles. LABPARK fournit des usines pilotes de génie chimique, de bioprocédés et biotechnologie, ainsi que de traitement de l'eau et de l'environnement pour les universités, les instituts de recherche et les entreprises. Nos solutions vous aident à modéliser précisément l'équipement réel, à valider le comportement de diffusion et de conduction, à former les opérateurs et les étudiants, et à accélérer la recherche sur les processus avec des conceptions adaptées à la géométrie, une instrumentation fiable et des conseils techniques pratiques.