Aerodynamic study of time-trial helmets in cycling racing using CFD analysis

The aerodynamic drag of three different time-trial cycling helmets was analyzed numerically for two different cyclist head positions. Computational Fluid Dynamics (CFD) methods were used to investigate the detailed airflow patterns around the cyclist for a constant velocity of 15 m/s without wind. The CFD simulations have focused on the aerodynamic drag effects in terms of wall shear stress maps and pressure coefficient distributions on the cyclist/helmet system. For a given head position, the helmet shape, by itself, obtained a weak effect on a cyclist’s aerodynamic performance (<1.5%). However, by varying head position, a cyclist significantly influences aerodynamic performance; the maximum difference between both positions being about 6.4%. CFD results have also shown that both helmet shape and head position significantly influence drag forces, pressure and wall shear stress distributions on the whole cyclist’s body due to the change in the near-wake behavior and in location of corresponding separation and attachment areas around the cyclist.     

Aerodynamic investigation of the thermo-dependent flow structure in the wake of a cyclist

The main purpose of this study was to assess the influence of the environmental temperature on both the aerodynamic flow evolving around the bicycle and cycling power output. The CFD method was used to investigate the detailed flow field around the cyclist/bicycle system for a constant speed of 11.1 m/s (40 km/h) without wind. In complement, a mathematical model was used to determine the temperature-dependent power output in the range [−10; 40 °C]. The numerical investigation gives valuable information about the turbulent flow field in the cyclist's wake which evolves accordingly the surrounding temperature. A major result of this study is that the areas of overpressure upstream of the cyclist and of underpressure downstream of him are less extensive for a temperature of 40 °C compared to −10 °C. The results suggest that the aerodynamic braking effect of the bicycle is minimized when the air temperature is high, as a lower air density results in a reduction in drag on the cyclist. This study showed that the power required to maintain a constant speed is reduced when the temperature is high, the reason being a lower aerodynamic resistance.     

Aerodynamic study of different cyclist positions: CFD analysis and full-scale wind-tunnel tests

Three different cyclist positions were evaluated with Computational Fluid Dynamics (CFD) and wind-tunnel experiments were used to provide reliable data to evaluate the accuracy of the CFD simulations. Specific features of this study are: (1) both steady Reynolds-averaged Navier–Stokes (RANS) and unsteady flow modelling, with more advanced turbulence modelling techniques (Large-Eddy Simulation – LES), were evaluated; (2) the boundary layer on the cyclist’s surface was resolved entirely with low-Reynolds number modelling, instead of modelling it with wall functions; (3) apart from drag measurements, also surface pressure measurements on the cyclist’s body were performed in the wind-tunnel experiment, which provided the basis for a more detailed evaluation of the predicted flow field by CFD. The results show that the simulated and measured drag areas differed about 11% (RANS) and 7% (LES), which is considered to be a close agreement in CFD studies. A fair agreement with wind-tunnel data was obtained for the predicted surface pressures, especially with LES. Despite the higher accuracy of LES, its much higher computational cost could make RANS more attractive for practical use in some situations. CFD is found to be a valuable tool to evaluate the drag of different cyclist positions and to investigate the influence of small adjustments in the cyclist’s position. A strong advantage of CFD is that detailed flow field information is obtained, which cannot easily be obtained from wind-tunnel tests. This detailed information allows more insight in the causes of the drag force and provides better guidance for position improvements.     

Optimal cycling time trial position models: Aerodynamics versus power output and metabolic energy

The aerodynamic drag of a cyclist in time trial (TT) position is strongly influenced by the torso angle. While decreasing the torso angle reduces the drag, it limits the physiological functioning of the cyclist. Therefore the aims of this study were to predict the optimal TT cycling position as function of the cycling speed and to determine at which speed the aerodynamic power losses start to dominate. Two models were developed to determine the optimal torso angle: a ‘Metabolic Energy Model’ and a ‘Power Output Model’. The Metabolic Energy Model minimised the required cycling energy expenditure, while the Power Output Model maximised the cyclists? power output. The input parameters were experimentally collected from 19 TT cyclists at different torso angle positions (0–24°). The results showed that for both models, the optimal torso angle depends strongly on the cycling speed, with decreasing torso angles at increasing speeds. The aerodynamic losses outweigh the power losses at cycling speeds above 46 km/h. However, a fully horizontal torso is not optimal. For speeds below 30 km/h, it is beneficial to ride in a more upright TT position. The two model outputs were not completely similar, due to the different model approaches. The Metabolic Energy Model could be applied for endurance events, while the Power Output Model is more suitable in sprinting or in variable conditions (wind, undulating course, etc.). It is suggested that despite some limitations, the models give valuable information about improving the cycling performance by optimising the TT cycling position.     

