an efficient C++ finite element environment
|
|
Go to the documentation of this file.
30 int main(
int argc,
char**argv) {
35 Float We_incr = (argc > 2) ? atof(argv[2]) : 0.1;
36 Float We_max = (argc > 3) ? atof(argv[3]) : 0.1;
37 Float delta_t0 = (argc > 4) ? atof(argv[4]) : 0.005;
38 string restart = (argc > 5) ? argv[5] :
"";
44 Float delta_t_min = 1e-5;
45 Float We_incr_min = 1e-5;
48 branch even (
"We",
"tau",
"u",
"p");
50 pb.
initial (omega, tau_h, uh, ph, restart);
51 dout << even (pb.
We, tau_h, uh, ph);
54 if (ok) pb.
We += We_incr;
55 derr <<
"# We = " << pb.
We <<
" delta_t = " << pb.
delta_t << endl;
56 field tau_h0 = tau_h, uh0 = uh, ph0 = ph;
57 ok = pb.
solve (tau_h, uh, ph);
59 dout << even (pb.
We, tau_h, uh, ph);
62 tau_h = tau_h0; uh = uh0; ph = ph0;
64 derr <<
"# solve failed: decreases We_incr and retry..." << endl;
68 if (We_incr < We_incr_min)
break;
70 derr <<
"# solve failed: decreases delta_t and retry..." << endl;
The Oldroyd problem by the theta-scheme – class header.
see the catchmark page for the full documentation
see the field page for the full documentation
bool solve(field &tau_h, field &uh, field &ph)
rheolef - reference manual
void initial(const geo &omega, field &tau_h, field &uh, field &ph, string restart)
see the environment page for the full documentation
This file is part of Rheolef.
The Oldroyd problem on the driven cavity benchmark – boundary conditions.
odiststream derr(cerr)
see the diststream page for the full documentation
see the Float page for the full documentation
see the branch page for the full documentation
int main(int argc, char **argv)
odiststream dout(cout)
see the diststream page for the full documentation
see the geo page for the full documentation