Simulation code for the Master's Thesis: Mathematical Modeling of Microplastics Spread in the Marine Environment
As the PDE system is stiff, the simulation code normalizes the system. For this we define the referenced parameters,
and
Inserting these values into the system, yields
This normalized system then gets solved. As the reference values, the peak concentration of the Gaussian is used.
To solve the complete system implicitly, FiPy requires the terms to be inside Implicit... Terms, otherwise
terms would be treated explicitly. When solving implicitly, the terms must also be added to the
complete equation and if the term is negative, the sign must be inside the Implicit... term to have correct behavior.