• Journal of KSNVE
• /
• v.7 no.5
• /
• pp.773-780
• /
• 1997

Using Generalized Inverse Method presented by Udwadia and Kalaba in 1992, we can obtain equations to exactly describe the motion of constrained systems. When the differential equations are numerically integrated by any numerical integration scheme, the numerical results are generally found to veer away from satisfying constraint equations. Thus, this paper deals with the numerical integration of the differential equations describing constrained systems. Based on Baumgarte method, we propose numerical methods for reducing the errors in the satisfaction of the constraints.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.45 no.3
• /
• pp.202-211
• /
• 2017

In this study, numerical simulations were performed using the ${\gamma}-Re_{\theta}$ transition model of KFLOW for the transitional flow over the KARI-11-180 airfoil. Numerical results of KFLOW were compared with experimental data and two other numerical results of XFoil and MSES. Fully turbulence model was predicted high skin friction drag than transition model because fully turbulence model could not solve the transitional flow. Numerical predictions using the ${\gamma}-Re_{\theta}$ model of KFLOW show a good agreement with experimental data and other numerical results. Present numerical results were confirmed the state of drag bucket due to dramatic changing of the transition location on the airfoil surface.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.5
• /
• pp.418-428
• /
• 2007

In the calculation of wave propagation by Boussinesq equation model, it is very common to experience numerical noises generated from nonlinear interaction and breaking wave occurrence, and the numerical solution is rapidly diverged unless the noises are properly controlled. A comparative study was here undertaken for the characteristics of three different lowpass filters (FFT filter, Gaussian filter and Shapiro filter) which are all designed to be applied to the interim results of numerical calculation. The numerical results obtained with application of respective filter techniques were compared with the results of an existing hydraulic experiment for the aspects of noise suppression, conservation of main signal and altering time. The results show that the Shapiro filter can be best applied with optimal choices of its element number, pass number and filtering tune interval. The combination of the number of filter element off, pass number of 50 or less, and application interval of 100 to 200 time steps generally showed good performance in both accuracy and efficiency of the numerical calculation.

• International Journal of Fluid Machinery and Systems
• /
• v.4 no.1
• /
• pp.199-208
• /
• 2011

This work presents numerical results, using OpenFOAM, of the flow in the swirl flow generator test rig developed at Politehnica University of Timisoara, Romania. The work shows results computed by solving the unsteady Reynolds Averaged Navier Stokes equations. The unsteady method couples the rotating and stationary parts using a sliding grid interface based on a GGI formulation. Turbulence is modeled using the standard k-${\varepsilon}$ model, and block structured wall function ICEM-Hexa meshes are used. The numerical results are validated against experimental LDV results, and against design velocity profiles. The investigation shows that OpenFOAM gives results that are comparable to the experimental and design profiles. The unsteady pressure fluctuations at four different positions in the draft tube is recorded. A Fourier analysis of the numerical results is compared whit that of the experimental values. The amplitude and frequency predicted by the numerical simulation are comparable to those given by the experimental results, though slightly over estimated.

• Steel and Composite Structures
• /
• v.12 no.1
• /
• pp.53-72
• /
• 2012

In this study, the physically-based failure models for matrix and fibers in compression and tension loading are introduced. For the 3D stress based fiber kinking model a modification is proposed for calculation of the fiber misalignment angle. All of these models are implemented into the finite element code by using the advantage of damage variable and the numerical results are discussed. To investigate the matrix failure model, purely in-plane transverse compression experiments are carried out on the specimens made by Glass/Epoxy to obtain the fracture surface angle and then a comparison is made with the calculated numerical results. Furthermore, shear failure of $({\pm}45)_s$ model is investigated and the obtained numerical results are discussed and compared with available experimental results. Some experiments are also carried out on the woven laminated composites to investigate the fracture pattern in the matrix failure mode and shown that the presented matrix failure model can be used for the woven composites. Finally, the obtained numerical results for stress based fiber kinking model and improved ones (strain based model) are discussed and compared with each other and with the available results. The results show that these models can predict the kink band angle approximately.

• Journal of computational fluids engineering
• /
• v.20 no.2
• /
• pp.46-51
• /
• 2015

