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.

Otimização da co-pirólise de gasóleo pesado com embalagens de PEAD pós-uso utilizando planejamento fatorial

A co-pirólise é uma rota promissora, uma vez que minimiza o impacto ambiental causado pela disposição do plástico de maneira inadequada, evita seu acúmulo em lixões e permite um melhor aproveitamento de um recurso natural não-renovável, o petróleo, matéria prima importante para a geração de energia e obtenção de produtos químicos. O presente trabalho teve como objetivo a definição das condições experimentais mais propícias à obtenção de líquidos pirolíticos com alta fração de óleo diesel, resultantes da co-pirólise de polietileno de alta densidade (PEAD) pós-consumo com gasóleo pesado tilizando-se catalisador de FCC (Fluid Catalytic Cracking). Como instrumento de otimização das condições experimentais, optou-se pela Metodologia Planejamento Fatorial. Foi também estudado o efeito das condições experimentais, como: a temperatura de reação, a relação gasóleo/polietileno e a quantidade de catalisador no meio reacional. As amostras de polietileno, gasóleo e catalisador foram submetidas à co-pirólise em sistema de leito fixo, sob fluxo constante de nitrogênio, variando-se a temperatura entre 450 C a 550 C, a quantidade de PEAD no meio reacional foi de 0,2 a 0,6 g, e a quantidade de catalisador foi de zero a 0,06 g, mantendo-se fixa a quantidade de gasóleo em 2 g. Foram efetuadas as caracterizações física e química do gasóleo, polietileno pós-uso e do catalisador. Como resultado, obteve-se a produção de 87% de fração de óleo diesel em duas condições diferentes: (a) 550 0C de temperatura sem catalisador; (b) 500 0C de temperatura e 25% de catalisador FCC. Em ambos os casos, a quantidade de gasóleo pesado e PEAD foram constantes (2 g Gasóleo; 0,2 g PEAD), assim com o tempo de reação de 15 minutos. A fração de óleo diesel obtida neste estudo alcançou o poder calorífico de 44,0 MJ/Kg que é similar ao óleo diesel comercial. É importante ressaltar que em ambos os casos nenhum resíduo foi produzido, sendo uma rota ambientalmente importante para reciclagem de embalagens plásticas contaminadas com óleo lubrificante originárias de postos de serviço, visando à recuperação de ambos conteúdo energético e orgânico dos resíduos de embalagens plásticas pós-uso

Characteristics of blood flow in microchannels and relevant impact on modelling blood behaviour in biochip separators

This paper reports a perspective investigation of computational modelling of blood fluid in microchannel devices as a preparation for future research on fluid-structure interaction (FSI) in biofluid mechanics. The investigation is carried out through two aspects, respectively on physical behaviours of blood flow in microchannels and appropriate methodology for modelling. The physics of blood flow is targeted to the challenges for describing blood flow in microchannels, including rheology of blood fluid, suspension features of red blood cells (RBCs), laminar hydrodynamic influence and effect of surface roughness. The analysis shows that due to the hyperelastic property of RBC and its comparable dimension with microchannels, blood fluid shows complex behaviours of two phase flow. The trajectory and migration of RBCs require accurate description of RBC deformation and interaction with plasma. Following on a discussion of modelling approaches, i.e. Eulerian method and Lagrangian method, the main stream modelling methods for multiphase flow are reviewed and their suitability to blood flow is analysed. It is concluded that the key issue for blood flow modelling is how to describe the suspended blood cells, modelled by Lagrangian method, and couple them with the based flow, modelled by Eulerian method. The multiphase flow methods are thereby classified based on the number of points required for describing a particle, as follows: (i) single-point particle methods, (ii) mutli-point particle methods, (iii) functional particle methods, and (iv) fluid particle methods. While single-point particle methods concentrate on particle dynamic movement, multipoint and functional particle methods can take into account particle mechanics and thus offer more detailed information for individual particles. Fluid particle methods provide good compromise between two phases, but require additional information for particle mechanics. For furthermore detailed description, we suggest to investigate the possibility using two domain coupling method, in which particles and base flow are modelled by two separated solvers. It is expected that this paper could clarify relevant issues in numerical modelling of blood flow in microchannels and induce some considerations for modelling blood flow using multiphase flow methods. © 2012 IEEE.

The impact of sand slugs against beams and plates: Coupled discrete particle/finite element simulations

