(* Modeling the reaction drain-pipe, see Figs.14.2 and 14.6*)
Dx=20
Dy=20
Rex=5
Rey=5
alfx=2
alfy=1
V1=Dx (((((Rex/x)))^12)-2 ((((Rex/x)))^6))+Dy (((((Rey/y)))^12)-2 ((((Rey/y)))^6))
V2=50 (E^(-0.1 (((x+y-(Rex/0.6)))^2)))
V3=V1+V2
Plot3D[V3,{x,4.35,10.},{y,4.35,10.}]
ContourPlot[V3,{x,4.4,10.},{y,4.2,10.},PlotPoints->60]
Plot3D[V3,{x,4.35,10.},{y,4.35,10.},ViewPoint->{6,6,3}]
Dx=20
Dy=30
Rex=5
Rey=5
alfx=2
alfy=1
V1=Dx (((((Rex/x)))^12)-2 ((((Rex/x)))^6))+Dy (((((Rey/y)))^12)-2 ((((Rey/y)))^6))
V2=50 (E^(-0.1 (((x+y-(Rex/0.6)))^2)))
V3=V1+V2
Plot3D[V3,{x,4.35,10.},{y,4.35,10.}]
ContourPlot[V3,{x,4.4,10.},{y,4.2,10.},PlotPoints->60]
Dx=30
Dy=20
Rex=5
Rey=5
alfx=2
alfy=1
V1=Dx (((((Rex/x)))^12)-2 ((((Rex/x)))^6))+Dy (((((Rey/y)))^12)-2 ((((Rey/y)))^6))
V2=50 (E^(-0.1 (((x+y-(Rex/0.6)))^2)))
V3=V1+V2
Plot3D[V3,{x,4.35,10.},{y,4.35,10.}]
ContourPlot[V3,{x,4.4,10.},{y,4.2,10.},PlotPoints->60]
Plot3D[V3,{x,4.35,10.},{y,4.35,10.},ViewPoint->{6,6,3}]
Plot3D[V3,{x,4.35,10.},{y,4.35,10.},ViewPoint->{6,10,3}]
Dx=30
(*Dy=25*)
Dy=30
Rex=5
Rey=5
alfx=2
alfy=1
V1=Dx (((((Rex/x)))^30)-2 ((((Rex/x)))^16))+Dy (((((Rey/y)))^8)-2 ((((Rey/y)))^2))
V2=50 (E^(-0.2 (((x+y-(Rex/0.6)))^2)))
V3=V1+V2
Plot3D[V3,{x,4.35,10.},{y,4.35,10.}]
ContourPlot[V3,{x,4.4,10.},{y,4.2,10.},PlotPoints->60]
Plot3D[V3,{x,4.35,10.},{y,4.35,10.},ViewPoint->{10,6,3}]
Plot3D[V3,{x,4.35,10.},{y,4.35,10.},ViewVector->{{20,20,3},{5,5,3}}]
Plot3D[V3,{x,4.35,10.},{y,4.35,10.},ViewPoint->{2,6,3}]
Plot3D[V3,{x,4.35,10.},{y,4.35,10.},ViewPoint->{6,10,3.5}]
Plot3D[V3,{x,4.35,10.},{y,4.35,10.},ViewPoint->{6,10,3}]
Plot3D[V3,{x,4.35,10.},{y,4.35,10.},ViewPoint->{6,10,2.5}]
,ViewVector->{{5,-5,5},{5,0,-5}}