The models of stress corrosion and pressure solution established by Yasuhara et al. were introduced into the 2D FEM code of thermo-hydro-mechanical-migratory coupling analysis for dual-porosity medium developed by the authors. Aiming at a hypothetical model for geological disposal of nuclear waste in an unsaturated rock mass from which there is a nuclide leak, two computation conditions were designed. Then the corresponding two-dimensional numerical simulation for the coupled thermo-hydro-mechanical-migratory processes were carried out, and the states of temperatures, rates and magnitudes of aperture closure, pore and fracture pressures, flow velocities, nuclide concentrations and stresses in the rock mass were investigated. The results show: the aperture closure rates caused by stress corrosion are almost six orders higher than those caused by pressure solution, and the two kinds of closure rates climb up and then decline, furthermore tend towards stability; when the effects of stress corrosion and pressure solution are considered, the negative fracture pressures in near field rise very highly; the fracture aperture and porosity are decreases in the case 1, so the relative permeability coefficients reduce, therefore the nuclide concentrations in pore and fracture in this case are higher than those in case 2.

A stabilized Incompressible Smoothed Particle Hydrodynamics (ISPH) is proposed to simulate free surface flow problems. In the ISPH, pressure is evaluated by solving pressure Poisson equation using a semi-implicit algorithm based on the projection method. Even if the pressure is evaluated implicitly, the unrealistic pressure fluctuations cannot be eliminated. In order to overcome this problem, there are several improvements. One is small compressibility approach, and the other is introduction of two kinds of pressure Poisson equation related to velocity divergence-free and density invariance conditions, respectively. In this paper, a stabilized formulation, which was originally proposed in the framework of Moving Particle Semi-implicit (MPS) method, is applied to ISPH in order to relax the density invariance condition. This formulation leads to a new pressure Poisson equation with a relaxation coefficient, which can be estimated by a preanalysis calculation. The efficiency of the proposed formulation is tested by a couple of numerical examples of dam-breaking problem, and its effects are discussed by using several resolution models with different particle initial distances. Also, the effect of eddy viscosity is briefly discussed in this paper.

Pulsatile flow of blood in constricted narrow arteries under periodic body acceleration is analyzed, modeling blood as non-Newtonian fluid models with yield stress such as (i) Herschel-Bulkley fluid model and (ii) Casson fluid model. The expressions for various flow quantities obtained by Sankar and Ismail (2010) for Herschel-Bulkley fluid model and Nagarani and Sarojamma (2008), in an improved form, for Casson fluid model are used to compute the data for comparing these fluid models. It is found that the plug core radius and wall shear stress are lower for H-B fluid model than those of the Casson fluid model. It is also noted that the plug flow velocity and flow rate are considerably higher for H-B fluid than those of the Casson fluid model. The estimates of the mean velocity and mean flow rate are considerably higher for H-B fluid model than those of the Casson fluid model.

The spreading of oil in an open ocean may cause serious damage to a marine environmental system. Thus, an accurate prediction of oil spill is very important to minimize coastal damage due to unexpected oil spill accident. The movement of oil may be represented with a numerical model that solves an advection-diffusion-reaction equation with a proper numerical scheme. In this study, the spilled oil dispersion model has been established in consideration of tide and tidal currents simultaneously. The velocity components in the advection-diffusion-reaction equation are obtained from the shallow-water equations. The accuracy of the model is verified by applying it to a simple but significant problem. The results produced by the model agree with corresponding analytical solutions and field-observed data. The model is then applied to predict the spreading of an oil spill in a real coastal environment.

The natural convection in an inclined porous square cavity is investigated numerically. The left wall is assumed to have spatial sinusoidal temperature variations about a constant mean value, while the right wall is cooled. The horizontal walls are considered adiabatic. A finite difference method is used to solve numerically the nondimensional governing equations. The effects of the inclination angle of the cavity, the amplitude and wave numbers of the heated sidewall temperature variation on the natural convection in the cavity are studied. The maximum average Nusselt number occurs at different wave number. It also found that the inclination could influence the Nusselt number.

