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

<b><i>* Cahn-Hilliard simulations *</i></b><br>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.<br>-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-&gt; C system with initial noise condition.Entering:<i>python3 CH_subcritical_2D_numba_10.py 100 </i><i>--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 </i><i>--source_noise run_5101_noise_source.npz -IC run_5101_sRest.hdf5</i>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.<br><b>* Brusselator Simulations *</b>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. <br>