Minor cleanup in Nutils participants of two-scale-heat-conduction tutorial

Author: IshaanDesai

two-scale-heat-conduction/macro-nutils/macro.py

@@ -15,8 +15,6 @@ def main():
     """
     is_coupled_case = True  # If False, single-physics problem is solved
 
-    # Original case from Bastidas et al.
-    # topo, geom = mesh.rectilinear([np.linspace(0, 1.0, 9), np.linspace(0, 0.5, 5)])
     topo, geom = mesh.rectilinear([np.linspace(0, 1.0, 5), np.linspace(0, 0.5, 4)])
 
     ns = function.Namespace(fallback_length=2)

two-scale-heat-conduction/micro-nutils/micro.py

@@ -20,7 +20,6 @@ def __init__(self, sim_id):
         self._sim_id = sim_id
 
         # Initial parameters
-        # self._nelems = 10  # Elements in one direction (original case from Bastidas et al.)
         self._nelems = 6  # Elements in one direction
 
         self._ref_level = 3  # Number of levels of mesh refinement
@@ -72,17 +71,6 @@ def initialize(self):
         # Initialize phase field once more on refined topology
         solphi = self._get_analytical_phasefield(self._topo, self._ns, self._degree_phi, self._ns.lam, self._r_initial)
 
-        target_porosity = 1 - math.pi * self._r_initial ** 2
-        print("Target amount of void space = {}".format(target_porosity))
-
-        # Solve Allen-Cahn equation till we reach target porosity value
-        # dt_initial = 1E-3
-        # psi = 0
-        # while psi < target_porosity:
-        #     print("Solving Allen-Cahn equation to achieve initial target grain structure")
-        #     solphi = self._solve_allen_cahn(self._topo, solphi, 0.5, dt_initial)
-        #     psi = self._get_avg_porosity(self._topo, solphi)
-
         self._solphi = solphi  # Save solution of phi
         psi = self._get_avg_porosity(self._topo, solphi)
         self._psi_nm1 = psi  # Average porosity value of last time step