• Journal of Korean Society of Steel Construction
• /
• v.12 no.2 s.45
• /
• pp.221-230
• /
• 2000

Checking the stability of anisotropic plates with free edges, it is impossible that buckling loads and modes are found via existing classical methods about various loads and boundary conditions. For solving this problems. finite difference method(FDM) is used to analyze the buckling behaviors for arbitrary boundary conditions. Using FDM, it is difficult to treat the fictitious points on free edges. So, this paper analyzes buckling behaviors of analytic models with one edge free and the other edges clamped and with opposite two edges free and other two edges clamped. The various buckling loads and mode characteristics through numerical results are given for buckling behaviors of anisotropic plates on free edges.

• Journal of Korean Society of Steel Construction
• /
• v.12 no.3 s.46
• /
• pp.303-315
• /
• 2000

The main object of this study is to analyze mechanical behaviors as anisotropic three-dimensional body under various static loads. This paper presents the applicability of the finite difference method to three dimensional problem of anisotropic body. The finite difference method as applied here is generalized to anisotropic three-dimensional problem of elastic body where the governing differential equations of equilibrium of such bodies are expressed in terms of the displacement u, v, and w in the coordinates axes x, y and z, care being taken to modify the finite difference expressions to satisfy the appropriate boundary conditions. By adopting a new three dimensional finite difference modelling including elimination of pivotal difference points in the case of free boundary condition, the three dimensional problem of anisotropic body was successfully completed. Several numerical results show quick convergence and numerical validity of finite difference technique in three dimensional problem.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.3
• /
• pp.80-89
• /
• 1992

This paper deals with the open boundary problems, and two numerical schemes are used for the implementation of open boundary condition. One is to add the artificial damping term to dynamic free-surface boundary condition. Determination of suitable damping coefficient and the damping cone is the most important in this scheme. The other scheme is a modified Orlanski's method. This will be useful for the problems with unidirectional waves. A few typical free-surface wave problems are modeled for the numerical test. Method of solution is fundamental source-distribution method and the fully nonlinear boundary conditions are applied. The computed results are compared with those of others for the proof of practicality of these schemes.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.12
• /
• pp.1165-1172
• /
• 2009

The laminar boundary layer along a streamwise corner formed by two flat plates intersecting at right angle is measured by using Particle Image Velocimetry(PIV) technique. The free stream velocity ranges from 2.96m/s to 3.0m/s. The angle of incidence of the corner is set to 1.2 degree providing slightly favourable pressure gradient to ensure a laminar flow in the corner region. A round shape leading edge is used and the length of the model is about 1000mm. In the bisector plane, the measurement data show separation type velocity profiles having an inflection point which is a typical characteristic of laminar corner boundary layers. As the distance away from the bisector plane increases, velocity profiles are found to change into the Blasius profile. The change completes around half length of the boundary layer thickness in the bisector plane away from the bisector plane along the plate. In the bisector plane, the growth characteristic of the boundary layer thickness and the approximate similarity of velocity profiles are confirmed from the measurement data.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.98-111
• /
• 1993

Two-dimensional nonlinear free-surface wave problems are analyzed with consideration of viscosity. Navier-Stokes equation and continuity equation are solved by the application of Finite Analytic Method, and MAC scheme is used far the treatment of free surface. Surface tension effect is also considered and laminar flow is assumed. The free-surface waves in shallow water, the flows around a vortex-pair with free surface and the wave ahead of a rectangular body are simulated to test the present numerical scheme. In the shallow water problem, viscous effect due to the friction on the bottom is observed. In the second problem, the approach of a vortex-pair to the free surface is simulated to examine the interaction of vortex-pair with the free surface. In the third problem, the wave ahead of a semi-infinite floating body is simulated.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.6
• /
• pp.1445-1451
• /
• 1994

To calculate the wave pressure acting on coastal structures under the design wave condition, it is often necessary to numerically reproduce the big standing wave profiles close to wave breaking condition. For this, the governing equation and all nonlinear terms occurring in boundary conditions should be effectively considered in the numerical wave profile. In particular, the velocity square term in the free surface boundary condition is very important. A boundary element method is applied here to calculate the standing wave profile with the velocity square term fully treated by Newton iterative method. In order to check the validity of the method, the numerical wave profiles are compared to ones calculated by the perturbation method, the Fourier approximation method and the hydraulic experiment.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.8 no.1
• /
• pp.37-43
• /
• 1996

