• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.5
• /
• pp.10-17
• /
• 2003

Dynamic stability is critically required to stabilize an airship which is statically unstable. Numerical computations were performed in order to support and confirm the foced oscillation wind tunnel tests. To analyze the low-speed flow filed around the airship, a low-Mach number preconditioned method was applied. Using two computation methods, variations of the dynamic damping coefficients were examined. Numerical results show that it is dynamically stable for three directional moments, but unstable for normal or side force. It is revealed that the damping coefficients are more sensitive to the direction of the angular rate than the angle of attack or the magnitude og angular rate.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.9 no.3
• /
• pp.140-154
• /
• 1997

On November 15, 1994 at 03:17 local time, an earthquake of surface magnitude (M$_{s}$) 7.0 occurred on the northern coast of the Mindoro in Philippine. A major tsunami was generated by this earthquake, extremely large tsunami waves engulfed the Mindoro and the Verde islands. This tsunami caused tremendous casualities and damage. The tsunami propagated to the Luzon island and felt at the Batangas after 10 minutes. The present paper intends to understand the propagation and inundation this tsunami with the aid of numerical computation model and computer graphic aided video animation.n.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.8
• /
• pp.1429-1437
• /
• 1992

A new finite element formulation based on the relaxation type hereditary integral is presented for a time-domain analysis of isotropic, linear viscoelastic problems. The semi-discrete variational approximation and elastic-viscoelastic correspondence principle are used in the theoretical development of the proposed method. In a time-stepping procedure of final, linear algebraic system equations, only a small additional computation for past history is required since the equivalent stiffness matrix is constant. The viscoelasticity matrices are derived and the stress computation algorithm is given in matrix form. The effect of time increment and Gauss point numbers on the numerical accuracy is examined. Two dimensional numerical examples of plane strain and plane stress are solved and compared with the analytical solutions to demonstrate the versatility and accuracy of the present method.

• Journal of Korea Water Resources Association
• /
• v.36 no.1
• /
• pp.87-104
• /
• 2003

A numerical model of surface runoff and river flow has been applied to small- and large-size catchment areas in order to investigate the physical characteristics of river flow during flood period. Several refinements are made on the existing model SIRG-RS for the ways of rainfall input through surface runoff, river junction treatment and the computation of river flow on steep slope. For the computation of frictional forces, employed is the power law of friction factor which is a function of Reynolds number and relative roughness height. The empirical equation of friction factor is developed using recent field data as well as laboratory data. The refined model has been applied to small-size catchment area as well as large-size catchment area, and the computation results are found in good agreement with the observations in both cases.

• KSII Transactions on Internet and Information Systems (TIIS)
• /
• v.14 no.6
• /
• pp.2422-2443
• /
• 2020

This paper studies a single-user Mobile Edge Computing (MEC) system where mobile device (MD) includes an application consisting of multiple computation components or tasks with dependencies. MD can offload part of each computation-intensive latency-sensitive task to the AP integrated with MEC server. In order to accomplish the application faultlessly, we calculate out the optimal task offloading strategy in a time-division manner for a predetermined execution order under the constraints of limited computation and communication resources. The problem is formulated as an optimization problem that can minimize the energy consumption of mobile device while satisfying the constraints of computation tasks and mobile device resources. The optimization problem is equivalently transformed into solving a nonlinear equation with a linear inequality constraint by leveraging the Lagrange Multiplier method. And the proposed dual Bi-Section Search algorithm Bi-JOTD can efficiently solve the nonlinear equation. In the outer Bi-Section Search, the proposed algorithm searches for the optimal Lagrangian multiplier variable between the lower and upper boundaries. The inner Bi-Section Search achieves the Lagrangian multiplier vector corresponding to a given variable receiving from the outer layer. Numerical results demonstrate that the proposed algorithm has significant performance improvement than other baselines. The novel scheme not only reduces the difficulty of problem solving, but also obtains less energy consumption and better performance.

• Journal of the Korean Chemical Society
• /
• v.23 no.3
• /
• pp.125-131
• /
• 1979

