Finites differences python code to solve CH equation with a source term and Comsol routine to solve Brusselator equation in radial domains.

* Cahn-Hilliard simulations *
Finite difference code implementing the modified Cahn Hilliard equation with a forward Euler scheme and the possibility to parallelize the solver using the numba python library.

-Additional files containing the source term (run_5101_noise_source) and the initial condition (run_5101_sRest.hdf5) are provided to reproduce the coupling with A+B-> C system with initial noise condition.
Entering:
python3 CH_subcritical_2D_numba_10.py 100 --sigma 0.1 -B -1.0 -Lx 60 -Nx 14 -dt 2e4 --CH_norm -cpu 3 --interpolator --source_file CH2470_PS_SO_source_hres_sigma_0.1_flatten_interpolator.npz --source_noise run_5101_noise_source.npz -IC run_5101_sRest.hdf5
In a bash terminal should be able to reproduce run in Fig 2a when storing all the provided file in the local folder.
The comsol routine CH2470_PS_SO allow to generate different source term adapting the size of the simulation domain an the duration of the simulation.

* Brusselator Simulations *
The comsol routine BR0140_PS_A_B is written to reproduce the simulations at DY=20 and variable DX. Additional parameters can be explored by simple modifications in the parameters menu of the same routine.