Modelling and validation of wet steam flow in a low pressure steam turbine

Results of numerical investigations of the wet steam flow in a three stage low pressure steam turbine test rig are presented. The test rig is a scale model of a modern steam turbine design and provides flow measurements over a range of operating conditions which are used for detailed comparisons with the numerical results. For the numerical analysis a modern CFD code with user defined models for specific wet steam modelling is used. The effect of different theoretical models for nucleation and droplet growth are examined. It is shown that heterogeneous condensation is highly dependent on steam quality and, in this model turbine with high quality steam, a homogeneous theory appears to be the best choice. The homogeneous theory gives good agreement between the test rig traverse measurements and the numerical results. The differences in the droplet size distribution of the three stage turbine are shown for different loads and modelling assumptions. The different droplet growth models can influence the droplet size by a factor of two. An estimate of the influence of unsteady effects is made by means of an unsteady two-dimensional simulation. The unsteady modelling leads to a shift of nucleation into the next blade row. For the investigated three stage turbine the influence due to wake chopping on the condensation process is weak but to confirm this conclusion further investigations are needed in complete three dimensions and on turbines with more stages. Copyright © 2011 by ASME.

Effects and numerical analysis of argon gas flow on the oxygen concentration in Czochralski silicon single crystal growth

Argon gas, as a protective environment and carrier of latent heat, has an important effect on the temperature distribution in crystals and melts. Numeric simulation is a potent tool for solving engineering problems. In this paper, the relationship between argon gas flow and oxygen concentration in silicon crystals was studied systematically. A flowing stream of argon gas is described by numeric simulation for the first time. Therefore, the results of experiments can be explained, and the optimum argon flow with the lowest oxygen concentration can be achieved. (C) 2002 Elsevier Science B.V. All rights reserved.

Operational Optimization of GaN Thin Film Growth Employing Numerical Simulation in a Showerhead MOCVD Reactor

A detailed reaction-tran sport model was studied in a showerhead reactor for metal organic chemical vapor deposition of GaN film by using computational fluid dynamics simulation. It was found that flat flow lines without swirl are crucial to improve the uniformity of the film growth, and thin temperature gradient above the suscptor can increase the film deposition rate. By above-mentioned research, we can employ higher h (the distance from the susceptor to the inlet), P (operational pressure) and the rate of susceptor rotation to improve the film growth.

The Deposition and Burning Characteristics During Slagging Co-Firing Coal and Wood: Modeling and Numerical Simulation

Numerical analysis was used to study the deposition and burning characteristics of combining co-combustion with slagging combustion technologies in this paper. The pyrolysis and burning kinetic models of different fuels were implanted into the WBSF-PCC2 (wall burning and slag flow in pulverized co-combustion) computation code, and then the slagging and co-combustion characteristicsespecially the wall burning mechanism of different solid fuels and their effects on the whole burning behavior in the cylindrical combustor at different mixing ratios under the condition of keeping the heat input samewere simulated numerically. The results showed that adding wood powder at 25% mass fraction can increase the temperature at the initial stage of combustion, which is helpful to utilize the front space of the combustor. Adding wood powder at a 25% mass fraction can increase the reaction rate at the initial combustion stage; also, the coal ignitability is improved, and the burnout efficiency is enhanced by about 5% of suspension and deposition particles, which is helpful for coal particles to burn entirely and for combustion devices to minimize their dimensions or sizes. The results also showed that adding wood powder at a proper ratio is helpful to keep the combustion stability, not only because of the enhancement for the burning characteristics, but also because the running slag layer structure can be changed more continuously, which is very important for avoiding the abnormal slag accumulation in the slagging combustor. The theoretic analysis in this paper proves that unification of co-combustion and slagging combustion technologies is feasible, though more comprehensive and rigorous research is needed.

Effects and numerical analysis of argon gas flow on the oxygen concentration in Czochralski silicon single crystal growth

Argon gas, as a protective environment and carrier of latent heat, has an important effect on the temperature distribution in crystals and melts. Numeric simulation is a potent tool for solving engineering problems. In this paper, the relationship between argon gas flow and oxygen concentration in silicon crystals was studied systematically. A flowing stream of argon gas is described by numeric simulation for the first time. Therefore, the results of experiments can be explained, and the optimum argon flow with the lowest oxygen concentration can be achieved. (C) 2002 Elsevier Science B.V. All rights reserved.

Numerical simulation of a bubble rising in shear-thinning fluids

