Computational and experimental study on performance of sails of a yacht

It is important to understand flow characteristics and performances of sails for both sailors and designers who want to have efficient thrust of yacht. In this paper the viscous flows around sail-like rigid wings, which are similar to main and jib sails of a 30 feet sloop, are calculated using a CFD tool. Lift, drag and thrust forces are estimated for various conditions of gap distance between the two sails and the center of effort of the sail system are obtained. Wind tunnel experiments are also carried out to measure aerodynamic forces acting on the sail system and to validate the computation. It is found that the combination of two sails produces the lift force larger than the sum of that produced separately by each sail and the gap distance between the two sails is an important factor to determine total lift and thrust.     

Analysis of Influencing Factors on Lift Coefficients of Autonomous Sailboat Double Sail Propulsion System Based on Vortex Panel Method

Sail is the core part of autonomous sailboat and wing sail is a new type of sail. Wing sail generates not only propulsion but also lateral force and heeling moment. The latter two will affect the navigation status and bring resistance. Double sail can effectively reduce the center of wind pressure and heeling moment. In order to study the effect of distance between two sails, airfoil and attack angle on the total lift coefficient of double sail propulsion system, pressure coefficient distribution and lift coefficient calculation model have been established based on vortex panel method. By using the basic finite solution, the fluid dynamic forces on the two-dimensional sails are computed.The results show that, the distance in the range of 0 to 1 time chord length, when using the same airfoil in the fore and aft sail, the total lift coefficient of the double sail increases with the increase of distance, finally reaches a stable value in the range of one to three times chord length. Lift coefficients of thicker airfoils are more sensitive to the change of distance. The thicker the airfoil, the longer distance is required of the total lift coefficient toward stable.When different airfoils are adopted in fore and aft sail, the total lift coefficient increases with the increase of the thickness of aft sail. The smaller the thickness difference is, the more sensitive to the distance change the lift coefficient is. The thinner the fore sail is, the lower the influence will be on the lift coefficient of aft sail.     

Sloshing in a three-dimensional rectangular tank: Numerical simulation and experimental validation

Pressure variations and three-dimensional effects on liquid sloshing loads in a moving partially filled rectangular tank have been carried out numerically and experimentally. A numerical algorithm based on the volume of fluid (VOF) technique is used to study the non-linear behavior and damping characteristics of liquid sloshing. A moving coordinate system is used to include the non-linearity and avoid the complex boundary conditions of moving walls. The numerical model solves the complete Navier–Stokes equations in primitive variables by using of the finite difference approximations. In order to mitigate a series of discrete impacts, the signal computed is averaged over several time steps. In order to assess the accuracy of the method used, computations are compared with the experimental results. Several configurations of both baffled and unbaffled tanks are studied. Comparisons show good agreement for both impact and non- impact type slosh loads in the cases investigated.     

An experimental method for determining the frequency-dependent added mass and added mass moment of inertia for a floating body in heave and pitch motions

Most of the existing relevant materials have been obtained from experiments, in which evaluating the added mass at the resonant frequency corresponding to the peak of a frequency-response curve obtained from the “forced” vibration analysis is the most popular technique. In this paper, a simple experimental method was presented where the “free” vibration responses instead of the “forced” ones were used to determine the values of m_ah and I_ap. The main part of the experimental system is composed of a floating body (model) and a spring–shaft shaker. The “free” vibration of this main part was induced by imposing on it an initial displacement (and/or an initial velocity), and from the time histories of displacements information such as the “damped” natural frequencies, damping ratios, sectional added mass coefficients (C_V and C_P) were obtained. Since the displacements of the spring–shaft shaker are “translational” and those of the floating body due to pitch motions are “angular”, a technique for the transformation between the associated parameters of the two components of the main part was presented.     

On the identification of ship coupled heave–pitch motions using neural networks

This paper outlines a procedure for the derivation of the differential equations describing the free response of a heaving and pitching ship from its stationary response to random waves. The coupled heave–pitch motion of a ship in random seas is modelled as a multi-dimensional Markov process. The partial differential equation describing the transition probability density function, known as the Fokker-Planck equation, for this process is derived. The Fokker-Planck equation is used to derive the random decrement equations for the coupled heave–pitch motion. The parameters in these equations are then identified using a neural network approach. The method is validated using numerical simulations and experimental results. The experimental data was obtained using an icebreaker ship model heaving and pitching in random waves. It is shown that the method produces good results when the system is lightly damped. An extension for using this method to identify couple heave–pitch motion in realistic seas is suggested.     

