• Journal of Korean Society for Atmospheric Environment
• /
• v.9 no.E
• /
• pp.347-357
• /
• 1993

The Lagrangian grid-free numerical method, the discrete vortex method, was applied to solve the Navier-Stokes euqations. This method avoids the introduction of numerical viscosity swamping the real physical viscosity at high Reynolds number, unlike Eulerian method, e.g. finite difference and element methods. The boundary integral equation method for the potential flow solution was included to make the discrete vortex method more feasible for complex geometries. The fast adaptive multipole expansion method was incorporated to reduce the computational time from $O(N^2)$ to O(N) for the computations of vortex-vortex interactions. The test problems were air flow around one circular cylinder and two circular cylinders in tandem with various gaps. The numerical results were in excellent gareement with the experimental and other computational results. The applicabilty of the method was discussed with the indoor and the outdoor air pollution problems, especially the contaminant transport in the recirculation regions.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.6
• /
• pp.744-752
• /
• 2003

A level set method is combined with the volume-of-fluid method so that the coupled method can not only calculate an interfacial curvature accurately but also can achieve mass conservation well. The coupled level set and volume-of-fluid(CLSVOF) method is efficiently implemented by employing an interface reconstruction algorithm which is based on the explicit relationship between the interface configuration and the fluid volume function. The CLSVOF method is applied for numerical simulation of droplet impact on solid surfaces with variable contact angles. The numerical results are found to preserve mass conservation and to be in good agreement with the data reported in the literature. Also, the present method proved to be applicable to the complex phenomena such as breakup and rebound of a droplet.

• Geomechanics and Engineering
• /
• v.22 no.1
• /
• pp.97-108
• /
• 2020

The time-consuming and less objectivity are the main problems of conventional micromechanical parameters calibration method of Particle Flow Code simulations. Thus this study aims to address these two limitation of the conventional "trial-and-error" method. A new calibration method for the linear parallel bond model (CM-LPBM) is proposed. First, numerical simulations are conducted based on the results of the uniaxial compression tests on limestone. The macroscopic response of the numerical model agrees well with the results of the uniaxial compression tests. To reduce the number of the independent micromechanical parameters, numerical simulations are then carried out. Based on the results of the orthogonal experiments and the multi-factor variance analysis, main micromechanical parameters affecting the macro parameters of rocks are proposed. The macro-micro parameter functions are ultimately established using multiple linear regression, and the iteration correction formulas of the micromechanical parameters are obtained. To further verify the validity of the proposed method, a case study is carried out. The error between the macro mechanical response and the numerical results is less than 5%. Hence the calibration method, i.e., the CM-LPBM, is reliable for obtaining the micromechanical parameters quickly and accurately, providing reference for the calibration of micromechanical parameters.

• Journal of applied mathematics & informatics
• /
• v.29 no.1_2
• /
• pp.49-53
• /
• 2011

In this paper, we adopt a numerical method to establish the smallest set to contain all Ger$\v{s}$gorin discs of a given complex matrix and its some similar matrices. With the smallest set, a new estimation for all eigenvalues of the matrix is obtained.

• Journal of Navigation and Port Research
• /
• v.32 no.7
• /
• pp.543-552
• /
• 2008

The diffraction of waves by three bottom fixed vertical circular cylinders is investigated by using the boundary element method. This method has been successfully applied to the isolated vertical circular cylinder and now is used to study the interaction between waves and multiple vertical cylinders. In this paper, a numerical analysis by the boundary element method is developed by the linear potential theory. The numerical analysis by the boundary element method is based on Green's second theorem and introduced to an integral equation for the fluid velocity potential around the vertical circular cylinders. To verify this method, the results obtained in present study are compared with the results computed by the multiple scattering method. The results of the comparisons show strong agreement. Also in this paper, several numerical examples are given to illustrate the effects of various parameters on the wave exciting force such are the separation distance, the wave number and the incident wave angle. This numerical computation method might be used broadly for the design of various offshore structures to be constructed in the future.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1998.04a
• /
• pp.568-574
• /
• 1998

