The effect of surface hardening on the elastic shakedown of elliptical contacts

The effect of varying both the aspect ratio and the coefficient of friction of contacts with elliptical geometry on their elastic shakedown performance has been examined theoretically for surfaces with two types of subsurface hardness or strength profiles. In stepwise hardening the hard layer is of uniform strength while in linear hardening its strength reduces from a maximum at the surface to that of the core at the base of the hardened layer. The shakedown load is expressed as the ratio of the maximum Hertzian pressure to the strength of the core material. As the depth of hardening, expressed as a multiple of the elliptical semi-axis, is increased so the potential shakedown load increases from a level that is appropriate to a uniform half-space of unhardened material to a value reflecting the hardness of the surface and near-surface material. In a step-hardened material, the shakedown limit for a surface 'pummelled' by the passage of a sequence of such loads reaches a cut-off or plateau value, which cannot be exceeded by further increases in hardening depth irrespective of the value of the friction coefficient. For a linear-hardened material the corresponding plateau is approached asymptotically. The work confirms earlier results on the upper bounds on shakedown of both point and line contacts and provides numerical values of shakedown loads for intermediate geometries. In general, the case depth required to achieve a given shakedown limit reduces in moving from a transversely moving nominal line load to an axisymmetric point load.

The contact of rough surfaces carrying pressure sensitive boundary layers

When two rough surfaces are loaded together it is well known that the area of true contact is very much smaller then the geometric area and that, consequently, local contact pressures are very much greater than the nominal value. If the asperities on each surface can be thought of as possessing smooth summits and each of the solids is elastically isotropic then the pressure distribution will consist of a series of small, but severe, Hertzian patches. However, if one of both of the surfaces in question is protected by a boundary layer then both the number and dimensions of these patches, and the form of the pressure distribution within them, will be modified. Recent experimental evidence from studies using both Atomic Force Microscopy and micro-tribometry suggests that boundary films produced by the action of commercial anti-wear additives, such as ZDTP, exhibit mechanical properties, which are affected by local values of pressure. These changes bring about further modifications to local conditions. These effects have been explored in a numerical model of rough surface contact and the implications for the mechanisms of surface distress and wear are discussed. © 2000 Elsevier Science B.V. All rights reserved.

Friction of sliding surfaces carrying adsorbed lubricant layers

In concentrated contacts the behaviour of lubricants is much modified by the high local pressures: changes can arise both from molecular ordering within the very thin film lubricant layers present at the interface as well as from the deposition on the component surfaces of more solid-like polymeric boundary layers. These 'third bodies' separating the solid surfaces may have rheological or mechanical properties very different from those observed in the bulk. Classical elasto-hydrodynamic theory considers the entrapped lubricant to exhibit a piezo-viscous behaviour while the conventional picture of more solid boundary lubricant layers views their shear strength r as being linearly dependent on local pressure p, so that T = TO + ap where TO and a are constants. If TO is relatively small, then the coefficient of friction \i = T Ip ~ a and so Amonton's laws are recovered. However, the properties of adsorbed or deposited surface films, or indeed other third bodies such as debris layers, may be more complex than this. A preliminary study has looked quantitatively at the influence of the pressure dependence of the shear strength of any surface layer on the overall friction coefficient of a contact which is made up of an array of asperities whose height varies in a Gaussian manner. Individual contact points may be elastic or plastic. The analysis results in plots of coefficient of friction versus the service or load parameter PIH&NRa where P is the nominal pressure on the contact, HS the hardness of the deforming surface, N the asperity density, R the mean radius of curvature of the asperities, and a is the standard deviation of their height distribution. In principle, any variation oft withp can be incorporated into the model; however, in this initial study we have used data on colloidal suspensions from the group at the Ecole Centrale de Lyon as well as examining the effect of functional relationships of somewhat greater complexity than a simple linear form. Results of the analysis indicate that variations in fj. are possible as the load is varied which depend on the statistical spread of behaviour at individual asperity contacts. The value of this analysis is that it attempts to combine the behaviour of films on the molecular scale with the topography of real engineering surfaces and so give an indication of the effects at the full-size or macro-scale that can be achieved by chemical or molecular surface engineering.

