Simulation of Coal Combustion by AUSM Turbulence-Chemistry Char Combustion Model and a Full Two-Fluid Model

An algebraic unified second-order moment (AUSM) turbulence-chemistry model of char combustion is introduced in this paper, to calculate the effect of particle temperature fluctuation on char combustion. The AUSM model is used to simulate gas-particle flows, in coal combustion in a pulverized coal combustor, together with a full two-fluid model for reacting gas-particle flows and coal combustion, including the sub-models as the k-epsilon-k(p) two-phase turbulence niodel, the EBU-Arrhenius volatile and CO combustion model, and the six-flux radiation model. A new method for calculating particle mass flow rate is also used in this model to correct particle outflow rate and mass flow rate for inside sections, which can obey the principle of mass conservation for the particle phase and can also speed up the iterating convergence of the computation procedure effectively. The simulation results indicate that, the AUSM char combustion model is more preferable to the old char combustion model, since the later totally eliminate the influence of particle temperature fluctuation on char combustion rate.

Two-dimensional fluid model simulation of bell jar top inductively coupled plasma

In the present paper, argon (Ar) plasmas in a bell jar inductively coupled plasma (ICP) source are systematically studied over pressures from 5 to 20 mtorr and power inputs from 0.2 to 0.5 kW. In this study, both a two-dimensional (2-D) fluid model simulation and global model calculation are compared, The 2-D fluid model simulation with a self-consistent power deposition is developed to describe the Ar plasma behavior as well as predict the plasma parameter distributions, Finally, a quantitative comparison between the global model and the fluid model is made to test their validity.

The interaction of two pairs of solitary waves in a two-fluid system

This paper deals with the interaction of solitary waves in a two-fluid system which consistsof two superimposed incompressible inviscid fluids with a free surface and a horizontal rigidbottom. Under the assumption of shallow water wave, we first derive the basic equationssuitable for the model considered, a generalized form of the Boussinesq equations, then usingthe PLK method and the reductive perturbation method, obtain the second-order approximatesolution for the head-on collision between two pairs of interface and surface solitary waves,and give their maximum amplitudes during the collision and the nonuniform phase shiftsafter the collision which lead to the distortion of the wave profiles.

Higher-order Boussinesq-type equations for interfacial waves in a two-fluid system

Interfacial waves propagating along the interface between a three-dimensional two-fluid system with a rigid upper boundary and an uneven bottom are considered. There is a light fluid layer overlying a heavier one in the system, and a small density difference exists between the two layers. A set of higher-order Boussinesq-type equations in terms of the depth-averaged velocities accounting for stronger nonlinearity are derived. When the small parameter measuring frequency dispersion keeping up to lower-order and full nonlinearity are considered, the equations include the Choi and Camassa's results (1999). The enhanced equations in terms of the depth-averaged velocities are obtained by applying the enhancement technique introduced by Madsen et al. (1991) and Schaffer and Madsen (1995a). It is noted that the equations derived from the present study include, as special cases, those obtained by Madsen and Schaffer (1998). By comparison with the dispersion relation of the linear Stokes waves, we found that the dispersion relation is more improved than Choi and Camassa's (1999) results, and the applicable scope of water depth is deeper.

Shock-Induced Near-Wall Two-Phase Flow Structure Over a Micron-Sized Particles Bed

The present paper studies numerical modelling of near-wall two-phase flows induced by a normal shock wave moving at a constant speed, over a micronsized particles bed. In this two-fluid model, the possibility of particle trajectory intersection is considered and a full Lagrangian formulation of the dispersed phase is introduced. The finiteness of the Reynolds and Mach numbers of the flow around a particle as well as the fineness of the particle sizes are taken into account in describing the interactions between the carrier- and dispersed- phases. For the small mass-loading ratio case, the numerical simulation of flow structure of the two phases is implemented and the profiles of the particle number density are obtained under the constant-flux condition on the wall. The effects of the shock Mach number and the particle size and material density on particle entrainment motion are discussed in detail.The obtained results indicate that interphase non-equilibrium in the velocity and temperature is a common feature for this type of flows and a local particle accumulation zone may form near the envelope of the particle trajectory family.

