Control of two-dimensional turbulent wall jet on a coanda surface

The Coanda effect that is the tendency of a fluid jet to stay attached to an adjacent curved surface that is very well shaped has been employed to improve the performance of various devices. The main objective of this paper is to investigate ways of keeping the flow attached to a larger length of a Coanda surface. There are considered two possibilities: one passive, which uses a slot that connects the low pressure and high pressure points on the Coanda surface and an active one, based on the principle of synthetic jet, created through an orifice located near the point of detachment of the jet. Reynolds averaged Navier-Stokes simulations (RANS) with shear stress transport k-ω (SST model) of Menter have been used to compute the two-dimensional turbulent wall jet. The numerical results are presented for the two methods considered. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A model for diffusion of mulucomponent information

The epidemic model of diffusion of news (or disease) is generalized to describe the diffusion of a multi‐component information. The multivaluedness of information in our model arises due to the large number (k) of constituent components or items of the information in question. When the different components of information are assumed to bear no hierarchy, the master equation of the model contains an intractably large number of variables (2 ^k ). The dynamics of the model, however, displays some simplifying features, one of which is the conservation of homogeneity of distribution of population over the different information vectors (in the sense defined in the text). The homogenized version of the model is found to be numerically tractable. The growth curves for large k continue to display sigmoid shapes, but with large ‘saturation times’. The dependence of ‘saturation time’ (i.e. the time required for spread of all the information in almost the entire population) on various parameters of the model, for uniform initial distributions, is numerically investigated. The ‘saturation time’ is found to vary inversely with the intensity of interaction (ß) and the size of population (N), as expected. An important numerical feature that emerges is that the ‘saturation time’ seems to be in linear proportion to the number of information items (k).     

On multiserver feedback retrial queue with finite buffer

This paper deals with a multiserver feedback retrial queueing system with finite waiting position and constant retrial rate. This system is analyzed as a quasi-birth-and-death (QBD) process and the necessary and sufficient condition for stability of the system is investigated. Some important system performance measures are obtained using matrix geometric method. The effect of various parameters on the system performance measures are illustrated numerically. Finally, the algorithmic development of the full busy period for the model under consideration is discussed.     

Assessment of turbulence modeling for gas flow in two-dimensional convergent–divergent rocket nozzle

In the present study, the turbulent gas flow dynamics in a two-dimensional convergent–divergent rocket nozzle is numerically predicted and the associated physical phenomena are investigated for various operating conditions. The nozzle is assumed to have impermeable and adiabatic walls with a flow straightener in the upstream side and is connected to a plenum surrounding the nozzle geometry and extended in the downstream direction. In this integrated component model, the inlet flow is assumed a two-dimensional, steady, compressible, turbulent and subsonic. The physics based mathematical model of the considered flow consists of conservation of mass, momentum and energy equations subject to appropriate boundary conditions as defined by the physical problem stated above. The system of the governing equations with turbulent effects is solved numerically using different turbulence models to demonstrate their numerical accuracy in predicting the characteristics of turbulent gas flow in such complex geometry. The performance of the different turbulence models adopted has been assessed by comparing the obtained results of the static wall pressure and the shock position with the available experimental and numerical data. The dimensionless shear stress at the nozzle wall and the separation point are also computed and the flow field is illustrated. The various implemented turbulence models have shown different behavior of the turbulent characteristics. However, the shear-stress transport (SST) k–ω model exhibits the best overall agreement with the experimental measurements. In general, the proposed numerical procedure applied in the present paper shows good capability in predicting the physical phenomena and the flow characteristics encountered in such kinds of complex turbulent flow.     

A computational model of a transcritical R744 ejector based on a homogeneous real fluid approach