Velocity distributions for low pressure turbines

A parametric set of velocity distributions has been investigated using a flat plate experiment. Three different diffusion factors and peak velocity locations were tested. These were designed to mimic the suction surfaces of Low Pressure (LP) turbine blades. Unsteady wakes, inherent in real turbomachinery flows, were generated using a moving bar mechanism. A turbulence grid generated a freestream turbulence level that is believed to be typical of LP turbines. Measurements were taken across a Reynolds number range of 50,000-220,000 at three reduced frequencies (0.314, 0.628, 0.942). Boundary layer traverses were performed at the nominal trailing edge using a Laser Doppler Anemometry system and hot-films were used to examine the boundary layer behaviour along the surface. For every velocity distribution tested, the boundary layer separated in the diffusing flow downstream of the peak velocity. The loss production is dominated by the mixing in the reattachment process, mixing in the turbulent boundary layer downstream of reattachment and the effects of the unsteady interaction between the wakes and the boundary layer. A sensitive balance governs the optimal location of peak velocity on the surface. Moving the velocity peak forwards on the blade was found to be increasingly beneficial when bubblegenerated losses are high, i.e. at low Reynolds number, at low reduced frequency and at high levels of diffusion. Copyright © 2008 by ASME.

The internal structure of igniting turbulent sprays as revealed by complex chemistry DNS

A parametric study of spark ignition in a uniform monodisperse turbulent spray is performed with complex chemistry three-dimensional Direct Numerical Simulations in order to improve the understanding of the structure of the ignition kernel. The heat produced by the kernel increases with the amount of fuel evaporated inside the spark volume. Moreover, the heat sink by evaporation is initially higher than the heat release and can have a negative effect on ignition. With the sprays investigated, heat release occurs over a large range of mixture fractions, being high within the nominal flammability limits and finite but low below the lean flammability limit. The burning of very lean regions is attributed to the diffusion of heat and species from regions of high heat release, and from the spark, to lean regions. Two modes of spray ignition are reported. With a relatively dilute spray, nominally flammable material exists only near the droplets. Reaction zones are created locally near the droplets and have a non-premixed character. They spread from droplet to droplet through a very lean interdroplet spacing. With a dense spray, the hot spark region is rich due to substantial evaporation but the cold region remains lean. In between, a large surface of flammable material is generated by evaporation. Ignition occurs there and a large reaction zone propagates from the rich burned region to the cold lean region. This flame is wrinkled due to the stratified mixture fraction field and evaporative cooling. In the dilute spray, the reaction front curvature pdf contains high values associated with single droplet combustion, while in the dense spray, the curvature is lower and closer to the curvature associated with gaseous fuel ignition kernels. © 2011 The Combustion Institute.

Multiplexed model predictive control

This paper proposes a form of MPC in which the control variables are moved asynchronously. This contrasts with most MIMO control schemes, which assume that all variables are updated simultaneously. MPC outperforms other control strategies through its ability to deal with constraints. This requires on-line optimization, hence computational complexity can become an issue when applying MPC to complex systems with fast response times. The Multiplexed MPC (MMPC) scheme described in this paper solves the MPC problem for each subsystem sequentially, and updates subsystem controls as soon as the solution is available, thus distributing the control moves over a complete update cycle. The resulting computational speed-up allows faster response to disturbances, which may result in improved performance, despite finding sub-optimal solutions to the original problem. This paper describes nominal and robust MMPC, states some stability results, and demonstrates the effectiveness of MMPC through two examples. © 2011 Elsevier Ltd. All rights reserved.

Multi-objective optimization for multiphase orbital rendezvous missions

