Yes - that should work.

Set ifuservp = .true. and then, in uservp()

      if (ifield.eq.1) then ! fluid properties

            udiff = param(2)

            utrans = param(1)

      else                            ! thermal properties

            udiff = your diffusivity

            utrans = your rhoCp

      endif

Hi Paul,

Yep, there's some melting going on.  There's a porous structure that does not melt and wax trapped in the porous structure.  The wax melts and re-solidifies.  There really nothing complex happening with the geometry.  A static mesh should do the job.  The material properties are the challenge.  I think I could model the material behavior as having a non-linear specific heat capacity as a function of temperature.

Adam

Hi Adam,

Can you be more specific?  Is it melting?  

In principle, it should be possible as that was one of the original applications for Nek5000.  In practice, however, I would rate that as still a research area (at least from the Nek perspective).  I guess that melting is much easier than freezing, however.

Paul

---

From: Nek5000-users <nek5000-users-bounces@lists.mcs.anl.gov> on behalf of nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov>
Sent: Thursday, August 3, 2017 8:29:09 AM
To: nek5000-users@lists.mcs.anl.gov
Subject: [Nek5000-users] Conjugate heat transfer with complex properties

 

Hello,

I have a conjugate heat transfer problem where the solid material is a
composite that undergoes a phase transformation.  Is it possible to
model a material like this in Nek5000?

Thanks,

Adam Lange

_______________________________________________
Nek5000-users mailing list
Nek5000-users@lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users

_______________________________________________
Nek5000-users mailing list
Nek5000-users@lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users