Merge pull request #560 from circlexiang/master

Add 2D Riemann problem case

Author: pkuLmq

README.md

@@ -27,6 +27,13 @@ The neural network models used in the tutorial examples can be found at– [AIS
 Detailed guide for installation and tutorials is available on [our documentation website](https://deepflame.deepmodeling.com).
 
 ## Features
+New in v1.6 (2025/5/30):
+- Add a new solver, dfSteadyFoam, a steady-state compressible flow solver. It supports turbulence and uses the flexible PIMPLE algorithm to efficiently compute steady solutions
+- Extend energy model to support sensible enthalpy (hs) in addition to absolute enthalpy (ha) and internal energy (ea). Temperature can now be computed independently of Cantera via in-code Newton iteration, improving flexibility and removing external dependencies
+- Improve install.sh by displaying system architecture (e.g., x86_64) after build for clearer platform information
+- Switch to sphinx-book-theme and update documentation build environment (Ubuntu 24.04, Python 3.12) for better appearance and compatibility
+- Add 2D Riemann problem example case
+
 New in v1.5 (2025/1/8):
 - Provide a new boundary condition, totalFlowRateAdvectiveDiffusion (adopted from OpenFOAM v7), which accounts for diffusion effects at the boundary
 - Develop a new solver, `dfBuoyancyFoam`(adopted from fireFoam in OpenFOAM v7), a transient, compressible solver designed to model turbulent reacting flows incorporating buoyancy effects

examples/dfHighSpeedFoam/twoD_RiemannProblem/0_orig/N2

@@ -0,0 +1,46 @@
+/*--------------------------------*- 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      N2;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0.79;
+
+boundaryField
+{
+	inlet
+	{
+	    type            zeroGradient;
+	}
+	top
+	{
+		type           	zeroGradient;
+	}
+	bottom
+    {
+		type           	zeroGradient;
+    }
+	outlet
+	{
+	    type            zeroGradient;
+	}
+	frontAndBack
+	{
+	    type            empty;
+	}
+}
+
+// ************************************************************************* //

examples/dfHighSpeedFoam/twoD_RiemannProblem/0_orig/O2

@@ -0,0 +1,46 @@
+/*--------------------------------*- 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.21;
+
+boundaryField
+{
+	inlet
+	{
+	    type            zeroGradient;
+	}
+	top
+	{
+		type           	zeroGradient;
+	}
+	bottom
+    {
+		type           	zeroGradient;
+    }
+	outlet
+	{
+	    type            zeroGradient;
+	}
+	frontAndBack
+	{
+	    type            empty;
+	}
+}
+
+// ************************************************************************* //

examples/dfHighSpeedFoam/twoD_RiemannProblem/0_orig/T

@@ -0,0 +1,46 @@
+/*--------------------------------*- 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      T;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 1 0 0 0];
+
+internalField   uniform 348.4; // to be overwritten 
+
+boundaryField
+{
+    inlet
+    {
+        type            zeroGradient;
+    }
+    outlet
+    {
+        type            zeroGradient;
+    }
+    bottom
+    {
+        type           	zeroGradient;
+    }
+    top
+    {
+        type           	zeroGradient;
+    }
+    frontAndBack
+    {
+        type            empty;
+    }
+}
+
+// ************************************************************************* //

examples/dfHighSpeedFoam/twoD_RiemannProblem/0_orig/U

@@ -0,0 +1,46 @@
+/*--------------------------------*- 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       volVectorField;
+    location    "0";
+    object      U;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+	inlet
+	{
+	    type            zeroGradient;
+	}
+	top
+	{
+		type           	zeroGradient;
+	}
+	bottom
+    {
+		type           	zeroGradient;
+    }
+	outlet
+	{
+	    type            zeroGradient;
+	}
+	frontAndBack
+	{
+	    type            empty;
+	}
+}
+
+// ************************************************************************* //

examples/dfHighSpeedFoam/twoD_RiemannProblem/0_orig/Ydefault

@@ -0,0 +1,46 @@
+/*--------------------------------*- 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      Ydefault;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0.0;
+
+boundaryField
+{
+	inlet
+	{
+	    type            zeroGradient;
+	}
+	top
+	{
+		type           	zeroGradient;
+	}
+	bottom
+    {
+		type           	zeroGradient;
+    }
+	outlet
+	{
+	    type            zeroGradient;
+	}
+	frontAndBack
+	{
+	    type            empty;
+	}
+}
+
+// ************************************************************************* //

examples/dfHighSpeedFoam/twoD_RiemannProblem/0_orig/p

@@ -0,0 +1,46 @@
+/*--------------------------------*- 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      p;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -2 0 0 0 0];
+
+internalField   uniform 150000; // to be overwritten 
+
+boundaryField
+{
+    inlet
+    {
+        type            zeroGradient;
+    }
+    outlet
+    {
+        type            zeroGradient;
+    }
+    bottom
+    {
+        type           	zeroGradient;
+    }
+    top
+    {
+        type           	zeroGradient;
+    }
+    frontAndBack
+    {
+        type            empty;
+    }
+}
+
+// ************************************************************************* //

examples/dfHighSpeedFoam/twoD_RiemannProblem/Allclean

@@ -0,0 +1,11 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # Run from this directory
+
+echo "Cleaning log.*"
+rm log.*
+echo "Cleaning processor*"
+rm -r processor*
+echo "Cleaning polyMesh/"
+rm -r constant/polyMesh
+echo "Cleaning 0/"
+rm -r 0

examples/dfHighSpeedFoam/twoD_RiemannProblem/Allrun

@@ -0,0 +1,13 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # Run from this directory
+
+# Source tutorial run functions
+. $WM_PROJECT_DIR/bin/tools/RunFunctions
+
+application=dfHighSpeedFoam
+
+cp -r 0_orig/ 0/
+runApplication blockMesh
+runApplication setFields
+runApplication decomposePar
+runApplication mpirun -np 50 --allow-run-as-root $application -parallel