I am a relatively new user to
I am a relatively new user to the OpenFOAM community, and would appreciate some feedback/direction on a new solver I will be working on in OpenFOAM as part of my PhD.
From the documentation I have seen to date, there does not appear to be an implementation of the Phase Field multiphase flow model in OpenFOAM at this time. I have been reviewing the existing solvers and just about all the documentation I can get my hands on, but I am finding that there is still not a good step-by-step tutorial out there on development of a new solver in OpenFOAM. Below is a brief outline of the application I am trying to develop: //*********************************************// Phase Field Model Most similar to level-set method, but with the order parameter governed by the Cahn-Hilliard equation which uses chemical potential gradients and the minimization of free energy to govern interface movement). Parameters: phi - typically representing mole fraction velocity field pressure field Solver Pseudocode: Step 1: Update phi: 1) Calculate chemical potential field based on the leplacian of phi and an energy of mixing term (linear equation which is a function of phi). 2) Solve the Cahn-Hilliard equation implicitely for the new phi (involves temporal term, convective term, and diffusive term with the leplacian of the chemical potential). In reality, this step could be solved in one stage, but unfortunately there is a forth order derivative involved for which the documentation has lead me to believe is not supported with OpenFOAM's grid structure (cannot find second neighbor for 5-point descretization). So instead, I plan to calculate the chemical potential, then take the leplacian to get the final results. Step 2: Solve the momentum conservation equation, with an additional source term incorporating the momentum forcing from the chemical potential gradients (similar to curvature * surface tension * surface normal vector). Step 3: Perform pressure velocity correction Step 4: Calculate residuals for phi and velocity, and confirm a divergence free velocity vield. Step 5: If all residuals are below the specified tolerance, move to the next timestep. otherwise, go back to step 1, but using the new velocity field. //*********************************************// In regards to the questions and feedback: 1) The most obvious, Is there an existing solver in OpenFOAM (or one that you know of that is not in the standard release) that could be easily adapted to meet the needs of this application? 2) Has anyone come across a comprehensive outline describing the addition of a new solver from start to finish (i.e. adding new scalar fields, defining them in the FoamX configuration files so that they can be accessed in FoamX, outlining the function of each #include statement, and describing when each is required. etc.) 3) I am a bit concerned about my intended method for solving the cahn hilliard equation via a calculated chemical potential field (cheat for 4th order derivative). I know that it will use additional resources and could introduce additional numerical instability, but I have not been able to find a viable alternative as of yet... I would appreciate some input on any complications that you might foresee, or any methods which might prove to be better suited. Thank you in advance for any help provided. |
As an additional note, the mod
As an additional note, the model will eventually need to include easily edited functions describing such properties as the mixing energy (used in the calculation of the chemical poential), and two-phase properties as a function of phi (density, viscosity, mobility).
There are also a number of constants required which need to be incorporated into the dictionary files for ease of modification. Thanks. Adam |
Hi Adam,
I have been toying
Hi Adam,
I have been toying with modeling a Cahn-Hilliard solution in openfoam as well. I'm not sure that I can be of too much help at this point, but my initial thought is that the biharmonic operator should be fairly stable out of the box, so to speak. This might be helpful as well: http://portal.acm.org/citation.cfm?id=1036150 If you get things working at all please do share. I plan to start working with this model soon, so I'll post back here as well if I get anywhere. Best, Evan |
Adam,
Still haven't had a
Adam,
Still haven't had a chance to work too much on this, but thought I'd pass along a neat opensource C code that might come in handy: http://gururajan.mp.googlepages.com/phase-field thanks to Mogadalai Pandurangan Gururajan. Evan |
Hi,
Has anybody got a OF solver for the Cahn-Hilliard equation? regards Bjorn |
Phase-Field solver built?
Hello adona058!
I have just read you post which was posted in March, well, a half year ago. I just want to inform myself if you had any succeed to build a phase field solver for openfoam? cheers, Claus |
Hi all,
I am also considering Phase-filed method for my case. I just wonder if anyone is successful in implementing that method in OF? Cheers, Duong |
I am also very interested on whether anyone has already implemented that in openfoam?
|
no news about a phase-field formulation in OF?
|
Hi,
I am pretty far on implementing the Cahn-Hilliard model into the two-phase solver interFoam. The new solver compiles and it also runs. However I have to fix some problem with the alpha limiter. The problem is the 4th order diffusion term most certainly. When this is fixed I am going to validate the solver and share it with the public. If you guys are still interested just let me know. Regards, Sebastian. |
Hi Sebastian,
How far are you in two phase solver? I am working on same but do not know how to implement 4th order derivative. Could you please share some details? Regards, Hrushi |
Quote:
I've got the same problem. Leaving this 4th order diffusion term out of the equation works and the results seem realistic. However I want to implement the whole equation so I am still working on it. I've build my solver upon the interFoam code and I used the equation from fipy documentaion: http://www.ctcms.nist.gov/fipy/examp...2DCoupled.html I introduced a second equation for Psi for the diffusion term in the alpha1 equation which is solved at first. However the 4th order diffusion term gives me weird alphas maybe because of the explicit calculation of Psi. I am pretty much stuck at this point. I can upload the solver if you're interested. Which equations are you using and how far are you? Regards, Sebastian. |
Do you get non physical values of alpha1? Because that is where I am stuck too. I am solving the diffusion equation and not the Navier-Stokes couple.
BTW, the Psi term does not look physical with second order derivative of free energy. Is there a documentation for this term? Regards, Hrushi |
Quote:
I think the second order derivative of the free energy is correct because you have to write f'(phi) in canonical form. Thats how you end up with the equation for Psi. Also this equation works in my code without the energy gradient term. Regards, Sebastian. |
Hi Sebastian,
I am interested in Phase Field Modeling in OpenFOAM. Can you upload your solver so that we can learn from it? I am a beginner and what i have learnt upto now is: -interFoam is based on VOF whereas the phase field modeling needs to be created by replacing that VOF technique with Cahn-Hilliard equation. Thanking you in advance. Yours Anil Kunwar |
Quote:
Hello Anil, the VOF or EULER methods have a different background. These methods differ in the treatment of the whole multiphase system. The VOF method uses an additional equation (alpha = phase fraction) apart from 1 set of navier stokes equations for all phases to distinguish between multiple phases whereas the EULER method uses 3N navier stokes equations for N phases. However neither can be used to simulate phase separation because the diffusion mechanisms in the navier stokes and alpha eqn are missing. The Cahn Hilliard model doesn't replace the alpha eqn of the VOF model. It just adds a diffusion term to the alpha eqn and replaces the CSF (continuum surface force) model in the navier stokes eqn, since diffusion is a primary transport mechanism during the phase separation. Are you familiar with OF programming? I am asking because I am not done with that solver. I would like to release a version which works :) Kind regards, Sebastian |
Hi Sebastian,
Quote:
I can read OpenFOAM codes and understand the things; that is i can do part of a big code. Also, i have used FiPy in my simulation works. I can understand the Cahn-Hilliard equation solved in FiPy way. The only thing that i am not ready with OpenFOAM programming is I have not built a complete code with myself alone. If you need my input, I am ready to work with you as per your suggestions and assist you for the completion of your work. You can write to me at : anilkunwar318@yahoo.com Yours Anil Kunwar |
Just to avoid confusion:
- why are we talking about Eulerian multiphase models here (meaning by this, models that are based on some sort of volume/ensemble averaging)? I think they have nothing to do with diffuse interphase models - I think C-H eq DO replace the alpha eqn of the VOF even if it is based on a totally different formalism and theory. Especially if the diffuse interphase model is used only as a interphase tracking scheme, without caring about the thermodynamical interpretation and phase separation, I don't see any problem with that (but I also don't see any major advantages in this case). However in the C-H/N-S system, also the N-S equations should be modified accordingly with a new stress term in principle however it's not wrong the idea of modifying the alpha eqn, but, as it is written in this post somewhere before, I think the main problem will be the fourth order derivative appearing in the equation Quote:
|
Hi gigilentini8,
If i have to simulate a heterogeneous system, bubble in liquid tin : can i use this C-H equation to develop an OpenFOAM solver for both purposes: -diffuse interface modeling -thermodynamic evaluation Actually, i like to simulate how the bubbles merge in stagnant liquid? Have you written any openFOAM codes regarding phase field models? Yours Anil Kunwar |
Anil,
you don't need the C-H model for your purpose. The normal VOF model can simulate the agglomeration of the bubbles also. Regards |
Hi Sebastian,
I recently came across this literature http://math.arizona.edu/~meissen/doc...interfaces.pdf to learn that Allen-Cahn and Cahn-Hilliard equations are the canonical forms of reaction-diffusion systems. Apart from bubble evolution in molten metals(250 degree centigrade); I am also interested in how a bubble prevents the growth of Intermetallic Compounds for solid Cu/Sn/Cu joint (150 degree centigrade). The main thing is : "At working temperature of 150 degree centigrade, i just have solids. I want to know whether this exclusively implements (1) phase field model (PFM) (2) Some other solid mechanics solver". Yours Anil Kunwar |
Quote:
unless you want to model phase separation or bubbles that, depending on the thermodynamic condition, have mixed miscibile/immiscible behavior (not usually the case in normal conditions) |
Anil,
the biggest problem will be to find a solver which can handle this kind of complex physics. You need a multiphase solver (> 2 phases) which can handle energy transport and phase change! There's no official solver in OF from what I know. And you probably won't find an unofficial one also. Clearly the solver I am working on is not suitable for that since it is based on interFoam, which is limited to 2 phases. I suggest you start looking for a solidification solver which includes phase change an energy transport. This would be more suitable for your case. Hope this helps. Regards, Sebastian |
Quote:
Have you released any information in conference/paper? Please provide the citation link Thanx :) |
Hey Andrew,
no, I have not released any paper yet. Are you interested in the CH solver? Regards, Sebastian |
Cahn_Hilliard solver
Dear nero235,
Im interested in the solver of the cahn-hilliard equation, since im working on a final year project to dimonstrate the equation in openfoam, how can I get the solver to try it? Thank you in advance, Maher |
nero235,
I am interested to hear how your experience with the Cahn Hilliard model in OpenFOAM concluded. Were you able to formulate a stable model, and if so, how did you treat the fourth order terms? Thanks, Joe |
Quote:
I haven't done much on this topic lately since I had to work on another problem first. I have been in contact with Holger Marschall though who is developing another solver for the Cahn-Hilliard model in OF. He has much more experience in this field and therefore I am waiting for his release of the solver, which should be soon. However there is a solver from A. Donaldson which was written for an old OF version. I've begun porting it to a 2.2 version. Check out the git repo for that: https://github.com/11101011. Regards, Sebastian |
Hi Sebastian,
I got the phaseFieldFoam code from the link (https://github.com/11101011). While trying to compile the original code in the OpenFoam version 30, i am getting a lot of errors. can you help me out rectifying those errors or do you have the latest version of the code for the OpenFOAM version 30 ? regards |
phase field confusion
i have some confusion with phase field and hope you solve my confusion.I don't know whether phase field is suitable for the calculation and simulation of multi-phase turbulence. Most of the literature is VOF and level set. I hope you can provide some help~
|
FWIW, I managed to compile Sebastian's port of Adam Donaldson's code with openFoam 6, blindly adapting a few lines. The fork is here. Adam's most relevant paper describing the code is probably this one, btw.
With that I won't say that I'm sure that the computation is correct, only that the code seems to run plausibly. For the cases I'm currently interested into, i.e. microfluidics, it may be not the right approach. At the sizes of my system, surface tension sigma and contact angle theta seem to have no effect. Maybe my port missed something, maybe it's a limitation of the modelling. There are another couple of research codes out there which I still want to try, before resigning to VOF. Feedback appreciated! https://www.cfd-online.com/Forums/mo...r-marker-1.pnghttps://www.cfd-online.com/Forums/mo...r-marker-1.png |
Quote:
i'm also interested in Holger Marschall's solver, have you got it? |
Computational modelling using Phase Filed method in OpenFOAM
Hello All,
Has anyone developed a solver using phase field method to solve C-H equation and diffusion equation for solidification of a binary alloy? Thanks |
Phase Field modeling in OpenFOAM
Hello Foamers,
Have you been able to develop Phase-Field solver using OpenFOAM. I am interested in the solver that solves the Cahn-hilliard 4th order diffusion equation. Please let me know. Thanks |
Quote:
Have you been able to develop the PhaseField solver using OpenFOAM ? |
Quote:
|
All times are GMT -4. The time now is 22:50. |