• Environmental Engineering Research
• /
• 제25권4호
• /
• pp.605-613
• /
• 2020

This study is to investigate the possibility of using activated carbon prepared from Iraqi date-pits (ADP) which are produced from palm trees (Phoenix dactylifera L.) as low-cost reactive material in the permeable reactive barrier (PRB) for treating lead (Pb⁺²) from the contaminated groundwater, and then compare the results experimentally with other common reactive materials such as commercial activated carbon (CAC), zeolite pellets (ZP). Factors influencing sorption such as contact time, initial pH of the solution, sorbent dosage, agitation speed, and initial lead concentration has been studied. Two isotherm models were used for the description of sorption data (Langmuir and Freundlich). The maximum lead sorption capacities were measured for ADP, CAC, and ZP and were found to be 24.5, 12.125, and 4.45 mg/g, respectively. The kinetic data were analyzed using various kinetic models particularly pseudo-first-order, pseudo-second-order, and intraparticle diffusion. COMSOL Multiphysics 3.5a depend on finite element procedure was applied to formulate transmit of lead (Pb⁺²) in the two-dimensional numerical (2D) model under an equilibrium condition. The numerical solution shows that the contaminant plume is hindered by PRB.

• Coupled systems mechanics
• /
• 제7권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
• /
• 제19권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.

• 지구물리와물리탐사
• /
• 제21권1호
• /
• pp.1-7
• /
• 2018

전기탐사 2차원 모델링에서는 다수의 파수영역 전위를 계산하고 이를 푸리에 역변환하여 공간영역 전위를 계산한다. 푸리에 역변환은 여러 개의 서로 다른 파수에서의 파수영역 전위를 사용하여 수치적으로 얻어진다. 적분의 정확도를 향상시키기 위하여 파수의 크기에 따라 적분 구간을 지수 근사와 대수 근사 구간으로 분할하는 방법이 널리 사용되고 있다. 푸리에 역변환에는 크게 구간 적분법과 가우스 적분법이 사용되고 있다. 그러나 이들 방법은 송수신 간격을 고려하지 못하므로 송수신 간격에 따른 오차를 피할 수 없다. 특히 송수신 간격이 매우 작거나 클 경우 오차가 급격하게 증가하는 문제점을 가지고 있다. 이 연구에서는 송수신 간격을 고려하여 가우스 좌표값 및 가중값을 적용하는 새로운 수치 적분법을 개발하였다. 반무한 공간에 대한 수치 실험 결과, 개발된 수치 적분법은 송수신 간격에 관계없이 0.4% 이하의 정밀도를 나타내었다.

• Structural Engineering and Mechanics
• /
• 제32권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
• /
• 제26권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.

• 상하수도학회지
• /
• 제21권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.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2004년도 제22회 춘계학술대회논문집
• /
• 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년도 제32회 KOSCO SYMPOSIUM 논문집
• /
• 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
• /
• 제5권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.