<< Click to Display Table of Contents >> Navigation:  Sample Problems > Usage > Functions > vector_functions

{ VECTOR_FUNCTIONS.PDE  

 

 This example illustrates the vector functions

 

 VECTOR

 MAGNITUDE

 DOT

 CROSS

 NORMAL

 TANGENTIAL

 

}

title

"vector functions"

select

 elevationgrid=500

{no variables}

definitions

 u= exp(-x^2+ y)               { A scalar potential, perhaps }

 f= grad(u)                   { F = grad(u)  is a vector }

 df= div(f)                   { Divergence of F is a scalar}

 cf= curl(f)                   { Curl of F is a new vector }

 vx= -sin(y)   vy= 2*sin(x)   { vector components }

 v= vector(vx,vy)             { Another vector }

 mv= magnitude(v)             { Magnitude of v }

 cv= curl(v)

 ccv= curl(curl(v))

 tvv = v*v { v*v is a tensor }

 divtx = 2*vx*dx(vx)+vx*dy(vy)+vy*dy(vx)   {x-component of div(t) }

 divty = vx*dx(vy)+vy*dx(vx)+2*vy*dy(vy)   {y-component of div(t) }

 divt = vector(divtx,divty)

{no equations}

{plot domain -- required}

boundaries

region 1

  start "Outer" (-1,0)

  line to (1,0) to (1,1) to (-1,1) to close

feature

  start "inner" (-1/2,1/2) line to (1/2,1/2)

plots

vector(f)

elevation(normal(f)) on "Outer"

elevation(tangential(f)) on "inner"

contour(df)   as "Div F"

contour(mv)   as "Magnitude V"

contour(dot(v,vector(x,0)))

contour(zcomp(cross(f,v)))

contour(zcomp(cv))   as "Curl V"

vector(ccv)   as "Curl Curl V"

vector(div(v*v)) as "Div(V*V) inline"

vector(divt) as "Div(V*V) expanded"

vector(div(tvv)) as "Div(V*V) tensor parameter"

end