Numerical and experimental study of a flapping foil generator

The energy extraction performance of a flapping foil generator is studied through experiment and numerical simulation. A practical flapping foil generator has been proposed. The heave motion and the pitch motion of the foil are adjusted through a crankshaft-like structure. The heave and pitch motions of the foil are transferred to the rotational motion of the main shaft. A pair of gears is adopted to increase the pitch angle. A prototype with pitch amplitude θ₀ = 60^∘ has been built and the experiment is carried out in a tunnel. The overall performance of the mechanism has been analysed. Good agreement of numerical results and experiment data has been found. Further simulations with larger pitch amplitudes are carried out. It is found that higher efficiency can be achieved with larger pitch amplitude at medium frequency.     

Computations of self-propulsion free to sink and trim and of motions in head waves of the KRISO Container Ship (KCS) model

Two computations of the KCS model with motions are presented. Self-propulsion in model scale free to sink and trim are studied with the rotating discretized propeller from the Hamburg Model Basin (HSVA) at Fr = 0.26. This case is particularly complex to simulate due to the close proximity of the propeller to the rudder. The second case involves pitch and heave in regular head waves. Computations were performed with CFDShip-Iowa version 4.5, a RANS/DES CFD code designed for ship hydrodynamics. The self-propulsion computations were carried out following the procedure described in Carrica et al. [1], in which a speed controller is used to find the propeller rotational speed that results in the specified ship velocity. The rate of revolutions n, sinkage, trim, thrust and torque coefficients K_T, K_Q and resistance coefficient CT_(SP) are thus obtained. Comparisons between CFD and EFD show that the rate of revolutions n, thrust and torque coefficients K_T and K_Q have higher prediction accuracies than sinkage and trim. For the simulation of pitch and heave in head waves, the geometry includes KCS hull and rudder under three conditions with two Froude numbers and three wave length and amplitude combinations. 0th and 1st harmonic amplitudes and 1st harmonic phase are computed for total resistance coefficient C_T, heave motion z and pitch angle θ. Comparisons between CFD and EFD show that pitch and heave are much better predicted than the resistance. In both cases comparisons with simulations by other authors presented at the G2010 CFD Workshop [2] using different CFD methodologies are included.     

Calorimetric study of the thermodynamics of formation of MgSO4 and NaSO4− ion pairs at the ionic strength of seawater

The formation of the ion pairs MgSO₄ and NaSO₄^? was investigated calorimetrically at 0.75 M ionic strength, 25°C, 1 atm. Simultaneous determinations of enthalpy changes, ΔH₁⁰, and stability constants, K₁, were not possible, and values of K₁ determined independently had to be introduced for the calculation of ΔH₁⁰. The values of ΔH₁⁰ obtained were 1–3 kJ mol^?1 for MgSO₄ and 0 kJ mol^?1 for NaSO₄^?.     

A new formula for saturated water steam pressure within the temperature range −25 to 220°C

Instead of approximation formula ln(E(t)/E(0)) = [(a ? bt)t/(c + T)] commonly used at present for representing dependence of pressure of saturated streams of liquid water E upon temperature we suggested new approximation formula of greater accuracy in the form ln(E(t)/E(0)) = [(A ? Bt + Ct ²)t/T], where t and T are temperature in °C and K respectively. For this formula with parameters A = 19.846, B = 8.97 × 10^?3, C = 1.248 × 10^?5 and E(0) = 6.1121 GPa with ITS-90 temperature scale and for temperature range from 0°C to 110°C relative difference of approximation applying six parameter formula by W. Wagner and A. Pruß 2002, developed for positive temperatures, is less than 0.005%, that is approximately 15 times less than accuracy obtained with the firs formula. Increase of temperature range results in relative difference increasing, but for even temperature range from 0°C to 220°C it does not higher than 0.1%. For negative temperatures relative difference between our formula and a formula of D. M. Murphy and T. Koop, 2005, is less than 0.1% for temperatures higher than ?25°C. This paper also presents values of coefficients for approximation of Goff and Grach formula recommended by IMO. The procedure of finding dew point T _d for known water steam pressure e _n based on our formula adds up to solving an algebraic equation of a third degree, which coefficients are presented in this paper. For simplifying this procedure this paper also includes approximation ratio applying a coefficient A noted above, in the form T _d (e _n ) = \(\frac{{AT_0 }}{{A - \varepsilon }}\) + 0.0866?² + 0.0116?^10/3, where ? = ln(e _n /E(T ₀)). Error of dew point recovery in this ratio is less than 0.005 K within the range from 0 to 50°C.     

