Pavel

From KratosWiki
(Difference between revisions)
Jump to: navigation, search
(Introduction)
(Convection Diffusion module)
Line 1: Line 1:
== <span style="color:#FF0000">  Convection Diffusion module  </span> ==
 
 
 
=== <span style="color:#0000FF"> Introduction </span>===
 
=== <span style="color:#0000FF"> Introduction </span>===
  
The numerical solution of convection-diffusion transport problems arises in many important applications in science and engineering. These problems occur in many applications such as in the transport of air and ground water pollutants, oil reservoir flow, in the modeling of semiconductors, and so forth.  
+
The numerical solution of convection-diffusion transport problems arises in many important applications in science and engineering. These problems occur in many applications such as in the transport of air and ground water pollutants, oil reservoir flow, in the modeling of semiconductors, and so forth. This paper describes the Convection Diffusion Applications for solving this equation.
 +
 
 +
Theory
 +
 
 +
Under the assumption of incompressibility, the governing equations are given by
 +
 
 +
[[Image:placa.jpg|thumb|right|300px|Problem description]]
 +
<math>
 +
  \rho C \frac{\partial T}{\partial t} + \rho C v \cdot \nabla T + \nabla
 +
  \cdot q = 0
 +
</math>
 +
(1)
 +
 
  
Some examples are shown below:
 
  
{| class="wikitable" width="100%" style="text-align:left; background:#d0d9dd; border:0px solid #e1eaee; font-size:100%; -moz-border-radius-topleft:0px; -moz-border-radius-bottomleft:0px; padding:0px 0px 0px 0px;" valign="top"
+
<math>
!
+
  q = -\kappa \nabla \cdot T
!
+
</math>
!
+
(2)
!
+
|-style="background:#F1FAFF;"
+
| [[Image:shape.jpg|200px]]
+
| Heat transfer: convection-diffusion
+
ADVERTISMENT STYLE no numerical details!!!
+
|-style="background:#F1FAFF;"
+
| [[Image:shape.jpg|200px]]
+
| Brief description of what the model means
+
|-style="background:#F1FAFF;"
+
| [[Image:shape.jpg|200px]]
+
| Brief description of what the model means
+
|-style="background:#F1FAFF;"
+
| [[Image:shape.jpg|200px]]
+
| Brief description of what the model means
+
|}
+
  
 +
In the context of mass difussion within a fluid, (1) is is the mass
 +
conservation equation and (2) is a contitutive law proposed by Fourier. The
 +
notation is standard: <math>\rho</math> is the density, <math>C</math> the heat capacity, <math>\kappa </math>
 +
the thermal conductivity, <math>T</math> is the temperature, <math>v</math> is the velocity field
 +
and <math>q</math> is the diffusive flux per unit fluid density.
  
  
 +
Remark: systen can be decoupled since we can plug (2) into (1) and solve the
 +
scalar equation
  
Description of the underlying theory and schematic list of the problems this application can solve.
+
<math>
 +
  \rho C \frac{\partial T}{\partial t} + \rho Cv \cdot \nabla T -
 +
  \nabla \cdot (\kappa \nabla \cdot T) = 0
 +
</math>
 +
(3)
  
In this application the Navier Stokes equations are solved bla bla..
 
  
=== <span style="color:#0000FF"> Structure </span> ===
+
Problem statement
  
====  Fluid types ====
+
Let us consider the transport by convection and diffusion in an open set
* '''Incompressible''' fluid
+
<math>\Omega</math> <math>\subset \Re^d</math>(d=2 or 3) \ with piecewise smooth boundary <math>\Gamma</math>,
aaa
+
such that <math>\Gamma = \Gamma_d \cup \Gamma_N</math>. The unit outward normal vector to
* '''Compressible''' fluid
+
<math>\Gamma</math> is denoted <math>n</math>. The convection-diffusion initial-boundary value
 +
problem can be stated as follows: given a divergence-free velocity field a,
 +
the diffusion tensor <math>\kappa</math> and adequate initial and boundary conditions, find <math>T</math> :
 +
<math>\bar{\Omega} \times \left[ 0, T \left] \rightarrow \Re \right. \right.</math> such
 +
that
  
'''Constitutive laws'''
+
<math>
* '''Newtonian'''
+
\rho C \frac{\partial T}{\partial t} + \rho Cv \cdot \nabla T -
* '''No-Newtonian'''
+
  \nabla \cdot (\kappa \nabla T) = 0 in \Omega \times (0, T)
** Bingham plastics
+
</math>
** Variable yield model
+
(4)
  
