TITLE "1-d polymer jet"

COORDINATES cartesian1 

select
CHANGELIM=0.01
REMATRIX=On 


VARIABLES       
  V
TXX
TYY

definitions
M=0.001
R=0.00122
L=83000
De=0.01
Re=40
Ug=10
Li=1


E=(2*R)/(De*L*Re*M)
qt=((0.65*R*L)/V^0.5)*( (Ug-V)^2)*((Re*(Ug-V))/V^(0.5))^(-0.81)

Initial Value
V=1*Exp(-2*x)
TXX=13000*Exp(-20*x)
TYY=0.5

EQUATIONS       

V:  dx(V) = ( (-E*(TXX-TYY))/((De*V^2)+qt))/(1-((E*(TXX+2*TYY+3))/V^2))
TXX: dx(TXX) = (1/V)*((2*dx(V)*TXX)+(2*dx(V))-(TXX/De))
TYY: dx(TYY) = (1/V)*((-dx(V)*TYY)-(dx(V))-(TYY/De))


BOUNDARIES      
region 1
start(0) 
point value (V)=1
point value (TXX)=13000
point value (TYY)=0.000001
line to (Li)


PLOTS           
  elevation(V) from (0) to (Li) 
elevation(TXX) from (0) to (Li) 
elevation(TYY) from (0) to (Li) 
END