scholarly journals Mathematical Modeling of Reactive Extraction of Solute from Slab Solid Material

2020 ◽  
Vol 20 (2) ◽  
pp. 458
Author(s):  
Indah Hartati ◽  
Hary Sulistyo ◽  
Wahyudi Budi Sediawan ◽  
Muhammad Mufti Azis ◽  
Moh Fahrurrozi

Reactive extraction is gaining higher attention due its wide application in various solute separation processes. Here, a mathematical model of reactive extraction in slab has been proposed. The model was developed by considering simultaneous processes of active compound intra particle diffusion, second order elemental reaction of solute-active compound, and intra-particle product diffusion. The obtained partial differential equations (PDEs) were solved using Finite Difference Approximation (FDA) method by using realistic parameters. Concentration profile as well as product yield were evaluated as a function of time. As a result, the model proposed here may serve as a basis design for reactive extraction unit. Sensitivity analyses was conducted to inspect the influence of slab thickness, diffusivity and reaction rate constant to the product yield. Eventually, model validation was conducted by comparing the simulation results with analytical solutions for special cases. Validation results showed that the model gave good agreement with the analytical solution.

2020 ◽  
Vol 0 (0) ◽  
Author(s):  
Indah Hartati ◽  
Wahyudi Budi Sediawan ◽  
Hary Sulistyo ◽  
Muhammad Mufti Azis ◽  
Moh Fahrurrozi

AbstractHydrotropes have been largely explored as reactive extraction agent for lignin separation. In this paper, a mathematical model of hydrotropic-reactive extraction of sugarcane bagasse lignin was proposed and validated by experimental data from literature. The mathematical model was developed by assuming the particle is in slab shape, and by considering simultaneous processes of hydrotrope intra particle diffusion, second order reaction of lignin-hydrotrope, and intra-particle soluble delignification product diffusion. The proposed model results in a set of partial differential equations which were then solved by explicit finite difference approximation method. The mathematical model parameters were determined by fitting the model to the hydrotropic reactive extraction experimental data reported by Ansari and Gaikar (2014). Simulations show that the mathematical model of the hydrotropic-reactive extraction were well fitted to the experimental data with the obtained hydrotrope effective diffusivity (DeA) of 5.0 × 10−11 m2/s, effective diffusivity of soluble lignin product (DeC) of 9.0 × 10−12 m2/s and reaction rate constant (kr) of 1.78 × 10−10 m3/(g.s). It was also observed that the reaction was first order to the hydrotrope (n = 1), and one half order to the lignin (m = 0.5). Meanwhile the pseudo-stoichiometric mass ratio of hydrotrope to lignin was 6.4 g hydrotrope/g lignin.


1973 ◽  
Vol 40 (1) ◽  
pp. 7-12 ◽  
Author(s):  
W. H. Yang ◽  
C. H. Lu

A set of three nonlinear partial-differential equations is derived for general finite deformations of a thin membrane. The material that composes the membrane is assumed to be hyperelastic. Its mechanical property is represented by the neo-Hookean strain-energy function. The equations reduce to special cases known in the literature. A fast convergent algorithm is developed. The numerical solutions to the finite-difference approximation of the differential equations are computed iteratively with a trivial initial iterant. As an example, the problem of inflating a rectangular membrane with fixed edges by a uniform pressure applied on one side is presented. The solutions and their convergence are displayed and discussed.


2012 ◽  
Vol 12 (1) ◽  
pp. 193-225 ◽  
Author(s):  
N. Anders Petersson ◽  
Björn Sjögreen

AbstractWe develop a stable finite difference approximation of the three-dimensional viscoelastic wave equation. The material model is a super-imposition of N standard linear solid mechanisms, which commonly is used in seismology to model a material with constant quality factor Q. The proposed scheme discretizes the governing equations in second order displacement formulation using 3N memory variables, making it significantly more memory efficient than the commonly used first order velocity-stress formulation. The new scheme is a generalization of our energy conserving finite difference scheme for the elastic wave equation in second order formulation [SIAM J. Numer. Anal., 45 (2007), pp. 1902-1936]. Our main result is a proof that the proposed discretization is energy stable, even in the case of variable material properties. The proof relies on the summation-by-parts property of the discretization. The new scheme is implemented with grid refinement with hanging nodes on the interface. Numerical experiments verify the accuracy and stability of the new scheme. Semi-analytical solutions for a half-space problem and the LOH.3 layer over half-space problem are used to demonstrate how the number of viscoelastic mechanisms and the grid resolution influence the accuracy. We find that three standard linear solid mechanisms usually are sufficient to make the modeling error smaller than the discretization error.


1966 ◽  
Vol 6 (03) ◽  
pp. 217-227 ◽  
Author(s):  
Hubert J. Morel-Seytoux

Abstract The influence of pattern geometry on assisted oil recovery for a particular displacement mechanism is the object of investigation in this paper. The displacement is assumed to be of unit mobility ratio and piston-like. Fluids are assumed incompressible and gravity and capillary effects are neglected. With these assumptions it is possible to calculate by analytical methods the quantities of interest to the reservoir engineer for a great variety of patterns. Specifically, this paper presentsvery briefly, the methods and mathematical derivations required to obtain the results of engineering concern, andtypical results in the form of graphs or formulae that can be used readily without prior study of the methods. Results of this work provide checks for solutions obtained from programmed numerical techniques. They also reveal the effect of pattern geometry and, even though the assumptions of piston-like displacement and of unit mobility ratio are restrictive, they can nevertheless be used for rather crude but quick, cheap estimates. These estimates can be refined to account for non-unit mobility ratio and two-phase flow by correlating analytical results in the case M=1 and the numerical results for non-Piston, non-unit mobility ratio displacements. In an earlier paper1 it was also shown that from the knowledge of closed form solutions for unit mobility ratio, quantities called "scale factors" could be readily calculated, increasing considerably the flexibility of the numerical techniques. Many new closed form solutions are given in this paper. INTRODUCTION BACKGROUND Pattern geometry is a major factor in making water-flood recovery predictions. For this reason many numerical schemes have been devised to predict oil recovery in either regular patterns or arbitrary configurations. The numerical solutions, based on the method of finite difference approximation, are subject to errors often difficult to evaluate. An estimate of the error is possible by comparison with exact solutions. To provide a variety of checks on numerical solutions, a thorough study of the unit mobility ratio displacement process was undertaken. To calculate quantities of interest to the reservoir engineer (oil recovery, sweep efficiency, etc.), it is necessary to first know the pressure distribution in the pattern. Then analytical procedures are used to calculate, in order of increasing difficulty: injectivity, breakthrough areal sweep efficiency, normalized oil recovery and water-oil ratio as a function of normalized PV injected. BACKGROUND Pattern geometry is a major factor in making water-flood recovery predictions. For this reason many numerical schemes have been devised to predict oil recovery in either regular patterns or arbitrary configurations. The numerical solutions, based on the method of finite difference approximation, are subject to errors often difficult to evaluate. An estimate of the error is possible by comparison with exact solutions. To provide a variety of checks on numerical solutions, a thorough study of the unit mobility ratio displacement process was undertaken. To calculate quantities of interest to the reservoir engineer (oil recovery, sweep efficiency, etc.), it is necessary to first know the pressure distribution in the pattern. Then analytical procedures are used to calculate, in order of increasing difficulty: injectivity, breakthrough areal sweep efficiency, normalized oil recovery and water-oil ratio as a function of normalized PV injected.


Sign in / Sign up

Export Citation Format

Share Document