Difference between palpation and optoelectronics recording of scapular motion

The aim of this study is to determine the errors of scapular localisation due to skin relative to bone motion with an optoelectronic tracking system. We compared three-dimensional (3D) scapular positions obtained with skin markers to those obtained through palpation of three scapular anatomical landmarks. The scapular kinematics of nine subjects were collected. Static positions of the scapula were recorded with the right arm elevated at 0°, 40°, 80°, 120° and 160° in the sagittal plane. Palpation and subsequent digitisation of anatomical landmarks on scapula and thorax were done at the same positions. Scapular 3D orientation was also computed during 10 repeated movements of arm elevation between 0° and 180°. Significant differences in scapular kinematics were seen between static positions and palpation when considering anterior/posterior tilt and upward/downward rotation at angles over 120° of humeral elevation and only at 120° for internal/external rotation. There was no significant difference between positions computed during static positions and during the movement for the three scapular orientations. A rotation correction model is presented in order to reduce the errors between static position and palpation measurement.     

Aerodynamics of cyclist posture, bicycle and helmet characteristics in time trial stage

The present work is focused on the aerodynamic study of different parameters, including both the posture of a cyclist's upper limbs and the saddle position, in time trial (TT) stages. The aerodynamic influence of a TT helmet large visor is also quantified as a function of the helmet inclination. Experiments conducted in a wind tunnel on nine professional cyclists provided drag force and frontal area measurements to determine the drag force coefficient. Data statistical analysis clearly shows that the hands positioning on shifters and the elbows joined together are significantly reducing the cyclist drag force. Concerning the saddle position, the drag force is shown to be significantly increased (about 3%) when the saddle is raised. The usual helmet inclination appears to be the inclination value minimizing the drag force. Moreover, the addition of a large visor on the helmet is shown to provide a drag coefficient reduction as a function of the helmet inclination. Present results indicate that variations in the TT cyclist posture, the saddle position and the helmet visor can produce a significant gain in time (up to 2.2%) during stages.     

Computational fluid dynamics analysis of drag and convective heat transfer of individual body segments for different cyclist positions

This study aims at investigating drag and convective heat transfer for cyclists at a high spatial resolution. Such an increased spatial resolution, when combined with flow-field data, can increase insight in drag reduction mechanisms and in the thermo-physiological response of cyclists related to heat stress and hygrothermal performance of clothing. Computational fluid dynamics (steady Reynolds-averaged Navier-Stokes) is used to evaluate the drag and convective heat transfer of 19 body segments of a cyclist for three different cyclist positions. The influence of wind speed on the drag is analysed, indicating a pronounced Reynolds number dependency on the drag, where more streamlined positions show a dependency up to higher Reynolds numbers. The drag and convective heat transfer coefficient (CHTC) of the body segments and the entire cyclist are compared for all positions at racing speeds, showing high drag values for the head, legs and arms and high CHTCs for the legs, arms, hands and feet. The drag areas of individual body segments differ markedly for different cyclist positions whereas the convective heat losses of the body segments are found to be less sensitive to the position. CHTC-wind speed correlations are derived, in which the power-law exponent does not differ significantly for the individual body segments for all positions, where an average value of 0.84 is found. Similar CFD studies can be performed to assess drag and CHTCs at a higher spatial resolution for applications in other sport disciplines, bicycle equipment design or to assess convective moisture transfer.     

Individual flight styles in ski jumping: results obtained during Olympic Games competitions

