• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.8
• /
• pp.693-698
• /
• 2007

A numerical method for accurate simulations of rotor aerodynamics is proposed. The numerical diffusion in the typically coarse grids away from the rotor blades is improved by implying a fourth-order of interpolation of local characteristic variables of the flow in the reconstruction stage of MUSCL approach in the framework of a finite volume formulation. In addition, different slope limiters are applied to the different characteristic fields, such as compressive limiters to linear characteristic fields to reduce the numerical dissipation whereas, diffusive limiters to nonlinear characteristic fields to increase numerical stability. Various exemplary problems related to the rotor aerodynamics applications are tested and the numerical results show a significant improvement in wake capturing capability. However, rotor aeroacoustic calculations show no meaningful difference over traditional MUSCL approach.

• Structural Engineering and Mechanics
• /
• v.65 no.1
• /
• pp.121-128
• /
• 2018

Three structure-dependent integration methods with no numerical dissipation have been successfully developed for time integration. Although these three integration methods generally have the same numerical properties, such as unconditional stability, second-order accuracy, explicit formulation, no overshoot and no numerical damping, there still exist some different numerical properties. It is found that TLM can only have unconditional stability for linear elastic and stiffness softening systems for zero viscous damping while for nonzero viscous damping it only has unconditional stability for linear elastic systems. Whereas, both CEM and CRM can have unconditional stability for linear elastic and stiffness softening systems for both zero and nonzero viscous damping. However, the most significantly different property among the three integration methods is a weak instability. In fact, both CRM and TLM have a weak instability, which will lead to an adverse overshoot or even a numerical instability in the high frequency responses to nonzero initial conditions. Whereas, CEM possesses no such an adverse weak instability. As a result, the performance of CEM is much better than for CRM and TLM. Notice that a weak instability property of CRM and TLM might severely limit its practical applications.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.15 no.1
• /
• pp.1-8
• /
• 1990

In this paper, the transient electromagneti wave scattering at dielectric cylinder is studied by new numerical analysis method. Basic formulation of boundary integral equation (BIE) for numerical method is started weighted residual technique. BIE is made to two simultaneous equation at surface inner and outside point of dielectric cylinder in extended boundary condition (EBC) and surface boundary condition (SBC). Numerical method is used Boundary element method (BEM) that is two form, one is direct method and the other is indirect method, so that this method that transformes operator inversion martics is used numerical analysis. A good agreement of this numerical solution and the other results is obtained.

• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.50-57
• /
• 2008

The prediction of hydrate pellet decomposition characteristics is required to design the regasification process of GTS (gas to solid) technology, which is considered as an economic alternative for LNG technology to transport natural gas produced from small and stranded gas wells. Mathematical model based on the conservation principles, the phase equilibrium relation, equation of gas state and phase change kinetics was set up and numerical solution procedure employing volume averaged fixed grid formulation and extended enthalpy method are implemented. Initially, porous methane hydrate pellet is at uniform temperature and pressure within hydrate stable region. The pressure starts to decrease with a fixed rate down to the final pressure and is kept constant afterwards while the bounding surface of pellet is heated by convection. The predicted convective heat and mass transfer accompanied by the decomposed gas flow through hydrate/ice solid matrix is reported focused on the comparison of spherical and cylindrical pellets having the same effective radius.

• 한국연소학회:학술대회논문집
• /
• 2003.05a
• /
• pp.57-62
• /
• 2003

A general procedure of obtaining reliable one-step kinetics model for hydrocarbon mixture from the fully detailed chemistry is described iin this study. One-step theoretical formulation of the induction parameter model IPM uses a theoretical reconstruction of the induction time database obtained from a detailed kinetics library. Non-dimensional induction time calculations is compared with that of detailed kinetics. The IPM was latter implemented to fluid dynamics code and applied for the numerical simulation of detonation wave propagation. The numerical results including the numerical smoked-foil record show the all the details of the detonation wave propagation characteristics at the cost around 1/100 of the detailed kinetics calculation.

• Computational Structural Engineering
• /
• v.6 no.2
• /
• pp.103-114
• /
• 1993

An advanced time-domain finite element formulation is presented for the displacement and stress analysis of isotropic, linear viscoelastic problems of a hereditary-type constitutive law. The semidiscrete finite element method with linear time-stepping scheme and an elastic-viscoelastic correspondence principle are used in the theoretical development. An efficient treatment of hereditary integral is introduced to improve the numerical accuracy, to reduce the computation time, and to avoid the use of large memory storage. Two-dimensional numerical examples of plane strain and plane stress are solved under the assumption on the material property of being elastic in dilatation and like three-element Voigt model in distorsion, and compared with the analytical solutions and the past numerical results to show the versatility and efficiency of the proposed method.

• Steel and Composite Structures
• /
• v.3 no.2
• /
• pp.75-95
• /
• 2003

A numerical model is presented to simulate the mechanical behaviour of composite steel and concrete columns taking into account the interaction between the hollow steel section and the concrete core. The model, based on displacement finite element methods with an Updated Lagrangian formulation, allows for geometrical and material non linearities combined with heating over all or a part of the section and column length. Comparisons of numerical calculations made using the model with 33 fire resistance tests show that the model is able to predict the fire resistance, expressed in minutes of fire exposure, of composite columns with a good accuracy.

• Structural Engineering and Mechanics
• /
• v.65 no.5
• /
• pp.513-521
• /
• 2018

In this paper, the resonance response of spar-type floating platform in coupled heave and pitch motion is investigated using a CPU time-effective numerical method. A coupled nonlinear 2-DOF equation of motion is derived based on the potential wave theory and the rigid-body hydrodynamics. The transient responses are solved by the fourth-order Runge-Kutta (RK4) method and transformed to the frequency responses by the digital Fourier transform (DFT), and the first-order approximation of heave response is analytically derived. Through the numerical experiments, the theoretical derivation and the numerical formulation are verified from the comparison with the commercial software AQWA. And, the frequencies of resonance arising from the nonlinear coupling between heave and pitch motions are investigated and justified from the comparison with the analytically derived first-order approximation of heave response.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.3 no.4
• /
• pp.143-148
• /
• 1983

A general numerical procedure for the time dependent analysis of concrete structures due to the effects of creep, shrinkage, aging of concrete and temperature variations is presented. Creep is represented by a superposition integral, and is considered as dependent on concrete age and temperature. An efficient numerical formulation of creep is possible through the use of a special exponential form of specific creep function which forms the kernel function of the superposition integral. The accuracy of this method is demonstrated by a numerical example.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.465-467
• /
• 2004

Finite element analysis of welding processes, which entail phase evolution, heat transfer and deformation, is considered in this paper. Attention focuses on numerical implementation of the thermo-elastic-plastic constitutive equation proposed by Leblond et al in consideration of the transformation plasticity. Based upon the multiplicative decomposition of deformation gradient, hyperelastic formulation is employed for efficient numerical integration, and the algorithmic consistent moduli for elastic-plastic deformations including transformation plasticity are obtained in the closed form. The convergence behavior of the present implementation is demonstrated via a couple of numerical examples. Several locking phenomena removed by Solid-shell element.