code: def _update_physics(self):
          
          phase = self.project.phases()[self.settings['phase']]
          physics = self.project.find_physics(phase=phase)
          for item in self.settings['sources']:
              # Regenerate models with new guess
              quantity = self.settings['quantity']
              # Put quantity on phase so physics finds it when regenerating
              phase[quantity] = self[quantity]
              # Regenerate models, on either phase or physics
              phase.regenerate_models(propnames=item)
              for phys in physics:
                  phys.regenerate_models(propnames=item)
docs: r
          Update physics using the current value of 'quantity'
  
          Notes
          -----
          The algorithm directly writes the value of 'quantity' into the phase.
          This method was implemented relaxing one of the OpenPNM rules of
          algorithms not being able to write into phases.