Numerical Study of Coal Gasification Using Eulerian-Eulerian Multiphase Model

Author(s):  
Shaoping Shi ◽  
Christopher Guenther ◽  
Stefano Orsino

Gasification converts the carbon-containing material into a synthesis gas (syngas) which can be used as a fuel to generate electricity or used as a basic chemical building block for a large number of uses in the petrochemical and refining industries. Based on the mode of conveyance of the fuel and the gasifying medium, gasification can be classified into fixed or moving bed, fluidized bed, and entrained flow reactors. Entrained flow gasifiers normally feature dilute flow with small particle size and can be successfully modeled with the Discrete Phase Method (DPM). For the other types, the Eulerian-Eulerian (E-E) or the so called two-fluid multiphase model is a more appropriate approach. The E-E model treats the solid phase as a distinct interpenetrating granular “fluid” and it is the most general-purposed multi-fluid model. This approach provides transient, three-dimensional, detailed information inside the reactor which would otherwise be unobtainable through experiments due to the large scale, high pressure and/or temperature. In this paper, a transient, three-dimensional model of the Power Systems Development Facility (PSDF) transport gasifier will be presented to illustrate how Computational Fluid Dynamics (CFD) can be used for large-scale complicated geometry with detailed physics and chemistry. In the model, eleven species are included in the gas phase while four pseudo-species are assumed in the solid phase. A total of sixteen reactions, both homogeneous (involving only gas phase species) and heterogeneous (involving species in both gas and solid phases), are used to model the coal gasification chemistry. Computational results have been validated against PSDF experimental data from lignite to bituminous coals under both air and oxygen blown conditions. The PSDF gasifier geometry was meshed with about 70,000, hexahedra-dominated cells. A total of six cases with different coal, feed gas, and/or operation conditions have been performed. The predicted and measured temperature profiles along the gasifier and gas compositions at the outlet agreed fairly well.

Energies ◽  
2021 ◽  
Vol 14 (5) ◽  
pp. 1233
Author(s):  
Umair Jamil Ur Rahman ◽  
Artur Krzysztof Pozarlik ◽  
Thomas Tourneur ◽  
Axel de Broqueville ◽  
Juray De Wilde ◽  
...  

In this paper, an intensified spray-drying process in a novel Radial Multizone Dryer (RMD) is analyzed by means of CFD. A three-dimensional Eulerian–Lagrangian multiphase model is applied to investigate the effect of solids outlet location, relative hot/cold airflow ratio, and droplet size on heat and mass transfer characteristics, G-acceleration, residence time, and separation efficiency of the product. The results indicate that the temperature pattern in the dryer is dependent on the solids outlet location. A stable, symmetric spray behavior with maximum evaporation in the hot zone is observed when the solids outlet is placed at the periphery of the vortex chamber. The maximum product separation efficiency (85 wt %) is obtained by applying high G-acceleration (at relative hot/cold ratio of 0.75) and narrow droplet size distribution (45–70 µm). The separation of different sized particles with distinct drying times is also observed. Smaller particles (<32 µm) leave the reactor via the gas outlet, while the majority of big particles leave it via the solids outlet, thus depicting in situ particle separation. The results revealed the feasibility and benefits of a multizone drying operation and that the RMD can be an attractive solution for spray drying technology.


2018 ◽  
Vol 619 ◽  
pp. A82
Author(s):  
Man Zhang ◽  
Yu Fen Zhou ◽  
Xue Shang Feng ◽  
Bo Li ◽  
Ming Xiong

In this paper, we have used a three-dimensional numerical magnetohydrodynamics model to study the reconnection process between magnetic cloud and heliospheric current sheet. Within a steady-state heliospheric model that gives a reasonable large-scale structure of the solar wind near solar minimum, we injected a spherical plasmoid to mimic a magnetic cloud. When the magnetic cloud moves to the heliospheric current sheet, the dynamic process causes the current sheet to become gradually thinner and the magnetic reconnection begin. The numerical simulation can reproduce the basic characteristics of the magnetic reconnection, such as the correlated/anticorrelated signatures in V and B passing a reconnection exhaust. Depending on the initial magnetic helicity of the cloud, magnetic reconnection occurs at points along the boundary of the two systems where antiparallel field lines are forced together. We find the magnetic filed and velocity in the MC have a effect on the reconnection rate, and the magnitude of velocity can also effect the beginning time of reconnection. These results are helpful in understanding and identifying the dynamic process occurring between the magnetic cloud and the heliospheric current sheet.


