Robust stabilization of a nonlinear cement mill model

Plugging is well known to be a major cause of instability in industrial cement mills. A simple nonlinear model able to simulate the plugging phenomenon is presented. It is shown how a nonlinear robust controller can be designed in order to fully prevent the mill from plugging.

Using mesoscopic models to design strong and tough biomimetic polymer networks

Using computational modeling, we investigate the mechanical properties of polymeric materials composed of coiled chains, or "globules", which encompass a folded secondary structure and are cross-linked by labile bonds to form a macroscopic network. In the presence of an applied force, the globules can unfold into linear chains and thereby dissipate energy as the network is deformed; the latter attribute can contribute to the toughness of the material. Our goal is to determine how to tailor the labile intra- and intermolecular bonds within the network to produce material exhibiting both toughness and strength. Herein, we use the lattice spring model (LSM) to simulate the globules and the cross-linked network. We also utilize our modified Hierarchical Bell model (MHBM) to simulate the rupture and reforming of N parallel bonds. By applying a tensile deformation, we demonstrate that the mechanical properties of the system are sensitive to the values of N in and N out, the respective values of N for the intra- and intermolecular bonds. We find that the strength of the material is mainly controlled by the value of N out, with the higher value of N out providing a stronger material. We also find that, if N in is smaller than N out, the globules can unfold under the tensile load before the sample fractures and, in this manner, can increase the ductility of the sample. Our results provide effective strategies for exploiting relatively weak, labile interactions (e.g., hydrogen bonding or the thiol/disulfide exchange reaction) in both the intra- and intermolecular bonds to tailor the macroscopic performance of the materials. © 2011 American Chemical Society.

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.

Towards robust design of axial compressors with uncertainty quantification

Operational uncertainties such as throttle excursions, varying inlet conditions and geometry changes lead to variability in compressor performance. In this work, the main operational uncertainties inherent in a transonic axial compressor are quantified to deter- mine their effect on performance. These uncertainties include the effects of inlet distortion, metal expansion, ow leakages and blade roughness. A 3D, validated RANS model of the compressor is utilized to simulate these uncertainties and quantify their effect on polytropic efficiency and pressure ratio. To propagate them, stochastic collocation and sparse pseudospectral approximations are used. We demonstrate that lower-order approximations are sufficient as these uncertainties are inherently linear. Results for epistemic uncertainties in the form of meshing methodologies are also presented. Finally, the uncertainties considered are ranked in order of their effect on efficiency loss. © 2012 AIAA.

The sensitivity of vibration characteristics of reinforced concrete beams under incremental static and cyclic loading

The use of changes in vibration properties for global damage detection and monitoring of existing concrete structures has received great research attention in the last three decades. To track changes in vibration properties experimentally, structures have been artificially damaged by a variety of scenarios. However, this procedure does not represent realistically the whole design-life degradation of concrete structures. This paper presents experimental work on a set of damaged reinforced concrete beams due to different loading regimes to assess the sensitivity of vibration characteristics. Of the total set, three beams were subject to incremental static loading up to failure to simulate overloading, and two beams subject to 15 million loading cycles with varying amplitudes to produce an accelerated whole-life degradation scenario. To assess the vibration behaviour in both cases, swept sine and harmonic excitations were conducted at every damage level. The results show that resonant frequencies are not sensitive enough to damage due to cyclic loading, whereas cosh spectral and root mean square distances are more sensitive, yet more scattered. In addition, changes in non-linearity follow a softening trend for beams under incremental static loading, whilst they are significantly inconsistent for beams under cyclic loading. Amongst all examined characteristics, changes in modal stiffness are found to be most sensitive to damage and least scattered, but modal stiffness is tedious to compute due mainly to the difficulty of constructing restoring force surfaces from field measurements. © (2013) Trans Tech Publications.

Experimental study of surface roughness noise

