Utilisateur:Ellande/Brouillon

Ellande

Outils

Discussion

Brouillon 1

Brouillon 2

Brouillon 3

Brouillon 4

Brouillon 5

Brouillon 6

Brouillon 7

Brouillon 8

Brouillon 9

Bac à sable

Détails des calculs de dérivées particulaires

Méthode 1

\frac{d (ρ \vec{v})}{d t} = \frac{\partial (ρ \vec{v})}{\partial t} + (\vec{v} \cdot \vec{g r a d}) (ρ \vec{v}) = \frac{\partial (ρ \vec{v})}{\partial t} + (\vec{v} \cdot \vec{\nabla}) (ρ \vec{v})

\frac{d (ρ \vec{v})}{d t} = ρ \frac{\partial \vec{v}}{\partial t} + \frac{\partial ρ}{\partial t} \vec{v} + ρ v_{x} \frac{\partial \vec{v}}{\partial x} + \frac{\partial ρ}{\partial x} v_{x} \vec{v} + ρ v_{y} \frac{\partial \vec{v}}{\partial y} + \frac{\partial ρ}{\partial y} v_{y} \vec{v} + ρ v_{z} \frac{\partial \vec{v}}{\partial z} + \frac{\partial ρ}{\partial z} v_{z} \vec{v}

\frac{d (ρ \vec{v})}{d t} = ρ \frac{\partial \vec{v}}{\partial t} + \frac{\partial ρ}{\partial t} \vec{v} + ρ (v_{x} \frac{\partial \vec{v}}{\partial x} + v_{y} \frac{\partial \vec{v}}{\partial y} + v_{z} \frac{\partial \vec{v}}{\partial z}) + (\frac{\partial ρ}{\partial x} v_{x} + \frac{\partial ρ}{\partial y} v_{y} + \frac{\partial ρ}{\partial z} v_{z}) \vec{v}

\frac{d (ρ \vec{v})}{d t} = ρ \frac{\partial \vec{v}}{\partial t} + \frac{\partial ρ}{\partial t} \vec{v} + ρ (\vec{v} \cdot \vec{g r a d}) \vec{v} + (\vec{v} \cdot \vec{g r a d} ρ) \vec{v}

\frac{d (ρ \vec{v})}{d t} + ρ \vec{v} div \vec{v} = ρ (\frac{\partial \vec{v}}{\partial t} + (\vec{v} \cdot \vec{g r a d}) \vec{v}) + (\frac{\partial ρ}{\partial t} + \vec{v} \cdot \vec{g r a d} ρ + ρ div \vec{v}) \vec{v}

\frac{d (ρ \vec{v})}{d t} + ρ \vec{v} div \vec{v} = ρ (\frac{\partial \vec{v}}{\partial t} + (\vec{v} \cdot \vec{g r a d}) \vec{v}) + \underset{= 0}{\underset{⏟}{(\frac{\partial ρ}{\partial t} + div (ρ \vec{v}))}} \vec{v}

\frac{d (ρ \vec{v})}{d t} + ρ \vec{v} div \vec{v} = ρ \frac{d \vec{v}}{d t}

Méthode 2

\frac{d (ρ \vec{v})}{d t} = ρ \frac{d \vec{v}}{d t} + \frac{d ρ}{d t} \vec{v}

\frac{d (ρ \vec{v})}{d t} + ρ \vec{v} div \vec{v} = ρ \frac{d \vec{v}}{d t} + \underset{= 0}{\underset{⏟}{(\frac{d ρ}{d t} + ρ div \vec{v})}} \vec{v}

Méthode 3

\frac{d (ρ \vec{v})}{d t} = ρ (\frac{\partial \vec{v}}{\partial t} + v_{x} \frac{\partial \vec{v}}{\partial x} + v_{y} \frac{\partial \vec{v}}{\partial y} + v_{z} \frac{\partial \vec{v}}{\partial z}) + (\frac{\partial ρ}{\partial t} + \frac{\partial ρ}{\partial x} v_{x} + \frac{\partial ρ}{\partial y} v_{y} + \frac{\partial ρ}{\partial z} v_{z}) \vec{v}

