One two d amr

Allwmake

@@ -12,6 +12,8 @@ wmake src/dfChemistryModel
 wmake src/TurbulenceModels/compressible
 wmake src/regionModels/surfaceFilmModels
 src/lagrangian/Allwmake $targetType $*
+wmake src/dynamicMesh
+wmake src/dynamicFvMesh
 
 wmake applications/solvers/df0DFoam
 wmake applications/solvers/dfLowMachFoam

applications/solvers/dfHighSpeedFoam/dfHighSpeedFoam.C

@@ -60,6 +60,7 @@ int main(int argc, char *argv[])
     double time_monitor_flow;
     double time_monitor_chem;
     double time_monitor_Y;
+    double time_monitor_AMR;
     clock_t start, end;
 
     turbulence->validate();
@@ -69,6 +70,7 @@ int main(int argc, char *argv[])
     #include "readFluxScheme.H"
 
     dimensionedScalar v_zero("v_zero", dimVolume/dimTime, 0.0);
+    dimensionedScalar refCri("refCri", dimensionSet(1, -4, 0, 0, 0), 0.0);
 
     // Courant numbers used to adjust the time-step
     scalar CoNum = 0.0;
@@ -80,13 +82,28 @@ int main(int argc, char *argv[])
     {
         #include "readTimeControls.H"
 
+        //used for AMR
+        refCri = max(mag(fvc::grad(rho)));
+        tmp<volScalarField> tmagGradrho = mag(fvc::grad(rho));
+        volScalarField normalisedGradrho
+        (
+            "normalisedGradrho",
+            tmagGradrho()/refCri
+        );
+        normalisedGradrho.writeOpt() = IOobject::AUTO_WRITE;
+        tmagGradrho.clear();
+
         if (!LTS)
         {
             #include "setDeltaT.H"
             runTime++;
 
             // Do any mesh changes
+            start = std::clock();
             mesh.update();
+            end = std::clock();
+            time_monitor_AMR += double(end - start) / double(CLOCKS_PER_SEC);
+
         }
 
         // --- Directed interpolation of primitive fields onto faces
@@ -225,6 +242,7 @@ int main(int argc, char *argv[])
         Info<< "MonitorTime_chem = " << time_monitor_chem << " s" << nl << endl;
         Info<< "MonitorTime_Y = " << time_monitor_Y << " s" << nl << endl;
         Info<< "MonitorTime_flow = " << time_monitor_flow << " s" << nl << endl;
+        Info<< "MonitorTime_AMR = " << time_monitor_AMR << " s" << nl << endl;
 
         Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
             << "  ClockTime = " << runTime.elapsedClockTime() << " s"

examples/dfHighSpeedFoam/det_1D_H2_01mm/constant/CanteraTorchProperties → examples/dfHighSpeedFoam/oneD_detonationH2/constant/CanteraTorchProperties

@@ -21,8 +21,8 @@ CanteraMechanismFile "H2_Li.xml";
 transportModel "UnityLewis";
 loadbalancing
 {
-        active  false;
-        //log   true;
+        active  true;
+        log   true;
 }
 odeCoeffs
 {}
@@ -32,9 +32,4 @@ zeroDReactor
 {
     //cvorcp "UV";
 }
-loadbalancing
-{
-        active  true;
-        //log   true;
-}
 // ************************************************************************* //

examples/dfHighSpeedFoam/oneD_detonationH2/constant/dynamicMeshDict

@@ -0,0 +1,68 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  7
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      dynamicMeshDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dynamicFvMesh   dfDynamicRefineFvMesh;
+
+multiCritRefinementControls
+{
+    enableMultiCritRefinementControl  false;
+}
+
+dfDynamicRefineFvMeshCoeffs
+{
+    // How often to refine
+    refineInterval  1; 
+
+    // Field to be refinement on
+    field           normalisedGradrho;
+
+    // Refine field in between lower..upper
+    lowerRefineLevel 0.2;
+    upperRefineLevel 1;
+
+    // If value < unrefineLevel unrefine
+    unrefineLevel   0.2;
+
+    // Have slower than 2:1 refinement
+    nBufferLayers   3;
+
+    // Refine cells only up to maxRefinement levels
+    maxRefinement   3;
+
+    // Stop refinement if maxCells reached
+    maxCells        1000000;
+
+    // Flux field and corresponding velocity field. Fluxes on changed
+    // faces get recalculated by interpolating the velocity. Use 'none'
+    // on surfaceScalarFields that do not need to be reinterpolated.
+    correctFluxes
+    (
+        (phi none)
+        //(nHatf none) 
+        //(rhoPhi none)
+        //(alphaPhi0.water none)
+        //(ghf none)
+        (neg none)
+        (pos none)
+    );
+
+    // Write the refinement level as a volScalarField
+    dumpLevel       true;
+}
+
+
+// ************************************************************************* //

examples/dfHighSpeedFoam/det_1D_H2_01mm/constant/thermophysicalProperties → examples/dfHighSpeedFoam/oneD_detonationH2/constant/thermophysicalProperties

@@ -14,5 +14,5 @@ FoamFile
     object      thermophysicalProperties;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-inviscid      true;
+inviscid      false;
 // ************************************************************************* //

examples/dfHighSpeedFoam/det_1D_H2_01mm/system/blockMeshDict → examples/dfHighSpeedFoam/oneD_detonationH2/system/blockMeshDict

@@ -30,7 +30,7 @@ vertices
 
 blocks
 (
-    hex (0 1 2 3 4 5 6 7) (5000 1 1) simpleGrading (1 1 1)
+    hex (0 1 2 3 4 5 6 7) (625 1 1) simpleGrading (1 1 1)
 );
 
 edges

examples/dfHighSpeedFoam/det_1D_H2_01mm/system/controlDict → examples/dfHighSpeedFoam/oneD_detonationH2/system/controlDict

@@ -14,7 +14,9 @@ FoamFile
     object      controlDict;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-libs ("libfieldFunctionObjects.so");
+libs (
+      "libdfDynamicFvMesh.so"
+      );
 
 application     dfHighSpeedFoam;
 
@@ -24,7 +26,7 @@ startTime       0;
 
 stopAt          endTime;
 
-endTime         2.4e-4;
+endTime         2.2e-4;
 
 deltaT          1e-09;
 
@@ -48,7 +50,7 @@ runTimeModifiable true;
 
 adjustTimeStep  yes;
 
-maxCo           0.02;
+maxCo           0.1;
 
 maxDeltaT       1;
 
@@ -57,8 +59,8 @@ functions
     minMax
     {
       type          fieldMinMax;
-      libs          ("libfieldFunctionObjects.so");
-      writeControl  timeStep; //writeTime;
+      libs          ("libdfFieldFunctionObjects.so");
+      writeControl  writeTime; //timeStep;
       fields        (p);
     }
 }