From the physics point of view, the jump length in ski jumping depends on: the in-run velocity v(0), the velocity perpendicular to the ramp v(p0) due to the athlete's jumping force, the lift and drag forces acting during take-off and during the flight, and the weight of the athlete and his equipment. The aerodynamic forces are a function of the flight position and of the equipment features. They are a predominant performance factor and can largely be influenced by the athlete. The field study conducted during the Olympic Games competitions 2002 at Park City (elevation: 2000 m) showed an impressive ability of the Olympic medallists to reproduce their flight style and remarkable differences between different athletes have been found. The aerodynamic forces are proportional to the air density. Elite athletes are able to adapt their flight style to thin air conditions in order to maximise jump length and to keep the flight stable. The effects of flight position variations on the performance have been analysed by means of a computer model which is based on the equations of motion and on wind tunnel data corresponding to the flight positions found in the field. Athletes have to solve extremely difficult optimisation problems within fractions of a second. The computer simulation can be used as a reliable starting point for the improvement of training methods and gives an insight into the "implicit" knowledge of physics that the ski jumping athlete must have available for a good performance.     

Lumbar spine angles and intervertebral disc characteristics with end-range positions in three planes of motion in healthy people using upright MRI

Understanding changes in lumbar spine (LS) angles and intervertebral disc (IVD) behavior in end-range positions in healthy subjects can provide a basis for developing more specific LS models and comparing people with spine pathology. The purposes of this study are to quantify 3D LS angles and changes in IVD characteristics with end-range positions in 3 planes of motion using upright MRI in healthy people, and to determine which intervertebral segments contribute most in each plane of movement. Thirteen people (average age = 24.4 years, range 18–51 years; 9 females; BMI = 22.4 ± 1.8 kg/m²) with no history of low back pain were scanned in an upright MRI in standing, sitting flexion, sitting axial rotation (left, right), prone on elbows, prone extension, and standing lateral bending (left, right). Global and local intervertebral LS angles were measured. Anterior-posterior length of the IVD and location of the nucleus pulposus was measured. For the sagittal plane, lower LS segments contribute most to change in position, and the location of the nucleus pulposus migrated from a more posterior position in sitting flexion to a more anterior position in end-range extension. For lateral bending, the upper LS contributes most to end-range positions. Small degrees of intervertebral rotation (1–2°) across all levels were observed for axial plane positions. There were no systematic changes in IVD characteristics for axial or coronal plane positions.     

A computer simulation of free-range exercise in the laboratory.

We present a technique for simulating dynamic field (free-range) exercise, using a novel computer-controlled cycle ergometer. This modified cycle ergometer takes into account the effect of friction and aerodynamic drag forces on a 70-kg cyclist in a racing position. It also affords the ability to select different gear ratios. We have used this technique to simulate a known competition cycle route in Cape Town, South Africa. In an attempt to analyze the input stimulus, in this case the generated power output of each cyclist, eight subjects cycled for 40 min at a self-selected, comfortable pace on the first part of the simulated route. Our results indicate that this exercise input excites the musculocardiorespiratory system over a wide range of power outputs, both in terms of amplitude and frequency. This stimulus profile thereby complies with the fundamental requirement for nonlinear (physiological) systems analysis and identification. Through a computer simulation, we have devised a laboratory exercise protocol that not only is physiologically real but also overcomes the artificiality of most traditional laboratory exercise protocols.     

Scapular kinematics during unloaded and maximal loaded isokinetic concentric and eccentric shoulder flexion and extension movements

Characterization of scapular kinematics under demanding load conditions might aid to distinguish between physiological and clinically relevant alterations. Previous investigations focused only on submaximal external load situations. How scapular movement changes with maximal load remains unclear. Therefore, the present study aimed to evaluate 3D scapular kinematics during unloaded and maximal loaded shoulder flexion and extension. Twelve asymptomatic individuals performed shoulder flexion and extension movements under unloaded and maximal concentric and eccentric loaded isokinetic conditions. 3D scapular kinematics assessed with a motion capture system was analyzed for 20° intervals of humeral positions from 20° to 120° flexion. Repeated measures ANOVAs were used to evaluate kinematic differences between load conditions for scapular position angles, scapulohumeral rhythm and scapular motion extent. Increased scapular upward rotation was seen during shoulder flexion and extension as well as decreased posterior tilt and external rotation during eccentric and concentric arm descents of maximal loaded compared to unloaded conditions. Load effects were further seen for the scapulohumeral rhythm with greater scapular involvement at lower humeral positions and increased scapular motion extent under maximal loaded shoulder movements. With maximal load applied to the arm physiological scapular movement pattern are induced that may imply both impingement sparing and causing mechanisms.     

