louisgag
...
updated 6 years ago
The FreeCAD files: forum.freecad.org/viewtopic.php?t=85139
The main goal and difference between this aircraft and other cyclogyros is that it has a single drive shaft and still allows high agility.
The dream behind this device would be to make a human powered VTOL implementation of it, as you can see in the related patent (number DE102022114599). This German patent is also available for licensing opportunities or acquisition.
The rotor has 4 blades, a radius of 84.25 mm and is designed to rotate at up to 5300 rpm.
A prototype rotor with the same design has been built and tested and is presented in the paper "Aerodynamic Performance of a Hovering Cycloidal Rotor: A CFD Study with Experimental Validation," which will be soon published by the Journal of the American Helicopter Society.
This aircraft was designed and tested for the project "A Novel and Simple Aircraft Requiring Minimal Power to Hover" funded by the Alexander von Humboldt Foundation Research Fellowships for Postdoctoral Researchers and conducted at the Institute of Aerodynamics and Gas Dynamics at the University of Stuttgart. This project would have never been possible without the support of many colleagues at the institute. More info: http://louisgagnon.com/research/AvH_cyclo.html
The 3D model and animation were done in FreeCAD with its Assembly4 and Render workbenches, using many tricks and Python to make things run smoothly.
To facilitate code improvement by developers and the community, a repository (doi.org/10.18419/darus-3340) provides all the necessary files to either reconstruct the case or run it as is.
The rotor has 4 blades, a radius of 84.25 mm and rotates at 1695 RPM. The maximum pitch angles of the upper and lower paths are 35° and 38°, respectively, and they occur at positions 100° and 277°.
The case was run using OpenFOAM (pimpleFoam) on the HLRS-Hawk supercomputer.
A more detailed description of the case is available in the (currently preprint) paper "Parametric Aerodynamic Study of a 3D Printed Cycloidal Rotor Prototype Designed for the Best Figure of Merit" which was submitted to the Journal of the American Helicopter Society: researchgate.net/publication/369473318_Parametric_Aerodynamic_Study_of_a_3D_Printed_Cycloidal_Rotor_Prototype_Designed_for_the_Best_Figure_of_Merit
Project Info: “A Novel and Simple Aircraft Requiring Minimal Power to Hover” project funded by the Alexander von Humboldt Foundation Research Fellowships for Postdoctoral Researchers: http://louisgagnon.com/research/AvH_cyclo.html and conducted at the Institute of Aerodynamics and Gas Dynamics
To facilitate code improvement by developers and the community, a repository (doi.org/10.18419/darus-3340) provides all the necessary files to either reconstruct the case or run it as is.
The rotor has 4 blades, a radius of 84.25 mm and rotates at 1695 RPM. The maximum pitch angles of the upper and lower paths are 35° and 38°, respectively, and they occur at positions 100° and 277°.
The case was run using OpenFOAM (pimpleFoam) on the HLRS-Hawk supercomputer.
A more detailed description of the case is available in the (currently preprint) paper "Parametric Aerodynamic Study of a 3D Printed Cycloidal Rotor Prototype Designed for the Best Figure of Merit" which was submitted to the Journal of the American Helicopter Society: researchgate.net/publication/369473318_Parametric_Aerodynamic_Study_of_a_3D_Printed_Cycloidal_Rotor_Prototype_Designed_for_the_Best_Figure_of_Merit
Project Info: “A Novel and Simple Aircraft Requiring Minimal Power to Hover” project funded by the Alexander von Humboldt Foundation Research Fellowships for Postdoctoral Researchers: http://louisgagnon.com/research/AvH_cyclo.html and conducted at the Institute of Aerodynamics and Gas Dynamics
Tip: watch full screen with 1440p resolution (gear at bottom, quality menu) and loop (right-click menu).
This video shows the streamlines (velocity lines rendered as LIC) and pressure distribution in the middle of the cycloidal rotor blade for both a URANS (above) and a DDES (below) case using the same mesh and timestep.
Both cases are ran in 3D using OpenFOAM (pimpleFoam) on the HLRS-Hawk computers.
The rotor has 4 blades, a radius of 84.25 mm and turns at 5300 RPM. The maximum pitch angles of the upper and lower paths are 35° and 38°, respectively, and they occur at position 100° and 277°, respectively.
Typical dynamic stall behavior can be observed: blade pressure and suction sides switch between the upper and lower portions of the rotation and dynamic stall occurs and the end of the upper course.
NOTE: This is an improved version of my "Pressure plot and LIC-Streamlines in the middle of a Cyclorotor blade" video.
Project Info: “A Novel and Simple Aircraft Requiring Minimal Power to Hover” project funded through the Alexander von Humboldt Foundation Research Fellowships for Postdoctoral Researchers: http://louisgagnon.com/research/AvH_cyclo.html
Rotor has a radius of 84.25 mm and 4 blades.
The cycloidal rotor is represented by a cylinder having the radius and span length as the rotor.
Project Info: “A Novel and Simple Aircraft Requiring Minimal Power to Hover” project funded through the Alexander von Humboldt Foundation Research Fellowships for Postdoctoral Researchers: http://louisgagnon.com/research/AvH_cyclo.html
This video shows the streamlines and pressure distribution in the middle of the cyclorotor blade for both a URANS (above) and a DDES (below) case using the same mesh.
Both cases are ran in 3D using OpenFOAM (pimpleFoam) on the HLRS-Hawk computers.
The rotor has 4 blades, a radius of 84.25 mm, and turns at 5300 RPM.
This animation show the interpolation interfaces of the background mesh (light yellow) and component blade mesh (here dark pink).
More info on GitHub: github.com/louisgag/openFoam-Overset-SpeedUp
The mesh shown is for the overset twin version of this AMI-case: youtube.com/watch?v=q0pafX63_x0
Project Info: “A Novel and Simple Aircraft Requiring Minimal Power to Hover” project funded through the Alexander von Humboldt Foundation Research Fellowships for Postdoctoral Researchers: http://louisgagnon.com/research/AvH_cyclo.html
The rendering shows the lambda2 isosurfaces colored by distance.
The simulation is described in an upcoming paper named "Investigation of the Reynolds Number on the Performance of a Cycloidal Rotor" as: "a NACA 0012 airfoil with a chord length of c = 0.15 m [with] an initial angle of attack m (AoA) of 10° and [...] oscillating with an angular velocity of Ω = 46.6667 rad/s. The free stream velocity is U∞ = 35 m/s with a turbulence intensity of Tu = 0.08%. The kinematic viscosity is ν = 38.889 · 10−6 m^2/s ."
Accompanying paper: doi.org/10.1115/1.4056844 and bachelor thesis: https://stg.ibs-bw.de/aDISWeb/app?service=direct/0/Home/$DirectLink&sp=SOPAC02&sp=SAKSWB-IdNr1763925978
Project Info: “A Novel and Simple Aircraft Requiring Minimal Power to Hover” project funded through the Alexander von Humboldt Foundation Research Fellowships for Postdoctoral Researchers: http://louisgagnon.com/research/AvH_cyclo.html
This videos accompanied the presentation of the master thesis research done by Ms. Doudou Huang at IAG on the effect of dynamically morphing blades on the efficiency of a cycloidal rotor. Presentation: http://louisgagnon.com/articlesScientifiques/HuangGagnon2021_EASN_morphingCyclorotor.pdf
Project Info: “A Novel and Simple Aircraft Requiring Minimal Power to Hover” project funded through the Alexander von Humboldt Foundation Research Fellowships for Postdoctoral Researchers: http://louisgagnon.com/research/AvH_cyclo.html
The CFD model is a small progression based on the original work by Shawn Cogan: https://stg.ibs-bw.de/aDISWeb/app?service=direct/0/Home/$DirectLink&sp=SOPAC02&sp=SAKSWB-IdNr1742829899 which was presented at the 10th EASN Conference (https://easnconference.eu/sites/default/files/10th_EASN_Virtual_Conference-Prefinal_Agenda.pdf) and the video will be shown during the 11th EASN Conference: https://easnconference.eu/sites/default/files/11th_EASN_Virtual_Conference-Preliminary_Agenda_4.pdf
Project Info: “A Novel and Simple Aircraft Requiring Minimal Power to Hover” project funded through the Alexander von Humboldt Foundation Research Fellowships for Postdoctoral Researchers: http://louisgagnon.com/research/AvH_cyclo.html
Original blog entry: http://louisgagnon.com/scBlog/rotor3DCFD.html
What is shown:
00:08 - Relative atmospheric pressure on the blades, the effects of the major dynamic stall event are mostly seen shortly after 270°. Most of the thrust is generated by pressure and one can see that the bottom passage has the strongest contribution.
00:45 - The coefficient of friction on the blade, considering that the blade rotates and oscillates, identifies in yellow positions where the relative velocity between blade and fluid goes to zero and indicates likelihood for detached flow. Once again, the position shortly after 270° shows the strongest action, because of dynamic stall.
01:08 - The velocity LIC lines in the rotating rotor reference frame (rotating but not locally oscillating) on the first layer of cells above the blade give another, more global, view of the dynamic stall action and the vortices that are created by the blade.
01:33 - The lambda2 contour surfaces show these vortices, such has the horseshoe vortex which detaches from the blade around angle 360° and the tip vortices of the blades. The lambda2 surfaces in the first 1 cm above the blade are cut out from the rendering to allow better visualization of the blade's coefficient of friction.
02:04 - The same lambda2 surfaces are now colored by the vorticity projected on the non-oscillating orbital blade path to highlight the direction of the vortices, in particular the tip vortices which have a different direction of rotation on the top passage as in the bottom one, because of the combination direction of motion of the blade and the thrust it generates.
02:36 - The sensitivity of the lambda2 surfaces is extended to identify the large turbulent structures that leave the upper blades and travel through the rotor. The pressure bubbles shown on the rearwards half the blades correlate with the magnitude of the generated thrust and also indicate that the large eddies coming from the upper blade have little influence on the lower passage. The LIC lines indicate the path of the flow inside the vortices identified by the lambda2 surfaces and the color indicates the radial distance from the center of the rotor, to ease understanding of motion of the turbulent structures traveling through the rotor.
03:16 - The arrows now allow identifying the inflow and wake generated by the rotor. Their size is scaled according to the square root of the velocity vector magnitude and they are shown for any velocity magnitude above 5 cm/s, thus allowing to identifying the zone of influence of the rotor. Although the arrows are filtered out near the blades, some local irregularities are seen and correspond to the flow response to the passage of a blade. Because of the small aspect ratio and the absence of endplates or side disks, the inflow is highly three-dimensional and the flow leaving the rotor forms, contrary to what is usually seen in 2D simulations, a rather wide wake. The arrows also highlight the acceleration of the flow through the rotor and also the marked difference of velocity between the upper and lower passages of the blades.
04:10 - The side view of the wake further pictures that air flows into the rotor from all sides and also in the middle of the rotor, thus between the two blade passages, while the wake remains rather bidimensional. As expected the flow is also curved through the rotor, thus bending the wake towards the positive x-direction.
Geometry: Blade chord is 15 cm and aspect ratio is 2; rotor radius is 50 cm. The motion function is given at the beginning of the video.
CFD model: Turbulence is implemented using the unsteady Reynolds-averaged Navier-Stokes method with a k-ω SST turbulence model with a fully resolved wall using a first grid cell thickness of 6.6 μm. The v2012 pimpleFoam solver is used and the grid has pill-shaped moving mesh zones which contain the blades and which are coupled to a rotating background mesh through the AMI interface method. No symmetry is used and the complete mesh has 20.5 million cells. Second order schemes in both time and space are used.
Computational aspects: The simulation was run on the Hawk supercomputer of the High Performance Computing Center Stuttgart (HLRS) using 2048 processes on 32 nodes, taking approximately 6 seconds per timestep, each covering 0.5 degrees of rotation, thus leading to a maximum Courant number varying between 160 and 200.
Acknowledgments: This simulation was conducted at the Institute of Aerodynamics and Gas Dynamics of the University of Stuttgart for the "A Novel and Simple Aircraft Requiring Minimal Power to Hover" and is financed by the Humboldt Research Fellowship for postdoctoral researchers. Free open-source software used: FreeCAD and naca4gen and GNU Octave for the blade profile, OpenFOAM for the simulation, and ParaView and Python for post-processing.
Made with FreeCAD's sketcher panel
Project Info: “A Novel and Simple Aircraft Requiring Minimal Power to Hover” project funded through the Alexander von Humboldt Foundation Research Fellowships for Postdoctoral Researchers: http://louisgagnon.com/research/AvH_cyclo.html
Project Info: “A Novel and Simple Aircraft Requiring Minimal Power to Hover” project funded through the Alexander von Humboldt Foundation Research Fellowships for Postdoctoral Researchers: http://louisgagnon.com/research/AvH_cyclo.html
Project Info: “A Novel and Simple Aircraft Requiring Minimal Power to Hover” project funded through the Alexander von Humboldt Foundation Research Fellowships for Postdoctoral Researchers: http://louisgagnon.com/research/AvH_cyclo.html
This is a companion video for the following article: Gagnon, L., Masarati, P., « Autonomous Untethered Flight of Multibody Dynamics Rotorcraft With Cycloidal Rotors », DETC2018-85152, ASME IDETC/CIE 2018, Québec City, Canada, August 26-29, 2018
Perfect static position match is still not attained, possibly due to the combination of objectives between energy use, cabin stability, and path following capacity. Other possibilities: integrator part of the PID, or SAS parameters. Model optimized using a different path than the one displayed here. Each bar represents one meter.
More info on the research project: http://louisgagnon.com and http://www.researchgate.net/publication/310625903_Flight_scenarios_study_of_two_cycloidal_rotor_aircraft_concepts
More info on the research project: http://louisgagnon.com and http://www.researchgate.net/publication/310625903_Flight_scenarios_study_of_two_cycloidal_rotor_aircraft_concepts
checkout http://www.electricsnowmobile.mcgill.ca for more!
checkout http://www.electricsnowmobile.mcgill.ca for more!
checkout http://www.electricsnowmobile.mcgill.ca for more!


