• Journal of the Korean Society of Groundwater Environment
• /
• v.6 no.4
• /
• pp.188-193
• /
• 1999

Transport of pollutants in aquifer largely depends on groundwater flow which is governed by aquifer hydraulic parameters. Determination of these parameters and associated groundwater modeling become essential for adequate remediation of contaminate groundwater. The objective of this paper is to analyze groundwater flow and determine the optimum hydraulic parameters by performing groundwater modeling based on sensitivity analysis for unconfined sandy gavel aquifer constructed in a laboratory scale under various boundary condition. Results revealed that the simulated drawdown was lower than the observed drawdown irrespective of boundary conditions. and specific yield (S$_{y}$) had less effect on the grondwater flow than permeability (K) in the aquifer. Water balance analysis showed that the measured drawdown in neighboring observation wells during pumping was higher than either simulated or recovered water table. The indicated that a difference might exist in the water tables between aquifer and wells. The difference was investigated by time domain reflectometry (TDR) measurements on water contents in the region of water table and capillary fringe, and explained by a delayed response of water table during gravitational drainage as the water table was lowered as a result of pumping.g.

• International Journal of Vascular Biomedical Engineering
• /
• v.2 no.1
• /
• pp.17-24
• /
• 2004

Tumor angiogenesis was simulated using a two-dimensional computational model. The equation that governed angiogenesis comprised a tumor angiogenesis factor (TAF) conservation equation in time and space, which was solved numerically using the Galerkin finite element method. The time derivative in the equation was approximated by a forward Euler scheme. A stochastic process model was used to simulate vessel formation and vessel elongation towards a paracrine site, i.e., tumor-secreted basic fibroblast growth factor (bFGF). In this study, we assumed a two-dimensional model that represented a thin (1.0 mm) slice of the tumor. The growth of the tumor over time was modeled according to the dynamic value of bFGF secreted within the tumor. The data used for the model were based on a previously reported model of a brain tumor in which four distinct stages (namely multicellular spherical, first detectable lesion, diagnosis, and death of the virtual patient) were modeled. In our study, computation was not continued beyond the 'diagnosis' time point to avoid the computational complexity of analyzing numerous vascular branches. The numerical solutions revealed that no bFGF remained within the region in which vessels developed, owing to the uptake of bFGF by endothelial cells. Consequently, a sharp, declining gradient of bFGF existed near the surface of the tumor. The vascular architecture developed numerous branches close to the tumor surface (the brush-border effect). Asymmetrical tumor growth was associated with a greater degree of branching at the tumor surface.

• Journal of Energy Engineering
• /
• v.2 no.3
• /
• pp.268-275
• /
• 1993

A series of parametric investigations are performed in order to resolve the flow characteristic of a Quebec city stoker incinerator. The parameters considered in this study are five internal configurations of the Quebec city stoker itself and its modified ones, primary air velocity, the injection velocity and angle of the secondary air, and the reduction of the stoker exit area. A control-volume based finite-difference method by Patankar together with the power-law scheme is employed for discretization. The resolution of the pressure-velocity coupling is made by the use of SIMPLEC algorithm. The standard, two equation, k-$\varepsilon$ model is incorporated for the closure of turbulence. The size of recirculation region, turbulent viscosity, the mass fraction of the secondary air and pressure drop are calculated in order to analyze the characteristics of flow field. The results are physically acceptable and discussed in detail. The flow field of the Quebec city stoker shows the strong recirculation zone together with the high turbulence intensity over the upper part of the incinerator.

• Advances in environmental research
• /
• v.3 no.1
• /
• pp.45-69
• /
• 2014

