def setup_SimInfo() -> ch.SimInfo:
info = ch.SimInfo()
info.t0 = 0.0
info.nx = 150
info.ny = 150
info.dx = 1. / info.nx
info.dy = 1. / info.ny
info.bc = 'neumann'
return info
# ***************************************************************
# ********* POLYMER PARAMETERS *********
Xmin = 0.055
Xmax = 0.5
N = np.mean([200, 2000])
L_repeat = (10**-9) * np.mean([20, 80]) # meters
n_repeat = 15
L_omega = n_repeat * L_repeat
L_kuhn = (10**-9) * np.mean([0.5, 3.0]) # meters
# **************************************
# *********** INPUTS ***********
# Setup simulation info object
info = ch.SimInfo()
info.t0 = 0.0
info.nx = 128
info.ny = 128
info.dx = 1. / info.nx
info.dy = 1. / info.ny
info.bc = 'neumann'
info.rhs_type = 'ch_thermal_no_diffusion'
# Set up grid for spatial-field quantities
nx = int(info.nx)
xx, yy = np.meshgrid(np.arange(0, 1, 1 / info.nx),
np.arange(0, 1, 1 / info.nx))
# Setup CH parameter info object
chparams = ch.CHparamsVector(info.nx, info.ny)