• Proceedings of the Korea Water Resources Association Conference
• /
• 2009.05a
• /
• pp.516-523
• /
• 2009

The artificial compressibility (AC) method for the incompressible Navier-Stokes equations in the generalized curvilinear coordinates using the primitive form is implemented. The main advantage of the AC approach is that the resulting system of equations resembles the system of compressible N-S equations and can thus be integrated in time using standard, well-established time-marching methods. The errors, which are the odd-even oscillation, for pressure field in using the artificial compressibility can be eliminated by using the $4^{th}$ order artificial dissipation term which is explicitly included. Even though this paper focuses exclusively on 2D laminar flows to validate and assess the performance of this solver, this numerical method is general enough so that it can be readily extended to carry out 3D URANS simulation of engineering flows. This algorithm yields practically identical velocity profiles and primary vortex and secondary vortices that are in excellent overall agreement with the results of the vorticity-stream function formulation (Ghia et al., 1982). However, the grid resolution have to be required to be large enough to express the various vortices.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2004.10a
• /
• pp.700-703
• /
• 2004

A reliability analysis and design procedure based on the design of experiment (DOE) is combined with the response surface method (RSM) for numerical efficiency. The procedure established is based on a 3$^n$ full factorial DOE for numerical quadrature using explicit formula of optimum levels and weights derived for general distributions. The full factorial moment method (FFMM) shows good performance in terms of accuracy and ability to treat non-normally distributed random variables. But, the FFMM becomes very inefficient because the number of function evaluation required increases exponentially as the number of random variables considered increases. To enhance the efficiency, the response surface moment method (RSMM) is proposed. In RSMM, experiments only with high probability are conducted and the rest of data are complemented by a quadratic response surface approximation without mixed terms. The response surface is updated by conducting experiments one by one until the value of failure probability is converged. It is calculated using the Pearson system and the four statistical moments obtained from the experimental data. A measure for checking the relative importance of an experimental point is proposed and named as influence index. During the update of response surface, mixed terms can be added into the formulation.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.6
• /
• pp.1163-1170
• /
• 1992

A computer program was developed to analyze the two-dimensional unsteady incompressible viscous flow behind a rectangular bluff body between two parallel plates. The Peaceman-Rachford alternating direction implicit numerical method and Wachspress parameter were adopted to solve the governing equations in vorticity-transport and stream function formulation. The steady state flow and the vortex flow behind a rectangular bluff body in a chemical were investigated for Reynolds numbers of 200 and 500. The vortex shedding was generated by a physical pertubation numerically imposed at the center of the flow field for a short time. It was observed that the perturbed flow became periodic after a transient period.

• Journal of Mechanical Science and Technology
• /
• v.16 no.6
• /
• pp.810-818
• /
• 2002

When the contribution of lightweight components to the total energy of a system is small, tole inertia effects are sometimes ignored by replacing them to massless links. For example, a revolute-spherical massless link generates two kinematic constraint equations between adjacent bodies and allows four relative degrees of freedom. In this paper, to implement a massless link systematically in a computer program using the velocity transformation technique, the velocity transformation matrix of massless links is derived and numerically implemented. The velocity transformation matrix for a revolute-spherical massless link and a revolute-universal massless link are appeared as a 6$\times$4 matrix and a 6$\times$3 matrix, respectively. A massless link model in a suspension composite joint transmitting external forces is also developed and the numerical efficiency of the proposed model is compared to a conventional multibody model. For a massless link transmitting external forces, forces acting on links are resolved and transmitted to the attached points with a quasi-static assumption. Numerical examples are presented to verify the formulation.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.27 no.2
• /
• pp.139-147
• /
• 2017

Acoustic performance of splitter silencers was investigated by using 3-dimensional commercial software and experiments. Flow resistivity of sound absorbing material was indirectly estimated by using an impedance tube setup and a curve fitting method. In addition the acoustic impedance of perforated plate was determined by an empirical formulation. Such properties have been used as input parameters in the commercial software. The prediction for a splitter silencer with 1000 mm length was compared with the experimental result. The numerical method is then applied to identify the effects of number of splitters, length of splitters, absorptive material density, and porosity of a perforated plate on the performance of the splitter silencers. As the number and length of splitter increases, the acoustic performance significantly increases. Although the increase of density of absorptive material also increase the acoustic performance, a change in the density over a certain level hardly affect it. The increase of porosity will enhance the performance especially at higher frequencies.

• Computers and Concrete
• /
• v.7 no.4
• /
• pp.365-386
• /
• 2010

