[Nek5000-users] Surface source term

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

Hi Paul,

             yes, I use also "iel = gllel(ieg)".

I have just checked the energy balance and it is fine.

Thanks a lot for your answer.

Regards,
Andreas.
________________________________________
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: Friday, March 09, 2012 6:10 PM
To: nek5000-users at lists.mcs.anl.gov
Subject: Re: [Nek5000-users] Surface source term

Hi Andreas,

That looks right to me...  Did you also use:

       iel = gllel(ieg)

where ieg is the argument passed into the call line?

Also, do you get the expected energy balance in a test case?

Best,

Paul



On Fri, 9 Mar 2012, nek5000-users at lists.mcs.anl.gov wrote:

> Dear NEKs,
>
>                I am trying to add a surface source term in the energy equation on a certain surface of my 2D domain.
>
> The way I did was to put into userq the following expression:
>
> qvol = qref*area(ix,1,ifc,iel)/bm1(ix,iy,iz,iel)
>
> The case seems to work fine but I would like to ask you just a confirmation that the way I did is correct (maybe there are some mathematical/physical reasons why that way is wrong).
>
> Thanks in advance for your help.
>
> Regards,
> Andreas.
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users