A multi-objective optimization approach was proposed for multiphase orbital rendezvous missions and validated by application to a representative numerical problem. By comparing the Pareto fronts obtained using the proposed method, the relationships between the three objectives considered are revealed, and the influences of other mission parameters, such as the sensors' field of view, can also be analyzed effectively. For multiphase orbital rendezvous missions, the tradeoff relationships between the total velocity increment and the trajectory robustness index as well as between the total velocity increment and the time of flight are obvious and clear. However, the tradeoff relationship between the time of flight and the trajectory robustness index is weak, especially for the four- and five-phase missions examined. The proposed approach could be used to reorganize a stable rendezvous profile for an engineering rendezvous mission, when there is a failure that prevents the completion of the nominal mission.

Open cycle thorium–uranium-fuelled nuclear energy systems

The potential for countries that currently have a nominal nuclear energy infrastructure to adopt thorium–uranium-fuelled nuclear energy systems, using a once-through ‘open’ nuclear fuel cycle, has been suggested by the International Atomic Energy Agency. This review paper highlights generation II, III and III+ nuclear energy technologies that could potentially adopt an open thorium–uranium fuel cycle and qualitatively highlights the main differences between the open thorium–uranium and open uranium fuel cycles.

Slow peaking and low-gain designs for global stabilization of nonlinear systems

This paper presents an analysis of the slow-peaking phenomenon, a pitfall of low-gain designs that imposes basic limitations to large regions of attraction in nonlinear control systems. The phenomenon is best understood on a chain of integrators perturbed by a vector field up(x, u) that satisfies p(x, 0) = 0. Because small controls (or low-gain designs) are sufficient to stabilize the unperturbed chain of integrators, it may seem that smaller controls, which attenuate the perturbation up(x, u) in a large compact set, can be employed to achieve larger regions of attraction. This intuition is false, however, and peaking may cause a loss of global controllability unless severe growth restrictions are imposed on p(x, u). These growth restrictions are expressed as a higher order condition with respect to a particular weighted dilation related to the peaking exponents of the nominal system. When this higher order condition is satisfied, an explicit control law is derived that achieves global asymptotic stability of x = 0. This stabilization result is extended to more general cascade nonlinear systems in which the perturbation p(x, v) v, v = (ξ, u) T, contains the state ξ and the control u of a stabilizable subsystem ξ = a(ξ, u). As an illustration, a control law is derived that achieves global stabilization of the frictionless ball-and-beam model.

Fertile-free annular fuel for plutonium recycling

Up to 50% increase in the power density of the existing pressurized water reactor (PWR)-type reactors can be achieved by the use of internally and externally cooled annular fuel geometry. As a result, the accumulated stock-piles of Pu, especially if incorporated infertile-free inert matrix, can be burnt at a substantially higher rate as compared with the conventional mixed oxide-fueled reactors operating at standard power density. In this work, we explore the basic feasibility of a PWR core fully loaded with Pu incorporated infertile-free fuel of annular internally and externally cooled geometry and operating at 150% of nominal power density. We evaluate basic burnable poison designs, fuel management strategies, and reactivity feedback coefficients. The three-dimensional full core neutronic analysis performed with Studsvik Core Management System showed that the design of such a Pu-loaded annular fuel core is feasible but significantly more challenging than the Pu fertile-free core with solid fuel pins operating at nominal power density. The main difficulty arises from the fact that the annular fuel core requires at least 50% higher initial Pu loading in order to maintain the standard fuel cycle length of 18 months. Such a high Pu loading results in hardening of the neutron spectrum and consequent reduction in reactivity worth of all reactivity control mechanisms and, in some cases, positive moderator temperature coefficient (MTC). The use of isotopically enriched Gd and Er burnable poisons was found to be beneficial with respect to maximizing Pu burnup and reducing power peaking factors. Overall, the annular fertile-free Pu-loaded high-power-density core appears to be feasible, although it still has relatively high power peaking and potential for slightly positive MTC at beginning of cycle. However, we estimate that limiting the power density to 140% of the nominal case would assure acceptable core power peaking and negative MTC at all times during the cycle.

Active control of high frequency vibration in uncertain structures

