Based on the potential flow assumption, this paper outlines the formulation of the boundary-value problem governing the steady free surface flow around a trimaran advancing at a constant forward. A Rankine panel metho...Based on the potential flow assumption, this paper outlines the formulation of the boundary-value problem governing the steady free surface flow around a trimaran advancing at a constant forward. A Rankine panel method was developed for the solution of the problem. The side-hulls were treated as lifting bodies and circulation was introduced into the flow field to implement the asymmetrical flow around each side hull. The linearized free surface condition was used. The appropriate treatment for transom-stern ship was also discussed. Additionally, a test was performed to quantify the resistance and especially the interference between the main and side hulls. The numerical results are in good argeement with experimental data.展开更多
Using the axial symmetry results of marker and cell (MAC) method as initial value in this paper, two numerical calculating methods are presented for the late wavemaking response induced by explosion in harbour. One of...Using the axial symmetry results of marker and cell (MAC) method as initial value in this paper, two numerical calculating methods are presented for the late wavemaking response induced by explosion in harbour. One of the methods is the superposition method of the vibration mode based on fluid slosh in container. Another one is the joining method of the MAC results with the shallow wave theory calculation in time domain. As a practical example, it is conducted to the numerical calculation about 1000 ton TNT equivalent explosion within touch of water surface. The results show that it can be rationally described with the methods to the wavemaking progress and character. The numerical results are identical with the observed scene on the spot experiment. The methods are simple and applicable in the engineering design.展开更多
Panel methods for the calculation of wavemaking resistance result in a linear equation system for the unknown singularities.The coefficient matrix is full but not well conditioned.In this paper an incomplete LU decomp...Panel methods for the calculation of wavemaking resistance result in a linear equation system for the unknown singularities.The coefficient matrix is full but not well conditioned.In this paper an incomplete LU decomposition (ILU) method and a combined multigrid ILU method are used to solve the linear system.Systematic computations using the ILU method have shown that the CPU time can be reduced to 30% to 40% of that using an incomplete Gaussian elimination method. In the proposed multigrid ILU method an averaged restriction and a piecewise constant prolongation are used.The construction of the coefficient matrix at coarse levels is based on geometrical considerations.It turns out that the condition of the relative consistency is fulfilled.Comparison computations have shown that nearly the same results were obtained.However,due to additional CPU time needed for the execution of the matrix vector products in the restriction and the prolongation proceses of the multigrid method,a further reduction of the total CPU time could not be reailized.展开更多
Under the' assumption of linearization of the free-surface condition, making use of Green's function method and the convolution theorem, analytic solutions of perturbation velocity potentials which correspond ...Under the' assumption of linearization of the free-surface condition, making use of Green's function method and the convolution theorem, analytic solutions of perturbation velocity potentials which correspond to three dimensional unsteady thickness problem and lifting problem caused respectively by arbitrary motions of a body and a hydrofoil beneath the water surface can be achieved in the closed form, In general, the whole perturbation velocity potential consists of three terms, namely φ=φ1+φ2+φ3 , where φ1 denotes the induced velocity potential of the surface singularity distribution in an unbounded fluid, φ2 denotes its mirror image and φ3 denotes that of wave formation which includes the memory effect of the action of the singularity distribution. Utilizing the polynomial expansion of sin[(t-τ)] , the similarity between φ2 and φ3 is discovered and thus a simpler differential relation between them is obtained. Applying this relation, the amount of work in calculation of φ3 which is the most time-consuming one will be reduced significantly. It is favorable not only for dealing with unsteady wave- making problems but also for solving the steady ones in virtue of evading a major difficulty which has to be encountered during the evaluation of an improper inte- gral containing a singularity in the Green's function. The limitation of this new technique turns out to be its slower convergence as the Froude number is lower.展开更多
文摘Based on the potential flow assumption, this paper outlines the formulation of the boundary-value problem governing the steady free surface flow around a trimaran advancing at a constant forward. A Rankine panel method was developed for the solution of the problem. The side-hulls were treated as lifting bodies and circulation was introduced into the flow field to implement the asymmetrical flow around each side hull. The linearized free surface condition was used. The appropriate treatment for transom-stern ship was also discussed. Additionally, a test was performed to quantify the resistance and especially the interference between the main and side hulls. The numerical results are in good argeement with experimental data.
基金Subsidized subject financially supported by the National Natural Science Foundation of China
文摘Using the axial symmetry results of marker and cell (MAC) method as initial value in this paper, two numerical calculating methods are presented for the late wavemaking response induced by explosion in harbour. One of the methods is the superposition method of the vibration mode based on fluid slosh in container. Another one is the joining method of the MAC results with the shallow wave theory calculation in time domain. As a practical example, it is conducted to the numerical calculation about 1000 ton TNT equivalent explosion within touch of water surface. The results show that it can be rationally described with the methods to the wavemaking progress and character. The numerical results are identical with the observed scene on the spot experiment. The methods are simple and applicable in the engineering design.
文摘Panel methods for the calculation of wavemaking resistance result in a linear equation system for the unknown singularities.The coefficient matrix is full but not well conditioned.In this paper an incomplete LU decomposition (ILU) method and a combined multigrid ILU method are used to solve the linear system.Systematic computations using the ILU method have shown that the CPU time can be reduced to 30% to 40% of that using an incomplete Gaussian elimination method. In the proposed multigrid ILU method an averaged restriction and a piecewise constant prolongation are used.The construction of the coefficient matrix at coarse levels is based on geometrical considerations.It turns out that the condition of the relative consistency is fulfilled.Comparison computations have shown that nearly the same results were obtained.However,due to additional CPU time needed for the execution of the matrix vector products in the restriction and the prolongation proceses of the multigrid method,a further reduction of the total CPU time could not be reailized.
文摘Under the' assumption of linearization of the free-surface condition, making use of Green's function method and the convolution theorem, analytic solutions of perturbation velocity potentials which correspond to three dimensional unsteady thickness problem and lifting problem caused respectively by arbitrary motions of a body and a hydrofoil beneath the water surface can be achieved in the closed form, In general, the whole perturbation velocity potential consists of three terms, namely φ=φ1+φ2+φ3 , where φ1 denotes the induced velocity potential of the surface singularity distribution in an unbounded fluid, φ2 denotes its mirror image and φ3 denotes that of wave formation which includes the memory effect of the action of the singularity distribution. Utilizing the polynomial expansion of sin[(t-τ)] , the similarity between φ2 and φ3 is discovered and thus a simpler differential relation between them is obtained. Applying this relation, the amount of work in calculation of φ3 which is the most time-consuming one will be reduced significantly. It is favorable not only for dealing with unsteady wave- making problems but also for solving the steady ones in virtue of evading a major difficulty which has to be encountered during the evaluation of an improper inte- gral containing a singularity in the Green's function. The limitation of this new technique turns out to be its slower convergence as the Froude number is lower.