[Nek5000-users] Surface source term

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

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
>