Effects of different movement directions on electromyography recorded from the shoulder muscles while passing the target positions

PurposeWe compared electromyography (EMG) recorded from the shoulder joint muscles in the same position for different movement directions.MethodsFifteen healthy subjects participated. They performed shoulder elevation from 0° to 120°, shoulder depression from 120° to 0°, shoulder horizontal adduction from ?15° to 105°, and shoulder horizontal abduction from 105° to ?15°. The target positions were 90° shoulder elevation in the 0°, 30°, 60°, and 90° planes (0°, 30°, 60°, and 90° positions). EMG signals were recorded from the supraspinatus (SSP) muscle by fine-wire electrodes. EMG signals from the infraspinatus (ISP), anterior deltoid, middle deltoid, and posterior deltoid muscles were recorded using active surface electrodes.ResultsDuring elevation and horizontal abduction, the SSP showed significantly higher activity than that shown during depression and during horizontal adduction in the 0°, 30°, and 60° positions. During elevation, the ISP showed significantly higher activity than during depression and during horizontal adduction in the 90° position. During horizontal abduction, the ISP showed significantly higher activity than during depression in the 90° position.ConclusionsWhen the movement tasks were performed in different movement directions at the same speed, each muscle showed characteristic activity.     

Nucleic Acid Structure Analysis: A Users Guide to a Collection of New Analysis Programs

Abstract

Common nomenclature describing the geometry of nucleic acid structures was established at a 1988 EMBO Workshop on DNA Curvature and Bending (Diekmann, S. (1988) J. Mol. Biol. 208, 787–791; Diekmann, S. (1989) The EMBO Journal 8, 1–4; Sarma, RH. (1988) J. Biomol. Structure & Dynamics 6, 391–395; Dickerson, R.E. (1989) J. Biomol. Structured Dynamics 6, 627–634; Dickerson, RE. et al. (1989)Nuc. Acids Res. 17, 1979–1803). We have subsequently developed and incorporated sophisticated mathematics in a computer program to calculate the parameters described by the guidelines. The program calculates all the local parameters relating complementary bases and neighboring base and base pairs in both Cartesian and helical coordinate frames. In addition, the main mathematical property requested by the EMBO guidelines—that the magnitude of the parameters be independent of strand or direction of measurement—is accomplished without the use of a midway coordinate frame for the rotational parameters. The mathematics preserve the physical intuition used in defining the parameters; in particular, the rotational parameters are true rotations based on a simple physical model (rotation at constant angular velocity for a unit amount of time), not Euler angles or angles between vectors and planes as is the case with other approaches. As a result, the mathematical equations are symmetric with the property that a 5° tilt is the same as a 5° roll or a 5° twist, except that the rotations take place about different axes. In other approaches, a 5° tilt can mean a different amount of net rotation from a 5° roll or a 5° twist. In addition, a great deal of flexibility is built into the program so that the user has control over the analysis, including the input format, the coordinate frame used for the base pairing relationship, the point about which the rotations are performed, and which geometric relationships are analyzed. While there is a great deal of flexibility, the program is easy to use. Interactive queries and user accessible files make the options in the program very convenient to tailor to individual needs. In addition, there is also a program that calculates bond lengths, valence angles, and torsion angles along the nucleic acid backbone, and within the sugar and base rings. Another program ‘learns’ the identities of the bond lengths, valence angles, and torsion angles that the user would like to determine. This last program is especially useful for calculating the hydrogen bonds between atoms in complementary strands as well as the unusual hydrogen bonds found in recently determined nucleic acid NMR structures or within protein/nucleic acid complexes.     

Particle-Image Velocimetry and Force Measurements of Leading-Edge Serrations on Owl-Based Wing Models