Lattice Boltzmann method is implemented to study hydrodynamically and thermally developing steady laminar flows in a channel. Numerical simulation of two-dimensional convective heat transfer problem is conducted using two-dimensional, nine directional D2Q9 thermal lattice Boltzmann arrangements. The velocity and temperature profiles in the developing region predicted by Lattice Boltzmann method are compared against those obtained by ANSYS-FLUENT. Velocity and temperature profiles as well as the skin friction and the Nusselt numbers agree very well with those predicted by the self-similar solutions of fully developed flows. It is clearly shown here that thermal lattice Boltzmann method is an effective computational fluid dynamics (CFD) tool to study nonisothermal flow problems.

The dispersion coefficient tensor including off-diagonal components was introduced in the flow with secondary currents, which is called skewed shear flow dispersion (SSFD) coefficient tensor, in this paper. To observe the detailed effect of cross-dispersion terms in SSFD model on solute dispersion, mathematical analysis of eigenvalue problem with respect to the equation with SSFD coefficient tensor was performed. The analysis results show the several differences of SSFD model compared to CSFD (conventional shear flow dispersion) model: the oblique direction of principal dispersion with respect to the streamline, the increase of peak concentration, and the change in the eccentricity of elliptical tracer cloud. SSFD coefficient tensor in a streamwise curvilinear coordinate system of curved channel was transformed to those components of fixed Cartesian coordinate system, and 2D numerical model with finite element method was established in the Eulerian-Cartesian coordinate. Through this process, the transformation equation using the depth-averaged velocity field was derived. Several numerical tests were performed to assure the results obtained in the mathematical analysis and to show the applicability of the derived transformation equation on the flow with continuously changing flow direction.

Offshore structures are occasionally located at a relatively deep water region, the outside of breakwater. In this case, these structures may be damaged by the supposition of incident and reflected waves from a vertical breakwater. To prevent the damage, the reflected waves are controlled by installing porous structures at the face of the vertical breakwater. In this study, numerical experiments are carried out to identify the characteristics of wave reflection from the porous structures installing in front of a vertical or slit caisson.

The problem of thermal diffusion and diffusion thermo effects on thermosolutal Marangoni convection flow of an electrically conducting fluid over a permeable surface is investigated. Using appropriate similarity transformations, the governing system of partial differential equation is transformed to a set of nonlinear ordinary differential equations, then solved numerically using the Runge-Kutta-Fehlberg method. The effects of thermal diffusion and diffusion thermo, magnetic field parameter, thermosolutal surface tension ratio, and suction/injection parameter on the flow field, heat transfer characteristic, and concentration are thoroughly examined. Numerical results are obtained for temperature and concentration profiles as well as the local Nusselt and Sherwood numbers are presented graphically and analyzed. It is found that these governing parameters affect the variations of the temperature and concentration and also the local Nusselt and Sherwood numbers.

The flow and heat transfer characteristics of incompressible viscous flow over a nonlinearly stretching sheet with the presence of viscous dissipation is investigated numerically. The similarity transformation reduces the time-independent boundary layer equations for momentum and thermal energy into a set of coupled ordinary differential equations. The obtained equations, including nonlinear equation for the velocity field $f$ and differential equation by variable coefficient for the temperature field $\theta$, are solved numerically by using the fourth order of Runge-Kutta integration scheme accompanied by shooting technique with Newton-Raphson iteration method. The effect of various values of Prandtl number, Eckert number and nonlinear stretching parameter are studied. The results presented graphically show some behaviors such as decrease in dimensionless temperature $\theta$ due to increase in Pr number, and curve relocations are observed when heat dissipation is considered.