The groundwater has been a major source of water supply throughout the ages. Around 50% of the rural as well as urban population in the developing countries like India depends on groundwater for drinking. The groundwater is also an important source in the agriculture and industrial sector. In many parts of the world, groundwater resources are under increasing threat from growing demands, wasteful use and contamination. A good planning and management practices are needed to face this challenge. A key to the management of groundwater is the ability to model the movement of fluids and contaminants in the subsurface environment. It is obvious that the contaminant source activities cannot be completely eliminated and perhaps our water bodies will continue to serve as receptors of vast quantities of waste. In such a scenario, the goal of water quality protection efforts must necessarily be the control and management of these sources to ensure that released pollutants will be sufficiently attenuated within the region of interest and the quality of water at points of withdrawal is not impaired. In order to understand the behaviour of contaminant transport through different types of media, several researchers are carrying out experimental investigations through laboratory and field studies. Many of them are working on the analytical and numerical studies to simulate the movement of contaminants in soil and groundwater of the contaminant transport. With the advent of high power computers especially, a numerical modelling has gained popularity and is indeed of particular relevance in this regard. This paper provides the state of the art of contaminant transport and reviews the allied research works carried out through experimental investigation or using the analytical solution and numerical method. The review involves the investigation in respect of both, saturated and unsaturated, porous media.

• Journal of the Society of Naval Architects of Korea
• /
• v.54 no.1
• /
• pp.43-48
• /
• 2017

In recent years, the demand for ships and offshore platforms that can navigate and operate through the Arctic Ocean has been rapidly increasing due to global warming and large reservoirs of oil and natural gas in the area. Winterization design is one of the key issues to consider in the robust structural safety design and building of ships that operate in the Arctic and Sub-Arctic regions. However, international regulations for winterization design in Arctic condition regulated that only those ships and offshore platforms with a Polar Class designation and/or an alternative standard. In order to cope with the rising demand for operating in the Arctic region, existing and new Arctic vessels with a Polar Class designation are lacking to cover for adequate winterization design with HSE philosophy. Existing ships and offshore platform was not designed based on reliable data based on numerical and experiment studies. There are only designed as a performance and functional purposes. It is very important to obtain of reliable data and provide of design guidance of the anti-icing structures by taking the effects of low temperature into consideration. Therefore, the main objective of this paper reconsiders anti-icing design of aluminum helideck using the heating cable. To evaluate of reliable data and recommend of anti-icing design method, various types of analysis and methods can be applied in general. In the present study, finite element method carried out the thermal analysis with cold chamber testing for performance and capacity of heating cables.

• 제어로봇시스템학회:학술대회논문집
• /
• 2005.06a
• /
• pp.700-705
• /
• 2005

The purpose of this project is the derivation and development of techniques for the new estimation of robustness for the systems having uncertainties. The basic ideas to analyze the system which is the originally nonlinear is Lyapunov direct theorems. The nonlinear systems have various forms of terms inside the system equations and this investigation is confined in the form of bounded uncertainties. Bounded means the uncertainties are with same positive/negative range. The number of uncertainties will be the degree of freedoms in the calculation of the stability region. This is so called the robustness bounds. This proposition adopts the theoretical analysis of the Lyapunov direct methods, that is, the sign properties of the Lyapunov function derivative integrated along finite intervals of time, in place of the original method of the sign properties of the time derivative of the Lyapunov function itself. This is the new sufficient criteria to relax the stability condition and is used to generate techniques for the robust design of control systems with structured perturbations. Using this relaxing stability conditions, the selection of Lyapunov candidate function is of various forms. In this paper, the quadratic form is selected. this generated techniques has been demonstrated by recent research interest in the area of robust control design and confirms that estimation of robustness bounds will be improved upon those obtained by results of the original Lyapunov method. In this paper, the symbolic algebraic procedures are utilized and the calculating errors are reduced in the numerical procedures. The application of numerical procedures can prove the improvements in estimations of robustness for one-and more structured perturbations. The applicable systems is assumed to be linear with time-varying with nonlinear bounded perturbations. This new techniques will be extended to other nonlinear systems with various forms of uncertainties, especially in the nonlinear case of the unstructured perturbations and also with various control method.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.19 no.1
• /
• pp.9-17
• /
• 2015