Laboratory studies of interleaving

This is a review of laboratory experiments on double-diffusive interleaving. Several configurations are discussed and compared, including thermohaline fronts, heated sidewalls, melting ice-blocks, point heat sources, and double-diffusive plumes, using both the heat-salt and sugar-salt systems to generate property anomalies. Two parallels emerge. The vertical scale of intrusions in most configurations is proportional to the ‘natural’ scale given by the property anomaly divided by the density gradient. The speed of advance of intrusions into undisturbed water is proportional to the product of the buoyancy frequency and intrusion thickness, with a constant of proportionality of order 0.005, but which depends on the vertical gradient density ratio and on the experimental configuration.We end with a long list of questions. Further progress will require new or more complex experimental configurations, more quantitative observations, and close comparisons with intrusion theories, both linear and finite-amplitude.     

奥运青岛帆船赛场气象服务关键技术分析

帆船比赛是与气象关系最密切的比赛项目。以青岛8，9月份的天气为分析、预报的对象，在深入分析及参照悉尼2000年奥运气象服务的预报技术组成、实施方法与技术要求的基础上，分析讨论了青岛奥运帆船比赛气象服务体系中关键性环节的技术要求与技术构成。主要研究问题包括：（1）GIS青岛气象资料预报服务系统；（2）青岛地区站点的布设及资料观测系统；（3）青岛夏秋Nowcasting集成预报系统。     

Physical modeling of untrenched submarine pipeline instability

Wave-induced instability of untrenched pipeline on sandy seabed is a ‘wave–soil–pipeline’ coupling dynamic problem. To explore the mechanism of the pipeline instability, the hydrodynamic loading with U-shaped oscillatory flow tunnel is adopted, which is quite different from the previous experiment system. Based on dimensional analysis, the critical conditions for pipeline instability are investigated by altering pipeline submerged weight, diameter, soil parameters, etc. Based on the experimental results, different linear relationships between Froude number (Fr) and non-dimensional pipeline weight (G) are obtained for two constraint conditions. Moreover, the effects of loading history on the pipeline stability are also studied. Unlike previous experiments, sand scouring during the process of pipe’s losing stability is detected in the present experiments. In addition, the experiment results are compared with the previous experiments, based on Wake II model for the calculation of wave-induced forces upon pipeline. It shows that the results of two kinds of experiments are comparable, but the present experiments provide better physical insight of the wave–soil–pipeline coupling effects.     

Simulation of wave propagation over a submerged bar using the VOF method with a two-equation k– turbulence modeling

A two-equation k– turbulence model is used in this paper to simulate the propagation of cnoidal waves over a submerged bar, where the free surface is handled by the volume-of-fluid (VOF) method. Using a VOF partial-cell variable and a donor–acceptor method, the model is capable of treating irregular boundaries, including arbitrary bottom topography and internal obstacles, where the no-slip condition is satisfied. The model also allows the viscous sublayer to be modeled by a wall function approximation implemented in the grid nodes that are immediately adjacent to a wall boundary. The numerical model applied to the propagation of cnoidal waves over a submerged bar can produce results that are in general agreement with some laboratory measurements. Some remarks arising from the comparison between the computational and experimental results are presented.     

Residual ultimate strength of steel plates with longitudinal cracks under axial compression–experiments

The main objective of the present study is to examine the residual ultimate strength characteristics of steel plates with cracking damages under axial compressive actions through experimental investigations. The present study is a sequel of the author's previous paper (Paik, J.K., Satish Kumar, Y.V., Lee, J.M., 2005. Ultimate strength of cracked plate elements under axial compression or tension. Thin-Walled Structures, 43, 237–272). In contrast to the previous paper which dealt with transverse cracks located in the direction normal to the axial loading direction, the present paper is concerned with longitudinal cracks which are located in parallel to the axial loading direction. Similar to the previous paper, the orientation/location and size of cracks inside the plates are varied for the present experimental investigations. The details of experimental results are documented. The database and insights developed from the present work will be useful for cracking damage-tolerant design of steel-plated structures and also for condition assessment or health monitoring of aging steel-plated structures.     

Numerical simulation of the hydrodynamic behaviour of submerged plane nets in current

The aim of this paper is to investigate the shape and tension distribution of fishing nets in current. A numerical model is developed, based on lumped mass method to simplify the net. The motion equation is set up for each lumped mass. The Runge–Kutta–Verner fifth-order and sixth-order method is used to solve these simultaneous equations, and then the displacement and tension of each lumped mass are obtained. In order to verify the validity of the numerical method, model tests have been carried out. The results by the numerical simulation agree well with the experimental data.     

