[Nek5000-users] nelgv in uservp of conj_ht

[nek5000-users at lists.mcs.anl.gov]

Hi Sandeep,

The temperature should be continuous if everything is working correctly.

One way things can go wrong is if the initial temperature is not continuous, because subsequent
fields are computed as updates of the original.   All corrections are continuous:

    T^n = T^n-1  + dT

dT is continuous.   T^n will be continuous only if T^n-1 is continuous.   useric() generally ensures
continuity of temperature and velocity via the dsavg operation before the start of the calculation.

Even if the user prescribes a discontinuous initial condition, useric() invokes dsavg to force continuity.

How do you know that the function is discontinuous ?

Paul

________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov [nek5000-users-bounces at lists.mcs.anl.gov] on behalf of nek5000-users at lists.mcs.anl.gov [nek5000-users at lists.mcs.anl.gov]
Sent: Tuesday, January 24, 2017 8:27 AM
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] nelgv in uservp of conj_ht

Hello All,
I have a question in uservp of conjugate heat transfer. I am trying to modify the conjugate heat transfer problem (with additional passive scalars) and use for solving MHD problem with vector potential approach.

I tested the conjugate heat transfer case changing the IF condition in uservp in following ways.

Case 1:
        if (ieg .gt. nelgv) then  ! properties in the solid
            udiff   = 0.1*param(8) ! conductivity
            utrans  = 1.0
         endif

Case 2:
if (y .ge. 1.5D0) then
udiff= 0.1*param(8)
elseif (y. le. 0.5D0) then
udiff= 0.1*param(8)
else
udiff   = param(8)
endif

I got a small difference in the temperature profile at the interface between solid and fluid.
Could you please let me know, how do the above two cases affect the conditions imposed at the interface (what is the mathematical treatment and how does it change).
Thank you in advance.

Yours sincerely,
Sandeep
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170124/e49fb4ed/attachment.html>