# Propeller disc theory in Star CCM+

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

 May 4, 2012, 09:59 Propeller disc theory in Star CCM+ #1 New Member   Christian Rodriguez Join Date: Jan 2012 Posts: 9 Rep Power: 5 Hey guys, I am currently doing my master thesis in rotor aerodynamics for a small scaled helicopter. The first step before I include the main rotors in my mesh is to make a simplified verison of the model. My idea is to create a infinite thin disc over my helicopter body to simulate the flow of the rotor blade. The model will be simulated in a box which acts as a windtunnel (with a certain freestream velocity). Is there a easy way to model a disc in Star CCM+? I was wondering about this due to the fact that adding the disc will create a discontinuity in velocity and pressure around the rotor. Great forum btw //Christian

 May 6, 2012, 09:02 #2 Senior Member   Join Date: Oct 2009 Location: Germany Posts: 637 Rep Power: 12 I'm not too sure about an infinite thin disk. Maybe it works with a fan interface, but I haven't tried that. Another option would be to create a disk with a finite thickness. Create a region from that and put a momentum source. That should work for sure (we did something similar for wind park simulations to simplify the turbines). __________________ We do three types of jobs here: GOOD, FAST AND CHEAP You may choose any two!

 May 7, 2012, 08:24 #3 Member   Ryan Coe Join Date: Jun 2010 Location: Albuquerque, NM Posts: 98 Rep Power: 7 I've used the finite thickness method described by abdul099 to model marine propellers. You can use a field function to set the strength of the momentum source based on some known propeller curve (ie advance ratio vs thrust and torque output). __________________ Ryan

December 13, 2014, 22:07
#4
Member

Join Date: Nov 2014
Posts: 79
Rep Power: 2
I am new to CCM+ and have some problems implementing an actuator disc to simulate a propellor/wind turbine. Need some help here please as i have spent very long hours but not getting anywhere.

I have seen many of the given tutorials and I thought the following tutorials are relevant to what I wanna achieve.

Quote:
 1. Motion \ Marine Resistance Prediction: KCS Hull with a Rudder and 2. Motion \ Moving Reference Frames: Marine Propellor in Open Water and 3. Incompressible Flow \ Porous Resistance: Isotropic Media
My current situation
Quote:
 I have managed to simulate a a solid disc inside a wind tunnel test section. I used Mesh Models described in 1. KCS_Geometry Tutorial, SURFACE REMESHER | TRIMMER | PRISM LAYER MESHER. Reference Values | Base size was set as the my Disc's diameter size. Other than that, Mesh Default mesh settings was DEFAULT as my whole mesh got worse after I started playing around with the settings. All the sides of the wind tunnels were set as WALLS as they should. At this stage, the disc was set to WALLS on the inlet, outlet and sides. Solver was IMPLICIT UNSTEADY. Inlet Region was 10 m/s with turbulent intensity 2% and turbulent length scale 0.01. actuatordisc.jpg These are my results (pardon poor pics, I can't find a straighforward way in CCM+ to output pictures nicely). Velocity Scene. velocity.jpg Results don't look well especially at the wall of the wind tunnel and the disc, I would expect lower velocity due to the no-slip condition. I still am abit lost on how I should tweak my mesh settings. Any ideas please. Vorticity Scene. vorticity.png EDIT: My base size was previously set to the size of my disc. I guess that was too big, which let to large mesh grids. I now set the base to 1/2 of my diameter. Mesh looks better now, but my solution doesn't converge. I guess it must be something to do with the Courant number. So, I set my time step from 0.02 to 0.01 and I have convergence now. Any wisdow on how to how low my base size or timestep, or is it a trial and error thing? I wondering how low should I go to get a fine solution at the walls.

