Add coded function object to calculate mass flux

Author: staneuski

etc/caseDicts/functions/fields/compressibilityMachNo

@@ -2,27 +2,37 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     | Website:  https://openfoam.org
-    \\  /    A nd           | Version:  dev
+    \\  /    A nd           | Version:  dev (20241102)
      \\/     M anipulation  |
+-------------------------------------------------------------------------------
+Description
+    Caches and writes Mach number based on compressibility
+
 \*---------------------------------------------------------------------------*/
 
 type            coded;
+libs            ( "libutilityFunctionObjects.so" );
 
-libs            ("libutilityFunctionObjects.so");
+executeControl  writeTime;
+writeControl    writeTime;
 
-name            compressibilityMachNo;
+codeFields
+#{
+    fields.append("psi");
+    fields.append("U");
+#};
 
-writeControl    writeTime;
+codeExecute
+#{
+    const volScalarField& psi = mesh().lookupObject<volScalarField>("psi");
+    const volVectorField& U = mesh().lookupObject<volVectorField>("UMean");
+
+    store("Ma", mag(U)*sqrt(psi));
+#};
 
 codeWrite
 #{
-    const volVectorField& U(mesh().lookupObject<volVectorField>("U"));
-    const volScalarField& psi(mesh().lookupObject<volScalarField>("psi"));
-
-    volScalarField::New(
-        "Ma",
-        mag(U) * sqrt(psi)
-    )->write();
+    return writeObject("Ma");
 #};
 
 // ************************************************************************* //

etc/caseDicts/functions/fields/massFlux

@@ -0,0 +1,43 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev (20241102)
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+Description
+    Caches and writes mass flux
+
+\*---------------------------------------------------------------------------*/
+
+type            coded;
+libs            ( "libutilityFunctionObjects.so" );
+
+executeControl  writeTime;
+writeControl    writeTime;
+
+codeInclude
+#{
+    #include "fvcFlux.H"
+#};
+
+codeFields
+#{
+    fields.append("rho");
+    fields.append("U");
+#};
+
+codeExecute
+#{
+    const volScalarField& rho = mesh().lookupObject<volScalarField>("rho");
+    const volVectorField& U = mesh().lookupObject<volVectorField>("U");
+
+    store("varphi", linearInterpolate(rho*U) & mesh().Sf());
+#};
+
+codeWrite
+#{
+    return writeObject("varphi");
+#};
+
+// ************************************************************************* //