Turbulent particle dispersion in an electrostatic precipitator

The behaviour of charged particles in turbulent gas flow in electrostatic precipitators (ESPs) is crucial information to optimise precipitator efficiency. This paper describes a strongly coupled calculation procedure for the rigorous computation of particle dynamics during ESP taking into account the statistical particle size distribution. The turbulent gas flow and the particle motion under electrostatic forces are calculated by using the commercial computational fluid dynamics (CFD) package FLUENT linked to a finite volume solver for the electric field and ion charge. Particle charge is determined from both local electrical conditions and the cell residence time which the particle has experienced through its path. Particle charge density and the particle velocity are averaged in a control volume to use Lagrangian information of the particle motion in calculating the gas and electric fields. The turbulent particulate transport and the effects of particulate space charge on the electrical current flow are investigated. The calculated results for poly-dispersed particles are compared with those for mono-dispersed particles, and significant differences are demonstrated.     

Turbulent diffusion of mass in circular pipe flow

An implicit finite difference scheme was used to solve the convective-diffusion equation to predict the steady-state transport of a conservative, neutrally bouyant tracer injected along the centreline into a fully developed turbulent pipe flow. Three different distributions for the radial mass diffusivity have been compared with two independent sets of experimental data. The results indicate that the distribution based on the turbulent kinematic eddy viscosity predicted by a k?l model produces the closest agreement between the numerical model predictions and the experimentally observed tracer distribution.     

An assessment of turbulence models applied to the simulation of a two-dimensional submerged jet

This paper is concerned with the investigation of the performance of different turbulence models in the numerical prediction of transient flow caused by a confined submerged jet. Several widely used models, i.e., the standard k–ε, RNG k–ε, low Reynolds number k–ε models and the differential Reynolds stress model, as included in CFD codes, were compared with each other for a two-dimensional, incompressible, turbulent jet flow and with reported experimental data. A flapping oscillation was predicted regardless of the model used. A chosen Strouhal (St) number definition brought the fundamental frequencies from both the experiments and computations into close proximity. However, different turbulence models have exhibited quite different behaviours in terms of the frequency and regularity of the oscillation and in terms of the scale and duration of the vortices generated. All versions of the k–ε model yielded regular oscillations, which agree with experimental observations. On the other hand, the Reynolds stress (RS) model produced a complex pattern but a slower dissipation of vortices. In addition, some aspects of gridding and inflow representation are also discussed.     

Mathematical modeling of coupled turbulent flow and solidification in a single belt caster with electromagnetic brake

A mathematical model has been developed to simulate turbulent fluid flow and solidification in the presence of a DC magnetic field in an extended nozzle for metal delivery to a single belt caster. This paper reports on predicted effects of DC magnetic field conditions in modifying flows and solidification behavior in the metal delivery system. It is shown that the application of a DC magnetic brake to the proposed system can result in a reasonably uniform feeding of melt onto the cooled moving belt. This, in turn, optimises the rate of even shell growth along the chilled substrate. In order to account for the effects of turbulence, a revised low-Reynolds k–ε turbulent model was employed. A Darcy-porosity approach was used to simulate fluid flow within the mushy solidification region. Simulations were carried out for plain carbon steel strip casting. The fully coupled transport equations were numerically solved using the finite volume method. The computed flow patterns were compared with those reported in the literature. The performance of the magnetic flow control device proposed in this work is evaluated and compared with flow modifications obtained by inserting a ceramic filter within the reservoir.     

Multi-objective analysis of dam release flows in rivers downstream from hydropower reservoirs

This paper assesses dam releases from hydropower reservoirs in order to optimize power production and fish habitat protection. A multi-objective programming model includes output from 2-D hydraulic simulation for habitat assessment to optimize power production and fish habitat suitability as a Pareto set. To identify the optimal Pareto set two different approaches are used and compared: ε-constraint methods and non-dominant-sorting genetic algorithm (NSGA II). To formulate the ecological objective the river habitat quality is quantified by the weighted usable area (WUA). The relation between the WUA and the river flow-rate is obtained by using a 2D hydraulic model in which the hydraulic characteristics of river current – depth and velocity – are calculated by a finite difference numerical integration of two-dimensional shallow water equations on a boundary fitted non orthogonal curvilinear mesh. This approach allows the integration of motion equations on geometrically complex domains as those representing the morphology of natural watercourses. The performance of the proposed methodology is analyzed in a case study of a stretch of the Piave river downstream of the dam of the Pieve di Cadore reservoir (Belluno, Italy).     

On the structure of some group codes

In this paper, we study the following problem: Which characteristics does a codeC possess when the syntactic monoidsyn(C ^*) of the star closureC ^* ofC is a group? For a codeC, if the syntactic monoidsyn(C ^*) is a group, then we callC a group code. This definition of a group code is different from the one in [1] (see [1], 46–47). Schützenberger had characterized the structure of finite group codes and had proved thatC is a finite group code if and only ifC is a full uniform code (see [5], [8]). Fork-prefix andk-suffix codes withk≥2,k-infix,k-outfix,p-infix,s-infix, right semaphore codes and left semaphore codes, etc., we obtain similar results. It is proved that the above mentioned codes are group codes if and only if they are uniform codes.     

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.     

