• Journal of applied mathematics & informatics
• /
• v.40 no.3_4
• /
• pp.603-618
• /
• 2022

In this paper, a numerical solution of the generalized diffusion equation with a delay has been obtained by a numerical technique based on the Galerkin finite element method by applying the cubic B-spline basis functions. The time discretization process is carried out using the forward Euler method. The numerical scheme is required to preserve the delay-independent asymptotic stability with an additional restriction on time and spatial step sizes. Both the theoretical and computational rates of convergence of the numerical method have been examined and found to be in agreement. As it can be observed from the numerical results given in tables and graphs, the proposed method approximates the exact solution very well. The accuracy of the numerical scheme is confirmed by computing L₂ and L_∞ error norms.

• Asian Journal of Atmospheric Environment
• /
• v.2 no.1
• /
• pp.1-13
• /
• 2008

A numerical simulation method has been developed to predict atmospheric flow and stack gas diffusion using a calculation domain of several km around a stack under complex terrain conditions containing buildings. The turbulence closure technique using a modified k-$\varepsilon$-type model under a non hydrostatic assumption was used for the flow calculation, and some of the calculation grids near the ground were treated as buildings using a terrain-following coordinate system. Stack gas diffusion was predicted using the Lagrangian particle model, that is, the stack gas was represented by the trajectories of released particles. The numerical model was applied separately to the flow and stack gas diffusion around a cubical building and to a two-dimensional ridge in this study, before being applied to an actual terrain containing buildings in our next study. The calculated flow and stack gas diffusion results were compared with those obtained by wind tunnel experiments, and the features of flow and stack gas diffusion, such as the increase in turbulent kinetic energy and the plume spreads of the stack gas behind the building and ridge, were reproduced by both calculations and wind tunnel experiments. Furthermore, the calculated profiles of the mean velocity, turbulent kinetic energy and concentration of the stack gas around the cubical building and the ridge showed good agreement with those of wind tunnel experiments.

• Environmental Engineering Research
• /
• v.23 no.4
• /
• pp.442-448
• /
• 2018

This paper presents the numerical simulation of advection-diffusion mechanism of BOD concentration which was used as an indicator of waste only in one flow-direction of waste stabilization ponds (1-dimension (1-D)). This model was represented in partial differential equation order 2. The purpose of this paper was to determine the simulation of the model 1-D of wastewater transport phenomena based advection-diffusion mechanism and did validate the model. Numerical methods which was used for the solution of this model is finite difference method with Forward Time Central Space scheme. The simulation results which was obtained would be compared with field observation data as a validation model. Collection of field data was carried out in the Wastewater Treatment Plant Sewon, Bantul, D.I. Yogyakarta. The results of numerical simulations were indicate that the advection-diffusion mechanism takes place continuously over time. Then validation of the model was state that there was a difference between the calculation results with the field data, with a correlation value of 0.998.

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.63-72
• /
• 1998

The numerical simulations of flowfield and pollutant dispersion over two-dimensional hills of various shapes are described. The Reynolds-averaged Wavier-Stokes equations and concentration diffusion equation based on the gradient diffusion theory have been applied to the atmospheric shear flow over the bell-shaped hills which are basic components of the complex terrain. The flow characteristics such as velocity profiles of the geophysical boundary layer, speed-up phenomena, mean pollutant concentration profiles are compared with experimental data to validate the present numerical procedure and it has been found that the present numerical results agree well with experiments and other numerical data. It has been also found that the distributions of ground level concentration are strongly influenced by the source location and height.

• Journal of Ocean Engineering and Technology
• /
• v.24 no.4
• /
• pp.78-85
• /
• 2010

This study investigates reducing characteristics for the spreading of dredged soil and the diffusion of contaminant by silt protector curtain through three dimensional numerical experiment. The numerical medel is modified by combining the sediment transport characteristics for cohesive sediment into the previously developed model. Several numerical experiments have been given in order to investigate the reducing effect of silt protector using two dimensional numerical channel model under various parameters such as upstream flow velocity, depth of silt curtain and the position of dumped materials. Through the evaluation of several simulation results, we knew that the careful design has to be given in the determination of depth and position of silt protector.