A mathematical model of the compressible transonic single- and two-phase flow of a real fluid is discussed in this paper. The model was originally developed to simulate a refrigerant flow through a heat pump ejector. In the proposed approach, a temperature-based energy equation is replaced with an enthalpy-based formulation, in which the specific enthalpy, instead of the temperature, is an independent variable. A thermodynamic and mechanical equilibrium between gaseous and liquid phases is assumed for the two-phase flow. Consequently, real fluid properties, such as the density, the dynamic viscosity and the diffusion coefficient, are defined as functions of the pressure and the specific enthalpy. The energy equation formulation is implemented in commercial CFD software using subroutines that were developed in-house. The formulations was tested extensively for a single-phase flow of the R141b refrigerant, and for a two-phase flow of the R744 fluid (carbon dioxide) that occurred in a 3-D model of the ejector motive nozzle. In the model validation procedure, a satisfactory comparison between the experimental and computational results of the primary and secondary mass flow rates was obtained for both flow regimes. In addition, in the case of the R744 flow, the pressure distribution along the centre line of the ejector was accurately predicted as well. Furthermore, the results also shows that geometry modelling and measurement accuracy play an important in the final numerical results. As a result of the reasonable computational times, this method can be effectively used for the design of ejectors and also in geometric optimisation computations.     

Dynamic stability of glass-reinforced plastic shells

The problem of the dynamic stability of circular-cylindrical glass-reinforced plastic shells subjected to external transverse pressure is examined in the nonlinear formulation. After the Lagrange equations have been constructed, the problem reduces to the integration of a system of ordinary differential equations with aperiodic coefficients. The integration has been carried out numerically on a computer for various loading rates and shell parameters. Analogous problems for isotropic metal shells were examined in [1–4]. A review of the subject may be found in [5].Mekhanika Polimerov, Vol. 4, No. 1, pp. 109–115, 1968     

Parallel Strategies Applied to FE Simulations of Multiphasic Problems

While the theoretical background of various porous media models is well understood, it is still a demanding task to deal with these models numerically. In this contribution, a triphasic model is presented, which is capable of describing partially saturated soils. In quasi‐static conditions, this model results in the primary variables solid displacement, pore‐liquid pressure and pore‐gas pressure. For a stable numerical implementation, Taylor‐Hood elements are required, which need quadratic ansatz functions for the displacement and linear ansatz functions for the pressure terms. Looking at numerical simulations in 2‐d, challenging finite element calculations have already been realized in combination with adaptivity in time and space [1]. Nevertheless, new strategies have to be considered for a realization of applications of the model in 3‐d in order to handle the huge amount of unknowns arising from the discretization with Taylor‐Hood elements. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Second derivative estimation using harmonic analysis

Simulation sensitivity analysis is an important problem for simulation practitioners analyzing complex systems. The significance of this problem has resulted in the development of various gradient estimators that can be used to address this issue. Although higher derivative estimators have been discussed concurrently, less attention has been given to assess the efficiency and feasibility of computing such estimators. In this paper, two second derivative estimators are presented. The first estimators, called the HFD estimators, combine harmonic gradient estimators with finite differences second derivative estimators. The resulting hybrid estimators requireO(p) fewer simulation runs to implement compared to the straightforward finite differences approach, wherep is the number of input parameters in the simulation model. The second estimators, called the HA estimators, incorporate harmonic analysis directly, requiring one or two simulation runs to implement, depending on whether a control variate simulation run is made. Expressions for the bias and the variance of the HFD and the HA estimators (with and without variance reduction techniques) are derived. Optimal mean squared error convergence rates are also discussed. In particular, the convergence rates for both these estimators are shown to be the same, though the computational performance of the HFD estimators is better than that for the HA estimators on anM/M/1 queue simulation model. Computational results for the HFD estimators on an (s, S) inventory system simulation model are also included.     

Peripheral layer viscosity effects on periodic blood flow through rigid tube of thin diameter

A mathematical model has been presented for periodic blood flow in a rigid circular tube of thin diameter. Blood is presented as a 3-layered fluid by considering core fluid as a casson fluid which is covered by a thin layer of Newtonian fluid (plasma). The energy integral method has been used to obtain the unsteady pressure gradients as suggested by Elkouh [2]. The results obtained for velocity profiles have been compared with the experimental results of Bugliarello and Sevilla (Biorheology 7 (1970), 85). The effects of various parameters on wall shearing stress has also been brought out and discussed.     

Simulation of Multiphase Flows with Strong Shocks and Density Variations

