1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
|
#ifndef VISCOUS_STRESS_H
#define VISCOUS_STRESS_H
#include <map>
#include <string>
#include <vector>
#include <fstream>
#include "MatrixLibrary.h"
#include "ATC_TypeDefs.h"
namespace ATC {
/**
* @class ViscousStress
* @brief Base class that defines interface for a constitutive law
* @brief that computes viscous stresses given all field and gradient information.
*/
class ViscousStress
{
public:
ViscousStress() {};
virtual ~ViscousStress() {};
virtual void initialize(void){};
//* Returns parameter values, (Nothing uses this).
virtual void parameters(std::map<std::string,double> & /* parameters */) {}
//* Computes viscous stress given a strain rate tensor.
//* Units: mvv/L^3 (i.e. for units Real: g/(mol ps^2 A^2) )
virtual void viscous_stress(const FIELD_MATS &fields,
const GRAD_FIELD_MATS &gradFields,
DENS_MAT_VEC &stress)=0;
virtual void viscosity(const FIELD_MATS & /* fields */,
DENS_MAT & /* coefs */) const
{throw ATC_Error("ViscousStress::viscosity: unimplemented function");}
//* Returns the derivative of the stress tensor for a given strain-rate tensor.
virtual void tangent(const MATRIX & /* F */, MATRIX & /* C */) const
{throw ATC_Error("ViscousStress::tangent: unimplemented function");}
};
/**
* @class ViscousStressConstant
* @brief Class for computing stress for a constant viscosity material
* assuming divergence-free flow
*/
class ViscousStressConstant : public ViscousStress
{
public:
ViscousStressConstant():viscosity_(0.){};
ViscousStressConstant(std::fstream &matfile);
ViscousStressConstant(double viscosity)
: viscosity_(viscosity) {};
void viscous_stress(const FIELD_MATS &fields,
const GRAD_FIELD_MATS &gradFields,
DENS_MAT_VEC &flux);
void viscosity(const FIELD_MATS & fields,
DENS_MAT & coefs) const;
protected:
double viscosity_;
void set_tangent();
};
}
#endif
|
