• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.2
• /
• pp.159-173
• /
• 1991

This paper describes a potential-based panel method formulated for the analysis of a super-cavitating two-dimensional hydrofoil. The method employs normal dipoles and sources distributed on the foil and cavity surfaces to represent the potential flow around the cavitating hydrofoil. The kinematic boundary condition on the wetted portion of the foil surface is satisfied by requiring that the total potential vanish in the fictitious inner flow region of the foil, and the dynamic boundary condition on the cavity surface is satisfied by requiring thats the potential vary linearly, i.e., the tangential velocity be constant. Green's theorem then results in a potential-based integral equation rather than the usual velocity-based formulation of Hess & Smith type. With the singularities distributed on the exact hydrofoil surface, the pressure distributions are predicted with improved accuracy compared to those of the linearized lilting surface theory, especially near the leading edge. The theory then predicts the cavity shape and cavitation number for an assumed cavity length. To improve the accuracy, the sources and dipoles on the cavity surface are moved to the newly computed cavity surface, where the boundary conditions are satisfied again. This iteration process is repeated until the results are converged. Characteristics of iteration and discretization of the present numerical method are much faster and more stable than the existing nonlinear theories. The theory shows good correlations with the existing theories and experimental results for the super-cavitating flow. In the region of small angles of attack, the present prediction shows and excellent comparison with the Geurst's linear theory. For the long cavity, the method recovers the trends of the Wu's nonlinear theory. In the intermediate regions of the short super-cavitation, the method compares very well with the experimental results of Parkin and also those of Silberman.

• Journal of computational fluids engineering
• /
• v.13 no.2
• /
• pp.27-41
• /
• 2008

An unstructured hybrid mesh flow solver has been developed for the simulation of three-dimensional steady and unsteady incompressible flow fields. The incompressible Navier-Stokes equations with an artificial compressibility method were discretized by using a node-based finite-volume method. For the unsteady time-accurate computation, a dual-time stepping method was adopted to satisfy a divergence-free flow field at each physical time step. An implicit time integration method with local time stepping was implemented to accelerate the convergence in the pseudo-time sub-iteration procedure. The one-equation Spalart-Allmaras turbulence model has been adopted to solve high-Reynolds number flow fields. The flow solver was parallelized to minimize the CPU time and to overcome the computational overhead. This method has been applied to calculate steady and unsteady flow fields around submarine configurations and a 3-D infinite cylinder. Validations were made by comparing the predicted results with those of experiments or other numerical results. It was demonstrated that the present method is efficient and robust for the prediction of steady and unsteady incompressible flow fields.

• The Journal of the Korea institute of electronic communication sciences
• /
• v.10 no.4
• /
• pp.507-512
• /
• 2015

In the wireless localization systems using range measurements, iteration method-based approximated solutions have been used. Also, linear closed-form solutions have been investigated in the light of local minimum problem and computational load. However, each closed-form solution has individual error factors that cause usage limit of the solutions. In this paper, a fusion method integrating two representative closed-form solutions is presented. The presented method cancels the error factors of each solution out. Weights for integrating the standalone solutions are determined using the error factors-based fuzzy method. The performance of the proposed method is verified using some simulation results.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.11a
• /
• pp.1022-1027
• /
• 2007

This paper presents a method to analyze the unbalance response of a high speed polygon mirror scanner motor supported by sintered metal bearing and flexible structures by using the finite element method and the mode superposition method considering the asymmetry of the gyroscopic effect and sintered metal bearing. The eigenvalues and eigenvectors are calculated by solving the eigenvalue problem and the adjoint eigenvalue problem by using the restarted Arnoldi iteration method. The decoupled equations of motion can be obtained from global finite element motion equations by using the orthogonal relation between the right eigenvectors and left eigenvectors. The decoupled equations of motion are used to analyze the unbalance response of a high speed polygon mirror scanner motor. The validity of the proposed method is verified by comparing the simulated unbalance response with the experimental results.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.26 no.1
• /
• pp.122-128
• /
• 1989

In this paper we propose a new method for solving a set of nonlinear algebraic equations encountered in the DC and transient analyses of electronic circuits. This method will be called Quadratic Newton-Raphson Method (QNRM), since it is based on the Newton-Raphson Method (NRM) but effectively takes into accoujnt the second order derivative terms in the Taylor series expansion of the nonlinear algebraic equations. The second order terms are approximated by linear terms using a carefully estimated solution at each iteration. Preliminary simulation results show that the QNRM saves the overall computational time significantly in the DC and transient analysis, compared with the conventional NRM.

• Journal of Korean Society of Steel Construction
• /
• v.8 no.3 s.28
• /
• pp.21-34
• /
• 1996

