Modeling delta function as source
Hi Foamers,
I would like to know if it is possible to include a delta function as a source term in OpenFOAM. For instance, I would like to implement, div(a*gradient(phi)) = a_0*(c+(d*delta(xx_c))) where x_c is the position corresponding to the nodal point where delta function is singular and the position changes at every time step. Please let me know if we can implement this in OpenFOAM. Thanks Regards Vishal 
Hi Vishal,
I was wondering if you found an answer to this? thanks 
I assume, that you mean the dirac distribution by delta(xx0).
You can approximate delta by a field that has the value 1/V in the grid cell containing x0 and 0 elsewhere, resulting in: integral(f(x)delta(xx0))dV=f(x0). For a more accurate solution you may need to implement a jump boundary condition and split your mesh on each "side" of x0. 
You can create a volScalarField filled with 0s. Then loop over the cells and check which one of them satisfy the delta condition, so that you can change their values to 1

All times are GMT 4. The time now is 13:22. 