The impact of a slug of dry sand particles against a metallic sandwich beam or circular sandwich plate is analysed in order to aid the design of sandwich panels for shock mitigation. The sand particles interact via a combined linear-spring-and-dashpot law whereas the face sheets and compressible core of the sandwich beam and plate are treated as rate-sensitive, elastic-plastic solids. The majority of the calculations are performed in two dimensions and entail the transverse impact of end-clamped monolithic and sandwich beams, with plane strain conditions imposed. The sand slug is of rectangular shape and comprises a random loose packing of identical, circular cylindrical particles. These calculations reveal that loading due to the sand is primarily inertial in nature with negligible fluid-structure interaction: the momentum transmitted to the beam is approximately equal to that of the incoming sand slug. For a slug of given incoming momentum, the dynamic deflection of the beam increases with decreasing duration of sand-loading until the impulsive limit is attained. Sandwich beams with thick, strong cores significantly outperform monolithic beams of equal areal mass. This performance enhancement is traced to the "sandwich effect" whereby the sandwich beams have a higher bending strength than that of the monolithic beams. Three-dimensional (3D) calculations are also performed such that the sand slug has the shape of a circular cylindrical column of finite height, and contains spherical sand particles. The 3D slug impacts a circular monolithic plate or sandwich plate and we show that sandwich plates with thick strong cores again outperform monolithic plates of equal areal mass. Finally, we demonstrate that impact by sand particles is equivalent to impact by a crushable foam projectile. The calculations on the equivalent projectile are significantly less intensive computationally, yet give predictions to within 5% of the full discrete particle calculations for the monolithic and sandwich beams and plates. These foam projectile calculations suggest that metallic foam projectiles can be used to simulate the loading by sand particles within a laboratory setting. © 2013 Elsevier Ltd.

Wetness loss prediction for a low pressure steam turbine using computational fluid dynamics

Two-phase computational fluid dynamics modelling is used to investigate the magnitude of different contributions to the wet steam losses in a three-stage model low pressure steam turbine. The thermodynamic losses (due to irreversible heat transfer across a finite temperature difference) and the kinematic relaxation losses (due to the frictional drag of the drops) are evaluated directly from the computational fluid dynamics simulation using a concept based on entropy production rates. The braking losses (due to the impact of large drops on the rotor) are investigated by a separate numerical prediction. The simulations show that in the present case, the dominant effect is the thermodynamic loss that accounts for over 90% of the wetness losses and that both the thermodynamic and the kinematic relaxation losses depend on the droplet diameter. The numerical results are brought into context with the well-known Baumann correlation, and a comparison with available measurement data in the literature is given. The ability of the numerical approach to predict the main wetness losses is confirmed, which permits the use of computational fluid dynamics for further studies on wetness loss correlations. © IMechE 2013 Reprints and permissions: sagepub.co.uk/journalsPermissions.nav.

Mantle-derived gaseous components in ore-forming fluids of the Xiangshan uranium

The Xiangshan U deposit, the largest hydrothermal U deposit in China, is hosted in late Jurassic felsic volcanic rocks although the U mineralization post dates the volcanics by at least 20 Ma. The mineralization coincides with intrusion of local mantle-derived mafic dykes formed during Cretaceous crustal extension in South China. Ore-forming fluids are rich in CO2, and U in the fluid is thought to have been dissolved in the form of UO2 (CO3)22− and UO2 (CO3) 34− complexes. This paper provides He and Ar isotope data of fluid inclusions in pyrites and C isotope data of calcites associated with U mineralization (pitchblende) in the Xiangshan U deposit. He isotopic compositions range between 0.1 and 2.0Ra (where Ra is the 3He/4He ratio of air=1.39×10−6) and correlates with 40Ar/36Ar; although there is potential for significant 3He production via 6Li(n,α)3H(β)3He reactions in a U deposit (due to abundant neutrons), nucleogenic production cannot account for either the 3He concentration in these fluids, nor the correlations between He and Ar isotopic compositions. It is more likely that the high 3He/4He ratios represent trapped mantle-derived gases. A mantle origin for the volatiles of Xiangshan is consistent with the δ13C values of calcites, which vary from −3.5‰ to −7.7‰, overlapping the range of mantle CO2. The He, Ar and CO2 characteristics of the ore-forming fluids responsible for the deposit are consistent with mixing between 3He- and CO2-rich mantle-derived fluids and CO2-poor meteoric fluids. These fluids were likely produced during Cretaceous extension and dyke intrusion which permitted mantle-derived CO2 to migrate upward and remobilize U from the acid volcanic source rocks, resulting in the formation of the U deposit. Subsequent decay of U within the fluid inclusions has reduced the 3He/4He ratio, and variations in U/3He result in the range in 3He/4He observed with U/3He ratios in the range 5–17×103 likely corresponding to U concentrations in the fluids b0.2 ppm.

A volume of fluid numerical model for wave impacts at coastal structures

