23 Feb
2016
23 Feb
'16
2:40 p.m.
On 23/02/16 14:37, William Booker wrote:
So if I'm understanding this correctly, for the flux F , testfunction phi and a boundary \Omega_h,
I'll have \int phi F .n ds(\Omega_h) on the boundary to impose the boundary conditions weakly.
So if i include known information as the numbers and unknown information as the variables, then this would fulfil my Dirichlet bcs. In the case of non linear shallow water, then it would be
+(hu*phi)*ds(1) + (hu*phi)*ds(2) # boundary flux for height +(0.5*g*1+ hu*hu/1)*xi*ds(1)+(0.5*0.12*0.12 + hu*hu/0.12)*xi*ds(2) #boundary flux for momentum hl = 1, hr = 0.12
As the height and not the momentum are known on the boundary.
Yes, that logic looks correct. Lawrence