A turbulent boundary-layer flow over a rough wall generates a dipole sound field as the near-field hydrodynamic disturbances in the turbulent boundary-layer scatter into radiated sound at small surface irregularities. In this paper, phased microphone arrays are applied to the experimental study of surface roughness noise. The radiated sound from two rough plates and one smooth plate in an open jet is measured at three streamwise locations, and the beamforming source maps demonstrate the dipole directivity. Higher source strengths can be observed in the rough plates than the smooth plate, and the rough plates also enhance the trailing-edge noise. A prediction scheme in previous theoretical work is used to describe the strength of a distribution of incoherent dipoles over the rigid plate and to simulate the sound detected by the microphone array. Source maps of measurement and simulation exhibit encouraging similarities in both source pattern and source strength, which confirms the dipole nature and the predicted magnitude of roughness noise. The simulations underestimate the streamwise gradient of the source strengths and overestimate the source strengths at the highest frequency. © 2007 by Yu Liu and Ann P. Dowling.

Oxygen carrier dispersion in inert packed beds to improve performance in chemical looping combustion

Various packed beds of copper-based oxygen carriers (CuO on Al2O3) were tested over 100 cycles of low temperature (673K) Chemical Looping Combustion (CLC) with H2 as the fuel gas. The oxygen carriers were uniformly mixed with alumina (Al2O3) in order to investigate the level of separation necessary to prevent agglomeration. It was found that a mass ratio of 1:6 oxygen carrier to alumina gave the best performance in terms of stable, repeating hydrogen breakthrough curves over 100 cycles. In order to quantify the average separation achieved in the mixed packed beds, two sphere-packing models were developed. The hexagonal close-packing model assumed a uniform spherical packing structure, and based the separation calculations on a hypergeometric probability distribution. The more computationally intensive full-scale model used discrete element modelling to simulate random packing arrangements governed by gravity and contact dynamics. Both models predicted that average 'nearest neighbour' particle separation drops to near zero for oxygen carrier mass fractions of x≥0.25. For the packed bed systems studied, agglomeration was observed when the mass fraction of oxygen carrier was above this threshold. © 2013 Elsevier B.V.

Large eddy simulations in low-pressure turbines: Effect of wakes at elevated free-stream turbulence

The transition of a separated shear layer over a flat plate, in the presence of periodic wakes and elevated free-stream turbulence (FST), is numerically investigated using Large Eddy Simulation (LES). The upper wall of the test section is inviscid and specifically contoured to impose a streamwise pressure distribution over the flat plate to simulate the suction surface of a low-pressure turbine (LPT) blade. Two different distributions representative of a 'high-lift' and an 'ultra high-lift' turbine blade are examined. Results obtained from the current LES compare favourably with the extensive experimental data previously obtained for these configurations. The LES results are then used to further investigate the flow physics involved in the transition process.In line with experimental experience, the benefit of wakes and FST obtained by suppressing the separation bubble, is more pronounced in 'ultra high-lift' design when compared to the 'high-lift' design. Stronger 'Klebanoff streaks' are formed in the presence of wakes when compared to the streaks due to FST alone. These streaks promoted much early transition. The weak Klebanoff streaks due to FST continued to trigger transition in between the wake passing cycles.The experimental inference regarding the origin of Klebanoff streaks at the leading edge has been confirmed by the current simulations. While the wake convects at local free-stream velocity, its impression in the boundary layer in the form of streaks convects much slowly. The 'part-span' Kelvin-Helmholtz structures, which were observed in the experiments when the wake passes over the separation bubble, are also captured. The non-phase averaged space-time plots manifest that reattachment is a localized process across the span unlike the impression of global reattachment portrayed by phase averaging. © 2013 Elsevier Inc.

Adapting data processing to compare model and experiment accurately: A discrete element model and magnetic resonance measurements of a 3d cylindrical fluidized bed