Direct numerical simulation of turbulent non-Newtonian flow using OpenFOAM

Understanding transition and turbulence in the flow of shear-thinning non-Newtonian fluids remains substantially unresolved and additional research is required to develop better computational methods for wall-bounded turbulent flows of these fluids. Previous DNS studies of shear-thinning fluids mainly use purpose-built codes and simple geometries such as pipes and channels. However in practical application, the geometry of mixing vessels, pumps and other process equipment is far more complex, and more flexible computational methods are required. In this paper a general-purpose DNS approach for shear-thinning fluids is undertaken using the OpenFOAM CFD library. DNS of turbulent Newtonian and non-Newtonian flow in a pipe flow are conducted and the accuracy and efficiency of OpenFOAM are assessed against a validated high-order spectral element-Fourier DNS code – Semtex. The results show that OpenFOAM predicts the flow of shear-thinning fluids to be a little more transitional than the predictions from Semtex, with lower radial and azimuthal turbulence intensities and higher axial intensity. Despite this, the first and second order turbulence statistics differ by at most 16%, and usually much less. An assessment of the parallel scaling of OpenFOAM indicates that OpenFOAM scales very well for the CPUs from 8 to 512, but the intranode scalability is poor for less than 8CPUs. The present work shows that OpenFOAM can be used for DNS of shear-thinning fluids in the simple case of pipe flow, and suggests that more complex flows, where flow separation is often important, are likely to be simulated with accuracies that are acceptably good for engineering application.     

Pipe roughness identification of water distribution networks: The full turbulent case

This paper proposes a technique to identify individual pipe roughness parameters in a water distribution network by means of the inversion of the steady-state hydraulic network equations. By enabling the reconstruction of these hydraulic friction parameters to be reliable, this technique improves the conventional model’s accuracy and thereby promises to enhance model-based leakage detection and localization. As it is the case in so-called fireflow tests, this methodology is founded on the premise to measure the pressure distributed at a subset of nodes in the network’s graph while assuming the nodal consumption to be known. Beside of the proposed problem formulation, which is restricted to only allow turbulent flow in each of the network’s pipes initially, developed algorithms are presented and evaluated using simulation examples.     

A mass formula for self-dual permutation codes

For a module V over a finite semisimple algebra A we give the total number of self-dual codes in V. This enables us to obtain a mass formula for self-dual codes in permutation representations of finite groups over finite fields of coprime characteristic.     

Constacyclic codes over finite commutative semi-simple rings

Finite commutative semi-simple rings are direct sum of finite fields. In this study, we investigate the algebraic structure of λ-constacyclic codes over such finite semi-simple rings. Among others, necessary and sufficient conditions for the existence of self-dual, LCD, and Hermitian dual-containing λ-constacyclic codes over finite semi-simple rings are provided. Using the CSS and Hermitian constructions, quantum MDS codes over finite semi-simple rings are constructed.     

CFD simulation of precession in sudden pipe expansion flows with low inlet swirl

In this paper, three-dimensional, time-dependent calculations are carried out using the finite volume CFD code CFX4 and the VLES approach with standard k–ε model to simulate the turbulent swirl flow in an axisymmetric sudden expansion with an expansion ratio of 5.0 for a Reynolds number of 10⁵. This flow is unstable over the entire swirl number range considered between 0 and 0.48, and a large-scale coherent structure is found to precess about the centerline. Compared with the unswirled case, inclusion of a slight inlet swirl (swirl number below 0.23) can reduce the precession speed, cause the precession to be against the mean swirl and suppress the flapping motion. Several modes of precession are predicted as the swirl intensity increases, in which the precession, as well as the spiral structure, reverses direction. Accompanying the transition between different modes, abrupt changes in precession frequency are also experienced. Grid sensitivity and comparison with smaller expansion ratio data are also discussed.     

CFD of turbulent transport of particles behind a backward-facing step using a new model—k-ε-Sp

The k-ε-S_p model, describing two-dimensional gas–solid two-phase turbulent flow, has been developed. In this model, the diffusion flux and slip velocity of solid particles are introduced to represent the particle motion in two-phase flow. Based on this model, the gas–solid two-phase turbulent flow behind a vertical backward-facing step is simulated numerically and the turbulent transport velocities of solid particles with high density behind the step are predicted. The numerical simulation is validated by comparing the results of the numerical calculation with two other two-phase turbulent flow models (k-ε-A_p, k-ε-k_p) by Laslandes and the experimental measurements. This model, not only has the same virtues of predicting the longitudinal transport of the solid particles as the present practical two-phase flow models, but also can predict the lateral transport of the solid particles correctly.     

Non-invertible-element constacyclic codes over finite PIRs

