
[Sponsors] 
June 21, 2011, 05:55 
Prescribed mesh motion does not correspond to actual motion

#1 
New Member
Richard Barrett
Join Date: Jun 2011
Posts: 14
Rep Power: 6 
Hi, I am attempting to simulate the compression of a gas within a pistoncylinder apparatus using Ansys CFX. I am using CEL expressions to specify the mesh motion, coupled together with the 'Specified Displacement' option. I am applying the mesh motion to a subdomain of the main fluid domain and I am finding that the results of the simulation do not match the prescribed motion that I have applied through the CEL expressions, i.e. the CFX simulation lags behind the actual.
I was just wondering if it is correct to apply this type of mesh motion in this manner and if so, why do my results lag behind the prescribed motion. Thanks in advance Regards Richard 

June 21, 2011, 08:34 

#2 
Senior Member
Join Date: Apr 2009
Posts: 530
Rep Power: 13 
This means you've set something up wrong. Could you post the expressions used and the relevant CCL for the subdomain and boundary conditions?


June 21, 2011, 09:26 

#3 
New Member
Richard Barrett
Join Date: Jun 2011
Posts: 14
Rep Power: 6 
Thanks a million for your reply and apologies for the length of my reply.
I am attempting to model the compression of two opposed pistons in one cylinder, hence there is a line of symmetry down the centre of the combustion chamber and the model is axisymmetric about the central axis. Initially the cylinder is 178 mm and the stroke length is 167.81 mm. Included below are the CEL expressions used in this simulation. Note that anything to do with 'Function' are step functions used to define different phases of the compression, i.e. acceleration, CV and deceleration stages and I have omitted these CEL expressions. Also included is an image of the desired displacement function that I wish to apply as the mesh motion. This is applied to a subdomain which comprises of the full fluid domain and the function drops off linearly to zero at the mid chamber location (as there is zero mesh motion at this location). The BC's used in this simulation are all either wall or symmetry BC's, all with 'unspecified' mesh motion as mesh motion is dependent on the motion of the subdomain. Disp_Plot_CFX.PNG Disp = OverallFunction*sDisp TDC = 167.81e03 [m] TotalTime = 116.6e03 [s] aAccel = 2079.05 [m/s^2] aDecel = 23493.3 [m/s^2] sDisp = sDispAccel+sDispCV+sDispDecel+sDispZero sDispAccel = (0.5*aAccel*t*t)*Function1 sDispCV = ((velCV*t)(0.5*aAccel*tAccel*tAccel))*Function2 sDispDecel = (((ttimeDecel)*(velCV))+(0.5*aDecel*(ttimeDecel)*(ttimeDecel))+((velCV*timeDecel)(0.5*aAccel*tAccel*tAccel)))*Function3 sDispZero = TDC*Function4 tAccel = 6.07055e03 [s] tCV = 9.992233e03 [s] tTotal = 16.6e03 [s] timeDecel = tAccel+tCV velCV = 12.621 [m/s] 

June 21, 2011, 09:42 

#4 
Senior Member
Join Date: Apr 2009
Posts: 530
Rep Power: 13 
So 'Disp' is the function applied to the subdomain, with 'OverallFunction' providing the linear weighting? Since this appears to only depend on t (and x/y/z to get the linear weighting I assume), then the results should be fairly easy to debug. In your results file you can get the "Total Mesh Displacement" at one end of the cylinder. Given t and x/y/z at that location you can evaluate your displacement expression. Are you saying these two don't match?


June 21, 2011, 10:07 

#5 
New Member
Richard Barrett
Join Date: Jun 2011
Posts: 14
Rep Power: 6 
Thanks for your reply and yes, 'Disp' is the desired displacement function and applied as the mesh motion to the subdomain. As is shown in the graph above, the applied total displacement is 167.81 mm at the end of the cylinder exposed to the maximum displacement. From the results file, the 'Total Mesh Displacement' at this end of the cylinder is only 91.72 mm.


June 21, 2011, 13:48 

#6 
Senior Member
Join Date: Apr 2009
Posts: 530
Rep Power: 13 
What is OverallFunction equal to? If it depends on x/y/z did you account for the fact that x/y/z will change during the solution as the mesh moves?


June 21, 2011, 15:26 

#7 
New Member
Richard Barrett
Join Date: Jun 2011
Posts: 14
Rep Power: 6 
The two variables that every CEL expression depends on are 't' and 'aaxis'.... Does the 'aaxis' variable change as mesh deformation occurs?


June 21, 2011, 16:35 

#8 
Senior Member
Join Date: Apr 2009
Posts: 530
Rep Power: 13 
aaxis is usually just equal to z, unless you have some local coord frame defined. So if the mesh moves in the z direction then yes, aaxis will change.


June 30, 2011, 13:22 

#9 
New Member
Richard Barrett
Join Date: Jun 2011
Posts: 14
Rep Power: 6 
Thanks a million for your help.


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Negative Volume during Mesh Motion Analysis  giov_ingr  FLUENT  2  December 13, 2013 07:09 
Dynamic moving mesh  PeiYing Hsieh (Hsieh)  OpenFOAM Running, Solving & CFD  64  June 7, 2012 10:04 
a mesh motion question  womo  CFX  5  March 16, 2011 11:54 
external flow with snappyHexMesh  chelvistero  OpenFOAM  11  January 15, 2010 20:43 
mesh motion  samad87  FLUENT  0  August 6, 2009 03:15 