The first stages in the development of a new design tool, to be used by coastal engineers to improve the efficiency, analysis, design, management and operation of a wide range of coastal and harbour structures, are described. The tool is based on a two-dimensional numerical model, NEWMOTICS-2D, using the volume of fluid (VOF) method, which permits the rapid calculation of wave hydrodynamics at impermeable natural and man-made structures. The critical hydrodynamic flow processes and forces are identified together with the equations that describe these key processes. The different possible numerical approaches for the solution of these equations, and the types of numerical models currently available, are examined and assessed. Preliminary tests of the model, using comparisons with results from a series of hydraulic model test cases, are described. The results of these tests demonstrate that the VOF approach is particularly appropriate for the simulation of the dynamics of waves at coastal structures because of its flexibility in representing the complex free surfaces encountered during wave impact and breaking. The further programme of work, required to develop the existing model into a tool for use in routine engineering design, is outlined.

Turbulence modelling and its impact on CFD predictions for cooling of electronic components

This paper will discuss Computational Fluid Dynamics (CFD) results from an investigation into the accuracy of several turbulence models to predict air cooling for electronic packages and systems. Also new transitional turbulence models will be proposed with emphasis on hybrid techniques that use the k-ε model at an appropriate distance away from the wall and suitable models, with wall functions, near wall regions. A major proportion of heat emitted from electronic packages can be extracted by air cooling. This flow of air throughout an electronic system and the heat extracted is highly dependent on the nature of turbulence present in the flow. The use of CFD for such investigations is fast becoming a powerful and almost essential tool for the design, development and optimization of engineering applications. However turbulence models remain a key issue when tackling such flow phenomena. The reliability of CFD analysis depends heavily on the turbulence model employed together with the wall functions implemented. In order to resolve the abrupt fluctuations experienced by the turbulent energy and other parameters located at near wall regions and shear layers a particularly fine computational mesh is necessary which inevitably increases the computer storage and run-time requirements. The PHYSICA Finite Volume code was used for this investigation. With the exception of the k-ε and k-ω models which are available as standard within PHYSICA, all other turbulence models mentioned were implemented via the source code by the authors. The LVEL, LVEL CAP, Wolfshtein, k-ε, k-ω, SST and kε/kl models are described and compared with experimental data.

Investigation into the performance of turbulence models for fluid flow and heat transfer phenomena in electronic applications

Heat is extracted away from an electronic package by convection, conduction, and/or radiation. The amount of heat extracted by forced convection using air is highly dependent on the characteristics of the airflow around the package which includes its velocity and direction. Turbulence in the air is also important and is required to be modeled accurately in thermal design codes that use computational fluid dynamics (CFD). During air cooling the flow can be classified as laminar, transitional, or turbulent. In electronics systems, the flow around the packages is usually in the transition region, which lies between laminar and turbulent flow. This requires a low-Reynolds number numerical model to fully capture the impact of turbulence on the fluid flow calculations. This paper provides comparisons between a number of turbulence models with experimental data. These models included the distance from the nearest wall and the local velocity (LVEL), Wolfshtein, Norris and Reynolds, k-ε, k-ω, shear-stress transport (SST), and kε/kl models. Results show that in terms of the fluid flow calculations most of the models capture the difficult wake recirculation region behind the package reasonably well, although for packages whose heights cause a high degree of recirculation behind the package the SST model appears to struggle. The paper also demonstrates the sensitivity of the models to changes in the mesh density; this study is aimed specifically at thermal design engineers as mesh independent simulations are rarely conducted in an industrial environment.

Predicting the likely impact of aircraft post crash fire on aircraft evacuation using fire and evacuation simulation

The SMARTFIRE Computational Fluid Dynamics (CFD) fire field model has successfully reproduced the observed characteristics including measured temperatures, species concentrations and time to flashover for a post-crash fire experiment conducted by the FAA within their C-133 cabin test facility. In this test only one exit was open in order to provide ventilation for the developing cabin fire. In real post-crash fires, many exits are likely to be open as passangers attempt to evacuate. In this paper, the likely impacts on evacuation of a post-crash fire in which various exiting combinations are available are investigated. The fire scenario, investigated using the SMARTFIRE software, is based on the C-133 experiment but with a fully furnished cabin and with four different exit availability options. The fire data is imported into the airEXODUS evacuation simulation software and the resulting evacuations examined. The combined fire and evacuation analysis reveals that even though the aircraft configuration is predicted to comfortably satisfy the evacuation certification requirement, when fire is included, a number of casualties result, even from the certification compliant exit configuration.

Emission lines of [Cl II] in the optical spectra of gaseous nebulae