Discrete element modeling is being used increasingly to simulate flow in fluidized beds. These models require complex measurement techniques to provide validation for the approximations inherent in the model. This paper introduces the idea of modeling the experiment to ensure that the validation is accurate. Specifically, a 3D, cylindrical gas-fluidized bed was simulated using a discrete element model (DEM) for particle motion coupled with computational fluid dynamics (CFD) to describe the flow of gas. The results for time-averaged, axial velocity during bubbling fluidization were compared with those from magnetic resonance (MR) experiments made on the bed. The DEM-CFD data were postprocessed with various methods to produce time-averaged velocity maps for comparison with the MR results, including a method which closely matched the pulse sequence and data processing procedure used in the MR experiments. The DEM-CFD results processed with the MR-type time-averaging closely matched experimental MR results, validating the DEM-CFD model. Analysis of different averaging procedures confirmed that MR time-averages of dynamic systems correspond to particle-weighted averaging, rather than frame-weighted averaging, and also demonstrated that the use of Gaussian slices in MR imaging of dynamic systems is valid. © 2013 American Chemical Society.

Mitigation of lake eutrophication: Loosen nitrogen control and focus on phosphorus abatement

w Traditionally, nitrogen control is generally considered an important component of reducing lake eutrophication and cyanobacteria blooms. However, this viewpoint is refuted recently by researchers in China and North America. In the present paper, the traditional viewpoint of nitrogen control is pointed out to lack a scientific basis: the N/P hypothesis is just a subjective assumption; bottle bioassay experiments fail to simulate the natural process of nitrogen fixation. Our multi-year comparative research in more than 40 Yangtze lakes indicates that phosphorus is the key factor determining phytoplankton growth regardless of nitrogen concentrations and that total phytoplankton biomass is determined by total phosphorus and not by total nitrogen concentrations. These results imply that, in the field, nitrogen control will not decrease phytoplankton biomass. This finding is supported by a long-term whole-lake experiment from North America. These outcomes can be generalized in terms that a reduction in nitrogen loading may not decrease the biomass of total phytoplankton as it can stimulate blooms of nitrogen-fixing cyanobacteria. To mitigate eutrophication, it is not nitrogen but phosphorus that should be reduced, unless nitrogen concentrations are too high to induce direct toxic impacts on human beings or other organisms. Finally, details are provided on how to reduce controls on nitrogen and how to mitigate eutrophication. (C) 2009 National Natural Science Foundation of China and Chinese Academy of Sciences. Published by Elsevier Limited and Science in China Press. All rights reserved.

Plankton Community Succession in Artificial Systems Subjected to Cyanobacterial Blooms Removal using Chitosan-Modified Soils

Using artificial systems to simulate natural lake environments with cyanobacterial blooms, we investigated plankton community succession by polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE) fingerprinting and morphological method. With this approach, we explored potential ecological effects of a newly developed cyanobacterial blooms removal method using chitosan-modified soils. Results of PCR-DGGE and morphological identification showed that plankton communities in the four test systems were nearly identical at the beginning of the experiment. After applying the newly developed and standard removal methods, there was a shift in community composition, but neither chemical conditions nor plankton succession were significantly affected by the cyanobacteria removal process. The planted Vallisneria natans successfully recovered after cyanobacteria removal, whereas that in the box without removal process did not. Additionally, canonical correspondence analysis indicated that other than for zooplankton abundance, total phosphorus was the most important environmental predictor of planktonic composition. The present study and others suggest that dealing with cyanobacteria removal using chitosan-modified soils can play an important role in controlling cyanobacterial blooms in eutrophicated freshwater systems.

On the application of LES to seal geometries

Five Large Eddy Simulation (LES) and hybrid RANS-NLES (Reynolds-Averaged Navier-Stokes-Numerical-LES) methods are used to simulate flow through a labyrinth seal geometry and are contrasted with RANS solutions. Results show that LES and RANS-NLES is capable of accurately predicting flow behaviour of two seal flows with a scatter of less than 5 %. RANS solutions show the potential to perform poorly for the turbulence models tested. LES and hybrid RANS-NLES are found to be consistent and in agreement with measurements, providing a flexible numerical platform for design investigations. It also allows greater flow physics insights. © Springer Science+Business Media Dordrecht 2013.

