Calculated fvSource on the fly

friedenhe · Jun 23, 2024 · 7b3f157 · 7b3f157
1 parent 77be385
commit 7b3f157
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/src/adjoint/DASolver/DAHeatTransferFoam/DAHeatTransferFoam.C b/src/adjoint/DASolver/DAHeatTransferFoam/DAHeatTransferFoam.C
@@ -113,6 +113,12 @@ label DAHeatTransferFoam::solvePrimal(
             Info << "Time = " << runTime.timeName() << nl << endl;
         }
 
+        if (hasFvSource_)
+        {
+            volScalarField& fvSource = fvSourcePtr_();
+            daFvSourcePtr_->calcFvSource(fvSource);
+        }
+
         fvScalarMatrix TEqn(
             fvm::laplacian(k, T)
             + fvSource);