The motion of a single bubble rising freely in quiescent non-Newtonian viscous fluids was investigated experimentally and computationally. The non-Newtonian effects in the flow of viscous inelastic fluids are modeled by the Carreau theological model. An improved level set approach for computing the incompressible two-phase flow with deformable free interface is used. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The simulation results demonstrate that the algorithm is robust for shear-thinning liquids with large density (rho(1)/rho(g) up to 10(3)) and high viscosity (eta(1)/eta(g) up to 10(4)). The comparison of the experimental measurements of terminal bubble shape and velocity with the computational results is satisfactory. It is shown that the local change in viscosity around a bubble greatly depends on the bubble shape and the zero-shear viscosity of non-Newtonian shear-thinning liquids. The shear-rate distribution and velocity fields are used to elucidate the formation of a region of large viscosity at the rear of a bubble as a result of the rather stagnant flow behind the bubble. The numerical results provide the basis for further investigations, such as the numerical simulation of viscoelastic fluids. (C) 2010 Elsevier B.V. All rights reserved.

Direct Numerical Simulation of Shock/Turbulent Boundary Layer Interaction in a Supersonic Compression Ramp

A direct numerical simulation of the shock/turbulent boundary layer interaction flow in a supersonic 24-degree compression ramp is conducted with the free stream Mach number 2.9. The blow-and-suction disturbance in the upstream wall boundary is used to trigger the transition. Both the mean wall pressure and the velocity profiles agree with those of the experimental data, which validates the simulation. The turbulent kinetic energy budget in the separation region is analyzed. Results show that the turbulent production term increases fast in the separation region, while the turbulent dissipation term reaches its peak in the near-wall region. The turbulent transport term contributes to the balance of the turbulent conduction and turbulent dissipation. Based on the analysis of instantaneous pressure in the downstream region of the mean shock and that in the separation bubble, the authors suggest that the low frequency oscillation of the shock is not caused by the upstream turbulent disturbance, but rather the instability of separation bubble.

FLOW-PIPE-SEEPAGE COUPLING ANALYSIS OF SPANNING INITIATION OF A PARTIALLY-EMBEDDED PIPELINE

The initiation of pipeline spanning involves the coupling between the flow over the pipeline and the seepage-flow in the soil underneath the pipeline. The pipeline spanning initiation is experimentally observed and discussed in this article. It is qualitatively indicated that the pressure-drop induced soil seepage failure is the predominant cause for pipeline spanning initiation. A flow-pipe-seepage sequential coupling Finite Element Method (FEM) model is proposed to simulate the coupling between the water flow-field and the soil seepage-field. A critical hydraulic gradient is obtained for oblique seepage failure of the sand in the direction tangent to the pipe. Parametric study is performed to investigate the effects of inflow velocity, pipe embedment on the pressure-drop, and the effects of soil internal friction angle and pipe embedment-to-diameter ratio on the critical flow velocity for pipeline spanning initiation. It is indicated that the dimensionless critical flow velocity changes approximately linearly with the soil internal friction angle for the submarine pipeline partially-embedded in a sandy seabed.

Numerical Simulation of a Low-power Hydrazine Arcjet Thruster

A modeling study is conducted to investigate the plasma flow and heat transfer characteristics of low-power (kW class) arc-heated thrusters (arcjets) with 2:1 hydrogen/nitrogen to simulate decomposed hydrazine as the propellant. The all-speed SIMPLE algorithm is employed to solve the governing equations, which take into account the effects of compressibility, the Lorentz force and Joule heating, as well as the temperature- and pressure-dependence of the gas properties. Typical computed results about the temperature, velocity and Mach number distributions within arcjet thruster are presented for the case with arc current of 9 A and inlet stagnant pressure of 3.3×105 Pa to show the flow and heat transfer characteristics. It is found that the propellant is heated mainly in the near-cathode and constrictor region, with the highest plasma temperature appearing near the cathode tip, and the flow transition from the subsonic to supersonic regime occurs within the constrictor region. The effect of gas viscosity on the plasma flow within arcjet thruster is examined by an additional numerical test using artificially reduced values of gas viscosity. The test results show that the gas viscosity appreciably affects the plasma flow and the performance of the arcjet thruster for the cases with the hydrazine or hydrogen as the propellant. The integrated axial Lorentz force in the thruster nozzle is also calculated and compared with the thrust force of the arcjet thruster. It is found that the integrated axial Lorentz force is much smaller than the thrust force for the low-power arcjet thruster. Modeling results for the NASA 1-kW class arcjet thruster with simulated hydrazine as the propellant are found to be reasonably consistent with available experimental data.

Measurement of neutron flow with the external target experimental device on CSR