First principles interatomic potential for tungsten based on Gaussian process regression

An accurate description of atomic interactions, such as that provided by first principles quantum mechanics, is fundamental to realistic prediction of the properties that govern plasticity, fracture or crack propagation in metals. However, the computational complexity associated with modern schemes explicitly based on quantum mechanics limits their applications to systems of a few hundreds of atoms at most. This thesis investigates the application of the Gaussian Approximation Potential (GAP) scheme to atomistic modelling of tungsten - a bcc transition metal which exhibits a brittle-to-ductile transition and whose plasticity behaviour is controlled by the properties of $\frac{1}{2} \langle 111 \rangle$ screw dislocations. We apply Gaussian process regression to interpolate the quantum-mechanical (QM) potential energy surface from a set of points in atomic configuration space. Our training data is based on QM information that is computed directly using density functional theory (DFT). To perform the fitting, we represent atomic environments using a set of rotationally, permutationally and reflection invariant parameters which act as the independent variables in our equations of non-parametric, non-linear regression. We develop a protocol for generating GAP models capable of describing lattice defects in metals by building a series of interatomic potentials for tungsten. We then demonstrate that a GAP potential based on a Smooth Overlap of Atomic Positions (SOAP) covariance function provides a description of the $\frac{1}{2} \langle 111 \rangle$ screw dislocation that is in agreement with the DFT model. We use this potential to simulate the mobility of $\frac{1}{2} \langle 111 \rangle$ screw dislocations by computing the Peierls barrier and model dislocation-vacancy interactions to QM accuracy in a system containing more than 100,000 atoms.

Weighted least squares with expectation-maximization algorithm for burst detection in U.K. water distribution systems

Flow measurement data at the district meter area (DMA) level has the potential for burst detection in the water distribution systems. This work investigates using a polynomial function fitted to the historic flow measurements based on a weighted least-squares method for automatic burst detection in the U.K. water distribution networks. This approach, when used in conjunction with an expectationmaximization (EM) algorithm, can automatically select useful data from the historic flow measurements, which may contain normal and abnormal operating conditions in the distribution network, e.g., water burst. Thus, the model can estimate the normal water flow (nonburst condition), and hence the burst size on the water distribution system can be calculated from the difference between the measured flow and the estimated flow. The distinguishing feature of this method is that the burst detection is fully unsupervised, and the burst events that have occurred in the historic data do not affect the procedure and bias the burst detection algorithm. Experimental validation of the method has been carried out using a series of flushing events that simulate burst conditions to confirm that the simulated burst sizes are capable of being estimated correctly. This method was also applied to eight DMAs with known real burst events, and the results of burst detections are shown to relate to the water company's records of pipeline reparation work. © 2014 American Society of Civil Engineers.

The influence of boundary conditions on tip leakage flow

Most of the current understanding of tip leakage flows has been derived from detailed cascade experiments. However, the cascade model is inherently approximate since it is difficult to simulate the boundary conditions present in a real machine, particularly the secondary flows convecting from the upstream stator row and the relative motion of the casing and blade. This problem is further complicated when considering the high pressure turbine rotors of aero engines, where the high Mach numbers must also be matched in order to correctly model the aerodynamics and heat transfer. More realistic tests can be performed on high-speed turbines, but the experimental fidelity and resolution achievable in such set-ups is limited. In order to examine the differences between cascade models and real-engine behavior, the influence of boundary conditions on the tip leakage flow in an unshrouded high pressure turbine rotor is investigated using RANS calculations. This study examines the influence of the rotor inlet condition and relative casing motion. A baseline calculation with a simplified inlet condition and no relative endwall motion exhibits similar behavior to cascade studies. Only minor changes to the leakage flow are induced by introducing either a more realistic inlet condition or relative casing motion. However when both of these conditions are applied simultaneously the pattern of leakage flow is very different, with ingestion of flow over much of the early suction surface. The paper explores the physical processes driving this change and the impact on leakage losses and modeling requirements. Copyright © 2013 by ASME.