• Journal of Mechanical Science and Technology
• /
• v.16 no.3
• /
• pp.395-402
• /
• 2002

This paper shows how to formulate the transient analysis of 2-dimensional Hagen-Poiseuille flow using smoothed particle hydrodynamics (SPH). Treatments of viscosity, particle approximation and boundary conditions are explained. Numerical tests are calculated to examine effects caused by the number of particles, the number of particles per smoothing length, artificial viscosity and time increments for 2-dimensional Hagen-Poiseuille flow. Artificial viscosity for reducing the numerical instability directly affects the velocity of the flow, though effects of the other parameters do not produce as much effect as artificial viscosity. Numerical solutions using SPH show close agreement with the exact ones for the model flow, but SPH parameter must be chosen carefully Numerical solutions indicate that SPH is also an effective method for the analysis of 2-dimensional Hagen-Poiseuille flow.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.04a
• /
• pp.153-158
• /
• 2005

Numerical simulations of two-dimensional steady incompressible lid-driven flow in a square cavity are presented to verify the validity of a new solution code(PowerCFD) with unstructured grids. The code uses the non-staggered(collocated) grid approach which is very popular for incompressible flow analysis because of its numerical efficiency on the curvilinear or unstructured grids. Solutions are obtained for configurations with a Reynolds number as high as 10,000 with both rectangular and hybrid types of unstructured grid mesh. Interesting features of the flow are presented in detail and comparisons are made with benchmark solutions found in the literature. It is found that the code is capable of producing accurately the nature of the lid-driven cavity flow at high Reynolds numbers.

• The KSFM Journal of Fluid Machinery
• /
• v.9 no.4 s.37
• /
• pp.27-35
• /
• 2006

A Navier-Stokes code has been developed to simulate the flow through a cross flow fan. It is based on an unstructured finite volume method and uses moving grid technique to model the rotation of the fan. A low Reynolds number turbulence model is used to calculate eddy viscosity. The basic algorithm is SIMPLE. Numerical simulations over a wide range of flow rate are carried out to validate the code. Comparison of all numerical solutions with experimental data confirms the validity of the present code. Present numerical solutions show a noticeable improvement over a previous numerical method which is based on a model of body force to simulate the rotation of the impeller.

• Honam Mathematical Journal
• /
• v.43 no.3
• /
• pp.373-383
• /
• 2021

In this study, a numerical method deals with the Chebyshev wavelet collocation and Adomian decomposition methods are proposed for solving Korteweg-de Vries equation. Integration of the Chebyshev wavelets operational matrices is derived. This problem is reduced to a system of non-linear algebraic equations by using their operational matrix. Thus, it becomes easier to solve KdV problem. The error estimation for the Chebyshev wavelet collocation method and ADM is investigated. The proposed method's validity and accuracy are demonstrated by numerical results. When the exact and approximate solutions are compared, for non-linear or linear partial differential equations, the Chebyshev wavelet collocation method is shown to be acceptable, efficient and accurate.

• Journal of the Korean Mathematical Society
• /
• v.58 no.4
• /
• pp.835-847
• /
• 2021

We study behavior of numerical solutions for a nonlinear eigenvalue problem on ℝⁿ that is reduced from a dispersion managed nonlinear Schrödinger equation. The solution operator of the free Schrödinger equation in the eigenvalue problem is implemented via the finite difference scheme, and the primary nonlinear eigenvalue problem is numerically solved via Picard iteration. Through numerical simulations, the results known only theoretically, for example the number of eigenpairs for one dimensional problem, are verified. Furthermore several new characteristics of the eigenpairs, including the existence of eigenpairs inherent in zero average dispersion two dimensional problem, are observed and analyzed.

• Coupled systems mechanics
• /
• v.4 no.1
• /
• pp.67-98
• /
• 2015

Numerical prediction of dynamic behavior of fully coupled saturated porous media is of great importance in many engineering problems. Specifically, static and dynamic response of soils - porous media with pores filled with fluid, such as air, water, etc. - can only be modeled properly using fully coupled approaches. Modeling and simulation of static and dynamic behavior of soils require significant Verification and Validation (V&V) procedures in order to build credibility and increase confidence in numerical results. By definition, Verification is essentially a mathematics issue and it provides evidence that the model is solved correctly, while Validation, being a physics issue, provides evidence that the right model is solved. This paper focuses on Verification procedure for fully coupled modeling and simulation of porous media. Therefore, a complete Solution Verification suite has been developed consisting of analytical solutions for both static and dynamic problems of porous media, in time domain. Verification for fully coupled modeling and simulation of porous media has been performed through comparison of the numerical solutions with the analytical ones. Modeling and simulation is based on the so called, u-p-U formulation. Of particular interest are numerical dispersion effects which determine the level of numerical accuracy. These effects are investigated in detail, in an effort to suggest a compromise between numerical error and computational cost.

