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.
_______________________________________________
Mailing list: https://launchpad.net/~dolfin
Post to : [email protected]
Unsubscribe : https://launchpad.net/~dolfin
More help : https://help.launchpad.net/ListHelp