\frac{d (ρ \vec{v})}{d t} = ρ (\frac{\partial \vec{v}}{\partial t} + (\vec{v} \cdot \vec{\nabla}) \cdot \vec{v}) + \frac{d ρ}{d t} \vec{v}

\frac{d (ρ \vec{v})}{d t} + ρ \vec{v} div \vec{v} = ρ (\frac{\partial \vec{v}}{\partial t} + (\vec{v} \cdot \vec{\nabla}) \cdot \vec{v}) + ρ \vec{v} div \vec{v}

\frac{d (ρ \vec{v})}{d t} + ρ \vec{v} div \vec{v} = ρ \frac{\partial \vec{v}}{\partial t} + ρ (\vec{v} \cdot \vec{\nabla}) \cdot \vec{v}

Autre

$\vec{\nabla} \cdot (ρ \vec{v} \otimes \vec{v}) = d i v (ρ \vec{v}) \vec{v} + ρ (\vec{v} \cdot \vec{g r a d}) \vec{v}$

$\frac{\partial (ρ \vec{v})}{\partial t} + \vec{\nabla} \cdot (ρ \vec{v} \otimes \vec{v}) = - \vec{\nabla} p + \vec{\nabla} \cdot \overline{\overline{τ}} + ρ \vec{f}$

$\vec{\nabla} \cdot (ρ \vec{v} \otimes \vec{v}) = \vec{\nabla} \cdot ((\begin{matrix} ρ v_{x} \\ ρ v_{y} \\ ρ v_{z} \end{matrix}) \times (\begin{matrix} v_{x} & v_{y} & v_{z} \end{matrix})) = \vec{\nabla} \cdot (\begin{matrix} ρ v_{x} v_{x} & ρ v_{x} v_{y} & ρ v_{x} v_{z} \\ ρ v_{y} v_{x} & ρ v_{y} v_{y} & ρ v_{y} v_{z} \\ ρ v_{z} v_{x} & ρ v_{z} v_{y} & ρ v_{z} v_{z} \end{matrix})$

$\vec{\nabla} \cdot (ρ \vec{v} \otimes \vec{v}) = (\begin{matrix} \frac{\partial (ρ v_{x} v_{x})}{\partial x} + \frac{\partial (ρ v_{x} v_{y})}{\partial y} + \frac{\partial (ρ v_{x} v_{z})}{\partial z} \\ \frac{\partial (ρ v_{y} v_{x})}{\partial x} + \frac{\partial (ρ v_{y} v_{y})}{\partial y} + \frac{\partial (ρ v_{y} v_{z})}{\partial z} \\ \frac{\partial (ρ v_{z} v_{x})}{\partial x} + \frac{\partial (ρ v_{z} v_{y})}{\partial y} + \frac{\partial (ρ v_{z} v_{z})}{\partial z} \end{matrix})$

$\vec{\nabla} \cdot (ρ \vec{v} \otimes \vec{v}) = (\begin{matrix} \frac{\partial ρ}{\partial x} v_{x} v_{x} + ρ \frac{\partial v_{x}}{\partial x} v_{x} + ρ \frac{\partial v_{x}}{\partial x} v_{x} + \frac{\partial ρ}{\partial y} v_{x} v_{y} + ρ \frac{\partial v_{x}}{\partial y} v_{y} + ρ \frac{\partial v_{y}}{\partial y} v_{x} + \frac{\partial ρ}{\partial z} v_{x} v_{z} + ρ \frac{\partial v_{x}}{\partial z} v_{z} + ρ \frac{\partial v_{z}}{\partial z} v_{x} \\ \frac{\partial ρ}{\partial x} v_{y} v_{x} + ρ \frac{\partial v_{y}}{\partial x} v_{x} + ρ \frac{\partial v_{x}}{\partial x} v_{y} + \frac{\partial ρ}{\partial y} v_{y} v_{y} + ρ \frac{\partial v_{y}}{\partial y} v_{y} + ρ \frac{\partial v_{y}}{\partial y} v_{y} + \frac{\partial ρ}{\partial z} v_{y} v_{z} + ρ \frac{\partial v_{y}}{\partial z} v_{z} + ρ \frac{\partial v_{z}}{\partial v_{z}} v_{y} \\ \frac{\partial ρ}{\partial x} v_{z} v_{x} + ρ \frac{\partial v_{z}}{\partial x} v_{x} + ρ \frac{\partial v_{x}}{\partial x} v_{z} + \frac{\partial ρ}{\partial y} v_{z} v_{y} + ρ \frac{\partial v_{z}}{\partial y} v_{y} + ρ \frac{\partial v_{y}}{\partial y} v_{z} + \frac{\partial ρ}{\partial z} v_{z} v_{z} + ρ \frac{\partial v_{z}}{\partial z} v_{z} + ρ \frac{\partial v_{z}}{\partial z} v_{z} \end{matrix})$