Recent R-matrix calculations of electron impact excitation rates among the 3s(2)3p(4) levels of Cl II are used to derive the nebular emission-line intensity ratios R-1=I(6161.8 Angstrom)/I(8578.7 Angstrom) and R-2=I(6161.8 Angstrom)/I(9123.6 Angstrom) as a function of electron temperature (T-e) and density (N-e). The ratios are found to be very sensitive to changes in T-e but not N-e for densities lower than 10(5) cm(-3). Hence, they should, in principle, provide excellent optical T-e diagnostics for planetary nebulae. The observed values of R-1 and R-2 for the planetary nebulae NGC 6741 and IC 5117, measured from spectra obtained with the Hamilton echelle spectrograph on the 3 m Shane Telescope, imply temperatures in excellent agreement with those derived from other diagnostic lines formed in the same region of the nebula as [Cl II]. This provides some observational support for the accuracy of the [Cl II] line ratio calculations and hence the atomic data on which they are based. The [Cl II] 8578.7 and 9123.6 Angstrom lines are identified for the first time (to our knowledge) in a high-resolution spectrum of the symbiotic star RR Telescopii, obtained with the University College London Echelle Spectrograph on the 3.9 m Anglo- Australian Telescope. However, the 6161.8 Angstrom feature is unfortunately too weak to be identified in the RR Telescopii observations, consistent with its predicted line strength.

Emission-line ratios for [N II] in gaseous nebulae and a comparison between theory and observation

R-matrix calculations of electron impact excitation rates among the 2s(2)2p(2) P-3, D-1, S-1, and 2s2p(3) S-5 levels of N II are presented. These results are used in conjunction with other recent calculations of electron impact excitation rates and Einstein A-coefficients for N II to derive the emission-line ratio: ratio diagrams and where (R-1, R-2) (R-1, R-3), where R-1 = I(5756.2 Angstrom)/I(6549.9 + 6585.2 Angstrom), R-2 = I(2143.5 Angstrom)/I(6549.9 + 6585.2 Angstrom), and R-3 = I(2139.7 Angstrom)/I(6549.9 + 658.2 Angstrom), for a range of electron temperatures (T-e = 5000-20,000 K) and electron densities (N-e = 10(2)-10(7) cm(-3)) appropriate to gaseous nebulae. These diagrams should, in principle, allow the simultaneous determination of T-e and N-e from measurements of the [N II] lines in a spectrum. Plasma parameters deduced for a sample of gaseous nebulae, using observational data obtained from ground-based telescopes plus the International Ultraviolet Explorer and Hubble Space Telescope satellites, are found to show generally excellent internal consistency and to be in good agreement with the values of T-e and N-e estimated from other line ratios. These results provide observational support for the accuracy of the theoretical ratios and hence the atomic data adopted in their derivation. Theoretical ratios are also presented for the infrared line pair R-4 = I(122 mum)/I(205 mum), and the usefulness of R-4 as an electron density diagnostic is briefly discussed.

Effective collision strengths for electron impact excitation of C1III

Effective collision strengths for electron-impact excitation of the phosphorus-like ion Cl III are presented for all fine- structure transitions among the levels arising from the lowest 23 LS states. The collisional cross sections are computed in the multichannel close-coupling R-matrix approximation, where sophisticated configuration-interaction wave functions are used to represent the target states. The 23 LS states are formed from the basis configurations 3s(2)3p(3). 3s3p(4). 3s(2)3p(2)3d, and 3s(2)3p(2)4s, and correspond to 49 fine- structure levels, leading to a total possible 1176 fine- structure transitions. The effective collision strengths. obtained by averaging the electron collision strengths over a Maxwellian distribution of electron velocities. are tabulated in this paper for all 1176 transitions and for electron temperatures in the ranges T(K) = 7500-25.000 and log T(K) = 4.4-5.3. The former range encompasses the temperatures of particular importance for application to gaseous nebulae. while the latter range is more applicable to the study of solar and laboratory-type plasmas. (C) 2001 Academic Press.

Value-oriented impact assessment: the economics of a new approach to impact assessment

Environmental Impact Assessment has gained a prominent position as a tool to evaluate the environmental effects of economic activities. However, all approaches proposed so far use a burden-oriented logic. They concentrate on the different environmental impacts in order to ascertain the overall environmental damage caused by economic activity. This paper argues that such a burden-oriented view is (a) hampered by a series of methodological shortcomings which hinders its widespread use in practice; and (b) is analytically incomplete. The paper proposes a value-oriented approach to impact assessment. For this purpose an economic analysis of the optimal use of environmental and social resources is conducted from both a burden-oriented and a value-oriented standpoint. The basic logic of a value-oriented impact assessment is explained, as well as the resulting economic conditions for an optimal use of resources. In addition, it is shown that value- and burden-oriented approaches are complementary to achieve optimality. Finally, the paper discusses the conditions under which the use of burden- or value-oriented impact assessments is appropriate, respectively.