• Wind and Structures
• /
• v.34 no.1
• /
• pp.59-72
• /
• 2022

In this work the numerical results of the flow around a 5:1 rectangular cylinder at Reynolds numbers 3 000 and 40 000, zero angle of attack and smooth incoming flow condition are presented. Implicit Large Eddy Simulations (ILES) have been performed with a high-order accurate spatial scheme and an implicit high-order accurate time integration method. The spatial approximation is based on a discontinuous Galerkin (dG) method, while the time integration exploits a linearly-implicit Rosenbrock-type Runge-Kutta scheme. The aim of this work is to show the feasibility of high-fidelity flow simulations with a moderate number of DOFs and large time step sizes. Moreover, the effect of different parameters, i.e., dimension of the computational domain, mesh type, grid resolution, boundary conditions, time step size and polynomial approximation, on the results accuracy is investigated. Our best dG result at Re=3 000 perfectly agrees with a reference DNS obtained using Nek5000 and about 40 times more degrees of freedom. The Re=40 000 computations, which are strongly under-resolved, show a reasonable correspondence with the experimental data of Mannini et al. (2017) and the LES of Zhang and Xu (2020).

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.12
• /
• pp.3080-3087
• /
• 2000

We obtain the mathematical model of the hard disk drive actuator system the system response data of the finite element analysis or experimental results. The model is base on the Rayleigh-Ritz method to approximate the dynamic response of the actuator system. The basic idea is to use the curve-fit technique to obtain the approximation coefficients. It allows the dynamic analysis of the actuator system without resort to the repetitive finite element modeling work. Even though the dynamic characteristics of the system of the system are affected somewhat by the structural modification and the change of the material properties, we can use the modified size and dynamic properties of the actuator system in the mathematical model to some extent. In this study, we express the mathematical model of the simplified rectangular plate first and then proceed to the actual hard disk drive actuator system.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.4
• /
• pp.507-514
• /
• 2003

Most of CAE analyses for injection molding have been based on the Mele Shaw's approximation: two-dimensional flow analysis. in some cases, that approximation causes significant errors due to loss of the geometrical information as well as simplification of the flow characteristics in the thickness direction. Although injection molding analysis software using three-dimensional solid elements has been developed recently, such as Moldflow Flow3D, it does not contain a deformation analysis function yet. The present work covers three-dimensional deformation analysis or injection molded plastic parts using solid elements. A numerical scheme for deformation analysis has bun proposed from the results of injection molding analysis using Moldflow Flow3D. The accuracy of the proposed approach has been verified through a numerical analysis of rectangular plates with various thicknesses in comparison with the classical shell-based approach. In addition, the reliability of the approach has also been proved through an industrial example. an optical plastic lens, in comparison of real experiments.

• Nuclear Engineering and Technology
• /
• v.48 no.1
• /
• pp.43-54
• /
• 2016

In the present paper, development of the three-dimensional (3D) computational code based on Galerkin finite element method (GFEM) for solving the multigroup forward/adjoint diffusion equation in both rectangular and hexagonal geometries is reported. Linear approximation of shape functions in the GFEM with unstructured tetrahedron elements is used in the calculation. Both criticality and fixed source calculations may be performed using the developed GFEM-3D computational code. An acceptable level of accuracy at a low computational cost is the main advantage of applying the unstructured tetrahedron elements. The unstructured tetrahedron elements generated with Gambit software are used in the GFEM-3D computational code through a developed interface. The forward/adjoint multiplication factor, forward/adjoint flux distribution, and power distribution in the reactor core are calculated using the power iteration method. Criticality calculations are benchmarked against the valid solution of the neutron diffusion equation for International Atomic Energy Agency (IAEA)-3D and Water-Water Energetic Reactor (VVER)-1000 reactor cores. In addition, validation of the calculations against the $P_1$ approximation of the transport theory is investigated in relation to the liquid metal fast breeder reactor benchmark problem. The neutron fixed source calculations are benchmarked through a comparison with the results obtained from similar computational codes. Finally, an analysis of the sensitivity of calculations to the number of elements is performed.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2004.05b
• /
• pp.557-562
• /
• 2004

A surface-piercing barrier model is presented for understanding morphological development in the sheltered region and investigating the main factors causing the severe accumulation. Surface-piercing structures like vertical barriers, surface docks and floating breakwaters are recently favored from the point of view of a marine scenario since they do not in general partition the natural sea. The numerical solutions are compared with experimental data on wave profiles and morphological change rates within a rectangular harbor of a constant depth protected by surface-piercing thin breakwaters as a simplified problem. Our numerical study involves several modules: 1) wave dynamics analyzed by a plane-wave approximation, 2) suspended sediment transport combined with sediment erosion-deposition model, and 3) concurrent morphological changes. Scattering waves are solved by using a plane wave method without inclusion of evanescent modes. Evanescent modes are only considered in predicting the reflection ratio against the vertical barrier and energy losses due to vortex shedding from the lower edge of plate are taken into account. A new relationship to relate the near-bed concentration to the depth-mean concentration is presented by analyzing the vertical structure of concentration. The numerical solutions were also compared with experimental data on morphological changes within a rectangular harbor of constant water depth. Through the numerical experiments, the vortex-induced flow appears to be not ignorable in predicting the morphological changes although the immersion depth of a plate is not deep.