A numerical method is presented for computation of eigenvector derivatives used an iterative procedure with guaranteed convergence. An approach for treating the singularity in calculating the eigenvector derivatives is presented, in which a shift in each eigenvalue is introduced to avoid the singularity. If the shift is selected properly, the proposed method can give very satisfactory results after only one iteration. A criterion for choosing an adequate shift, dependent on computer hardware is suggested ; it is directly dependent on the eigenvalue magnitudes and the number of bits per numeral of the computer. Another merit of this method is that eigenvector derivatives with repeated eigenvalues can be easily obtained if the new eigenvectors are calculated. These new eigenvectors lie "adjacent" to the m (number of repeated eigenvalues) distinct eigenvectors, which appear when the design parameter varies. As an example to demonstrate the efficiency of the proposed method in the case of distinct eigenvalues, a cantilever plate is considered. The results are compared with those of Nelson's method which can find the exact eigenvector derivatives. For the case of repeated eigenvalues, a cantilever beam is considered. The results are compared with those of Dailey's method which also can find the exact eigenvector derivatives. The design parameter of the cantilever plate is its thickness, and that of the cantilever beam its height.

• Structural Engineering and Mechanics
• /
• v.46 no.5
• /
• pp.735-753
• /
• 2013

The wind load is always the dominant load of cooling tower due to its large size, complex geometry and thin-wall structure. At present, when computing the wind-induced response of the large-scale cooling tower, the wind pressure distribution is obtained based on code regulations, wind tunnel test or computational fluid dynamic (CFD) analysis, and then is imposed on the tower structure. However, such method fails to consider the change of the wind load with the deformation of cooling tower, which may result in error of the wind load. In this paper, the analysis of the large cooling tower based on the iterative method for wind pressure is studied, in which the advantages of CFD and finite element method (FEM) are combined in order to improve the accuracy. The comparative study of the results obtained from the code regulations and iterative method is conducted. The results show that with the increase of the mean wind speed, the difference between the methods becomes bigger. On the other hand, based on the design of experiment (DOE), an approximate model is built for the optimal design of the large-scale cooling tower by a two-level optimization strategy, which makes use of code-based design method and the proposed iterative method. The results of the numerical example demonstrate the feasibility and efficiency of the proposed method.

• Structural Engineering and Mechanics
• /
• v.30 no.2
• /
• pp.165-176
• /
• 2008

In this paper, a new iterative method for solving vehicle-bridge interaction problems is proposed. Iterative methods have advantages over the non-iterative methods in that it is not necessary to update the system matrix for a given wheel location, and the method can be applied for a new type of car or bridge with few or no modifications. In the proposed method, the necessity of system matrices update is eliminated using the equivalent interaction force acting on the bridge, which is obtained iteratively. Ballast stiffness is included in the interaction forces and the geometric compatibility at the contact points are used as convergence criteria. The bridge is considered as an elastic Bernoulli-Euler beam with surface irregularity and ballast stiffness. The moving vehicle is modeled as a multi-axle mass-spring-damper system having many degrees of freedom depending on the number of axles. The pitching effect, which is the interaction effect between the rear and front wheels when a vehicle begins to enter or leave the bridge, is also considered in the formulation including extended ground boundaries having surface irregularity and ballast stiffness. The applicability of the proposed method is illustrated in the numerical studies.

• Computational Structural Engineering
• /
• v.11 no.1
• /
• pp.179-190
• /
• 1998

A geometrically nonlinear finite element formulation of spatial cable networks is presented using two cable elements. Firstly, derivation procedures of tangent stiffness and mass matrices for the space truss element and the elastic catenary cable element are summarized. The load incremental method based on Newton-Raphson iteration method and the dynamic relaxation method are presented in order to determine the initial static state of cable nets subjected to self-weights and support motions. Furthermore, static non-linear analysis of cable structures under additional live loads are performed based on the initial configuration. Challenging example problems are presented and discussed in order to demonstrate the feasibility of the present finite element method and investigate static nonlinear behaviors of cable nets.

• Computational Structural Engineering
• /
• v.2 no.3
• /
• pp.97-103
• /
• 1989

In this study, a usefulness of the iterative solution procedures is reviewed for the elasto-plastic large deflection analysis of imperfect plates by finite element method. Three typical solution techniques such as simple incremental(SI) method, Newton-Raphson(NR) method and modified Newton-Raphson (mNR) method are compared. It is concluded that for thin plates which are given rise to the large deflection, iteration for the convergence of the unbalance force should be performed and in this case mNR method is more useful than NR method since the computing time of the former becomes to be a half of the latter, in which the accuracy of the result remains same. For thick plates or thin plates with large initial deflection, however, the use of SI method is quite better since the unbalance force may be negligible.