• Geomechanics and Engineering
• /
• v.28 no.5
• /
• pp.437-454
• /
• 2022

The hydro-mechanical (HM) two-way sequential coupling in the TOUGH2-FLAC3D linking algorithm is validated completely and successfully in both M to H and H to M directions, which are initiated by mechanical surface loading for geomechanical validation and hydrological groundwater pumping for hydrogeological validation, respectively. For such complete and successful validation, a TOUGH2-FLAC3D linked numerical model is developed first by adopting the TOUGH2-FLAC3D linking algorithm, which uses the two-way (fixed-stress split) sequential coupling scheme and the implicit backward time stepping method. Two geomechanical and two hydrogeological validation problems are then simulated using the linked numerical model together with basic validation strategies and prerequisites. The second geomechanical and second hydrogeological validation problems are also associated with the Mandel effect and the Noordbergum and Rhade effects, respectively, which are three phenomenally well-known but numerically challenging HM effects. Finally, sequentially coupled numerical solutions are compared with either analytical solutions (verification) or fully coupled numerical solutions (benchmarking). In all the four validation problems, they show almost perfect to extremely or very good agreement. In addition, the second geomechanical validation problem clearly displays the Mandel effect and suggests a proper or minimum geometrical ratio of the height to the width for the rectangular domain to maximize agreement between the numerical and analytical solutions. In the meantime, the second hydrogeological validation problem clearly displays the Noordbergum and Rhade effects and implies that the HM two-way sequential coupling scheme used in the linked numerical model is as rigorous as the HM two-way full coupling scheme used in a fully coupled numerical model.

• The Pure and Applied Mathematics
• /
• v.24 no.3
• /
• pp.147-153
• /
• 2017

In this paper, we propose a numerical method to solve fuzzy differential equations. Numerical experiments show that when the step size is small, the new method has significantly good approximate solutions of fuzzy differential equation. Graphical representation of fuzzy solutions in three-dimension is also provided as a reference of visual convergence of the solution sequence.

• 제어로봇시스템학회:학술대회논문집
• /
• 1994.10a
• /
• pp.721-724
• /
• 1994

This paper proposes a continuous hitting by a flexible link hammer. This hammer system is used only the first mode of vibration for a desired hitting. The input of the hammer driver for a continuous hitting is obtained from numerical solutions of two sets of non-linear simultaneous equations which satisfy the hitting conditions. Being too complicated, these numerical calculations are not useful for online processing. Therefore, the multi-layered neural networks are applied to the generation of the input patterns of the hammer driver. The trained network outputs agree well to the numerical solutions.

• Water for future
• /
• v.28 no.6
• /
• pp.159-168
• /
• 1995

The run-up and the evolution of solitary waves on steep beaches are investigated by using a two-dimensional boundary integral equation model. The model is first used to compute the run-up heights of solitary waves on a relatively mind slope. The model is verified by comparing the computed numerical solutions with available experimental data, other numerical solutions and approximated analytical solutions. The agreement between the present numerical solutions and the other data is found to be excellent. The model is then applied to the calculation of run-up heights on very steep slopes. As far as the maximum run-up of solitary waves is concerned, the boundary integral equation model provides reasonable and reliable solutions. Finally, the evolution on steep beaches is also examined and the obtained wave heights are compared with those calculated from the Green's law.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.21 no.2
• /
• pp.89-107
• /
• 2017

In this paper, we have proposed a modified Marker-And-Cell (MAC) method to investigate the problem of an unsteady 2-D incompressible flow with heat and mass transfer at low, moderate, and high Reynolds numbers with no-slip and slip boundary conditions. We have used this method to solve the governing equations along with the boundary conditions and thereby to compute the flow variables, viz. u-velocity, v-velocity, P, T, and C. We have used the staggered grid approach of this method to discretize the governing equations of the problem. A modified MAC algorithm was proposed and used to compute the numerical solutions of the flow variables for Reynolds numbers Re = 10, 500, and 50000 in consonance with low, moderate, and high Reynolds numbers. We have also used appropriate Prandtl (Pr) and Schmidt (Sc) numbers in consistence with relevancy of the physical problem considered. We have executed this modified MAC algorithm with the aid of a computer program developed and run in C compiler. We have also computed numerical solutions of local Nusselt (Nu) and Sherwood (Sh) numbers along the horizontal line through the geometric center at low, moderate, and high Reynolds numbers for fixed Pr = 6.62 and Sc = 340 for two grid systems at time t = 0.0001s. Our numerical solutions for u and v velocities along the vertical and horizontal line through the geometric center of the square cavity for Re = 100 has been compared with benchmark solutions available in the literature and it has been found that they are in good agreement. The present numerical results indicate that, as we move along the horizontal line through the geometric center of the domain, we observed that, the heat and mass transfer decreases up to the geometric center. It, then, increases symmetrically.

