• Wind and Structures
• /
• v.5 no.2_3_4
• /
• pp.369-378
• /
• 2002

The application of Large Eddy Simulation (LES) in a curvilinear coordinate system to the flow around a square cylinder is presented. In order to obtain sufficient resolution near the side of the cylinder, we use an O-type grid. Even with a curvilinear coordinate system, it is difficult to avoid the numerical oscillation arising in high-Reynolds-number flows past a bluff body, without using an extremely fine grid used. An upwind scheme has the effect of removing the numerical oscillations, but, it is accompanied by numerical dissipation that is a kind of an additional sub-grid scale effect. Firstly, we investigate the effect of numerical dissipation on the computational results in a case where turbulent dissipation is removed in order to clarify the differences between the effect of numerical dissipation. Next, the applicability and the limitations of the present method, which combine the dynamic SGS model with acceptable numerical dissipation, are discussed.

• Journal of Industrial Technology
• /
• v.10
• /
• pp.69-72
• /
• 1990

In this paper the numerical method for the study of wave reflection from and transmission through submerged permeable breakwaters using the boundary element method is developed. The numerical analysis technique is based on the wave pressure function instead of velocity potential because it is difficult to define the velocity potential in the each region arising the energy dissipation. Also, the non-linear energy dissipation within the submerged porous structure is simulated by introducing the linear dissipation coefficient and the tag mass coefficient equivalent to the non-linear energy dissipation. For the validity of this analysis technique, the numerical results obtained by the present boundary element method are compared with those obtained by the other computation method. Good agreements are obtained and so the validity of the present numerical analysis technique is proved.

• Structural Engineering and Mechanics
• /
• v.85 no.2
• /
• pp.207-216
• /
• 2023

Three time-discontinuous Galerkin quadrature element methods (TDGQEMs) are developed for structural dynamic problems. The weak-form time-discontinuous Galerkin (TDG) statements, which are capable of capturing possible displacement and/or velocity discontinuities, are employed to formulate the three types of quadrature elements, i.e., single-field, single-field/least-squares and two-field. Gauss-Lobatto quadrature rule and the differential quadrature analog are used to turn the weak-form TDG statements into a system of algebraic equations. The stability, accuracy and numerical dissipation and dispersion properties of the formulated elements are examined. It is found that all the elements are unconditionally stable, the order of accuracy is equal to two times the element order minus one or two times the element order, and the high-order elements possess desired high numerical dissipation in the high-frequency domain and low numerical dissipation and dispersion in the low-frequency domain. Three fundamental numerical examples are investigated to demonstrate the effectiveness and high accuracy of the elements, as compared with the commonly used time integration schemes.

• Earthquakes and Structures
• /
• v.7 no.2
• /
• pp.159-178
• /
• 2014

Although the family methods with unconditional stability and numerical dissipation have been developed for structural dynamics they all are implicit methods and thus an iterative procedure is generally involved for each time step. In this work, a new family method is proposed. It involves no nonlinear iterations in addition to unconditional stability and favorable numerical dissipation, which can be continuously controlled. In particular, it can have a zero damping ratio. The most important improvement of this family method is that it involves no nonlinear iterations for each time step and thus it can save many computationally efforts when compared to the currently available dissipative implicit integration methods.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2003.03a
• /
• pp.275-285
• /
• 2003

Pinching is an important property of reinforced concrete member which characterizes its cyclic behavior. In the present study, numerical studies were performed to investigate the characteristics and mechanisms of pinching behavior and the energy dissipation capacity of flexure-dominated reinforced concrete members. By analyzing existing experimental studies and numerical results, it was found that energy dissipation capacity of a member is directly related to energy dissipated by re-bars rather than concrete that is a brittle material, and that it is not related to magnitude of axial compressive force applied to the member. Therefore, for a member with specific arrangement and amount of re-bars, the energy dissipation capacity remains uniform regardless of the flexural strength that is changed by the magnitude of axial force applied. Due to the uniformness of energy dissipation capacity pinching appears in axial compression member. The flexural pinching that is not related to shear force becomes conspicuous as the flexural strength increases relatively to the uniform energy dissipation capacity. Based on the findings, a practical method for estimating energy dissipation capacity and damping modification factor was developed and verified with existing experiments.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.10a
• /
• pp.11-19
• /
• 2001

An adaptive nonlinear artificial dissipation model is presented for performing aeroacoustic computations by the high-order and high-resolution numerical schemes based on the central finite differences. An effective formalism of it is devised by combining a selective background smoothing term and a well-established nonlinear shock-capturing term which is for the temporal accuracy as well as the numerical stability. A conservative form of the selective background smoothing term is presented to keep accurate phase speeds of the propagating nonlinear waves. The nonlinear shock-capturing term that has been modeled by the second-order derivative term is combined with it to improve the resolution of discontinuities and stabilize the strong nonlinear waves. It is shown that the improved artificial dissipation model with an adaptive control constant which is independent of problem types reproduces the correct profiles and speeds of nonlinear waves, suppresses numerical oscillations near discontinuity and avoids unnecessary damping on the smooth linear acoustic waves. The feasibility and performance of the adaptive nonlinear artificial dissipation model are investigated by the applications to actual computational aeroacoustics problems.

• Structural Engineering and Mechanics
• /
• v.89 no.4
• /
• pp.411-420
• /
• 2024

The purpose of this research was to investigate the cyclic behavior of the D-shaped dampers (DSD). Similarly, at first, the numerical model was calibrated using the experimental sample. Then, parametric studies were conducted in order to investigate the effect of the radius and thickness of the damper on energy dissipation, effective and elastic stiffness, ultimate strength, and equivalent viscous damping ratio (EVDR). An analytical equation for the elastic stiffness of the DSD was also proposed, which showed good agreement with experimental results. Additionally, approximate equations were introduced to calculate the elastic and effective stiffness, ultimate strength, and energy dissipation. These equations were presented according to the curve fitting technique and based on numerical results. The results indicated that reducing the radius and increasing the thickness led to increased energy dissipation, effective stiffness, and ultimate strength of the damper. On the other hand, increasing the radius and thickness resulted in an increase in EVDR. Moreover, the ratio of effective stiffness to elastic stiffness also played a crucial role in increasing the EVDR. The thickness and radius of the damper were evaluated as the most effective dimensions for reducing energy dissipation and EVDR.

• Bulletin of the Korean Mathematical Society
• /
• v.54 no.1
• /
• pp.145-158
• /
• 2017

We present a finite difference method for solving the Ohta-Kawasaki model, representing a model of mesoscopic phase separation for the block copolymer. The numerical methods for solving the Ohta-Kawasaki model need to inherit the mass conservation and energy dissipation properties. We prove these characteristic properties and solvability and unconditionally gradient stability of the scheme by using Hessian matrices of a discrete functional. We present numerical results that validate the mass conservation, and energy dissipation, and unconditional stability of the method.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.1
• /
• pp.29-45
• /
• 2013

Central schemes offer a simple and versatile approach for computing approximate solutions of nonlinear systems of hyperbolic conservation laws. However, there are large numerical dissipation in case of contact discontinuity. We study semi-discrete central upwind scheme by changing flux functions to reduce the numerical dissipation and we perform numerical computations for various problems in case of contact discontinuity.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.30 no.4
• /
• pp.143-152
• /
• 2018

Numerical convergence and stability of TVD schemes have been applied in the FUNWAVE-TVD model were compared. The fourth order accurate MUSCL-TVD scheme using minmod limiter suggested by Yamamoto and Daiguji (1993), the fourth order accurate MUSCL-TVD scheme using van-Leer limiter suggested by Erduran et al. (2005) and the second order accurate MUSCL-TVD scheme using van-Leer limiter in Zhou et al. (2001) were compared. Comparisons of the numerical scheme were conducted with experimental data of Vincent and Briggs irregular wave experiments. In comparison with the fourth order accurate scheme using van-Leer limiter, the fourth order accurate scheme using minmod limiter is less dissipative but required lower CFL condition for stable numerical solution. On the other hand, the scheme using van-Leer limiter required smaller resolution spatial grid due to numerical dissipation, but relatively higher CFL condition can be used compared to the scheme using minmod limiter. In the breaking wave experiments which were conducted using high resolution spatial grid to reduce numerical dissipation, the characteristic of the schemes can be clearly observed. Numerical instabilities and blow-up of the numerical solutions were found in the irregular wave breaking simulation with the scheme using minmod limiter. However, the simulation can be completed with the scheme using van-Leer limiter, but required low CFL condition. Good agreements with the observed data were also observed in the results using van-Leer limiter.