• Journal of computational fluids engineering
• /
• v.10 no.2
• /
• pp.60-68
• /
• 2005

For study on the unsteady wall interference effect, flows around a forced oscillating airfoil in closed test-section wind tunnels have been numerically investigated by solving compressible Navier-Stokes equations. The numerical scheme is based on a node-based finite-volume method with the Roe's flux-difference splitting and an implicit time-integration method coupled with dual time-step sub-iteration. The Spalart-Allmaras one-equation model is employed for the turbulence effect. The computed results of the oscillating airfoil having a thin wake showed that the lift curve slope is increased and the magnitude of hysteresis loop is reduced by the interference effects. Since the vortex around the airfoil is generated and convected downstream faster than the free-air condition, the phase of lift, drag and pitching moment coefficients was shifted. The pressure on the test section wall shows harmonic terms having the oscillating frequency contained in the wail effect.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.11
• /
• pp.1075-1080
• /
• 2010

For the aeroelastic analysis of a wing with friction damping, coupled time integration method was used to obtain time responses in the subsonic and transonic regions. To take into account aerodynamic nonlinearity induced by shock wave on the lifting surface, transonic small disturbance equation with in-phase periodic boundary condition was used for unsteady aerodynamic calculation. For 2-DOF airfoil system with displace-dependent friction dampers, the effects of normal load slope and Mach number on flutter boundary were investigated.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2001.11b
• /
• pp.983-988
• /
• 2001

Dynamic behaviors and stability of an optical disk drive coupled with an automatic ball balancer(ABB) are analyzed by a theoretical approach. The feeding system is modeled a rigid body with six degree-of-freedom. Using Lagrange's equation, we derive the nonlinear equations of motion for a non-autonomous system with respect to the rectangular coordinate. To investigate the dynamic stability of the system in the neighborhood of the equilibrium positions, the monodromy matrix technique is applied to the perturbed equations. On the other hand, time responses are computed by the Runge-Kutta method. We also investigate the effects of the damping coefficient and the position of ABB on the dynamic behaviors of the system.

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

The propagation properties of various pulse signal types(square pulse, Gaussian pulse, trapezoid pulse, RF pulse) on coupled microstrip lines are investigated. Numerical integration technique which has its accuracy and is easily simulated, is used to obtain the time domain response of pulse signals. Frequency-dependent characteristics of coupled microstrip line is obtained using Jansen's approximate equation. The propagation properties of pulse signal on coupled microstrip lines is analyzed regarding to its geometric structure (relative permittivity ${varepsilon}_r$ substrate height h, strip width w of the microstrip line) and pulse width ${\tau}$ of signal pulse. The simulation results show that space between two lines is very significant parameter in pulse distortion in comparison of any other parameters. The results of this paper are compatible to the trade-off determination of relative permittivity, substrate height, strip width and pulse width of signal pulse when a design of MIC and MMIC is necessary.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.9
• /
• pp.1905-1913
• /
• 2002

Using the finite element method, this study investigates the dynamic time responses of a flexible spinning disk of which axis of rotation is misaligned with the axis of symmetry. The misalignment between the axes of symmetry and rotation is one of the major vibration sources in optical disk drives such as CD-ROM, CD-R, CD-RW and DVD drives. Based upon the Kirchhoff plate theory and the von-Karman strain theory, three coupled equations of motion for the misaligned disk are obtained: two of the equations are for the in-plane motion while the other is for the out-of-plane motion. After transforming these equations into two weak forms for the in-plane and out-of-plane motions, the weak forms are discretized by using newly defined annular sector finite elements. Applying the generalized-$\alpha$ time integration method to the discretized equations, the time responses and the displacement distributions are computed and then the effects of the misalign ment on the responses and the distributions are analyzed. The computation results show that the misalignment has an influence on the magnitudes of the in-plane displacements and it results in the amplitude modulation or the beat phenomenon in the time responses of the out-of-plane displacement.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.24 no.2
• /
• pp.87-93
• /
• 2020

Seismic responses due to the dynamic coupling between a primary structure and secondary system connected to a structure are analyzed in this study. The seismic responses are compared based on dynamic coupling criteria and according to the error level in the natural frequency, with the recent criteria being reliant on the error level in the spectral displacement response. The acceleration responses and relative displacement responses of a primary structure and a secondary system for a coupled model and two different decoupled models of two degrees-of-freedom system are calculated by means of the time integration method. Errors in seismic responses of the uncoupled models are reduced with the recent criteria. As the natural frequency of the secondary system increases, error in the natural frequency decreases, but seismic responses of uncoupled models can be underestimated compared to that of coupled model. Results in this paper can help determine dynamic coupling and predict uncoupled models' response conservatism.

• Interaction and multiscale mechanics
• /
• v.1 no.1
• /
• pp.1-31
• /
• 2008

The present paper presents multiscale modelling via coupling of the discrete and finite element methods. Theoretical formulation of the discrete element method using spherical or cylindrical particles has been briefly reviewed. Basic equations of the finite element method using the explicit time integration have been given. The micr-macro transition for the discrete element method has been discussed. Theoretical formulations for macroscopic stress and strain tensors have been given. Determination of macroscopic constitutive properties using dimensionless micro-macro relationships has been proposed. The formulation of the multiscale DEM/FEM model employing the DEM and FEM in different subdomains of the same body has been presented. The coupling allows the use of partially overlapping DEM and FEM subdomains. The overlap zone in the two coupling algorithms is introduced in order to provide a smooth transition from one discretization method to the other. Coupling between the DEM and FEM subdomains is provided by additional kinematic constraints imposed by means of either the Lagrange multipliers or penalty function method. The coupled DEM/FEM formulation has been implemented in the authors' own numerical program. Good performance of the numerical algorithms has been demonstrated in a number of examples.

• Ocean Systems Engineering
• /
• v.4 no.2
• /
• pp.83-97
• /
• 2014

Wave propagation in a three-dimensional (3D) fully nonlinear numerical wave tank (NWT) is studied based on velocity potential theory. The governing Laplace equation with fully nonlinear boundary conditions on the moving free surface is solved using the indirect desingularized boundary integral equation method (DBIEM). The fourth-order predictor-corrector Adams-Bashforth-Moulton scheme (ABM4) and mixed Eulerian-Lagrangian (MEL) method are used for the time-stepping integration of the free surface boundary conditions. A smoothing algorithm, B-spline, is applied to eliminate the possible saw-tooth instabilities. The artificial wave speed employed in MTF (multi-transmitting formula) approach is investigated for fully nonlinear wave problem. The numerical results from incorporating the damping zone (DZ), MTF and MTF coupled DZ (MTF+DZ) methods as radiation condition are compared with analytical solution. An effective MTF+DZ method is finally adopted to simulate the 3D linear wave, second-order wave and irregular wave propagation. It is shown that the MTF+DZ method can be used for simulating fully nonlinear wave propagation very efficiently.

• Wind and Structures
• /
• v.21 no.5
• /
• pp.461-487
• /
• 2015

In recent years, the Graphics Processing Unit (GPU) has become a competitive computing technology in comparison with the standard Central Processing Unit (CPU) technology due to reduced unit cost, energy and computing time. This paper describes the derivation and implementation of GPU-based algorithms for the analysis of wind loading uncertainty on high-rise systems, in line with the research field of probability-based wind engineering. The study begins by presenting an application of the GPU technology to basic linear algebra problems to demonstrate advantages and limitations. Subsequently, Monte-Carlo integration and synthetic generation of wind turbulence are examined. Finally, the GPU architecture is used for the dynamic analysis of three high-rise structural systems under uncertain wind loads. In the first example the fragility analysis of a single degree-of-freedom structure is illustrated. Since fragility analysis employs sampling-based Monte Carlo simulation, it is feasible to distribute the evaluation of different random parameters among different GPU threads and to compute the results in parallel. In the second case the fragility analysis is carried out on a continuum structure, i.e., a tall building, in which double integration is required to evaluate the generalized turbulent wind load and the dynamic response in the frequency domain. The third example examines the computation of the generalized coupled wind load and response on a tall building in both along-wind and cross-wind directions. It is concluded that the GPU can perform computational tasks on average 10 times faster than the CPU.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.20 no.2
• /
• pp.7-13
• /
• 2019

This study aims to develop a FE Model to simulate dissimilar friction stir welding and to address its potential for fundamental analysis and practical applications. The FE model is based on Coupled Eulerian-Lagrangian approach. Multiphysics systems are calculated using explicit time integration algorithm, and heat generations by friction and inelastic heat conversion as well as heat transfer through the bottom surface are included. Using the developed model, friction stir welding between an Al6061T6 plate and an AZ61 plate were simulated. Three simulations are carried out varying the welding parameters. The model is capable of predicting the temperature and plastic strain fields and the distribution of void. The simulation results showed that temperature was generally greater in Mg plates and that, as a rotation speed increase, not the maximum temperature of Mg plate increased, but did the temperature of Al plate. In addition, the model could predict flash defects, however, the prediction of void near the welding tool was not satisfactory. Since the model includes the complex physics closely occurring during FSW, the model possibly analyze a lot of phenomena hard to discovered by experiments. However, practical applications may be limited due to huge simulation time.