Hi!

Q_T = [0,L] x [0,H] x [0,T]

1- We have the following variational formulation:

\int_{Q_T} dV/dt \psi dx dy dt = - D \int_{Q_T} \nabla V . \nabla \psi dx dy dt + D \int_0^T \int_0^L \psi(x,0,t) (\gamma VC - h Vf) dx dt - u0 \int_{Q_T} y \dV/dx \psi dx dy dt - \sigma_2 \int_{Q_T} V \psi dx dy dt

is it’s following formulation with FF++ is correct

problem system2(V,VH,solver=UMFPACK)=
int2d(Th)(V*VH/dt)+int2d(Th)(D*(dx(V)*dx(VH)+dy(V)*dy(VH)))
+int2d(Th)(u0*y*dx(V)*VH)+int2d(Th)(sigma2*V*VH)
-int2d(Th)(oldV*VH/dt)
+int1d(Th,1)(D*h*Vf*VH)-int1d(Th,1)(D*gamma*VC*VH);


​2- what is the mathematical problem corresponding to the FF++ formulation:

problem extend(Vk,Vt) = int2d(Th)(dy(Vk)*dy(Vt)) + on(1,Vk=V);

Regards

mesh Th = square(N,N,[L*x,H*y]);