Fixed the wallHeatFlux BC and added xyz output for deformedFFDs (#834)

* Made writeDeformedFFDs also write xyz files.

* Fixed an issue for the fixedWallHeatFlux BC.

Author: friedenhe

dafoam/mphys/mphys_dafoam.py

@@ -480,6 +480,7 @@ def solve_linear(self, d_outputs, d_residuals, mode):
                                 )
                             else:
                                 self.DVGeo.writeTecplot("deformedFFDs_%d.dat" % self.solution_counter)
+                                self.DVGeo.writePlot3d("deformedFFDs_%d.xyz" % self.solution_counter)
 
                         # write the deformed constraints for post-processing
                         if DASolver.getOption("writeDeformedConstraints"):

src/adjoint/DAMisc/fixedWallHeatFlux/fixedWallHeatFluxFvPatchScalarField.C

@@ -111,6 +111,11 @@ void fixedWallHeatFluxFvPatchScalarField::write(Ostream& os) const
 {
     fvPatchScalarField::write(os);
     writeEntry("value", os);
+    os.writeEntry("Pr", Pr_);
+    os.writeEntry("Prt", Prt_);
+    os.writeEntry("nu", nu_);
+    os.writeEntry("Cp", Cp_);
+    os.writeEntry("heatFlux", heatFlux_);
 }
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //