Modified Weak Galerkin of Heat equation in 2D

Dear all experts, Respected Francois Sir, and Respected Hecht sir, I have tried to implement modified weak Galerkin in the heat equation 2D FF++ code using the poisson equation modified weak Galerkin code . I am wondering why they both not giving good results sir. Can you please check my code sir??.

  1. Heat equation in 2D with u=0 in boundary. Code>
    Heat_WGM_Sine.edp (2.4 KB)
  2. Heat equation in 2D with neumann derivative of u=0 in boundary. Code>
    MWG_Heat_2D_Cosine.edp (2.3 KB)

Thanks in advance!.

You miss exp(t) in the last term of definition of f, and *(nTonEdge-1) in the two intalledges involving the mean
Heat_WGM_Sine.edp (2.4 KB)
For Neumann BC, do not put the factor (nTonEdge-1) in the two intalledges involving the mean, but put it in the intalledges involving the jumps.

1 Like

Thank you so much sir. If you get time kindly look the post of integro-partial differential equation code in FF++. I have stored the solution of previous steps and computed the approximate for integral but results not coming good sir.