This paper presents the physical formulation of a 1D material model suitable for seismic applications. It is written within the framework of thermodynamics with internal variables that is, especially, very efficient for the phenomenological representation of material behaviors at macroscale: those of the representative elementary volume. The model can reproduce the main characteristics observed for concrete, that is nonsymetric loading rate-dependent (viscoelasticity) behavior with appearance of permanent deformations and local hysteresis (continuum plasticity), stiffness degradation (continuum damage), cracking due to displacement localization (discrete plasticity or damage). The parameters have a clear physical meaning and can thus be easily identified. Although this point is not detailed in the paper, this material model is developed to be implemented in a finite element computer program. Therefore, for the benefit of the robustness of the numerical implementation, (i) linear state equations (no local iteration required) are defined whenever possible and (ii) the conditions in which the presented model can enter the generalized standard materials class - whose elements benefit from good global and local stability properties - are clearly established. To illustrate the capabilities of this model - among them for Earthquake Engineering applications - results of some numerical applications are presented.

• Structural Engineering and Mechanics
• /
• v.45 no.1
• /
• pp.95-110
• /
• 2013

In the recent decade, practical of wavelet technique is being utilized in various domain of science. Particularly, engineers are interested to the wavelet solution method in the time series analysis. Fundamentally, seismic responses of structures against time history loading such as an earthquake, illustrates optimum capability of systems. In this paper, a procedure using particularly discrete Haar wavelet basis functions is introduced, to solve dynamic equation of motion. In the proposed approach, a straightforward formulation in a fluent manner is derived from the approximation of the displacements. For this purpose, Haar operational matrix is derived and applied in the dynamic analysis. It's free-scaled matrix converts differential equation of motion to the algebraic equations. It is shown that accuracy of dynamic responses relies on, access of load in the first step, before piecewise analysis added to the technique of equation solver in the last step for large scale of wavelet. To demonstrate the effectiveness of this scheme, improved formulations are extended to the linear and nonlinear structural dynamic analysis. The validity and effectiveness of the developed method is verified with three examples. The results were compared with those from the numerical methods such as Duhamel integration, Runge-Kutta and Wilson-${\theta}$ method.

• Steel and Composite Structures
• /
• v.19 no.4
• /
• pp.1011-1033
• /
• 2015

In this article, large amplitude bending behaviour of laminated composite flat panel under combined effect of moisture, temperature and mechanical loading is investigated. The laminated composite panel model has been developed mathematically by introducing the geometrical nonlinearity in Green-Lagrange sense in the framework of higher-order shear deformation theory. The present study includes the degraded composite material properties at elevated temperature and moisture concentration. In order to achieve any general case, all the nonlinear higher order terms have been included in the present formulation and the material property variations are introduced through the micromechanical model. The nonlinear governing equation is obtained using the variational principle and discretised using finite element steps. The convergence behaviour of the present numerical model has been checked. The present proposed model has been validated by comparing the responses with those available published results. Some new numerical examples have been solved to show the effect of various parameters on the bending behaviour of laminated composite flat panel under hygro-thermo-mechanical loading.

• Structural Engineering and Mechanics
• /
• v.51 no.6
• /
• pp.923-937
• /
• 2014

The rational finite element method is different from the standard finite element method, which is constructed using basic solutions of the governing differential equations as interpolation functions in the elements. Therefore, it is superior to the isoparametric approach because of its obvious physical meaning and accuracy; it has successfully been applied to the isotropic elasticity problem. In this paper, the formulation of rational finite elements for plane orthotropic elasticity problems is deduced. This method is formulated directly in the physical domain with full consideration of the requirements of the patch test. Based on the number of element nodes and the interpolation functions, different approaches are applied with complete polynomial interpolation functions. Then, two special stiffness matrixes of elements with four and five nodes are deduced as a representative application. In addition, some typical numerical examples are considered to evaluate the performance of the elements. The numerical results demonstrate that the present method has a high level of accuracy and is an effective technique for solving plane orthotropic elasticity problems.

• Transactions of Materials Processing
• /
• v.25 no.1
• /
• pp.49-55
• /
• 2016

Electrohydraulic forming (EHF) is a high-speed forming process that uses an electric arc discharge in water. Shock waves resulting from the electric arc discharge are propagated to the blank through water and the blank moves toward the die. Advantages of EHF include improved formability due to the high-speed process and reduction of the bouncing effect. In the current study, a numerical simulation of EHF was developed using LS-DYNA. In the simulation, the model for the electric arc was assumed as an adiabatic gas expansion and an Arbitrary Lagrange-Eulerian (ALE) multi material formulation was used to describe the interaction between the electric arc and the water. In order to model the Fluid-Structure Interaction (FSI), a coupling mechanism was used. The blank of Al 1100-O was simulated using shell elements. The results of the simulation showed that the blank was deformed due to the pressure propagation of water and the bouncing effect did not affect the formability of blank.