What I was thinking is this: A = Gas, B=Liquid
Equation 1 for A: dCA/dt=d/dz(D*dCA/dz) - ka*CA*CB + kb*CA*CB Equation 2 for B: dCB/dt=d/dz(D*dCB/dz) - kb*CA*CB + ka*CA*CB I am not sure yet of the order of the reaction rates or if I will need to add a third component (since the reaction is in fact A + B <-> C). Can you guide me in setting up such a system? Salo On Tue, Jan 22, 2013 at 10:26 AM, Daniel Wheeler <[email protected]>wrote: > 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 ] > > -- Salomon Turgman Cohen Postdoctoral Associate Cornell University (919) 341-9650
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