examples/dfHighSpeedFoam/oneD_detonationH2/system/decomposeParDict

@@ -0,0 +1,32 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  7
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      decomposeParDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+numberOfSubdomains  4;
+
+method          scotch;
+
+constraints
+{
+    dfRefinementHistory
+    {
+        //- Decompose cells such that all cell originating from single cell
+        //  end up on same processor
+        type    dfRefinementHistory;
+    }
+}
+
+// ************************************************************************* //

examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/AR

@@ -0,0 +1,42 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  7
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      AR;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0.8861; 
+
+
+boundaryField
+{
+    Left
+    {
+        type            zeroGradient;
+    }
+
+    Right
+    {
+        type            zeroGradient;
+    }
+
+    empty
+    {
+        type            empty;
+    }
+}
+
+
+// ************************************************************************* //

examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/H2

@@ -0,0 +1,41 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  7
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      H2;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0.01266;
+
+
+boundaryField
+{
+    Left
+    {
+        type            zeroGradient;
+    }
+
+    Right
+    {
+        type            zeroGradient;
+    }
+
+    empty
+    {
+        type            empty;
+    }
+}
+
+// ************************************************************************* //

examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/O2

@@ -0,0 +1,41 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  7
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      O2;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0.101266;
+
+
+boundaryField
+{
+    Left
+    {
+        type            zeroGradient;
+    }
+
+    Right
+    {
+        type            zeroGradient;
+    }
+
+    empty
+    {
+        type            empty;
+    }
+}
+
+// ************************************************************************* //