Unsteady thermoelastic analysis of a circular disk with a hot central core and subjected to adiabatic conditions

The unsteady thermoelastic analysis of a cooling circular disk or cylinder which is originally at uniform temperature is a classical problem of the theory of thermal stresses. More recent studies consider the case of composite structural configurations. The present paper deals with a situation which, apparently, has not been previously considered: unsteady thermal stresses caused by the presence of a hot, central nucleus. The temperature field is obtained in terms of a Fourier–Bessel expansion and then, radial and tangential stresses are evaluated analytically. The problem is of basic interest in mechanical and naval engineering systems.     

Distribution and abundance of the Pacific sardine (Sardinops sagax) in the Gulf of California and their relation with the environment

In 1989–90 the small pelagic fishery of the Gulf of California began to show a very marked decline in the catch of its main component, the Pacific sardine (Sardinops sagax). The catch plummeted from 292,000 t in 1988–89 to 7000 t in 1991–92 and 1992–93. This caused a serious economic crisis in the local fishery fleet and industry, and resulted in the loss of 3000 jobs. In 1993–94 the fishery showed signs of recovery as the abundance of the Pacific sardine began to recover. The catch improved to 128,000 t in 1993–94 and further to 215,000 t in 1996–97. In trying to understand this great variability, we proposed the hypothesis that the distribution and the abundance of the Pacific sardine of the Gulf of California is determined by the wind patterns (upwelling) and the sea surface temperature. The results of analyzing data from 25 cruises showed the period of low relative abundance between 1990 and 1993 and one of high abundance between 1993 and 1996. The range of the sardine's distribution expanded as its abundance increased and contracted when abundances were low. The relationship between the abundances of the sardine and environmental variables proved to nonlinear and bell-shaped. The adjusted pattern explained 78.8% of the variability of the sardine abundance. The highest abundance are produced by moderate upwelling (13–18 m³s⁻¹ per 10 m of coastline) and sea surface temperatures of between 19°C and 25°C.     

Experimental research on the statistical characteristics of radar signal spikes at small grazing angles

The results of research into the characteristics of spikes in a radar signal scattered by the sea surface at small grazing angles and horizontal polarization of radiation are reported. Studies were carried out at wind speeds of 4.4–10.3 m s^–1. It is shown that the experimental distribution of the amplitudes of the radar signal spikes is described well by an exponential function. The distribution of the amplitudes of spikes, normalized to their average value for the measurement cycle, is of universal character. It does not depend on the observational conditions. The relationship between the average level of signal and the wind speed agrees with the results of computation on the model suggested.Translated by Mikhail M. Trufanov.     

Turbulent free-surface flow around a Wigley hull using the slightly compressible flow formulation

The slightly compressible flow formulation is applied to the free-surface, three-dimensional turbulent flow around a Wigley hull. Two turbulence models (large eddy simulation and Baldwin–Lomax) are used and compared. The simulation conditions are the ones for which experimental and numerical results exist. The computational grid is built using an algebraic grid generator with the model fixed in space. The codes use the interface-capturing technique for computing the free-surface displacements and the Beam and Warming scheme for marching in time the numerical model. The results compare well with the experimental data available.     

Sailboat propeller drag

All but the smallest classes of modern keelboats are fitted with inboard engines and consequently, when making way under sail, the craft experience parasitic drag due to trailing propellers and associated appendages. The variety of screw configurations used on sailing boats includes fixed-blade, feathering, and folding set-ups, with blades numbering two or three. Although the magnitude of the resultant drag is thought to have a significant influence on sailing performance, the published literature having regard to this problem is sparse. Here, the aim was to evaluate the drag effect of fixed-blade propellers of types commonly used on sailing craft. The results of towing tank tests on full-scale propellers are presented for the locked shaft condition; these are presented along with reconfigured data from the few previously published sources. For the case in which the propeller is allowed to rotate, tests were conducted on a typical screw with a range of braking torques being applied. It was hypothesised that the performance coefficients of the Wageningen B-Screw Series could be used to characterise adequately the types of screw of interest and that these could be extrapolated to enable prediction of the drag of a freewheeling propeller; an assessment of this formed part of the investigation.     

Cross-sectional areas of mid-ocean ridge axes bounding the Easter and Juan Fernandez microplates

