• 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 the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.2
• /
• pp.197-215
• /
• 2015

Modeling water flow in variably saturated, porous media is important in many branches of science and engineering. Highly nonlinear relationships between water content and hydraulic conductivity and soil-water pressure result in very steep wetting fronts causing numerical problems. These include poor efficiency when modeling water infiltration into very dry porous media, and numerical oscillation near a steep wetting front. A one-dimensional finite element formulation is developed for the numerical simulation of variably saturated flow systems. First order backward Euler implicit and second order Crank-Nicolson time discretization schemes are adopted as a solution strategy in this formulation based on Picard and Newton iterative techniques. Five examples are used to investigate the numerical performance of two approaches and the different factors are highlighted that can affect their convergence and efficiency. The first test case deals with sharp moisture front that infiltrates into the soil column. It shows the capability of providing a mass-conservative behavior. Saturated conditions are not developed in the second test case. Involving of dry initial condition and steep wetting front are the main numerical complexity of the third test example. Fourth test case is a rapid infiltration of water from the surface, followed by a period of redistribution of the water due to the dynamic boundary condition. The last one-dimensional test case involves flow into a layered soil with variable initial conditions. The numerical results indicate that the Crank-Nicolson scheme is inefficient compared to fully implicit backward Euler scheme for the layered soil problem but offers same accuracy for the other homogeneous soil cases.

• Geophysics and Geophysical Exploration
• /
• v.21 no.1
• /
• pp.1-7
• /
• 2018

In the two-dimensional (2D) modeling of electrical method, the potential in the space domain is reconstructed with the calculated potentials in the wavenumber domain using inverse Fourier transform. The inverse Fourier integral is numerically evaluated using the transformed potential at different wavenumbers. In order to improve the precision of the integration, either the logarithmic or exponential approximation has been used depending on the size of wavenumber. Two numerical methods have been generally used to evaluate the integral; interval integration and Gaussian quadrature. However, both methods do not consider the distance from the current source. Thus the resulting potential in the space domain shows some error. Especially when the distance from the current source is very small or large, the error increases abruptly and the evaluated potential becomes extremely unstable. In this study, we developed a new method to calculate the integral accurately by introducing the distance from the current source to the rescaled Gauss abscissa and weight. The numerical tests for homogeneous half-space model show that the developed method can yield the error level lower than 0.4 percent over the various distances from the current source.

• Structural Engineering and Mechanics
• /
• v.32 no.2
• /
• pp.251-267
• /
• 2009

Reinforced concrete (RC) structures consist of two different materials: concrete and steel bar. The stress transfer behaviour between the two materials through bond plays an important role in the load-carrying capacity of RC structures, especially when they subject to lateral load such as blast and seismic load. Therefore, bond and slip between concrete and reinforcement bar will affect the response of RC structures under such loads. However, in most numerical analyses of blast-induced structural responses, the perfect bond between concrete and steel bar is often assumed. The main reason is that it is very difficult to model bond slip in the commercial finite element software, especially in hydrodynamic codes. In the present study, a one-dimensional slide line contact model in LS-DYNA for modeling sliding of rebar along a string of concrete nodes is creatively used to model the bond slip between concrete and steel bars in RC structures. In order to model the bond slip accurately, a new approach to define the parameters of the one-dimensional slide line model from common pullout test data is proposed. Reliability and accuracy of the proposed approach and the one-dimensional slide line in modelling the bond slip between concrete and steel bar are demonstrated through comparison of numerical results and experimental data. A case study is then carried out to investigate the bond slip effect on numerical analysis of blast-induced responses of a RC column. Parametric studies are also conducted to investigate the effect of bond shear modulus, maximum elastic slip strain, and damage curve exponential coefficient on blast-induced response of RC columns. Finally, recommendations are given for modelling the bond slip in numerical analysis of blast-induced responses of RC columns.

• Computers and Concrete
• /
• v.26 no.2
• /
• pp.107-114
• /
• 2020

In this study, a two dimensional model of receding contact problem has been analyzed using finite element method (FEM) based software ANSYS and ABAQUS. For this aim finite element modeling of elastic layer and two homogeneous, isotropic and symmetrical elastic quarter planes pressed by means of a rigid circular punch has been presented. Mass forces and friction are neglected in the solution. Since the problem is examined for the plane state, the thickness along the z-axis direction is taken as a unit. In order to check the accuracy of the present models, the obtained results are compared with the available results of the open literature as well as the results of two software are compared using Root Mean Square Error (RMSE) and good agreements are found. Numerical analyses are performed considering different values of the external load, rigid circular radius, quarter planes span length and material properties. The contact lengths and contact stresses of these values are examined, and their results are presented. Consequently, it is concluded that the considered non-dimensional quantities have noteworthy influence on the contact lengths and contact stress distributions, additionally if FEM analysis is used correctly, it can be an efficient alternative method to the analytical solutions that need time.

• Journal of Korean Society of Water and Wastewater
• /
• v.21 no.1
• /
• pp.101-112
• /
• 2007

