|February 15, 2016, 03:36||
Cahn Hilliard implementation (Biharmonic operator)
Join Date: Feb 2016
Location: Santa Barbara, CA, USA
Posts: 2Rep Power: 0
I'm interested in modeling viscoelastic phase separations in simple flows. As a first step in the process, I want to implement a simple diffusive phase separation such as the Cahn Hilliard equation using OpenFOAM.
There are at least several threads with people asking for advice for solving the Cahn Hilliard equation, and as of last active discourse (~2012) there was no publicly distributed resolution to the main problem (an implicit treatment of the fourth order derivative in the problem)
A link to one of the prior threads on this topic:
The grand challenge is that, for a stable implementation, a fourth order spatial derivative (the biharmonic operator) must be calculated implicitly at some stage. One common way to accomplish this is to solve simultaneously for two variables such that the model assumes the following structure:
x = laplacian(y)
dx/dt = laplacian(x + y)
I'm new to the environment, but so far it seems that in OpenFOAM the solver only solves a single field at a time; can I trick it into solving for these two variables simultaneously? Is there another way that I should be thinking about the problem?
Any help or insight into this challenge would be greatly appreciated!
|biharmonic, cahn hilliard, phase field, phase separations|
|Thread||Thread Starter||Forum||Replies||Last Post|
|ANSYS Workbench on "Uncertified" Linux Distros||hsr||CFX||262||September 7, 2016 21:05|
|Create a biharmonic operator||Pascal_doran||OpenFOAM Programming & Development||12||May 22, 2016 06:02|
|FLUENT installation on UBUNTU 12.04 (LTS)||teymourj||FLUENT||0||August 26, 2013 20:27|
|cahn hilliard equation finite volume method||dvdromnu||OpenFOAM Programming & Development||0||December 14, 2009 03:09|
|Cahn Hilliard couple with Navier Stokes||Khan||Main CFD Forum||0||September 2, 2006 02:19|