| Author | Message | ||
     Kevin Ellwood (ikevin) Member Username: ikevin Post Number: 13 Registered: 06-2003 |
Hello All I have a reaction diffusion problem in which there are several different material regions. Each region has its own kinetic parameters and diffusivities and I have a first level of the model working fine. Now I wish to modify this model such that the diffusivities are governed by: dD/dt = function of species I have handled this for a single material by introducin the following artificial diffusion equation into flexpde: dD/dt = div(k*grad(D)) + function of species where k is a small number (not too small :-) ). I have tested it and it works. My problem is how do I handle the different initial conditions in the case of different materials? Specifying "D" as being an unkown promotes it to a nodal quatity, which also causes some issues. This is really an element quatity that should be computed at quass points within each element. Are there any strategies/tricks for handling this? Thanks for any comments Kevin | ||
| Henadzi Filipenka (hfil) Junior Member Username: hfil Post Number: 3 Registered: 05-2004 |
Please see my paper at : www.sciteclibrary.ru "On the problem of crystal metallic lattice..." Best regards!Henadzi Filipenka |