Abstract: Formulation for drying and pyrolysis process in packed beds at slow heating rates is presented. Drying of biomass particles bed is described by mass diffusion equation and local moisture-vapour-equilibrium relations. In gasifiers, volatilization rate during pyrolysis of biomass is modeled by using apparent kinetic rate expression, while product compositions at slow heating rates is modeled using empirical fitted mass ratios (i.e., CO/CO2, ME/CO2, H2O/CO2) in terms of pyrolysis temperature. The drying module is validated fairly with available chemical kinetics scheme and found that the testing zone in gasifier bed constituted of relatively smaller particles having high airflow with high isothermal temperature expedite the drying process. Further, volatile releases more quickly within the shorter zone height at high temperatures (isothermal). Both, moisture loss and volatile release profiles are found to be sensitive to temperature, although the influence of initial moisture content on volatile release profile is not so sensitive.
Abstract: This paper sets to demonstrate a modeling of electrokinetic mixing employing electroosmotic stationary and time-dependent microchannel using alternate zeta patches on the lower surface of the micromixer in a lab on chip microfluidic device. Electroosmotic flow is amplified using different 2D and 3D model designs with alternate and geometric zeta potential values such as 25, 50, and 100 mV, respectively, to achieve high concentration mixing in the electrokinetically-driven microfluidic system. The enhancement of electrokinetic mixing is studied using Finite Element Modeling, and simulation workflow is accomplished with defined integral steps. It can be observed that the presence of alternate zeta patches can help inducing microvortex flows inside the channel, which in turn can improve mixing efficiency. Fluid flow and concentration fields are simulated by solving Navier-Stokes equation (implying Helmholtz-Smoluchowski slip velocity boundary condition) and Convection-Diffusion equation. The effect of the magnitude of zeta potential, the number of alternate zeta patches, etc. are analysed thoroughly. 2D simulation reveals that there is a cumulative increase in concentration mixing, whereas 3D simulation differs slightly with low zeta potential as that of the 2D model within the T-shaped micromixer for concentration 1 mol/m3 and 0 mol/m3, respectively. Moreover, 2D model results were compared with those of 3D to indicate the importance of the 3D model in a microfluidic design process.
Abstract: In this paper, we present a simple effective numerical geometric method to estimate the divergence of a vector field over a curved surface. The conservation law is an important principle in physics and mathematics. However, many well-known numerical methods for solving diffusion equations do not obey conservation laws. Our presented method in this paper combines the divergence theorem with a generalized finite difference method and obeys the conservation law on discrete closed surfaces. We use the similar method to solve the Cahn-Hilliard equations on evolving spherical surfaces and observe stability results in our numerical simulations.
Abstract: The diffusion-reaction equations are important Partial Differential Equations in mathematical biology, material science, physics, and so on. However, finding efficient numerical methods for diffusion-reaction systems on curved surfaces is still an important and difficult problem. The purpose of this paper is to present a convergent geometric method for solving the reaction-diffusion equations on closed surfaces by an O(r)-LTL configuration method. The O(r)-LTL configuration method combining the local tangential lifting technique and configuration equations is an effective method to estimate differential quantities on curved surfaces. Since estimating the Laplace-Beltrami operator is an important task for solving the reaction-diffusion equations on surfaces, we use the local tangential lifting method and a generalized finite difference method to approximate the Laplace-Beltrami operators and we solve this reaction-diffusion system on closed surfaces. Our method is not only conceptually simple, but also easy to implement.
Abstract: Stochastic modeling concerns the use of probability
to model real-world situations in which uncertainty is present.
Therefore, the purpose of stochastic modeling is to estimate the
probability of outcomes within a forecast, i.e. to be able to predict
what conditions or decisions might happen under different situations.
In the present study, we present a model of a stochastic diffusion
process based on the bi-Weibull distribution function (its trend
is proportional to the bi-Weibull probability density function). In
general, the Weibull distribution has the ability to assume the
characteristics of many different types of distributions. This has
made it very popular among engineers and quality practitioners, who
have considered it the most commonly used distribution for studying
problems such as modeling reliability data, accelerated life testing,
and maintainability modeling and analysis. In this work, we start
by obtaining the probabilistic characteristics of this model, as the
explicit expression of the process, its trends, and its distribution by
transforming the diffusion process in a Wiener process as shown in
the Ricciaardi theorem. Then, we develop the statistical inference of
this model using the maximum likelihood methodology. Finally, we
analyse with simulated data the computational problems associated
with the parameters, an issue of great importance in its application to
real data with the use of the convergence analysis methods. Overall,
the use of a stochastic model reflects only a pragmatic decision on
the part of the modeler. According to the data that is available and
the universe of models known to the modeler, this model represents
the best currently available description of the phenomenon under
consideration.
Abstract: The modelling of physical phenomena, such as the
earth’s free oscillations, the vibration of strings, the interaction of
atomic particles, or the steady state flow in a bar give rise to Sturm-
Liouville (SL) eigenvalue problems. The boundary applications of
some systems like the convection-diffusion equation, electromagnetic
and heat transfer problems requires the combination of Dirichlet and
Neumann boundary conditions. Hence, the incorporation of Robin
boundary condition in the analyses of Sturm-Liouville problem. This
paper deals with the computation of the eigenvalues and
eigenfunction of generalized Sturm-Liouville problems with Robin
boundary condition using the finite element method. Numerical
solution of classical Sturm–Liouville problem is presented. The
results show an agreement with the exact solution. High results
precision is achieved with higher number of elements.
Abstract: In this paper, a backward semi-Lagrangian scheme
combined with the second-order backward difference formula
is designed to calculate the numerical solutions of nonlinear
advection-diffusion equations. The primary aims of this paper are
to remove any iteration process and to get an efficient algorithm
with the convergence order of accuracy 2 in time. In order to achieve
these objects, we use the second-order central finite difference and the
B-spline approximations of degree 2 and 3 in order to approximate
the diffusion term and the spatial discretization, respectively. For the
temporal discretization, the second order backward difference formula
is applied. To calculate the numerical solution of the starting point
of the characteristic curves, we use the error correction methodology
developed by the authors recently. The proposed algorithm turns out
to be completely iteration free, which resolves the main weakness
of the conventional backward semi-Lagrangian method. Also, the
adaptability of the proposed method is indicated by numerical
simulations for Burgers’ equations. Throughout these numerical
simulations, it is shown that the numerical results is in good
agreement with the analytic solution and the present scheme offer
better accuracy in comparison with other existing numerical schemes.
Abstract: This paper presents a complete dynamic modeling
of a membrane distillation process. The model contains two
consistent dynamic models. A 2D advection-diffusion equation
for modeling the whole process and a modified heat equation
for modeling the membrane itself. The complete model describes
the temperature diffusion phenomenon across the feed, membrane,
permeate containers and boundary layers of the membrane. It gives
an online and complete temperature profile for each point in the
domain. It explains heat conduction and convection mechanisms that
take place inside the process in terms of mathematical parameters, and
justify process behavior during transient and steady state phases. The
process is monitored for any sudden change in the performance at any
instance of time. In addition, it assists maintaining production rates
as desired, and gives recommendations during membrane fabrication
stages. System performance and parameters can be optimized
and controlled using this complete dynamic model. Evolution of
membrane boundary temperature with time, vapor mass transfer along
the process, and temperature difference between membrane boundary
layers are depicted and included. Simulations were performed over
the complete model with real membrane specifications. The plots
show consistency between 2D advection-diffusion model and the
expected behavior of the systems as well as literature. Evolution
of heat inside the membrane starting from transient response till
reaching steady state response for fixed and varying times is
illustrated.
Abstract: Carrier scatterings in the inversion channel of MOSFET dominates the carrier mobility and hence drain current. This paper presents an analytical model of the subthreshold drain current incorporating the effective electron mobility model of the pocket implanted nano scale n-MOSFET. The model is developed by assuming two linear pocket profiles at the source and drain edges at the surface and by using the conventional drift-diffusion equation. Effective electron mobility model includes three scattering mechanisms, such as, Coulomb, phonon and surface roughness scatterings as well as ballistic phenomena in the pocket implanted n-MOSFET. The model is simulated for various pocket profile and device parameters as well as for various bias conditions. Simulation results show that the subthreshold drain current data matches the experimental data already published in the literature.
Abstract: Numerical investigation on the generality of nanoparticle velocity equation had been done on the previous published work. The three dimensional governing equations (continuity, momentum and energy) were solved using finite volume method (FVM). Parametric study of thermal performance between pure water-cooled and nanofluid-cooled are evaluated for volume fraction in the range of 1% to 4%, and nanofluid type of gamma-Al2O3 at Reynolds number range of 67.41 to 286.77. The nanofluid is modeled using single and two phase approach. Three different existing Brownian motion velocities are applied in comparing the generality of the equation for a wide parametric condition. Deviation in between the Brownian motion velocity is identified to be due to the different means of mean free path and constant value used in diffusion equation.
Abstract: For a quick and accurate calculation of spatial neutron
distribution in nuclear power reactors 3D nodal codes are usually
used aiming at solving the neutron diffusion equation for a given
reactor core geometry and material composition. These codes use a
second order polynomial to represent the transverse leakage term. In
this work, a nodal method based on the well known nodal expansion
method (NEM), developed at COPPE, making use of this polynomial
expansion was modified to treat the transverse leakage term for the
external surfaces of peripheral reflector nodes.
The proposed method was implemented into a computational
system which, besides solving the diffusion equation, also solves the
burnup equations governing the gradual changes in material
compositions of the core due to fuel depletion. Results confirm the
effectiveness of this modified treatment of peripheral nodes for
practical purposes in PWR reactors.
Abstract: In this paper, the difference between the Alternating
Direction Method (ADM) and the Non-Splitting Method (NSM) is
investigated, while both methods applied to the simulations for 2-D
multimaterial radiation diffusion issues. Although the ADM have the
same accuracy orders with the NSM on the uniform meshes, the
accuracy of ADM will decrease on the distorted meshes or the
boundary of domain. Numerical experiments are carried out to
confirm the theoretical predication.
Abstract: Calcium [Ca2+] dynamics is studied as a potential form
of neuron excitability that can control many irregular processes like
metabolism, secretion etc. Ca2+ ion enters presynaptic terminal and
increases the synaptic strength and thus triggers the neurotransmitter
release. The modeling and analysis of calcium dynamics in neuron
cell becomes necessary for deeper understanding of the processes
involved. A mathematical model has been developed for cylindrical
shaped neuron cell by incorporating physiological parameters like
buffer, diffusion coefficient, and association rate. Appropriate initial
and boundary conditions have been framed. The closed form solution
has been developed in terms of modified Bessel function. A computer
program has been developed in MATLAB 7.11 for the whole
approach.
Abstract: Mathematical and computational modeling of calcium
signalling in nerve cells has produced considerable insights into how
the cells contracts with other cells under the variation of biophysical
and physiological parameters. The modeling of calcium signaling in
astrocytes has become more sophisticated. The modeling effort has
provided insight to understand the cell contraction. Main objective
of this work is to study the effect of voltage gated (Operated)
calcium channel (VOC) on calcium profile in the form of advection
diffusion equation. A mathematical model is developed in the form
of advection diffusion equation for the calcium profile. The model
incorporates the important physiological parameter like diffusion
coefficient etc. Appropriate boundary conditions have been framed.
Finite volume method is employed to solve the problem. A program
has been developed using in MATLAB 7.5 for the entire problem
and simulated on an AMD-Turion 32-bite machine to compute the
numerical results.
Abstract: The leaching rate of 137Cs from spent mix bead (anion and cation) exchange resins in a cement-bentonite matrix has been studied. Transport phenomena involved in the leaching of a radioactive material from a cement-bentonite matrix are investigated using three methods based on theoretical equations. These are: the diffusion equation for a plane source an equation for diffusion coupled to a firstorder equation and an empirical method employing a polynomial equation. The results presented in this paper are from a 25-year mortar and concrete testing project that will influence the design choices for radioactive waste packaging for a future Serbian radioactive waste disposal center.
Abstract: In the present work, the effective moisture diffusivity
and activation energy were calculated using an infinite series solution
of Fick-s diffusion equation. The results showed that increasing
drying temperature accelerated the drying process. All drying
experiments had only falling rate period. The average effective
moisture diffusivity values varied from 2.807x10-10 to 6.977x10-10m2
s_1 over the temperature and velocity range. The temperature
dependence of the effective moisture diffusivity for the thin layer
drying of the ginger slices was satisfactorily described by an
Arrhenius-type relationship with activation energy values of 19.313-
22.722 kJ.mol-1 within 40–70 °C and 0.8-3 ms-1 temperature range.
Abstract: Multiport diffusers are the effective engineering
devices installed at the modern marine outfalls for the steady
discharge of effluent streams from the coastal plants, such as
municipal sewage treatment, thermal power generation and seawater
desalination. A mathematical model using a two-dimensional
advection-diffusion equation based on a flat seabed and incorporating
the effect of a coastal tidal current is developed to calculate the
compounded concentration following discharges of desalination
brine from a sea outfall with multiport diffusers. The analytical
solutions are computed graphically to illustrate the merging of
multiple brine plumes in shallow coastal waters, and further
approximation will be made to the maximum shoreline's
concentration to formulate dilution of a multiport diffuser discharge.
Abstract: In this paper, a new time discontinuous expanded mixed finite element method is proposed and analyzed for two-order convection-dominated diffusion problem. The proofs of the stability of the proposed scheme and the uniqueness of the discrete solution are given. Moreover, the error estimates of the scalar unknown, its gradient and its flux in the L1( ¯ J,L2( )-norm are obtained.
Abstract: In this study, the density dependent nonlinear reactiondiffusion
equation, which arises in the insect dispersal models, is
solved using the combined application of differential quadrature
method(DQM) and implicit Euler method. The polynomial based
DQM is used to discretize the spatial derivatives of the problem. The
resulting time-dependent nonlinear system of ordinary differential
equations(ODE-s) is solved by using implicit Euler method. The
computations are carried out for a Cauchy problem defined by a onedimensional
density dependent nonlinear reaction-diffusion equation
which has an exact solution. The DQM solution is found to be in a
very good agreement with the exact solution in terms of maximum
absolute error. The DQM solution exhibits superior accuracy at large
time levels tending to steady-state. Furthermore, using an implicit
method in the solution procedure leads to stable solutions and larger
time steps could be used.
Abstract: In this paper we study some numerical methods to solve a model one-dimensional convection–diffusion equation. The semi-discretisation of the space variable results into a system of ordinary differential equations and the solution of the latter involves the evaluation of a matrix exponent. Since the calculation of this term is computationally expensive, we study some methods based on Krylov subspace and on Restrictive Taylor series approximation respectively. We also consider the Chebyshev Pseudospectral collocation method to do the spatial discretisation and we present the numerical solution obtained by these methods.