Breaking waves were generated in a 2-D flume. A piston-type wavemaker was operated in accordance with signals which consist of elementary harmonics with appropriate phase differences. These phase differences were estimated by using a linear wave theory so that wave crests were to be concentrated at the same position. The stroke of wavemaker was controlled to create plunging-type breaking waves. The signal with small amplitude could not generate breaking waves. In the case of moderate amplitudes, various breaking waves could be obtained. Most of breaking waves were spilling type. Only when the wavemaker was operated with appropriate amplitude, plunging-type breaking waves were generated. The parameters of breaking waves are the wave steepness and the frequency bandwidth. If the central frequency was low, breaking waves were not generated. Based on experimental data, we found that the wave height of breaking inception was H = 0.0113 gT$^2$. We also made computations by using a mixed Euler-Lagrangian scheme under the assumption of potential flow. The numerical results show good agreements with tank measurements.

• Journal of Korea Water Resources Association
• /
• v.52 no.10
• /
• pp.697-706
• /
• 2019

A multiphase flow modeling approach equipped with a hybrid turbulence modeling method is applied to compute the gravity currents in a rectangular channel. The present multiphase solver considers the dense fluid, the less-dense ambient fluid and the air above free surface as three phases with separate flow equations for each phase. The turbulent effect is simulated by the IDDES (improved delayed detach eddy simulation), a hybrid RANS/LES, approach which resolves the turbulent flow away from the wall in the LES mode and models the near wall flow in RANS mode on moderately fine computational meshes. The numerical results show that the present model can successfully reproduce the gravity currents in terms of the propagation speed of the current heads and the emergence of large-scale Kelvin-Helmholtz type interfacial billows and their three dimensional break down into smaller turbulent structures, even on the relatively coarse mesh for wall-modeled RANS computation with low-Reynolds number turbulence model. The present solutions reveal that the modeling approach can capture the large-scale three dimensional behaviors of gravity current head accompanied by the lobe-and-cleft instability at affordable computational resources, which is comparable to the LES results obtained on much fine meshes. It demonstrates that the multiphase modeling method using the hybrid turbulence model can be a promising engineering solver for predicting the physical behaviors of gravity currents in natural environmental configurations.

• Computational Structural Engineering
• /
• v.8 no.3
• /
• pp.51-57
• /
• 1995

A simplified stress analysis of perforated plates was carried out using homogenization technique. Homogenization technique, which introduced miroscale expansion in the standard finite element method, reconstructed the plate with regularly placed holes into a set of macroscale and microscale models. The microscale model helped compute homogenized material constants of the unit cell, which were used to compute macroscale displacements in the macroscale model. Also it was possible to compute the stress field of the plate using the microscale model. It was found that reasonable equivalent material constants were computed and that the required degrees of freedom was drastically reduced when homogenization technique was employed in the stress analyses. The microscale modeling in the homogenization technique provided a useful concept of pre- and post-processing in the stress analysis of perforated plates.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.9 no.2
• /
• pp.97-104
• /
• 2005

A simple but efficient grid generation technique by using the modified compressible form of stream function has been formulated. Transformation of a physical plane to a streamline plane, the Von Mises Transformation, has been widely used to solve the differential equations governing flow phenomena, however, limitation arises in low velocity region of boundary layer, mixing layer and wake region where the relatively large grid spacing is inevitable. Modified Von Mises Transformation with simple mathematical adjustment for the stream function is suggested and applied to solve the confined coaxial turbulent jet mixing with simple $\kappa-\epsilon$ turbulence model. Comparison with several experimental data of axial mean velocity, turbulent kinetic energy, and Reynolds shear stress distribution shows quite good agreement in the mixing layer except in the centerline where the turbulent kinetic energy distributions were somewhat under estimated. This formulation is strongly suggested to be utilized specially for free turbulent mixing layers in axisymmetric flow conditions such as the investigation of mixing behavior, jet noise production and reduction for Turbofan engines.