This study has used RMA2 model and RMA4 model, which are depth-averaged two-dimensional flow and water-quality prediction models, to analyze the variation of the water-quality by the gates operation in the Seonakdong River. Sensitivity analysis is performed to get the Manning coefficient and the coefficient of eddy viscosity for RMA2 model, and to get the diffusion coefficient for RMA4 model. Since the numerical simulation using RMA2 and RMA4 models did not consider tributary pollutant load except for that of Joman River, it could make a little difference from the natural phenomenon. Nevertheless, the numerical simulation shows that the discharge of $30m^3/s$, which is the continuous inflow from the Daedong-gate, can make it possible to achieve the target water-quality (BOD 4.3mg/L) of Nakbon-N watershed about 10 days later if the Daejeo-gate could remain opened in connection with the Noksan-gate operation.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.825-830
• /
• 2004

A numerical analysis based on two-dimensional and three-dimensional incompressible Navier-Stokes equations has been carried out for double-circular-arc (DCA) compressor cascades. Two types of double-circular-arc cascades were used in this analysis. The appropriate turbulence model for compressor analysis was selected among the conventional turbulence models such as Baldwin-Lomax, k-$\varepsilon$ and k-$\varepsilon$ models. The results of current study were compared with available experimental data at various incidence angles. The 2-D and 3-D computational codes based on SIMPLE/PWIM algorithm for collocated grid and hybrid scheme for the convective terms were the main features of numerical tools. As commonly known, turbulence modeling is very important for the prediction of cascade flows, which are extremely complex with separation and reattachment by adverse pressure gradient. For selection of turbulence model, 2-D analysis was performed. And then, k-$\varepsilon$ turbulence model with wall function chosen as the reasonable turbulence model for 3-D calculation was used to increase the efficiency of computation times. A reasonable result of 3-D flow pattern passing through the double-circular-arc cascade was obtained.

• 한국연소학회:학술대회논문집
• /
• 2006.04a
• /
• pp.49-54
• /
• 2006

The present study conducted a numerical modeling on the diesel SCR (selective catalytic reduction) system using ammonia as a reductant over vanadium-based catalysts $(V_2O_5-WO_3/TiO_2)$. Transient modeling for ammonia adsorption/desorption on the catalyst surface was firstly carried out, and then the SCR reaction was modeled considering for it. In the current catalytic reaction model, we extended the pure chemical kinetic model based on laboratory-scale powdered-phase catalyst experiments to the chemico-physical one applicable to realistic commercial SCR reactors. To simulate multi-dimensional heat and mass transfer phenomena, the SCR reactor was modeled in two dimensional, axisymmetric domain using porous medium approach. Also, since diesel engines operate in transient mode, the present study employed an unsteady model. In addition, throughout simulations using the developed code, effects of space velocity on the DeNOx performance were investigated.

• Earthquakes and Structures
• /
• v.5 no.1
• /
• pp.111-127
• /
• 2013

This paper describes effects of infill walls on behavior of RC frame with low strength, including numerical modeling, modal testing and finite-element model updating. For this purpose full scaled, one bay and one story RC frame is produced and tested for plane and brick in-filled conditions. Ambient-vibration testis applied to identify dynamic characteristics under natural excitations. Enhanced Frequency Domain Decomposition and Stochastic Subspace Identification methods are used to obtain experimental dynamic characteristics. A numerical modal analysis is performed on the developed two-dimensional finite element model of the frames using SAP2000 software to provide numerical frequencies and mode shapes. Dynamic characteristics obtained by numerical and experimental are compared with each other and finite element model of the frames are updated by changing some uncertain modeling parameters such as material properties and boundary conditions to reduce the differences between the results. At the end of the study, maximum differences in the natural frequencies are reduced on average from 34% to 9% and a good agreement is found between numerical and experimental dynamic characteristics after finite-element model updating. In addition, it is seen material properties are more effective parameters in the finite element model updating of plane frame. However, for brick in-filled frame changes in boundary conditions determine the model updating process.

• Geomechanics and Engineering
• /
• v.9 no.6
• /
• pp.775-791
• /
• 2015

Due to the extension of communication ways (metro, highways, railways), hence, to improve traffic flow imposes often the difficult crossing that generally drive to the construction of underground works (tunnel, water conveyance tunnel...) plays a major role in the redevelopment of urban areas. This study is focused on the assessment of the interaction response of parallel tunnels, so this study uses the results from the simulation of two tunnels to illustrate a few observations that may aid in practical designs. In this article, simultaneous drilling of highway's twin tunnels is simulated by means of Finite Element Method (FEM) implemented in Plaxis program. So the treated subject appears in a setting of geotechnical where one can be to construct several tunnels sometimes in a ground of weak mechanical characteristics. The objective of this study is to simulate numerically the interaction effects caused by construction of two parallels tunnels. This is an important factor in the study of the total answer of the problem interaction between parallels underground works. The importance of the effects transmitted is function of several parameters as the type of the works, and the mechanical characteristics (tunnel size, depth, and the relative position between two tunnels, lining thickness...). This article describes numerical analyses of two parallels tunnels interaction. This study will be applied to a real case of a section tunnel T4 of the highway East-West (Algeria); the study presented below comprises a series of numerical simulations of two tunnels using the computer program Plaxis which is used in the analyses is based on Finite Element Method.