Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Using dirac delta function as an initial condition

Please login with a confirmed email address before reporting spam

Hi,

I am trying to model the effect of gravity on brownian motion for which the general equation from Smoluchowski is

dw/dt = div . (qB^-2grad.w - KB^-1w).

Considering only the z direction then the simplified form is

dw/dt = D(d2w/dz2) + C(dw/dz) where D = kT/mB and C = (1-(p1/p2))g/B ( terms constituting D and C are known)

the boundary conditions would be
w-> dirac delta(z - z0) where z0 is the height of fluid measured from the bottom where z = 0.
D(dw/dz) + Cw = 0 at z = 0 for all time ( prevention of particles crossing z = 0 )

Currently i have tried using pde convection-diffusion equation (cdeq) (version 3.5a) and the initial value seems reasonable but questionable due to the negative u. Also, the solution does not progress much from the initial value solution as the expected result would be a curve shifting in time from from one extreme to the other in the shape of / similar to a gaussian distribution.

I'm not sure how to implement the dirac function as a initial concentration and i've tried using flc1hs(z-z0,1e-5) as the initial condition for the subdomain setting.

Any help would be welcome and thanks in advance.

Regards,
Chun Wei


0 Replies Last Post 20 set 2010, 13:48 GMT-4
COMSOL Moderator

Hello Tan Chun Wei

Your Discussion has gone 30 days without a reply. If you still need help with COMSOL and have an on-subscription license, please visit our Support Center for help.

If you do not hold an on-subscription license, you may find an answer in another Discussion or in the Knowledge Base.

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.