﻿ Sample Problems > Usage > Functions > vector_functions

vector_functions

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 }

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