Solving Diffusion PDE By Finite difference Method in fortran

cool2shiv
Messages
5
Reaction score
0
Hey,
I want to solve a parabolic PDE with boundry conditions by using FINITE DIFFERENCE METHOD in fortran. (diffusion) See the attachment for the problem

The problem is that there is a droplet on a leaf and it is diffusing in the leaf
the boundry conditions are
dc/dn= 0 at the upper surface of drop as well as the leaf
and
dc/dz = 0
for the bottom most layer
and the width is taken very large

Can anyone help me please in making the grid for using FDM
the Droplet is not hemispherical in shape.
 

Attachments

  • problem_ figure1.jpg
    problem_ figure1.jpg
    41.4 KB · Views: 720
Physics news on Phys.org
I told you earlier

a) Make a Mesh via matrix, arrays, selecting regions, and inputting appriopriate zeros e.g. a zeros array
b) Input your initial conditions into that array
c) discretize your equations into finite difference (conditions for the top, bottom, and central part of the body)
d) solve via time-marching, or where you solve for the next time interval from the previous.. since this is a transient problem
e) make sure your answer is reasonable
 
Thread 'Direction Fields and Isoclines'
I sketched the isoclines for $$ m=-1,0,1,2 $$. Since both $$ \frac{dy}{dx} $$ and $$ D_{y} \frac{dy}{dx} $$ are continuous on the square region R defined by $$ -4\leq x \leq 4, -4 \leq y \leq 4 $$ the existence and uniqueness theorem guarantees that if we pick a point in the interior that lies on an isocline there will be a unique differentiable function (solution) passing through that point. I understand that a solution exists but I unsure how to actually sketch it. For example, consider a...
Back
Top