共查询到20条相似文献，搜索用时 31 毫秒
1.
A fast multipole boundary element method for three-dimensional potential flow problems 总被引：2，自引：0，他引：2
A fast multipole methodology (FMM) is developed as a numerical approach to reduce the computational cost and memory requirements in solving large-scale problems. It is applied to the boundary element method (BEM) for three-dimensional potential flow problems. The algorithm based on mixed multipole expansion and numeric, al integration is implemented in combination with an iterative solver. Numerical examinations, on Dirichlet and Neumann problems, are carried out to demonstrate the capability and accuracy of the present method. It has been shown that the method has evident advantages in saving memory and computing time when used to solve huge-scale problems which may be prohibitive for the traditional BEM implementation. 相似文献
2.
Shape optimization of two-dimensional cavitators in supercavitating flows, using NSGA II algorithm 总被引：1，自引：0，他引：1
The reduction of energy consumption of high speed submersible bodies is an important challenge in hydrodynamic researches. In this paper, shape optimization of two-dimensional cavitators in supercavitating flows is studied. A two dimensional supercavitation potential flow passes a symmetric two dimensional cavitator, which is placed perpendicular to the flow in a channel of infinite length and immediately a cavity is formed behind the cavitator. This is because of the generation of a gas or vapor cavity between the body and the surrounding liquid due to the change in a high speed flow direction passing the cavitator. Drag force acting on this supercavitating body dictates the thrust requirements for the propulsion system, to maintain a required cavity at the operating speed. Therefore, any reduction in the drag force, by modifying the shape of the cavitator, will lead to decrease this force. This study concentrates on the optimization of two dimensional cavitators in order to decrease drag coefficient for a specified after body length and velocity in a potential flow. To achieve this goal a multi-objective optimization problem is defined to optimize cavitator shapes in supercavitating flow. The so-called NSGA II (Non-dominated Sorting Genetic Algorithm) algorithm is used as an optimization method. Design parameters and constraints are obtained according to supercavitating flow characteristics and cavitator modeling and objective functions are generated using Linear Regression Method. The obtained results are compared with other classic optimization methods, like the weighted sum method, for validation. 相似文献
3.
A. V. Glazunov 《Izvestiya Atmospheric and Oceanic Physics》2006,42(3):282-299
The neutrally stratified boundary layer over a smooth rough surface is consider. The turbulent flow is simulated using a finite-difference eddy-resolving model of the atmospheric boundary layer (ABL). The model includes different turbulence closure schemes and numerical approximations for advection components of the momentum balance equation. We investigate the quality of reproduction of spectral characteristics of the turbulent flow and the model’s capabilities to reproduce the observed profile of mean wind velocity near the rough surface. It is shown that the best result is obtained by coupling a numerical scheme of higher order of accuracy with a mixed closure scheme based on an adaptive estimation of the mixing length for subgrid-scale fluctuations. Here, we are able to reproduce the asymptotics of the fluctuation spectrum of the longitudinal component of wind velocity near the surface and within the boundary layer as well as the logarithmic profile of mean velocity near the surface. 相似文献
4.
In this study, the flow around the pod unit is analysed and the performance characteristics of the propeller on the pod are investigated. The main objective of the present work is to further improve the original numerical method developed before for the prediction of performance of podded propellers and to further validate the earlier developed numerical model with a specific emphasis on the hydrodynamic interaction amongst the propulsor components. While in the earlier numerical method, the axial induced velocities by pod and strut parts were included into the calculations on the propeller disc plane, in the present method the tangential induced velocities on the propeller disc plane are included in the calculations as well. The flow domain around the podded propeller is mainly divided into three parts; the axisymmetric pod part, the strut part and the propeller part. While the pod and strut parts are modelled by a low-order boundary element method (BEM), the propeller is represented by a vortex lattice method (VLM). Coupling of the BEM and the VLM is carried out in an iterative manner to incorporate the effect of the pod on the propeller, and vice versa. The present numerical method is applied to two different podded propellers with zero yaw angles in order to compare the results with those of experimental measurements. The present numerical method is also validated in the case of 15° of yaw angle for a podded propulsor. The effect of pod and strut on the propeller and vice versa are discussed. 相似文献
5.
Bottom ventilated cavitation has been proven as a very effective drag reduction technology for river ships and planning boats. The ability of this technology to withstand the sea wave impact usual for seagoing ships depends on the ship bottom shape and could be enhanced by some active flow control devices. Therefore, there is the need in numerical tools to estimate the effects of bottom changes and to design such devices. The fundamentals of active flow control for the ship bottom ventilated cavitation are considered here on the basis of a special model of cavitating flows. This model takes into account the air compressibility in the cavity, as well as the multi-frequency nature of the incoming flow in wavy seas and of the cavity response on perturbations by incoming flow. The numerical method corresponding to this model was developed and widely manifested with an example of a ship model tested in a towing tank at Froude numbers between 0.4 and 0.7.The impact of waves in head seas and following seas on cavities has been studied in the range of wavelengths from 0.45 to 1.2 of the model (or ship) length. An oscillating cavitator-spoiler was considered as the flow controlling devices in this study. The oscillation magnitude and the phase shift between cavitator oscillation and the incoming waves have been varied to determine the best flow control parameters. The main results of the provided computational analysis include oscillations of cavity surface, of the pressure in cavity and of the moment of hydrodynamic load on the cavitator. The major part of computations has been carried out for the flap oscillating at the frequency coinciding with the wave frequency, but the effect of a frequency shift is also analyzed. 相似文献
6.
Experiments and numerical methods are developed to investigate the water entry of a freefall wedge with a focus on the evolution of the pressure on the impact sides (the side contacting water) and the top side (the dry side on the top of the wedge), evolution of the global hydrodynamic loads, evolution of the air–water interface, and wedge motion. It is found that a typical water entry of a freefall wedge can be divided into slamming, transition, collapse and post-closure stages. A single-fluid numerical model is presented to simulate the first three stages. The results are compared to experiments and good agreements are obtained. A two-fluid BEM is proposed to investigate the influence of the air flow before the closure of the cavity created on the top of the wedge. It is found that for the closure of the 2D cavity, the air flow starts to play an important role just before closure but due to the short duration, the influence of air flow on the body velocity and configuration of the air–water interface is limited. 相似文献
7.
8.
Nonlinear Dynamics and Bifurcations of a Supercavitating Vehicle 总被引：1，自引：0，他引：1
In this effort, a numerical study of the bifurcation behavior of a supercavitating vehicle is conducted. The vehicle model is nonsmooth; this is a result of the planing force acting on the vehicle. With a focus on dive-plane dynamics, bifurcations with respect to a quasi-static variation of the cavitation number are studied. The system is found to exhibit rich and complex nonlinear dynamics including nonsmooth bifurcations such as the grazing bifurcation; smooth bifurcations such as Hopf bifurcations, cyclic fold bifurcations, and period-doubling bifurcations; and aperiodic behaviors such as transient chaotic motions and chaotic crises. The tailslap phenomenon of the supercavitating vehicle is identified as the consequence of a Hopf bifurcation followed by a grazing event. It is shown that the occurrences of these bifurcations can be delayed or triggered earlier by using dynamic linear feedback control laws employing washout filters. 相似文献
9.
A practical, low order and potential-based surface panel method is presented to predict the flow around a three-dimensional rectangular foil section including the effect of boundary layer. The method is based on a boundary-integral formulation, known as the “Morino formulation” and the boundary layer effect is taken into account through a complementary thin boundary layer model. The numerical approach used in the method presents a strongly convergent solution based on the iterative wake roll-up and contraction model including the boundary layer effect. The method is applied to a three-dimensional foil section for which the velocity distribution around the foil was measured using a 2D Laser Doppler Velocimetry system in a large cavitation tunnel. Comparison of the predicted velocity distributions both inside and outside of the boundary layer of the foil as well as the boundary layer shapes obtained from the numerical model show fairly good correlation with the measurements, indicating the robustness and practical worthiness of the proposed method. 相似文献
10.
Very Large Eddy Simulation of Cavitation from Inception to Sheet/Cloud Regimes by A Multiscale Model
《中国海洋工程》2021,(3)
The cavitating flow in different regimes has the intricate flow structure with multiple time and space scales. The present work develops a multiscale model by coupling the volume of fluid(VOF) method and a discrete bubble model(DBM), to simulate the cavitating flow in a convergent-divergent test section. The Schnerr-Sauer cavitation model is used to calculate the mass transfer rate to obtain the macroscale phase structure, and the simplified Rayleigh-Plesset equation is applied to simulate the growing and collapsing of discrete bubbles. An algorithm for bridging between the macroscale cavities and microscale bubbles is also developed to achieve the multiscale simulation. For the flow field, the very large eddy simulation(VLES) approach is applied. Conditions from inception to sheet/cloud cavitation regimes are taken into account and simulations are conducted. Compared with the experimental observations, it is shown that the cavitation inception, bubble clouds formation and glass cavity generation are all well represented, indicating that the proposed VOF-DBM model is a promising approach to accurately and comprehensively reveal the multiscale phase field induced by cavitation. 相似文献
11.
Ventilated cavitation which is acknowledged as an efficient drag reduction technology for underwater vehicle is characterised by the very disparate length and time scales, posing great difficulty in the application of this technology. A multiscale numerical approach which integrates a sub-grid air entrainment model into the two-fluid framework is proposed in this paper to resolve the complex flow field created by ventilated cavity. Simulations have been carried out for the partially ventilated cavity underneath flat plate, with special efforts putting on understanding the gas entrainment at the cavity tail and the bubble dispersion process downstream. The flow parameters including the void fraction, the bubble velocity and the bubble size distributions in and downstream of the ventilated cavity are fully investigated. Comparisons between the numerical results with the experimental data are in satisfactory agreement, demonstrating the potential of the proposed methodology. The ventilation rate effect on the cavity shape and bubbly flow parameters are further investigated, obtaining the law of bubble dispersion and the bubble size evolution. This research not only provide a useful method for the investigation on the multiscale multiphase flow, but also give insight on understanding the combined drag reduction mechanism resulted from large-scale cavity and microbubbles. 相似文献
12.
The successful designs of hulls for ships employing drag reduction by air bottom cavitation have been based on solutions of inverse problems of the theory of ideal incompressible fluid. However, prediction of the drag reduction ratio, the air demand by ventilated cavities and the cavity impact on the hull–propeller interaction is impossible in the framework of this theory because all mentioned characteristics depend on interaction of air cavities with the ship boundary layers. Because the known CFD tools are not fitted to ventilated cavitation at low Froude numbers, an analysis of this interaction requires a novel flow model. This model includes the incompressible air flow in the ventilated cavity, the compressible flow of a water–air mixture in the boundary layer on cavities and downstream of them and the curl-free incompressible outer water flow. The provided 2D computations employing this model allows for explanations of the earlier observed effects and for prediction of the air demand by ventilated cavities. The computed velocity profiles downstream of cavities are in the accordance with the available experimental data. 相似文献
13.
Analysis of body supercavitation in shallow water 总被引：1，自引：0，他引：1
Motion of a cavitating body in shallow waters undergoes a blockage effect. There are influences of the rigid boundary (the sea bottom) and the free boundary (sea surface) in shallow waters. As shown by computation carried out with the ideal fluid theory, the combination of these influences leads to an increase of cavitation number for a cavity of a fixed length and to 3D deformations of the cavity cross-sections, with a swelling of the down cavity part. 相似文献
14.
《Coastal Engineering》1999,36(2):111-146
A numerical model based upon a low Reynolds number turbulence closure is proposed to study Reynolds number variation in reciprocating oscillatory boundary layers. The model is used to compute the boundary layer for flow regimes ranging from smooth laminar to rough turbulent. Criteria for fully developed turbulence are derived for walls of the smooth and rough types. In particular, a new criterion to identify the rough turbulent regime is determined based on the time-averaged turbulence intensity. The reliability of the present model is assessed through comparisons with detailed experimental data collected by other investigators. The model globally improves upon standard high Reynolds number closures. Variation through the wave cycle of the main flow variables (ensemble-averaged velocity, shear stress, turbulent kinetic energy) is remarkably well-predicted for smooth walls. Predictions are satisfactory for rough walls as well. Yet, the turbulence level in the rough turbulent regime is overpredicted in the vicinity of the bed. 相似文献
15.
A numerical method, which combines the boundary element method (BEM) and the volume of the fluid method (VOF method), has been presented to solve wave–structure interactions; the intense wave motion at the proximity of the structure is modeled by the VOF method and the rest of the fluid region is modeled by the BEM. The combined method can considerably reduce the time-consuming VOF domain, and thus practically makes it possible to apply the VOF method for random wave calculations, in which long time computations are usually required to obtain statistically meaningful results, and therefore the use of the single-VOF model often becomes prohibitive in terms of computational time and storage memories. A VOF model CADMAS-SURF, which is based on SMAC scheme and had been constructed by a number of VOF researchers in coastal engineering in Japan, is used in the combined BEM–VOF model. The two-way coupling treatment, which enables us to deal with bidirectional wave propagations, which was originally given for the SOLA-VOF model by Yan et al. (2003a) and later improved by Kim et al. (2007), was modified for the SMAC scheme. The coupling treatments are described in detail in the paper. The validity of the combined BEM–VOF model was investigated by comparing the numerical results with the theoretical results for the propagations of Stokes 5th order waves and random waves. 相似文献
16.
To make a curvilinear motion in the horizontal plane is one of the most contents for realizing the maneuverability of the supercavitating vehicle. It is significant to achieve the controllability and maneuverability of the vehicle in three dimensions both theoretically and practically on research. Models of angle of attack, gravity and inertial force effects on the supercavity in the horizontal curvilinear motion are established, respectively. The supercavity is simulated based on these models in combination with Logvinovich model and the unsteady gas-leakage rate model at the given ventilation rate, and the effect of the inertial force on it is analyzed numerically. Results show that the maximum deviation of the center line of the cross section of supercavity towards the outward normal direction of its trajectory increases as the cavitation number or curvature radius decrease and always occur in the tail because of the increase of inertial effects along the axis of supercavity from the cavitator when other models and flow parameters are constant for the given trajectory curvature. For the variable curvature, the supercavity sheds due to its instability caused by the time-varying angle of attack. The deviation increases along the length of supercavity if the curvature remains the same sign. 相似文献
17.
To make a curvilinear motion in the horizontal plane is one of the most contents for realizing the maneuverability of the supercavitating vehicle. It is significant to achieve the controllability and maneuverability of the vehicle in three dimensions both theoretically and practically on research. Models of angle of attack, gravity and inertial force effects on the supercavity in the horizontal curvilinear motion are established, respectively. The supercavity is simulated based on these models in combination with Logvinovich model and the unsteady gas-leakage rate model at the given ventilation rate, and the effect of the inertial force on it is analyzed numerically. Results show that the maximum deviation of the center line of the cross section of supercavity towards the outward normal direction of its trajectory increases as the cavitation number or curvature radius decrease and always occur in the tail because of the increase of inertial effects along the axis of supercavity from the cavitator when other models and flow parameters are constant for the given trajectory curvature. For the variable curvature, the supercavity sheds due to its instability caused by the time-varying angle of attack. The deviation increases along the length of supercavity if the curvature remains the same sign. 相似文献
18.
A simple relationship has been developed between the wall coordinate y+ and Kolmogorov's length scale using direct numerical simulation (DNS) data for a steady boundary layer. This relationship is then utilized to modify two popular versions of low Reynolds number k–ε model. The modified models are used to analyse a transitional oscillatory boundary layer. A detailed comparison has been made by virtue of velocity profile, turbulent kinetic energy, Reynolds stress and wall shear stress with the available DNS data. It is observed that the low Reynolds number models used in the present study can predict the boundary layer properties in an excellent manner. 相似文献
19.
The present study focuses on the implementation of PISO algorithm to simulate cavitating flows. For simulation of unsteady behaviors of cavitation which have practical applications, the development of unsteady PISO algorithm based on the non-conservative approach is investigated. The effects of mixture compressibility are considered to improve accuracy of simulations. For multi-phase simulation, single-fluid Navier–Stokes equations, along with the volume fraction transport equation, are employed. The bubble dynamics model is utilized to simulate phase change. To prove capabilities of the developed PISO algorithm to simulate cavitating flows, unsteady simulation of cavitation around NACA0015 hydrofoil, a two dimensional flat plat, and a three dimensional circular disk are performed. The frequency of flow, pressure distributions, cavitating vortex shedding, and cavity characteristics are analyzed to discern results accuracy. To investigate accuracy of results, comparisons with available published experimental data are made and good agreement is achieved. 相似文献
20.
Gilbert R.W. Zedler E.A. Grilli S.T. Street R.L. 《Oceanic Engineering, IEEE Journal of》2007,32(1):236-248
In this paper, we report on the use of a numerical wave tank (NWT), based on fully nonlinear potential flow (FNPF) equations, in driving simulations of flow and sediment transport around partially buried obstacles. The suspended sediment transport is modeled in the near-field in a Navier-Stokes (NS) model using an immersed-boundary method and an attached sediment transport simulation module. Turbulence is represented by large eddy simulation (LES). The NWT is based on a higher order boundary element method (BEM), with an explicit second-order time stepping. Hence, only the NWT boundary is discretized. The solution for the velocity potential and its derivatives along the boundary is obtained in the BEM, which subsequently provides a solution at any required internal point within the domain. At initial time, the NS-LES model domain is initialized with the 3-D velocity field provided by the NWT and driven for later time by the pressure gradient field obtained in the NWT. Incident wave fields, as specified in the NWT to drive sediment transport, can be arbitrary. Applications are presented here for single frequency waves, such as produced by a harmonic piston wavemaker in the laboratory, and modulated frequency wave groups. The feasibility of coupling the irrotational flow and NS solutions is demonstrated. 相似文献