====  Kinematical approaches ====
 
  
* '''Eulerian'''
+
<math>
** With '''free surface''' (level set)
+
T(x, 0) = T_0 (x) on \Omega
** Without free surface
+
</math>
* Lagrangian '''PFEM''' (implicitly with free surface)
+
  
==== Solution strategy  ====
+
<math>
 +
T = T_D on \Gamma_D \times (0, T)
 +
</math>
  
* '''Fractional step'''
+
<math>
* '''Monolithic'''
+
k (\nabla T) \cdot n = q on \Gamma_N \times (0, T)
Different solvers are availables (LINK TO SOLVER SECTION!!!!)
+
</math>
  
  
In both cases a Newton Raphson residual based strategy is used for linearizing the problem.
 
  
==== Elements ====
+
Space discretization method
Linear triangular elements in 2D and linear tetrahedra elements in 3D.
+
Multiplying Eq.(4) by a test function <math>W</math> and intehrating on the whole domain
 +
<math>\Omega</math> the equation reads
  
Stabilization techniques availables:
 
* ASGS
 
* OSS
 
  
==== Boundary conditions ====
+
<math>\int_\Omega \rho C \frac{\partial T}{\partial t} W d V + \int_\Omega \rho Cv \cdot
 +
  \nabla T W d V = \int_\Omega (\nabla \cdot (k \nabla T)) W d V
 +
</math>
 +
(5)
  
* Velocity boundary condition: Inlet of water
+
Integratin by parts the right term of Eq.(5) leads to
* Pressure boundary condition: Pressure can be imposed strongly or weakly...
+
* Wall boundary condition:
+
** Slip/no slip boundary condition
+
** Wall law
+
* Flag variable?????
+
  
==== Initial conditions ====
+
<math> \int_\Omega \rho C \frac{\partial T}{\partial t} W d V + \int_\Omega v \cdot \nabla T W d
Initial condition both in velocity and pressure can be set.
+
  V = - \int_\Omega k \nabla T \cdot \nabla W d V</math> (6)
  
==== Turbulence models ====
+
Finite element discretization
The user can chose wether to use or not a turbulence model.
+
Those available in kratos are:
+
* Smagorinsky-Lily
+
* Spalart-Allmaras
+
  
==== HPC ====
+
The temperature is discretized in the standard finite element method manner as
The code can be run in shared or distributed memory:
+
<math>
* OpenMP:
+
T = \sum N_i T_i
* MPI:
+
</math>(7)
 +
where <math>N_i</math> are the nodal shape functions. Substituting the finite element
 +
approximation (7) into the variational equation () and choosing a Galerling formulation (<math>W_i=N_i</math>) leads to the following equation:  
  
==== Problem parameters ====
+
<math> \rho C M \frac{\partial T}{\partial t} + \rho C S T
 +
= - \kappa L T </math> (7)
  
==== Others relevand aspects ====
+
Time discretization method
 +
Consider a first-order BDF (that is, the Euler implicit scheme)
  
=== <span style="color:#0000FF"> Benchmarking </span>===
 
  
=== <span style="color:#0000FF"> Tutorials </span>===
+
<math>
 +
\frac{\partial T}{\partial t}=\frac{(T^{n + 1} - T^n)}{\Delta t}
 +
</math>
 +
or a second-order BDF
  
=== <span style="color:#0000FF"> Contact people </span>===
 
  
=== <span style="color:#0000FF"> Akcnowledgements </span>===
+
<math>
 +
\frac{\partial T}{\partial t}= \frac{1}{2} \frac{(3 \times T^{n + 1} - 4 \times T^n + T^{n-1})}{\Delta t}
 +
</math>
  
 
== <span style="color:#FF0000"> Computational Structural Mechanics module </span> ==
 
== <span style="color:#FF0000"> Computational Structural Mechanics module </span> ==

Revision as of 14:21, 16 July 2013

Contents

Introduction

The numerical solution of convection-diffusion transport problems arises in many important applications in science and engineering. These problems occur in many applications such as in the transport of air and ground water pollutants, oil reservoir flow, in the modeling of semiconductors, and so forth. This paper describes the Convection Diffusion Applications for solving this equation.

Theory

Under the assumption of incompressibility, the governing equations are given by

Problem description


  \rho C \frac{\partial T}{\partial t} + \rho C v \cdot \nabla T + \nabla
  \cdot q = 0 
(1)



  q = -\kappa \nabla \cdot T 
(2)

In the context of mass difussion within a fluid, (1) is is the mass conservation equation and (2) is a contitutive law proposed by Fourier. The notation is standard: ρ is the density, C the heat capacity, κ the thermal conductivity, T is the temperature, v is the velocity field and q is the diffusive flux per unit fluid density.