Quote:
 Now I want to implement a a pressure drop across the disc. Based on what I saw from the tutorial 3. Incompressible Flow \ Porous Resistance: Isotropic Media, I believe I need to use define a Porous region for my disc porous.jpg In my mind, I have a self-determined thrust coefficient (C_T) of my disc and then I can apply T = C_T * 0.5* 1/2 * rho * Vinf^2 * A. Since I am assuming uniform thrust across my disc, my thrust per unit area would be t = T/A = C_T * 0.5* 1/2 * rho * Vinf^2. But now, I am a bit lost on how I should implement this field function. spent fustrating hours and hours and getting nowhere. Also, when I tried to regenerate the mesh under this case, an error occurs. mesherror.jpg I believe my error occured because the way the region has been defined. To define 2 regions (1 static region - flow outside ) and (1 rotating region - disc), I had to apply "ASSIGN PARTS TO REGION" to the PARTS on (1) Disc and (2) Wind Tunnel separately. However, this resulted in the disc REGION being inside the Wind Tunnel REGION, and which I suspect, caused this error. how should I go about this defining my disc in this manner? I couldn't find any specific tutorial that teaches this.

your guidance and help in part or whole is really really appreciated and would really make this whole work more enjoyable.

thanks!

Last edited by hwsv07; December 14, 2014 at 08:33.

December 14, 2014, 16:15
#5
Member

Join Date: Nov 2014
Posts: 79
Rep Power: 2
In my mind, to create an Actuator Disc simulation in a wind tunnel, the aim is to create 2 regions:
a) Rotating Region (representing the Disc)
b) Fluid Region (representing the wind coming onto the Disc)

Quote:
 I have decided to adopt a Rigid Body Frame just like in the tutorial Motion \ Rigid Body Motion: Rotating Fan. This means, the disc will rotate. So I created a DISC in a TEST SECTION in 3D-CAD. setup.jpg a) Since they are essentially 2 PARTS (disc and wind tunnel box) intersecting each other, I need to SUBTRACT the volume of the DISC from the TEST SECTION. b) Since I need to specify interaction between the rotating disc and the rest of the fluid, I need to use the INTERSECT operation and obtain the intersecting surfaces of the disc (i.e. the upstream surface, downstream surface and the sides).
Quote:
 With PARTS: SUBTRACT and INTERSECT created, I proceeded to assign regions (i.e. 2 regions)to these parts. I also created INTERFACE relationship on those disc surfaces.
So I thought I was right up to this point. Then I wanted to proceed to generate the VOLUME MESH but got this error:

Code:
`Trimmer cannot currently accommodate multiple regions in one continuum.`

I really do not understand why.

December 15, 2014, 17:23
#6
Member

Join Date: Nov 2014
Posts: 79
Rep Power: 2
Ok so I made a bit of progress after much fiddling around over and over.

Quote:
 I managed to get to past the problem with defining separate regions for Rotating Fluid and Static Fluid. I also implemented a rotating disc based on a Moving Reference Frame (MRF) approach and specified some rpm value for my disc.
Quote:
 DiscInlet_Velocity.jpg Here you can see the velocity of the fluid at the disc upstream surface. My disc is rotating clockwise, so I thought the velocity field should be flipped the other way around.
Quote:
 Pressure Drop.jpg Across the disc, the expected pressure drop is observed.

Moving ahead, I have now specified a porous region for the disc.

Under porosity option, there were a few options of inputing porosity: Constant, Field Function etc...

I tried an extreme case of porosity 0 as a constant value. So I expected the flow to just pass through my disc as if it was non-existent. But this did not happen and there was absolutely no change in my thrust coefficient.

How should I specify a specifc thrust for my disc?

Thank you!

 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 TAREK GANAT STAR-CCM+ 1 May 21, 2013 22:15 clo OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 33 September 26, 2012 04:04 Kart OpenFOAM Meshing & Mesh Conversion 1 February 4, 2010 05:38 maurizio CD-adapco 3 October 16, 2007 05:17 michele OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 2 July 15, 2005 04:15

All times are GMT -4. The time now is 16:21.