With the construction of the neutron detection wall at the external target position on Heavy Ion Research Facility in Lanzhou-Cooling Storage Ring (HIRFL-CSR), it will be possible to detect high energy neutron. A BUU model is applied to simulate the flow in both symmetric (Ni+Ni, Pb+Pb) and asymmetric(Pb+Ni) systems. It is shown that at above several hundreds MeV/u, the flow signals are very obvious and depend clearly on the centrality of the collisions. Based on the products in the forward angle less than 20 degrees, the simulation also reveals that the determination of the reaction plane and the selection of the impact parameter, both of which are essential in the flow measurement, are well implemented. The double event and its influence on the determination of the neutron flow are also simulated.

Enhancements of the simulation method on the edge effect in resin transfer molding processes

In resin transfer molding processes, small clearances exist between the fiber preform and the mold edges, which result in a preferential resin flow in the edge channel and then disrupt the flow patterns during the mold filling stage. A mathematical model including the effect of cavity thickness on resin flow was developed for flow behavior involving the interface between an edge channel and a porous medium. According to the mathematical analysis of momentum equations in a fully developed rectangular duct and formulations of the equivalent edge permeability, comparing with three-dimensional Navier-Stokes equations, the governing equations were modified in the edge channel. The volume of fluid (VOF) method was applied to track the flow front. A simple case is numerically simulated using the modified governing equations. The effects of edge channel width and cavity thickness on flow front and inlet pressure are analyzed, and the evolution characteristics of simulated results are in agreement with the experimental results. (c) 2007 Elsevier B.V. All rights reserved

Numerical Simulation of Reactive Extrusion Processes for Activated Anionic Polymerization

In order to deal with the complicated relationships among the variables of the reactive extrusion process for activated anionic polymerization, a three-dimensional equivalent model of closely intermeshing co-rotating twin screw extruders was established. Then the numerical computation expressions of the monomer concentration, the monomer conversion, the average molecular weight and the fluid viscosity were deduced, and the numerical simulation of the reactive extrusion process of Styrene was carried out. At last, our simulated results were compared with Michaeli's simulated results and experimental results. (C) 2007 Elsevier B.V. All rights reserved

Computer simulation of reactive extrusion processes for free radical polymerization

The reactive extrusion for polymerization is an integrated polymer processing technology. A new semi-implicit iterative algorithm was proposed to deal with the complicated relationships among the chemical reaction, the macromolecular structure and the chemorheological property. Then the numerical computation expressions of the average molecular weight, the monomer conversion, and the initiator concentration were deduced, and the computer simulation of the reactive extrusion process for free radical polymerization was carried out, on basis of which reactive processing conditions can be optimized.

Preparation and simulation of an additional catalyst layer for direct methanol fuel cell

An additional anode catalyst layer with PtRu/C was hot pressed between two Nafion (R) 112 membranes and a conventional direct methanol fuel cell (DMFC) cathode/membrane/anode assembly with the above membranes as separator was fabricated. The additional catalyst layer formed an assistant cell with the cathode to prevent methanol crossover. A simple one-dimensional mathematical model was presented to describe the performance of this new type of membrane electrode assembly system. As seen from both experimental result and model analysis, the additional catalyst layer can not only effectively prevent the methanol crossover, but also generate electrical power with the crossover methanol. The percentage of output power of the assistant cell to the total power analyzed by the model is about 40% under usual condition, which is much higher than that from experimental result, indicating the potential of the development in the DMFC designing. It was also discovered that the electrical power generated from the assistant cell with crossover methanol could take higher percentage in total electrical power when the main DMFC current density became lower.

Mesoscopic simulation of the impregnating process of unidirectional fibrous preform in resin transfer molding

The resin transfer molding has gained popularity in the preparation of fiber-reinforced polymer-matrix composites because of its high efficiency and low pollution. The non-uniform inter-tow and intra-tow flows are regarded as the reason of void formation in RTM. According to the process characteristics, the axisymmetric model was developed to study the interaction between the flow in the inter-tow space and that in the intra-tow space. The flow behavior inside the fiber tows was formulated using Brinkman's equation, while that in the open space around the fiber tows was formulated by Stokes' equation. The volume of fluid (VOF) method was applied to track the flow front, and the effects of filling velocity, resin viscosity, inter-tow dimension and intra-tow permeability on fluid pressure and flow front were analyzed. The results show that the flow front difference between the inter-tow and intra-tow becomes larger with the decrease of intra-tow permeability, as well as the increase of filling velocity and inter-tow dimension.