I find this rotate line in my sidewall boundary conditions, which I had more-or-less forgotten about over the last two years. Here is an example:
- Code: Select all
Boundary Condition 1
Name = "sidewall y=yMin"
Target Boundaries(1) = 1
Normal-Tangential Velocity = Logical True
Normal-Tangential Velocity Rotate = Logical False
Slip Coefficient 1 = Real 10.0
Slip Coefficient 2 = Real 0.00001
Slip Coefficient 3 = Real 0.00001
External Pressure = Variable depth
Real MATC "1.0 * rhoi * gravity * tx"
If I remove it then I get ice leaking out of the bottom of my side wall boundaries. But I don't know what it does and I can't find it in the documentation... I guess the rotate command somehow ensures that the normal at the bottom of the side walls doesn't point off at 45 degrees below the horizontal, which I guess would allow the ice leakage. But how exactly? Does it force the whole side wall boundary to have the same normal? Does it only effect the edges of the boundary? I need to understand how this works before I apply it to more realistic cases with less idealised geometry.