The flow of leaked hydrogen gas in tunnel structures is simulated through a free, open source computational fluid dynamics (CFD) code for incompressible thermal convection flow. A one-fifth scale experimental model of a real tunnel is the target model to be simulated. To achieve this, studies on the effects of different boundary conditions, in particular, the wind speed, are carried out on smaller tunnel structures with the same hydrogen inlet boundary conditions. The results suggest a threshold/limiting value of air speed through tunnel. The target model computed with the most suitable boundary conditions shows some agreement with the experimental ones. A method to compute the buoyancy factor used in the code is also presented.

This paper aims at analyzing three different solutions suggested for traffic congestion relief in Port Louis, the busiest city of Mauritius. It evaluates the impact of the three alternatives which are the use of Light Rail Transit (LRT) as an alternative mode of transport, the construction of a Ring road around Port Louis, and the upgrading of the current bus network into a Bus Rapid Transit (BRT) system. The impact of these three solutions has been evaluated by performing Traffic Cellular Automata (TCA) simulations. Our studies reveal that the Ring road will lead to more congestion while introducing the LRT or upgrading the current bus network will reduce congestion significantly.

Radiant floor heating has received increasing attention due to its diverse advantages, especially the energy saving as compared to the conventional dwelling heating system. This paper presents a numerical investigation of airflow and heat transfer in the slot-vented room with the radiant floor heating unit. Combination of fluid convection and thermal radiation has been implemented through the thermal boundary conditions. Spatial distributions of indoor air temperature and velocity, as well as the heat transfer rates along the radiant floor and the outer wall, have been presented and analyzed covering the domains from complete natural convection to forced convection dominated flows. The numerical results demonstrate that the levels of average temperature in the room with lateral slot-ventilation are higher than those without slot-ventilation, but lower than those in the room with ceiling slot-ventilation. Overall, the slot-ventilation room with radiant floor heating unit could offer better indoor air quality through increasing the indoor air temperature and fresh air exchanging rate simultaneously. Concerning the airborne pollutant transports and moisture condensations, the performance of radiant floor heating unit will be further optimized in our future researches.

This paper investigates the determinants of intraindustry trade (IIT), horizontal IIT (HIIT), and Vertical IIT (VIIT) in the automobile industry in Portugal. The trade in this sector between Portugal and the European Union (EU-27) was examined, between 1995 and 2008, using a dynamic panel data. We apply the GMM system to solve the problems of serial correlation and the endogeneity of some explanatory variables. The findings are consistent with the literature. The difference between per capita incomes and factor endowments present a positive sign. These results are according to Heckscher-Ohlin predictions. The economic dimension has a positive impact on trade. A negative effect of the distance on bilateral trade was expected and the results confirm this, underlining the importance of neighbour partnerships for all trade.

When the permeability field of a given porous medium domain is heterogeneous by the existence of randomly distributed fractures such that numerical investigation becomes cumbersome, another level of upscaling may be required. That is such complex permeability field could be relaxed (i.e., smoothed) by constructing an effective permeability field. The effective permeability field is an approximation to the real permeability field that preserves certain quantities and provides an overall acceptable description of the flow field. In this work, the effective permeability for a fractured rock system is obtained for different coarsening scenarios starting from very coarse mesh all the way towards the fine mesh simulation. In all these scenarios, the effective permeability as well as the pressure at each cell is obtained. The total flux at the exit boundary is calculated in all these cases, and very good agreement is obtained.

The purpose of this study is to model the flow movement in an idealized dam-break configuration. One-dimensional and two-dimensional motion of a shallow flow over a rigid inclined bed is considered. The resulting shallow water equations are solved by finite volumes using the Roe and HLL schemes. At first, the one-dimensional model is considered in the development process. With conservative finite volume method, splitting is applied to manage the combination of hyperbolic term and source term of the shallow water equation and then to promote 1D to 2D. The simulations are validated by the comparison with flume experiments. Unsteady dam-break flow movement is found to be reasonably well captured by the model. The proposed concept could be further developed to the numerical calculation of non-Newtonian fluid or multilayers fluid flow.

