CEL Code for Sinusoidal Movement of Piston

ghorrocks · March 19, 2015, 07:17

I have no idea what you are modelling so I cannot say whether it is correct. You also have not shown how you implemented anything.

The equations you describe are the Bingham model, but as I said I would try to use the existing built in Bingham model if possible - but if you want to make yield stress a function of something I suspect you may have to define the model yourself.

gbrajtm · March 19, 2015, 07:29

I can mail you my simulation if you provide the email id.

I am trying simulate a magneto rheological fluid damper to study the performance by obtaining Force developed in the damper.

The small thin gap is the annular gap of the piston where the fluid get magnetised under motion with certain velocity.

ghorrocks · March 19, 2015, 17:36

I do not have time to go through your simulation in detail.

The forum works best when you identify a specific problem you are having and post a question about that problem. Make sure you include enough information so we can understand what you are talking about and its context.

gbrajtm · March 31, 2015, 06:50

Hi Glenn Horrocks,

Greetings,Hope you doing well.

I have added viscosity model as function of magnetic flux density using CEL, But i have not considered the piston location while calculating the Dynamic viscosity.
This the Expression i have used.

if(Shear Strain Rate >= 15[s^-1], (5000[kg m^-1 s^-2]/Shear Strain Rate)+0.11[kg m^-1 s^-1], 55[kg m^-1 s^-1])

I am finding difficulties in considering the Piston location as it moves up and down.

Do you have any solution for this.

Kind Regards
Guru

ghorrocks · March 31, 2015, 08:00

I would expect convergence difficulties with a viscosity function like that. It has a big step in viscosity at SSR=15[s^-1]. Are you sure it is realistic?

Looking back on this thread you state that you are defining the piston motion. So why can't you simply use the variable you defined for the piston motion as an input to this equation as well?

gbrajtm · March 31, 2015, 08:14

yes, the value provided are not realistic. The values are different, for just your ference i have added.

How to add the piston motion equation as the input to the Viscosity expression.

ghorrocks · March 31, 2015, 18:23

You just include the x, y or z displacement of the piston in the CEL expression - it is not that hard. You just need to be careful about getting the units right.

gbrajtm · April 1, 2015, 03:30

Sorry for troubling you. I didn't get how to add.

My displacement expression for Piston motion in X direction is: 0.0425[m]*sin(t/1 [s])

And dynamic viscosity is

if(Shear Strain Rate >= 15[s^-1], (5000[kg m^-1 s^-2]/Shear Strain Rate)+0.11[kg m^-1 s^-1], 55[kg m^-1 s^-1])

ghorrocks · April 1, 2015, 08:03

How about this:

PistonX = 0.0425[m]*sin(t/1 [s])
Tolerance = 0.1 [m]
CloseToPiston = if(abs(x-PistonX)<Tolerance,1,0)
Viscosity = if(Shear Strain Rate >= 15[s^-1], (5000[kg m^-1 s^-2]/Shear Strain Rate)+0.11[kg m^-1 s^-1]+CloseToPiston*10[Pa s], 55[kg m^-1 s^-1] + CloseToPiston*10[Pa s])

You might need to fix the units up but hopefully you get the idea.

This function adds 10 Pa.s to the viscosity when you are within 0.1 metres of the current piston location.

gbrajtm · April 1, 2015, 08:46

Thank you very much i will simulate and let you.

Thanks

gbrajtm · April 11, 2015, 05:39

Hi,
Glenn Horrocks,
The expression below are used in the simulation.

St = 10000*tanh(Shear Strain Rate/0.001[s^-1])
pistonx = 0.0425[m]*cos(1.33*t/1 [s])
Vel = 0.06[m/s]*sin(1.33*t/1 [s])
Tolerance = 0.01 [m]
CloseToPiston = if(abs(x-PistonX)<Tolerance,1,0)
Visco = if(Shear Strain Rate >= 0.001[s^-1], (St[kg m^-1 s^-2]/Shear Strain Rate)+0.11[kg m^-1 s^-1]+CloseToPiston*10[Pa s], 100[kg m^-1 s^-1] + CloseToPiston*10[Pa s])

I am getting error(shows as Overflow and negative volume element) in solver if i use total time 0.5s and time step 0.01s. if i change it 0.1s and 0.005 time steps, simulation goes fine. any suggestions on this.

Please find the attachment for details of my simulation

Second is can you please elaborate. CloseToPiston = if(abs(x-PistonX)<Tolerance,1,0)

Thanks in advance.

ghorrocks · April 12, 2015, 08:52

Fix the negative volume element issue first. It might the cause of the overflow error as well.

Rerun the simulation generating transient results files on every time step, or at least the time steps just before it crashes. Make sure you include the mesh in the transient results files. Then use the post processor to look at the mesh and work out where the mesh is turning in side out. Then you know what you need to fix.

