CFD Analysis of the Survivability of a Square Rigged Sailing Vessel

 Emory Sullivan
 8 months ago
 Views:
Transcription
1 Engineering Applications of Computational Fluid Mechanics ISSN: (Print) X (Online) Journal homepage: CFD Analysis of the Survivability of a Square Rigged Sailing Vessel William C. Lasher & Logan S. Flaherty To cite this article: William C. Lasher & Logan S. Flaherty (2009) CFD Analysis of the Survivability of a SquareRigged Sailing Vessel, Engineering Applications of Computational Fluid Mechanics, 3:1, 7183, DOI: / To link to this article: Copyright 2009 Taylor and Francis Group LLC Published online: 19 Nov Submit your article to this journal Article views: 343 View related articles Citing articles: 1 View citing articles Full Terms & Conditions of access and use can be found at
2 Engineering Applications of Computational Fluid Mechanics Vol. 3, No. 1, pp (2009) CFD ANALYSIS OF THE SURVIVABILITY OF A SQUARERIGGED SAILING VESSEL William C. Lasher* and Logan S. Flaherty Penn State Erie, The Behrend College, School of Engineering, 5101 Jordan Road, Erie, PA , USA * (Corresponding Author) ABSTRACT: A CFDbased model has been developed for predicting the aerodynamic forces on the rig and sails of the U.S. Brig Niagara. Wind tunnel tests and fullscale experiments were performed to validate the model. The model was then used to predict the heel angle for different wind speeds and sail configurations, with a focus on determining whether the ship would be able to survive a storm. The results show that the Captain s practice of not allowing the ship to heel greater than 10 degrees is prudent, and that one of the sail configurations examined should be safe for any reasonably expected wind conditions on the Great Lakes. These results are useful for providing guidance to the Niagara s officers regarding survivability of the ship. Keywords: sailing, stability, survivability, squarerigged 1. INTRODUCTION For several hundred years wooden traditionallyrigged sailing ships were the dominant form of ocean transportation. Even though fuel power replaced sail as the primary mode of propulsion long ago, there are still an estimated 300 traditionallyrigged sailing vessels worldwide. These vessels, which are also known as tall ships or squarerigged ships, are used today primarily for historical reasons. One example of such a vessel is the U.S. Brig Niagara (Fig. 1). Niagara is a reconstruction of a sailing warship built for the War of 1812 between the U.S. and England. She has a dual mission of historic interpretation as well as skills preservation. Within this historic interpretation mission there are two main tracks. The ship is interpreted as representing a warship from an important battle in the War of 1812, a story of naval history. The ship also represents seafaring as an industrial process. The other, and dominant, aspect of the mission is a skills preservation program carried out through the active sailing of the ship on 60 to 75 days each summer. Niagara sails today with a crew of up to forty people, half professional, half trainee. As a comparison, during the time that the Niagara was commissioned as a U.S. Naval Vessel, she carried a crew of 155. Due to the difference in both size and experience, today s crew is less able to respond to sudden changes in weather conditions, such as squalls that can quickly appear on the Great Lakes. As a result, they are very concerned about survivability that is, how much sail she can carry and still remains upright during a storm. One might think that a question such as what it would take to capsize a squarerigged sailing ship would have been answered in the historical record; after all, these vessels were around for hundreds of years and sailed through almost every imaginable condition. There is, in fact, a small body of literature describing practical seamanship in sailing vessels, but it is largely anecdotal and the information is spotty (for example, Harland and Myers, 1984). In many cases the information only applies to specific vessels; in some cases it is simply wrong. Fig. 1 Niagara under partial sail. Received: 20 May 2008; Revised: 31 Jul. 2008; Accepted: 7 Aug
3 There is a significant body of published work (both experimental and computational) on sail aerodynamics, most of which applies to sailboat racing. For example, there are several wind tunnels that specialize in sail aerodynamics research. Most of these tunnels have capabilities that are important to studying sail aerodynamics, such as the tunnel at the Wolfson Unit at the Southampton University (Deakin, 1991). This tunnel has been involved in wind tunnel tests on sailing vessels since the late 1960 s, and is capable of producing an atmospheric boundary layer, which will be described later. Deakin investigated the wind heeling characteristics of sailing vessels in both steady winds and in gusts a problem similar to that investigated here, but under different conditions. Flay (1996) described the development of a twisted flow wind tunnel at the Yacht Research Unit at the University of Auckland. Twisted flow occurs due to the combination of the atmospheric boundary layer and the motion of the boat relative to the wind. This tunnel has been used for a number of studies for grandprix racing yachts, such as those described by Fallow (1996). Ranzenbach and Mairs (1997) implemented advanced blockage correction techniques at the Glenn L. Martin Wind Tunnel at the University of Maryland and used these techniques for a variety of studies, such as the investigation of forces on downwind (sailing in the direction of the wind) sails (Ranzenbach and Mairs, 1999). Blockage correction is especially important when the flow over the sails is separated, which occurs when a vessel is sailing downwind, or in conditions that are examined in this paper. Computational techniques have also been used for sail aerodynamics research for some time. Milgram (1968) first employed panel techniques to sail aerodynamics with some success, and Charvet and Huberson (1992) used vortex particle methods on this problem. These inviscid techniques work reasonably well for upwind sailing (sailing towards the wind), where the sails are flat and separation is at a minimum; however, they cannot be applied to downwind sails due to large areas of separated flow (Richards, 1997). Hedges, Richards and Mallinson (1996) were among the first to used ReynoldsAveraged Navier Stokes (RANS) simulation for downwind sail aerodynamics. Lasher and Richards (2007) compared the results from RANS simulations to wind tunnel data for spinnakers (a special type of downwind sail) in an untwisted atmospheric boundary layer, and in general the lift and drag coefficients were predicted with reasonable accuracy. Sail forces have also been found through fullscale testing. One of the most noted publications on fullscale testing is by Gerritsma, Kerwin and Moeyes (1975), who used equilibrium concepts and data from model tests of the hull to estimate the sail forces. Recognizing the limitations of this approach, Milgram, Peters and Eckhouse (1993) built a sailing dynamometer, which is a specially constructed boat where the mast and rigging are attached to a rigid frame that is structurally separate from the hull. This allows for direct measurement of the sail forces without having to estimate the hull hydrodynamic forces. Other sailing dynamometers have been built as reported by Masuyama and Fukasawa (1997) and by Hochkirch and Brandt (1999). While a sailing dynamometer is a major improvement over the forcebalance approach, there are still significant problems determining accurate force coefficients. Due to problems associated with accurately measuring boat speed, wind speed near the boat (the sails influence the incoming free stream), sail shape, and the inherent fluctuations in both wind speed and direction, there is a tremendous amount of scatter in the data. For example, Milgram, Peters and Eckhouse (1993) showed drag coefficient vs. lift coefficient for upwind sailing. At a lift coefficient of 1.0, the mean experimental drag coefficient is near 0.19, with a standard deviation of and a range from 0.14 to While the stability limits of squarerigged vessels can be found through onthewater experimentation, such an approach is risky, as it may result in the loss of the ship. The value of a sciencebased approach to predicting these limits, such as those described above, is clear. As a result of this need, two teams of researchers began investigating the stability of two different squarerigged vessels (Lasher et al., 2007; Miles et al., 2007). Lasher et al. (2007) developed a CFDbased model using Fluent 6.3 (Fluent, 2006) for predicting the stability of the Niagara, and Miles et al. (2007) began an experimental program to record wind and stability data during sailing conditions on the Pride of Baltimore II. The present work is an outgrowth of the computational work of Lasher et al. (2007), and is part of a larger program on squarerigged vessel stability. The two teams mentioned above have begun collaborating on a joint computational/experimental program to develop and validate models of squarerigged ship stability, and use these models to provide guidance to naval architects for design purposes, 72
4 as well as ship operators for decision making purposes. This paper focuses on one aspect of the stability question survival in a storm. During heavy weather the ship s captain may decide to go hoveto, which is a state where just enough sail is hoisted to keep the ship pointed in the desired direction, and forward motion is minimized. In the following sections an explanation of sailing aerodynamics and terminology will first be presented. The mathematical model used to predict stability will then be explained, including a discussion of model validation. The particular cases (or hoveto configurations) that were analyzed will be described, followed by a presentation of results and discussion. 2. A PRIMER ON SAILING AND SQUARE RIG TERMINOLOGY When a ship is sailing, the wind as seen from the ship s reference is a vector combination of the true (or actual) wind and the wind generated by the forward motion of the ship. This wind is known as the apparent wind, and the angle that the wind makes with the bow (or front) of the ship is known as the apparent wind angle (AWA). The sails act as airfoils, and the apparent wind generates aerodynamic lift and drag forces as shown in Fig. 2. These forces can be resolved into drive and side forces. The drive force propels the ship forward, and the side force has two effects it makes the ship slide sideways, and also causes a heeling moment, which makes the ship heel, or tip sideways. The side force is resisted by an equal and opposite hydrodynamic force on the ship s hull, and the heeling moment is resisted by a righting moment, which is generated by a shift in the ship s center of buoyancy. The main challenge in assessing stability is to determine the heeling moment due to the sails and rig (the rig consists of spars, such as the masts and yard arms; and rigging, which are lines that hold the masts up). The righting moment of the hull as a function of heel angle (or tipping angle) can be determined by the naval architect through hydrostatic analysis of the ship s hull shape and an onthewater inclining experiment, from which the ship s vertical center of gravity can be found. The actual heel angle can then be determined by setting the heeling moment equal to the righting moment. The sailor s primary goal is to select the combination of sails and sail settings (angle of attack of the sails) to produce the desired effect, which is normally to maximize the drive force and minimize the side force and heeling moment. This is a complex process which is made even more complicated in Niagara s case by the large number of sails and possible combinations. Fig. 3 shows a sail plan for Niagara, along with the names of the different sails. The sails shown in white are square sails; those shown in grey are known as fore and aft sails. To the novice the names of these sails appear to be a foreign language, but to an experienced sailor they have very specific meanings. For simplicity, references to specific sail names will generally be followed by the letter shown in the legend for Fig. 3. A. Spanker H. Main Topgallant Staysail B. Mainsail I. Fore Topsail C. Main Topsail J. Main Topmast Staysail D. Main Topgallant K. Foresail E. Main Royal L. Foretopmast Staysail F. Fore Royal M. Jib G. Fore Topgallant N. Flying jib Fig. 3 Niagara s sail plan. Fig. 2 Aerodynamic forces acting on a sailing vessel. AWA is the apparent wind angle. The angle of attack of each of these sails to the wind can be adjusted on the ship by various lines. When sailing the angle of the sail is referenced to 73
5 the centerline of the ship rather than to the wind, and on squarerigged vessels this angle is known as the bracing angle (the term bracing angle technically only applies to the squarerigged sails, but for simplicity the term will be used for all sails). For example, if a sail is set so that the chord is parallel to the centerline of the ship, then the bracing angle would be zero; if it were eased out square to the ship its bracing angle would be 90. In Fig. 2 the bracing angle of the sail is 35 degrees. 3. MATHEMATICAL MODEL 3.1 Overview of the aerodynamic model Due to geometric complexity, the computation of the heeling moment was broken into two parts the moment due to the sails, and the moment due to the rig (spars and rigging). The moment due to the sails was determined through CFD computation, and the moment due to the spars and rigging was determined using a simple drag coefficient approach. Interaction between the sails, spars and rigging was neglected, and the two moments were summed to find the total heeling moment. This moment was then used in conjunction with the righting moment curve to find the heel angle. All analyses were done using a reference wind speed of 40 knots (20.6 m/sec), and the resulting heeling moments were then scaled to different wind speeds assuming that the moments varied with the square of the wind velocity. While wind speeds of 40 knots are not common on the Great Lakes, they do occur, especially during thunderstorms. It is also known that ships sail in an atmospheric boundary layer, so a velocity profile appropriate for the open ocean was used (Kerwin, 1978): V( y) = 1.96 ln(1000y+ 1) (1) In Eq. (1), y is the vertical distance from the water surface measured in meters and V the wind velocity in meters per second. This produces the reference wind speed of 40 knots at the head of the mainmast on Niagara, which is where the anemometer is located. Since the wind seen by the anemometer can be influenced by the presence of the ship and sails, the CFD analysis was used to determine what velocity would be measured by the anemometer when Eq. (1) was used as an inlet boundary condition. In the present analysis this was very close to the reference velocity (typically within 1%), since all of the sails analyzed were low on the rig and far from the anemometer. The forces on each of the spars and rigging were calculated using a drag coefficient of 1.13 (Kerwin, 1978). The lateral components of these forces were then multiplied by the distance to the center of lateral resistance of underwater portion of the ship s hull, which was assumed to be 40% of the ship s draft below the waterline. The value of 40% is used because this represents the location of the center of lift on an ellipticallyloaded airfoil. While the hull is not acting strictly as an airfoil this is a reasonable approximation, and any error introduced by this assumption will be small due to the large height of the rig compared to the draft of the ship. The calculated heeling moment was then added to the heeling moment due to the sails as predicted by CFD. 3.2 CFD model for sails Flow equations The equations governing incompressible turbulent flow are the wellknown ReynoldsAveraged Navier Stokes (RANS) equations, written in Cartesian tensor notation as: Du 1 u i p = + ν i uu ' ' Dt ρ x x x i j i j j and (2) u i = 0 (3) x i Equation (2) is derived from the momentum equation, where the velocities u i are the timeaveraged velocities and u i ' are the fluctuating velocities caused by the turbulence, so the equation represents a timeaveraged balance of momentum. Eq. (3) is the continuity equation, representing conservation of mass. The above equations do not form a closed set (i.e., there are not enough equations for the number of unknowns) due to the last term in Eq. (2), which is called the Reynolds stress. Various models for this term have been proposed. Lasher and Sonnenmeier (2008) performed simulations using six different turbulence models on a series of sail models that had been previously tested in a wind tunnel, and compared the simulations to experimental data. They showed that the Realizable kε model (developed by Shih et al., 1995) was the best of the six models tested in terms of accuracy and stability. This model was therefore used for the present work. The Realizable kε model is based on the Boussinesq eddyviscosity hypothesis, which 74
6 relates the Reynolds stress to the mean flow strain rate tensor s ij, the turbulent kinetic energy k, and a turbulent viscosity v t by: ' ' 2 uu = kδ + 2ν s i j 3 ij t ij and 1 u u = j s i + ij 2 x x j i (4) (5) If the turbulent kinetic energy k and turbulent viscosity v t can be defined as a function of the mean flow variables u i and p, then the set of Eqs. (2) through (5) can be solved. In the Realizable kε model the turbulent viscosity is determined from the turbulent kinetic energy k and turbulent dissipation ε from 2 ν = c k t μ ε (6) The turbulence quantities k and ε are determined from modeled differential Eqs. (7) and (8): Dk ν = ν+ t k + G ε (7) Dt x σ x j k j Dε ν ε 2 = ν+ t ε + DSε D (8) Dt x σ x 1 2 k + νε j ε j G in Eq. (7) is the generation rate of turbulent kinetic energy, and is given by u u =ν j u G i + i t x x x j i j (9) The term S in Eq. (8) is based on the magnitude of the mean vorticity and is given by S = 2Ω Ω ij ij (10) where Ω ij is the mean rateofrotation tensor given by 1 u u Ω = i j ij 2 x x j i (11) The constants in this model generally take on the following values (Fluent, 2006): σ k = 1.00, σ ε = 1.20, D 2 = The coefficient D 1 in Eq. (8) is a function of S, k and ε, and the coefficient c μ in Eq. (6) is a function of the turbulent kinetic energy and dissipation, as well as the mean strain rate s ij Domain and boundary conditions The ship was placed in the center of a rectangular domain measuring 305 m 183 m (with the 305 m dimension in the streamwise direction), and a vertical height of 137 m (Fig. 4). The overall length of the ship is approximately 60 m, and the top of the mainmast is approximately 37 m above the water. To verify that the domain was large enough to not influence the results, a sample simulation was performed using a domain that was 33% smaller in each direction. The resulting predicted heeling moment was 2% larger, indicating that the domain size is appropriate. Z Fig. 4 Y X Computational domain shown in isometric perspective. Wind is in the positive xdirection. Equation (1) was used as a boundary condition on velocity at the domain inlet. Turbulence boundary conditions at the inlet were determined by comparing Eq. (1) to a generalized model for atmospheric boundary layers (Cook, 1985): ( ) u * y V y = ln κ y0 (12) In Eq. (12), u * is the friction velocity, κ is von Karman s constant (taken as 0.41), and y 0 is the surface roughness. Eqs. (1) and (12) produce identical velocity profiles but are in different forms. Equation (1) assumes an average surface roughness of m that is independent of wind velocity, whereas Eq. (12) allows for surface roughness that varies with wind speed. At the reference velocity selected for this paper, the friction velocity determined by comparing the two equations is m/s. The turbulent kinetic energy for an atmospheric boundary layer can be found from Eq. (13) 2 * k = 5.5u (13) 75
7 and the turbulent dissipation rate ε can be found from Eq. (14). u ε = * (14) κ( y+ y ) 0 Equations (13) and (14) were determined by a balance of the standard kε model using Eq. (12) as the velocity profile (Richards and Hoxey, 1993). A wall boundary condition with wall functions was used for the water s surface. With this boundary condition one specifies the velocity at the wall (zero in this case) as well as the turbulent shear stress, which determines the nearwall dissipation rate and kinetic energy. The shear stress on the water s surface τ w is given by Eq. (15): 2 τ w = ρu * (15) In Eq. 15, ρ is the density of the air. The surface roughness in Fluent is defined by the product of two constants: a roughness height K s and a roughness constant C s. A comparison of Fluent s roughness model with Eq. (12) shows that the product of the two constants must be approximately 10 times the surface roughness y 0 used in Eq. (12) to produce the same velocity profile, so the default constant of 0.5 was used for C s along with a surface roughness of 0.02 m for K s. A moving wall boundary condition was also used for the upper surface, with the wall velocity set to the velocity computed from Eq. (1). A symmetry condition was applied to the sides of the domain, and a uniform reference pressure of zero was applied to the outlet. Since the intent of being hoveto is to produce no motion, the velocity of the ship was set to zero. All of the sail shapes used in the CFD model were fixed (that is, it was assumed that the sails were rigid, with the only variable being the bracing angle). The sail shape used for the square sails was based on photographs taken from the deck. Since these sails are the largest and high in the rig, they are the most critical in terms of matching predicted performance with actual performance. Fore and aft sail shapes were based both on photographic evidence as well as commonly accepted values for the actual camber. The crosssection shape used for all sails was a circular arc Solution procedure The transport equations described above were solved using the finite volume method with the SIMPLE technique for determining pressure (Patankar, 1980) as implemented in Fluent (Fluent, 2006). The convective terms were calculated using a QUICKtype scheme (Leonard, 1979), which is secondorder accurate as implemented in Fluent for unstructured grids. Both the normalized residuals and the aerodynamic lift and drag coefficients were monitored, and the solution was iterated until either the normalized residuals were less than , or the force coefficients changed by less than 0.5% after several time steps. The number of time steps required for convergence depended on the angle of attack of the sails. The solutions took between 2 to 6 hours on a 3 GHz Pentium Processor. Two grids were used to assess numerical error due to grid spacing. The coarse grid used an initial grid spacing of 0.6 m on the sail surfaces, increasing at a geometric expansion ratio of 1.08, with a maximum grid spacing of 9 m. The fine grid used initial and maximum grid spacings of 0.3 m and 4.5 m, respectively, with the same expansion ratio. It was not possible to obtain a converged solution using the fine grid for all of the angles of attack due to the unsteady nature of separated flow, so the coarse grid was used for all subsequent computations to provide a consistent basis for comparison. Using Richardson extrapolation (Roache, 1998) the error in predicted force coefficients on the coarse grid was estimated to be approximately 5%. A fundamental problem in computing separated flows such as those seen in the present work is that they are inherently unsteady. Technically, the simulations should be transient with the results timeaveraged to fully capture largescale coherent structures such as vortex shedding. This would require a much finer mesh than is practical for this type of problem, especially given the uncertainty in the computed rigging drag, the actual sail shape, and the actual wind velocity profile (which can vary greatly, especially during storm conditions). In addition, the simulations technically should be performed using an unsteady turbulence model, which is currently not available in Fluent. The grid parameters used in the present work probably underresolve the flow; however, a detailed study on the simulation of separated flow over sails by Lasher and Sonnenmeier (2008) showed that these grid parameters produced results which reasonably agreed with wind tunnel tests. When using wall functions the distance of the first point from the sail surface is important. This distance is measured using the variable y + =u τ y/υ, where u τ is the shear velocity at the sail surface 76
8 (τ wall /ρ) 0.5. It is generally recommended that y + for the first grid points should be no less than 30 but as close to 30 as possible. In the present case most of the first points had y + values between 30 and 120, with a few less than 30 and some as high as 200 (depending on the sail shape and angle of attack). Given the highly complex nature of this flow this is a reasonable range of values for y +, as it is not practical to have a constant value. Further, Lasher and Richards (2007) showed that the predicted force coefficients for sail aerodynamics are not sensitive to either the nearwall grid spacing or the type of wall function used. 3.3 Heel calculation The righting moment for the ship (Fig. 5) was determined by the naval architect based on hydrostatics and the results of an inclining experiment. A 5thorder polynomial was fit to the data. The maximum righting moment for Niagara occurs at 30 degrees, and the ship would capsize if the heel exceeds this amount. As mentioned above, the total heeling moment was calculated by adding the heeling moment due to the sails to that due to the spars and rigging. Since all calculations were done at a reference velocity of 40 knots, the heeling moment was scaled to other wind speeds by assuming that it was proportional to the square of the velocity. creates a downward force and additional heeling moment. The magnitude of this additional moment is uncertain; however, Deakin (1991) showed in wind tunnel tests that the heeling moment for sailing yachts decreases by cos 1.3 (θ). At 30 degrees heel the difference between using cos(θ) and cos 1.3 (θ) is about 4%, so cos(θ) was used to make the model conservative. The reduced heeling moment was set equal to the righting moment, and the resulting nonlinear equation was solved for the heel angle using NewtonRhapson. 4. MODEL VALIDATION 4.1 Wind tunnel tests Since the square sails are a major contributor to both driving force and heeling moment, it is especially important to know how accurate the CFD simulations are for these sails. To assess the accuracy of the simulations, a rigid scale model of a topsail (C or I) was built and tested in a wind tunnel using the procedure described by Lasher et al. (2003). Lift and drag coefficients (the force coefficients perpendicular and parallel to the free stream) were measured at angles of attack (defined as the angle between the incoming free stream and the chord of the sail) ranging from 20 to 90 in 10 increments. CFD analysis was also performed for this sail using the conditions in the wind tunnel. The experimental and predicted force coefficients are shown in Fig. 6. Fig. 5 Niagara s righting moment determined from inclining experiments and hydrostatics. All calculations were performed with the ship at zero degrees heel, so it was necessary to correct the heeling moment due to the reduction in projected area as the ship heels. In the present work the heeling moment was reduced by cos(θ), where θ is the heel angle. It can be argued that the heeling moment should decrease by cos 2 (θ), since both the projected area and the heeling arm are reduced; however, as the ship heels the sails become tilted horizontally to the free stream, causing the flow to be deflected upwards, which Fig. 6 Comparison of wind tunnel data and results from CFD simulation for a single rigid topsail. The predicted force coefficients are always larger than the experimental values. Over the range of angles tested, the average predicted lift and drag are 11% and 15% greater than the experimentally measured lift and drag, respectively. At 50, the predicted lift and drag are 24% and 28% higher than the measured lift and drag, respectively. These results suggest that the heeling moment predicted by the model will be higher than what is 77
9 experienced by the ship, depending on the influence of the other factors. This means that the model will most likely provide a conservative estimate of the heel angle; that is, the actual heel angle is likely to be less than that predicted by the model. 4.2 Fullscale validation It is extremely difficult to obtain accurate heel measurements on the ship to compare with predictions from the model. One must have a relatively long period of steady wind conditions with low waves to allow the ship to reach equilibrium. The wind speed must also be high enough to generate a measurable heel angle (on the order of 3 5 ) and/or consistent ship speed, yet not so high that the master becomes concerned with the ship s and crew s safety. Finally, it is difficult to measure the actual sail trim, so this must be estimated, introducing additional error. Due to both Niagara s schedule and the wind conditions on the Great Lakes, only a few measurements were taken that are useful for validation; however, these measurements do provide some illuminating information. Three different cases for validating the predicted heel angle are presented. In each of these cases the wind speed was measured using the anemometer and the heel angle was measured using an inclinometer. Sail trim was estimated by observation. CFD simulations were performed to match the recorded conditions and sail trim. In the first case the ship was hoveto with the foretopmast staysail (L) aback (set backwards) and single reefed (reduced area) topsails (C, I) with the fore topsail aback (see Fig. 7). The wind was 17 knots on the beam (i.e., at an AWA of 90 degrees). The heel angle was measured to be 1.8 degrees using the dinghy davits on the photograph and the water surface as a reference. In the second case the ship was sailing at about 60 degrees AWA in knots of wind. The ship was carrying the jib (M), foretopmast staysail (L), maintopmast staysail (J), both topsails (C, I), and a single reef in the spanker (A). The measured heel angle was approximately 4 degrees. In the third case the ship had been sailing close to the wind when a front came through. She was carrying both topsails (C, I) with a single reef, the jib (M), foretopmast staysail (L), and maintopmast staysail (J). The sails were trimmed appropriately for the apparent wind direction when the master had the ship turned away from the wind. When the wind was on the beam, the heel angle was approximately 5 degrees, and the wind speed was approximately knots (due to the rapidly changing situation there is some uncertainty as to what the actual wind speed was at the time the heel angle was recorded). Fig. 7 Niagara practices going hoveto in moderate winds. Notice the two topsails (square sails) are pointed in opposite directions; the forwardmost topsail is backed. The measured and predicted heel angles for all three cases are shown in Fig. 8. Error bars of ± 0.5 are shown on the measured values; these represent an estimate of the error in measuring the heel angle due to the resolution of the inclinometer as well as the fact that the ship may have had an initial unmeasured heel. In all three cases the predicted heel angle is greater than the measured heel angle by degrees. For cases 2 and 3, the predicted heel angle is approximately 15% higher than what was measured. The heel angle for case 1 is too small to calculate a meaningful percent error. Fig. 8 Comparison of predicted and measured heel angles. 78
10 While these results cannot be considered definitive, they do indicate that the model is reasonably accurate. The overprediction of heel angle is consistent with what was found comparing the wind tunnel data with CFD simulations, and suggests that the other assumptions made in the model are not contributing significantly to the total error. As mentioned in the introduction, a large scale program is currently under way to measure and record wind conditions and heel angle on both Niagara and the Pride of Baltimore II. These data will be used for further validation of the mathematical model. 5. DESCRIPTION OF HOVETO CONFIGURATIONS A ship will go hoveto when the wind is too high to maneuver safely. This process consists of setting a small number of sails, sometimes reversed (or backed), with the rudder positioned to keep the ship pointed in the desired direction relative to the wind. Different hoveto configurations can be used in different conditions, as explained below. The goal of being hoveto is to minimize drift of the ship while riding out the storm with minimum strain on the vessel. The Captain of the Niagara was most interested in four specific hoveto configurations, which are summarized in Table 1 and shown in Fig. 9. A brief note on each of these configurations is provided below: Configuration 1 Since the wind is on the beam and waves will cause the ship to roll, this is only used in light to moderate weather for man overboard recovery, fire drills, abandon ship drills, etc. The Captain is interested in knowing the upper wind limit of this tactic. Configuration 2 Niagara has resorted to this on several occasions but not in wind velocities of over 30 knots, or wave heights greater than 5 feet. This is historically the preferred method of heavingto in a square rigged ship, and is considered to be safe to ride out in a major storm. The Captain wants to know whether this configuration would be safe for Niagara in very heavy winds. Configuration 3 This has been tried on Niagara but found to be uncomfortable, primarily because it was undercanvassed (i.e., did not provide enough force) in the winds in which it was tested (25 knots). The Captain wants to know whether this approach might work in higher winds. Configuration 4 This would be for extreme conditions and used in combination with a sea anchor, which is a drogue that reduces drift of the ship and helps keep her headed in the proper direction. 45 degrees AWA is the most likely orientation, but 70 degrees is the worst case (i.e., maximum heeling moment), so both of these apparent wind angles were analyzed. Table 1 Summary of hoveto sail configurations. Configuration AWA, degrees Sails Sail bracing angle, degrees (a) 70 (b) 45 Fore topsail (I) Fore topgallant (G) Main topsail (C) Main topgallant (D) Spanker (A) Jib (M) Fore topmast staysail (L) Fore topsail (I), double reef Main topsail (C), double reef Fore topmast staysail (L) Spanker (A), double reef Fore topmast staysail (L) Main topmast staysail (J) Backed (30) Backed (30) Backed (15) 5 Backed (15) 20 Spanker (A), triple reef 5 79
11 Configuration 1 Y Z Y X Z X Configuration 2 Z Y X Y Z X Configuration 3 Z Y X Configuration 4a (Configuration 4b is not shown; same geometry at different apparent wind angle) Z Y X Fig. 9 The four hoveto sail configurations. (The apparent wind is in the direction of the postitive axis.) 6. RESULTS AND DISCUSSION The forces and moments calculated by the CFD simulations for the sails and the force coefficient calculations for the spars and rigging are summarized in Table 2. All of the results in Table 2 are computed at a reference velocity of 40 knots. A number of conclusions can be drawn from this table: As expected, the total heeling moment decreases with decreasing sail area. 80
12 With the exception of configuration 2, the driving force for all cases is quite small compared to the side force. This is generally the desired effect, as the goal in heavingto is to minimize forward motion. The large lateral area of the underwater portion of the hull creates a significant reactive force, which will restrict sideways motion due to the side force. The larger relative driving force for configuration 2 is due to the fact that the two topsails (C, I) are set in the normal driving position, and these sails are large compared to the backed fore topmast staysail (L). The contribution to the total heeling moment from the spars and rigging increases from 17% for configuration 1 to 87% for configuration 4(a). This means that the accuracy of the model for calculating the heeling moment due to the spars and rigging becomes more important as sail is reduced; for configurations 4(a) and (b) the ultimate stability of the vessel is almost entirely determined by windage on the rig. This suggests that further fullscale experimental verification of the rig model is critical for determining the safety of the ship in these configurations. The heeling moment due to the spars and rigging is not the same for all configurations because the apparent wind angle is different, and this changes the direction of the drag vector. In addition, there are differences in the exposed area of the spars, as some of the spars have different bracing angles and some of them are lowered, depending on the configuration. Table 2 Forces and heeling moments at a reference velocity of 40 knots. Force or moment Configuration number (a) 4(b) Driving force, m 2,331 28,263 7,339 1,139 1,321 Side force, m 178,583 63,117 57,806 11,133 11,947 Heel moment from sails, Nm 2,830, , ,774 77,862 83,764 Heel moment from rig, Nm 584, , , , ,074 Total heeling moment, Nm 3,415,409 1,481,501 1, , ,838 Fig. 10 Heel angle vs. wind speed for all four configurations. The heeling moments from Table 2 were scaled by the square of the wind velocity (i.e., it was assumed that the heeling moment coefficient was constant), and this scaled moment was used to calculate the heel angle as described earlier. The predicted heeling angles vs. wind speed are shown in Fig. 10. It should be noted that the validation discussed earlier indicates that the model overpredicts the heel angle by approximately 15%, which provides a small safety factor. Since the heeling moment varies with the square of the wind speed, this reduces the underprediction of wind speed that would cause a particular heel angle to about 7%. For example, the wind speed predicted by the model to produce a heel angle of 30 degrees for configuration 1 is about 39 knots; a 15% overprediction of the heel angle means that the actual wind speed required to produce this heel angle would be approximately 42 knots a difference on the order of the accuracy of the anemometer. Also, since the validation was done using wind tunnel testing of a sail only (no rig), and the fullscale validation was performed using configurations with a large amount of sail area, there is a higher uncertainty in the results that have smaller sail area configurations. This uncertainty needs to be further explored by fullscale onthewater testing. 81
13 In addition to the general trends, the captain of Niagara is interested in two specific wind velocities those which create heel angles of 10 and 30 degrees. Ten degrees is the comfortable limit of heel due to the importance of safety, the captain does not want to exceed this heel angle as a matter of practice. At 30 degrees the ship will capsize and be lost, so the captain does not want to get even close to that limit. An examination of Fig. 10 shows the wisdom of not wanting to exceed 10 degrees of heel. For example, for configuration 1 the wind speed required to produce a heel angle of 10 degrees is approximately 25 knots, and the maximum wind speed for survivability of the vessel is 39 knots. In wind conditions with mean wind speeds of 25 knots it is not uncommon to have temporary gusts in excess of 40 knots. While it will take time for the ship to respond to a gust of short duration, the possibility that the ship will capsize must be avoided. Wind speeds in excess of 60 knots are infrequent on the Great Lakes, so it is doubtful that the Captain would have to be concerned about configuration 4 while sailing in home waters. The only concern in this case would be if the ship were to be on an extended ocean voyage, which has not been attempted in her history. Configuration 4 can be considered safe for essentially any weather condition that the ship is likely to experience. It was previously mentioned that configuration 3 was uncomfortable because it was undercanvassed in winds of 25 knots. Figure 10 shows that this configuration should produce a heel angle of approximately 4 degrees at 25 knots, and is better suited for wind speeds of knots, thus supporting the Captain s original observations. Finally, the Captain has noted that the wind speeds that would capsize the ship are somewhat lower than what he expected. Niagara does not have a very deep draft due to the fact that it had to be moved over a sandbar to get from the harbor where it was built into Lake Erie. As a result of this low draft the ship does not have the same stability as other ships of that era with similar size and mission. The low relative stability of Niagara makes the results of this work even more important. 7. CONCLUSIONS A CFDbased model for predicting stability for the U.S. Brig Niagara was developed. Results from the model were shown to be in good agreement with fullscale experimental data. Some of the more important results from this analysis are: The model overpredicts the heel angle of the ship by about 15%. Although the data set is small and further verification is warranted, these results are consistent with a comparison between the CFD simulations of the sails and wind tunnel tests of a small, rigid sail. The results provide useful guidance for the crew that confirms some of their standard practices. These results will help the Captain make decisions such as whether is it better, when caught on a lee shore (a shore closely downwind from the ship), to try to sail away from the shore and risk capsize, or to head directly to the shore and risk beaching. The relative importance of the heeling moment created by the rig for the smaller sail area configurations indicates that experimental verification of the rig model is warranted, so that the accuracy of the predictions for these cases can be determined. A program is currently in progress to further develop and validate this model, and will be an important step towards providing naval architects and sailors with methodologies for designing and operating traditionallyrigged vessels safely. ACKNOWLEDGEMENTS This work was supported by grants from the Behrend College Undergraduate Student Research Programs. The authors also wish to express their thanks to sailmaker David Bierig for originally suggesting the line of inquiry, and to Niagara Captains Walter Rybka and Wesley Heerssen for their assistance in collecting full scale experimental data. REFERENCES 1. Charvet T, Huberson S (1992). Numerical Calculation of the Flow around Sails. European Journal of Mechanics, B/Fluids 11(5): Cook NJ (1985). The Designer's Guide to Wind Loading of Building Structures, Part 1. Butterworth. 3. Deakin B (1991). Model Test Techniques Developed to Investigate the Wind Heeling Characteristics of Sailing Vessels and their Response to Gusts. Proceedings of the Tenth 82
14 Chesapeake Sailing Yacht Symposium, February 1991, Annapolis, MD, Fallow JB (1996). America s Cup Sail Design. Journal of Wind Engineering and Industrial Aerodynamics 63: Flay RGJ (1996). A Twisted Flow Wind Tunnel for Testing Yacht Sails. Journal of Wind Engineering and Industrial Aerodynamics 63: Fluent (2006). Fluent 6.3 User s Guide. Fluent, Inc., Lebanon, NH. 7. Gerritsma J, Kerwin JE, Moeyes G (1975). Determination of Sail Forces Based on Full Scale Measurements and Model Tests. MIT Dept of Ocean Engineering, Boston, MA Report No Harland J, Myers M (1984). Seamanship in the Age of Sail. Naval Institute Press, Annapolis, MD. 9. Hedges KL, Richards PJ, Mallinson GD (1996). Computer Modelling of Downwind Sails. Journal of Wind Engineering and Industrial Aerodynamics 63: Hochkirch K, Brandt H (1999). Fullscale Hydrodynamics Force Measurement on the Berlin Sailing Dynamometer. Proceedings of the Fourteenth Chesapeake Sailing Yacht Symposium, January 30, 1999, Annapolis, MD, Kerwin JE (1978). A Velocity Prediction Program for Ocean Racing Yachts. Report No , Massachusetts Institute of Technology. 12. Lasher WC, Musho TD, McKee KC, Rybka, W (2007). An Aerodynamic Analysis of the U.S. Brig Niagara. Proceedings of the Eighteenth Chesapeake Sailing Yacht Symposium, March 2 3, 2007, Annapolis, MD, Lasher WC, Richards PJ (2007). Validation of Reynoldsaveraged NavierStokes Simulations for International American s Cup Class Spinnaker Force Coefficients in Atmospheric Boundary Layer. Journal of Ship Research 51(1): Lasher WC, Sonnenmeier JR (2008). A Practical Analysis of RANS Simulations for Spinnaker Aerodynamics. Journal of Wind Engineering and Industrial Aerodynamics 96(2): Lasher WC, Sonnenmeier JR, Forsman DR, Zhang C, White K (2003). Experimental Force Coefficients for a Parametric Series of Spinnakers. Proceedings of the Sixteenth Chesapeake Sailing Yacht Symposium, March 21 22, 2003, Annapolis, MD, Leonard BP (1979). A Stable and Accurate Convective Modelling Procedure Based on Quadratic Upstream Interpolation. Computer Methods in Applied Mechanics and Engineering 19: Masuyama Y, Fukasawa T (1997). Full Scale Measurement of Sail Force Coefficients and the Validation of Numerical Calculation Procedure. Proceedings of the Thirteenth Chesapeake Sailing Yacht Symposium, January 24 25, 1997, Annapolis, MD, Miles JC, Johnson B, Womak J, Franzen I (2007). SNAME s Stabiility Letter Improvement Project (SLIP) for Passenger Sailing Vessels. Proceedings of the Eighteenth Chesapeake Sailing Yacht Symposium, March 2 3, 2007, Annapolis, MD, Milgram JH (1968). The Analytical Design of Yacht Sails. Transactions of the Society of Naval Architects and Marine Engineers. 20. Milgram JH, Peters DB, Eckhouse DN (1993). Modeling IACC Sail Forces by Combining Measurements with CFD. Proceedings of the Eleventh Chesapeake Sailing Yacht Symposium, Annapolis, MD, Patankar SV (1980). Numerical Heat Transfer and Fluid Flow. Hemisphere, Washington, D.C. 22. Ranzenbach R, Mairs C (1997). Experimental Determination of Sail Performance and Blockage Corrections. Proceedings of the Thirteenth Chesapeake Sailing Yacht Symposium, Annapolis, MD, 1997, Ranzenbach R, Mairs C (1999). Wind Tunnel Testing of Offwind Sails. Proceedings of the Fourteenth Chesapeake Sailing Yacht Symposium, Annapolis, MD, 1999, Richards PJ (1997). The Effect of Wind Profile and Twist on Downwind Sail Performance. Journal of Wind Engineering and Industrial Aerodynamics 67&68: Richards PJ, Hoxey RP (1993). Appropriate Boundary Conditions for Computational Wind Engineering Using the kε Turbulence Model. Journal of Wind Engineering and Industrial Aerodynamics 46&47: Roache PJ (1998). Verification of Codes and Calculations. AIAA Journal 36(5): Shih TH, Liou W, Shabbir A, Yang Z, Zhu J (1995). A New kε EddyViscosity Model for High Reynolds Number Turbulent Flows Model Development and Validation. Computers and Fluids 24(3):
Computational Analysis of the S Airfoil Aerodynamic Performance
Computational Analysis of the 2453S Airfoil Aerodynamic Performance Luis VelazquezAraque and Jiří Nožička 2 Department of Mechanical Engineering National University of Táchira, San Cristóbal 5, Venezuela
More informationANALYSIS OF AERODYNAMIC CHARACTERISTICS OF A SUPERCRITICAL AIRFOIL FOR LOW SPEED AIRCRAFT
ANALYSIS OF AERODYNAMIC CHARACTERISTICS OF A SUPERCRITICAL AIRFOIL FOR LOW SPEED AIRCRAFT P.Sethunathan 1, M.Niventhran 2, V.Siva 2, R.Sadhan Kumar 2 1 Asst.Professor, Department of Aeronautical Engineering,
More informationAIRFOIL PROFILE OPTIMIZATION OF AN AIR SUCTION EQUIPMENT WITH AN AIR DUCT
THERMAL SCIENCE, Year 2015, Vol. 19, No. 4, pp. 12171222 1217 AIRFOIL PROFILE OPTIMIZATION OF AN AIR SUCTION EQUIPMENT WITH AN AIR DUCT by Li QIU a,b, Rui WANG a,*, XiaoDong CHEN b, and DePeng WANG
More informationDatabase of Sail Shapes vs. Sail Performance and Validation of Numerical Calculation for Upwind Condition
THE 18 th CHESAPEAKE SAILING YACHT SYMPOSIUM ANNAPOLIS, MARYLAND, MARCH 2007 Database of Sail Shapes vs. Sail Performance and Validation of Numerical Calculation for Upwind Condition Yutaka Masuyama, Kanazawa
More information1. A tendency to roll or heel when turning (a known and typically constant disturbance) 2. Motion induced by surface waves of certain frequencies.
Department of Mechanical Engineering Massachusetts Institute of Technology 2.14 Analysis and Design of Feedback Control Systems Fall 2004 October 21, 2004 Case Study on Ship Roll Control Problem Statement:
More informationCentre for Offshore Renewable Energy Engineering, School of Energy, Environment and Agrifood, Cranfield University, Cranfield, MK43 0AL, UK 2
Fluid Structure Interaction Modelling of A Novel 10MW VerticalAxis Wind Turbine Rotor Based on Computational Fluid Dynamics and Finite Element Analysis Lin Wang 1*, Athanasios Kolios 1, PierreLuc Delafin
More informationA new methodology to evaluate the flying shape of a sail
International conference on Graphics Engineering June 19 th 21 st, 2013 Madrid,Spain INGEGRAF ADM AIP PRIMECA A new methodology to evaluate the flying shape of a sail T. Ingrassia (a), V. Nigrelli (a),
More informationWind tunnel effects on wingtip vortices
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition 47 January 2010, Orlando, Florida AIAA 2010325 Wind tunnel effects on wingtip vortices Xin Huang 1, Hirofumi
More informationComputational Investigation of Airfoils with Miniature Trailing Edge Control Surfaces
AIAA245 Computational Investigation of Airfoils with Miniature Trailing Edge Control Surfaces HakTae Lee, Ilan M. Kroo Stanford University, Stanford, CA 9435 Abstract Miniature trailing edge effectors
More informationSTABILITY OF MULTIHULLS Author: Jean Sans
STABILITY OF MULTIHULLS Author: Jean Sans (Translation of a paper dated 10/05/2006 by Simon Forbes) Introduction: The capsize of Multihulls requires a more exhaustive analysis than monohulls, even those
More informationS0300A6MAN010 CHAPTER 2 STABILITY
CHAPTER 2 STABILITY 21 INTRODUCTION This chapter discusses the stability of intact ships and how basic stability calculations are made. Definitions of the state of equilibrium and the quality of stability
More informationTHE 21 st CHESAPEAKE SAILING YACHT SYMPOSIUM ANNAPOLIS, MARYLAND, MARCH 2013
THE 21 st CHESAPEAKE SAILING YACHT SYMPOSIUM ANNAPOLIS, MARYLAND, MARCH 2013 A wind tunnel study of the interaction between two sailing yachts P.J. Richards, D.J. Le Pelley, D. Jowett, J. Little, O. Detlefsen
More informationNUMERICAL SIMULATION OF THE TURBULENT FLOW PAST UPWIND YACHT SAILS
NUMERICAL SIMULATION OF THE TURBULENT FLOW PAST UPWIND YACHT SAILS S. J. COLLIE, M. G. GERRITSEN, AND M. J. O SULLIVAN Abstract. The ßow past upwind yacht sails is simulated using the commercial computational
More informationNumerical study of aeroelasticity of sails
Numerical study of aeroelasticity of sails Alessandro Leone, Luciano Teresi* SMFM@DiS, Università Roma Tre, Roma, Italy * Corresponding author: Dip. of Studies on Structures, via Corrado Segre 6, 00146
More informationOPTIMIZATION OF SINGLE STAGE AXIAL FLOW COMPRESSOR FOR DIFFERENT ROTATIONAL SPEED USING CFD
http:// OPTIMIZATION OF SINGLE STAGE AXIAL FLOW COMPRESSOR FOR DIFFERENT ROTATIONAL SPEED USING CFD Anand Kumar S malipatil 1, Anantharaja M.H 2 1,2 Department of Thermal Power Engineering, VTURO Gulbarga,
More informationWind tunnel test and numerical simulation of wind pressure on a highrise building
Journal of Chongqing University (English Edition) [ISSN 16718224] Vol. 9 No. 1 March 2010 Article ID: 16718224(2010)01004707 To cite this article: AL ZOUBI Feras, LI Zhengliang, WEI Qike, SUN Yi.
More informationFORCE AND PRESSURE INVESTIGATION OF MODERN ASYMMETRIC SPINNAKERS
FORCE AND PRESSURE INVESTIGATION OF MODERN ASYMMETRIC SPINNAKERS I M Viola, R G J Flay, Yacht Research Unit, Mechanical Engineering Department, The University of Auckland, New Zealand SUMMARY An innovative
More informationTHE 22nd CHESAPEAKE SAILING YACHT SYMPOSIUM ANNAPOLIS, MARYLAND, MARCH 2016
THE 22nd CHESAPEAKE SAILING YACHT SYMPOSIUM ANNAPOLIS, MARYLAND, MARCH 2016 A Comparison of RANS and LES for Upwind Sailing Aerodynamics Stefano Nava, University of Auckland, Auckland, New Zealand John
More informationINCLINOMETER DEVICE FOR SHIP STABILITY EVALUATION
Proceedings of COBEM 2009 Copyright 2009 by ABCM 20th International Congress of Mechanical Engineering November 1520, 2009, Gramado, RS, Brazil INCLINOMETER DEVICE FOR SHIP STABILITY EVALUATION Helena
More informationNumerical Study of GiromillType Wind Turbines with Symmetrical and Nonsymmetrical Airfoils
European International Journal of Science and Technology Vol. 2 No. 8 October 2013 Numerical Study of GiromillType Wind Turbines with Symmetrical and Nonsymmetrical Airfoils Prathamesh Deshpande and
More informationInvestigations on the Aerodynamic Forces of 2D Square Lattice Tower Section Using CFD
J. Energy Power Sources Vol. 1, No. 5, 2014, pp. 270277 Received: August 14, 2014, Published: November 30, 2014 Journal of Energy and Power Sources www.ethanpublishing.com Investigations on the Aerodynamic
More informationWind resource assessment over a complex terrain covered by forest using CFD simulations of neutral atmospheric boundary layer with OpenFOAM
Wind resource assessment over a complex terrain covered by forest using CFD simulations of neutral atmospheric boundary layer with OpenFOAM Nikolaos Stergiannis nstergiannis.com nikolaos.stergiannis@vub.ac.be
More informationEFFECT OF CORNER CUTOFFS ON FLOW CHARACTERISTICS AROUND A SQUARE CYLINDER
EFFECT OF CORNER CUTOFFS ON FLOW CHARACTERISTICS AROUND A SQUARE CYLINDER Yoichi Yamagishi 1, Shigeo Kimura 1, Makoto Oki 2 and Chisa Hatayama 3 ABSTRACT It is known that for a square cylinder subjected
More informationAbstract. 1 Introduction
Developments in modelling ship rudderpropeller interaction A.F. Molland & S.R. Turnock Department of Ship Science, University of Southampton, Highfield, Southampton, S017 IBJ, Hampshire, UK Abstract A
More informationWIND LOADS / MOORING & FISH TAILING. Arjen Koop, Senior Project Manager Offshore Rogier Eggers, Project Manager Ships
WIND LOADS / MOORING & FISH TAILING Arjen Koop, Senior Project Manager Offshore Rogier Eggers, Project Manager Ships OVERVIEW Wind Loads Wind shielding Fish tailing? 2 WIND LOADS FOR OFFSHORE MARS TLP
More informationAIR FLOW DISTORTION OVER MERCHANT SHIPS.
AIR FLOW DISTORTION OVER MERCHANT SHIPS. M. J. Yelland, B. I. Moat and P. K. Taylor April 2001 Extended Abstract Anemometers on voluntary observing ships (VOS) are usually sited above the bridge in a region
More informationLift for a Finite Wing. all real wings are finite in span (airfoils are considered as infinite in the span)
Lift for a Finite Wing all real wings are finite in span (airfoils are considered as infinite in the span) The lift coefficient differs from that of an airfoil because there are strong vortices produced
More informationDetermination of the wind pressure distribution on the facade of the triangularly shaped highrise building structure
Determination of the wind pressure distribution on the facade of the triangularly shaped highrise building structure Norbert Jendzelovsky 1,*, Roland Antal 1 and Lenka Konecna 1 1 STU in Bratislava, Faculty
More informationEXPERIMENTAL AND NUMERICAL STUDY OF A TWO ELEMENT WING WITH GURNEY FLAP
25 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES EXPERIMENTAL AND NUMERICAL STUDY OF A TWO ELEMENT WING WITH GURNEY FLAP F.M. Catalano PhD.( catalano@sc.usp.br ) *, G. L. Brand * * Aerodynamic
More informationA COMPARATIVE STUDY OF MIX FLOW PUMP IMPELLER CFD ANALYSIS AND EXPERIMENTAL DATA OF SUBMERSIBLE PUMP
IMPACT: International Journal of Research in Engineering & Technology (IMPACT: IJRET) ISSN 23218843 Vol. 1, Issue 3, Aug 2013, 5764 Impact Journals A COMPARATIVE STUDY OF MIX FLOW PUMP IMPELLER CFD ANALYSIS
More informationDevelopment of TEU Type Mega Container Carrier
Development of 8 700 TEU Type Mega Container Carrier SAKAGUCHI Katsunori : P. E. Jp, Manager, Ship & Offshore Basic Design Department, IHI Marine United Inc. TOYODA Masanobu : P. E, Jp, Ship & Offshore
More informationStatic Extended Trailing Edge for Lift Enhancement: Experimental and Computational Studies
Static Extended Trailing Edge for Lift Enhancement: Experimental and Computational Studies T. Liu, J. Montefort, W. Liou Western Michigan University Kalamazoo, MI 49008 and Q. Shams NASA Langley Research
More informationImproving pedestrian wind comfort around a high rise building in an urban area: a case study by CFD simulations and onsite measurements
Improving pedestrian wind comfort around a high rise building in an urban area: a case study by CFD simulations and onsite measurements WD Janssen, B Blocken, and T van Hooff Unit Building Physics and
More informationEXPERIMENTAL STUDY OF WIND PRESSURES ON IRREGULAR PLAN SHAPE BUILDINGS
BBAA VI International Colloquium on: Bluff Bodies Aerodynamics & Applications Milano, Italy, July, 224 8 EXPERIMENTAL STUDY OF WIND PRESSURES ON IRREGULAR PLAN SHAPE BUILDINGS J. A. Amin and A. K. Ahuja
More informationWINDINDUCED LOADS OVER DOUBLE CANTILEVER BRIDGES UNDER CONSTRUCTION
WINDINDUCED LOADS OVER DOUBLE CANTILEVER BRIDGES UNDER CONSTRUCTION S. Pindado, J. Meseguer, J. M. Perales, A. SanzAndres and A. Martinez Key words: Wind loads, bridge construction, yawing moment. Abstract.
More informationAir Craft Winglet Design and Performance: Cant Angle Effect
Journal of Robotics and Mechanical Engineering Research Air Craft Winglet Design and Performance: Cant Angle Effect Eslam Said Abdelghany 1, Essam E Khalil 2*, Osama E Abdellatif 3 and Gamal elhariry 4
More informationa highrise (78 m) building in the city of Antwerp. CFD simulations are performed for the building with and without secondskin facade concept impleme
The Seventh International Colloquium on Bluff Body Aerodynamics and Applications (BBAA7) Shanghai, China; September 26, 2012 CFD analysis of wind comfort on highrise building balconies: validation and
More informationTHEORY OF WINGS AND WIND TUNNEL TESTING OF A NACA 2415 AIRFOIL. By Mehrdad Ghods
THEORY OF WINGS AND WIND TUNNEL TESTING OF A NACA 2415 AIRFOIL By Mehrdad Ghods Technical Communication for Engineers The University of British Columbia July 23, 2001 ABSTRACT Theory of Wings and Wind
More informationNumerical simulation and analysis of aerodynamic drag on a subsonic train in evacuated tube transportation
Journal of Modern Transportation Volume 20, Number 1, March 2012, Page 4448 Journal homepage: jmt.swjtu.edu.cn DOI: 10.1007/BF03325776 1 Numerical simulation and analysis of aerodynamic drag on a subsonic
More informationBlade Design and Performance Analysis of Wind Turbine
International Journal of ChemTech Research CODEN( USA): IJCRGG ISSN : 09744290 Vol.5, No.2, pp 10541061, AprilJune 2013 ICGSEE2013[14 th 16 th March 2013] International Conference on Global Scenario
More informationTHE BRIDGE COLLAPSED IN NOVEMBER 1940 AFTER 4 MONTHS OF ITS OPENING TO TRAFFIC!
OUTLINE TACOMA NARROWS BRIDGE FLOW REGIME PAST A CYLINDER VORTEX SHEDDING MODES OF VORTEX SHEDDING PARALLEL & OBLIQUE FLOW PAST A SPHERE AND A CUBE SUMMARY TACOMA NARROWS BRIDGE, USA THE BRIDGE COLLAPSED
More informationNUMERICAL INVESTIGATION OF FLOW OVER MULTI ELEMENT AIRFOILS WITH LIFTENHANCING TABS
28 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES NUMERICAL INVESTIGATION OF FLOW OVER MULTI ELEMENT AIRFOILS WITH LIFTENHANCING TABS Zhenhui Zhang, Dong Li National Key Laboratory of Science
More informationEdit this text for your title
Edit this text for your title MEK 4450 Marine Operations Edit this text for your subtitle Presenter name, location, date etc. Kværner ASA / DNV, Fall 2013 Lesson 2/3 Lift phases Load out Transportation
More informationON THE EFFECT OF LIFT FORCES IN BUBBLE PLUMES
Ninth International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 1012 December 2012 ON THE EFFECT OF LIFT FORCES IN BUBBLE PLUMES Jan Erik OLSEN* and Mihaela POPESCU
More informationAerodynamic Modification of CFR Formula SAE Race Car
Aerodynamic Modification of CFR Formula SAE Race Car Brandon M. Verhun, Trevor D. Haight, and Thomas A. Mahank Department of Mechanical Engineering Saginaw Valley State University University Center, MI
More informationThe OTSS System for Drift and Response Prediction of Damaged Ships
The OTSS System for Drift and Response Prediction of Damaged Ships Shoichi Hara 1, Kunihiro Hoshino 1,Kazuhiro Yukawa 1, Jun Hasegawa 1 Katsuji Tanizawa 1, Michio Ueno 1, Kenji Yamakawa 1 1 National Maritime
More informationLargeeddy simulation of a turbulent buoyant helium plume
Center for Turbulence Research Annual Research Briefs 8 45 Largeeddy simulation of a turbulent buoyant helium plume By G. Blanquart AND H. Pitsch. Motivation and objectives The numerical simulation of
More informationInvestigation on Divergent Exit Curvature Effect on Nozzle Pressure Ratio of Supersonic Convergent Divergent Nozzle
RESEARCH ARTICLE OPEN ACCESS Investigation on Divergent Exit Curvature Effect on Nozzle Pressure Ratio of Supersonic Convergent Divergent Nozzle Shyamshankar.M.B*, Sankar.V** *(Department of Aeronautical
More informationOPTIMIZING THE LENGTH OF AIR SUPPLY DUCT IN CROSS CONNECTIONS OF GOTTHARD BASE TUNNEL. Rehan Yousaf 1, Oliver Scherer 1
OPTIMIZING THE LENGTH OF AIR SUPPLY DUCT IN CROSS CONNECTIONS OF GOTTHARD BASE TUNNEL Rehan Yousaf 1, Oliver Scherer 1 1 Pöyry Infra Ltd, Zürich, Switzerland ABSTRACT Gotthard Base Tunnel with its 57 km
More informationEFFECT OF GURNEY FLAPS AND WINGLETS ON THE PERFORMANCE OF THE HAWT
Chapter6 EFFECT OF GURNEY FLAPS AND WINGLETS ON THE PERFORMANCE OF THE HAWT 6.1 Introduction The gurney flap (wicker bill) was a small flat tab projecting from the trailing edge of a wing. Typically it
More informationAn Experimental Investigation of Asymmetric Spinnaker Aerodynamics Using Pressure and Sail Shape Measurements
An Experimental Investigation of Asymmetric Spinnaker Aerodynamics Using Pressure and Sail Shape Measurements D. Motta, Yacht Research Unit, University of Auckland, New Zealand, dmot267@aucklanduni.ac.nz
More informationNUMERICAL SIMULATION OF THE FLOW AROUND A CCLASS CATAMARAN Filipe Santos Carvalho
NUMERICAL SIMULATION OF THE FLOW AROUND A CCLASS CATAMARAN Filipe Santos Carvalho filipe.s.carvalho@tecnico.ulisboa.pt Instituto Superior Técnico, Universidade de Lisboa, Portugal December 2016 ABSTRACT
More informationNumerical and Experimental Investigations of Lift and Drag Performances of NACA 0015 Wind Turbine Airfoil
International Journal of Materials, Mechanics and Manufacturing, Vol. 3, No., February 2 Numerical and Experimental Investigations of Lift and Drag Performances of NACA Wind Turbine Airfoil İzzet Şahin
More informationHigh Sweptback Delta Wing Flow
Advanced Materials Research Submitted: 20140625 ISSN: 16628985, Vol. 1016, pp 377382 Accepted: 20140625 doi:10.4028/www.scientific.net/amr.1016.377 Online: 20140828 2014 Trans Tech Publications,
More informationDevelopment of Biomimicry Wind Louver Surface Design
International Proceedings of Chemical, Biological and Environmental Engineering, V0l. 93 (2016) DOI: 10.7763/IPCBEE. 2016. V93. 6 Development of Biomimicry Wind Louver Surface Design Jaepil Choi 1, Donghwa
More informationPredicting the Suction Gas Superheating in Reciprocating Compressors
Purdue University Purdue epubs International Compressor Engineering Conference School of Mechanical Engineering 2014 Predicting the Suction Gas Superheating in Reciprocating Compressors Jonatas Ferreira
More informationAgood tennis player knows instinctively how hard to hit a ball and at what angle to get the ball over the. Ball Trajectories
42 Ball Trajectories Factors Influencing the Flight of the Ball Nathalie Tauziat, France By Rod Cross Introduction Agood tennis player knows instinctively how hard to hit a ball and at what angle to get
More informationAdjustment of the kω SST turbulence model for prediction of airfoil characteristics near stall
Journal of Physics: Conference Series PAPER OPEN ACCESS Adjustment of the ω SST turbulence model for prediction of airfoil characteristics near stall To cite this article: A A Matyusheno and A V Garbaru
More informationIMAGEBASED STUDY OF BREAKING AND BROKEN WAVE CHARACTERISTICS IN FRONT OF THE SEAWALL
IMAGEBASED STUDY OF BREAKING AND BROKEN WAVE CHARACTERISTICS IN FRONT OF THE SEAWALL Weijie Liu 1 and Yoshimitsu Tajima 1 This study aims to study the breaking and broken wave characteristics in front
More informationAPPLICATION OF COMPUTATIONAL FLUID DYNAMICS (CFD) IN WIND ANALYSIS OF TALL BUILDINGS. Damith Mohotti, Priyan Mendis, Tuan Ngo
APPLICATION OF COMPUTATIONAL FLUID DYNAMICS (CFD) IN WIND ANALYSIS OF TALL BUILDINGS Damith Mohotti, Priyan Mendis, Tuan Ngo Department of Infrastructures Engineering, The University of Melbourne, Victoria,
More informationNumerical simulation of topographic effects on wind flow fields over complex terrain
The Eighth AsiaPacific Conference on Wind Engineering, December 10 14, 2013, Chennai, India Numerical simulation of topographic effects on wind flow fields over complex terrain B.W. Yan 1, Q.S. LI 2,
More informationImproved Aerodynamic Characteristics of Aerofoil Shaped Fuselage than that of the Conventional Cylindrical Shaped Fuselage
International Journal of Scientific & Engineering Research Volume 4, Issue 1, January213 1 Improved Aerodynamic Characteristics of Aerofoil Shaped Fuselage than that of the Conventional Cylindrical Shaped
More informationValidating CFD Simulation Results: Wind flow around a surface mounted cube in a turbulent channel flow
Validating CFD Simulation Results: Wind flow around a surface mounted cube in a turbulent channel flow ISLAM ABOHELA 1, NEVEEN HAMZA 1, STEVEN DUDEK 1 1 School of Architecture, Planning and Landscape,
More informationInternational Journal of Technical Research and Applications eissn: , Volume 4, Issue 3 (MayJune, 2016), PP.
DESIGN AND ANALYSIS OF FEED CHECK VALVE AS CONTROL VALVE USING CFD SOFTWARE R.Nikhil M.Tech Student Industrial & Production Engineering National Institute of Engineering Mysuru, Karnataka, India 570008
More informationWind loading on cladding and glazed façades
Wind loading on cladding and glazed façades Dr Mauro OVEREND Lecturer in Architectural Structures School of the Built Environment, University of Nottingham UK Kenneth ZAMMIT PhD Student School of the Built
More informationKinematics of Vorticity
Kinematics of Vorticity Vorticity Ω Ω= V 2 circumferentially averaged angular velocity of the fluid particles Sum of rotation rates of perpendicular fluid lines Nonzero vorticity doesn t imply spin.ω=0.
More informationCASE STUDY FOR USE WITH SECTION B
GCE A level 135/01B PHYSICS ASSESSMENT UNIT PH5 A.M. THURSDAY, 0 June 013 CASE STUDY FOR USE WITH SECTION B Examination copy To be given out at the start of the examination. The prerelease copy must
More informationAvailable online at ScienceDirect. Procedia Engineering 126 (2015 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 126 (2015 ) 431 435 7th International Conference on Fluid Mechanics, ICFM7 Selfpropulsion of a flexible plunging foil near
More informationOFFSHORE RACING CONGRESS
World Leader in Rating Technology OFFSHORE RACING CONGRESS ORC Speed Guide Explanation 1. INTRODUCTION The ORC Speed Guide is a customcalculated manual for improving performance for an individual boat.
More informationScales of Atmospheric Motion Scale Length Scale (m) Time Scale (sec) Systems/Importance Molecular (neglected)
Supplement Wind, Fetch and Waves Scales of Atmospheric Motion Scale Length Scale (m) Time Scale (sec) Systems/Importance Molecular 107  102 101 (neglected) Coriolis not important Turbulent 102 10
More informationComparative Stability Analysis of a Frigate According to the Different Navy Rules in Waves
Comparative Stability Analysis of a Frigate According to the Different Navy Rules in Waves ABSTRACT Emre Kahramano lu, Technical University, emrek@yildiz.edu.tr Hüseyin Y lmaz,, hyilmaz@yildiz.edu.tr Burak
More informationSAILS AERODYNAMICS. A seatrip around sail shapes. V.G. Chapin, Marine 2017, Nantes
SAILS AERODYNAMICS A seatrip around sail shapes V.G. Chapin, Marine 2017, Nantes «Ceux qui pensent que c est impossible sont priés de ne pas déranger ceux qui essaient.» OUTLINE Sail: Curry (1925) Sails
More informationAn Investigation into the Capsizing Accident of a Pusher Tug Boat
An Investigation into the Capsizing Accident of a Pusher Tug Boat Harukuni Taguchi, National Maritime Research Institute (NMRI) taguchi@nmri.go.jp Tomihiro Haraguchi, National Maritime Research Institute
More informationPerformance of GOE387 Airfoil Using CFD
International Journal of Aerospace Sciences 2017, 5(1): 17 DOI: 10.5923/j.aerospace.20170501.01 Performance of GOE387 Airfoil Using CFD Mohamed A. Fouad Kandil *, Abdelrady Okasha Elnady Mechatronics
More informationWind Flow Validation Summary
IBHS Research Center Validation of Wind Capabilities The Insurance Institute for Business & Home Safety (IBHS) Research Center fullscale test facility provides opportunities to simulate natural wind conditions
More informationA Numerical Study of Thickness Effect of the Symmetric NACA 4Digit Airfoils on Self Starting Capability of a 1kW HType Vertical Axis Wind Turbine
International Journal of Mechanical Engineering and Applications 2015; 3(31): 716 Published online February 13, 2015 (http://www.sciencepublishinggroup.com/j/ijmea) doi: 10.11648/j.ijmea.s.2015030301.12
More informationConventional Ship Testing
Conventional Ship Testing Experimental Methods in Marine Hydrodynamics Lecture in week 34 Chapter 6 in the lecture notes 1 Conventional Ship Testing  Topics: Resistance tests Propeller open water tests
More informationWind loads investigations of HAWT with wind tunnel tests and site measurements
loads investigations of HAWT with wind tunnel tests and site measurements Shigeto HIRAI, Senior Researcher, Nagasaki R&D Center, Technical Headquarters, MITSUBISHI HEAVY INDSUTRIES, LTD, Fukahori, Nagasaki,
More informationNational Maritime Center
National Maritime Center Providing Credentials to Mariners (Sample Examination) Page 1 of 8 Choose the best answer to the following Multiple Choice Questions. 1. In illustration D001SL, what is the edge
More informationAnalysis of Wind Turbine Blade
International Journal of Current Engineering and Technology EISSN 2277 4106, PISSN 23475161 2014 INPRESSCO, All Rights Reserved Available at http://inpressco.com/category/ijcet Research Article Analysis
More informationAvailable online at Procedia Engineering 200 (2010) (2009) In situ drag measurements of sports balls
Available online at www.sciencedirect.com Procedia Engineering 200 (2010) (2009) 2437 2442 000 000 Procedia Engineering www.elsevier.com/locate/procedia 8 th Conference of the International Sports Engineering
More informationOPTIMIZATION OF A WAVE CANCELLATION MULTIHULL SHIP USING CFD TOOLS
OPTIMIZATION OF A AVE CANCELLATION MULTIHULL SHIP USING CFD TOOLS C. Yang, R. Löhner and O. Soto School of Computational Sciences, George Mason University Fairfax VA 0304444, USA ABSTRACT A simple CFD
More informationIMPROVING THE DESIGN OF SAILS USING CFD AND OPTIMIZATION ALGORITHMS
High Performance Yacht Design Conference Auckland, 46 December, IMPROVING THE DESIGN OF SAILS USING CFD AND OPTIMIZATION ALGORITHMS Sriram Shankaran 1, ssriram@stanford.edu Tyler Doyle, tyler@stanford.edu
More informationNUMERICAL SIMULATION OF CROSSFLOW AROUND FOUR CIRCULAR CYLINDERS INLINE SQUARE CONFIGURATION NEAR A PLANE WALL
NUMERICAL SIMULATION OF CROSSFLOW AROUND FOUR CIRCULAR CYLINDERS INLINE SQUARE CONFIGURATION NEAR A PLANE WALL A. Grummy Wailanduw 1, Triyogi Yuwono 2 and Wawan Aries Widodo 2 1 Student of Doctoral Program
More informationNUMERICAL AND EXPERIMENTAL INVESTIGATION OF THE FLOWFIELD IN A BLOWDOWN WIND TUNNEL
ICAS 2002 CONGRESS NUERICAL AND EXPERIENTAL INVESTIGATION OF THE FLOWFIELD IN A BLOWDOWN WIND TUNNEL Ali Rebaine, ahmood Khalid, Cabot Broughton and Fred Ellis Institute for Aerospace Research, National
More informationPower efficiency and aerodynamic forces measurements on the Dettwilerwind turbine
Institut für Fluiddynamik ETH Zentrum, ML H 33 CH8092 Zürich P rof. Dr. Thomas Rösgen Sonneggstrasse 3 Telefon +4144632 2646 Fax +4144632 1147 roesgen@ifd.mavt.ethz.ch www.ifd.mavt.ethz.ch Power efficiency
More informationIMO REVISION OF THE INTACT STABILITY CODE. Proposal of methodology of direct assessment for stability under dead ship condition. Submitted by Japan
INTERNATIONAL MARITIME ORGANIZATION E IMO SUBCOMMITTEE ON STABILITY AND LOAD LINES AND ON FISHING VESSELS SAFETY 49th session Agenda item 5 SLF 49/5/5 19 May 2006 Original: ENGLISH REVISION OF THE INTACT
More informationOptimizing the shape and speed performance of the boats
Optimizing the shape and speed performance of the boats A report submitted to Modelling Week and Study Group Meeting on Industrial Problems2013 by AWL Pubudu Thilan 1, Baikuntha Nath Sahu 2, Darshan Lal
More informationCOURSE OBJECTIVES CHAPTER 9
COURSE OBJECTIVES CHAPTER 9 9. SHIP MANEUVERABILITY 1. Be qualitatively familiar with the 3 broad requirements for ship maneuverability: a. Controls fixed straightline stability b. Response c. Slow speed
More informationComputational evaluation of pedestrian wind comfort and wind safety around a highrise building in an urban area
International Environmental Modelling and Software Society (iemss) 7th Intl. Congress on Env. Modelling and Software, San Diego, CA, USA, Daniel P. Ames, Nigel W.T. Quinn and Andrea E. Rizzoli (Eds.) http://www.iemss.org/society/index.php/iemss2014proceedings
More informationJOURNAL PUBLICATIONS
1 JOURNAL PUBLICATIONS 71. Lee, T., Mageed, A., Siddiqui, B. and Ko, L.S., (2016) Impact of ground proximity on aerodynamic properties of an unsteady NACA 0012 airfoil, submitted to Journal of Aerospace
More informationReview of Equivalent Neutral Winds and Stress
Review of Equivalent Neutral Winds and Stress Mark A. Bourassa Center for OceanAtmospheric Prediction Studies, Geophysical Fluid Dynamics Institute & Department of Earth, Ocean and Atmospheric Science
More informationAre Advanced Wind Flow Models More Accurate? A Test of Four Models
Are Advanced Wind Flow Models More Accurate? A Test of Four Models Philippe Beaucage, PhD Senior Research Scientist Michael C. Brower, PhD Chief Technical Officer Brazil Wind Power Conference 2012 Albany
More informationClarification of Behavior of Huge Tsunami Action on Bridges  Hydraulic Model Experiment and Simulation Technology 
Clarification of Behavior of Huge Tsunami Action on Bridges  Hydraulic Model Experiment and Simulation Technology  21 TOSHIMITSU SUZUKI *1 RIKUMA SHIJO *2 KAORU YOKOYAMA *3 SYUNICHI IKESUE *4 HIROFUMI
More informationOptimization of a WingSail shape for a small boat
STAR Global Conference 2014 Vienna, March 1719 Optimization of a WingSail shape for a small boat G. Lombardi, F. Cartoni, M. Maganzi Dept. of Civil and Industrial Engineering of Pisa Aeronautical Section
More informationCFD SIMULATION STUDY OF AIR FLOW AROUND THE AIRFOIL USING THE MAGNUS EFFECT
Magnus effect, simulation, air flow Patryk SOKOŁOWSKI *, Jacek CZARNIGOWSKI **, Paweł MAGRYTA *** CFD SIMULATION STUDY OF AIR FLOW AROUND THE AIRFOIL USING THE MAGNUS EFFECT Abstract The article presents
More informationNew Technologies applied to the design and optimization of tunnel ventilation systems
New Technologies applied to the design and optimization of tunnel ventilation systems Ing. Justo Suárez Area Manager  ZITRON Dot. Ing. Massimiliano Bringiotti Managing Director  GEOTUNNEL Ing. Ana Belén
More informationAerodynamic investigation of Winglets on Wind Turbine Blades using CFD
RisøR1543(EN) Aerodynamic investigation of Winglets on Wind Turbine Blades using CFD Jeppe Johansen and Niels N. Sørensen Risø National Laboratory Roskilde Denmark February 26 Author: Jeppe Johansen
More informationWrite important assumptions used in derivation of Bernoulli s equation. Apart from an airplane wing, give an example based on Bernoulli s principle
HW#3 Sum07 #1. Answer in 4 to 5 lines in the space provided for each question: (a) A tank partially filled with water has a balloon well below the free surface and anchored to the bottom by a string. The
More informationThe Effect of Von Karman Vortex Street on Building Ventilation
The Effect of Von Karman Vortex Street on Building Ventilation P.Praveen Kumar Abstract This paper deals with the utilisation of the von Karman vortex street principle to maximise air flow into buildings.
More information