Carmen Code
 All Classes Files Functions Variables Macros Pages
Functions
Source.cpp File Reference

Computes the source terms of the system. More...

#include "Carmen.h"
Include dependency graph for Source.cpp:

Functions

Vector Source (Cell &UserCell)
 Returns the source term in the cell UserCell. More...
 

Detailed Description

Computes the source terms of the system.

Function Documentation

Vector Source ( Cell UserCell)

Returns the source term in the cell UserCell.

Parameters
UserCellCell value
Returns
Vector

Gravity vector

24 {
25  // --- Local variables ---
26 
27  Vector Force(Dimension);
28  Vector Result(QuantityNb);
29  Result.setZero();
30 
32  Vector V(3);
33  real Gx=0., Gy=0., Gz=0.,rho=0.;
34  for(int i=1;i<=3;i++)
35  V.setValue(i,UserCell.average(i+1));
36  rho = UserCell.density();
37  Gz = 0.2;
38  Result.setValue(2,rho*Gx);
39  Result.setValue(3,rho*Gy);
40  Result.setValue(4,rho*Gz);
41  Result.setValue(5,rho*(Gx*V.value(1) + Gy*V.value(2) + Gz*V.value(3)));
42 
43  Result.setZero();
44  return Result;
45 
46 }
int QuantityNb
Definition: Parameters.cpp:171
Standard class for every vector in Carmen.
Definition: Vector.h:29
real density() const
Returns the cell-average density.
Definition: Cell.h:1262
int Dimension
Definition: Parameters.cpp:74
real average(const int QuantityNo) const
Returns the component no. QuantityNo of the cell-average value.
Definition: Cell.h:1128
#define real
Definition: PreProcessor.h:31

Here is the caller graph for this function: