- Thank you received: 0
LDG
- rhebergens
- Topic Author
- Offline
- Junior Member
Less
More
6 years 4 months ago #646
by rhebergens
LDG was created by rhebergens
Hello,
I was wondering what would be the best way to implement LDG for Stokes (see for example doi.org/10.1137/S0036142900380121 )
In LDG an auxiliary variable \sigma is introduced such that
\sigma = grad(u)
NGSolve has the VectorL2 finite element, which is good for vector unknowns, but \sigma is a matrix. Could you recommend a good way to treat matrix unknowns?
Thanks!
Sander
I was wondering what would be the best way to implement LDG for Stokes (see for example doi.org/10.1137/S0036142900380121 )
In LDG an auxiliary variable \sigma is introduced such that
\sigma = grad(u)
NGSolve has the VectorL2 finite element, which is good for vector unknowns, but \sigma is a matrix. Could you recommend a good way to treat matrix unknowns?
Thanks!
Sander
- rhebergens
- Topic Author
- Offline
- Junior Member
Less
More
- Thank you received: 0
6 years 4 months ago #647
by rhebergens
Replied by rhebergens on topic LDG
Hello,
I implemented it now as follows which seems to be working and still keeps the code compact:
M = L2(mesh, order = order)
X = FESpace([M, M, M, M])
L11, L12, L21, L22 = X.TrialFunction()
L = CoefficientFunction(( (L11, L12), (L21, L22) ))
LT = CoefficientFunction(( (L11, L21), (L12, L22) ))
Thanks,
Sander
I implemented it now as follows which seems to be working and still keeps the code compact:
M = L2(mesh, order = order)
X = FESpace([M, M, M, M])
L11, L12, L21, L22 = X.TrialFunction()
L = CoefficientFunction(( (L11, L12), (L21, L22) ))
LT = CoefficientFunction(( (L11, L21), (L12, L22) ))
Thanks,
Sander
- rhebergens
- Topic Author
- Offline
- Junior Member
Less
More
- Thank you received: 0
6 years 4 months ago #650
by rhebergens
Replied by rhebergens on topic LDG
Hi Christoph,
Thanks - I'll give it a go.
Sander
Thanks - I'll give it a go.
Sander
Time to create page: 0.101 seconds