Remark: systen can be decoupled since we can plug (2) into (1) and solve the scalar equation


  \rho C \frac{\partial T}{\partial t} + \rho Cv \cdot \nabla T -
  \nabla \cdot (\kappa \nabla \cdot T) = 0 
(3)


Problem statement

Let us consider the transport by convection and diffusion in an open set Ω \subset \Re^d(d=2 or 3) \ with piecewise smooth boundary Γ, such that \Gamma = \Gamma_d \cup \Gamma_N. The unit outward normal vector to Γ is denoted n. The convection-diffusion initial-boundary value problem can be stated as follows: given a divergence-free velocity field a, the diffusion tensor κ and adequate initial and boundary conditions, find T : \bar{\Omega} \times \left[ 0, T \left] \rightarrow \Re \right. \right. such that


\rho C \frac{\partial T}{\partial t} + \rho Cv \cdot \nabla T -
   \nabla \cdot (\kappa \nabla T) = 0 in \Omega \times (0, T) 
(4)


T(x,0) = T0(x)onΩ


T = T_D on \Gamma_D \times (0, T)


k (\nabla T) \cdot n = q on \Gamma_N \times (0, T)


Space discretization method Multiplying Eq.(4) by a test function W and intehrating on the whole domain Ω the equation reads


\int_\Omega \rho C \frac{\partial T}{\partial t} W d V + \int_\Omega \rho Cv \cdot
   \nabla T W d V = \int_\Omega (\nabla \cdot (k \nabla T)) W d V 
(5)

Integratin by parts the right term of Eq.(5) leads to

 \int_\Omega \rho C \frac{\partial T}{\partial t} W d V + \int_\Omega v \cdot \nabla T W d
   V = - \int_\Omega k \nabla T \cdot \nabla W d V (6)

Finite element discretization

The temperature is discretized in the standard finite element method manner as 
 T = \sum N_i T_i 
(7) where Ni are the nodal shape functions. Substituting the finite element approximation (7) into the variational equation () and choosing a Galerling formulation (Wi = Ni) leads to the following equation:

 \rho C M \frac{\partial T}{\partial t} + \rho C S T 
 = - \kappa L T (7)

Time discretization method Consider a first-order BDF (that is, the Euler implicit scheme)



\frac{\partial T}{\partial t}=\frac{(T^{n + 1} - T^n)}{\Delta t}

or a second-order BDF



\frac{\partial T}{\partial t}= \frac{1}{2} \frac{(3 \times T^{n + 1} - 4 \times T^n + T^{n-1})}{\Delta t}

Computational Structural Mechanics module

Introduction

Examples showing the class of problems that the code can solve (2-4 examples)

Description of the underlying theory and schematic list of the problems this application can solve.


The Computational Structural Mechanics module (CSM) is....

Application Structure

Analysis Type

The available solutions strategies are:

  • Static
  • Dynamic
  • Relaxed dynamic

With this module you can solve both linear and non linear problems. In case of non linear problems several methods are available:

  • Newton-Raphson
  • Newton Raphson with line search
  • Arch lenght

Different solvers are availables (LINK TO SOLVER SECTION!!!!)


Elements

  • Frame Elements:
    • Euler-Bernoulli beam short explanation
    • Crisfield truss short explanation
  • 2D elements
    • Linear triangular element:
  • Shell elements:
    • Isotropic shell: (change the name with the usual one!!!!)
    • Ansotropic shell: (change the name with the usual one!!!!)
    • EBST shell: (change the name with the usual one!!!!)
  • Membrane element:
  • Solid elements:
    • Linear tetrahedral element:
Dimension Element Type Kratos name Geometry Nonlinearity Material Type
1D Frame LinearBeamElement Line Isotropic
1D Truss CrisfieldTrussElement Line Large Displacement Isotropic
2D Solid TotalLagrangian 2D Geometries Large Displacement Isotropic
3D Solid TotalLagrangian 3D Geometries Large Displacement Isotropic
Shell ShellIsotropic 3D Triangle Large Displacement Isotropic
Shell ShellAnisotropic 3D Triangle Large Displacement Orthotropic

Boundary Conditions

Boundary conditions can be set fixing displacements and rotations degrees of freedom.

Loads

  • Self weight
  • Punctual force
  • Moment
  • Face pressure (sign convenction!!!!)
  • Distributed load

Constitutive laws

The following constitutive laws are available:

  • Linear elastic:
  • ...

HPC

The code can be run in shared or distributed memory:

  • OpenMP:
  • MPI:

Problem parameters

...

Others relevand aspects

...

Benchmarking

Here validation and verification examples should be inserted