A new method for solving nonlinear Volterra-Fredholm-Hammerstein (VFH) integral equations is presented. This method is based on reformulation of VFH to the simple form of Fredholm integral equations and hence converts it to optimal control problem. The existence and uniqueness of proposed method are achieved. Numerical results are given at the end of this paper.

The numerical simulation of aeroacoustic phenomena requires high-order accurate numerical schemes with low dispersion and low dissipation errors. A technique has recently been devised in a Computational Fluid Dynamics framework which enables optimal parameters to be chosen so as to better control the grade and balance of dispersion and dissipation in numerical schemes (Appadu and Dauhoo, 2011; Appadu, 2012a; Appadu, 2012b; Appadu, 2012c). This technique has been baptised as the Minimized Integrated Exponential Error for Low Dispersion and Low Dissipation (MIEELDLD) and has successfully been applied to numerical schemes discretising the 1-D, 2-D, and 3-D advection equations. In this paper, we extend the technique of MIEELDLD to the field of computational aeroacoustics and have been able to construct high-order methods with Low Dispersion and Low Dissipation properties which approximate the 1-D linear advection equation. Modifications to the spatial discretization schemes designed by Tam and Webb (1993), Lockard et al. (1995), Zingg et al. (1996), Zhuang and Chen (2002), and Bogey and Bailly (2004) have been obtained, and also a modification to the temporal scheme developed by Tam et al. (1993) has been obtained. These novel methods obtained using MIEELDLD have in general better dispersive properties as compared to the existing optimised methods.

Numerical simulation is performed to investigate the laminar force convection of Al₂O₃/water nanofluid in a flow channel with discrete heat sources. The heat sources are placed on the bottom wall of channel which produce much thermal energy that must be evacuated from the system. The remaining surfaces of channel are kept adiabatic to exchange energy between nanofluid and heat sources. In the present study the effects of Reynolds number ($\text{Re}=\mathrm{50,\; 100,\; 200,\; 400}$, and 1000), particle volume fraction ($\varphi =0$ (distilled water), 1 and 4%) on the average heat transfer coefficient (h), pressure drop ($\Delta P$), and wall temperature ($T_w$) are evaluated. The use of nanofluid can produce an asymmetric velocity along the height of the channel. The results show a maximum value 38% increase in average heat transfer coefficient and 68% increase in pressure drop for all the considered cases when compared to basefluid (i.e., water). It is also observed that the wall temperature decreases remarkably as Re and ϕ increase. Finally, thermal-hydraulic performance (η) is evaluated and it is seen that best performance can be obtained for $\text{Re}=1000$ and $\varphi =4$%.

The newly developed algorithm called the grid-by-grid inversion method is a very convenient method for converting an existing computer code for Newtonian flow simulations to that for viscoelastic flow simulations. In this method, the hyperbolic constitutive equation is split such that the term for the convective transport of stress tensor is treated as a source which is updated iteratively. This allows the stress tensors at each grid point to be expressed in terms of velocity gradient tensor at the same location, and the set of stress tensor components is found after inverting a small matrix at each grid point. To corroborate the robustness and accuracy of the grid-by-grid inversion method, we apply it to the 4 : 1 axisymmetric contraction problem. This algorithm is found to be robust and yields accurate results as compared with other finite volume methods. Any commercial CFD packages for Newtonian flow simulations can be easily converted to those for viscoelastic fluids exploiting the grid-by-grid inversion method.

Understanding platoon dispersion is critical for the coordination of traffic signal control in an urban traffic network. Assuming that platoon speed follows a truncated normal distribution, ranging from minimum speed to maximum speed, this paper develops a piecewise density function that describes platoon dispersion characteristics as the platoon moves from an upstream to a downstream intersection. Based on this density function, the expected number of cars in the platoon that pass the downstream intersection, and the expected number of cars in the platoon that do not pass the downstream point are calculated. To facilitate coordination in a traffic signal control system, dispersion models for the front and the rear of the platoon are also derived. Finally, a numeric computation for the coordination of successive signals is presented to illustrate the validity of the proposed model.

