[Nek5000-users] Laplacian as a forcing function and stability issues

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Thu Nov 4 17:52:50 CDT 2010


Hi Paul,

Thanks for offering to look into this. I have e-mailed you the files.

Regards
Shriram


On 4 November 2010 16:30, <nek5000-users at lists.mcs.anl.gov> wrote:

>
> Hi Shriram,
>
> Do you want to send .rea/.usr/SIZE in a tarfile?
>
> I can take a look... I've been having good luck w/ dyn. smag.
> of late.
>
> Paul
>
>
>
> On Thu, 4 Nov 2010, nek5000-users at lists.mcs.anl.gov wrote:
>
>  Hi,
>>
>> I am trying to implement a LES model which requires a term like c . Sij*
>> to
>> be added to the existing unfiltered viscous term. I am computing Sij* by
>> filtering the velocity and calling comp_gije() twice. This is then
>> included
>> in the forcing function. This seems to work for few steps but the solution
>> starts to blow up afterwards. I am wondering if this would be due to
>> stability issues since the laplacian here is not calculated by any of the
>> splitting methods. I would like to know if it would be possible to add
>> this
>> term to the Navier Stokes formulation that is solved each step. Any other
>> ideas/ thoughts would be helpful. Thanks.
>>
>> Regards
>> Shriram
>>
>>  _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20101104/89112224/attachment.html>


More information about the Nek5000-users mailing list