In this paper, we have developed numerical model for particulated flow through narrow slit using Eulerian-Granular method. Commercial software (FLUENT) was utilized as simulation tool and main focus was to identify the effect from various numerical options for modeling of solid particles as continuos phase in granular flow. Gidaspow model was chosen as basic model for solid viscosity and drag model. And lun-et-al model was used as solid pressure and radial distribution model, respectively. Several other model options in FLUENT were tested considering the cross related effect. Mass flow rate of the particulate through the slit was measured to compare. Due to the high volume density of the stacked particulates above the slit, effect from various numerical options were not significant. The numerical results from basic model were also compared with experimental results and showed very good agreement. The effects from the characteristics of particles such as diameter, angle of internal friction, and collision coefficient were also analyzed for future design of velocity resistance layer in solar thermal absorber. Angle of internal friction was found to be the dominat variable for the particle mass flow rate considerably. More defined 3D model along with energy equation for complete solar thermal absorber design is currently underway.

• International Journal of Concrete Structures and Materials
• /
• v.10 no.4
• /
• pp.461-478
• /
• 2016

This paper aims to provide guidelines for the numerical modeling of reinforced concrete (RC) frame elements in order to assess the seismic performance of structures. Several types of numerical models RC frame elements are available in nonlinear structural analysis packages. Since these numerical models are formulated based on different assumption and theories, the models accuracy, computing time, and applicability vary, which poses a great difficulty to practicing engineering and limits their confidence in the analysis resultants. In this study, the applicability of four representative numerical models of RC frame elements is evaluated through comparison with experimental results of four-storey bare frame available from European Laboratory for Structural Assessment. The accuracy of a numerical model is evaluated according to the top displacement, interstorey drift, Maximum storey shear, damage pattern and energy dissipation capacity of the frame structure. The results obtained allow a better understanding of the characteristics and potentialities of all procedures, helping the user to choose the best approach to perform nonlinear analysis.

• Coupled systems mechanics
• /
• v.7 no.1
• /
• pp.95-109
• /
• 2018

The relevance of turbulent mixing in estuarine numerical models for stratified two-layer shallow water flows is analysed in this paper. A one-dimensional numerical model was developed for this purpose by extending an immiscible two-layer model with an additional source term, which accounts for turbulent mixing effects, namely the entrainment of fluid from the lower to the upper layer. The entrainment rate is quantified by an empirical equation as a function of the bulk Richardson number. A finite volume method based on an approximated Roe solver was used to solve the governing coupled system of partial differential equations. A comparison of numerical results with and without entrainment is presented to illustrate the influence of entrainment on both the salt-water intrusion length and lower layer dynamics. Furthermore, one example is given to demonstrate how entrainment terms may help to stabilize the numerical scheme and prevent a possible loss of hyperbolicity. Finally, the model with entrainment is validated by comparing the numerical results to field measurements.

• Journal of Korea Water Resources Association
• /
• v.31 no.2
• /
• pp.189-198
• /
• 1998

Five dam-break floodwave models are t재 field data sets. The models included FLDWAV, SMPDBK, HEC-1, Tr66, and HEC Dimensionless Graph. The field data sets documented the disasters at Teton dam, and Yeunchun dam. The FLDWAV results are uesd to test the sensitivity of the floodwave to variations in Manning's roughness coefficient, breach size, and bottom slope. The HEC-1 analysis includes testing the sensitivity of the results to model parameters. The TR66 model and FLDWAV, with channel routing by TR66 in both cases. SMPDBK and the Dimensionless Graph procedure are applied without particular difficulties being encountered in both real world cases. It is necessary to analyze numerical limit of existing numerical models and then to apply the relatively accurate numerical model in real basin. It is found that FLDWAV model is superior in numerical accuracy and stability to any other model. This study will contribute to improve defect of numerical models and develop more precise numerical model for a efficient and rapid dam breach flood disaster predict.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2002.10a
• /
• pp.191-196
• /
• 2002

The paper presents a numerical and experimental investigation on dynamic behaviors of a towed low tension cable. In the numerical study, an implicit finite difference algorithm is employed for three-dimensional cable equations. Fluid and geometric non-linearity and bending stiffness are considered and solved by Newton-Raphson iteration. Block tri-diagonal matrix method is applied for the fast calculation of the huge size of matrices. In order to verify the numerical results and to see real physical phenomena, an experiment is carried out for a 6m cable in a deep and long towing tank. The cable is towed in two different ways; one is towed at a constant speed and the other is towed at a constant speed with top end horizontal oscillations. Cable tension and shear forces are measured at the top end. Numerical and experimental results are compared with good agreements in most cases but with some differences in a few cases. The differences are due to drag coefficients caused by vortex shedding. In the numerical modeling, non-uniform element length needs to be employed to cope with the sharp variation of tension and shear forces at near top end.