• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.16 no.3
• /
• pp.193-204
• /
• 2012

This paper is comparing numerical schemes for a differential equation with convection and fourth-order diffusion. Our model equation is $h_t+(h^2-h^3)_x=-(h^3h_{xxx})_x$, which arises in the context of thin film flow driven the competing effects of an induced surface tension gradient and gravity. These films arise in thin coating flows and are of great technical and scientific interest. Here we focus on the several numerical methods to apply the model equation and the comparison and analysis of the numerical results. The convection terms are treated with well known WENO methods and the diffusion term is treated implicitly. The diffusion and convection schemes are combined using a fractional step-splitting method.

• 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.

• Steel and Composite Structures
• /
• v.4 no.1
• /
• pp.23-36
• /
• 2004

The inelastic lateral-distortional buckling of doubly-symmetric hot-rolled I-section beam-columns subjected to a concentric axial force and uniform bending with elastic restraint which produce single curvature is investigated in this paper. The numerical model adopted in this paper is an energy-based method which leads to the incremental and iterative solution of a fourth-order eigenproblem, with very rapid solutions being obtained. The elastic restraint considered in this paper is full restraint against translation, but torsional restraint is permitted at the tension flange. Hitherto, a numerical method to analyse the elastic and inelastic lateral-distortional buckling of restrained or unrestrained beam-columns is unavailable. The prediction of the inelastic lateral-distortional buckling load obtained in this study is compared with the inelastic lateral-distortional buckling of restrained beams and the inelastic lateral-torsional buckling solution, by suppressing the out-of-plane web distortion, is published elsewhere and they agree reasonable well. The method is then extended to the lateral-distortional buckling of continuously restrained doubly symmetric I-sections to illustrate the effect of web distortion.

• 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.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.6
• /
• pp.447-457
• /
• 2018

In this study, a two-dimensional fully nonlinear transient wave numerical tank was developed using a desingularized indirect boundary integral equation method. The desingularized indirect boundary integral equation method is simpler and faster than the conventional boundary element method because special treatment is not required to compute the boundary integral. Numerical simulations were carried out in the time domain using the fourth order Runge-Kutta method. A mixed Eulerian-Lagrangian approach was adapted to reconstruct the free surface at each time step. A numerical damping zone was used to minimize the reflective wave in the downstream region. The interpolating method of a Gaussian radial basis function-type artificial neural network was used to calculate the gradient of the free surface elevation without element connectivity. The desingularized indirect boundary integral equation using an isolated point source and radial basis function has no need for information about the element connectivity and is a meshless method that is numerically more flexible. In order to validate the accuracy of the numerical wave tank based on the desingularized indirect boundary integral equation method and meshless technique, several numerical simulations were carried out. First, a comparison with numerical results according to the type of desingularized source was carried out and confirmed that continuous line sources can be replaced by simply isolated sources. In addition, a propagation simulation of a $2^{nd}$-order Stokes wave was carried out and compared with an analytical solution. Finally, simulations of propagating waves in shallow water and propagating waves over a submerged bar were also carried and compared with published data.

• Wind and Structures
• /
• v.23 no.5
• /
• pp.421-447
• /
• 2016

Wind-induced and earthquake-induced excitations on tall structures can be effectively controlled by Tuned Liquid Damper (TLD). This work presents a numerical simulation procedure to study the performance of tuned liquid tank- structure system through ${\sigma}$-transformation based fluid-structure coupled solver. For this, a 'C' based computational code is developed. Structural equations are coupled with fluid equations in order to achieve the transfer of sloshing forces to structure for damping. Structural equations are solved by fourth order Runge-Kutta method while fluid equations are solved using finite difference based sigma transformed algorithm. Code is validated with previously published results. The minimum displacement of structure is observed when the resonance condition of the coupled system is satisfied through proper tuning of TLD. Since real-time excitations are random in nature, the performance study of TLD under random excitation is also carried out in which the Bretschneider spectrum is used to generate the random input wave.

• Kyungpook Mathematical Journal
• /
• v.61 no.3
• /
• pp.591-612
• /
• 2021

In this paper, a uniformly convergent numerical scheme is designed for solving singularly perturbed reaction-diffusion problems. The problem is converted to an equivalent weak form and then a Galerkin finite element method is used on a piecewise uniform Shishkin mesh with linear basis functions. The convergence of the developed scheme is proved and it is shown to be almost fourth order uniformly convergent in the maximum norm. To exhibit the applicability of the scheme, model examples are considered and solved for different values of a singular perturbation parameter ε and mesh elements. The proposed scheme approximates the exact solution very well.

• 한국가스학회:학술대회논문집
• /
• 1997.09a
• /
• pp.14-19
• /
• 1997

A numerical method for determining the temperature vartiation in a natural gas transmission line is presented. By considering an element of the gas pipeline and assuming radially lumped heat transfer at steady-state conditions, the energy equation is developed. The integration of the developed nonlinear differential equation is done numerically using the fourth order Runge-Kutta scheme. The results of the present study have been compared with the results of Coulter equations, and show a fairly good agreement.

• Structural Engineering and Mechanics
• /
• v.67 no.6
• /
• pp.603-616
• /
• 2018

In this study, the acceleration vector in each time step is assumed to be a mth order time polynomial. By using the initial conditions, satisfying the equation of motion at both ends of the time step and minimizing the square of the residual vector, the m+3 unknown coefficients are determined. The order of accuracy for this approach is m+1, and it has a very low dispersion error. Moreover, the period error of the new technique is almost zero, and it is considerably smaller than the members of the Newmark method. The proposed scheme has an appropriate domain of stability, which is greater than that of the central difference and linear acceleration techniques. The numerical tests highlight the improved performance of the new algorithm over the fourth-order Runge-Kutta, central difference, linear and average acceleration methods.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.5 no.4
• /
• pp.350-356
• /
• 1993

We carry out a numerical calculation to understand the nonlinear wave deformation around breakwaters using the Boussinesq equation, which is weakly nonlinear and weakly dispersive shallow water equation. A numerical method based on a finite element scheme and fourth order Runge-Kutta algorithm is employed to investigate the diffraction of incident waves by the breakwater. As a computational model, two-dimensional wave flume is treated. The breakwaters is perpendicular to the side wall of a channel. From the numerical results, the wave deformations according to the change of the length and the thickness of breakwaters are investigated. We also investigate the effect of the nonlinearity by comparing the results with the linear solutions.