$\vec{\nabla} \cdot (ρ \vec{v} \otimes \vec{v}) = (\begin{matrix} (ρ d i v \vec{v} + \vec{g r a d} ρ \cdot \vec{v}) v_{x} + ρ \vec{g r a d} v_{x} \cdot \vec{v} \\ (ρ d i v \vec{v} + \vec{g r a d} ρ \cdot \vec{v}) v_{y} + ρ \vec{g r a d} v_{y} \cdot \vec{v} \\ (ρ d i v \vec{v} + \vec{g r a d} ρ \cdot \vec{v}) v_{z} + ρ \vec{g r a d} v_{z} \cdot \vec{v} \end{matrix})$

$\vec{\nabla} \cdot (ρ \vec{v} \otimes \vec{v}) = (\begin{matrix} d i v (ρ \vec{v}) v_{x} + ρ \vec{g r a d} v_{x} \cdot \vec{v} \\ d i v (ρ \vec{v}) v_{y} + ρ \vec{g r a d} v_{y} \cdot \vec{v} \\ d i v (ρ \vec{v}) v_{z} + ρ \vec{g r a d} v_{z} \cdot \vec{v} \end{matrix})$

$\vec{\nabla} \cdot (ρ \vec{v} \otimes \vec{v}) = d i v (ρ \vec{v}) \vec{v} + ρ (\vec{v} \cdot \vec{g r a d}) \vec{v}$

$\vec{\nabla} \cdot (ρ \vec{v} \otimes \vec{v}) = - \frac{\partial ρ}{\partial t} \vec{v} + ρ (\vec{v} \cdot \vec{g r a d}) \vec{v}$

$\frac{\partial (ρ \vec{v})}{\partial t} + \vec{\nabla} \cdot (ρ \vec{v} \otimes \vec{v}) = ρ \frac{\partial \vec{v}}{\partial t} + \frac{\partial ρ}{\partial t} \vec{v} - \frac{\partial ρ}{\partial t} \vec{v} + ρ (\vec{v} \cdot \vec{g r a d}) \vec{v}$

$\frac{\partial (ρ \vec{v})}{\partial t} + \vec{\nabla} \cdot (ρ \vec{v} \otimes \vec{v}) = ρ \frac{\partial \vec{v}}{\partial t} + ρ (\vec{v} \cdot \vec{g r a d}) \vec{v}$

$(\vec{v} \cdot \vec{\nabla}) ρ \vec{v} = (\begin{matrix} (\vec{v} \cdot \vec{\nabla}) ρ v_{x} \\ (\vec{v} \cdot \vec{\nabla}) ρ v_{y} \\ (\vec{v} \cdot \vec{\nabla}) ρ v_{z} \end{matrix}) = (\begin{matrix} v_{x} \frac{\partial ρ v_{x}}{\partial x} + v_{y} \frac{\partial ρ v_{x}}{\partial y} + v_{z} \frac{\partial ρ v_{x}}{\partial z} \\ v_{x} \frac{\partial ρ v_{y}}{\partial x} + v_{y} \frac{\partial ρ v_{y}}{\partial y} + v_{z} \frac{\partial ρ v_{y}}{\partial z} \\ v_{x} \frac{\partial ρ v_{z}}{\partial x} + v_{y} \frac{\partial ρ v_{z}}{\partial y} + v_{z} \frac{\partial ρ v_{z}}{\partial z} \end{matrix})$

