[shinigamirem]

Dear OpenFOAMERS,

I am very new to OpenFOAM and in the last months i have established a simulation of Fick's diffussion equation. My interest is to simulate the ion concentration in a solid electrolyte. For this i used the laplacianFoam and changed the temperature to a concentration with a constant diffusion coefficient in the transport properties. My mesh is so far one block created with blockMesh.

Now i want to further improve my simulation and i am not sure how to proceed. My new improvements i want to implement are:

• including the electrodes which is a different material (solid) as new blocks
• having two different diffusion coefficients for the two materials. Later on i want to make them concentration dependent
• including a migration term into the simulation from the Nernst-Planck equation
• limit the concentration as negative concentration makes no sense
• including electron diffusion depending on ions

This should be a preview where i want to go, so the recommendations can be based on them.

The limit of concentration is what i am working on right now, which i wanna solve through including the fvOption. Therefore i would need a thermophysical model in OpenFOAM. Hard to choose for a beginner, which is why i reaching out for help. I am still not very settled with OpenFOAM yet, which is why i have not a good overview on the possibilities and resources which may fit best for me.

I hope someone can help and give advice and/or recommendations

Thank you very much!