dolfin team mailing list archive

["B. Emek Abali" <question152565@xxxxxxxxxxxxxxxxxxxxx>]

Question #152565 on DOLFIN changed:
https://answers.launchpad.net/dolfin/+question/152565

B. Emek Abali proposed the following answer:
The point source is a mathematical trick of trying the make the
analytical solution possible, FEniCS can do it already better, just use
a Gauss bell-shape

pointload = Expression("A*exp(-B*(pow(x[0]-v*t, 2)+pow(x[1]-v*t, 2) )
)")

A is amplitude
has its maximum for time t at x-vt
so it is also moving in x,y plane
decaying due to minus sign as strong as the B

on the right hand side contract it with surface normal and test function, in a time loop don't forget to update it
pointload.t = ..

-- 
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.