Tenseur des contraintes

Force exercée sur une surface $S$

En ajoutant les forces orientées dans la même direction, la résultante de l'ensemble des forces sur une surface élémentaire d'orientation quelconque s'exprime :

$\vec{d F_{S}} = \overline{\overline{τ}} \times \vec{d S} = (\begin{matrix} σ_{x x} & τ_{x y} & τ_{x z} \\ τ_{y x} & σ_{y y} & τ_{y z} \\ τ_{z x} & τ_{z y} & σ_{z z} \end{matrix}) \times (\begin{matrix} d S_{x} \\ d S_{y} \\ d S_{z} \end{matrix}) = (\begin{matrix} σ_{x x} d S_{x} + τ_{x y} d S_{y} + τ_{x z} d S_{z} \\ τ_{y x} d S_{x} + σ_{y y} d S_{y} + τ_{y z} d S_{z} \\ τ_{z x} d S_{x} + τ_{z y} d S_{y} + σ_{z z} d S_{z} \end{matrix})$ .

Force exercée sur un élément de volume $d V$

Les forces surfaciques qui s'appliquent sur les faces d'un élément de volume $d V = d x d y d z$ sont modélisées par le tenseur des contraintes. L'illustration ci-contre permet de comprendre comment se décomposent ces forces.

Sur chaque face, par convention, le vecteur surface est orienté vers l'extérieur du volume. Par exemple, la face la plus proche de nous sur l'illustration a un vecteur surface $\vec{d S} = d S_{x} \vec{e_{x}} = d y d z \vec{e_{x}}$ . La force qui s'exerce sur cette face peut être décomposée en 3 forces :

une force normale à la surface $\vec{d F_{1}} = σ_{x x} d y d z \vec{e_{x}}$ ;
deux forces dans le plan de la surface :
- $\vec{d F_{2}} = τ_{y x} d y d z \vec{e_{y}}$ dans la direction de l'axe $(O y)$ ;
- $\vec{d F_{2}} = τ_{z x} d y d z \vec{e_{z}}$ dans la direction de l'axe $(O y)$ .

Les indices associés à chaque contrainte indique, dans l'ordre, la direction de la force et la face sur laquelle la force s'applique. Les contraintes situées sur la diagonale correspondent à des forces de pression ce qui justifie que l'on leur affecte un nom différent. Les autres correspondent à des contraintes de cisaillement dues à la viscosité dans le cas de la mécanique des fluides.

La résultante des forces qui s'exerce sur l'élément de volume peut s'écrire :

$\vec{d F_{V}} = \sum_{i} \vec{d F_{i}} = (\begin{matrix} [σ_{x x} (x + d x) - σ_{x x} (x)] d y d z + [τ_{x y} (y + d y) - τ_{x y} (y)] d x d z + [τ_{x z} (z + d z) - τ_{x z} (z)] d x d y \\ [τ_{y x} (x + d x) - τ_{y x} (x)] d y d z + [σ_{y y} (y + d y) - σ_{y y} (y)] d x d z + [τ_{y z} (z + d z) - τ_{y z} (z)] d x d y \\ [τ_{z x} (x + d x) - τ_{z x} (x)] d y d z + [τ_{z y} (y + d y) - τ_{z y} (y)] d x d z + [σ_{z z} (z + d z) - σ_{z z} (z)] d x d y \end{matrix})$ ,

