• Journal of Advanced Marine Engineering and Technology
• /
• v.28 no.5
• /
• pp.818-826
• /
• 2004

The spray-induced mixing characteristics and thermal decomposition of aqueous urea solution into ammonia have been studied to design optimum sizes and geometries of the mixing chamber in SCR(Selective Catalytic Reduction) system. The cold flow tests about the urea-injection nozzle were performed to clarify the parameters of spray mixing characteristics such as mean diameter and velocity of drops and spray width determined from the interactions between incoming air and injected drops. Discrete particle model in Fluent code was adopted to simulate spray-induced mixing process and the experimental results on the spray characteristics were used as input data of numerical calculations. The simulation results on the spray-induced mixing were verified by comparing the spray width extracted from the digital images with the simulated Particle tracks of injected drops. The single kinetic model was adopted to predict thermal decomposition of urea solution into ammonia and solved simultaneously along with the verified spray model. The hot air generator was designed to match the flow rate and temperature of the exhaust gas of the real engines The measured ammonia productions in the hot air generator were compared with the numerical predictions and the comparison results showed good agreements. Finally, we concluded that the design capabilities for sizing optimum mixing chamber were established.

• Journal of Soil and Groundwater Environment
• /
• v.7 no.4
• /
• pp.10-16
• /
• 2002

The multi-aquifer system separated by semipervious leaky beds was analyzed. The finite difference scheme of the Crank-Nicolson method is applied to obtain the solution for this system. The solution of this scheme was compared with the analytical solution for two-layer aquifer systems with one-dimensional steady state. The results showed a good agreement between analytical and numerical solution for two-layer aquifer systems. So, the numerical scheme can be extended to multi-aquifer system. When the pumping is tried for single or multi aquifer, the computation of the groundwater heads was possible for each aquifer in the multi-aquifer with two-dimensional system. So, it might be helpful for the effective groundwater management.

• Geomechanics and Engineering
• /
• v.10 no.2
• /
• pp.189-205
• /
• 2016

Soils are subjected to additional stresses due to the loads transferred by the foundations of the buildings. The distribution of stress in soil has great importance in geotechnical engineering projects such as stress, settlement and liquefaction analyses. The purpose of this study is to examine the shear stresses on horizontal plane below the rectangular foundations subjected to biaxial bending on an elastic soil. In this study, closed-form analytical solutions for shear stresses in x and y directions were obtained from Boussinesq's stress equations. The expressions of analytical solutions were simplified by defining the shear stress influence values ($I_1$, $I_2$, $I_3$), and solution charts were presented for obtaining these values. For some special loading conditions, the expressions for shear stresses in the soil below the corners of a rectangular foundation were also given. In addition, a computer program was developed to calculate the shear stress increment at any point below the rectangular foundations. A numerical example for illustrating the use of the presented solution charts was given and, finally, shear stress isobars were obtained for the same example by a developed computer program. The shear stress expressions obtained in this work can be used to determine monotonic and cyclic behavior of soils below rectangular foundations subjected to biaxial bending.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.5
• /
• pp.439-446
• /
• 2012

This paper presents an implicit moving least squares(MLS) difference method for improving the solution accuracy of 1-D free boundary problems, which implicitly updates the topology change of moving interface. The conventional MLS difference method explicitly updates the moving interface; it requires no iterative solution procedure but results in the loss of accuracy. However, the newly developed implicit scheme makes the total system nonlinear involving iterative solution procedure, but numerical verification show that it dramatically elevates the solution accuracy with moderate computation increase. Through numerical experiments for melting problems having moving singularity, it is verified that the proposed method can achieve the second order accuracy.

• Advances in materials Research
• /
• v.6 no.4
• /
• pp.317-328
• /
• 2017

In this paper, an improved theoretical solution for interfacial stress analysis is presented for simply supported concrete beam bonded with a sandwich FGM plate. Interfacial stress analysis is presented for simply supported concrete beam bonded with a sandwich plate. This improved solution is intended for application to beams made of all kinds of materials bonded with a thin plate, while all existing solutions have been developed focusing on the strengthening of reinforced concrete beams, which allowed the omission of certain terms. It is shown that both the normal and shear stresses at the interface are influenced by the material and geometry parameters of the composite beam. A numerical parametric study was performed for different simulated cases to assess the effect of several parameters. Numerical comparisons between the existing solutions and the present new solution enable a clear appreciation of the effects of various parameters. The results of this study indicated that the FGM sandwich panel strengthening systems are effective in enhancing flexural behavior of the strengthened RC beams.

• Coupled systems mechanics
• /
• v.11 no.2
• /
• pp.167-198
• /
• 2022

In this paper we deal with classical instability problems of heterogeneous Euler beam under conservative loading. It is chosen as the model problem to systematically present several possible solution methods from simplest deterministic to more complex stochastic approach, both of which that can handle more complex engineering problems. We first present classical analytic solution along with rigorous definition of the classical Euler buckling problem starting from homogeneous beam with either simplified linearized theory or the most general geometrically exact beam theory. We then present the numerical solution to this problem by using reduced model constructed by discrete approximation based upon the weak form of the instability problem featuring von Karman (virtual) strain combined with the finite element method. We explain how such numerical approach can easily be adapted to solving instability problems much more complex than classical Euler's beam and in particular for heterogeneous beam, where analytic solution is not readily available. We finally present the stochastic approach making use of the Duffing oscillator, as the corresponding reduced model for heterogeneous Euler's beam within the dynamics framework. We show that such an approach allows computing probability density function quantifying all possible solutions to this instability problem. We conclude that increased computational cost of the stochastic framework is more than compensated by its ability to take into account beam material heterogeneities described in terms of fast oscillating stochastic process, which is typical of time evolution of internal variables describing plasticity and damage.

• Geomechanics and Engineering
• /
• v.28 no.3
• /
• pp.225-237
• /
• 2022

Settlement evaluation is important for shallow tunnels in big cities to estimate the settlement that occurs due to the excavation of twin tunnels. The majority of earlier research on analytical solutions, on the other hand, concentrated on calculating the settlement for a single tunnel. This research introduces a procedure to evaluate the settlement induced by the excavation of twin tunnels (two parallel tunnels). In this study, a series of numerical analysis were performed to validate the analytical solution results. Two geological conditions were considered to derive the settlement depending on each case. The analytical and numerical methods were compared, which involved considering many sections and conducting a parametric study; the results have good agreement. Moreover, a comparison of the 3D flat model and 2D (FEM) with the analytical solution shows that in the fill soil, the maximum settlement values were obtained by the analytical solution. In contrast, the values obtained by the analytical solution in the rock is more conservative than those in the fill. Finally, this method was shown to be appropriate for twin tunnels dug side by side by utilizing finite element analysis 3D and 2D (PLAXIS 3D and PLAXIS 2D) to verify the analytical equations. Eventually, it will be possible to use this approach to predict settlement troughs over twin tunnels.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2023.05a
• /
• pp.171-171
• /
• 2023

This study investigates the capability of Physics-Informed Neural Networks (PINNs) for solving the solution of partial differential equations. Particularly, the 1D Saint-Venant Equations (SVEs) were considered, which describe the movement of water in a domain with shallow depth compared to its horizontal extent, and are widely adopted in hydrodynamics, river, and coastal engineering. The core contribution of this work is to combine the robustness of neural networks with the physical constraints of the SVEs. The PINNs method utilized a neural network to approximate the solutions of SVEs, while also enforcing the underlying physical principles of the equations. This allows for a more effective and reliable solution, especially in areas with complex geometry and varying bathymetry. To validate the robustness of the PINNs method, numerical experiments were conducted on several benchmark problems. The results show that the PINNs could be achieved high accuracy when compared with the solution from the numerical solution. Overall, this study demonstrates the potential of using PINNs and highlights the benefits of integrating neural network and physics information for improved efficiency and accuracy in solving SVEs.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.136-141
• /
• 2001

In this study, we analysed tree surface flow by using the experimental and numerical method with a different surfactant concentration. We compared numerical solution with experimental results for one-dimensional model. The result shows that in general the tree surface velocity can well be reproduced by the one-dimensional model for various surfactant concentration.

• Water for future
• /
• v.17 no.3
• /
• pp.211-219
• /
• 1984

A 4-level baroclinic numerical model is designed by using the vorticity equation and Omega equation. Block-Cyclic-Reduction method is applied to the solution of the Helmholtz defferential equation, which is proved to be better than the Relaxation method from the composite viewpoint of accuracy, stability and economy. It was investigated whether the model explains the physical process influenced by voricity and temperature advection. It was also examined if the model atmosphere describes the general circulation. This examination is similar to Phillips(1956). The result of this numerical experiment shows that the model explains qualitatively the Quasi-Geostrophic theory for the development of Baroclinic wave, as throughly described in Holton(1972).