High-resolution Particle-Image Velocimetry （PIV） and time-resolved force measurements were performed to analyze the impact of the comb-like structure on the leading edge of barn owl wings on the flow field and overall aerodynamic performance. The Reynolds number was varied in the range of 40,000 to 120,000 and the range of angle of attack was 0° to 6° for the PIV and -15° to ＋20° for the force measurements to cover the full flight envelope of the owl. As a reference, a wind-tunnel model which possessed a geometry based on the shape of a typical barn owl wing without any owl-specific adaptations was built, and measurements were performed in the aforementioned Reynolds number and angle of attack： range. This clean wing model shows a separation bubble in the distal part of the wing at higher angles of attack. Two types of comb-like structures, i.e., artificial serrations, were manufactured to model the owl＇s leading edge with respect to its length, thickness, and material properties. The artificial structures were able to reduce the size of the separation region and additionally cause a more uniform size of the vortical structures shed by the separation bubble within the Reynolds number range investigated, resulting in stable gliding flight independent of the flight velocity. However, due to increased drag coefficients in conjunction with similar lift coefficients, the overall aerodynamic performance, i.e., lift-to-drag ratio is reduced for the serrated models. Nevertheless, especially at lower Reynolds numbers the stabilizing effect of the uniform vortex size outperforms the lower aerodynamic performance.     

Functionally interpretable local coordinate systems for the upper extremity using inertial & magnetic measurement systems

BackgroundInertial Magnetic Measurement Systems (IMMS) are becoming increasingly popular by allowing for measurements outside the motion laboratory. The latest models enable long term, accurate measurement of segment motion in terms of joint angles, if initial segment orientations can accurately be determined. The standard procedure for definition of segmental orientation is based on the measurement of positions of bony landmarks (BLM). However, IMMS do not deliver position information, so an alternative method to establish IMMS based, anatomically understandable segment orientations is proposed.MethodsFor five subjects, IMMS recordings were collected in a standard anatomical position for definition of static axes, and during a series of standardized motions for the estimation of kinematic axes of rotation. For all axes, the intra- and inter-individual dispersion was estimated. Subsequently, local coordinate systems (LCS) were constructed on the basis of the combination of IMMS axes with the lowest dispersion and compared with BLM based LCS.FindingsThe repeatability of the method appeared to be high; for every segment at least two axes could be determined with a dispersion of at most 3.8°. Comparison of IMMS based with BLM based LCS yielded compatible results for the thorax, but less compatible results for the humerus, forearm and hand, where differences in orientation rose to 17.2°.InterpretationAlthough different from the ‘gold standard’ BLM based LCS, IMMS based LCS can be constructed repeatable, enabling the estimation of segment orientations outside the laboratory.ConclusionsA procedure for the definition of local reference frames using IMMS is proposed.     

Dosimetric effect of uncorrected rotations in lung SBRT with stereotactic imaging guidance

PurposeTo evaluate the dosimetric impact of uncorrected rotations on the planning target volume (PTV) coverage for early stage non-small cell lung cancer patients treated with stereotactic body radiotherapy using Brainlab ExacTrac image guidance.MethodsTwenty-two patients were retrospectively selected. Two scenarios of uncorrected rotations were simulated with magnitude of 1°, 2°, 3° and 5°: (1) rotation around the treatment isocenter; and (2) roll and yaw rotations around a setup isocenter. The D₉₅ of PTV from recalculated dose on the rotated CT was compared to that from the clinical plan. A logistic regression model was used to predict the probability of dose differences between recalculated and original plans that are less than 2% based on the rotation angle, PTV volume, and distance between the treatment and setup isocenter.ResultsLogistic regression model showed the uncorrected isocentric rotations of up to 2.5° in all directions have negligible dosimetric impact. For non-isocentric rotations, a rotational error of 2° may cause significant under-dose of the PTV. Statistically significant (p < 0.05) parameters in the logistic regression model were angle for isocentric rotations, angle and distance for non-isocentric roll rotations, and angle, distance and the PTV volume for non-isocentric yaw rotations.ConclusionsThe severity of the dose deviations due to uncorrected rotations depends on the type and magnitude of the rotation, the volume of the PTV, and the distance between the treatment and setup isocenter, which should be taken into consideration when making clinical judgment of whether the rotational error could be ignored.     