Oil-Water Two-Phase Flow Inside T-Junction

The oil/water two-phase flow inside T-junctions was numerically simulated with a 3-D two-fluid model, and the turbulence was described using the mixture k - epsilon model. Some experiments of oil/water flow inside a single T-junction were conducted in the laboratory. The results show that the separating performance of T-junction largely depends oil the inlet volumetric fraction and flow patterns. A reasonable agreement is reached between the numerical simulation and the experiments for both the oil fraction distribution and the separation efficiency.

Study of drag reduction by gas injection for power-law fluid flow in horizontal stratified and slug flow regimes

In this work, the drag reduction by gas injection for power-law fluid flow in stratified and slug flow regimes has been studied. Experimentswere conducted to measure the pressure gradient within air/CMC solutions in a horizontal Plexiglas pipe that had a diameter of 50mm and a length of 30 m. The drag reduction ratio in stratified flow regime was predicted using the two-fluid model. The results showed that the drag reduction should occur over the large range of the liquid holdup when the flow behaviour index remained at the low value. Furthermore, for turbulent gas-laminar liquid stratified flow, the drag reduction by gas injection for Newtonian fluid was more effective than that for shear-shinning fluid, when the dimensionless liquid height remained in the area of high value. The pressure gradient model for a gas/Newtonian liquid slug flow was extended to liquids possessing the Ostwald–de Waele power law model. The proposed model was validated against 340 experimental data point over a wide range of operating conditions, fluid characteristics and pipe diameters. The dimensionless pressure drop predicted was well inside the 20% deviation region for most of the experimental data. These results substantiated the general validity of the model presented for gas/non-Newtonian two-phase slug flows.

Characteristics of a CW HF chemical-laser calculated from a simplified two-dimensional model

A two-dimensional simplified model of an HF chemical laser is introduced. Using an implicit finite difference scheme, the solution of two adjacent parallel streams with diffusion mixing and chemical reaction is generated. A contour of mixing and reaction boundary is obtained without presupposition. The distribution of the HF(v) concentrations, gas temperature and the optical small signal gain (alpha sub V, J) on the flowing plane (X, Y) are presented. Compared with the solution solved directly from a set of Navier-Stokes equations, the results of these two methods agree with each other qualitatively. The influences of the different velocity, temperature (T sub 0) and composition of the two streams on the small signal gain after the nozzle exit are investigated. It is interesting that for larger J with a fixed v, the peaks of alpha sub v-T sub 0 profiles move towards higher T sub 0. The computing method is simple and only a short computing time is needed.

An Improved CE/SE Scheme for Numerical Simulation of Gaseous and Two-Phase Detonations

A new structure of solution elements and conservation elements based on rectangular mesh was pro- posed and an improved space-time conservation element and solution element (CE/SE) scheme with sec- ond-order accuracy was constructed. Furthermore, the application of improved CE/SE scheme was extended to detonation simulation. Three models were used for chemical reaction in gaseous detonation. And a two-fluid model was used for two-phase (gas–droplet) detonation. Shock reflections were simu- lated by the improved CE/SE scheme and the numerical results were compared with those obtained by other different numerical schemes. Gaseous and gas–droplet planar detonations were simulated and the numerical results were carefully compared with the experimental data and theoretical results based on C–J theory. Mach reflection of a cellular detonation was also simulated, and the numerical cellular pat- terns were compared with experimental ones. Comparisons show that the improved CE/SE scheme is clear in physical concept, easy to be implemented and high accurate for above-mentioned problems.

Two-component model for the interstellar scattering in direction to the Vela pulsar

