GeneralTemplate

From KratosWiki
(Difference between revisions)
Jump to: navigation, search
(Created page with "== <span style="color:#FF0000"> Computational Fluid Dynamics module </span> == === <span style="color:#0000FF"> Introduction </span>=== Examples showing the class of problem...")
 
 
(3 intermediate revisions by one user not shown)
Line 1: Line 1:
== <span style="color:#FF0000">  Computational Fluid Dynamics module </span> ==
+
== <span style="color:#FF0000">  Computational Fluid Dynamics module (CFDM) </span> ==
  
=== <span style="color:#0000FF"> Introduction </span>===
+
=== <span style="color:#0000FF"> CFDM Introduction </span>===
  
  
Line 33: Line 33:
 
In this application the Navier Stokes equations are solved bla bla..
 
In this application the Navier Stokes equations are solved bla bla..
  
=== <span style="color:#0000FF"> Structure </span> ===
+
=== <span style="color:#0000FF"> CFDM Structure </span> ===
  
====  Fluid types ====
+
====  CFDM Fluid types ====
 
* '''Incompressible''' fluid
 
* '''Incompressible''' fluid
 
aaa
 
aaa
Line 46: Line 46:
 
** Variable yield model
 
** Variable yield model
  
====  Kinematical approaches ====
+
====  CFDM Kinematical approaches ====
  
 
* '''Eulerian'''  
 
* '''Eulerian'''  
Line 53: Line 53:
 
* Lagrangian '''PFEM''' (implicitly with free surface)
 
* Lagrangian '''PFEM''' (implicitly with free surface)
  
==== Solution strategy  ====
+
==== CFDM Solution strategy  ====
  
 
* '''Fractional step'''
 
* '''Fractional step'''
Line 62: Line 62:
 
In both cases a Newton Raphson residual based strategy is used for linearizing the problem.
 
In both cases a Newton Raphson residual based strategy is used for linearizing the problem.
  
==== Elements ====
+
==== CFDM Elements ====
 
Linear triangular elements in 2D and linear tetrahedra elements in 3D.
 
Linear triangular elements in 2D and linear tetrahedra elements in 3D.
  
Line 69: Line 69:
 
* OSS
 
* OSS
  
==== Boundary conditions ====
+
==== CFDM Boundary conditions ====
  
 
* Velocity boundary condition: Inlet of water
 
* Velocity boundary condition: Inlet of water
Line 78: Line 78:
 
* Flag variable?????
 
* Flag variable?????
  
==== Initial conditions ====
+
==== CFDM Initial conditions ====
 
Initial condition both in velocity and pressure can be set.
 
Initial condition both in velocity and pressure can be set.
  
==== Turbulence models ====
+
==== CFDM Turbulence models ====
 
The user can chose wether to use or not a turbulence model.
 
The user can chose wether to use or not a turbulence model.
 
Those available in kratos are:
 
Those available in kratos are:
Line 87: Line 87:
 
* Spalart-Allmaras
 
* Spalart-Allmaras
  
==== HPC ====
+
==== CFDM HPC ====
 
The code can be run in shared or distributed memory:
 
The code can be run in shared or distributed memory:
 
* OpenMP:  
 
* OpenMP:  
 
* MPI:
 
* MPI:
  
==== Problem parameters ====
+
==== CFDM Problem parameters ====
  
==== Others relevand aspects ====
+
==== CFDM Others relevand aspects ====
  
=== <span style="color:#0000FF"> Benchmarking </span>===
+
=== <span style="color:#0000FF"> CFDM Benchmarking </span>===
  
=== <span style="color:#0000FF"> Tutorials </span>===
+
=== <span style="color:#0000FF"> CFDM Tutorials </span>===
  
=== <span style="color:#0000FF"> Contact people </span>===
+
=== <span style="color:#0000FF"> CFDM Contact people </span>===
  
=== <span style="color:#0000FF"> Akcnowledgements </span>===
+
=== <span style="color:#0000FF"> CFDM Akcnowledgements </span>===
  
== <span style="color:#FF0000"> Computational Structural Mechanics module </span> ==
 
  
=== <span style="color:#0000FF"> Introduction </span> ===
+
 
 +
== <span style="color:#FF0000"> Computational Structural Mechanics module (CSMM) </span> ==
 +
 
 +
=== <span style="color:#0000FF"> CSMM Introduction </span> ===
  
 
Examples showing the class of problems that the code can solve (2-4 examples)
 
Examples showing the class of problems that the code can solve (2-4 examples)
Line 115: Line 117:
 
The Computational Structural Mechanics module (CSM) is....
 
The Computational Structural Mechanics module (CSM) is....
  
=== <span style="color:#0000FF"> Application Structure </span>===
+
=== <span style="color:#0000FF"> CSMM Application Structure </span>===
  
==== Analysis Type ====
+
==== CSMM Analysis Type ====
 
The available solutions strategies are:
 
The available solutions strategies are:
 
* '''Static'''
 
* '''Static'''
Line 132: Line 134:
  
  
==== Elements ====
+
==== CSMM Elements ====
 
* '''Frame Elements''':  
 
* '''Frame Elements''':  
 
** '''Euler-Bernoulli beam'''  short explanation
 
** '''Euler-Bernoulli beam'''  short explanation
Line 199: Line 201:
 
|}
 
