CEL Code for Sinusoidal Movement of Piston

gbrajtm · January 12, 2015, 07:48

Hi,

How to write CEL expression for Sinusoidal movement of piston in cylinder.

How make different fluid regions in the global domain using CCL expressions.

Kind Regard's
Guru

JuPa · January 12, 2015, 12:16

Your equation will have to be a function of time which means transient runs.

gbrajtm · January 13, 2015, 00:42

Thanks for your time,

I know the equation for motion, ie
Stroke S=S(max) COS (wt)
Velocity U=U(max) SIN(wt)

But i dont know how to implement it in CFX,
Presently i am using Periodic motion option in rigid body.

Could please help in this,

Kind regard's
Guru

ghorrocks · January 18, 2015, 06:51

You do not need rigid bodies to model a piston in a cylinder, just moving mesh.

Also: If you are modelling a standard conrod driven piston (like an IC engine) then your equation of motion is wrong. The conrod introduces higher order terms into the motion which are significant in almost all cases.

gbrajtm · January 18, 2015, 23:22

Hi Glenn,

Thanks for ur respond,

I am simulating a MR damper, where the fluid near piston is non newtonian and rest region is newtonian fluid.
The piston moves similar to shock absorber.
I am facing problem in specifying the Fluid domain and piston motion.

Please let me know of you have any suggestions.

Kind regards
Guru

ghorrocks · January 18, 2015, 23:29

Please do not use jargon. What is a MR damper? .... I just looked it up on google, is it a Magnetorheological damper?

What are you trying to learn with this simulation? What are you trying to do? The type of model you do depends on what you are trying to achieve.

gbrajtm · January 19, 2015, 00:41

Yes its magnetorheological fluid damper.
Trying to get force velocity distribution.

ghorrocks · January 19, 2015, 17:19

What do you mean by force velocity distribution? How does that fit in with your original request of how to put a sinusoidal motion on the piston?

gbrajtm · January 20, 2015, 00:02

I wanted to study the force developed(pressure) at different velocity of the piston, when it is under reciprocatory action.

The viscosity of fluid going to change near the piston, when it is exposed to magnetic field.

ghorrocks · January 20, 2015, 02:03

OK, thanks.

In that case you definitely do not want to use the rigid body solver. Use mesh motion and put your equation of motion on the moving piston. Put unspecified motion on the side boundaries and stationary boundaries on the top face. This will be the simplest way to implement this simulation.

gbrajtm · January 20, 2015, 02:33

Presently i am using periodic displacement in the mesh motion. But if i wanted use the expression where exactly i supposed to put it.

Kind regards
Guru

ghorrocks · January 20, 2015, 05:25

There used to be example tutorials of pure mesh motion but they don't seem to be in V15. Have a look on the ANSYS community website, there might be some examples there.

gbrajtm · January 20, 2015, 05:30

Thank you very much for your time.

Kind Regards
Guru

gbrajtm · February 23, 2015, 08:08

Hi Glenn Horrocks,

Hope you doing well.

I need one more help, my problem statement is below.

I have magnetorheological damper, Where the area near to the piston (piston has holes to pass fluid) will become non newtonian, when it gets magnetised fluid either side of it will be remain newtonian. I wanted to move the piston with certain stroke and velocity under the said condition.

I don't know how to write an CCL command language for the fluid property.

could you please help.

ghorrocks · February 24, 2015, 02:05

The CEL expressions to do material properties are similar to any other expressions (for example your piston motion equations). Have a look in the tutorials for examples of CEL expressions.

gbrajtm · March 19, 2015, 00:34

Hi Glenn,

I got the expression where the Bingham model is defined in terms of shear rate. How can i incorporate it into fix?

Equations are

μ= τ_y/γ + μp γ ≥ 15
μ = μs γ < 15

μs ≥ 500 * μp

ghorrocks · March 19, 2015, 02:29

The Bingham non-Newtonian model is already a built in model in CFX. There is no need to write the equations to define it. Just use the existing model.

The built in model also has the advantage that it is properly integrated into the solver, so derivatives are correctly calculated for linearisation to assist convergence and error checking is done. You won't get this if you define your own viscosity function.

gbrajtm · March 19, 2015, 06:57

In same computational domain i have two types of fluid properties. In annular gap its not newtonian fluid due to viscosity and in rest other region it is newtonian. so viscosity is defined in terms of yield stress(function of magnetic flux density).
Hence wherever the shear rate(exceeds the specified value) changes that is non newtonian region and rest will be newtonian, this is according to the expression i have mentioned earlier.

Please have look on the model, which i have attached.

Regard's
Guru

ghorrocks · March 19, 2015, 07:06

The function you describe is a simple CEL expression. A couple of if functions and some other basic operators should do it. Have a look in the CFX reference guide for the available functions and CEL syntax. The tutorials give examples of several CEL expressions.

The model you propose sounds complex - a non-Newtonian fluid linked to a magnetic flux density field. If you are going to do a simulation as complex as this you will need to be able to do basic stuff like CEL expressions yourself.

gbrajtm · March 19, 2015, 07:12

Yes its little complex, I have obtained the yield stress value from the magnetic flux density.

Is the way i am doing is correct? may i know what you think.?

January 12, 2015, 07:48	CEL Code for Sinusoidal Movement of Piston	#1
gbrajtm Member Dr Gurubasavaraju Join Date: Dec 2014 Location: Bengaluru India Posts: 78 Rep Power: 11	Hi, How to write CEL expression for Sinusoidal movement of piston in cylinder. How make different fluid regions in the global domain using CCL expressions. Kind Regard's Guru

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

January 12, 2015, 12:16		#2
JuPa Senior Member Mr CFD Join Date: Jun 2012 Location: Britain Posts: 361 Rep Power: 15	Your equation will have to be a function of time which means transient runs.

