An Example with Explicit Velocity Vectors 

Although completely general solutions are not forthcoming at this point, we can still evaluate some commonly used expressions that are often needed.  To demonstrate this, we first need to define the velocity vector v(x,y,z).  I use completely arbitrary definitions for the velocity components in terms of x, y, z and the unknown function Δ(x).  These functions were inspired by the Maple code of professor Davis as seen in ex1241.mw although they do not correspond exactly. 

 

> vx := unapply(vinf*`ϕ`(y/(x^2+C*y)), x, y, z); 1
 

(Typesetting:-mprintslash)([vx := proc (x, y, z) options operator, arrow; vinf*`ϕ`(y/(x^2+C*y)) end proc], [proc (x, y, z) options operator, arrow; vinf*`ϕ`(y/(x^2+C*y)) end proc]) 

> vy := unapply((D[1](vx))(x, y, z), x, y, z); 1
 

(Typesetting:-mprintslash)([vy := proc (x, y, z) options operator, arrow; -2*vinf*(D(`ϕ`))(y/(x^2+C*y))*y*x/(x^2+C*y)^2 end proc], [proc (x, y, z) options operator, arrow; -2*vinf*(D(`ϕ`... 

> vz := unapply(10, x, y, z); 1
 

(Typesetting:-mprintslash)([vz := proc (x, y, z) options operator, arrow; 10 end proc], [proc (x, y, z) options operator, arrow; 10 end proc]) 

We can see that the components of v have now been assigned 

> v(x, y, z); 1
 

(Typesetting:-mprintslash)([Vector[column]([[vinf*`ϕ`(y/(x^2+C*y))], [-2*vinf*(D(`ϕ`))(y/(x^2+C*y))*y*x/(x^2+C*y)^2], [10]], [ 

 

Now we can see what some of the terms will look like at an arbitrary point (x,y,z) 

 

> eval(`τxy`(x, y, z)); 1
 

-mu*(4*vinf*((`@@`(D, 2))(`ϕ`))(y/(x^2+C*y))*y^2*x^2/(x^2+C*y)^4+8*vinf*(D(`ϕ`))(y/(x^2+C*y))*y*x^2/(x^2+C*y)^3-2*vinf*(D(`ϕ`))(y/(x^2+C*y))*y/(x^2+C*y)^2+vinf*(D(`ϕ`))(y/(...
-mu*(4*vinf*((`@@`(D, 2))(`ϕ`))(y/(x^2+C*y))*y^2*x^2/(x^2+C*y)^4+8*vinf*(D(`ϕ`))(y/(x^2+C*y))*y*x^2/(x^2+C*y)^3-2*vinf*(D(`ϕ`))(y/(x^2+C*y))*y/(x^2+C*y)^2+vinf*(D(`ϕ`))(y/(...
-mu*(4*vinf*((`@@`(D, 2))(`ϕ`))(y/(x^2+C*y))*y^2*x^2/(x^2+C*y)^4+8*vinf*(D(`ϕ`))(y/(x^2+C*y))*y*x^2/(x^2+C*y)^3-2*vinf*(D(`ϕ`))(y/(x^2+C*y))*y/(x^2+C*y)^2+vinf*(D(`ϕ`))(y/(...
 

> `τzy`(x, y, z); 1
 

0 

> `τxx`(x, y, z); 1
 

4*mu*vinf*(D(`ϕ`))(y/(x^2+C*y))*y*x/(x^2+C*y)^2+(2/3*mu-kappa)*(D[1](vinf)*`ϕ`(y/(x^2+C*y))+vinf*D[1](`ϕ`(y/(x^2+C*y)))-2*D[2](vinf)*(D(`ϕ`))(y/(x^2+C*y))*y*x/(x^2+C*y)^2-2...
4*mu*vinf*(D(`ϕ`))(y/(x^2+C*y))*y*x/(x^2+C*y)^2+(2/3*mu-kappa)*(D[1](vinf)*`ϕ`(y/(x^2+C*y))+vinf*D[1](`ϕ`(y/(x^2+C*y)))-2*D[2](vinf)*(D(`ϕ`))(y/(x^2+C*y))*y*x/(x^2+C*y)^2-2...
4*mu*vinf*(D(`ϕ`))(y/(x^2+C*y))*y*x/(x^2+C*y)^2+(2/3*mu-kappa)*(D[1](vinf)*`ϕ`(y/(x^2+C*y))+vinf*D[1](`ϕ`(y/(x^2+C*y)))-2*D[2](vinf)*(D(`ϕ`))(y/(x^2+C*y))*y*x/(x^2+C*y)^2-2...
4*mu*vinf*(D(`ϕ`))(y/(x^2+C*y))*y*x/(x^2+C*y)^2+(2/3*mu-kappa)*(D[1](vinf)*`ϕ`(y/(x^2+C*y))+vinf*D[1](`ϕ`(y/(x^2+C*y)))-2*D[2](vinf)*(D(`ϕ`))(y/(x^2+C*y))*y*x/(x^2+C*y)^2-2...
4*mu*vinf*(D(`ϕ`))(y/(x^2+C*y))*y*x/(x^2+C*y)^2+(2/3*mu-kappa)*(D[1](vinf)*`ϕ`(y/(x^2+C*y))+vinf*D[1](`ϕ`(y/(x^2+C*y)))-2*D[2](vinf)*(D(`ϕ`))(y/(x^2+C*y))*y*x/(x^2+C*y)^2-2...
 

>
 

 

These results, although lacking full generality, are still useful in many applications, especially since the results can be easily mapped to other coordinate systems as already demonstrated.