摆动尾鳍水动力性能的试验和数值研究

鱼类能够在水下高速度、低噪音、高效率地游动。鱼类出色的推进性能通过其摆动尾鳍实现。这种摆动尾鳍推进方式已经用在了水下无人航行器上。因此研究摆动尾鳍的水动力性能是非常有意义的。对摆动尾鳍的推进水动力性能进行了详尽的研究。设计、装配了一套仿尾鳍推进系统,并对其进行了相应的水动力试验。在试验中研究了运动参数对摆动尾鳍水动力性能的影响。与此同时,采用基于雷诺平均N-S方程的数值方法对摆动尾鳍的水动力性能进行了研究。在数值计算中采用了k-ωSST湍流模型和有限体积法。数值计算结果和水动力试验结果进行了比较。对尾鳍表面的压力分布和流场中的尾涡结构进行了分析。水动力试验和数值计算都表明摆动尾鳍可以产生推进力和较高的推进效率。     

Physico-chemical study of Dead Sea waters: II. Density measurements and equation of state of Dead Sea waters at 1 atm

The densities (p) of artificial solutions of Dead Sea waters have been measured at 20°, 25°, 30°, 35° and 40°C using a vibrating tube densimeter. The resulting relative densities (p — p_o, where p_o is the density of water) have been used to determine an equation of state for Dead Sea waters. The average deviation between experimental values and those calculated from the obtained equation was 0.000020 g ml^?1. A thermal expansion coefficient and the coefficients characterizing the influence of the changes in salt or ionic concentrations on the density of Dead Sea waters were calculated, and they were shown to be temperature and concentration dependent. The densities of Dead Sea waters were found to be very sensitive to any changes in ionic composition. The partial molal volumes of salt components were calculated and discussed.     

Oceanic processes of upper ocean heat content associated with two types of ENSO

The spatiotemporal variability of equatorial Pacific upper ocean heat content (HC) and subsurface heat during two types of El Niño-Southern Oscillation (ENSO), namely eastern and central Pacific (EP and CP) types, is investigated using subsurface ocean heat budget analysis. Results show that HC tendencies during both types of ENSO are mainly controlled by oceanic heat advection beneath the mixed layer to the thermocline, and the role of net surface heat flux can be neglected. The most important three terms are the zonal and vertical advections of anomalous heat by climatological currents (QU ₀ T′, QW ₀ T′) and zonal advection of climatological heat by anomalous current (QU′T ₀). The large contribution of QU ₀ T′ extends from west to east along the equatorial Pacific. The considerable contribution of QU′T ₀ is confined to the east of 160°W, and that of the QW ₀ T′ is observed in the central Pacific between 180°E and 120°W. In particular, a major contribution of QW ₀ T′ is also observed in the far eastern Pacific east of 100°W during EP ENSO. There is also a small contribution from meridional advection of climatological heat by anomalous current (QV′T ₀). In contrast, the meridional advection of anomalous heat by climatological currents (QV ₀ T′) and vertical advection of climatological heat by anomalous current (QW′T ₀) are two damping factors in the HC tendency, with the former dominating. Differences in spatial distribution of the heat advection associated with the two types of ENSO are also presented. We define a warm water heat index (WWH) as integrated heat content above 26 kg m^?3 potential density (26σ _? ) isopycnal depth within 130°E–80°W and 5°S–5°N. Further examination suggests that the recharge–discharge of WWH is involved in both types of El Niño, though with some differences. First, it takes about 42 (55) months for the evolution of a recharge–discharge cycle during an EP (CP) ENSO. Second, the EP El Niño event peaks during the discharge phase, 7–8 months after the recharge time. The CP El Niño peaks during the recharge phase, 4–5 months before the recharge time. The locations of HC anomalies in the El Niño mature phase relative to those at recharged time explain why the EP and CP El Niño peak in different stages of the recharge–discharge process.     