This paper discusses a rigorous treatment of the refractive scintillation of pulsar PSR B0833-45 caused by a two-component interstellar scattering medium. It is assumed that the interstellar scattering medium is composed of a thin screen ISM and an extended interstellar medium. We consider that the scattering of the thin screen concentrates in a thin layer presented by a delta function distribution and that the scattering density of the extended irregular medium satisfies the Gaussian distribution. We investigate and develop equations for the flux density structure function corresponding to this two-component ISM geometry in the scattering density distribution and compare our result with that of the Vela pulsar observations. We conclude that the refractive scintillation caused by this two-component ISM scattering gives a more satisfactory explanation for the observed flux density variation of the Vela pulsar than does the single extended medium model. The level of refractive scintillation is strongly sensitive to the distribution of scattering material along the line of sight. The logarithmic slope of the structure function is sensitive to thin screen location and is relatively insensitive to the scattering strength of the thin screen medium. Therefore, the proposed model can be applied to interpret the structure function of flux density observed in pulsar PSR B0833-45. The result suggests that the medium consists of a discontinuous distribution of plasma turbulence embedded in the Vela supernova remnant. Thus our work provides some insight into the distribution of the scattering along the line of sight to the Vela pulsar.

A new flow field and its two-dimension model for polymer electrolyte membrane fuel cells (PEMFCs)

A new flow field was designed to search flow fields fitting polymer electrolyte membrane fuel cells (PEMFCs) better due its extensible. There are many independent inlets and outlets in the new flow field. The new flow field we named NINO can extend to be more general when pressures at the inlet and outlet vary and some usual flow fields will be obtained. A new mathematical model whose view angle is obverse is used to describe the flow field.

Theoretical estimation of thermodynamic properties of the system PS/PPO on the basis of modified combining rule of Sanchez-Lacombe lattice fluid model

The three scaling parameters described in Sanchez-Lacombe lattice fluid theory (SLLFT), T*, P* and rho* of pure polystyrene (PS), pure poly(2,6-dimethyl-1,4-phenylene oxide) (PPO) and their mixtures are obtained by fitting corresponding experimental pressure volume-temperature data with equation-of-state of SLLFT. A modified combining rule in SLLFT used to match the volume per mer, v* of the PS/PPO mixtures was advanced and the enthalpy of mixing and Flory-Huggins (FH) interaction parameter were calculated using the new rule. It is found that the difference between the new rule and the old one presented by Sanchez and Lacombe is quite small in the calculation of the enthalpy of mixing and FH interaction parameter and the effect of volume-combining rule on the calculation of thermodynamic properties is much smaller than that of energy-combining rule. But the relative value of interaction parameter changes much due to the new volume-based combining rule. This effect can affect the position of phase diagram very much, which is reported elsewhere [Macromolecules 34 (2001) 6291]

Statistical thermodynamics of polydisperse polymer systems in the framework of lattice fluid model: Effect of molecular weight and its distribution on the spinodal in polymer solution

With the aid of thermodynamics of Gibbs, the expression of the spinodal was derived for the polydisperse polymer-solvent system in the framework of Sanchez-Lacombe Lattice Fluid Theory (SLLFT). For convenience, we considered that a model polydisperse polymer contains three sub-components. According to our calculation, the spinodal depends on both weight-average ((M) over bar (w)) and number-average ((M) over bar (n)) molecular weights of the polydisperse polymer, but the z-average molecular weight ((M) over bar (z)) dependence on the spinodal is invisible. The dependence of free volume on composition, temperature, molecular weight, and its distribution results in the effect of (M) over bar (n) on the spinodal. Moreover, it has been found that the effect of changing (M) over bar (w) on the spinodal is much bigger than that of changing (M) over bar (n) and the extrema of the spinodal increases with the rise of the weight-average molecular weight of the polymer in the solutions with upper critical solution temperature (UCST). However, the effect of polydispersity on the spinodal can be neglected for the polymer with a considerably high weight-average molecular weight. A more simple expression of the spinodal for the polydisperse polymer solution in the framework of SLLFT was also derived under the assumption of upsilon(*)=upsilon(1)(*)=upsilon(2)(*) and (1/r(1)(0))-(1/r(2i)(0))-->(1/r(1)(0)).