7 OutPut(
"Example: Poiseuille_CST_Oldroyd.h" << endl) ;
13 void ExactU1(
double x,
double y,
double *values)
15 values[0] = 4*y*(1-y);
21 void ExactU2(
double x,
double y,
double *values)
29 void ExactP(
double x,
double y,
double *values)
40 void BoundCondition(
int i,
double t, BoundCond &cond)
45 void U1BoundValue(
int BdComp,
double Param,
double &value)
51 case 1: value=4*Param*(1-Param);
55 case 3: value=4*Param*(1-Param);
57 default: cout <<
"wrong boundary part number" << endl;
62 void U2BoundValue(
int BdComp,
double Param,
double &value)
74 default: cout <<
"wrong boundary part number" << endl;
82 void LinCoeffs(
int n_points,
double *x,
double *y,
83 double **parameters,
double **coeffs)
89 for(i=0;i<n_points;i++)
102 void ExactS1(
double x,
double y,
double *values)
110 void ExactS2(
double x,
double y,
double *values)
119 void ExactS3(
double x,
double y,
double *values)
121 values[0] = (x*y) + 5.0;
130 void BoundCondition_CST(
int i,
double t, BoundCond &cond)
136 void S1BoundValue(
int BdComp,
double Param,
double &value)
140 case 0: value=Param + 3.0;
142 case 1: value=1 + 3.0;
144 case 2: value=(1-Param) + 3.0;
146 case 3: value=0 + 3.0;
148 default: cout <<
"wrong boundary part number" << endl;
156 void S2BoundValue(
int BdComp,
double Param,
double &value)
166 case 3: value=(1-Param);
168 default: cout <<
"wrong boundary part number" << endl;
177 void S3BoundValue(
int BdComp,
double Param,
double &value)
181 case 0: value=0 + 5.0;
183 case 1: value=Param + 5.0;
185 case 2: value=(1-Param) + 5.0;
187 case 3: value=0 + 5.0;
189 default: cout <<
"wrong boundary part number" << endl;
198 void LinCoeffs_CST(
int n_points,
double *X,
double *Y,
199 double **parameters,
double **coeffs)
203 double *coeff, *param, x, y;
204 double u1, u1x, u1y, u2, u2x, u2y, tau1, tau2, tau3, tau1x, tau1y, tau2x, tau2y, tau3x, tau3y;
206 for(i=0;i<n_points;i++)
232 coeff[1] = (u1*tau1x) + (u2*tau1y) + (nu*(tau1)) - (2*(tau1*u1x + tau2*u1y));
234 coeff[2] = (u1*tau2x) + (u2*tau2y) + (nu*tau2) - (tau1*u2x + tau3*u1y);
236 coeff[3] = (u1*tau3x) + (u2*tau3y) + (nu*(tau3)) - (2*(tau2*u2x + tau3*u2y));
double RE_NR
Definition: Database.h:313
static TParamDB * ParamDB
Definition: Database.h:1134