Stability of ion pairs from gypsum solubility degree of ion pair formation between the major constituents of seawater

The stability constants of the ion pairs NaSO₄^?, KSO₄^?, MgSO₄^?, CaSO₄, MgCl⁺ and CaCl⁺ were determined at 25°C and 0.7 M formal ionic strength, by measuring the solubility of gypsum (CaSO₄ · 2H₂O) in different media. The media used contained one or two of the following electrolytes: NaCl, KCl, MgCl₂, NaClO₄, Mg(ClO₄)₂, Na₂SO₄. Values for the stability constants are 1.22, 1.84, 12.3, 30.6, 0.48 and 1.20 M^?1, respectively, and the solubility product for gypsum is 2.87 · 10^?4M². The distribution of the main constituents of seawater was calculated using these results and the values of the carbonate and bicarbonate constants given by Dyrssen and Hansson (1972–1973). The solubility of gypsum in seawater as calculated and determined experimentally was 21.43 mM and 21.10 mM, respectively.     

Growth characteristics ofChattonella antiqua

Nutrient requirements of a red tide flagellate,Chattonella antiqua, were investigated in a laboratory culture experiment. Growth ofC. antiqua was supported by nitrate and ammonium, and by urea to a limited extent, but not by glycine, alanine and glutamate. Orthophosphate served as a good phosphorus source but glycerophosphate did not. Fe³⁺ (1µM) fully promoted the flagellate's growth in the presence of 80µM of EDTA. The addition of Mn²⁺ (0–20µM), Zn²⁺ (0–10µM) and Co²⁺ (0–0.4µM) did not show any effect. Among three vitamins tested, only B₁₂ (6 ng 1^?1) served as a growth factor. Glucose, acetate and glycolate did not improve growth in the light nor did they support growth in darkness. The minimum cell quotas for nitrogen, phosphorus, iron and B₁₂ were estimated to be 11 pmoles ce^?1, 1.0, ～0.09 and 1.1 fg cell^?1, respectively.     

Thickness of the mixed bottom layer in the Northern Atlantic

The thickness of the mixed bottom boundary layer (BBL) has been analyzed based on the CTD data at transoceanic sections in abyssal waters of the Northern Atlantic. The measurements were carried out at two transoceanic sections approximately along 48° N (ASV-99) and 5° N (AI-2000) in 1999 and 2000. These data, and the WOCE data obtained at four zonal sections (AR7E and AR12 along 57° N, AR01 along 24.5° S, and A06 along 7.5° N), were used for the calculation of the statistical characteristics of the BBL??s thickness H _B . The probability distribution function F(H _B ) was close to lognormal. The mean value ??H _B ?? at different latitudes was in the range from 30 to 60 m. The averaged BBL thickness = 46.1 m. The BBL??s thickness was about 1% of the ocean??s depth D; the ratio H _B /D was the minimum (0.8%) near the equator and increased up to 1.6% in the polar latitudes.     

斜向波浪作用下双层水平板式防波堤波浪荷载试验研究

In this study, systematic physical model tests were performed to investigate the wave forces on the twin-plate breakwater under irregular waves. Based on the experimental results, the effects of the relative plate width B/L,wave height Hs/D and incident angle θ0 on the wave forces were analyzed and discussed. The results showed that:(1) The envelopes of the total wave pressure were generally symmetrical along the direction of plate width under the incident angles(θ0) being 0°, 15°, 30°, 45° and 60°. In particular, the envelopes of wave pressure atθ0=30° were larger than all other cases.(2) The synchronous pressure distribution of the breakwater under oblique wave action was more complicated comparing to the normal incident waves.(3) Based on data analysis, an empirical formula was obtained to estimate the total vertical force of the twin-plate breakwater.This empirical formula can be a good reference for the design basis of engineering applications under specified wave conditions.     