• Journal of the Korean Statistical Society
• /
• v.35 no.2
• /
• pp.143-156
• /
• 2006

In the present paper, various solutions for generalized canonical correlation analysis (GCCA) are considered depending on the criteria and constraints. For the comparisons of some characteristics of the solutions, we provide with certain unifying stationary equations which might to also useful to obtain various generalized canonical correlation analysis solutions. In addition, we suggest an approach for the generalized canonical correlation analysis by exploiting the concept of maximum eccentricity originally de-signed to test the internal independence structure. The solutions, including new one, are compared through unifying stationary equations and by using some numerical illustrations. A type of iterative procedure for the GCCA solutions is suggested and some numerical examples are provided to illustrate several GCCA methods.

• Computational Structural Engineering : An International Journal
• /
• v.1 no.1
• /
• pp.49-57
• /
• 2001

This paper presents exact solutions for the modal stress-resultant distributions for vibrating rectangular Mindlin plates involving two opposite sides simply supported while the other two sides free. These exact stress-resultants of vibrating plates with free edges, hitherto unavailable, are very important because they serve as benchmark solutions for checking numerical solutions and methods. Using the exact solutions of a square plate, this paper highlights the problem of determining accurate stress-resultants, especially the transverse shear forces and twisting moments in thin plates, when employing the widely used numerical methods such as the Ritz method and the finite element method. Thus, this study shows that there is a need for researchers to develop refinements to the Ritz method and the finite element method for determining very accurate stress-resultants in vibrating plates with free edges.

• International Journal of Highway Engineering
• /
• v.18 no.3
• /
• pp.21-32
• /
• 2016

PURPOSES : In this paper, the analytical solutions suggested to simulate the behavior of rheological fluids were rigorously re-derived and investigated for fixed conditions to evaluate the applicability for the solutions on a mini-cone slump test of cement paste. The selected solutions with proper boundary conditions can be used as reference solutions to evaluate the performance of numerical simulation approaches, such as the discrete element method. METHODS : The slump, height, and spread radius for the given boundary and yield stress conditions that are determined by five different analytical solutions are compared. RESULTS : The analytical solution based on fluid mechanics for pure shear flow shows similar results to that for intermediate flow at low yield stresses. The fluid mechanics-based analytical solution resulted in a very similar trend to the geometry-based analytical solution. However, it showed a higher slump at high yield stress and lower slump at low yield stress ranges than the geometry-based analytical model. The analytical solution based on the mini-cone geometry was not significantly affected by the yield criteria, such as von Mises and Tresca. CONCLUSIONS : Even though differences among the analytical solutions in terms of slump and spread radius existed, the difference can be considered insignificant when the solutions were used as reference to evaluate the appropriateness of numerical approaches, such as the discrete element method.

• Structural Engineering and Mechanics
• /
• v.86 no.1
• /
• pp.69-83
• /
• 2023

Analytical solutions to problems are crucial because they provide high-quality comparison data for assessing the accuracy of numerical solutions. Benchmark analytical solutions for the vibrations of cracked functionally graded material (FGM) plates are not available in the literature because of the high level of complexity of such solutions. On the basis of first-order shear deformation plate theory (FSDT), this study proposes analytical series solutions for the vibrations of FGM rectangular plates with side or internal cracks parallel to an edge of the plates by using Fourier cosine series and the domain decomposition technique. The distributions of FGM properties along the thickness direction are assumed to follow a simple power law. The proposed analytical series solutions are validated by performing comprehensive convergence studies on the vibration frequencies of cracked square plates with various crack lengths and under various boundary condition combinations and by performing comparisons with published results based on various plate theories and the theory of three-dimensional elasticity. The results reveal that the proposed solutions are in excellent agreement with literature results obtained using the Ritz method on the basis of FSDT. The paper also presents tabulations of the first six nondimensional frequencies of cracked rectangular Al/Al2O3 FGM plates with various aspect ratios, thickness-to-width ratios, crack lengths, and FGM power law indices under six boundary condition combinations, the tabulated frequencies can serve as benchmark data for assessing the accuracy of numerical approaches based on FSDT.

• Journal of Korea Water Resources Association
• /
• v.30 no.6
• /
• pp.753-759
• /
• 1997

A new numerical model was doveloped to simulate density-dependent ground water flow and solute transport. Accuracy of a numerical model depends upon how well it simulates advection dominant situations because numerical oscillations can spoil solutions for these situations. Nonlinear oscillation-absorption finite element method. based on the variational principle, was employed. Unlike previous numerical models, this model can easily be expanded for more complex situations. Accuracy of the model is evaluated by comparing with analytical solutions and results of other numerical model.