The effects of a porous fence with a nonuniform porosity on flow fields are investigated numerically. First, an experiment with a non-uniform porous fence located in a wind tunnel is performed to obtain a reference data set. Then, a numerical model that utilizes the finite volume scheme with a weakly compressible-flow method to solve the continuity and momentum equations is developed. The numerical simulation is compared to experimental measurements for validation purposes. As a result, the numerical predictions show good agreements with the experimental data. Finally, the numerical investigations of the flow fields around porous fences with various combinations of upper and lower fence porosity are also presented. When the upper porosity is greater than the lower porosity, the Protection Index PI_0.1, PI_0.3 and PI_0.5, representing the adverse sheltering effect, decreases compared to that of the uniform porous fence. When the upper porosity is less than the lower porosity, the PI_0.5 increases and the variations of the PI_0.1 and PI_0.3, depend on the upper porosity, compared to that of the uniform porous fence. The results show that the porous fence with the upper fence porosity ${\epsilon }_U=0\%$ and the lower fence porosity ${\epsilon }_L=30\%$ gives the best sheltering effect among the porous fences in this study.

The mechanism of inclusion removal from liquid steel by gas bubbling and bubble attachment in the tundish is complex due to the great number of variables involved, and it is even more difficult to study because of the turbulent flow conditions. The main objective of this work is to analyze and improve the understanding of the alumina inclusion removal rate by bubble attachment and by gas bubbling fluid dynamics effects. The results show that the inclusion collection probability mainly depends on the attachment mechanism by collision. This parameter was determined by calculating the induction time, which is shorter when the rupture time and the formation time of a stable three phases contact (particle/liquid/gas) are ignored than when it is fully considered, affecting the attachment probability. In addition, to achieve acceptable inclusion removal, a smaller bubble diameter is required, such as 1 mm. This consideration is almost impossible to achieve during tundish operation; a more realistic bubble diameter around 10 mm is employed, resulting in a very inefficient inclusion removal process by bubble attachment. Nevertheless, in a real casting tundish the inclusion removal rate employing argon bubbling is efficient; is mainly due to the fluid flow pattern changes rather than bubble attachment. Consequently, it is imperative to consider the summation of both removal mechanisms to compute a better approximation of this important operation.

This work numerically studies the flow pattern of a magnetic fluid filled within an annulus whose inner cylinder is moving at a constant rotational speed, while the outer cylinder is stationary but under the influence of a nonuniform external magnetic field. The magnetic field consists of four basic configurations, that is, completely circular, semicircular, quarter circular, and alternately quarter circular. The strength of the external magnetic field is characterized using a reference Hartmann number. As the reference Hartmann number increases, the fluid elements need to overcome greater resistance to enter the region with magnetic field. Hence, there always exists an apparent recirculation cell within the region without externally applied magnetic field. The strength and size of the recirculation cell depend on the reference Hartmann number, the number and size of the discrete regions without external magnetic field. Only the shear stress on the moving cylinder always increases in magnitude with the reference Hartmann number and the span of the single external magnetic field region. Splitting and separating the external magnetic field may increase the magnitude of the shear stress on the moving inner cylinder but decrease that on the stationary outer cylinder. If the magnitude of the shear stress on the outer cylinder reduces beyond zero, a shear stress in the opposite sense will increase in magnitude with Hartmann number.

Semianalytical solutions are developed for turbulent hydrogen-air plume. We derived analytical expressions for plume centerline variables (radius, velocity, and density deficit) in terms of a single universal function, called plume function. By combining the obtained analytical expressions of centerline variables with empirical Gaussian expressions of the mean variables, we obtain semianalytical expressions for mean quantities of hydrogen-air plume (velocity, density deficit, and mass fraction).