|}
  
==== Boundary Conditions ====
+
==== CSMM Boundary Conditions ====
  
 
Boundary conditions can be set fixing displacements and rotations degrees of freedom.
 
Boundary conditions can be set fixing displacements and rotations degrees of freedom.
  
==== Loads ====
+
==== CSMM Loads ====
  
 
* Self weight  
 
* Self weight  
Line 211: Line 213:
 
* Distributed load
 
* Distributed load
  
==== Constitutive laws ====
+
==== CSMM Constitutive laws ====
  
 
The following constitutive laws are available:
 
The following constitutive laws are available:
Line 218: Line 220:
 
* ...
 
* ...
  
==== HPC ====
+
==== CSMM HPC ====
 
The code can be run in shared or distributed memory:
 
The code can be run in shared or distributed memory:
 
* OpenMP:  
 
* OpenMP:  
 
* MPI:
 
* MPI:
  
==== Problem parameters ====
+
==== CSMM Problem parameters ====
  
 
...
 
...
  
==== Others relevand aspects ====
+
==== CSMM Others relevand aspects ====
  
 
...
 
...
  
=== <span style="color:#0000FF"> Benchmarking </span>===
+
=== <span style="color:#0000FF"> CSMM Benchmarking </span>===
  
 
Here validation and verification examples should be inserted
 
Here validation and verification examples should be inserted
  
=== <span style="color:#0000FF"> Tutorials </span>===
+
=== <span style="color:#0000FF"> CSMM Tutorials </span>===
 +
 
 +
=== <span style="color:#0000FF"> CSMM Contact people </span>===
 +
 
 +
=== <span style="color:#0000FF">CSMM  Akcnowledgements </span>===
  
=== <span style="color:#0000FF"> Contact people </span>===
 
  
=== <span style="color:#0000FF"> Akcnowledgements </span>===
 
  
== <span style="color:#FF0000"> Convection Diffusion module </span>==
+
== <span style="color:#FF0000"> Convection Diffusion module (CDM)</span>==
  
=== <span style="color:#0000FF"> Introduction </span>===
+
=== <span style="color:#0000FF"> CDM 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. This paper describes the Convection Diffusion Applications for solving this equation.
 
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.
Line 355: Line 359:
 
</math>
 
</math>
  
=== <span style="color:#0000FF"> Structure </span>===
+
=== <span style="color:#0000FF"> CDM Structure </span>===
 
==== Analysis type ====
 
==== Analysis type ====
 
The available solution strategy is:
 
The available solution strategy is:
Line 361: Line 365:
 
With this module you can solve both linear and non linear problems.
 
With this module you can solve both linear and non linear problems.
  
==== Kinematical approaches ====
+
==== CDM Kinematical approaches ====
 
Eulerian and Lagrangian approach are available in order to solve the equation.
 
Eulerian and Lagrangian approach are available in order to solve the equation.
  
==== Solution strategies ====
+
==== CDM Solution strategies ====
  
==== Elements ====
+
==== CDM Elements ====
 
Linear triangular elements in 2D and linear tetrahedra elements in 3D. Both elements are stabilized with OSS.  
 
Linear triangular elements in 2D and linear tetrahedra elements in 3D. Both elements are stabilized with OSS.  
  
Line 373: Line 377:
 
ConvDiff3D
 
ConvDiff3D
  
==== Boundary conditions ====
+
==== CDM Boundary conditions ====
  
 
Dirichlet boundary condition:  
 
Dirichlet boundary condition:  
Line 386: Line 390:
 
</math>
 
</math>
  
