00001 /* 00002 volTensorField gradU = fvc::grad(U); 00003 volSymmTensorField D = symm(fvc::grad(U)); 00004 volTensorField Dprim = symm(fvc::grad(U - UMean)); 00005 00006 volScalarField prod = -((U - UMean)*(U - UMean)) && D; 00007 */ 00008 00009 /* 00010 volScalarField txx 00011 ( 00012 IOobject 00013 ( 00014 "txx", 00015 runTime.timeName(), 00016 mesh, 00017 IOobject::NO_READ, 00018 IOobject::NO_WRITE 00019 ), 00020 mesh, 00021 dimensionSet(0, 1, -1, 0, 0) 00022 ); 00023 txx =sqrt(Txx - (UMeanx*UMeanx)); 00024 txx.write(); 00025 00026 volScalarField tyy 00027 ( 00028 IOobject 00029 ( 00030 "tyy", 00031 runTime.timeName(), 00032 mesh, 00033 IOobject::NO_READ, 00034 IOobject::NO_WRITE 00035 ), 00036 mesh, 00037 dimensionSet(0, 1, -1, 0, 0) 00038 ); 00039 tyy = sqrt(Tyy - (UMeany*UMeany)); 00040 tyy.write(); 00041 00042 volScalarField tzz 00043 ( 00044 IOobject 00045 ( 00046 "tzz", 00047 runTime.timeName(), 00048 mesh, 00049 IOobject::NO_READ, 00050 IOobject::NO_WRITE 00051 ), 00052 mesh, 00053 dimensionSet(0, 1, -1, 0, 0) 00054 ); 00055 tzz = sqrt(Tzz - (UMeanz*UMeanz)); 00056 tzz.write(); 00057 00058 volScalarField txy 00059 ( 00060 IOobject 00061 ( 00062 "txy", 00063 runTime.timeName(), 00064 mesh, 00065 IOobject::NO_READ, 00066 IOobject::NO_WRITE 00067 ), 00068 mesh, 00069 dimensionSet(0, 2, -2, 0, 0) 00070 ); 00071 txy = Txy - (UMeanx*UMeany); 00072 txy.write(); 00073 */ 00074 00075 // ************************ vim: set sw=4 sts=4 et: ************************ //
