Initial & Boundary conditions

[vincent-reau]

Hi everyone, I'm having some issues with the boundary and initial conditions. I have 2 parameters Cs and Cg.

• x ∈ [0, 1]
• t ∈ [0, 1]

Initial conditions:

• Cs = 0, at x = 0 and t = 0
• Cg = 0, at x = 0 and t = 0

Boundary conditions:

• Cg = 1, at x = 0 and t !=0
• dCg/dx = 0, at x=L and t !=0

Here's my code. Is it well transcript or do i need to make some changes ?

def cs_func(x):
  return 0.0

def cg_func(x):
  return 0.0

ic_cs = dde.icbc.IC(geomtime, cs_func, lambda _, on_initial: on_initial, component=0)
ic_cg = dde.icbc.IC(geomtime, cg_func, lambda _, on_initial: on_initial, component=1)

def boundary_l(x,on_boundary):
  return on_boundary  and np.isclose(x[0], 0)  # définit le boundary x=0, left

def boundary_r(x,on_boundary):
  return on_boundary  and np.isclose(x[0], 1)  # définit le boundary x=L, right

bc_cg = dde.icbc.DirichletBC(
    geomtime,
    lambda x: 1,
    boundary_l,
    component = 1 
) # C_g(0,t)=1

bc_dcg = dde.icbc.OperatorBC(
    geomtime,
    lambda x, y, _ : dde.grad.jacobian(y, x, i=0, j=1), 
    boundary_r
) # dc_g/dz(L,t)=0

[lululxvi]

Looks OK.