CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > CFX

CEL Code for Sinusoidal Movement of Piston

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 12, 2015, 07:48
Default CEL Code for Sinusoidal Movement of Piston
  #1
Member
 
Dr Gurubasavaraju
Join Date: Dec 2014
Location: Bengaluru India
Posts: 77
Rep Power: 10
gbrajtm is on a distinguished road
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
gbrajtm is offline   Reply With Quote

Old   January 12, 2015, 12:16
Default
  #2
Senior Member
 
JuPa's Avatar
 
Mr CFD
Join Date: Jun 2012
Location: Britain
Posts: 361
Rep Power: 13
JuPa is on a distinguished road
Your equation will have to be a function of time which means transient runs.
JuPa is offline   Reply With Quote

Old   January 13, 2015, 00:42
Default
  #3
Member
 
Dr Gurubasavaraju
Join Date: Dec 2014
Location: Bengaluru India
Posts: 77
Rep Power: 10
gbrajtm is on a distinguished road
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
gbrajtm is offline   Reply With Quote

Old   January 18, 2015, 06:51
Default
  #4
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,261
Rep Power: 136
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
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.
ghorrocks is offline   Reply With Quote

Old   January 18, 2015, 23:22
Default
  #5
Member
 
Dr Gurubasavaraju
Join Date: Dec 2014
Location: Bengaluru India
Posts: 77
Rep Power: 10
gbrajtm is on a distinguished road
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
gbrajtm is offline   Reply With Quote

Old   January 18, 2015, 23:29
Default
  #6
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,261
Rep Power: 136
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
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.
ghorrocks is offline   Reply With Quote

Old   January 19, 2015, 00:41
Default
  #7
Member
 
Dr Gurubasavaraju
Join Date: Dec 2014
Location: Bengaluru India
Posts: 77
Rep Power: 10
gbrajtm is on a distinguished road
Yes its magnetorheological fluid damper.
Trying to get force velocity distribution.
gbrajtm is offline   Reply With Quote

Old   January 19, 2015, 17:19
Default
  #8
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,261
Rep Power: 136
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
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?
ghorrocks is offline   Reply With Quote

Old   January 20, 2015, 00:02
Default
  #9
Member
 
Dr Gurubasavaraju
Join Date: Dec 2014
Location: Bengaluru India
Posts: 77
Rep Power: 10
gbrajtm is on a distinguished road
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.
gbrajtm is offline   Reply With Quote

Old   January 20, 2015, 02:03
Default
  #10
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,261
Rep Power: 136
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
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.
ghorrocks is offline   Reply With Quote

Old   January 20, 2015, 02:33
Default
  #11
Member
 
Dr Gurubasavaraju
Join Date: Dec 2014
Location: Bengaluru India
Posts: 77
Rep Power: 10
gbrajtm is on a distinguished road
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
gbrajtm is offline   Reply With Quote

Old   January 20, 2015, 05:25
Default
  #12
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,261
Rep Power: 136
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
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.
ghorrocks is offline   Reply With Quote

Old   January 20, 2015, 05:30
Default
  #13
Member
 
Dr Gurubasavaraju
Join Date: Dec 2014
Location: Bengaluru India
Posts: 77
Rep Power: 10
gbrajtm is on a distinguished road
Thank you very much for your time.

Kind Regards
Guru
gbrajtm is offline   Reply With Quote

Old   February 23, 2015, 08:08
Default
  #14
Member
 
Dr Gurubasavaraju
Join Date: Dec 2014
Location: Bengaluru India
Posts: 77
Rep Power: 10
gbrajtm is on a distinguished road
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.
gbrajtm is offline   Reply With Quote

Old   February 24, 2015, 02:05
Default
  #15
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,261
Rep Power: 136
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
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.
ghorrocks is offline   Reply With Quote

Old   March 19, 2015, 00:34
Default
  #16
Member
 
Dr Gurubasavaraju
Join Date: Dec 2014
Location: Bengaluru India
Posts: 77
Rep Power: 10
gbrajtm is on a distinguished road
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
gbrajtm is offline   Reply With Quote

Old   March 19, 2015, 02:29
Default
  #17
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,261
Rep Power: 136
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
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.
ghorrocks is offline   Reply With Quote

Old   March 19, 2015, 06:57
Default
  #18
Member
 
Dr Gurubasavaraju
Join Date: Dec 2014
Location: Bengaluru India
Posts: 77
Rep Power: 10
gbrajtm is on a distinguished road
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
Attached Images
File Type: png Capture.PNG (42.6 KB, 43 views)
gbrajtm is offline   Reply With Quote

Old   March 19, 2015, 07:06
Default
  #19
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,261
Rep Power: 136
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
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.
ghorrocks is offline   Reply With Quote

Old   March 19, 2015, 07:12
Default
  #20
Member
 
Dr Gurubasavaraju
Join Date: Dec 2014
Location: Bengaluru India
Posts: 77
Rep Power: 10
gbrajtm is on a distinguished road
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.?
gbrajtm is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 09:32.