On May 19, 2013, at 12:23 PM, Michelangelo Formisano <[email protected]> wrote:
> I'm Michelangelo, an Italian post-doc, and I recently discovered FIPY. Welcome to FiPy. We hope it's useful to you. > I would like to develop a 3D thermal/geophysical model of asteroids by using > FIPY. My idea is to solve simultaneously the heat equation (with the source > term represented by the radionuclides) and the advection equation (which > controls the migration of metals through the silicatic matrix). I have > already developed a 1D finite difference method by using Python, but I'd like > to move on to making a 3D model and Fipy sounds great. > > Do you suggest me some examples to follow and a way to incorporate a > radiation boundary condition at the surface? Any of the diffusion examples in http://www.ctcms.nist.gov/fipy/examples/diffusion/index.html should get you started for both the heat equation and the diffusive part of metal migration. The convection examples in http://www.ctcms.nist.gov/fipy/examples/convection/index.html will show how to do the convective part. Many of the phase field examples in http://www.ctcms.nist.gov/fipy/examples/phase/index.html show one approach to dealing with phase transformation and coupling with solute transport. If you are interested in couplings between heat and species transport, e.g. Soret effect, I was involved in one such study: Mohanty et al. Diffusion under temperature gradient: A phase-field model study. Journal of Applied Physics (2009) vol. 106 (3) pp. 034912 and could probably dig up the FiPy code that was used. Radiation boundary conditions were recently discussed here: http://thread.gmane.org/gmane.comp.python.fipy/2942 _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