January 13, 2015, 00:42		#3
gbrajtm Member Dr Gurubasavaraju Join Date: Dec 2014 Location: Bengaluru India Posts: 78 Rep Power: 11	Thanks for your time, I know the equation for motion, ie Stroke S=S(max) COS (wt) Velocity U=U(max) SIN(wt) But i dont know how to implement it in CFX, Presently i am using Periodic motion option in rigid body. Could please help in this, Kind regard's Guru

January 18, 2015, 06:51		#4
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,854 Rep Power: 144	You do not need rigid bodies to model a piston in a cylinder, just moving mesh. Also: If you are modelling a standard conrod driven piston (like an IC engine) then your equation of motion is wrong. The conrod introduces higher order terms into the motion which are significant in almost all cases.

January 18, 2015, 23:22		#5
gbrajtm Member Dr Gurubasavaraju Join Date: Dec 2014 Location: Bengaluru India Posts: 78 Rep Power: 11	Hi Glenn, Thanks for ur respond, I am simulating a MR damper, where the fluid near piston is non newtonian and rest region is newtonian fluid. The piston moves similar to shock absorber. I am facing problem in specifying the Fluid domain and piston motion. Please let me know of you have any suggestions. Kind regards Guru

January 18, 2015, 23:29		#6
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,854 Rep Power: 144	Please do not use jargon. What is a MR damper? .... I just looked it up on google, is it a Magnetorheological damper? What are you trying to learn with this simulation? What are you trying to do? The type of model you do depends on what you are trying to achieve.

January 19, 2015, 00:41		#7
gbrajtm Member Dr Gurubasavaraju Join Date: Dec 2014 Location: Bengaluru India Posts: 78 Rep Power: 11	Yes its magnetorheological fluid damper. Trying to get force velocity distribution.

January 19, 2015, 17:19		#8
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,854 Rep Power: 144	What do you mean by force velocity distribution? How does that fit in with your original request of how to put a sinusoidal motion on the piston?

January 20, 2015, 00:02		#9
gbrajtm Member Dr Gurubasavaraju Join Date: Dec 2014 Location: Bengaluru India Posts: 78 Rep Power: 11	I wanted to study the force developed(pressure) at different velocity of the piston, when it is under reciprocatory action. The viscosity of fluid going to change near the piston, when it is exposed to magnetic field.

January 20, 2015, 02:03		#10
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,854 Rep Power: 144	OK, thanks. In that case you definitely do not want to use the rigid body solver. Use mesh motion and put your equation of motion on the moving piston. Put unspecified motion on the side boundaries and stationary boundaries on the top face. This will be the simplest way to implement this simulation.

January 20, 2015, 02:33		#11
gbrajtm Member Dr Gurubasavaraju Join Date: Dec 2014 Location: Bengaluru India Posts: 78 Rep Power: 11	Presently i am using periodic displacement in the mesh motion. But if i wanted use the expression where exactly i supposed to put it. Kind regards Guru

January 20, 2015, 05:25		#12
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,854 Rep Power: 144	There used to be example tutorials of pure mesh motion but they don't seem to be in V15. Have a look on the ANSYS community website, there might be some examples there.

January 20, 2015, 05:30		#13
gbrajtm Member Dr Gurubasavaraju Join Date: Dec 2014 Location: Bengaluru India Posts: 78 Rep Power: 11	Thank you very much for your time. Kind Regards Guru

February 23, 2015, 08:08		#14
gbrajtm Member Dr Gurubasavaraju Join Date: Dec 2014 Location: Bengaluru India Posts: 78 Rep Power: 11	Hi Glenn Horrocks, Hope you doing well. I need one more help, my problem statement is below. I have magnetorheological damper, Where the area near to the piston (piston has holes to pass fluid) will become non newtonian, when it gets magnetised fluid either side of it will be remain newtonian. I wanted to move the piston with certain stroke and velocity under the said condition. I don't know how to write an CCL command language for the fluid property. could you please help.

February 24, 2015, 02:05		#15
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,854 Rep Power: 144	The CEL expressions to do material properties are similar to any other expressions (for example your piston motion equations). Have a look in the tutorials for examples of CEL expressions.

March 19, 2015, 00:34		#16
gbrajtm Member Dr Gurubasavaraju Join Date: Dec 2014 Location: Bengaluru India Posts: 78 Rep Power: 11	Hi Glenn, I got the expression where the Bingham model is defined in terms of shear rate. How can i incorporate it into fix? Equations are μ= τ_y/γ + μp γ ≥ 15 μ = μs γ < 15 μs ≥ 500 * μp

March 19, 2015, 02:29		#17
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,854 Rep Power: 144	The Bingham non-Newtonian model is already a built in model in CFX. There is no need to write the equations to define it. Just use the existing model. The built in model also has the advantage that it is properly integrated into the solver, so derivatives are correctly calculated for linearisation to assist convergence and error checking is done. You won't get this if you define your own viscosity function.

March 19, 2015, 07:06		#19
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,854 Rep Power: 144	The function you describe is a simple CEL expression. A couple of if functions and some other basic operators should do it. Have a look in the CFX reference guide for the available functions and CEL syntax. The tutorials give examples of several CEL expressions. The model you propose sounds complex - a non-Newtonian fluid linked to a magnetic flux density field. If you are going to do a simulation as complex as this you will need to be able to do basic stuff like CEL expressions yourself.

March 19, 2015, 07:12		#20
gbrajtm Member Dr Gurubasavaraju Join Date: Dec 2014 Location: Bengaluru India Posts: 78 Rep Power: 11	Yes its little complex, I have obtained the yield stress value from the magnetic flux density. Is the way i am doing is correct? may i know what you think.?