Tutorials

Contact people

Akcnowledgements

Convection Diffusion module

Introduction

The numerical solution of convection-diffusion transport problems arises in many important applications in science and engineering. These problems occur in many applications such as in the transport of air and ground water pollutants, oil reservoir flow, in the modeling of semiconductors, and so forth. This paper describes the Convection Diffusion Applications for solving this equation.

Theory

Under the assumption of incompressibility, the governing equations are given by

Problem description


  \rho C \frac{\partial T}{\partial t} + \rho C v \cdot \nabla T + \nabla
  \cdot q = 0 
(1)



  q = -\kappa \nabla \cdot T 
(2)

In the context of mass difussion within a fluid, (1) is is the mass conservation equation and (2) is a contitutive law proposed by Fourier. The notation is standard: ρ is the density, C the heat capacity, κ the thermal conductivity, T is the temperature, v is the velocity field and q is the diffusive flux per unit fluid density.


Remark: systen can be decoupled since we can plug (2) into (1) and solve the scalar equation


  \rho C \frac{\partial T}{\partial t} + \rho Cv \cdot \nabla T -
  \nabla \cdot (\kappa \nabla \cdot T) = 0 
(3)


Problem statement

Let us consider the transport by convection and diffusion in an open set Ω \subset \Re^d(d=2 or 3) \ with piecewise smooth boundary Γ, such that \Gamma = \Gamma_d \cup \Gamma_N. The unit outward normal vector to Γ is denoted n. The convection-diffusion initial-boundary value problem can be stated as follows: given a divergence-free velocity field a, the diffusion tensor κ and adequate initial and boundary conditions, find T : \bar{\Omega} \times \left[ 0, T \left] \rightarrow \Re \right. \right. such that


\rho C \frac{\partial T}{\partial t} + \rho Cv \cdot \nabla T -
   \nabla \cdot (\kappa \nabla T) = 0 in \Omega \times (0, T) 
(4)


T(x,0) = T0(x)onΩ


T = T_D on \Gamma_D \times (0, T)


k (\nabla T) \cdot n = q on \Gamma_N \times (0, T)


Space discretization method Multiplying Eq.(4) by a test function W and intehrating on the whole domain Ω the equation reads


\int_\Omega \rho C \frac{\partial T}{\partial t} W d V + \int_\Omega \rho Cv \cdot
   \nabla T W d V = \int_\Omega (\nabla \cdot (k \nabla T)) W d V 
(5)

Integratin by parts the right term of Eq.(5) leads to

 \int_\Omega \rho C \frac{\partial T}{\partial t} W d V + \int_\Omega v \cdot \nabla T W d
   V = - \int_\Omega k \nabla T \cdot \nabla W d V (6)

Finite element discretization

The temperature is discretized in the standard finite element method manner as 
 T = \sum N_i T_i 
(7) where Ni are the nodal shape functions. Substituting the finite element approximation (7) into the variational equation () and choosing a Galerling formulation (Wi = Ni) leads to the following equation:

 \rho C M \frac{\partial T}{\partial t} + \rho C S T 
 = - \kappa L T (7)

Time discretization method Consider a first-order BDF (that is, the Euler implicit scheme)



\frac{\partial T}{\partial t}=\frac{(T^{n + 1} - T^n)}{\Delta t}

or a second-order BDF



\frac{\partial T}{\partial t}= \frac{1}{2} \frac{(3 \times T^{n + 1} - 4 \times T^n + T^{n-1})}{\Delta t}

Structure

Analysis type

The available solution strategy is:

   Dynamic

With this module you can solve both linear and non linear problems.

Kinematical approaches

Eulerian and Lagrangian approach are available in order to solve the equation.

Solution strategies

Elements

Linear triangular elements in 2D and linear tetrahedra elements in 3D. Both elements are stabilized with OSS.

ConvDiff2D

ConvDiff3D

Boundary conditions

Dirichlet boundary condition:


T = T_D on \Gamma_D \times (0, T)

Neumann boundary conditions: 
k (\nabla T) \cdot n = q on \Gamma_N \times (0, T)

Initial conditions

Initial condition in temperature can be set.

HPC

The code can be run in shared or distributed memory:

  • OpenMP:
  • MPI:

Problem parameters

The parameters involved in this problem are:

ρ : Density

C :heat capacity

κ: thermal conductivity

v : velocity field

q: diffusive flux per unit fluid density.

T a: ambient temperature.

σ: Stefen Boltzmann constant

e: emissivity

h: convection coefficient

Others relevand aspects

Benchmarking

Tutorials

Contact people

Akcnowledgements

Personal tools
Categories