Particulate aluminium fluxes in the eastern Atlantic

The atmospheric, primary down-column and sedimentary fluxes of particulate aluminium (Al_p) have been calculated for a number of regions in the Atlantic Ocean.The vertical down-column flux of Al_p from Atlantic surface waters exhibits a strong geographical variation, and its magnitude is influenced by supply mechanisms, which control the surface Al_p concentrations, and primary production, which affects the rate of down-column transport. Overall, the down-column transport of Al_p is greatest in the marginal regions of the Atlantic. In the eastern margins the highest surface water concentrations are found in the region lying between ～30°N and ～10°N, i.e. under the general path of the northeast trades. In this region there is excellent agreement between the dry (i.e. 1 cm^?1 s^?1 deposition velocity) atmospheric flux (～80 000 ng Al_p cm^?2 y^?1), the primary vertical down-column flux ($\text{? 70 000}\text{ng Al}{}_{\mathrm{p}}\text{cm}{}^{\mathrm{?2}}\text{y}{}^{\mathrm{?1}}$) and the sediment flux (～90 000 ng Al_p cm^?2 y^?1). In the regions to the north (i.e. ～40°N to ～30°N) and to the south (i.e. ～10°N to ～5°S) the primary down-column Al_p flux decreases to an average of ～19 000 μg cm^?2 y^?1, which makes a direct maximum contribution of ～20% of the sediment Al_p requirement. In the open-ocean South Atlantic the primary down-column flux of Al_p is ～3300 μg cm^?2 y^?1, this is similar to the dry (i.e. 1 cm^?1 s^?1 deposition velocity) atmospheric flux, and contributes ～20% of the Al_p required by the underlying deep-sea sediment.     

Ship propulsion in waves by actively controlled flapping foils

Flapping wings located beneath or to the side of the hull of the ship are investigated as unsteady thrusters, augmenting ship propulsion in waves. The main arrangement consists of horizontal wing(s) in vertical oscillatory motion which is induced by ship heave and pitch, while rotation about the wing pivot axis is actively controlled. In this work we investigate the energy extraction by the system operating in irregular wave conditions and its performance concerning direct conversion to propulsive thrust. More specifically, we consider operation of the flapping foil in waves characterised by a spectrum, corresponding to specific sea state, taking into account the coupling between the hull and the flapping foil dynamics. The effect of the wavy free surface is accounted for through the satisfaction of the corresponding boundary conditions and the consideration of the wave velocity on the formation of the incident flow. Numerical results concerning thrust and power coefficients are presented, indicating that significant thrust can be produced under general operating conditions. The present work can be exploited for the design and optimum control of such systems extracting energy from sea waves for augmenting marine propulsion in rough seas, with simultaneous reduction of ship responses offering also dynamic stabilisation.     

Phytoplankton production characteristics in the Eastern Atlantic and Atlantic sector of the Southern Ocean in October–November 2004

Production parameters of surface phytoplankton were measured along three transects: La Manche-Cape Town (I); Cape Town-54°S (II); 0°-49°W (along 54°S) (III). The Canary upwelling waters were most productive along transect I, where the surface chlorophyll a (Chl ₀) and the surface primary production (PP ₀) were as high as 4.3 mg/m³ and 173 mg C/m³ per day, respectively. Mosaic patterns in the distribution of these parameters were recorded in the northeastern regions of the South Subtropical Anticyclonic Gyre (Chl ₀ = 0.03–0.35 mg/m³; PP ₀ = 1.6–12.6 mg C/m³ per day). Along transect II, the average twofold southward increase in Chl ₀ (from 0.2 to 0.4 mg/m³) and the concurrent decline of the phytoplankton assimilation activity ( AN ₀) resulted in deviations from typical latitudinal changes inPP ₀. At most sites, PP ₀ values varied between 6 and 15 mg C/m³ per day. Negligible changes in Chl ₀ (0.36–0.85 mg/m³), PP ₀ (8–19 mg C/m³ per day), and AN ₀ (0.7–1.6 mg C/mg chl a per hour) were registered for the oceanic waters along transect III. Along all the transects, PP ₀ depended on Chl ₀ to a greater extent than AN ₀. The values of the latter parameter were largely determined by the water temperature and showed a slight correlation with the insolation. Along transect II, the integrated primary production (PP _int) and the layer-integrated chlorophyll a in the upper 200 m (Chl _0–200) generally varied from 180 to 360 mg C/m² per day and from 30 to 70 mg/m², respectively. In the Polar Front region, an increase in Chl _0–200, PP _int, Chl ₀, and PP ₀ up to respective values of 190 mg/m², 520 mg C/m² per day, 1.2 mg/m³, and 32 mg C/m³ per day was observed. A comparison of the water column (0–100 m) stability with the vertical distribution of the primary production and chlorophyll content along transect II implies that the thick (>100 m) upper mixed layer (UML) formed in response to the strong water cooling and wind forcing was largely responsible for the limited primary production in the Subantarctic and Antarctic regions. The large UML thickness resulted in an intense removal of plant cells from the photosynthetic layer and light starvation of a significant (up to 60%) part of the phytoplankton community.     