The system of extended Euler type hyperbolic equations is considered to describe a two-phase compressible flow. A numerical scheme for computing multi-component flows is then examined. The numerical approach is based on the mathematical model that considers interfaces between fluids as numerically diffused zones. The hyperbolic problem is tackled using a high resolution HLLC scheme on a fixed Eulerian mesh. The global set of conservative equations (mass, momentum and energy) for each phase is closed with a general two parameters equation of state for each constituent. The performance of various variants of a diffuse interface method is carefully verified against a comprehensive suite of numerical benchmark test cases in one and two space dimensions. The studied benchmark cases are divided into two categories: idealized tests for which exact solutions can be generated and tests for which the equivalent numerical results could be obtained using different approaches. The ability to simulate the Richtmyer-Meshkov instabilities, which are generated when a shock wave impacts an interface between two different fluids, is considered as a major challenge for the present numerical techniques. The study presents the effect of density ratio of constituent fluids on the resolution of an interface and the ability to simulate Richtmyer-Meshkov instabilities by various variants of diffuse interface methods. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Finite-difference analysis of natural convection flow of a viscous fluid in a porous channel with constant heat sources

The fully developed natural convection flow of a viscous fluid in a porous channel is modeled and studied numerically. The walls are kept at constant temperatures. The effects of various dimensionless parameters emerging in the model are studied graphically. It has been noted that the velocity and temperature both depend on the heat source and the free convection parameters.     

Numerical and analytical simulation of peristaltic flow of a Jeffrey-six constant fluid

This paper describes the Peristaltic flow of a Jeffrey-six constant fluid in an endoscope. The two-dimensional equation of Jeffrey-six constant fluid is simplified by making the assumptions of long wave length and low Reynolds number. The reduced momentum equations are solved with three methods, namely (i) Perturbation method, (ii) Homotopy analysis method, and (iii) shooting method. The comparison of the three solutions shows a very good agreement between the three results. The expressions for pressure rise and frictional forces per wave length have been also computed numerically. Finally, the pressure rise, frictional forces are plotted for different parameters of interest.     

Numerical investigation of pressure loss reduction in a power plant stack

The strengthened environmental laws require the power plants to reduce the emissions. Flue gas desulphurization and deNO_x involve adding chemicals to the flow stream, thereby resulting in increased mass flow. This problem could be overcome by reducing the pressure drop in the duct work and stack combination, so that a higher flow at reduced pressure drop can be handled by the existing fans. In this study, a power plant stack model of 1:40 was investigated numerically. The pressure reduction was achieved by introduction of baffles with various orientations and turning vanes at the inlet of the stack. The flows were modeled and analyzed using commercial computational fluid dynamics (CFD) software Fluent 6.2. The numerical results were validated with the experimental data. The 30° baffle without turning vanes was found to be the optimum baffle angle in terms of the pressure loss reduction. Variation of axial velocity, swirling component and turbulence kinetic energy along the axis of the stack was analyzed to understand the mechanism of the pressure loss reduction in a power plant stack. Guidelines for further pressure loss reduction were provided based on the insight gained from the simulation results.     

Computation of melting with natural convection inside a rectangular enclosure heated by discrete protruding heat sources

This paper presents the results of a numerical investigation of the heat transfer by natural convection during the melting of a phase change material (PCM, n-eicosane with melting point of 36 °C) contained in a rectangular enclosure. This latest is heated by three discrete protruding heat sources (simulating electronic components) placed on one of its vertical walls. The power generated by heat sources is dissipated in PCM. The advantage of using this cooling scheme is that the PCMs are able to absorb high amount of heat generated by the heat sources, without acting the fan during the charging process (melting of the PCM). The thermal behavior and thermal performance of the proposed PCM based-heat sink are numerically investigated by developing a mathematical model based on the mass, momentum and energy conservation equations. The obtained numerical results show the impact of various key parameters on the cooling capacity of the PCM-based heat sink. Correlations encompassing a wide range of parameters were developed in terms of the dimensionless secured operating time (time required by one of the electronic components before reaching its critical temperature, T_cr ∼ 75 °C) and the corresponding liquid fraction, using the asymptotic computational fluid dynamics (ACFD) technique.     

A novel constitutive law for Silicon under contact loading

