|
[Sponsors] |
Propeller simulation not matching experimental data |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
May 6, 2024, 15:26 |
|
#41 |
Senior Member
Join Date: Jun 2009
Posts: 1,875
Rep Power: 33 |
Q: what is the source of the values you are trying to match?
Q: are they absolutely free of errors? Experimental data has error bands are there other "fundamental checks" you can do that can increase the confidence on your results w/o having to look at someone else data? Recall in practice, there will be no benchmark, and conclusions must drawn anyways. What kind of back of the envelope thought problems can be used to infere if the current data is valid or not?
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
May 6, 2024, 17:33 |
Experimental data
|
#42 | |
Member
Kevin Gnanaraj
Join Date: Apr 2024
Posts: 34
Rep Power: 2 |
Quote:
It will be verified somewhat soon though, in the university's wind tunnel. My simulation is a part of the larger project. On a side-note, I made the flow compressible, and it is improving results a little, but not nearly enough to make up the difference. Also, were those questions meant to be rhetorical? Last edited by keg504; May 6, 2024 at 17:34. Reason: Fixed writing, added link |
||
May 6, 2024, 19:56 |
|
#43 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,862
Rep Power: 144 |
At Re=3E5 the flow is likely to be almost all turbulent for most airfoil shapes - unless this is a laminar flow airfoil. That means using a turbulence model for the entire this is likely to be OK (again, unless it is a laminar flow airfoil).
Propellers often have separations near the hub as the blade is often quite fat in that region. This usually does not make much of a difference the blade near the hub does not contribute much to propeller performance. To visualise the separations you will want to look at the wall shear on the blades. Where it goes to zero or reverses direction you have a separation.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
May 7, 2024, 10:21 |
|
#44 |
Senior Member
Join Date: Jun 2009
Posts: 1,875
Rep Power: 33 |
Not at all. I used to teach, and I would have likely asked you similar questions as you were making progress in your project.
The idea here is to help you find the solution by yourself and help solidify your modeling skills. You are making great progress.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
May 8, 2024, 18:28 |
Wall shear
|
#45 |
Member
Kevin Gnanaraj
Join Date: Apr 2024
Posts: 34
Rep Power: 2 |
Hi Glen, I've plotted the wall shear. It looks like the flow separates at a similar location for all rotational speeds, near the leading edge of the propeller. For almost all RPM values, the load distribution seems to be the same, with the highest pressures concentrated at the leading edge, and dropping off very quickly after that. Is there a way to tell whether the propeller is meant to be laminar or turbulent?
I can see that the contour plots show irregularities in the blade profile, so it seems most likely that my CAD model is the culprit here. Moving on to contra rotating coaxial rotors, would it still be acceptable to use frozen rotor as the interface between the rotors? From my understanding, FR is highly dependent on the rotor position, but I want a more general idea of the performance. Would I use a periodic interface, or stage, then? edit: I read up a bit more about pitch, so seems that mixing planes is not suitable at all. Here are the shear plots: 2273 RPM 2273_wall_shear.png 3727 RPM 3727_wall_shear.png 4471 RPM 4471_wall_shear.png Last edited by keg504; May 8, 2024 at 18:46. Reason: Correction after reading |
|
May 8, 2024, 19:45 |
|
#46 | |
Senior Member
Join Date: Jun 2009
Posts: 1,875
Rep Power: 33 |
Quote:
The interface between counter-rotating rotors should be either a "transient rotor-stator sliding interface" for a transient simulation, or a "mixing plane" for a steady-state simulation.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
||
May 8, 2024, 20:21 |
|
#47 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,862
Rep Power: 144 |
No - you are looking at the wrong wall shear.
If the flow is attached the wall shear will be in the direction of flow of the air over the surface. If the flow is separated it will be against the direction of the bulk air flow. So the wall shear direction you should be looking at is y direction on your images in post #45. You showed the Z direction which is not very meaningful for observing separations. Even better, draw vectors on the blade surface of wall shear using the X, Y and Z components. That should really clearly show any separation as a reversal of the wall shear. This approach has the advantage of working no matter the surface curvature is or the orientation of the blade. Looking at the Wall Shear Y component will only work in the current position of the blades.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
May 9, 2024, 07:01 |
Shear vectors
|
#48 |
Member
Kevin Gnanaraj
Join Date: Apr 2024
Posts: 34
Rep Power: 2 |
I did plot x and y shear contours, they both looked mostly negative, so I didn't think they were useful. But here are the shear vectors. I'm not quite sure how to interpret them, but it looks like the flow remains attached for the most part.
2273 RPM 2273_wall_shear_vect.png 3737 RPM 3727_wall_shear_vect.png 4471 RPM 4471_wall_shear_vect.png |
|
May 9, 2024, 07:08 |
|
#49 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,862
Rep Power: 144 |
The recirculations are usually a much lower velocity than the main flow. So you will have to scale the vector length up so you can see what is happening in the low flow regions.
The reason I ask this is because attached flows can be predicted quite accurately with normal turbulence models, but separated flows often require a more sophistocated model to get high accuracy. So working out whether you have significant separated regions determines which approach you need to look at.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
May 9, 2024, 07:54 |
Larger vectors
|
#50 |
Member
Kevin Gnanaraj
Join Date: Apr 2024
Posts: 34
Rep Power: 2 |
I scaled up the vectors to be as large as they would go in CFD-Post. Here are the results. It looks like the flows do remain attached, except near the hub, as you stated.
2273 RPM 2273_wall_shear_vect.png 3727 RPM 3727_wall_shear_vect.png 4471 RPM 4471_wall_shear_vect.png Last edited by keg504; May 9, 2024 at 07:54. Reason: Added information |
|
May 9, 2024, 08:02 |
|
#51 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,862
Rep Power: 144 |
Yes, that does look like there are no separations. That makes things considerably easier - it means a normal 2 equation turbulence model should be quite accurate.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
May 14, 2024, 10:54 |
Coaxial simulations
|
#52 |
Member
Kevin Gnanaraj
Join Date: Apr 2024
Posts: 34
Rep Power: 2 |
So I've been running simulations of the co-axial configuration. I've noticed a discontinuity in pressure across the mixing plane interface (I'm using stage average velocity with implicit averaging). Is that something to be worried about? The total pressure and velocity seem to remain continuous.
Also, I'm pretty sure the geometry is the problem in my case. I've been running a 60 in propeller that I have the manufacturer CAD for, and even without a proper mesh study (I'm using a similar number of elements), the difference is only ~3% under. Here you can see the total pressure and pressure plots for the coaxial simulation (the blank space in the middle is meant to represent the motor, it is a cylinder): Total pressure 3429_coax_final_tot_press.png Pressure 3429_coax_final_press.png |
|
May 14, 2024, 11:34 |
|
#53 | |
Senior Member
Join Date: Jun 2009
Posts: 1,875
Rep Power: 33 |
Quote:
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
||
May 14, 2024, 11:54 |
Stage average velocity
|
#54 |
Member
Kevin Gnanaraj
Join Date: Apr 2024
Posts: 34
Rep Power: 2 |
My reasoning is that even though I have turned on compressible flow, the flow isn't actually at Mach 0.3, so constant total pressure would compromise on the simulation if one of the two has to be reduced in relation to the other to maintain it. I do have it in mind to test with the other option, but time is not on my side right now.
I'm aware that the reasoning is incomplete and probably misunderstanding something, but I couldn't find a reference on which option to choose. edit: also, I was basing it on blade element theory, namely I'm attempting to approximate propeller velocity by using the stage average velocity Last edited by keg504; May 14, 2024 at 16:41. Reason: Added information on choice of mixing plane option |
|
May 14, 2024, 18:41 |
|
#55 | |
Senior Member
Join Date: Jun 2009
Posts: 1,875
Rep Power: 33 |
Quote:
The old default, pre-Ansys CFX-18R1 used to be stage average velocity and it was superseded by a more accurate model. Total pressure behavior should be irrelevant to compressible or incompressible fluid. If it were dependent on that difference, it would be a big problem since a user would need to know when to blend the models, i.e. which parameter determines which model to use.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
||
May 14, 2024, 18:52 |
|
#56 |
Member
Kevin Gnanaraj
Join Date: Apr 2024
Posts: 34
Rep Power: 2 |
Thanks for the clearing that up. I hope I have time to change it, but probably not right now, so I will use as-is for now.
Would you be so kind as to explain the actual difference between stage average velocity and constant total pressure, so that I know for future reference, please? |
|
May 23, 2024, 07:56 |
Summary
|
#57 |
Member
Kevin Gnanaraj
Join Date: Apr 2024
Posts: 34
Rep Power: 2 |
I changed it to constant total pressure, as you suggest, Opaque, and it helped a little, but not drastically, though I don't think it would, just would stand up more to scrutiny.
I think my problem has been largely solved, thank you ghorrocks and opaque for your help and patience. There might be optimisations that can be done, I'm not too worried though. As a summary for those might setup a similar case here it is: For a steady state simulation of static propeller thrust in an outdoor environment, use a stationary outer domain with air as an ideal gas, with an inlet set to an opening condition, and outlets set to also be opening. All openings have a reference pressure of 0. The inlet and outlet should be separate, otherwise the imbalances will not converge. The size of the air domain used is 8 * rotor diameter. The reference pressure is 1 atm. Compressible flow is also turned on. The rotor domain is rotating, using a frozen rotor interface with the air, and a mixing plane interface (using constant total pressure) with the other rotor. The rotor diameter used is a cylinder with diameter 1.1 * rotor diameter and a length of 0.4 * rotor diameter (the length can vary based on rotor spacing in the coaxial configuration). The domain dimensions come from "3D CFD Simulation and Experimental Validation of Small APC Slow Flyer Propeller Blade" (Kutty and Rajendran, 2017). If I've missed anything, let me know, and I'll edit this post . |
|
June 14, 2024, 10:12 |
Coaxial experimental data
|
#58 |
Member
Kevin Gnanaraj
Join Date: Apr 2024
Posts: 34
Rep Power: 2 |
The test rig for the coaxial configuration has been mostly set-up, save for automated data collection, but preliminary testing using multimeters and hand held stroboscopes are showing that my simulations roughly match the experimental data after the aforementioned nebulous adjustment factor. Better data will be collected, but right now, I would say the error is about 10%.
|
|
Tags |
floating point exception, mixing planes, propeller flow error |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[OpenFOAM] Paraview python script, creating data using only CLI, saving in csv/excel file | Ash Kot | ParaView | 1 | September 24, 2021 13:23 |
Error running AMI propeller simulation | luitzor | OpenFOAM Running, Solving & CFD | 0 | April 19, 2021 14:48 |
Pump CAD + experimental data for CFD verification study | bemism | Main CFD Forum | 0 | July 20, 2017 16:30 |
Data Produced From Fine Marine Cant Match with The Experimental Data | PeiSan | Fidelity CFD | 4 | August 23, 2014 06:33 |
How to compare the average velocity of the simulation with the Experimental data ? | nanavati | OpenFOAM Post-Processing | 2 | August 22, 2014 05:48 |