2000 ◽  
Author(s):  
M. Singh ◽  
P. K. Panigrahi ◽  
G. Biswas

Abstract A numerical study of rib augmented cooling of turbine blades is reported in this paper. The time-dependent velocity field around a pair of symmetrically placed ribs on the walls of a three-dimensional rectangular channel was studied by use of a modified version of Marker-And-Cell algorithm to solve the unsteady incompressible Navier-Stokes and energy equations. The flow structures are presented with the help of instantaneous velocity vector and vorticity fields, FFT and time averaged and rms values of components of velocity. The spanwise averaged Nusselt number is found to increase at the locations of reattachment. The numerical results are compared with available numerical and experimental results. The presence of ribs leads to complex flow fields with regions of flow separation before and after the ribs. Each interruption in the flow field due to the surface mounted rib enables the velocity distribution to be more homogeneous and a new boundary layer starts developing downstream of the rib. The heat transfer is primarily enhanced due to the decrease in the thermal resistance owing to the thinner boundary layers on the interrupted surfaces. Another reason for heat transfer enhancement can be attributed to the mixing induced by large-scale structures present downstream of the separation point.


2021 ◽  
Author(s):  
Vojtech Patocka ◽  
Nicola Tosi ◽  
Enrico Calzavarini

&lt;p&gt;We evaluate the equilibrium concentration of a thermally convecting suspension that is cooled from above and in which&lt;br&gt;solid crystals are self-consistently generated in the thermal boundary layer near the top. In a previous study (Patoc&amp;#780;ka et&lt;br&gt;al., 2020), we investigated the settling rate of solid particles suspended in a highly vigorous (Ra = 10&lt;sup&gt;8&lt;/sup&gt; , 10&lt;sup&gt;10&lt;/sup&gt;, and 10&lt;sup&gt;12&lt;/sup&gt; ),&lt;br&gt;finite Prandtl number (Pr = 10, 50) convection. In this follow-up study we additionally employ the model of crystal&lt;br&gt;generation and growth of Jarvis and Woods (1994), instead of using particles with a predefined size and density that are&lt;br&gt;uniformly injected into the carrier fluid.&lt;/p&gt;&lt;p&gt;We perform a series of numerical experiments of particle-laden thermal convection in 2D and 3D Cartesian geometry&lt;br&gt;using the freely available code CH4 (Calzavarini, 2019). Starting from a purely liquid phase, the solid fraction gradually&lt;br&gt;grows until an equilibrium is reached in which the generation of the solid phase balances the loss of crystals due to&lt;br&gt;sedimentation at the bottom of the fluid. For a range of predefined density contrasts of the solid phase with respect to&lt;br&gt;the density of the fluid (&amp;#961;&lt;sub&gt;p&lt;/sub&gt; /&amp;#961;&lt;sub&gt;f&lt;/sub&gt; = [0, 2]), we measure the time it takes to reach such equilibrium. Both this time and&lt;br&gt;the equilibrium concentration depend on the average settling rate of the particles and are thus non-trival to compute for&lt;br&gt;particle types that interact with the large-scale circulation of the fluid (see Patoc&amp;#780;ka et al., 2020).&lt;/p&gt;&lt;p&gt;We apply our results to the cooling of a large volume of magma, spanning from a large magma chamber up to a&lt;br&gt;global magma ocean. Preliminary results indicate that, as long as particle re-entrainment is not a dominant process, the&lt;br&gt;separation of crystals from the fluid is an efficient process. Fractional crystallization is thus expected and the suspended&lt;br&gt;solid fraction is typically small, prohibiting phenomena in which the feedback of crystals on the fluid begins to govern the&lt;br&gt;physics of the system (e.g. Sparks et al, 1993).&lt;/p&gt;&lt;p&gt;References&lt;br&gt;Patoc&amp;#780;ka V., Calzavarini E., and Tosi N.(2020). Settling of inertial particles in turbulent Rayleigh-Be&amp;#769;nard convection.&lt;br&gt;Physical Review Fluids, 26(4) 883-889.&lt;/p&gt;&lt;p&gt;Jarvis, R. A. and Woods, A. W.(1994). The nucleation, growth and settling of crystals from a turbulently convecting&lt;br&gt;fluid. J. Fluid. Mech, 273 83-107.&lt;/p&gt;&lt;p&gt;Sparks, R., Huppert, H., Koyaguchi, T. et al (1993). Origin of modal and rhythmic igneous layering by sedimentation in&lt;br&gt;a convecting magma chamber. Nature, 361, 246-249.&lt;/p&gt;&lt;p&gt;Calzavarini, E (2019). Eulerian&amp;#8211;Lagrangian fluid dynamics platform: The ch4-project. Software Impacts, 1, 100002.&lt;/p&gt;


Materials ◽  
2018 ◽  
Vol 11 (10) ◽  
pp. 1857 ◽  
Author(s):  
Swati Sharma

When certain polymers are heat-treated beyond their degradation temperature in the absence of oxygen, they pass through a semi-solid phase, followed by the loss of heteroatoms and the formation of a solid carbon material composed of a three-dimensional graphenic network, known as glassy (or glass-like) carbon. The thermochemical decomposition of polymers, or generally of any organic material, is defined as pyrolysis. Glassy carbon is used in various large-scale industrial applications and has proven its versatility in miniaturized devices. In this article, micro and nano-scale glassy carbon devices manufactured by (i) pyrolysis of specialized pre-patterned polymers and (ii) direct machining or etching of glassy carbon, with their respective applications, are reviewed. The prospects of the use of glassy carbon in the next-generation devices based on the material’s history and development, distinct features compared to other elemental carbon forms, and some large-scale processes that paved the way to the state-of-the-art, are evaluated. Selected support techniques such as the methods used for surface modification, and major characterization tools are briefly discussed. Barring historical aspects, this review mainly covers the advances in glassy carbon device research from the last five years (2013–2018). The goal is to provide a common platform to carbon material scientists, micro/nanomanufacturing experts, and microsystem engineers to stimulate glassy carbon device research.


Author(s):  
Ve´ronique Penin ◽  
Pascale Kulisa ◽  
Franc¸ois Bario

During the last few decades, the size and weight of turbo-machinery have been continuously reduced. However, by decreasing the distance between rows, rotor-stator interaction is strengthened. Two interactions now have the same magnitude: wake interaction and potential effect. Studying this effect is essential to understand rotor-stator interactions. Indeed, this phenomenon influences the whole flow, including the boundary layer of the upstream and downstream blades, ergo the stability of the flow and the efficiency of the machine. A large scale turbine cascade followed by a specially designed rotating cylinder system is used. Synchronised velocity LDA measurements on the vane profile show the flow and boundary layer behavior due to the moving bars. To help the general understanding and to corroborate our experimental results, numerical investigations are carried out with an unsteady three dimensional Navier-Stokes code. Moreover, the numerical study informs about the potential disturbance to the whole flow of the cascade.


Author(s):  
Rongbin Li ◽  
Mingzhuang Xie ◽  
Hui Jin ◽  
Liejin Guo ◽  
Fengqin Liu

AbstractThe three-dimensional (3-D) comprehensive mathematical model was developed to simulate the coal gasification process in an entrained flow gasifier with a swirl burner. The models employed or developed includes the coal devolatilization model, the char combustion and gasification model, the gas homogeneous reaction model, the random-trajectory model, gas turbulence model, and the P-1 radiation model. The solution of models was executed based on the computational fluid dynamics (CFD). By qualitatively comparing the results at different swirl number, the significant influences of swirl on characteristics of coal gasification such as flow distributions, gas temperature and product composition including hydrogen (H2), carbon monoxide (CO), etc., and on the performance of coal gasification such as averaged exit product composition, carbon conversion rate and cold gas efficiency, were in detail discussed. Especially, a proper swirl number (S ≤ 0.65) in favor of gasification was found for the investigated gasifier in this paper.


2015 ◽  
Vol 776 ◽  
pp. 161-199 ◽  
Author(s):  
Kyoungsik Chang ◽  
George Constantinescu

This numerical study investigates flow and turbulence structure through and around a circular array of solid circular cylinders of diameter $d$. The region containing the array of rigid cylinders resembles a porous circular cylinder of diameter $D$. The porous cylinder Reynolds number defined with the steady incoming flow velocity is $\mathit{Re}_{D}=10\,000$. Fully three-dimensional (3D) large eddy simulations (LES) are conducted to study the effects of the volume fraction of solids of the porous cylinder ($0.023<\text{SVF}<0.2$) and $d/D$ on the temporal variation and mean values of the drag/lift forces acting on the solid cylinders and on the porous cylinder. The effects of the bleeding flow through the circular porous cylinder on the wake structure and the influence of the SVF and $d/D$ on the onset of flow three-dimensionality within or downstream of the porous cylinder and transition to turbulence are discussed. Results are compared with experimental data, predictions of theoretical models available in the literature and also with the canonical case of a solid cylinder ($\text{SVF}=1,d/D=1$). Three-dimensional LES predict that large-scale wake billows are shed in the wake of the porous cylinder for $\text{SVF}>0.05$, similar to the von Karman vortex street observed for solid cylinders. As the SVF decreases, the length of the separated shear layers (SSLs) of the porous cylinder and the distance from the back of the porous cylinder at which wake billows form increase. For sufficiently low volume fractions of solids (e.g. $\text{SVF}=0.05$, 0.023), no wake billows are shed and the interactions among the wakes of the solid cylinders are weak. Even for $\text{SVF}=0.023$, SSLs containing large-scale turbulent eddies form on the two sides of the porous cylinder, but their ends cannot interact to generate wake billows. In both regimes, the force acting on some of the solid cylinders within the array is highly unsteady. As opposed to results obtained based on 2D simulations, no intermediate regime in which the force acting on the solid cylinders is close to steady is present. Interestingly, an energetic low frequency corresponding to a Strouhal number defined with the diameter of the porous cylinder of approximately 0.2 is present within the porous cylinder and near-wake regions not only for cases where wake billows are generated but also for cases where no wake billows form. In the latter cases, this frequency is due to an instability acting on the SSLs which induces in-phase large-scale undulatory deformations of the two SSLs. A combined drag parameter for the porous cylinder ${\it\Gamma}_{D}=\overline{C}_{d}\,aD/(1-\text{SVF})$ is introduced, where $aD$ is the non-dimensional frontal area per unit volume of the porous cylinder. This parameter characterizes by how much the velocity of the bleeding flow at the back of the porous cylinder is reduced compared with the incoming flow velocity for a given total drag force acting on the porous cylinder. Results from simulations conducted with different values of the SVF, $d/D$ and mean time-averaged solid cylinder streamwise drag parameter, $\overline{C}_{d}$, show that ${\it\Gamma}_{D}$ increases monotonically with increasing $aD$. Several ways of defining the spatial extent of the wake region in a less ambiguous way are proposed.


Author(s):  
Alexander P. Skuratov ◽  
Alexander V. Ivlev ◽  
Artem A. Pianykh

A three-dimensional mathematical model of the solidification process of a liquid metal is considered, taking into account the mobility of the boundaries at which the phase transition is carried out (Stefan boundary value problem). The algorithm of calculation is improved, allowing due to the use of the Dirac δ-function in determining the effective heat capacity to take into account the nonlinearity of the equation of unsteady thermal conductivity and the heat of the phase transition. A numerical study of heat transfer during solidification of lead-containing aluminum melt droplets in air and water is carried out. The influence of droplet size and melt overheating on the solidification dynamics of granules has been studied. An approximate ratio based on the square root law is proposed, taking into account the amount of overheating of the liquid phase and linking the thickness of the formed solid phase with the duration of the granulation process


Sign in / Sign up

Export Citation Format

Share Document