Quantitative evaluation of turbulent mixing in the Central Equatorial Pacific

Turbulent mixing in the central equatorial Pacific has been quantitatively evaluated by analyzing data from microstructure measurements and conductivity temperature depth profiler (CTD) observations in a meridionally and vertically large region. The result that strong turbulent mixing with dissipation rate ε (>O(10^?7) W kg^?1), continuing from sea-surface mixed layer to low Richardson number region below, in the area within 1° of the equator, shows that turbulent mixing has a close relationship to shear instability. ε > O(10^?7) W kg^?1 and turbulent diffusivity K _ρ  > O(10^?3) m² s^?1 were obtained from near-surface to 85 db at stations even southwardly beyond 3°S, where it is already far from the southern boundary (~2°S) of the Equatorial Undercurrent. Turbulence-induced heat flux and salinity flux were calculated, and both had their maxima in the equatorial upwelling region, though the former was downward and the latter was upward. Accordingly, vertical velocity in the upwelling region was estimated to be similar to the results derived by other methods. These fluxes and the vertical velocity suggest the critical importance of turbulent mixing in maintaining the well-mixed upper layer. Secondly, in the intermediate region (>500 db), turbulent eddies were investigated by applying Thorpe’s method to the CTD data. A large number of overturns were detected, with spatial-averaged K _ρ (700–1,000 db) being 3.3 × 10^?6 m² s^?1, and the corresponding K _ρ-max reaching to O(10^?4) m² s^?1 in the north (3°–13°N). The results suggest that, in the intermediate region, considerable turbulent mixing occurs and moderates the properties of the water masses.     

Evaluation of spatial distribution of turbulent mixing in the central Pacific

A long-term mean turbulent mixing in the depth range of 200–1000 m produced by breaking of internal waves across the middle and low latitudes (40°S–40°N) of the Pacific between 160°W and 140°W is examined by applying fine-scale parameterization depending on strain variance to 8-year (2005–2012) Argo float data. Results show that elevated turbulent dissipation rate (ε) is related to significant topographic regions, along the equator, and on the northern side of 20°N spanning to 24°N throughout the depth range. Two patterns of latitudinal variations of ε and the corresponding diffusivity (K_ρ) for different depth ranges are confirmed: One is for 200–450 m with significant larger ε and K_ρ, and the maximum values are obtained between 4°N and 6°N, where eddy kinetic energy also reaches its maximum; The other is for 350–1000 m with smaller ε and K_ρ, and the maximum values are obtained near the equator, and between 18°S and 12°S in the southern hemisphere, 20°N and 22°N in the northern hemisphere. Most elevated turbulent dissipation in the depth range of 350–1000 m relates to rough bottom roughness (correlation coefficient?=?0.63), excluding the equatorial area. In the temporal mean field, energy flux from surface wind stress to inertial motions is not significant enough to account for the relatively intensified turbulent mixing in the upper layer.     

Forces on a pitching plate: An experimental and numerical study

