Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD...
Transcript of Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD...
![Page 1: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/1.jpg)
Modélisation OD/1Davec
COMSOL Multiphysics 3.4
COMSOL FRANCEStephan Savarese
[email protected] 76 46 49 01
17 avril 2008
LMCS 2008
Printemps de la Recherche
EDF – Chatou
![Page 2: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/2.jpg)
1. Introduction – E.D.P
Ft
ud
t
ue aa =Γ⋅∇+
∂∂+
∂∂
2
2
( ) sourceheatQTkt
TCp =∇⋅∇−
∂∂ρ
( ) ( )cRcDt
c ∇⋅−=∇⋅∇−∂∂
u
( ) 00 =∇⋅∇− Vrεε
( ) 0=∇⋅∇− Vσ
Heat Flux:
Mass Convection-Diffusion:
DC Conductive:
Electrostatics:
Des équations...
![Page 3: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/3.jpg)
Ft
ud
t
ue aa =Γ⋅∇+
∂∂+
∂∂
2
2
( )( )( ) uuFuuIu T ∇⋅−=∇+∇−⋅∇+ ρηρ p
dt
d
0=×∇+∂∂
EBt
( ) SourceQp
t
p
c=
∇−⋅∇++∂∂
ρρ0
12
2
2
( ) 0=⋅∇ uρ
Acoustics:
Gauss’s Law:
Faraday’s Law :
Navier-Stokes:
ρ=⋅∇ D
[ ] Fσ =⋅∇Structural Analysis:
Modéliserpresque tout!
... pour modéliser...
![Page 4: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/4.jpg)
Générateur électrique Ecoulement diphasique
Ingénierie biomédicale Contacteur test
... la réalité (multi)physique!
![Page 5: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/5.jpg)
Analyse financière:
En 1D
Calcul de la valeur d’un indice européen avec l’équationde Black-Scholes
![Page 6: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/6.jpg)
Solitons : équation de Korteweg-De Vries
En 1D
Deux solitons à l’instant t=1,875s
![Page 7: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/7.jpg)
Vagues en eau peu profonde
En 1D
• Equations de Saint-Venant :
z = profondeur d’eauv = vitesse verticaleg = accélération de la pesanteur
ν = viscosité cinématique
0)()()(
0)(
=∆−∇+⋅⋅∇+∂
∂
=⋅∇+∂∂
zvzgzvzvt
zv
zvt
z
υCoupe de profil du fond et de la surface d’un lac.
![Page 8: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/8.jpg)
Equation de Boltzmann
En 1D
Résolution de l’équation de Boltzmann à deux termes: propriétés de transport (haut) etfonction de distribution de l’énergie (milieu)des électrons dans un plasma d’oxygène(bas).
![Page 9: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/9.jpg)
Perte de charge dans un réseau hydraulique
En 1D
Calcul de débit dans un réseau de tuyaux.
Forme faible:
![Page 10: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/10.jpg)
Réduction de modèles complexes
Couplage 0D / 2Daxi
Inducteur dans un Circuit Amplificateur- Densité de flux magnétique
Vin 1 0 1Vcc 4 0 15Rg 1 2 100Cin 2 3 10uR1 4 3 47kR2 3 0 10kRE 7 0 1kCout 5 6 10uRl 6 0 1kQ1 5 3 7 Q2N2222.LIB "femlab.lib" .END
inducteur
Tension des signaux d’entrée, de sortie et de l’inducteur
Densité de flux magnétique
![Page 11: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/11.jpg)
2. La Famille COMSOL au complet !
COMSOL Script
AcousticsModule
Chemical Engineering
Module
Structural MechModule
MEMS ModuleRF Module
Reaction Engineering Lab
Signals andSystems Lab
OptimizationLab
CAD ImportModule
Material Library
AC/DCModule
Heat TransferModule
Earth Science Module
COMSOL Multiphysics Base Package :• Acoustics• Electromagnetics• Mass Diffusion• Fluid Mechanics• Heat Transfer• Structural Mechanics• Classical Equations• Moving Mesh / ALE • All Solvers• Model Library
Une famille qui s’agrandit...
![Page 12: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/12.jpg)
COMSOL Multiphysics Base Package :• Acoustics• Electromagnetics• Mass Diffusion• Fluid Mechanics• Heat Transfer• Structural Mechanics• Classical Equations• Moving Mesh / ALE • All Solvers• Model Library
COMSOL & Autres ProgrammesMATLAB
CAD ImportModule
COMSOL Script
SolidWorks
3D CAD
2D CADdxf
iges, sat, etc.
Other FEAPrograms
ExperimentalData
Postprocess to grid.Interpolation Function.
SPICEModels
Reaction Engineering Lab
ChemKinReactions
direct link
Word, PowerPointtext file export
M-file or struct export
Un environnement ouvert
![Page 13: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/13.jpg)
Signals & Systems Lab
• New functionality for model calibration and parameter estimationusing nonlinear least-squares methods
– Implemented for multiple sets of data
• New examples
![Page 14: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/14.jpg)
Signals and Systems Lab
• Contains over 100 functions to support :– signal processing, – system simulation, – system identification, – design and analysis of control systems, – adaptive filtering, – non-stationary signal and statistics.
• Spectral tool GUI
• Transfer Function tool GUI
• Probability Density Function tool GUI
![Page 15: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/15.jpg)
FFT analysis of Genera Evolution
• Perform a frequency analysis of genera evolution
![Page 16: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/16.jpg)
Reaction Engineering Lab
• Calibrer et valider des modèles par comparaison avec les résultatsexpérimentaux : estimation de paramètres
• Modèles dépendants de l’espace dans Comsol Multiphysics.
• Modèles de systèmes de réactifs à partir deformules de réactions chimiques.
• Equations de bilan de matière et d’énergiepour les réacteurs idéaux.
![Page 17: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/17.jpg)
The Reaction Engineering Lab
• Reaction Engineering Lab sets up models of reacting systems fromchemical reaction formulas
• It solves material and energy balances for ideal reactors• Parameter estimation allows for calibration and validation of
models to experimental data• It sets up space dependent models in COMSOL Multiphysics
![Page 18: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/18.jpg)
The Reaction Engineering Lab
![Page 19: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/19.jpg)
Reaction Engineering Lab exports to COMSOL Multiphysics
• Transfer the kinetic model and physical properties of the reacting mixture from ideal reactors to space-dependent systems
• Move into detailed reactor analysis and design
( )i
ix Rdx
cud = ( ) iiii RccDt
c =+∇−⋅∇+∂∂
u
![Page 20: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/20.jpg)
Fluid Catalytic Cracking
k3k4
k5
k1
k2
H COheavy crude oil
Gas
LCOlight cycle oil
Coke
Gasolinek6
k7 k8
• Reaction Engineering Lab – Automatic generation of complex kinetics– Ideal reactor modeling
![Page 21: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/21.jpg)
Compression Ignition of Methane
• Reaction Engineering Lab – Import and solve combustion kinetics
![Page 22: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/22.jpg)
5.2 Contrôle Actif avec COMSOL 3.4
• Contrôle d’écoulements– Conditions aux limites (variables de couplage)– Corps solides dans un domaine fluide (ALE,
levelset, autres…)
• Analyse piézo-électrique enAcoustique, MEMS, …
• Interface d’import SPICE en électromagnétisme AC/DC
• Contrôle de systèmes: Signals & Systems Lab
![Page 23: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/23.jpg)
Feedback and Control with COMSOL …
You can simulate your control loops in COMSOL
Full 2D or 3D detail rather than Lumped Mass Models …
• On / Off Control• Proportional Control• Full PID Control
![Page 24: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/24.jpg)
The Problem … Compare Thermal Control Performance: On/Off vs Proportional
• Integration Coupling Variable – TTHERMOCOUPLE
• Feedback – Heatflux Boundary Cond - Peltier
Peltier Cooler: ON
Peltier Cooler: OFF
High Wattage IC Thermocouple
Heat Conducts on Board
Thermocouple Feels Heat
![Page 25: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/25.jpg)
On/Off vs. Proportional Control of Cooler• Feedback through Flux Boundary Condition at Peltier• Qon off = -q_peltier*(T_sense > T_threshold)• Qproportional = -q_peltier*(T_sense - T_threshold) / 10
![Page 26: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/26.jpg)
Related Problem … How about PID Control?
• Control top inlet gas flow rate by measuring concentration• Use new ODE Solver to evaluate time integral part of control
)()()(0
set
t
DsetIsetPin cct
kdtcckccku −∂∂+−+−= ∫
![Page 27: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/27.jpg)
Vous et COMSOL sur www.comsol.fr
• Version d’évaluation sur 2 semaineswww.comsol.fr/contact/evaluation.php
• Documentation, formation, etc– Bibliothèques de modèles : www.comsol.fr/showroom/
– Formations : www.comsol.fr/training/
– Livres : www.comsol.fr/stories/books/ – CD : www.comsol.fr/conference2007/cd/
– Conférence COMSOL 2008
• Support avec [email protected]– Étude gratuite de faisabilité
– Support technique
Vous et COMSOL
![Page 28: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/28.jpg)
• 250 applications• 190 présentations
• 40 modèles
Conférence COMSOL 2008
4-6 novembre 2008 à Hanovre
Conférences COMSOL 2007CD gratuit sur www.comsol.fr
• 300 participants• 98 présentations orales
+ de 150 rdv de support + de 540 places de minicours
En 2007 à Grenoble
Conférences COMSOL
![Page 29: Modélisation OD/1D avec COMSOL Multiphysics 3 2008... · COMSOL & Autres Programmes MATLAB CAD Import Module COMSOL Script SolidWorks 3D CAD 2D CAD dxf iges, sat, etc. ... ChemKin](https://reader033.fdocument.org/reader033/viewer/2022061610/5a7c8f127f8b9a563b8cf78c/html5/thumbnails/29.jpg)
Merci
de votre attention !
Préparez vos questions !!!
Comparaison d’un modèle de transfert thermique 3Davec un modèle réduit 1D.Prise en compte du rayonnement entre les facesinternes du tuyau par intégrales de convolution
Optimisation des Paramètres SPICE décrivant une Diode Semi-conductrice
Graphique en relief de la concentration de trous (porteurs de charge positive, figure de gauche) et superposition descaractéristiques de fonctionnement de la diode et du circuit équivalent (droite).