The active suppression of structural vibration is normally achieved by either feedforward or feedback control. In the absence of a suitable reference signal feedforward control cannot be employed and feedback control is the only viable approach. Conventional feedback control algorithms (e.g. LQR and LQG) are designed on the basis of a mathematical model of the system and ideally the performance of the system should be robust against uncertainties in this model. The aim of this paper is to numerically investigate the robustness of LQR and LQG algorithms by designing the controller for a nominal system, and then assessing (via Monte Carlo simulation) the effects of uncertainties in the system. The ultimate concern is with the control of high frequency vibrations, where the short wavelength of the structural deformation induces a high sensitivity to imperfection. It is found that standard algorithms such as LQR and LQG are generally unfeasible for this case. This leads to a consideration of design strategies for the robust active control of high frequency vibrations. The system chosen for the numerical simulation concerns two coupled plates, which are randomized by the addition of point masses at random locations.

All-chemical vapor deposited graphene/silicon nitride TFTs

All-chemical vapor deposited silicon nitride / monolayer graphene TFTs have been fabricated. Polychromatic Raman spectroscopy shows high quality monolayer graphene channels with uniform coverage and significant interfacial doping at the source-drain contacts. Nominal mobilities of approximately 1900 cm 2V-1s-1 have been measured opening up a potentially useful platform for analogue and RFR-based applications fabricated through allchemical vapor deposition processes. © The Electrochemical Society.

Predicting the response of structures to transient shock loading.

This work concerns the prediction of the response of an uncertain structure to a load of short duration. Assuming an ensemble of structures with small random variations about a nominal form, a mean impulse response can be found using only the modal density of the structure. The mean impulse response turns out to be the same as the response of an infinite structure: the response is calculated by taking into account the direct field only, without reflections. Considering the short duration of an impulsive loading, the approach is reasonable before the effect of the reverberant field becomes important. The convolution between the mean impulse response and the shock loading is solved in discrete time to calculate the response at the driving point and at remote points. Experimental and numerical examples are presented to validate the theory presented for simple structures such as beams, plates, and cylinders.

Inkjet printing of weakly elastic polymer solutions

Fluid assessment methods, requiring small volumes and avoiding the need for jetting, are particularly useful in the design of functional fluids for inkjet printing applications. With the increasing use of complex (rather than Newtonian) fluids for manufacturing, single frequency fluid characterisation cannot reliably predict good jetting behaviour, owing to the range of shearing and extensional flow rates involved. However, the scope of inkjet fluid assessments (beyond achievement of a nominal viscosity within the print head design specification) is usually focused on the final application rather than the jetting processes. The experimental demonstration of the clear insufficiency of such approaches shows that fluid jetting can readily discriminate between fluids assessed as having similar LVE characterisation (within a factor of 2) for typical commercial rheometer measurements at shearing rates reaching 104rads-1.Jetting behaviour of weakly elastic dilute linear polystyrene solutions, for molecular weights of 110-488. kDa, recorded using high speed video was compared with recent results from numerical modelling and capillary thinning studies of the same solutions.The jetting images show behaviour ranging from near-Newtonian to "beads-on-a-string". The inkjet printing behaviour does not correlate simply with the measured extensional relaxation times or Zimm times, but may be consistent with non-linear extensibility L and the production of fully extended polymer molecules in the thinning jet ligament.Fluid test methods allowing a more complete characterisation of NLVE parameters are needed to assess inkjet printing feasibility prior to directly jetting complex fluids. At the present time, directly jetting such fluids may prove to be the only alternative. © 2014 The Authors.

Output disturbance rejection using parallel model predictive control

The solution time of the online optimization problems inherent to Model Predictive Control (MPC) can become a critical limitation when working in embedded systems. One proposed approach to reduce the solution time is to split the optimization problem into a number of reduced order problems, solve such reduced order problems in parallel and selecting the solution which minimises a global cost function. This approach is known as Parallel MPC. The potential capabilities of disturbance rejection are introduced using a simulation example. The algorithm is implemented in a linearised model of a Boeing 747-200 under nominal flight conditions and with an induced wind disturbance. Under significant output disturbances Parallel MPC provides a significant improvement in performance when compared to Multiplexed MPC (MMPC) and Linear Quadratic Synchronous MPC (SMPC). © 2013 IEEE.