Dear colleagues; I have trouble creating a jump over an edge for my field, ψ, a P2 element in FreeFem. I am solving, where q and μ are known, and r is a multiplier: (see the eqn and geom attached)
I want a difference of ψ over the blue dashed line: i.e., ψ|below = ψ|above +2 π, for instance, while ∇ψ remains continuos on both sides (without cutting or any alterations to the domain).
What terms some I add to the
solve Mahdi( ψ , v ) = int2d(Th)( r*(dx( ψ )*dx(v) + dy( ψ )dy(v)) ) + int2d(Th)( ( μ x-rqx )dx(v) + ( μ y-rqy )*dy(v) )
Writing as: solve Mahdi( ψ , v )= int2d(Th)( r*(dx( ψ )*dx(v) + dy( ψ )dy(v)) ) + int2d(Th)( ( μ x-rqx )dx(v) + ( μ y-rqy )dy(v) ) + intalledges(Th)(label==1) v * (2 *pi)*jump( ψ );
Where the label for the blue dashed line is 1, it didn’t solve the problem.
- Badges
- Science topic
- Similar topics
- Mathematics