• Journal of Environmental Science International
• /
• v.11 no.12
• /
• pp.1321-1326
• /
• 2002

A numerical model is investigated to predict a behavior of the gaseous volatile organic compounds and a subsurface contamination caused by them in the unsaturated zone. Two dimensional advective-dispersion equation caused by a density difference and two dimensional diffusion equation are computed by a finite difference method in the numerical model. A laboratory experiment is also carried out to compare the results of the numerical model. The dimensions of the experimental plume are 1.2m in length, 0.5m in height, and 0.05m in thickness. In comparing the result of 2 methods used in the numerical model with the one of the experiment respectively, the one of the advective-dispersion equation shows better than the one the diffusion equation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.4
• /
• pp.264-271
• /
• 2009

Characteristics of a counter flowing diffusion flame, which is formulated by an oppositely-injected methane-jet flow in a narrow channel of a uniform air flow. The location of the flame fronts and the flame lengths were compared by changing the flow rates of fuel. To distinguish the effects of the narrow channel on the diffusion flame, a numerical simulation for an ideal two-dimensional flame was conducted. Overall trends of the flame behavior were similar in both numerical and experimental results. With the increase of the ratio of jet velocity to air velocity flame front moved farther upstream. It is thought that the flow re-direction in the channel suppresses fuel momentum more significantly due to the higher temperature and increased viscosity of burned gas. Actual flames in a narrow channel suffer heat loss to the ambient and it has finite length of diffusion flame in contrast to the numerical results of infinite flame length. Thus a convective heat loss was additionally employed in numerical simulation and closer results were obtained. These results can be used as basic data in development of a small combustor of a nonpremixed flame.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.37A no.11
• /
• pp.961-971
• /
• 2012

This article is concerned with effective numerical techniques for partial differential equation (PDE)-based image restoration. Numerical realizations of most PDE-based denoising models show a common drawback: loss of fine structures. In order to overcome the drawback, the article introduces a new time-stepping procedure, called the method of nonflat time evolution (MONTE), in which the timestep size is determined based on local image characteristics such as the curvature or the diffusion magnitude. The MONTE provides PDE-based restoration models with an effective mechanism for the equalization of the net diffusion over a wide range of image frequency components. It can be easily applied to diverse evolutionary PDE-based restoration models and their spatial and temporal discretizations. It has been numerically verified that the MONTE results in a significant reduction in numerical dissipation and preserves fine structures such as edges and textures satisfactorily, while it removes the noise with an improved efficiency. Various numerical results are shown to confirm the claim.

• Computers and Concrete
• /
• v.21 no.4
• /
• pp.471-484
• /
• 2018

Refined analysis depicting mass transportation and physicochemical reaction and reasonable computing load with acceptable DOFs are the two major challenges of numerical simulation for concrete durability. Mesoscopic numerical simulation for chloride diffusion considering binder, aggregate and interfacial transition zone is unable to be expended to the full structure due to huge number of DOFs. In this paper, a multiscale approach of combining both mesoscopic model including full-graded aggregate and equivalent macroscopic model was introduced. An equivalent conversion of chloride content at the Interfacial Transition Layer (ITL) connecting both models was considered. Feasibility and relative error were discussed by analytical deduction and numerical simulation. Case study clearly showed that larger analysis model in multiscale model expanded the diffusion space of chloride ion and decreased chloride content in front of rebar. Difference for single-scale simulation and multiscale approach was observed. Finally, this paper addressed some worth-noting conclusions about the chloride distribution and rebar corrosion regarding the configuration of rebar placement, rebar diameter, concrete cover and exposure period.

• Nuclear Engineering and Technology
• /
• v.16 no.1
• /
• pp.21-28
• /
• 1984

The reduction of the truncation error including numerical diffusion, has been one of the most important tasks in the development of numerical methods. The stream line method is used to cancel cross numerical diffusion and some of the non-diffusion type truncation error. The two-step stream line method which is the combination of the stream line method and finite difference methods is developed in this work for the solution of the govern ing equations of incompressible buoyant turbulent flow. This method is compared with the finite difference method. The predictions of both classes of numerical methods are compared with experimental findings. Truncation error analysis also has been performed in order to the compare truncation error of the stream line method with that of finite difference methods.