A method for calculation of two center overlap integrals for a pair of Slater type orbitals was developed by Mulliken et al. In this method the spherical polar coordinates for a pair of Slater type orbitals located at two different points are required to be transformed into a spheroidal coordinate set for calculation of two center overlap integrals. A new method, the expansion method for spherical harmonics, in which Slater type orbitals, located at two different points, are expressed in a common coordinate system has been applied for computation of two center overlap integrals. The new method for computation of two center overlap integrals is required to translate Slater type orbitals centered at two different points into the reference point for computation of two center overlap integrals. This work has been expanded the expansion method for spherical harmonics for computation of two center overlap integrals to $|3s{\g}$, $|5s{\g}$ and $|5s{\g}$. Master formulas for two center overlap integrals are derived for these orbitals, using the general expansion formulas. The numerical values of the two center overlap integrals evaluated for a hypothetical NO molecule are in agreement with those of the previous works.

• Structural Engineering and Mechanics
• /
• v.64 no.1
• /
• pp.81-92
• /
• 2017

Many materials in engineering exhibit different modulus in tension and compression, which are known as bi-modulus materials. Based on the bi-modulus elastic theory, a modified semi-analytical model, by introducing a stress function, is established in this paper to study the mechanical response of a bi-modulus cylinder placed in an axisymmetric temperature field. Meanwhile, a numerical procedure to calculate the temperature stresses in bi-modulus structures is developed. It is proved that the bi-modulus solution can be degenerated to the classical same modulus solution, and is in great accordance with the solutions calculated by the semi-analytical model proposed by Kamiya (1977) and the numerical solutions calculated both by the procedure complied in this paper and by the finite element software ABAQUS, which demonstrates that the semi-analytical model and the numerical procedure are accurate and reliable. The result shows that the modified semi-analytical model simplifies the calculation process and improves the speed of computation. And the numerical procedure simplifies the modeling process and can be extended to study the stress field of bi-modulus structures with complex geometry and boundary conditions. Besides, the necessity to introduce the bi-modulus theory is discussed and some suggestions for the qualitative analysis and the quantitative calculation of such structure are proposed.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.25 no.10
• /
• pp.562-567
• /
• 2013

1D-numerical analysis of the network algorithm with the orifice equation for the relationship between pressure difference and flowrate has been mostly used to analyse leakage flow at the gap. In this study, a 3D-numerical method applying momentum loss model to the gap region in the computational domain is represented to reflect effectively the effect of leakage flow by determining the proportion of pressure difference to air passage velocity. While the 3D-numerical method is verified through the computation of the two compartments model, the numerical analysis of the stack effect in a building stairway is performed. As the temperature of air outside drops, the pressure in the upper stairway and leakage flowrate through the gap in the door rise. The change of gap area does not have an effect on pressure in the stairway for the analysis conditions.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.8
• /
• pp.1182-1190
• /
• 2003

A parallelized FEM code based on domain decomposition method has been recently developed for large-scale computational fluid dynamics. A 4-step splitting finite element algorithm is adopted for unsteady flow computation of the incompressible Navier-Stokes equation, and Smagorinsky LES model is chosen for turbulent flow computation. Both METIS and MPI Libraries are used for domain partitioning and data communication between processors, respectively. Tiburon model of Hyundai Motor Company is chosen as the computational model at Re=7.5 $\times$ 10$^{5}$ , which is based on the car height. The calculation is carried out under both the wind tunnel condition and the road condition using IBM SP parallel architecture at KISTI Super Computing Center. Compared with the existing experimental data, both the velocity and pressure fields are predicted reasonably well and the drag coefficient is in good agreement. Furthermore, it is confirmed that the drag under the road condition is smaller than that under the wind-tunnel condition.

• International Journal of Control, Automation, and Systems
• /
• v.3 no.3
• /
• pp.430-443
• /
• 2005

A general construction algorithm of the Grassmann space parameters in linear systems - so-called, the Plucker matrix, 'L' in m-input, p-output, n-th order static output feedback systems and the Plucker matrix, $'L^{aug}'$ in augmented (m+d)-input, (p+d)-output, (n+d)-th order static output feedback systems - is presented for numerical checking of necessary conditions of complete static and complete minimum d-th order dynamic output feedback pole-assignments, respectively, and also for discernment of deterministic computation condition of their pole-assignable real solutions. Through the construction of L, it is shown that certain generically pole-assignable strictly proper mp > n system is actually none pole-assignable over any (real and complex) output feedbacks, by intrinsic rank deficiency of some submatrix of L. And it is also concretely illustrated that this none pole-assignable mp > n system by static output feedback can be arbitrary pole-assignable system via minimum d-th order dynamic output feedback, which is constructed by deterministic computation under full­rank of some submatrix of $L^{aug}$.