A systematic procedure for designing a direct-radiator-type loudspeaker has been developed, based on a numerical vibro-acoustic analysis and the Taguchi method. The finite-element model of the speaker cone has been used to calculate the vibration response of the cone excited by the voice coil. The vibration response of the speaker cone has been used as a boundary condition for the acoustic analysis, and the acoustic frequency characteristics of the loudspeaker have been calculated by the boundary element method. The numerical model has been confirmed by comparing the numerical results with experimental ones obtained in an anechoic chamber. Some design parameters contributing dominantly to the acoustic characteristics have been selected by using the Taguchi method, and the variations of the acoustic characteristics due to the changes of the parameter values have been examined using the numerical model.

• 한국방재학회:학술대회논문집
• /
• 2007.02a
• /
• pp.181-185
• /
• 2007

A free-surface correction(FSC) method is presented to solve the 3-D shallow water equations. Using the mode splitting process, FSC method can simulate shallow water flows under the hydrostatic assumption. For the hydrostatic pressure calculation, the momentum equations are firstly discretized using a semi-implicit scheme over the vertical direction leading to the tri-diagonal matrix systems. A semi-implicit scheme has been adopted to reduce the numerical instability caused by relatively small vertical length scale compare to horizontal one. and, as the free surface correction step the final horizontal velocity fields are corrected after the final surface elevations are obtained. Finally, the vertical final velocity fields can be calculated from the continuity equation. The numerical model is applied to the calculation of the simulation of flow fields in a rectangular open channel with the tidal influence. The comparisons with the analytical solutions show overall good agreements between the numerical results and analytical solutions.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.2
• /
• pp.206-218
• /
• 2014

The wave attenuation by floating breakwaters in high amplitude waves, which can lead to wave overtopping and breaking, is examined by numerical simulations. The governing equations, the Navier-Stokes equations and the continuity equation, are calculated in a fixed Cartesian grid system. The body boundaries are defined by the line segment connecting the points where the grid line and body surface meet. No-slip and divergence free conditions are satisfied at the body boundary cell. The nonlinear waves near the moving body is defined using the modified marker-density method. To verify the present numerical method, vortex induced vibration on an elastically mounted cylinder and free roll decay are numerically simulated and the results are compared with those reported in the literature. Using the present numerical method, the wave attenuations by three kinds of floating breakwaters are simulated numerically in a regular wave to compare the performance.

• Journal of information and communication convergence engineering
• /
• v.8 no.3
• /
• pp.312-316
• /
• 2010

Conformal mapping has been a familiar tool of science and engineering for generations. The methods of numerical mapping are usually classified into those which construct the map from standard domain such as the unit disk onto the 'problem domain', and those which construct the map in the reverse direction. We treat numerical conformal mapping from the unit disk onto the Jordan regions as the problem domain in this paper. The traditional standard methods of this type are based on Theodorsen integral equation. Wegmann's method is well known as a Newton-like efficient one for solving Theodorsen equation. An improved method for convergence by applying low frequency pass filter to the Wegmann's method was proposed. In this paper we propose an effective algorithm for numerical conformal mapping based on the improved method. This algorithm is able to determine the discrete numbers and initial values automatically in accordance with the given region and the required accuracy. This results come from analyzing the shape of given domain as seen in the Fourier Transform.

• Journal of computational fluids engineering
• /
• v.9 no.2
• /
• pp.43-51
• /
• 2004

Interaction between fluid flow and thermal radiation has received considerable attention due to its numerous applications in engineering field. In this case the thermofluid properties of radiating fluid vary with the variation of temperature field caused by absorption and emission of radiant heat. To analyze the radiation heat transfer in radiating fluid, the simultaneous solution of the radiative transfer equation (RTE) and the fluid dynamics equations is required. This means that the numerical procedure used for the RTE must be computationally efficient to permit its inclusion in the other submodels, and must be compatible with the other transport equations. The finite volume method (FVM) and the discrete ordinates method (DOM) are usually employed to simulate radiation problems in generalized coordinates. These two representative methods are examined and compared, especially in view of the numerical integration of the radiation intensity over solid angle. The FVM shows better accuracy than the DOM owing to less constraints of the selection of control angle.