==== Initial conditions ====
+
==== CDM Initial conditions ====
  
 
Initial condition in temperature can be set.
 
Initial condition in temperature can be set.
  
==== HPC ====
+
==== CDM HPC ====
 
The code can be run in shared or distributed memory:
 
The code can be run in shared or distributed memory:
 
* OpenMP:  
 
* OpenMP:  
 
* MPI:
 
* MPI:
  
==== Problem parameters ====
+
==== CDM Problem parameters ====
 
The parameters involved in this problem are:
 
The parameters involved in this problem are:
  
Line 416: Line 420:
 
<math>h</math>: convection coefficient
 
<math>h</math>: convection coefficient
  
==== Others relevand aspects ====
+
==== CDM Others relevand aspects ====
  
=== <span style="color:#0000FF"> Benchmarking </span>===
+
=== <span style="color:#0000FF"> CDM Benchmarking </span>===
  
=== <span style="color:#0000FF"> Tutorials </span>===
+
=== <span style="color:#0000FF"> CDM Tutorials </span>===
  
=== <span style="color:#0000FF"> Contact people </span>===
+
=== <span style="color:#0000FF"> CDM Contact people </span>===
  
=== <span style="color:#0000FF"> Akcnowledgements </span>===
+
=== <span style="color:#0000FF"> CDM Akcnowledgements </span>===

Latest revision as of 17:05, 30 July 2013

Contents

Computational Fluid Dynamics module (CFDM)

CFDM Introduction

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

Shape.jpg Brief description of wat the model means, eventually insert link to the benchmark section...or whatever...

ADVERTISMENT STYLE no numerical details!!!

Shape.jpg Brief description of what the model means
Shape.jpg Brief description of what the model means
Shape.jpg Brief description of what the model means



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

In this application the Navier Stokes equations are solved bla bla..

CFDM Structure

CFDM Fluid types

  • Incompressible fluid

aaa

  • Compressible fluid

Constitutive laws

  • Newtonian
  • No-Newtonian
    • Bingham plastics
    • Variable yield model

CFDM Kinematical approaches

  • Eulerian
    • With free surface (level set)
    • Without free surface
  • Lagrangian PFEM (implicitly with free surface)

CFDM Solution strategy

  • Fractional step
  • Monolithic

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


In both cases a Newton Raphson residual based strategy is used for linearizing the problem.

CFDM Elements

Linear triangular elements in 2D and linear tetrahedra elements in 3D.

Stabilization techniques availables:

  • ASGS
  • OSS

CFDM Boundary conditions

  • Velocity boundary condition: Inlet of water
  • Pressure boundary condition: Pressure can be imposed strongly or weakly...
  • Wall boundary condition:
    • Slip/no slip boundary condition
    • Wall law
  • Flag variable?????

CFDM Initial conditions

Initial condition both in velocity and pressure can be set.

CFDM Turbulence models

The user can chose wether to use or not a turbulence model. Those available in kratos are:

  • Smagorinsky-Lily
  • Spalart-Allmaras

CFDM HPC

The code can be run in shared or distributed memory:

  • OpenMP:
  • MPI:

CFDM Problem parameters

CFDM Others relevand aspects

CFDM Benchmarking

CFDM Tutorials

CFDM Contact people

CFDM Akcnowledgements

Computational Structural Mechanics module (CSMM)

CSMM 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....

CSMM Application Structure

CSMM 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!!!!)


CSMM 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

CSMM Boundary Conditions

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

CSMM Loads

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

CSMM Constitutive laws

The following constitutive laws are available:

  • Linear elastic:
  • ...

CSMM HPC

The code can be run in shared or distributed memory:

  • OpenMP:
  • MPI:

CSMM Problem parameters

...

CSMM Others relevand aspects

...

CSMM Benchmarking

Here validation and verification examples should be inserted

CSMM Tutorials

CSMM Contact people

CSMM Akcnowledgements

Convection Diffusion module (CDM)

CDM 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}

CDM Structure

Analysis type

The available solution strategy is:

   Dynamic

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

CDM Kinematical approaches

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

CDM Solution strategies

CDM Elements

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

ConvDiff2D

ConvDiff3D

CDM 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)

CDM Initial conditions

Initial condition in temperature can be set.

CDM HPC

The code can be run in shared or distributed memory:

  • OpenMP:
  • MPI:

CDM 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

CDM Others relevand aspects

CDM Benchmarking

CDM Tutorials

CDM Contact people

CDM Akcnowledgements

Personal tools
Categories