ERROR in the SWIRL.PDE example !

Log Out | Topics | Search Moderators | Register | Edit Profile

Author Message   Alexandr P. (alexp) New member Username: alexp Post Number: 1 Registered: 01-2006 Posted on Tuesday, January 24, 2006 - 03:17 am:    You write ... In two-dimensional steady-state axisymmetric form, this equation becomes three component equations, radial (vr), tangential (vt) and axial (vz): vr*dr(vr) - vr^2/r + vz*dz(vr) + dr(p) = visc*[div(grad(vr)) - vr/r^2] vr*dr(vt) + vr*vt/r + vz*dz(vt) = visc*[div(grad(vt)) - vt/r^2] vr*dr(vz) + vz*dz(vz) + dz(p) = visc*div(grad(vz)) but in script You write ... EQUATIONS ... { Tangential ("Swirling") Momentum equation } vt: dens*[vr*dr(vt) - vr*vt/r + vz*dz(vt)] - visc*[div(grad(vt))-vt/r^2]=F*(vt0-vt) ... In this expression before "vr*vt/r" there should be a sign "+".