A flat plate in pitching motion is considered as a fundamental source of locomotion in the general context of marine propulsion. The experimental as well as numerical investigation is carried out at a relatively small Reynold number of 2000 based on the plate length c and the inflow velocity U_∞. The plate oscillates sinusoidally in pitch about its 1/3 − c axis and the peak to peak amplitude of motion is 20°. The reduced frequency of oscillation k = πfc/U_∞ is considered as a key parameter and it may vary between 1 and 5. The underlying fluid-structure problem is numerically solved using a compact finite-differences Navier–Stokes solution procedure and the numerical solution is compared with Particle Image Velocimetry (PIV) measurements of the flow field around the pitching foil experimental device mounted in a water-channel. A good agreement is found between the numerical and experimental results and the threshold oscillation frequency beyond which the wake exhibits a reverse von Kármán street pattern is determined. Above threshold, the mean velocity in the wake exhibits jet-like profiles with velocity excess, which is generally considered as the footprint of thrust production. The forces exerted on the plate are extracted from the numerical simulation results and it is shown, that reliable predictions for possible thrust production can be inferred from a conventional experimental control volume analysis, only when besides the wake's mean flow the contributions from the velocity fluctuation and the pressure term are taken into account.     

Equilibrium and structural studies of silicon(IV) and aluminium(III) in aqueous solution. II. Formation constants for the monosilicate ions SiO(OH)3− and SiO2(OH)22−. A precision study at 25°C in a simplified seawater medium

The hydrolysis of silicic acid, Si(OH)₄, was studied in a simplified seawater medium (0.6 M Na(Cl)) at 25°C. The measurements were performed as potentiometric titrations (hydrogen electrode) in which OH^? was generated coulometrically. The total concentration of Si(OH)₄, B, and log[H⁺] were varied within the limits 0.00075 ? B ? 0.008 M and 2.5 ? -log[H⁺] ? 11.7, respectively. Within these ranges the formation of SiO(OH)₃^? and SiO₂(OH)₂^2? with formation constants log β_?11(Si(OH)₄ ? SiO(OH)₃^? + H⁺) = ?9.472 ±0.002 and log β_?21(Si(OH)₄ ? SiO₂(OH)₂^2? + 2H⁺) = ?22.07 ± 0.01 was established. With B > 0.003 M polysilicate complexes are formed, however, with $\text{-}\text{log[H}{}^{\mathrm{+}}\text{] ? 10.7}$ their formation does not significantly affect the evaluated formation constants. Data were analyzed with the least squares computer program LETAGROPVRID.     

Seawater: an explanation of differential isothermal compressibility measurements in terms of hydration and ion-water interactions

Hydration, ion-water interactions, and water structure effects in seawater were studied by determining differences (Δβ) between the compressibilities of test salt solutions and the compressibilities of reference solutions. The reference solutions were distilled water and seawater (35%0), and the test salt solutions were either 0.13 m or 0.26m with respect to one of the following test salts: LiCl, NaCl, KCl, CsCl, NaF, NaI, MgCl₂, CaCl₂, BaCl₂, Na₂SO₄, K₂SO₄, and MgSO₄. The compressibility measurements (to 900 bars) were carried out at 2°C and also at 15°C using a differential method in which a pressure increase or a temperature increase causes Δβ to become less negative. At 1 bar and 15°C, the Δβ (0.26 m, distilled water reference) values ranged from ?1.14 × 10^?6 bar^? for NaI to ?3.84 × 10^?6 bar^?1 for Na₂SO₄, and the Δβ (0.26 m, seawater reference) values ranged from ?1.30 × 10^?6 bar^?1 for NaCl to ?3.04 × 10^?6 bar^?1 for Na₂SO₄. The Δβ values were used to calculate hydration numbers. Entropy of transfer, excess hydrogen bond breaking (determined by NMR), and effective radii of ions are properties which can be used to describe the influence of ions on water structure. The extent to which these properties correlate with Δβ values depends upon whether the ion is an anion or a cation, and this correlation forms the thesis that anions alter water structure in a different way than do cations.