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)