• Communications for Statistical Applications and Methods
• /
• v.26 no.3
• /
• pp.305-313
• /
• 2019

Generalized canonical correlation analysis (GCCA) extends the canonical correlation analysis (CCA) to the case of more than two sets of variables and there have been many studies on how two-set canonical solutions can be generalized. In this paper, we derive certain stationary equations which can lead the higher-order solutions of several GCCA methods and suggest a type of iterative procedure to obtain the canonical coefficients. In addition, with some numerical examples we present the methods for graphical display, which are useful to interpret the GCCA results obtained.

• Geomechanics and Engineering
• /
• v.3 no.4
• /
• pp.291-321
• /
• 2011

The paper describes a simple numerical FLAC model that was developed to simulate the dynamic response of two instrumented reduced-scale model reinforced soil walls constructed on a 1-g shaking table. The models were 1 m high by 1.4 m wide by 2.4 m long and were constructed with a uniform size sand backfill, a polymeric geogrid reinforcement material with appropriately scaled stiffness, and a structural full-height rigid panel facing. The wall toe was constructed to simulate a perfectly hinged toe (i.e. toe allowed to rotate only) in one model and an idealized sliding toe (i.e. toe allowed to rotate and slide horizontally) in the other. Physical and numerical models were subjected to the same stepped amplitude sinusoidal base acceleration record. The material properties of the component materials (e.g. backfill and reinforcement) were determined from independent laboratory testing (reinforcement) and by back-fitting results of a numerical FLAC model for direct shear box testing to the corresponding physical test results. A simple elastic-plastic model with Mohr-Coulomb failure criterion for the sand was judged to give satisfactory agreement with measured wall results. The numerical results are also compared to closed-form solutions for reinforcement loads. In most cases predicted and closed-form solutions fall within the accuracy of measured loads based on ${\pm}1$ standard deviation applied to physical measurements. The paper summarizes important lessons learned and implications to the seismic design and performance of geosynthetic reinforced soil walls.

• Journal of Mechanical Science and Technology
• /
• v.15 no.6
• /
• pp.752-763
• /
• 2001

Although there are many different solution schemes proposed for multidimensional radiative transfer, reference solutions to benchmark these methods are very rare in the literature. In this paper we produced some accurate solutions for purely absorbing gray and nongray gases including H$_2$O and CO$_2$by using the discrete transfer method with sufficiently accurate T(sub)95 quadrature set. The spectral transmittances of the mixtures of H$_2$O and CO$_2$are estimated by using the narrow band model. The gray gas solutions are obtained for different absorption coefficients, and the nongray real gas solutions are obtained for different mixture fractions of H$_2$O and CO$_2$. The numerical solutions presented in this paper are proved to be sufficiently accurate as compared to the available exact solutions and they may be used as reference solutions in evaluating various solution schemes.

• Structural Engineering and Mechanics
• /
• v.64 no.1
• /
• pp.81-92
• /
• 2017

Many materials in engineering exhibit different modulus in tension and compression, which are known as bi-modulus materials. Based on the bi-modulus elastic theory, a modified semi-analytical model, by introducing a stress function, is established in this paper to study the mechanical response of a bi-modulus cylinder placed in an axisymmetric temperature field. Meanwhile, a numerical procedure to calculate the temperature stresses in bi-modulus structures is developed. It is proved that the bi-modulus solution can be degenerated to the classical same modulus solution, and is in great accordance with the solutions calculated by the semi-analytical model proposed by Kamiya (1977) and the numerical solutions calculated both by the procedure complied in this paper and by the finite element software ABAQUS, which demonstrates that the semi-analytical model and the numerical procedure are accurate and reliable. The result shows that the modified semi-analytical model simplifies the calculation process and improves the speed of computation. And the numerical procedure simplifies the modeling process and can be extended to study the stress field of bi-modulus structures with complex geometry and boundary conditions. Besides, the necessity to introduce the bi-modulus theory is discussed and some suggestions for the qualitative analysis and the quantitative calculation of such structure are proposed.