gbrajtm · April 20, 2015, 08:09

If want to consider the piston location along with other two co-ordinate, what changes i have to make to this expression.
CloseToPiston = if(abs(x-PistonX)<Tolerance,1,0)

ghorrocks · April 20, 2015, 08:40

Isn't that high school level maths? - http://www.purplemath.com/modules/distform.htm

gbrajtm · April 23, 2015, 03:16

Yes , I know the expressions, But didn't know how to make it compatible with ANSYS CEL.

My intention is the thin region should behave like Non newtonian as per the expression given, but rest region is should be newtonian.
Presently i have used Inside()@location for specifying the region.

if(Shear Strain Rate >= 0.001[s^-1], (40000[kg m^-1 s^-2]/Shear Strain Rate)+0.11[kg m^-1 s^-1], 55[kg m^-1 s^-1])*inside()@piston

is this correct?

But the simulation gives error for all turbulence models and runs without error for Laminar model.

ghorrocks · April 23, 2015, 06:53

Your function returns 0 if you are outside the region. You cannot set the viscosity to be zero.

gbrajtm · April 23, 2015, 07:13

thanks,
now I got it, so results are as not expected.
Closetopiston expression used also not service the goal.
is there any other method? where i can write expression depending on time in the same domain, while moving the piston in the domain, with specified CCL expression whether node is active region or not.

ghorrocks · April 23, 2015, 08:18

To make your latest CEL expression work you simply need to make it return a non-zero value when outside the region. A simple if statement will do that:

a = if(Shear Strain Rate >= 0.001[s^-1], (40000[kg m^-1 s^-2]/Shear Strain Rate)+0.11[kg m^-1 s^-1], 55[kg m^-1 s^-1])*inside()@piston
viscosity = if(inside@piston == 1,a, 1 [kg m^-1 s^-1])

Or something like that.

gbrajtm · May 6, 2015, 05:00

Hi Glenn Horrocks,
This equations works fine, even i am getting some results. Then i noticed that the @piston is 2D location as shown in the fig. But i am not sure whether " Inside" function take only the surface into account or the region below it along.
I was trying to add the 3D region as marked in the red. but it does not happening. There is no option for getting cells location which are in the thin region.

ghorrocks · May 6, 2015, 20:14

For a 3D region you need to define it in your mesh. So you will need to set this up in the meshing software.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Movement of piston in a cylinder with fluent?	Whitworth	ANSYS	5	March 4, 2015 00:50
Piston Simulation	helloworld922	Autodesk Simulation CFD	1	September 9, 2012 23:54
State of the art in CFD technology	Juan Carlos GARCIA SALAS	Main CFD Forum	39	November 1, 1999 15:34
public CFD Code development	Heinz Wilkening	Main CFD Forum	38	March 5, 1999 12:44
What kind of Cmmercial CFD code you feel well?	Lans	Main CFD Forum	13	October 27, 1998 11:20

March 19, 2015, 07:17		#21
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,854 Rep Power: 144	I have no idea what you are modelling so I cannot say whether it is correct. You also have not shown how you implemented anything. The equations you describe are the Bingham model, but as I said I would try to use the existing built in Bingham model if possible - but if you want to make yield stress a function of something I suspect you may have to define the model yourself.

March 19, 2015, 07:29		#22
gbrajtm Member Dr Gurubasavaraju Join Date: Dec 2014 Location: Bengaluru India Posts: 78 Rep Power: 11	I can mail you my simulation if you provide the email id. I am trying simulate a magneto rheological fluid damper to study the performance by obtaining Force developed in the damper. The small thin gap is the annular gap of the piston where the fluid get magnetised under motion with certain velocity.

March 19, 2015, 17:36		#23
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,854 Rep Power: 144	I do not have time to go through your simulation in detail. The forum works best when you identify a specific problem you are having and post a question about that problem. Make sure you include enough information so we can understand what you are talking about and its context.

March 31, 2015, 08:00		#25
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,854 Rep Power: 144	I would expect convergence difficulties with a viscosity function like that. It has a big step in viscosity at SSR=15[s^-1]. Are you sure it is realistic? Looking back on this thread you state that you are defining the piston motion. So why can't you simply use the variable you defined for the piston motion as an input to this equation as well?

March 31, 2015, 08:14		#26
gbrajtm Member Dr Gurubasavaraju Join Date: Dec 2014 Location: Bengaluru India Posts: 78 Rep Power: 11	yes, the value provided are not realistic. The values are different, for just your ference i have added. How to add the piston motion equation as the input to the Viscosity expression.

March 31, 2015, 18:23		#27
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,854 Rep Power: 144	You just include the x, y or z displacement of the piston in the CEL expression - it is not that hard. You just need to be careful about getting the units right.

