Abstract: In this paper, numerical solution for the generalized Rosenau-Burgers equation is considered and Crank-Nicolson finite difference scheme is proposed. Existence of the solutions for the difference scheme has been shown. Stability, convergence and priori error estimate of the scheme are proved. Numerical results demonstrate that the scheme is efficient and reliable.
Abstract: A conjugate heat transfer for steady two-dimensional
mixed convection with magnetic hydrodynamic (MHD) flow of an
incompressible quiescent fluid over an unsteady thermal forming
stretching sheet has been studied. A parameter, M, which is used to
represent the dominance of the magnetic effect has been presented in
governing equations. The similar transformation and an implicit
finite-difference method have been used to analyze the present
problem. The numerical solutions of the flow velocity distributions,
temperature profiles, the wall unknown values of f''(0) and '(θ (0) for
calculating the heat transfer of the similar boundary-layer flow are
carried out as functions of the unsteadiness parameter (S), the Prandtl
number (Pr), the space-dependent parameter (A) and
temperature-dependent parameter (B) for heat source/sink and the
magnetic parameter (M). The effects of these parameters have also
discussed. At the results, it will produce greater heat transfer effect
with a larger Pr and M, S, A, B will reduce heat transfer effects. At
last, conjugate heat transfer for the free convection with a larger G has
a good heat transfer effect better than a smaller G=0.
Abstract: This paper proposes a simple model of economic geography within the Dixit-Stiglitz-Iceberg framework that may be used to analyze migration patterns among three cities. The cost–benefit tradeoffs affecting incentives for three types of migration, including echelon migration, are discussed. This paper develops a tractable, heterogeneous-agent, general equilibrium model, where agents share constant human capital, and explores the relationship between the benefits of echelon migration and gross human capital. Using Chinese numerical solutions, we study the manifestation of echelon migration and how it responds to changes in transportation cost and elasticity of substitution. Numerical results demonstrate that (i) there are positive relationships between a migration-s benefit-and-wage ratio, (ii) there are positive relationships between gross human capital ratios and wage ratios as to origin and destination, and (iii) we identify 13 varieties of human capital convergence among cities. In particular, this model predicts population shock resulting from the processes of migration choice and echelon migration.
Abstract: The paper discusses a 3D numerical solution of the inverse boundary problem for a continuous casting process of alloy. The main goal of the analysis presented within the paper was to estimate heat fluxes along the external surface of the ingot. The verified information on these fluxes was crucial for a good design of a mould, effective cooling system and generally the whole caster. In the study an enthalpy-porosity technique implemented in Fluent package was used for modeling the solidification process. In this method, the phase change interface was determined on the basis of the liquid fraction approach. In inverse procedure the sensitivity analysis was applied for retrieving boundary conditions. A comparison of the measured and retrieved values showed a high accuracy of the computations. Additionally, the influence of the accuracy of measurements on the estimated heat fluxes was also investigated.
Abstract: The present paper develops and validates a numerical procedure for the calculation of turbulent combustive flow in converging and diverging ducts and throuh simulation of the heat transfer processes, the amount of production and spread of Nox pollutant has been measured. A marching integration solution procedure employing the TDMA is used to solve the discretized equations. The turbulence model is the Prandtl Mixing Length method. Modeling the combustion process is done by the use of Arrhenius and Eddy Dissipation method. Thermal mechanism has been utilized for modeling the process of forming the nitrogen oxides. Finite difference method and Genmix numerical code are used for numerical solution of equations. Our results indicate the important influence of the limiting diverging angle of diffuser on the coefficient of recovering of pressure. Moreover, due to the intense dependence of Nox pollutant to the maximum temperature in the domain with this feature, the Nox pollutant amount is also in maximum level.
Abstract: Recently, it is found that telegraph equation is more suitable than ordinary diffusion equation in modelling reaction diffusion for such branches of sciences. In this paper, a numerical solution for the one-dimensional hyperbolic telegraph equation by using the collocation method using the septic splines is proposed. The scheme works in a similar fashion as finite difference methods. Test problems are used to validate our scheme by calculate L2-norm and L∞-norm. The accuracy of the presented method is demonstrated by two test problems. The numerical results are found to be in good agreement with the exact solutions.
Abstract: In this paper the Laplace Decomposition method is developed to solve linear and nonlinear fractional integro- differential equations of Volterra type.The fractional derivative is described in the Caputo sense.The Laplace decomposition method is found to be fast and accurate.Illustrative examples are included to demonstrate the validity and applicability of presented technique and comparasion is made with exacting results.
Abstract: This is the second part of the paper. It, aside from the
core subroutine test reported previously, focuses on the simulation of
turbulence governed by the full STF Navier-Stokes equations on a
large scale. Law of the wall is found plausible in this study as a model
of the boundary layer dynamics. Model validations proceed to
include velocity profiles of a stationary turbulent Couette flow, pure
sloshing flow simulations, and the identification of water-surface
inclination due to fluid accelerations. Errors resulting from the
irrotational and hydrostatic assumptions are explored when studying
a wind-driven water circulation with no shakings. Illustrative
examples show that this numerical strategy works for the simulation
of sloshing-shear mixed flow in a 3-D rigid rectangular base tank.
Abstract: Nowadays, there is little information, concerning the
heat shield systems, and this information is not completely reliable to
use in so many cases. for example, the precise calculation cannot be
done for various materials. In addition, the real scale test has two
disadvantages: high cost and low flexibility, and for each case we
must perform a new test. Hence, using numerical modeling program
that calculates the surface recession rate and interior temperature
distribution is necessary. Also, numerical solution of governing
equation for non-charring material ablation is presented in order to
anticipate the recession rate and the heat response of non-charring
heat shields. the governing equation is nonlinear and the Newton-
Rafson method along with TDMA algorithm is used to solve this
nonlinear equation system. Using Newton- Rafson method for
solving the governing equation is one of the advantages of the
solving method because this method is simple and it can be easily
generalized to more difficult problems. The obtained results
compared with reliable sources in order to examine the accuracy of
compiling code.
Abstract: In the present study, a procedure was developed to
determine the optimum reaction rate constants in generalized
Arrhenius form and optimized through the Nelder-Mead method. For
this purpose, a comprehensive mathematical model of a fixed bed
reactor for dehydrogenation of heavy paraffins over Pt–Sn/Al2O3
catalyst was developed. Utilizing appropriate kinetic rate expressions
for the main dehydrogenation reaction as well as side reactions and
catalyst deactivation, a detailed model for the radial flow reactor was
obtained. The reactor model composed of a set of partial differential
equations (PDE), ordinary differential equations (ODE) as well as
algebraic equations all of which were solved numerically to
determine variations in components- concentrations in term of mole
percents as a function of time and reactor radius. It was demonstrated
that most significant variations observed at the entrance of the bed
and the initial olefin production obtained was rather high. The
aforementioned method utilized a direct-search optimization
algorithm along with the numerical solution of the governing
differential equations. The usefulness and validity of the method was
demonstrated by comparing the predicted values of the kinetic
constants using the proposed method with a series of experimental
values reported in the literature for different systems.
Abstract: This work focuses on analysis of classical heat transfer equation regularized with Maxwell-Cattaneo transfer law. Computer simulations are performed in MATLAB environment. Numerical experiments are first developed on classical Fourier equation, then Maxwell-Cattaneo law is considered. Corresponding equation is regularized with a balancing diffusion term to stabilize discretizing scheme with adjusted time and space numerical steps. Several cases including a convective term in model equations are discussed, and results are given. It is shown that limiting conditions on regularizing parameters have to be satisfied in convective case for Maxwell-Cattaneo regularization to give physically acceptable solutions. In all valid cases, uniform convergence to solution of initial heat equation with Fourier law is observed, even in nonlinear case.
Abstract: A mathematical model for the hydrodynamic
lubrication of parabolic slider bearings with couple stress lubricants
is presented. A numerical solution for the mathematical model using
finite element scheme is obtained using three nodes isoparametric
quadratic elements. Stiffness integrals obtained from the weak form
of the governing equations were solved using Gauss Quadrature to
obtain a finite number of stiffness matrices. The global system of
equations was obtained for the bearing and solved using Gauss Seidel
iterative scheme. The converged pressure solution was used to obtain
the load capacity of the bearing. Parametric studies were carried out
and it was shown that the effect of couple stresses and profile
parameter are to increase the load carrying capacity of the parabolic
slider bearing. Numerical experiments reveal that the magnitude of
the profile parameter at which maximum load is obtained increases
with decrease in couple stress parameter. The results are presented in
graphical form.
Abstract: In this paper, the local grid refinement is focused by
using a nested grid technique. The Cartesian grid numerical method is
developed for simulating unsteady, viscous, incompressible flows
with complex immersed boundaries. A finite volume method is used in
conjunction with a two-step fractional-step procedure. The key aspects
that need to be considered in developing such a nested grid solver are
imposition of interface conditions on the inter-block and accurate
discretization of the governing equation in cells that are with the
inter-block as a control surface. A new interpolation procedure is
presented which allows systematic development of a spatial
discretization scheme that preserves the spatial accuracy of the
underlying solver. The present nested grid method has been tested by
two numerical examples to examine its performance in the two
dimensional problems. The numerical examples include flow past a
circular cylinder symmetrically installed in a Channel and flow past
two circular cylinders with different diameters. From the numerical
experiments, the ability of the solver to simulate flows with
complicated immersed boundaries is demonstrated and the nested grid
approach can efficiently speed up the numerical solutions.
Abstract: Today, numerical simulation is a powerful tool to
solve various hydraulic engineering problems. The aim of this
research is numerical solutions of shallow water equations using
finite volume method for Simulations of dam break over wet and dry
bed. In order to solve Riemann problem, Roe-s approximate solver is
used. To evaluate numerical model, simulation was done in 1D and
2D states. In 1D state, two dam break test over dry bed (with and
without friction) were studied. The results showed that Structural
failure around the dam and damage to the downstream constructions
in bed without friction is more than friction bed. In 2D state, two
tests for wet and dry beds were done. Generally in wet bed case,
waves are propagated to canal sides but in dry bed it is not
significant. Therefore, damage to the storage facilities and
agricultural lands in wet bed case is more than in dry bed.
Abstract: An unstructured finite volume numerical model is
presented here for simulating shallow-water flows with wetting and
drying fronts. The model is based on the Green-s theorem in
combination with Chorin-s projection method. A 2nd-order upwind
scheme coupled with a Least Square technique is used to handle
convection terms. An Wetting and drying treatment is used in the
present model to ensures the total mass conservation. To test it-s
capacity and reliability, the present model is used to solve the
Parabolic Bowl problem. We compare our numerical solutions with
the corresponding analytical and existing standard numerical results.
Excellent agreements are found in all the cases.
Abstract: The optimal control problem for the viscoelastic melt
spinning process has not been reported yet in the literature. In this
study, an optimal control problem for a mathematical model of a
viscoelastic melt spinning process is considered. Maxwell-Oldroyd
model is used to describe the rheology of the polymeric material, the
fiber is made of. The extrusion velocity of the polymer at the spinneret
as well as the velocity and the temperature of the quench air and the
fiber length serve as control variables. A constrained optimization
problem is derived and the first–order optimality system is set up
to obtain the adjoint equations. Numerical solutions are carried out
using a steepest descent algorithm. A computer program in MATLAB
is developed for simulations.
Abstract: Based on the standard finite element method, a new
finite element method which is known as nonlocal finite element
method (NL-FEM) is numerically implemented in this article to
study the nonlocal effects for solving 1D nonlocal elastic problem.
An Eringen-type nonlocal elastic model is considered. In this model,
the constitutive stress-strain law is expressed interms of integral
equation which governs the nonlocal material behavior. The new
NL-FEM is adopted in such a way that the postulated nonlocal elastic
behavior of material is captured by a finite element endowed with a
set of (cross-stiffness) element itself by the other elements in mesh.
An example with their analytical solutions and the relevant numerical
findings for various load and boundary conditions are presented and
discussed in details. It is observed from the numerical solutions that
the torsional deformation angle decreases with increasing nonlocal
nanoscale parameter. It is also noted that the analytical solution fails
to capture the nonlocal effect in some cases where numerical
solutions handle those situation effectively which prove the
reliability and effectiveness of numerical techniques.
Abstract: In this study, the numerical solution of unsteady flow
between two concentric rotating spheres with suction and blowing at
their boundaries is presented. The spheres are rotating about a
common axis of rotation while their angular velocities are constant.
The Navier-Stokes equations are solved by employing the finite
difference method and implicit scheme. The resulting flow patterns
are presented for various values of the flow parameters including
rotational Reynolds number Re , and a blowing/suction Reynolds
number Rew . Viscous torques at the inner and the outer spheres are
calculated, too. It is seen that increasing the amount of suction and
blowing decrease the size of eddies generated in the annulus.
Abstract: This paper presents results obtained from the
numerical solution for the flow past an oscillating circular cylinder at
Reynolds number of 200. The frequency of oscillation was fixed to
the vortex shedding frequency from a fixed cylinder, f0, while the
amplitudes of oscillations were varied from to 1.1a, where a
represents the radius of the cylinder. The response of the flow
through the fluid forces acting on the surface of the cylinder are
investigated. The lock-on phenomenon is captured at low oscillation
amplitudes.
Abstract: The present work deals with analyses of the effects
of bearing curvature and non-Newtonian characteristics on the load capacity of an exponential rectangular squeeze film bearing using
Bingham fluids as lubricants. Bingham fluids are characterized by an
yield value and hence the formation of a “rigid" core in the region
between the plates is justified. The flow is confined to the region
between the core and the plates. The shape of the core has been
identified through numerical means. Further, numerical solutions for
the pressure distribution and load carrying capacity of the bearing
for various values of Bingham number and curvature parameter have
been obtained. The effects of bearing curvature and non-Newtonian
characteristics of the lubricant on the bearing performances have been
discussed.