On Mon, Jan 21, 2013 at 7:15 PM, Salomon Turgman Cohen <[email protected]>wrote:
> Hello list, > I have modeled two approximations of my system as 1D diffusion > equations. This is an approximation since in reality the system is > reactive, but I am sweeping the details of the reaction under the rug. Now > I would like to solve a two component system with reaction. The description > is the following: > > 1. I have a liquid and a gas (species A). The gas will be absorbed by the > liquid. > 2. Within the liquid the gas will react with the liquid (1:1 > stoichiometry) to form species B > 3. I have measured the diffusivity of species A and the diffusivity of > species B as a function of the concentration of species B > 4. The resulting model I am picturing would be two diffusion equations, > both with source terms to represent the reaction. > Can you write down the PDEs that you want to solve? Is the issue with writing down the equations of is the issue implementing the equations in FiPy? > How do I go about building a two component system in fipy. I found an > example elphf/diffusion/mesh1D.py which gives some idea, I just don't > understand enough about python yet. Will I need to create a new class to > accomplish this? > No. I assume it is just a couple set of PDEs. There are different ways to represent couple equations if FiPy, if you write down the equations we can get a better idea how to implement them. Cheers -- Daniel Wheeler
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
