Use of Different Numerical Solution Approaches for a Three-Phase Slurry Catalytic Reactor Model

Author(s):  
Mardonny Nazareno Barreira ◽  
Eduardo Coselli Vasco de Toledo ◽  
Rubens Maciel Filho ◽  
Maria das Graças Enrique

In this paper two different numerical approaches are used to find out the steady-state solution for multiphase chemical reactors. Such systems are distributed parameter systems and the steady-state formulation problem leads to a system of nonlinear algebraic equations. The problem is to solve the equations in a robust way so that reliable predictions can be made. Bearing this in mind in this work two methods are implement and their performances compared, to know the Orthogonal Collocation Method e the Finite Volume Method. The results show that in spite of good qualitative agreement in the predictions the latter approaches in more suitable to have a direct solution of the steady-state model of the reactor.

IEEE Access ◽  
2016 ◽  
Vol 4 ◽  
pp. 9479-9488 ◽  
Author(s):  
Yuntao Ju ◽  
Fuchao Ge ◽  
Wenchuan Wu ◽  
Yi Lin ◽  
Jing Wang

1989 ◽  
Vol 111 (2) ◽  
pp. 253-259 ◽  
Author(s):  
Y. S. Choi ◽  
S. T. Noah

The nonlinear, steady-state response of a displacement-forced link coupling with clearance with finite stiffness is determined. The solution procedure is derived from satisfying the boundary conditions at the contact points and then solving the resulting nonlinear algebraic equations by setting the duration of contact as a parameter. This direct approach to determining periodic solutions for systems with clearances with finite stiffness is substantially more efficient than numerical integration schemes. Results in terms of contact forces and durations of contact are pertinent to fatigue and wear studies. Parametric relations are presented for effects of the variation of damping, stiffness, exciting displacement, and gap length on the dynamic behavior of the link pair.


2009 ◽  
Vol 3 (2) ◽  
pp. 217-229 ◽  
Author(s):  
T. Zwinger ◽  
J. C. Moore

Abstract. We present steady state (diagnostic) and transient (prognostic) simulations of Midtre Lovénbreen, Svalbard performed with the thermo-mechanically coupled full-Stokes code Elmer. This glacier has an extensive data set of geophysical measurements available spanning several decades, that allow for constraints on model descriptions. Consistent with this data set, we included a simple model accounting for the formation of superimposed ice. Diagnostic results indicated that a dynamic adaptation of the free surface is necessary, to prevent non-physically high velocities in a region of under determined bedrock depths. Observations from ground penetrating radar of the basal thermal state agree very well with model predictions, while the dip angles of isochrones in radar data also match reasonably well with modelled isochrones, despite the numerical deficiencies of estimating ages with a steady state model. Prognostic runs for 53 years, using a constant accumulation/ablation pattern starting from the steady state solution obtained from the configuration of the 1977 DEM show that: 1 the unrealistic velocities in the under determined parts of the DEM quickly damp out; 2 the free surface evolution matches well measured elevation changes; 3 the retreat of the glacier under this scenario continues with the glacier tongue in a projection to 2030 being situated ≈500 m behind the position in 1977.


2009 ◽  
Vol 63 (2) ◽  
Author(s):  
Marcel Kotora ◽  
Zuzana Švandová ◽  
Jozef Markoš

AbstractNonequilibrium model for steady state simulation of catalytic distillation is presented. Mathematical model takes into account both mass and heat transfers across the gas liquid interface and through the liquid-solid (catalyst) interface. Equations describing the mentioned phenomena are based on the effective diffusivity approach. The resulting system of nonlinear algebraic equations was implemented in the FORTRAN programming language and solved by the BUNLSI (Ferraris & Tronconi, 1986) solver. The described model was verified using the experimental data obtained from a continuous distillation column equipped with catalytic packing. As an experimental model system, synthesis of propyl propionate from propan-1-ol and propionic acid was chosen. Comparison of experimental and simulation data is presented, and appropriateness of the developed model for other types of catalytic distillation processes is discussed.


1997 ◽  
Vol 15 (1) ◽  
pp. 53-64 ◽  
Author(s):  
N.A. Inogamov

The work begins with a presentation of results concerning periodic Rayleigh-Taylor (RT)and Richtmyer–Meshkov (RM) dynamics. The periodic RT flows are considered in 2D and 3D cases. The periodic RM flows are considered in a 2D case. Nonstationary and steadystate statements of RT and RM problems are given. Nonstationary statements of RT and RM problems are given in the2D case. Steady-state statements of RT and RM problems are considered in the 2D and 3D cases (RT) and in the 2D case (RM). Questions about turbulence are given in the second part of the work. In the nonstationary statement, the beginning of a transition from an equilibrium close to hydrostatics and behavior of a dynamic system when it is close to the steady-state solution are considered. To study the smooth evolution of the dynamic system in the RT case (at last during an initial of order of unity t ≃ 1 stage) we have to cut of f an inverse turbulent cascade incoming from the small scales (inverse means that the typical k decreases in time during the cascade). The cascade starts from some small scale k ≫ 1 (here g = 1) and appears in the scales k ≃ 1 after time t ≃ 1 after start. We have to consider exponentially smoothed initial data to cut it off. This means that amplitudes of initial spectrum of perturbations have to be exponentially decreased in the direction of high k. The steadystates have been studied in two ways. In both of them the spectral decomposition of a velocity field and an elevation of the contact surface has been used. In one of them, corresponding to the nonstationary situation, the ordinary differential system for time dependent spectral components has been integrated. The other solutions of the algebraic equations for constant components have been studied. These results have been compared with a computer simulation that has been done by large particles and artificial compressibility methods. High-order nonlinear terms corresponding to the interaction of a wide spectrum of harmonics are considered in both ways. The spatial structure of the many dimensional periodic solutions is studied. It is known that there are two important features in the flows under consideration: the bubble envelope and the jet “phase.” Spatial structures of the 2D and 3D flows are compared in the report. It is shown that the apices of the bubbles are isolated points and the bubble envelopes in both flows are qualitatively similar; the difference has a quantitative character. On the contrary, the jet “phases” differ qualitatively in the 2D and 3D cases. In the 2D case in the plane picture, the jets are separated from each other. In this picture they are isolated fingers. In the 3D case, near the bubble envelope after some time of development following after cos kx cos qy type initial perturbation, they are not the isolated fingers. They form a structure of the intersecting walls. We see, thus, that the 2D and 3D cases are topologically different. Most powerful ejections of dense fluid are in the points where these wall jets intersect each other.


2017 ◽  
Author(s):  
Bayu Triwibowo ◽  
Haniif Prasetiawan ◽  
Anwaruddin Hisyam ◽  
Mohammad Fariz Fauzan ◽  
Muhammad Habib Fahd Rizky

Sign in / Sign up

Export Citation Format

Share Document