Silicon (Si) remains the most important semiconductor material to date. The understanding of its deformation behavior under contact (indenter-) loading is crucial to improving technologically relevant abrasive machining techniques (lapping, sawing, grinding). While it has been long established that Si undergoes a series of stress driven phase transitions upon compression and subsequent pressure release, to the authors' knowledge, no material model is available that adequately captures this behavior. In particular, reverse transformation in unloading has received too little attention. A novel phenomenological, thermomechanical model based on experimental observations and MD predictions is presented in this work. It captures both the cd-Si → β-Si transition upon compression and the β-Si → a-Si transition upon rapid decompression, which are most relevant for indenter loading. To control inelasticity in unloading, the dissipation function was augmented by a kinematic constraint on the tensorial internal variable. In stress space, the transformation surfaces are hyperboloids of revolution aligned along the hydrostatic axis. The non-linear model was numerically implemented in a finite element code using an iterative implicit algorithm and successfully applied to simple loading cases. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Mixed-finite element and finite volume discretization for heavy brine simulations in groundwater

Recently, a new theory of high-concentration brine transport in groundwater has been developed. This approach is based on two nonlinear mass conservation equations, one for the fluid (flow equation) and one for the salt (transport equation), both having nonlinear diffusion terms. In this paper, we present and analyze a numerical technique for the solution of such a model. The approach is based on the mixed hybrid finite element method for the discretization of the diffusion terms in both the flow and transport equations, and a high-resolution TVD finite volume scheme for the convective term. This latter technique is coupled to the discretized diffusive flux by means of a time-splitting approach. A commonly used benchmark test (Elder problem) is used to verify the robustness and nonoscillatory behavior of the proposed scheme and to test the validity of two different formulations, one based on using pressure head ψ and concentration c as dependent variables, and one using pressure p and mass fraction ω as dependent variables. It is found that the latter formulation gives more accurate and reliable results, in particular, at large times. The numerical model is then compared against a semi-analytical solution and the results of a laboratory test. These tests are used to verify numerically the performance and robustness of the proposed numerical scheme when high-concentration gradients (i.e., the double nonlinearity) are present.     

A numerical model of the reacting multiphase flow in a pulp digester

This paper summarises the development, implementation and typical results from a new full-scale three-dimensional numerical model of the multiphase chemically reacting flow in a continuous digester. The model is based on continuity and momentum equations for wood chips and free liquid. A new sub-model describing the tangential stresses for the solid phase has been developed and incorporated into the digester model. The solution procedure that takes into account the high inter-phase friction terms and the high values of the solid pressure has been introduced. In order to illustrate the model’s performance, simulations have been carried on for the industrial digester. The results are in quantitative agreement with available field measurements and in qualitative agreement with operating personnel observations. The model utilises curvilinear body fitted coordinates and can be used to simulate the operation of various pulp digesters or any other chemical reactors working in a similar regime.     

Predictive Inference for the Elliptical Linear Model

This paper derives the prediction distribution of future responses from the linear model with errors having an elliptical distribution with known covariance parameters. For unknown covariance parameters, the marginal likelihood function of the parameters has been obtained and the prediction distribution has been modified by replacing the covariance parameters by their estimates obtained from the marginal likelihood function. It is observed that the prediction distribution with elliptical error has a multivariate Student'st-distribution with appropriate degrees of freedom. The results for some special cases such as the Intra-class correlation model, AR(1), MA(1), and ARMA(1,1) models have been obtained from the general results. As an application, theβ-expectation tolerance region has been constructed. An example has been added.     

New approach for the identification of power system economic dispatch parameters

It is well known that the coefficients of the input-output characteristics of the thermal steam-turbine model as well as the network model parameters have a great effect on the optimal economic operation of all thermal-electric power systems. Until today, these coefficients, the loss formula coefficients, theB-coefficients, and the active-reactive power loss models have been estimated using the well-known least-square estimation algorithm.In this paper, we present a new algorithm to estimate the power system parameters for economic dispatch calculation (EDC); this algorithm is based on the least absolute-value approximation (LAV)l ₁-norm. We compare the results obtained using the proposed algorithms with those obtained using the least-square error algorithm (LS). Optimal costs as well as overall network performance resulting from the implementation of each technique provide the basis of our conclusion.This work was supported by the Natural Science and Engineering Research Council of Canada, Grant No. A4146.