Allowed setting viscosity in primalBC (#249)

* Added the capability to change mu and nu in primalBC.

* Fixed a bug.

Author: friedenhe

src/adjoint/DAField/DAField.C

@@ -598,6 +598,64 @@ void DAField::setPrimalBoundaryConditions(const label printInfo)
             const IOMRFZoneListDF& MRF = db.lookupObject<IOMRFZoneListDF>("MRFProperties");
             scalar& omega = const_cast<scalar&>(MRF.getOmegaRef());
             omega = omegaNew;
+            
+            if (printInfo)
+            {
+                Info << "Setting MRF omega to " << omegaNew << endl;
+            }
+
+            continue;
+        }
+        else if (bcKey == "transport:nu")
+        {
+            // change the nu field
+            scalar nu = bcDict.getScalar("transport:nu");
+            volScalarField& nuField = const_cast<volScalarField&>(
+                db.lookupObject<volScalarField>("nu"));
+            forAll(nuField, cellI)
+            {
+                nuField[cellI] = nu;
+            }
+            forAll(nuField.boundaryField(), patchI)
+            {
+                forAll(nuField.boundaryField()[patchI], faceI)
+                {
+                    nuField.boundaryFieldRef()[patchI][faceI] = nu;
+                }
+            }
+            nuField.correctBoundaryConditions();
+
+            if (printInfo)
+            {
+                Info << "Setting transportProperties nu to " << nu << endl;
+            }
+
+            continue;
+        }
+        else if (bcKey == "thermo:mu")
+        {
+            // change the nu field
+            scalar mu = bcDict.getScalar("thermo:mu");
+            volScalarField& muField = const_cast<volScalarField&>(
+                db.lookupObject<volScalarField>("thermo:mu"));
+            forAll(muField, cellI)
+            {
+                muField[cellI] = mu;
+            }
+            forAll(muField.boundaryField(), patchI)
+            {
+                forAll(muField.boundaryField()[patchI], faceI)
+                {
+                    muField.boundaryFieldRef()[patchI][faceI] = mu;
+                }
+            }
+            muField.correctBoundaryConditions();
+
+            if (printInfo)
+            {
+                Info << "Setting thermalphysicalProperties mu to " << mu << endl;
+            }
+
             continue;
         }
 

tests/runTests_DARhoSimpleCFoam.py

@@ -53,6 +53,7 @@
         "k0": {"variable": "k", "patches": ["inout"], "value": [k0]},
         "omega0": {"variable": "omega", "patches": ["inout"], "value": [omega0]},
         "useWallFunction": True,
+        "thermo:mu": 1.8e-5
     },
     "primalVarBounds": {
         "UMax": 1000.0,

tests/runTests_DASimpleFoam.py

@@ -51,6 +51,7 @@
         "k0": {"variable": "k", "patches": ["inout"], "value": [k0]},
         "omega0": {"variable": "omega", "patches": ["inout"], "value": [omega0]},
         "useWallFunction": False,
+        "transport:nu": 1.5e-5
     },
     "fvSource": {
         "disk1": {