Hi Matt,
On Apr 16, 2007, at 3:01 PM, Matt Koch wrote:
Hello,
I was just wondering whether FiPy, in general, can solve a set of
PDEs including the following equations:
1) Navier-Stokes equation (transient, non-linear)
2) Heat Conduction (transient, non-linear, with phase change)
3) Species Diffusion (transient, non-linear, with phase change)
Yes to all the above. I am currently using fipy to solve a phase
field binary solid-liquid-vapor problem with full
Navier-Stokes. Examples for this problem will be publicly available
in the next month or two.
The boundary conditions would be as follows:
1) moving interface of liquid to gas, constant inflow
Should be okay, but I haven't tried it.
2) moving interface of liquid to solid, Stefan Problem, thermal
radiation
3) moving interface of liquid to solid, equivalent Stefan Problem
for species
Yes to the last two. I think the above are generally introduced
source terms that are applied to the "interface region" of the domain.
They would be implicitly handled with a phase field method, but would
require explicit construction for
the level set method. Note, that these are different from boundary
conditions in fipy that are applied to the external mesh.
I am not after exact examples or implementations, just a general
statement as to whether it might be possible to attack such a
problem successfully in FiPy at all?
Yes with caveats. fipy has efficiency limitations that we have only
partially addressed and no parallel as yet. Flow problems
with multiple species take a lot of memory and converge slowly so
this is definitely a consideration. I am solving my flow
problems in 2D with success (but not fast).
I am thinking, though, of capturing the liquid to solid interface
and the liquid to gas interface with two separate Level Sets?
I have never used the level set method in fipy with flow, I have used
it exclusively for electrodeposition problems. I may be
working with Andrew Reeve later this summer who is interested in
ground water flows on addressing some of the level set
and flow issues. The level set code as it stands only works on 1D and
2D grids and may need some tweaking for new
problems. Also, it is written entirely in python.
Hope this helps.
Cheers
Thanks and Regards,
Matt Koch
--
Daniel Wheeler