$\vec{d F_{V}} = (\begin{matrix} \frac{\partial σ_{x x}}{\partial x} d x d y d z + \frac{\partial τ_{x y}}{\partial y} d x d y d z + \frac{\partial τ_{x z}}{\partial z} d x d y d z \\ \frac{\partial τ_{y x}}{\partial y} d x d y d z + \frac{\partial σ_{y y}}{\partial y} d x d y d z + \frac{\partial τ_{y z}}{\partial y} d x d y d z \\ \frac{\partial τ_{z x}}{\partial y} d x d y d z + \frac{\partial τ_{z y}}{\partial y} d x d y d z + \frac{\partial σ_{z z}}{\partial y} d x d y d z \end{matrix})$

$\vec{d F_{V}} = \vec{d i v} (\begin{matrix} σ_{x x} & τ_{x y} & τ_{x z} \\ τ_{y x} & σ_{y y} & τ_{y z} \\ τ_{z x} & τ_{z y} & σ_{z z} \end{matrix}) d V = \vec{d i v} \overline{\overline{τ}} d V$

Décomposition

$\overline{\overline{σ}} = - P \overline{\overline{δ}} + \overline{\overline{τ}}$

$(\begin{matrix} σ_{x x} & τ_{x y} & τ_{x z} \\ τ_{y x} & σ_{y y} & τ_{y z} \\ τ_{z x} & τ_{z y} & σ_{z z} \end{matrix}) = - P (\begin{matrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{matrix}) + (\begin{matrix} τ_{x x} & τ_{x y} & τ_{x z} \\ τ_{y x} & τ_{y y} & τ_{y z} \\ τ_{z x} & τ_{z y} & τ_{z z} \end{matrix})$

$τ_{x x} + τ_{y y} + τ_{z z} = 0$

$τ_{i j} = μ (\frac{\partial v_{i}}{\partial x_{j}} + \frac{\partial v_{j}}{\partial x_{i}}) + μ^{'} d i v \vec{v} δ_{i j}$

$μ$ : viscosité dynamique

$μ^{'}$ : coefficient de seconde viscosité

tenseur de contraintes

Recherche ː travail des forces surfaciques

eq de l'énergie ; Tenseur des contraintes

Produit matriciel ; Matrice transposée

$\vec{div} \overline{\overline{τ}} = (\begin{matrix} \frac{\partial σ_{x x}}{\partial x} + \frac{\partial τ_{x y}}{\partial y} + \frac{\partial τ_{x z}}{\partial z} \\ \frac{\partial τ_{y x}}{\partial x} + \frac{\partial σ_{y y}}{\partial y} + \frac{\partial τ_{y z}}{\partial z} \\ \frac{\partial τ_{z x}}{\partial x} + \frac{\partial τ_{z y}}{\partial y} + \frac{\partial σ_{z z}}{\partial z} \end{matrix})$

Travail des forces de surface

forces dont la direction est selon l'axe y

$\begin{matrix} \sum W_{S_{y}} = & [(τ_{y x} v_{y})_{x + d x} - (τ_{y x} v_{y})_{x}] d y d z \\ + [(σ_{y y} v_{y})_{y + d y} - (σ_{y y} v_{y})_{y}] d x d z \\ + [(τ_{y z} v_{y})_{z + d z} - (τ_{y z} v_{y})_{z}] d x d y \end{matrix}$

$\sum W_{S_{y}} = (\frac{\partial (τ_{y x} v_{y})}{\partial x} + \frac{\partial (σ_{y y} v_{y})}{\partial y} + \frac{\partial (τ_{y z} v_{y})}{\partial z}) d V$

$\begin{matrix} \sum W_{S} = & (\frac{\partial (σ_{x x} v_{x})}{\partial x} + \frac{\partial (τ_{x y} v_{x})}{\partial y} + \frac{\partial (τ_{x z} v_{x})}{\partial z} \\ + \frac{\partial (τ_{y x} v_{y})}{\partial x} + \frac{\partial (σ_{y y} v_{y})}{\partial y} + \frac{\partial (τ_{y z} v_{y})}{\partial z} \\ + \frac{\partial (τ_{z x} v_{z})}{\partial x} + \frac{\partial (τ_{z y y} v_{z})}{\partial y} + \frac{\partial (σ_{z z} v_{z})}{\partial z}) d V \end{matrix}$

$\sum W_{S} = div (^{t} \overline{\overline{τ}} \times \vec{v}) d V = div (\begin{matrix} σ_{x x} v_{x} + τ_{y x} v_{y} + τ_{z x} v_{z} \\ τ_{x y} v_{x} + σ_{y y} v_{y} + τ_{z y} v_{z} \\ τ_{x z} v_{x} + τ_{y z} v_{y} + σ_{z z} v_{z} \end{matrix}) d V$

$\begin{matrix} \sum W_{S} = & (\frac{\partial σ_{x x}}{\partial x} v_{x} + σ_{x x} \frac{\partial v_{x}}{\partial x} + \frac{\partial τ_{x y}}{\partial y} v_{x} + τ_{x y} \frac{\partial v_{x}}{\partial y} + \frac{\partial τ_{x z}}{\partial z} v_{x} + τ_{x z} \frac{\partial v_{x}}{\partial z} \\ + \frac{\partial τ_{y x}}{\partial x} v_{y} + τ_{y x} \frac{\partial v_{y}}{\partial x} + \frac{\partial σ_{y y}}{\partial y} v_{y} + σ_{y y} \frac{\partial v_{y}}{\partial y} + \frac{\partial τ_{y z}}{\partial z} v_{y} + τ_{y z} \frac{\partial v_{y}}{\partial z} \\ + \frac{\partial τ_{z x}}{\partial x} v_{z} + τ_{z x} \frac{\partial v_{z}}{\partial x} + \frac{\partial τ_{z y y}}{\partial y} v_{z} + τ_{z y y} \frac{\partial v_{z}}{\partial y} + \frac{\partial σ_{z z}}{\partial z} v_{z} + σ_{z z} \frac{\partial v_{z}}{\partial z}) d V \end{matrix}$

$\begin{matrix} \sum W_{S} = \vec{v} \cdot \vec{div} \overline{\overline{𝝉}} + & (σ_{x x} \frac{\partial v_{x}}{\partial x} + τ_{x y} \frac{\partial v_{x}}{\partial y} + τ_{x z} \frac{\partial v_{x}}{\partial z} \\ + τ_{y x} \frac{\partial v_{y}}{\partial x} + σ_{y y} \frac{\partial v_{y}}{\partial y} + τ_{y z} \frac{\partial v_{y}}{\partial z} \\ + τ_{z x} \frac{\partial v_{z}}{\partial x} + τ_{z y y} \frac{\partial v_{z}}{\partial y} + σ_{z z} \frac{\partial v_{z}}{\partial z}) d V \end{matrix}$

Euler

En utilisant la dérivée particulaire :

\frac{d (ρ_{tot} 𝐯)}{d t} = \frac{\partial (ρ_{tot} 𝐯)}{\partial t} + (𝐯 \cdot 𝜵) (ρ_{tot} 𝐯) = ρ_{tot} \frac{\partial 𝐯}{\partial t} + ρ_{tot} (𝐯 \cdot 𝜵) 𝐯 + \frac{\partial ρ_{tot}}{\partial t} 𝐯 + 𝐯 (𝐯 \cdot 𝜵) ρ_{tot}

.

On ne tient compte que des deux premiers termes qui correspondent à la variation de quantité de mouvement par accélération locale $(ρ_{tot} \frac{\partial 𝐯}{\partial t})$ et par variation de masse locale $(\frac{\partial ρ_{tot}}{\partial t} 𝐯)$ . Les termes convectifs sont négligés.

La variation par accélération convective $ρ (𝐯 \cdot 𝜵) 𝐯$ est négligeable.
La variation par variation de masse convective $𝐯 (𝐯 \cdot 𝜵) 𝝆$ est négligeable.

La conservation de la quantité de mouvement s'écrit alors :

\frac{\partial (ρ_{tot} 𝐯)}{\partial t} = - \nabla P_{tot} . (3)