# Passive scalar is diffusive without mesh motion

 Register Blogs Members List Search Today's Posts Mark Forums Read

 July 23, 2015, 08:59 Passive scalar is diffusive without mesh motion #1 Super Moderator     Tobias Holzmann Join Date: Oct 2010 Location: Leoben (Austria) Posts: 1,947 Blog Entries: 6 Rep Power: 33 Hi all, I have a question to all of you. Maybe you are more familiar with that problem. The case is very simple. A 2d rectangular with an outlet velocity of 1m/s, p is fixed at the inlet with fixedValue; the walls are zeroGradient for U and p. The passive scalar is set to zeroGradient at each boundary. With setFields I generate an area with values = 1 for the passive scalar (everywhere else its zero). Case 1 has no mesh motion Case 2 has mesh motion due to the velocity field at the inlet (flux after the iteration = 0 at the inlet face). The position of the scalar is the same (as expected) but why is without mesh motion the scalar diffusive? The transport of the scalar is only done by time and convection (no diffusion term). In the mesh motion I added the mesh flux to the solver, but can this be the reason for less diffusion? But if I check the code: Code: ```fvm::ddt(S) +fvm::div(phi, S)``` There should be no diffusion in the case with no mesh motion. S is only transported by the flux, which is everywhere constant. :/ Here you can check the results: www.holzmann-cfd.de/cfd-online/passiveScalar.avi Any hint is appreciated. Thanks in advance, Tobi __________________ Keep foaming, Tobias Holzmann

 July 23, 2015, 09:46 #2 Senior Member   Join Date: Oct 2013 Posts: 382 Rep Power: 8 Are you sure there should be no diffusion at all? As far as I know, there is always some artificial, numerical diffusion (though I'm not sure if this requires that an actual laplacian is present in the equation). I have no idea why the diffusion is smaller with mesh motion, but I have no experience with dynamic meshes at all.

 July 23, 2015, 10:32 #3 Super Moderator     Tobias Holzmann Join Date: Oct 2010 Location: Leoben (Austria) Posts: 1,947 Blog Entries: 6 Rep Power: 33 Hi chriss, you are right. There is always Diffusion, hence there is only convection. But like you mentioned, I also don't know why I get the better result with mesh motion. Maybe it's due to the fact that the mesh flux reduces the convection term itself. (U-Umesh) which is smaller within the mesh motion case. __________________ Keep foaming, Tobias Holzmann

 July 23, 2015, 12:45 #4 Senior Member     Kyle Mooney Join Date: Jul 2009 Location: Amherst, MA USA - San Diego, CA USA Posts: 321 Rep Power: 11 I would say that the moving mesh case has less diffusion because there is near zero relative advection taking place. Like you said, with (U-Urel)~=0, there really isn't much math going on to even allow for numerical diffusion.

 July 31, 2015, 04:18 less diffusive scheme #5 Senior Member   Fabian Roesler Join Date: Mar 2009 Location: Germany Posts: 210 Rep Power: 11 Right, I agree with kmooney. The relative advection in the mesh motion case is smaller and so the artificial diffusion is smaller as well. What schemes did you apply and have you tried with a less diffusive scheme? Cheers Fabian

 July 31, 2015, 05:18 #6 Super Moderator     Tobias Holzmann Join Date: Oct 2010 Location: Leoben (Austria) Posts: 1,947 Blog Entries: 6 Rep Power: 33 Hi Fabian, at the beginning I was wondering why its less diffusive but the reason for that is exactly what we are talking about. (U - Urel) ~ 0. I also tried with different schemes, second order, 1st order, not much difference. Thanks for all feedbacks, __________________ Keep foaming, Tobias Holzmann

 August 5, 2015, 09:51 #7 Member   Hannes Join Date: Apr 2009 Location: Schleswig, Germany Posts: 37 Rep Power: 10 __________________ FluidEngineeringSolutions

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post alfaruk CFX 14 March 17, 2017 07:08 Ganesh FLUENT 13 January 22, 2014 05:11 Doginal CFX 2 January 12, 2014 07:21 DarrenC ANSYS Meshing & Geometry 11 August 5, 2013 06:42 rbarrett CFX 8 June 30, 2011 13:22

All times are GMT -4. The time now is 18:07.