共查询到20条相似文献，搜索用时 46 毫秒
1.
Further developments and applications of the 2D harmonic polynomial cell (HPC) method proposed by Shao and Faltinsen [22] are presented. First, a local potential flow solution coupled with the HPC method and based on the domain decomposition strategy is proposed to cope with singular potential flow characteristics at sharp corners fully submerged in a fluid. The results are verified by comparing them with the analytical added mass of a double-wedge in infinite fluid. The effect of the singular potential flow is not dominant for added mass and damping, but the error is non-negligible when calculating mean wave loads using direct pressure integration. Then, the double-layer nodes technique is used to simulate a thin free shear layer shed from lifting bodies, across which the velocity potential is discontinuous. The results are verified by comparing them with analytical results for steady and unsteady lifting problems of a flat plate in infinite fluid. The latter includes the Wagner problem and the Theodorsen functions. Satisfactory agreement with other numerical results is documented for steady linear flow past a foil and beneath the free surface. 相似文献
2.
Momentum and gravity effects during the constant velocity water entry of wedge-shaped sections 总被引：1，自引：0，他引：1
Computational fluid dynamics analysis was used to investigate the added mass momentum, flow momentum and gravity effects during the constant velocity water entry of wedge-shaped sections with deadrise angles from 5° to 45°. It is shown that the added mass continues to increase for a time after chine immersion and that added mass can be estimated in terms of a constant added mass coefficient and an effective wetted width. A momentum theory is presented in which the water entry force is explained as the sum of the rate of change of added mass momentum, which becomes zero at immersion to chine depth ratios greater than about three, and the rate of change of flow momentum, which continues at deep immersions. The effect of gravity on the water entry force is given as the hydrostatic force together with the force necessary to create the potential energy in the water pile up. Hydrodynamic forces are not significantly changed by the effect of gravity on the flow fields. 相似文献
3.
A computational model for catamaran wet deck slamming is developed on the basis of the variation of added mass as the hulls enter the water. In the case of a wave-piercing catamaran the bow cross section has a double arch cross section and slamming occurs when the arches fill. Residual air is entrained at the top of the arch due to bubble formation by turbulent mixing and this modifies the effect of the water added mass on the hull. The computational model therefore introduces a soft connection between the water added mass associated with the slam and the hull. The method has been evaluated by comparison with two-dimensional model drop tests in terms of the maximum forces and acceleration imposed on the hull, the variation of velocity during the slam event and the depth of penetration into the water. It is concluded that the added mass computation is adequate for slam modelling in global motions and loads calculations since it gives a good representation of the maximum total forces on the section and their duration. 相似文献
4.
A plate submerged at a certain depth underneath the sea surface has been proposed as a structure type for different purposes, including motion response reduction, wave control, and wave energy harvesting. In the present study, the three-dimensional wave radiation problem is investigated in the context of the linear potential theory for a submerged ring plate in isolation or attached to a floating column as an appendage. In the latter case, the ring plate is attached at a certain distance above the column bottom. The structure is assumed to undergo a heave motion. An analytical model is developed to solve the wave radiation problem via the eigenfunction expansion method in association with the region-matching technique. With the velocity potential being available, the hydrodynamic coefficients, such as added mass and radiation damping, are obtained through the direct pressure integration. An alternative solution of radiation damping has also been developed in this study, in which the radiation damping is related to the Kochin function in the wave radiation problem. After validating the present model, numerical analysis is performed in detail to assess the influence of various plate parameters, such as the plate size and submergence depth. It is noted that the additional added mass due to the attached ring plate is larger than that when the plate is in isolation. Meanwhile, the radiation damping of the column for the heave motion can vanish at a specific wave frequency by attaching a ring plate, corresponding to a condition that there exist no progressive waves in the exterior region. 相似文献
5.
The multiscale method is used to obtain asymptotic expansions up to the quantities of the third order for the elevations of the surface of the basin and the velocity potential of motion of liquid particles in the wave disturbances formed in the process of nonlinear interaction of periodic running waves of the first and second harmonics in a homogeneous ideal incompressible liquid of constant finite depth covered with broken ice. The dependences of the amplitude-phase structure of disturbances on the ice thickness, depth of the basin, and the parameters of interacting harmonics are investigated. We estimate the error of evaluation of the characteristics of the formed vertical displacement of the surface of the basin and nonlinear mass transfer introduced by neglecting the curvature of the wave profile in the expression for the velocity potential in deducing the kinematic and dynamic surface boundary conditions for nonlinear approximations. 相似文献
6.
Akitsugu Nadai Hiroshi Kuroiwa Masafumi Mizutori Shinnchi Sakai 《Journal of Oceanography》1999,55(1):13-30
The Communications Research Laboratory (CRL) has been developing high-frequency ocean surface radars (HFOSRs). The CRL dual-site HFOSR system can clarify the distribution of surface currents with a nominal range of 50 km. This paper presents a theoretical and experimental analysis of the measurement error of the current vector obtained by the CRL HFOSR system, using a comparison of instantaneous current vectors acquired by the HFOSR system and current meters moored at a depth of 2 m, taking account of the vertical current shear. The theoretical analysis shows that the probability distribution of the measurement error of the current vector forms concentric ellipses at a spatial scale that depends on the RMS measurement error of radial current velocity and with an aspect ratio that depends only on the azimuthal difference of the radar beams. When the azimuthal difference is a right angle, the measurement error of the current vector is at a minimum. A comparison between instantaneous current vectors measured by the CRL HFOSR system and moored current meters shows that the distribution of the difference vector between the radar current and the meter current agrees well with the theoretical measurement error of the current vector and that the RMS of difference vector length is about 10 cm s–1 while the azimuthal difference between two radar beams is between 45 and 135 degrees. The accuracy of current measurement by the dual-site HFOSR system is therefore considered to be less than 10 cm s–1 in this range of azimuthal difference. The theoretical analysis will be applicable for a wider range of the azimuthal difference of the radar beams. 相似文献
7.
N. G. Yakovlev 《Izvestiya Atmospheric and Oceanic Physics》2009,45(4):478-494
This paper presents the results of reconstructing the state of ice and snow covers on the Arctic Ocean from 1948 to 2002 obtained
with a couplod model of ocean circulation and sea-ice evolution. The area of the North Atlantic and Arctic Ocean north of
65° N, excluding Hudson Bay, is considered. The monthly mean ice areas and extents are analyzed. The trends of these areas
are calculated separately for the periods of 1970–1979, 1979–1990, and 1990–2002. A systematic slight underestimation by the
model is observed for the ice extent. This error is estimated to fit the error of 100 km in determining the position of the
ice edge (i.e., close to the model resolution). In summer the model fails to reproduce many observed polynias: by observational
data, the ice concentration in the central part of the Arctic Ocean constitutes around 0.8, while the model yields around
0.99. The average trend for the area of ice propagation in 1960–2002 is 13931 km2/year (or approximately 2% per decade); the trend of the ice area is 17643 km2/year (or approximately 3% per decade). This is almost three times lower than satellite data. The calculated data for ice
thickness in the late winter varies from 3.5 to 4.8 m, with a clear indication of periods of thick ice (the 1960s–1970s) and
relatively thin ice (the 1980s); 1995 is the starting point for quick ice-area reduction. The maximum ice accumulation is
in 1977 and 1988; here, the average trend is negative: −121 km3/year (or approximately 5.5% per decade). In 1996–2002, the average change in the ice thickness constituted +1.7 cm/year.
This speaks to the relatively fast disappearance of thin-ice fractions. This model also slightly underestimates the snow mass
with a trend of −2.5 km3/year (almost 0.35 mm of snow per year or 0.1 mm of liquid water per year). An analysis of the monthly mean ice-drift velocity
indicates the good quality of the model. Data on the average drift velocity and the results of comparisons between the calculated
and satellite data for individual months are presented. A comparison with observational data from 1990–1996 in the Fram Strait
shows that the model yields 3.28 m for the average ice thickness against the observed value of approximately 3.26 m. For the
same period, the model yields a monthly mean transport of 291.29 km3 as compared to the observed value of 237.17 km3. A comparison between the measured and calculated drift velocities in the Fram Strait indicates that the model value is around
9.78 cm/s, which is comparable to the measured value of 10.2 cm/s. The existing problems with describing the ice redistribution
by thickness gradations are illustrated by comparing data on ice thickness in the Fram Strait. 相似文献
8.
Correction method for full-depth current velocity with lowered acoustic Doppler current profiler (LADCP) 总被引：1，自引：0，他引：1
A new method is presented to process and correct full-depth current velocity data obtained from a lowered acoustic Doppler
current profiler (LADCP). The analysis shows that, except near the surface, the echo intensity of a reflected sound pulse
is closely correlated with the magnitude of the difference in vertical shear of velocity between downcast and upcast, indicating
an error in velocity shear. The present method features the use of echo intensity for the correction of velocity shear. The
correction values are determined as to fit LADCP velocity to shipboard ADCP (SADCP) and LADCP bottom-tracked velocities. The
method is as follows. Initially, a profile of velocity relative to the sea surface is obtained by integrating vertical shears
of velocity after low-quality data are rejected. Second, the relative velocity is fitted to the velocity at 100–800 dbar measured
by SADCP to obtain an “absolute” velocity profile. Third, the velocity shear is corrected using the relationship between the
errors in velocity shears and echo intensity, in order to adjust the velocity at sea bottom to the bottom-tracked velocity
measured by LADCP. Finally, the velocity profile is obtained from the SADCP-fitted velocity at depths less than 800 dbar and
the corrected velocity shear at depths greater than 800 dbar. This method is valid for a full-depth LADCP cast throughout
which the echo intensity is relatively high (greater than 75 dB in the present analysis). Although the processed velocity
may include errors of 1–2 cm s−1, this method produced qualitatively good current structures in the Northeast Pacific Basin that were consistent with the
deep current structures inferred from silicate distribution, and the averaged velocities were significantly different from
those calculated by the Visbeck (2002) method. 相似文献
9.
The horizontal movement of inertial particles in the intensive vortices, where the centrifugal force can be substantially higher than the gravity, is studied analytically. A similar problem was studied earlier for small (Stokes) particles at low Reynolds number, which allow one to be limited to the linear resistance law. It is shown that the previous results to a great extent can be extrapolated to the case of considerably heavier particles (e.g., water droplets with a diameter up to 1 mm at Reynolds numbers up to 103). The nonlinear nature of the resistance, i.e., its dependence on the particle velocity relative to the medium, should be taken into account for such particles. Some general laws are established for particle dynamics. In particular, their tangential velocity is close to the velocity of the medium, while the radial velocity is substantially lower (it is close on the order of magnitude to the geometric mean of the particle tangential velocity and the difference between the latter and the tangential velocity of the medium). The limits of applicability of the results are found, i.e., the restrictions to the size and mass/density of particles. 相似文献
10.
国内外在声学多普勒测速误差分析和修正算法方面已作了大量的研究工作.一般认为JANUS配置方式下,多普勒频移公式的简化误差可忽略.文中从理论上对多普勒测速的原理公式的简化过程进行分析,研究了大倾角状态下的测速公式的相对误差,认为在高速、大倾角工况下公式简化误差不可忽略,并提出了修正公式. 相似文献
11.
A methodology for the correction of systematic errors in a simplified atmospheric general‐circulation model is proposed. First, a method for estimating initial tendency model errors is developed, based on a 4‐dimensional variational assimilation of a long‐analysed dataset of observations in a simple quasi‐geostrophic baroclinic model. Then, a time variable potential vorticity source term is added as a forcing to the same model, in order to parameterize subgrid‐scale processes and unrepresented physical phenomena. This forcing term consists in a (large‐scale) flow dependent parametrization of the initial tendency model error computed by the variational assimilation. The flow dependency is given by an analogues technique which relies on the analysis dataset. Such empirical driving causes a substantial improvement of the model climatology, reducing its systematic error and improving its high frequency variability. Low‐frequency variability is also more realistic and the model shows a better reproduction of Euro‐Atlantic weather regimes. A link between the large‐scale flow and the model error is found only in the Euro‐Atlantic sector, other mechanisms being probably the origin of model error in other areas of the globe. 相似文献
12.
Toshihiko Teramoto 《Journal of Oceanography》1971,27(1):7-19
A formula to determine an effective electric conductivity of channel-bed from measurements of water velocity, cross-channel potential difference and electric current-density is presented. The effective conductivity defined in this paper is the conductivity averaged over depth from channel-bed surface to a depth of the order comparable to a breadth of the channel. Effective conductivities for the Izu Island region and the Straits of Tsugaru are estimated as 2.2×10–2 mhos/m and 6.0×10–2 mhos/m, which are considerably large and very large, respectively, in comparison to most of those determined for several regions near England. High conductivities for these regions are consistent with the distribution of underground conductivity anomaly deduced from that of geomagnetic field anomaly observed at the earth's surface in Japan. An influence of underground conductivity upon measurements of sea surface water velocity with towed electrodes is estimated for the regions. 相似文献
13.
In the present study,analyzed are the variation of added mass for a circular cylinder in the lock-in(synchronization) range of vortex-induced vibration(VIV) and the relationship between added mass and natural frequency.A theoretical minimum value of the added mass coefficient for a circular cylinder at lock-in is given.Developed are semi-empirical formulas for the added mass of a circular cylinder at lock-in as a function of flow speed and mass ratio.A comparison between experiments and numerical simulations shows that the semi-empirical formulas describing the variation of the added mass for a circular cylinder at lock-in are better than the ideal added mass.In addition,computation models such as the wake oscillator model using the present formulas can predict the amplitude response of a circular cylinder at lock-in more accurately than those using the ideal added mass. 相似文献
14.
The sway, heave and roll added masses of three uniform cylinders with semi-circular, rectangular and triangular cross-sectional shapes in shallow and narrow water are numerically analysed. The method is based on simulation of the potential flow induced by the cylinder's mode of motion. The effects of shallow and narrow water on added mass are analysed and presented. It is concluded that the shallow and narrow water effects on added mass depend on the different cross-section shapes of the cylinders. In particular, the water depth effect on sway added mass is stronger than that on heave added mass while the narrow water effect on sway is weaker than that on heave. The shallow water effect on added mass tends to weaken the narrow water effect. Lastly the effect of shallow and narrow water on added mass on a rectangular cylinder is the strongest while that on a triangular cylinder is the weakest. 相似文献
15.
Ocean prediction systems rely on an array of assumptions to optimize their data assimilation schemes. Many of these remain untested, especially at smaller scales, because sufficiently dense observations are very rare. A set of 295 drifters deployed in July 2012 in the north-eastern Gulf of Mexico provides a unique opportunity to test these systems down to scales previously unobtainable. In this study, background error covariance assumptions in the 3DVar assimilation process are perturbed to understand the effect on the solution relative to the withheld dense drifter data. Results show that the amplitude of the background error covariance is an important factor as expected, and a proposed new formulation provides added skill. In addition, the background error covariance time correlation is important to allow satellite observations to affect the results over a period longer than one daily assimilation cycle. The results show the new background error covariance formulations provide more accurate placement of frontal positions, directions of currents and velocity magnitudes. These conclusions have implications for the implementation of 3DVar systems as well as the analysis interval of 4DVar systems. 相似文献
16.
The highly accurate Boussinesq-type equations of Madsen et al. (Madsen, P.A., Bingham, H.B., Schäffer, H.A., 2003. Boussinesq-type formulations for fully nonlinear and extremely dispersive water waves: Derivation and analysis. Proc. R. Soc. Lond. A 459, 1075–1104; Madsen, P.A., Fuhrman, D.R., Wang, B., 2006. A Boussinesq-type method for fully nonlinear waves interacting with a rapidly varying bathymetry. Coast. Eng. 53, 487–504); Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945) are re-derived in a more general framework which establishes the correct relationship between the model in a velocity formulation and a velocity potential formulation. Although most work with this model has used the velocity formulation, the potential formulation is of interest because it reduces the computational effort by approximately a factor of two and facilitates a coupling to other potential flow solvers. A new shoaling enhancement operator is introduced to derive new models (in both formulations) with a velocity profile which is always consistent with the kinematic bottom boundary condition. The true behaviour of the velocity potential formulation with respect to linear shoaling is given for the first time, correcting errors made by Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945). An exact infinite series solution for the potential is obtained via a Taylor expansion about an arbitrary vertical position z = zˆ. For practical implementation however, the solution is expanded based on a slow variation of zˆ and terms are retained to first-order. With shoaling enhancement, the new models obtain a comparable accuracy in linear shoaling to the original velocity formulation. General consistency relations are also derived which are convenient for verifying that the differential operators satisfy a potential flow and/or conserve mass up to the order of truncation of the model. The performance of the new formulation is validated using computations of linear and nonlinear shoaling problems. The behaviour on a rapidly varying bathymetry is also checked using linear wave reflection from a shelf and Bragg scattering from an undulating bottom. Although the new models perform equally well for Bragg scattering they fail earlier than the existing model for reflection/transmission problems in very deep water. 相似文献
17.
Walter J. Senus 《Marine Geodesy》2013,36(3):263-288
The current accuracy of sea‐going and airborne gravity measurements is not bounded by the precision of the gravimeter but by the precision with which external parameters such as vehicle velocity (speed), azimuthal heading, and geographic position can be determined. Uncertainties in these parameters are summed up in the Eötvös correction in the reduction of the measured data. This work investigates the suitability of baseline navigation systems, in particular Loran‐C and Omega, to further reduce the uncertainty of the Eötvös correction. Emphasis is placed on the velocity measurement error. A new algorithm is developed which derives velocity based upon the change in hyperbolic (or circular) grid readings, as opposed to the standard change in geographic position technique. A comparative analysis shows the new algorithm to be as precise as the currently used conventional calculation. Further, this simplified technique is accompanied with a 20‐fold reduction in computational complexity. Application of the results presented in this paper to the Loran‐C and Omega systems shows a velocity determination capability of 0.1 knot over a six minute integration time, and 0.05 knot for a 15 minute integration time for Loran‐C. The minimum error attainable in Omega is 0.1 knot when determined over a 17 minute measurement time. Further precision can be gleaned by applying the calculated error as a correction to subsequent velocity calculations. 相似文献
18.
This paper deals with hydrodynamic forces of a single semisubmerged circular cylinder containing a concentric cylindrical hole constrained to move in a water domain of finite depth. The fluid domain is divided into inner and outer regions. The Laplace equations governing velocity potentials for the three regions are solved by separation of variables and expressed in terms of eigenfunctions of the resulting equations which satisfy appropriate boundary conditions. Continuity of pressure and velocity at the interface of the regions provides the necessary equations from which the velocity potentials, pressures and forces are obtained. Numerical results are plotted for added mass and damping coefficients for different draft-to-depth and radius-to-depth values and for various wave amplitudes. 相似文献
19.
The Norwegian wave-power buoy1 consists of a half-immersed floating sphere which is open to the sea at the bottom end. It is a two degree-of-freedom device involving the independent motion of the outer rigid sphere and the pressure across the internal free surface. A simpler model of the device is to represent the oscillatory flow through the bottom opening as another rigid body motion being that of the curved surface which would complete the sphere.The wave-induced forces on this surface and also the outer spherical surface due to independent oscillations of either are determined semi-analytically using a simple extension of the method of Havelock,2 recently simplified and generalised by Hulme,3 for the complete half-immersed sphere. In particular the 2 × 2 added mass and damping matrices are determined as a function of frequency and relative size of bottom opening to sphere radius. These quantities are essential in any theoretical analysis of a multi-degree of freedom wave-energy device. 相似文献
20.
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. 相似文献