| 36 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // |
| 37 | |
| 38 | int main() |
| 39 | { |
| 40 | scalarField eta(200); |
| 41 | scalarField C1mR(eta.size()); |
| 42 | |
| 43 | forAll(eta, i) |
| 44 | { |
| 45 | eta[i] = scalar(i)/10.0; |
| 46 | } |
| 47 | |
| 48 | scalar C1 = 1.42; |
| 49 | scalar eta0 = 4.38; |
| 50 | scalar beta = 0.012; |
| 51 | |
| 52 | C1mR = C1 - ((eta*(1.0 - eta/eta0))/(1.0 + beta*pow(eta, 3.0))); |
| 53 | |
| 54 | graph("C&1! - R", "C&1! - R", "@h!", eta, C1mR).write |
| 55 | ( |
| 56 | "C1mR", |
| 57 | "xmgr" |
| 58 | ); |
| 59 | |
| 60 | Info<< "end" << endl; |
| 61 | } |
| 62 | |
| 63 | |
| 64 | // ************************************************************************* // |