[Nek5000-users] Variable density flows

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Fri Oct 11 13:53:08 CDT 2013


Hi Paul and Praveen

div U=0    In the low Mach limit,  in the absence of external heat sources q' , time dependent variation of the background pressure dPo/dt, and reactions source terms 

and Mass conservation reduces to  conservation of density along particle path so D rho /Dt=  drho/dt + u. grad rho= 0


Ammar





On Oct 11, 2013, at 1:27 PM, nek5000-users at lists.mcs.anl.gov wrote:

> 
> Hi Praveen,
> 
> If you want to follow the model of the Guermond paper that you
> pointed out then I shouldn't think low Mach is necessary given
> that, in their model, div U = 0.
> 
> I've not looked into this in detail, however, so I suggest
> undertaking a bit of exploration and test.
> 
> To specify variable density in Nek, you set ifuservp=.true.
> and ifstrs=.true. in usrdat, and then in uservp() you set:
> 
>      if (ifield.eq.1) then
> 
>         utrans = mydensity
>         udiff  = my viscosity
> 
>      endif
> 
> You also need to set lx2=lx1-2, (same for ly2 and lz2 (if3d));
> and set lx1m=lx1,ly1m=ly1,lz1m=lz1 in the SIZE file.
> 
> Paul
> 
> 
> ----- Original Message -----
> From: nek5000-users at lists.mcs.anl.gov
> To: nek5000-users at lists.mcs.anl.gov
> Sent: Friday, October 11, 2013 11:54:38 AM
> Subject: Re: [Nek5000-users] Fwd: Variable density flows
> 
> 
> 
> 
> Hi Paul 
> 
> For variable density case do I need to use the low mach number version or some other way ? 
> 
> 
> Thanks 
> praveen 
> 
> 
> 
> On Fri, Oct 11, 2013 at 8:24 PM, < nek5000-users at lists.mcs.anl.gov > wrote: 
> 
> 
> 
> In principle, this should work with the Pn-Pn-2 STRESS formulation, but 
> I don't think it has been thoroughly tested. 
> 
> Regarding large viscosity variations, work needs to be done on the 
> preconditioners for the viscous terms. This hasn't been a traditional 
> focus for Nek5000, which has been aimed at high Reynolds number flows 
> where the viscous solves are essentially trivial (compared to the 
> pressure, say). 
> 
> Certainly, there are things that could be done to be effective in 
> this space. I would start by running some test problems with the 
> STRESS formulation to see if you get the accuracy you're after. If 
> so, then we can look into performance enhancement via better preconditioning 
> strategies for the case at hand. 
> 
> HTH, 
> 
> Paul 
> 
> 
> 
> ----- Original Message ----- 
> From: nek5000-users at lists.mcs.anl.gov 
> To: "nek5000" < nek5000-users at lists.mcs.anl.gov > 
> Sent: Friday, October 11, 2013 9:16:24 AM 
> Subject: [Nek5000-users] Fwd: Variable density flows 
> 
> 
> 
> Hello 
> 
> 
> 
> 
> 
> Can we solve variable density incompressible flows with nek5000 ? For e.g. see the model given by eq. 2.1 here 
> 
> 
> http://www.math.tamu.edu/~guermond/PUBLICATIONS/guermond_quartapelle_JCP_2000.pdf 
> 
> 
> 
> Also, can nek handle variable viscosity with large variations ? 
> 
> 
> Thanks 
> praveen 
> 
> _______________________________________________ 
> 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 
> 
> 
> _______________________________________________
> 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



More information about the Nek5000-users mailing list