|
PROJECT REPORT: VA170 |
|
Optimisation of a Motorbike Front Windshield and Fairing Using Sculptor |
Haakon Dahle Smith
22
October, 2004
This report presents the results of an optimisation carried out on the front windshield and fairing geometry of a Yamaha R1 motorbike to reduce drag. Sculptor was used to directly deform the CFD case file therefore removing the need to re-mesh.
The fairing geometry was controlled using four parameters, which were optimised using a response-surface method. Twenty-five runs were initially used to define the preliminary response surface and thereafter a total of four cases were run attempting to optimise the fairing for low drag.
The motorbike and rider were modelled in an upright position with rotating wheels and a rolling road. The model was run full scale at a speed of 27.8 m/s (100 kph)
The optimised front fairing / windshield geometry reduced drag by 3.75 N, 2.5 % of the total drag.
CONTENTS
This report follows on from previous studies carried out by ACFD on a Yamaha R1 motorbike, VA109 and VA27. The previous reports noted that although there is no clear requirement for downforce in motorbike racing the reduction of drag is vital in improving the performance of the bike.
A 2.3 million cell case was taken from VA109 and Sculptor was used to deform the fairing / windshield geometry in the case file directly. A response surface methodology was then used to optimise the deformation parameters.
Figure 1.1 identifies the front fairing and windshield (coloured in blue) that were modified as part of the study. Figure 1.2 shows the surface mesh on the bike.
Figure 1.1: CAD geometry of motorbike with fairing / windshield in blue
Figure 1.2: Surface mesh of bike
During the optimisation there were four parameters used that would alter different aspects of the fairing / windshield. The parameters were:
An overview of the deformations that were carried out is shown visually in Figures 1.3. The “net” surrounding the fairing and windshield, used to make the geometry deformations, can also be seen. An MPEG movie of the deformations is available on the web copy of this report. Once these parameters had been defined the deformations could be made interactively – directly to the case file without the need to re-mesh. Note that Figure 1.3 shows the limits used in the study.
|
Deformation |
Maximum Deflection |
Minimum Deflection |
|
Height of rear glass section |
||
|
Height of front / middle glass section |
||
|
Width of top section of fairing |
|
|
|
Height of top section of fairing |
|
|
Figure 1.3: The maximum and minimum deformations of the 4 parameters
The response-surface method required 25 geometries to be analysed to define the initial response-surface for the 4 parameters. A further 4 cases were subsequently recommended by the response-surface method in order to search for an optimum solution based on the performance criterion of lowest drag.
The man time used in setting up of all 30 runs in Sculptor was very small in proportion to the computational time taken to solve the case. The initial ASD volume allowing deformations to take place took 2-3 hours to setup, thereafter each change in the fairing geometry only took 2-3 minutes plus 10 minutes to save each case. A summary of all the drag and lift forces that were generated can be seen in Table 2.1 and are plotted in Figures 2.1 and 2.2.
Each case made use of the baseline data file as an initial flow solution and was therefore only required to be run-on 1000 more iterations over the previous 2500 of the baseline.
Table 2.1: Summary of forces over all runs
Figure 2.1: Variation of Lift Over All Runs
Figure 2.2: Variation of Drag Over All Runs
The run chosen for the best drag characteristics was RUN-it-04. It lost 3.76 N of drag (2.5% of total) whilst gaining 6 N of lift (28%). Figure 3.1 shows the difference in fairing / windshield geometry between the baseline (blue) and RUN-it-04 (red).
Figure 3.1: Difference in fairing / windshield geometry between baseline and RUN-it-04
Table 3.1 shows a summarised version of the full force table (Table 2.1) highlighting the differences in forces between the baseline and RUN-it-04.
Table 3.1: Summary of forces for baseline and RUN-it-04
Table 3.2 shows the difference in cooling between the baseline and RUN-it-04. This is represented by the drop in pressure, DP, across the radiator. It can be seen that the change is negligible.
Table 3.2: Summary of cooling flow for baseline and RUN-it-04
Figure 3.1 shows that the optimised geometry has a slightly lower rear section to the glass whilst the front and mid section of the glass has been raised more significantly. The width of the top section of the fairing has been narrowed and the height of the top section has also been lowered.
Figure 3.2 shows a comparison of contours of static pressure coefficient acting on the motorcycle and rider for the baseline and optimised cases. It can be seen on the fairing / windshield of RUN-it-04 that there is a larger area of lower pressure on the windshield (coloured in dark purple) and therefore a lower drag force. The corresponding area on the baseline case (coloured in blue) is at a higher pressure, meaning that there is a greater drag force on the component.
Figure 3.2: Banded contours of static pressure coefficient for baseline and RUN-it-04
The following conclusions can be drawn from the study:
Further work to be investigated could take the form of:
|
GENERAL DETAILS |
|
|
Case & Data file(s) |
Bike-upright-wheels-rotating.cas |
|
Mesh Type (if hybrid - areas of hex) |
Hybrid Two wake blocks |
|
Number of cells |
2.3e06 |
|
SOLVER CONTROLS |
|
|
Turbulence Model |
Spalart-Allmaras |
|
Near-wall Treatment |
Log-law-of the wall |
|
Discretization (Under-relaxation) Pressure Momentum Pressure-Velocity Coupling Turbulence |
2nd Order (0.2) QUICK (0.6) SIMPLE QUICK (0.7) |
|
Buoyancy forces |
Off |
|
MATERIALS |
|
|
Density Viscosity |
1.225 kg/m3 1.7894e-05 kg/ms |
|
BOUNDARY CONDITIONS |
|
|
Main Inlet |
Case A : 27.8m/s |
|
|
Turbulence intensity 0.8%; length scale 0.2m Dirn vector (1 ,0, 0) |
|
Ground |
As main inlet |
|
Front Wheel Case A : |
Axis X=0.00249; Z=-0.00462 w = -91.27 rad/s Contact Patch = (0.01278, -0.07687, 0) |
|
Rear Wheel Case A : |
Axis X=0.00603, Z=0.01052 w = -89.61 rad/s |
|
Radiator (porous media) Dirn vector 1 Dirn vector 2 Viscous resistance Inertial resistance |
(0.9377, -0.0045, 0.3474) (0.0007, 0.9999, 0.0112) (0, 0, 0) Not known |