7 OutPut(
"Example: Poiseuille_Coupled.h" << endl) ;
13 void ExactU1(
double x,
double y,
double *values)
15 values[0] = 4*y*(1-y);
23 void ExactU2(
double x,
double y,
double *values)
31 void ExactP(
double x,
double y,
double *values)
42 void BoundCondition(
int i,
double t, BoundCond &cond)
47 void U1BoundValue(
int BdComp,
double Param,
double &value)
53 case 1: value=4*Param*(1-Param);
57 case 3: value=4*Param*(1-Param);
59 default: cout <<
"wrong boundary part number" << endl;
66 void U2BoundValue(
int BdComp,
double Param,
double &value)
78 default: cout <<
"wrong boundary part number" << endl;
86 void LinCoeffs(
int n_points,
double *x,
double *y,
87 double **parameters,
double **coeffs)
93 for(i=0;i<n_points;i++)
106 void ExactS1(
double x,
double y,
double *values)
116 void ExactS2(
double x,
double y,
double *values)
126 void ExactS3(
double x,
double y,
double *values)
128 values[0] = (x*y) + 5.0;
137 void BoundCondition_CST(
int i,
double t, BoundCond &cond)
143 void S1BoundValue(
int BdComp,
double Param,
double &value)
147 case 0: value=Param + 3.0;
149 case 1: value=1 + 3.0;
151 case 2: value=(1-Param) + 3.0;
153 case 3: value=0 + 3.0;
155 default: cout <<
"wrong boundary part number" << endl;
163 void S2BoundValue(
int BdComp,
double Param,
double &value)
173 case 3: value=(1-Param);
175 default: cout <<
"wrong boundary part number" << endl;
184 void S3BoundValue(
int BdComp,
double Param,
double &value)
188 case 0: value=0 + 5.0;
190 case 1: value=Param + 5.0;
192 case 2: value=(1-Param) + 5.0;
194 case 3: value=0 + 5.0;
196 default: cout <<
"wrong boundary part number" << endl;
205 void LinCoeffs_CST(
int n_points,
double *X,
double *Y,
206 double **parameters,
double **coeffs)
210 double *coeff, *param, x, y;
211 double u1, u1x, u1y, u2, u2x, u2y, tau1, tau2, tau3, tau1x, tau1y, tau2x, tau2y, tau3x, tau3y;
213 for(i=0;i<n_points;i++)
219 param = parameters[i];
239 coeff[1] = (u1*tau1x) + (u2*tau1y) + ((1-(2.0*alpha))*nu*(tau1)) - (2*(tau1*u1x + tau2*u2x)) + (alpha*nu*((tau1*tau1) + (tau2*tau2)));
241 coeff[2] = (u1*tau2x) + (u2*tau2y) + ((1-(2.0*alpha))*nu*(tau2)) - (tau1*u1y + tau3*u2x) + (alpha*nu*((tau1*tau2) + (tau2*tau3)));
243 coeff[3] = (u1*tau3x) + (u2*tau3y) + ((1-(2.0*alpha))*nu*(tau3)) - (2*(tau2*u1y + tau3*u2y)) + (alpha*nu*((tau2*tau2) + (tau3*tau3)));
double RE_NR
Definition: Database.h:313
static TParamDB * ParamDB
Definition: Database.h:1134