• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.11
• /
• pp.1538-1551
• /
• 1997

Kinetic energy conservation for fixed and moving grids is examined in time-accurate finite element computation of fully unsteady inviscid flows. As numerical algorithms, fractional step method (FSM) and modified SIMPLE are used. To simulate the flow in moving grid system, arbitrary Lagrangian-Eulerian (ALE) method is adopted. In the present study, the energy conserving time integration rule for finite element algorithm is proposed and discussed schematically. It is shown that the discretization by Crank-Nicolson in time and Galerkin (central difference) in space must be used to ensure energy conservation. The developed code has been tested for a standing vortex in fixed or moving grid system, sloshing in a tank and propagation of a solitary wave, and has been shown to be a completely energy conserving algorithm.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.6
• /
• pp.1365-1381
• /
• 1990

A $C^{0}$ continuous displacement finite element method based on a higher-order shear deformation theory is employed in the prediction of the transient response of laminated composite plates subjected to low-velocity impact. A modified contact law was applied to calculate the contact force during impact. The discrete element chosen is a nine-noded quadrilateral with 5 degree-of-freedom per node. The Wilson-.theta. time integration algorithm is used for solving the time dependent equations of the impactor and the central difference method was adopted to perform time integration of the plate. Numerical results, including the contact force history, deflection, and velocity history, are presented. Comparisons of numerical results using a higher order theory and a first-order theory show that using a higher order theory provides more accurate results. Effects of boundary condition, impact velocity, and mass of the impactors are also discussed.d.

• Journal of Power System Engineering
• /
• v.11 no.1
• /
• pp.20-27
• /
• 2007

In this paper, finite element modeling was performed for vibration analysis of a rotor system installed in sunroof motor, and analysis process was developed for natural frequency and unbalance response analysis. At the same time, to reduce analysis modeling error caused by the difference between analysis and measured values, finite element model updating was conducted using an optimization algorithm, i.e. hybrid genetic algorithm and simulated annealing (HGASA) method. For this end experimental modal test was carried out and by using the measured frequency response function (FRF), model updating was performed considering both cases where core coil was removed and included. And acceptable result was obtained. Also, dynamic property coefficient of bush bearing which influences vibration response of the rotor system was estimated.

• Food Science of Animal Resources
• /
• v.33 no.1
• /
• pp.31-38
• /
• 2013

If the time-temperature indicator (TTI) experienced a different temperature than the accompanied packaged food, influenced by heat transfer between the TTI, package, and ambient air, TTI would incorrectly predict the food quality changes with temperature. Temperature distributions of a finite slab with different sizes, representing beef packaged with TTI, were estimated by the finite element method (FEM). The thermal properties of the beef and TTI, such as heat capacity, density, and heat conductivity, were estimated from the relevant equations using their chemical compositions. The FEM simulations were performed for three cases: different locations of TTIs on the beef, different thicknesses of beef, and non-isothermal conditions of ambient air. The TTIs were mounted in four different locations on the beef. There was little difference in temperature between four locations of the TTI on the package surface. As the thickness of the slab increased, the temperature of the TTI changed faster, followed by the corner surface, as well as middle and bottom parts, indicating the possible error for temperature agreement between the TTI and the slab. Consequently, it was found that any place on the package could be selected for TTI attachment, but the package size should carefully be determined within a tolerable error of temperature.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.6 no.3
• /
• pp.139-149
• /
• 2002

An explicit direct time integration method based solution algorithm is proposed to predict dynamic postbuckling response of thin plates. Based on the von Karman's plate equations and Marquerre's shallow shell theory, a rectangular plate finite element is formulated and utilized in this study. The element formulation takes into account geometrical nonlinearity and initial deflection of plates. The solution algorithm employs the central difference method. Using the computer program developed by the authors, dynamic postbuckling behavior of elastic thin plates under impact loading is investigated by considering the time variation of load and load duration. The efficiency of the proposed solution algorithm is examined through illustrative numerical examples.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.14 no.3
• /
• pp.189-200
• /
• 2010

A numerical method is proposed and analyzed to approximate a mathematical model of age-dependent population dynamics with spatial diffusion. The model takes a form of nonlinear and nonlocal system of integro-differential equations. A finite difference method along the characteristic age-time direction is considered and primal mixed finite elements are used in the spatial variable. A priori error estimates are derived for the relevant variables.

• Proceedings of the Korea Concrete Institute Conference
• /
• 1999.04a
• /
• pp.799-806
• /
• 1999

This study describes the results of experiment and numerical analysis for heating image by thermographic method when the size of void in concrete are changed. By comparing analytical solution by finite element method with measured image by thermography, the relationships between the surface temperature which can be confirmed by this method, the size of void and optimum time for detection of void and the difference of temperature are cleared.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.8 no.5
• /
• pp.486-487
• /
• 1997

This paper presents the pure capacitive lumped element equivalent circuits for several coplanar waveguide(CPW) discontinuities such as an open-end, an open-end with connected ground planes, a gap and an open-end CPW stub and gives their capacitive element values as a function of physical dimensions of the discontinuity and the frequency for a specific substrate. The capacitive element values are determined from the scattering parameters which are obtained using the finite-difference time-domain(FDTD) method. For an open-end, an open-end with connected ground planes and a gap, the numerical results of the FDTD are compared with the quasi-static results which are calculated using the three- dimensional finite difference method(3D-FDM).

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.22 no.7
• /
• pp.1459-1467
• /
• 1997

The current on the thin planar structure as an element of the transversely fed electromagnetically coupled(EMC) microstrip dipole array antenna is obtained by using the integral forms of the finite difference time domain(FDTD) method. This method was applied to calculating the optimum current distribution (Doplh-Tchebyscheff distribution) of each dipole element on the feed line as a function of their offset positions for the narrow main beam width and the side beam level below -20 dB. The current on each dipole substitutes for the electric and magnetic current densities on the virtual surface of the FDTD calculation to express the far field intensity, the calculation time and the computer memeory can be reduced to about 80% and 1.3 Mbyte, respectively. The calculated radiation patterns are compared to the measured values and these are in good agreement.

• Bulletin of the Korean Mathematical Society
• /
• v.61 no.1
• /
• pp.93-115
• /
• 2024

In this paper, the Gauge-Uzawa methods for the Darcy-Brinkman equations driven by temperature and salt concentration (DBTC) are proposed. The first order backward difference formula is adopted to approximate the time derivative term, and the linear term is treated implicitly, the nonlinear terms are treated semi-implicit. In each time step, the coupling elliptic problems of velocity, temperature and salt concentration are solved, and then the pressure is solved. The unconditional stability and error estimations of the first order semi-discrete scheme are derived, at the same time, the unconditional stability of the first order fully discrete scheme is obtained. Some numerical experiments verify the theoretical prediction and show the effectiveness of the proposed methods.