# Problem modelling potential equations with addtional variables

 Register Blogs Members List Search Today's Posts Mark Forums Read

 October 2, 2013, 07:39 Problem modelling potential equations with addtional variables #1 New Member   Tom Tranter Join Date: Mar 2013 Location: Leeds, UK Posts: 21 Rep Power: 6 Hi - I have a problem I'm hoping you can help with. Short Version: How do you solve a transport equation for a scalar additional variable that is not volumetric or specific e.g an electric potential. I know CFX has the electric potential functionality already but I need to define a new one to represent proton transport as well as electron transport. Converting units to get diffusivity in m^2/s incurs problems where additional variable is coupled to flow when it shouldn't be. Long Version: To give some context I am trying to model a fuel cell in CFX from scratch and this firstly requires solving two potential equations for electron transport through solid and porous domains and proton transport through a porous membrane. The governing equations are div(grad(PhiS)) + Rs = 0 and div(grad(PhiM)) + Rm = 0. On either side of the cell Rs=-Rm and has different magnitudes on either side. The equations are coupled to each other via the transfer current (Rs, Rm) source terms but at the moment not coupled to any fluid flow. Eventually the current transfer will be linked to species transport as hydrogen is split into protons and electrons on the anode side of the cell and recombines with oxygen on the cathode side of the cell to form water. I have created two additional variables, PhiS and PhiM which at first I specified as volumetric and have now swithed to specific but still get problems. S stands for solid and M stands for membrane. The units for the potentials are volts and the sources are in A/m^3. The fluid and solid models for the potentials in CFX are both poisson as the situation is not time-dependent. My first difficulty is that to solve the equations I must convert from Volts to density as if I am solving a species diffusion equation rather than an electric potential as CFX does not allow me to specify a general diffusivity in units of an electrical conductance (S/m) but rather requires the standard units of m^2/s. I have done this by multiplying the source terms by an expression called unitFactor = 1[kg V A^-1 s^-1]. Now when running the simulation to solve the equations I find that my potentials have become dependent on the fluid flow in the domain somehow when they should be independent. I have uploaded the ccl and an image of the problem which is a cross-section of the cell to dropbox. https://www.dropbox.com/sh/adc8n4805hvcwir/QhEwZEx-eg?m The values of the electric potential are set by boundary conditions at the top an bottom and should remain constant until the catalyst layers either side of the central membrane region where the transfer current source terms govern the membrane potential. I'm sure this general type of problem has been encountered before and I feel I'm close to a solution but may have done something wrong along the way. Thanks

 October 2, 2013, 17:48 #2 Super Moderator   Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 14,098 Rep Power: 109 If you use an unspecified additional variable it can be in any units you like and you define with with a CEL expression. Will this do what you want? I would also ask ANSYS support on this. They may well have an example of a similar simulation they can give you which will get you started. This is likely a better way to get started rather then doing it from scratch.

 October 3, 2013, 05:27 #3 New Member   Tom Tranter Join Date: Mar 2013 Location: Leeds, UK Posts: 21 Rep Power: 6 Hi Glenn, Thanks for the reply. Changing the variable to unspecified only leaves me with algebraic equation as the option to solve. How would I then implement the poisson equation I need to solve? I have contacted Ansys as you suggest but as this is for my PhD I would like to start from scratch to gain the best understanding. p.s for completeness I forgot to include the conductivity (sigma) in the governing equations div(sigma*(grad(PhiS))) + Rs = 0

 October 3, 2013, 05:33 #4 New Member   Tom Tranter Join Date: Mar 2013 Location: Leeds, UK Posts: 21 Rep Power: 6 Actually, forget I asked that last question. On closer inspection it looks like I had an error in one of my source terms and re-running the simulation now produces reasonable results. Sorry to waste your time.

 October 3, 2013, 08:20 #5 Super Moderator   Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 14,098 Rep Power: 109 Good to see you have some progress. But a general comment - a PhD does not mean you do things ignoring what others have done. You do not get a PhD by repeating what somebody else has done. pmtgt likes this.

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post bhh OpenFOAM Running, Solving & CFD 11 December 10, 2016 12:25 pete Site Help, Feedback & Discussions 20 June 16, 2015 13:52 subhsngh OpenFOAM Meshing & Mesh Conversion 0 July 6, 2013 16:30 aleh Main CFD Forum 0 February 3, 2009 20:59 Yang CFX 2 May 18, 2007 06:22

All times are GMT -4. The time now is 06:36.