We have calculated cross-sectional areas for the ridges bounding the Easter and Juan Fernandez microplates, 22°–28°S and 31°–35°S, obtaining accurate results where complete bathymetric data exist and estimates in other regions with partial bathymetric coverage and predicted bathymetry. We consider the reliability and usefulness of global predicted bathymetry in these calculations and the possible application of this dataset in other localities. The spreading rates on ridges bounding these microplates span the range from slow to superfast, allowing an investigation of ridge axis inflation over most of the rates active on Earth today. The across-axis areas of the Easter microplate ridge axes range from –29 km² to 7 km², while the Juan Fernandez ridge axis areas range from –27 km² to 8 km². Positive values correlate with regions usually interpreted as magmatically robust. Negative values arise from calculations in areas of propagating rift tips and deep grabens, such as Pito and Endeavor Deeps. Geochemical trends of Easter microplate axial basalts show decreasing MgO toward propagating rift tips and slight positive correlations between variables such as MgO vs. cross-sectional area, Na_8.0 vs. axial depth, and Na_8.0 vs. cross-sectional area. We document the decrease in the axial area approaching segment ends and propagating rift tips along both the West and East ridges of the microplates. On the Easter microplate both East and West ridge systems undergo large variations in spreading rate from >130 km Myr^–1 to <50 km Myr^–1. Inflation on these ridge segments is highly variable and only weakly correlated with spreading rate. On the Juan Fernandez microplate, West ridge spreading rates vary only between 115–140 km Myr^–1 and are systematically faster than on the East ridge, where rates vary between 10–35 km Myr^–1. Cross axis areas are systematically greater and significantly less variable on the faster spreading West ridge. Overall, compared to oceanic spreading centers bounding major plates with similar spreading rates, the axial areas are smaller on the microplate ridge systems, possibly because their rapidly changing configurations create a lag in the mantle response to the rigid plate boundary.     

Dynamic analysis of vehicle–deck interactions

This paper presents the review and studies at various levels of problems concerning the authors’ previous research on the dynamics of vehicle–deck interactions. The various levels of study include the dynamic structural behavior of vehicle–deck systems, vehicle vibrations, damping effects of vehicles on structural systems, dynamic interactions between tire and deck surface, and vehicle securing on decks during ship motions, etc. The study includes analytical, numerical and experimental analysis. Practical problems encountered by Ro–Ro ship designers are addressed by discussing those analysis. It is shown that influences from the dynamics of vehicle–deck interactions are relevant to a number of aspects of issues, such as the excitation frequency range, how detailed information of the structural system response is required, the structure characteristics, and positions and orientations of vehicles on decks, etc. The study contributes to the knowledge for the naval architect and vehicle engineer on how significant the dynamics of vehicle–deck interactions are when dealing with relevant problems.     

Downwind sail aerodynamics: A CFD investigation with high grid resolution

A computational fluid dynamics (CFD) code was applied to an America's Cup Class Yacht to investigate sailing performance in a downwind configuration. Apparent wind angles at 45°, 105° and 120° are reported, sailed with mainsail and asymmetrical spinnakers. Numerical results are in good agreement with wind tunnel data. A large mesh investigation was performed, ranging from 60,000 elements up to 37 million elements, which shows a converging trend to the experimental values with differences smaller than 3% in both lift and drag. The most commonly used turbulence models in sail applications were tested and the results are presented here in two meshes with 1 million elements and 6.5 millions, respectively. All turbulence models over-estimate forces more than solving the Navier–Stokes system without any additional equations, hence turbulence models do not increase solution accuracy according to these results.     

Observational Study on the Downward Solar Radiation at the Sea Surface in the Western Pacific

In situ observation of downward solar radiation in the Western Pacific were carried out with voluntary merchant ships for five years from autumn 1990 through autumn 1995. Daily means of the short wave radiation were computed from the observed solar radiation. Then, the effects of shadows of the ship's superstructures on the observed radiation were corrected if needed. A 5-year average of short wave radiation along the main sea-lanes in the Western Pacific was calculated based on the observed daily mean solar radiation. Maximum values of 270–280 Wm^–2 are found around 15°–20°N in May and June, while those of 290 Wm^–2 are observed south of 18°S in November and December along the lanes. Small annual variations are found in the equatorial region. Annual mean values at the equator are about 230 Wm^–2 between New Guinea and Indonesia, and 200 Wm^–2 east of New Guinea. The 5-year average of short wave radiation was compared with the climatologies given in previous studies. We have concluded that some of results of previous studies are significantly underestimated.