• Ocean Systems Engineering
• /
• v.2 no.3
• /
• pp.229-238
• /
• 2012

We discuss the bi - stability that is possibly exhibited by a liquid free surface in a parametrically - driven two-dimensional (2D) rectangular tank with finite liquid depth. Following the method of adaptive mode ordering, assuming two dominant modes and retaining polynomial nonlinearities up to third-order, a nonlinear finite-dimensional nonlinear modal system approximation is obtained. A "continuation method" of nonlinear dynamics is then used in order to elicit efficiently the instability boundary in parameters' space and to predict how steady surface elevation changes as the frequency and/or the amplitude of excitation are varied. Results are compared against those of the linear version of the system (that is a Mathieu-type model) and furthermore, against an intermediate model also derived with formal mode ordering, that is based on a second - order ordinary differential equation having nonlinearities due to products of elevation with elevation velocity or acceleration. The investigation verifies that, in parameters space, there must be a region, inside the quiescent region, where liquid surface instability is exhibited. There, behaviour depends on initial conditions and a wave form would be realised only if the free surface was substantially disturbed initially.

• International Journal of Aeronautical and Space Sciences
• /
• v.3 no.1
• /
• pp.95-104
• /
• 2002

Numerical solution of the full Navier-Stokes equation as well as the energy equation has been obtained for the unsteady natural convection in a rectangular enclosure. One side wall was maintained at very high temperature simulating fires. Especially the effect of surface radiation was taken into account. While the enclosed air was assumed to be transparent, the internal walls directly interacted one another through the surface radiation. Due to a significant temperature difference in the flow field, the equation of state was used instead of the Boussinesq approximation. It was found that the rapid heating of the adiabatic ceiling and floor by the incoming radiation from the hot wall made the evolution at thermo-fluid field highly unstable in the initial period. Therefore, the secondary cells brought about at the floor region greatly affected the heat transfer mechanism inside the enclosure. The heat transfer rate was augmented by the radiation, resulting in requiring less time for the flow to reach the steady state. At the steady state neglecting radiation two internal hydraulic jumps were clearly observed in upper/left as well as in lower/right comer. However, the hydraulic jump in the lower/right comer could not be observed for the case including radiation due to its high momentum flow over the bottom wall. Radiation resulted in a faster establishment of the steady state phenomena.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.13 no.4
• /
• pp.1-10
• /
• 1976

It is a well-known phenomenon that, in the case of vibrations of an elastic body in a fluid such as water, the presence of the surrounding fluid has the effect of lowering the natural frequencies of the vibration as compared with those in air or vacuum on account of the increased inertia, i.e. added mass. In this report, defining the mass increase factor as the ratio of added mass to vibration mass of the body in air, the author investigated the mass increased factor of an elastic plate vibrating in the fluid. It is assumed that the edges of the plate are simply supported, and that the surrounding fluid is an infinite ideal one. For the problem formulation the elliptical cylindrical coordinate system is adopted, so that a rectangular plate may be represented by a sheet degenerated from an elliptical cylinder. By virtue of the coordinate system adopted, plates which are chordwisely finite and lengthwisely contineous could directly be treated, but plates which are chordwisely finite in both directions could not be treated directly. For the latter, hence, plates which are chordwisely finite and lengthwisely semi-finite are investigated as an appropriate approximation. Some examples of the mass increase factor are numerically calculated for the fundamental mode and modes of zero or one nodal line in each direction with the range of the aspect ratio from 1 to 10 or more.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.19 no.8
• /
• pp.1595-1611
• /
• 1994

In this paper, the numerical wave propagation properties of the finite difference-time domain(FD-TD) method is investigated as a discrete model describing electromagnetic(EM) wave propagation phenomena. The leap-frog approximation of Maxwell's curl equations in time-space simulates EM wave propagation in terms of the numerical characteristic and the domain of dependence. A geometrical interpretation of the FD-TD numerical procedure is presented. The numerical dispersion error due to the leap-frog approximation and its dependence on the stability factor are illustrated. The FD-TD method using the leap-frog approximation is inherently a descriptive model. Thus, not only any physical picture about EM wave propagation phenomena can be drawn through this model, but also physical or engineering parameters in the frequency domain can be extracted from descriptive results. E-plane filter characteristics in the WR-28 rectangular waveguide and reflection property of an inductive iris in the WR-90 rectangluar waveguide extracted from simulation of the FD-TD model is included.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.3 no.2
• /
• pp.1-8
• /
• 1999

The present paper describes an approximation for estimation of earthquake-induced hydrodynamic loads in rigid storage tank which accelerated in horizontal direction. The storage tank is vertically cylindrical, and the sectional shape may be circular, rectangular or irregular. The solution for harmonic excitation is studied based on velocity potential theory, and then the time domain solution for earthquake is obtained by using design response spectrum. As a result, earthquake load is influenced primarily by the inertia force of high frequency effective mass of the storage tank, responding to the characteristics of design response spectrum, tank sectional shape, and the ratio of tank base length to depth. Earthquake-induced hydrodynamic loads in rigid storage tank can be effectively obtained by using the high frequency approximation method in case of quite large, or small ratio of the tank base length to water depth.