April 1, 2015, 03:30		#28
gbrajtm Member Dr Gurubasavaraju Join Date: Dec 2014 Location: Bengaluru India Posts: 78 Rep Power: 11	Sorry for troubling you. I didn't get how to add. My displacement expression for Piston motion in X direction is: 0.0425[m]*sin(t/1 [s]) And dynamic viscosity is if(Shear Strain Rate >= 15[s^-1], (5000[kg m^-1 s^-2]/Shear Strain Rate)+0.11[kg m^-1 s^-1], 55[kg m^-1 s^-1])

April 1, 2015, 08:03		#29
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,854 Rep Power: 144	How about this: PistonX = 0.0425[m]sin(t/1 [s]) Tolerance = 0.1 [m] CloseToPiston = if(abs(x-PistonX)<Tolerance,1,0) Viscosity = if(Shear Strain Rate >= 15[s^-1], (5000[kg m^-1 s^-2]/Shear Strain Rate)+0.11[kg m^-1 s^-1]+CloseToPiston10[Pa s], 55[kg m^-1 s^-1] + CloseToPiston*10[Pa s]) You might need to fix the units up but hopefully you get the idea. This function adds 10 Pa.s to the viscosity when you are within 0.1 metres of the current piston location.

April 1, 2015, 08:46		#30
gbrajtm Member Dr Gurubasavaraju Join Date: Dec 2014 Location: Bengaluru India Posts: 78 Rep Power: 11	Thank you very much i will simulate and let you. Thanks

April 12, 2015, 08:52		#32
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,854 Rep Power: 144	Fix the negative volume element issue first. It might the cause of the overflow error as well. Rerun the simulation generating transient results files on every time step, or at least the time steps just before it crashes. Make sure you include the mesh in the transient results files. Then use the post processor to look at the mesh and work out where the mesh is turning in side out. Then you know what you need to fix.

April 20, 2015, 08:09		#33
gbrajtm Member Dr Gurubasavaraju Join Date: Dec 2014 Location: Bengaluru India Posts: 78 Rep Power: 11	If want to consider the piston location along with other two co-ordinate, what changes i have to make to this expression. CloseToPiston = if(abs(x-PistonX)<Tolerance,1,0)

April 20, 2015, 08:40		#34
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,854 Rep Power: 144	Isn't that high school level maths? - http://www.purplemath.com/modules/distform.htm

April 23, 2015, 03:16		#35
gbrajtm Member Dr Gurubasavaraju Join Date: Dec 2014 Location: Bengaluru India Posts: 78 Rep Power: 11	Yes , I know the expressions, But didn't know how to make it compatible with ANSYS CEL. My intention is the thin region should behave like Non newtonian as per the expression given, but rest region is should be newtonian. Presently i have used Inside()@location for specifying the region. if(Shear Strain Rate >= 0.001[s^-1], (40000[kg m^-1 s^-2]/Shear Strain Rate)+0.11[kg m^-1 s^-1], 55[kg m^-1 s^-1])*inside()@piston is this correct? But the simulation gives error for all turbulence models and runs without error for Laminar model.

April 23, 2015, 06:53		#36
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,854 Rep Power: 144	Your function returns 0 if you are outside the region. You cannot set the viscosity to be zero.

April 23, 2015, 07:13		#37
gbrajtm Member Dr Gurubasavaraju Join Date: Dec 2014 Location: Bengaluru India Posts: 78 Rep Power: 11	thanks, now I got it, so results are as not expected. Closetopiston expression used also not service the goal. is there any other method? where i can write expression depending on time in the same domain, while moving the piston in the domain, with specified CCL expression whether node is active region or not.

April 23, 2015, 08:18		#38
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,854 Rep Power: 144	To make your latest CEL expression work you simply need to make it return a non-zero value when outside the region. A simple if statement will do that: a = if(Shear Strain Rate >= 0.001[s^-1], (40000[kg m^-1 s^-2]/Shear Strain Rate)+0.11[kg m^-1 s^-1], 55[kg m^-1 s^-1])*inside()@piston viscosity = if(inside@piston == 1,a, 1 [kg m^-1 s^-1]) Or something like that.

May 6, 2015, 05:00		#39
gbrajtm Member Dr Gurubasavaraju Join Date: Dec 2014 Location: Bengaluru India Posts: 78 Rep Power: 11	Hi Glenn Horrocks, This equations works fine, even i am getting some results. Then i noticed that the @piston is 2D location as shown in the fig. But i am not sure whether " Inside" function take only the surface into account or the region below it along. I was trying to add the 3D region as marked in the red. but it does not happening. There is no option for getting cells location which are in the thin region.

May 6, 2015, 20:14		#40
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,854 Rep Power: 144	For a 3D region you need to define it in your mesh. So you will need to set this up in the meshing software.