Forum Message

 

 

We have moved the forum to https://forum.ngsolve.org . This is an archived version of the topics until 05/05/23. All the topics were moved to the new forum and conversations can be continued there. This forum is just kept as legacy to not invalidate old links. If you want to continue a conversation just look for the topic in the new forum.

Notice

The forum is in read only mode.

How to handle Implicit Euler method when replacing domain and mesh

More
3 years 9 months ago #3016 by dong
How to handle Implicit Euler method when replacing domain and mesh was created by dong
I tried to implement the implicit Euler method for a time dependent problem.
[tex]\partial_t c+\nabla\cdot(cu-k\nabla c)=f[/tex].
Please see the attached file.

The code run well with the unit square domain.
Code:
mesh = Mesh(unit_square.GenerateMesh(maxh=0.5) )

When I replaced the unit square on line 32 (above) by a new domain (a unit square with two subdomains) by line 33.
Code:
mesh = Mesh(geo.GenerateMesh(maxh=0.5))
and line 83
Code:
a += IfPos(u*n,1,0)*B_ds_out*ds(skeleton=True,definedon=mesh.Boundaries("right|top|left|bottom"))
by line 84
Code:
a += IfPos(u*n,1,0)*B_ds_out*ds(skeleton=True,definedon=mesh.Boundaries("GammaS|GammaD"))
I think this would give a similar results. But the new one doesn't make any sense.

Could you please tell me what the mistakes I made? Do I need to change anything in the bilinear and linear forms? Any help would be appreciated.
Thank you so much.
Attachments:
Time to create page: 0.112 seconds

Login