Monte Carlo-based assessment of the trade-off between spatial resolution,field-of-view and scattered radiation in the variable resolution X-ray CT scanner

ObjectiveThe purpose of this work is to evaluate the impact of optimization of magnification on performance parameters of the variable resolution X-ray (VRX) CT scanner.MethodsA realistic model based on an actual VRX CT scanner was implemented in the GATE Monte Carlo simulation platform. To evaluate the influence of system magnification, spatial resolution, field-of-view (FOV) and scatter-to-primary ratio of the scanner were estimated for both fixed and optimum object magnification at each detector rotation angle. Comparison and inference between these performance parameters were performed angle by angle to determine appropriate object position at each opening half angle.ResultsOptimization of magnification resulted in a trade-off between spatial resolution and FOV of the scanner at opening half angles of 90°–12°, where the spatial resolution increased up to 50% and the scatter-to-primary ratio decreased from 4.8% to 3.8% at a detector angle of about 90° for the same FOV and X-ray energy spectrum. The disadvantage of magnification optimization at these angles is the significant reduction of the FOV (up to 50%). Moreover, magnification optimization was definitely beneficial for opening half angles below 12° improving the spatial resolution from 7.5 cy/mm to 20 cy/mm. Meanwhile, the FOV increased by more than 50% at these angles.ConclusionIt can be concluded that optimization of magnification is essential for opening half angles below 12°. For opening half angles between 90° and 12°, the VRX CT scanner magnification should be set according to the desired spatial resolution and FOV.     

Time/depth usage of Adélie penguins: an approach based on dive angles

Data on the swim speed, dive depth and feeding rates of three Adélie penguins (Pygoscelis adeliae) foraging in summer 1998/1999 in Adélie Land, Antarctica were collected using dorsally-mounted loggers, in tandem with oesophageal temperature sensors. Swim speed could be integrated, together with the rate of change of depth, to determine dive and return-to-surface angles. Overall, birds increased rates of change of depth during commuting phases so that dive angles were steeper in dives terminating at greater depths. Angles of descent and ascent during feeding dives were greater than during non-feeding dives. Variation in the descent angle over time of particular dives was generally less than 10°, but the angles of the ascent phases varied more widely. The importance of selecting the optimum descent and ascent angles with respect to prey exploitation, oxygen stores and time gained in the feeding area over the course of a dive by diving at a steeper angle is discussed.     

A method to perform spinal motion analysis from functional X-ray images

Identifying spinal instability is an important aim for proper surgical treatment. Analysis of functional X-ray images delivers measurements of the range of motion (RoM) and the center of rotation (CoR). In today's practice, CoR determination is often omitted, due to the lack of accurate methods. The aim of this work was to investigate the accuracy of a new analysis software (FXA?) based on an in vitro experiment.Six bovine spinal specimens (L3-4) were mounted in a robot (KR125, Kuka). CoRs were predefined by locking the robot actuator tool center point to the estimated position of the physiologic CoR and taking a baseline X-ray. Specimens were deflected to various RoM_preset flexion/extension angles about the CoR_preset. Lateral functional radiographs were acquired and specimen movements were recorded using an optical motion tracking system (Optotrak Certus). RoM and CoR errors were calculated from presets for both methods. Prior to the experiment, the FXA? software was verified with artificially generated images.For the artificial images, FXA? yielded a mean RoM-error of 0.01±0.03° (bias±standard deviation). In the experiment, RoM-error of the FXA?-software (deviation from presets) was 0.04±0.13°, and 0.10±0.16° for the Optotrak, respectively. Both correlated with 0.998 (p<0.001). For RoM<1.0°, FXA? determined CoR positions with a bias>20 mm. This bias progressively decreased from RoM=1° (bias=6.0 mm) to RoM=9° (bias<1.5 mm).Under the assumption that CoR location variances <5 mm are clinically irrelevant on the lumbar spine, the FXA? method can accurately determine CoRs for RoMs>1°. Utilizing FXA?, polysegmental RoMs, CoRs and implant migration measurements could be performed in daily practice.