In this paper we introduce the notion of λ-constacyclic codes over finite rings R for arbitrary element λ of R. We study the non-invertible-element constacyclic codes (NIE-constacyclic codes) over finite principal ideal rings (PIRs). We determine the algebraic structures of all NIE-constacyclic codes over finite chain rings, give the unique form of the sets of the defining polynomials and obtain their minimum Hamming distances. A general form of the duals of NIE-constacyclic codes over finite chain rings is also provided. In particular, we give a necessary and sufficient condition for the dual of an NIE-constacyclic code to be an NIE-constacyclic code. Using the Chinese Remainder Theorem, we study the NIE-constacyclic codes over finite PIRs. Furthermore, we construct some optimal NIE-constacyclic codes over finite PIRs in the sense that they achieve the maximum possible minimum Hamming distances for some given lengths and cardinalities.     

Modelling gas injection of a Peirce–Smith-converter

A commercial CFD-code PHOENICS was used to solve isothermal flow field of gas and liquid in a Peirce–Smith-converter. An Euler–Euler based algorithm was chosen for modelling fluid dynamics and evaluating controlling forces of a submerged gas injection. Predictions were made with a k–ε turbulence model in the body fitted coordinate system. The model has been verified with a 1/4 scale water model, and a parametric study with the mathematical model of submerged gas injection was made for the PS-process and the ladle injection processes. Limits of the modelling technique used were recognised, but calculated results indicate that the present model predicts the general flow field with reasonable accuracy. Predicted bubble distribution, pattern of the flow field and magnitude of flow velocities were used to evaluate scaling factors of physical models and general flow conditions of an industrial PS-converter.     

Linear closures of finite geometries

The interrelations between finite geometries (finite incidence structures) and linear codes over finite fields are discussed under some special fundamental aspects. For any incidence structure \({\mathcal{I}}\) block codes, block-difference codes and co-block codes over finite fields of characteristic p are discussed resp. introduced; correspondingly p-modular co-blocks are defined for \({\mathcal{I}}\). Orthogonality modulo p is introduced as a concept relating different geometries having the same point set. Conversely three types of block-tactical geometries may be derived from vector classes of fixed Hamming weight in a given linear code. These geometries are tactical configurations if the given code admits a transitive permutation group. A combination of both approaches leads to the concept of p-closure of a finite geometry and to the notions of p-closed, weakly p-closed and p-dense incidence structures. These geometric concepts are applied to simple or directed graphs via their natural “adjacency geometry”. Here the above mentioned code theoretic treatment leads to the concept of p-modular co-adjacent vertex sets. As instructive examples the Petersen graph, its complemetary graph and the Higman-Sims graph are considered.     

Constructive spherical codes near the Shannon bound

Shannon gave a lower bound in 1959 on the binary rate of spherical codes of given minimum Euclidean distance ρ. Using nonconstructive codes over a finite alphabet, we give a lower bound that is weaker but very close for small values of ρ. The construction is based on the Yaglom map combined with some finite sphere packings obtained from nonconstructive codes for the Euclidean metric. Concatenating geometric codes meeting the TVZ bound with a Lee metric BCH code over GF(p), we obtain spherical codes that are polynomial time constructible. Their parameters outperform those obtained by Lachaud and Stern (IEEE Trans Inf Theory 40(4):1140–1146, 1994). At very high rate they are above 98% of the Shannon bound.     

Self-Orthogonal and Self-Dual Codes Constructed via Combinatorial Designs and Diophantine Equations

Combinatorial designs have been widely used, in the construction of self-dual codes. Recently, new methods of constructing self-dual codes are established using orthogonal designs (ODs), generalized orthogonal designs (GODs), a set of four sequences and Diophantine equations over GF(p). These methods had led to the construction of many new self-dual codes over small finite fields and rings. In this paper, we used some methods to construct self-orthogonal and self dual codes over GF(p), for some primes p. The construction is achieved by using some special kinds of combinatorial designs like orthogonal designs and GODs. Moreover, we combine eight circulant matrices, a system of Diophantine equations over GF(p), and a recently discovered array to obtain a new construction method. Using this method new self-dual and self-orthogonal codes are obtained. Specifically, we obtain new self-dual codes [32,16,12] over GF(11) and GF(13) which improve the previously known distances.     

A 3D non-linear k–ε turbulent model for prediction of flow and mass transport in channel with vegetation

The results from a 3D non-linear k–ε turbulence model with vegetation are presented to investigate the flow structure, the velocity distribution and mass transport process in a straight compound open channel and a curved open channel. The 3D numerical model for calculating flow is set up in non-orthogonal curvilinear coordinates in order to calculate the complex boundary channel. The finite volume method is used to disperse the governing equations and the SIMPLEC algorithm is applied to acquire the coupling of velocity and pressure. The non-linear k–ε turbulent model has good useful value because of taking into account the anisotropy and not increasing the computational time. The water level of this model is determined from 2D Poisson equation derived from 2D depth-averaged momentum equations. For concentration simulation, an expression for dispersion through vegetation is derived in the present work for the mixing due to flow over vegetation. The simulated results are in good agreement with available experimental data, which indicates that the developed 3D model can predict the flow structure and mass transport in the open channel with vegetation.