The basic flow configuration is composed of a plane, double shear layer where relatively thin mid gas layer is sandwiched between air and fuel stream. The present study describes numerical investigations concerning the combustion enhancement according to a variation of mid layer thickness. In this case, the effect of heat release in turbulent mixing layers is important. For the numerical solution, a fully conservative unsteady $2^{nd}$ order time accurate sub-iteration method and $2^{nd}$ order TVD scheme are used with the finite volume method including k-${\omega}$ SST model. The results consists of three categories; single shear layer consists of fuel and air, inert gas sandwiched between fuel and air, cold fuel gas sandwiched between fuel and air. The numerical calculations has been carried out in case of 1, 2, 4 mm of mid layer thickness. The height of total gas stream is 4 cm. The combustion region is broadened in case of inert gas layer of 2, 4 mm thickness and cold fuel layer of 4 mm thickness compared with single shear layer.

• Interaction and multiscale mechanics
• /
• v.2 no.2
• /
• pp.177-187
• /
• 2009

Machine or structural members subjected to fatigue loading will have a crack initiated during early part of their life. Therefore analysis of members with cracks and other discontinuities is very important. Finite element method has enjoyed widespread use in engineering, but it is not convenient for crack problems as the region very close to crack tip is to be discretized with very fine mesh. However, as the body force method (BFM), requires only the boundary of the discontinuity (crack or hole) to be discretized it is easy versatile technique to analyze such problems. In the present work fundamental solution for concentrated load x + iy acting in the semi-infinite plate at an arbitrary point $z_0=x_0+iy_0$ is considered. These fundamental solutions are in complex form ${\phi}(z)$ and ${\psi}(z)$ (England 1971). These potentials are known as Melan potentials (Ramakrishna 1994). A crack in the semi-infinite plate as shown in Fig. 1 is considered. This crack is divided into number of divisions. By applying pair of body forces on a division, the resultant forces on the remaining 'N'divisions are to be found for which ${\phi}_1(z)$ and ${\psi}_1(z)$ are derived. Body force method is applied to calculate stress intensity factor for crack in semi-infinite plate. Also for the case of crack emanating from circular hole in semi-infinite plate radial stress, hoop stress and shear stress are calculated around the hole and crack. Convergent results are obtained by body force method. These results are compared with FEM results.

• Magazine of the Korea Concrete Institute
• /
• v.7 no.4
• /
• pp.149-156
• /
• 1995

The fracture process zone in concrete is a region ahead of a traction-free crack, in which two major mechanisms, microcracking and bridging, play important roles. The toughness due to bridging is dominant compared to toughness induced by microcracking, so that the bridging is dominani: mechanism governing the fracture process of concrete. Fracture mechanics does work for concrete provided that the fracture process zone is being considered, so that the development of model for the fracture process zone is most important to describe fracture phenomena in concrete. In this paper the bridging zone, which is a part of extended rnacrocrack with stresses transmitted by aggregates in concrete, is modelled by a Dugdale-Barenblatt type model with linear tension-softening curve. Two finite element techniques are shown for the analysis of progressive cracking in concrete based on the discrete crack approach: one with crack element, the other without crack element. The advantage of the technique with crack element is that it dees not need to update the mesh topology to follow the progressive cracking. Numerical results by the techniques are demonstrated.

• Journal of the Korea Concrete Institute
• /
• v.16 no.2 s.80
• /
• pp.155-162
• /
• 2004

An analytical method to predict the time dependent flexural behavior of composite girder is presented based on sectional analysis. The time dependent constitutive relation accounting for the early-age concrete properties including maturing of elastic modulus, creep and shrinkage is derived in an incremental format by the first order Taylor series expansion. The sectional analysis calculates the axial and curvature strains based on the force and moment equilibriums. The deflection curve of the girder approximated by the quadratic polynomial function is calculated by applying to the proper boundary conditions in the consecutive segments. Numerical applications are made for the 3-span double composite steel box girder which is a composite bridge girder filled with concrete at the bottom of the steel box in the negative moment region. The calculated results are compared with those by finite element analysis results. Close agreement is observed between the two approaches.