Optimal H2 order-one reduction by solving eigenproblems for polynomial equations

A method is given for solving an optimal H2 approximation problem for SISO linear time-invariant stable systems. The method, based on constructive algebra, guarantees that the global optimum is found; it does not involve any gradient-based search, and hence avoids the usual problems of local minima. We examine mostly the case when the model order is reduced by one, and when the original system has distinct poles. This case exhibits special structure which allows us to provide a complete solution. The problem is converted into linear algebra by exhibiting a finite-dimensional basis for a certain space, and can then be solved by eigenvalue calculations, following the methods developed by Stetter and Moeller. The use of Buchberger's algorithm is avoided by writing the first-order optimality conditions in a special form, from which a Groebner basis is immediately available. Compared with our previous work the method presented here has much smaller time and memory requirements, and can therefore be applied to systems of significantly higher McMillan degree. In addition, some hypotheses which were required in the previous work have been removed. Some examples are included.

Three-dimensional SBLI control for transonic airfoils

The application of shock control to transonic airfoils and wings has been demonstrated widely to have the potential to reduce wave drag. Most of the suggested control devices are two-dimensional, that is they are of uniform geometry in spanwise direction. Examples of such techniques include contour bumps and passive control. Recently it has been observed that a spanwise array of discrete three-dimensional controls can have similar benefits but also offer advantages in terms of installation complexity and drag. This paper describes research carried out in Cambridge into various three-dimensional devices, such as slots, grooves and bumps. In all cases the control device is applied to the interaction of a normal shock wave (M=1.3) with a turbulent boundary layer. Theoretical considerations are proposed to determine how such fundamental experiments can provide estimates of control performance on a transonic wing. The potential of each class of three-dimensional device for wave drag reduction on airfoils is discussed and surface bumps in particular are identified as offering potential drag savings for typical transonic wing applications under cruise conditions.

Shock / boundary-layer interaction control using three-dimensional bumps for transonic wings

Three-dimensional bumps have been developed and investigated, aiming at the two major objectives of shock-wave / boundary-layer interaction control, i.e. drag reduction and suppression of separation, simultaneously. An experimental investigation has been conducted for a default rounded bump in channel now at University of Cambridge and a computational study has been performed for a spanwise series of rounded bumps mounted on a transonic aerofoil at University of Stuttgart. Observed in both cases are wave drag reduction owing to A-shock structures produced by three-dimensional surface bumps and mild control effects on the boundary layer. The effects of rough surface and tall extension have been investigated as well as several geometric variations and multiple bump configurations. A double configuration of narrow rounded bumps has been found to best perform amongst the tested, considerably reducing wave drag through a well-established A-shock structure with little viscous penalty and thus achieving substantial overall drag reduction. Counter-rotating streamwise vortex pairs have been produced by some configurations as a result of local flow separation, but they have been observed to be confined in relatively narrow wake regions, expected to be beneficial in suppressing large-scale separation under off-design condition despite increase of viscous drag. On the whole a large potential of three-dimensional control with discrete rounded bumps has been demonstrated both experimentally and numerically, and experimental investigation of bumps fitted on a transonic aerofoil or wing is suggested toward practical application.

Multi-dimensional conditional moment closure modelling applied to a heavy-duty common-rail diesel engine

A multi-dimensional combustion code implementing the Conditional Moment Closure turbulent combustion model interfaced with a well-established RANS two- phase flow field solver has been employed to study a broad range of operating conditions for a heavy duty direct-injection common-rail Diesel engine. These conditions include different loads (25%, 50%, 75% and full load) and engine speeds (1250 and 1830 RPM) and, with respect to the fuel path, different injection timings and rail pressures. A total of nine cases have been simulated. Excellent agreement with experimental data has been found for the pressure traces and the heat release rates, without adjusting any model constants. The chemical mechanism used contains a detailed NOx sub-mechanism. The predicted emissions agree reasonably well with the experimental data considering the range of operating points and given no adjustments of any rate constants have been employed. In an effort to identify CPU cost reduction potential, various dimensionality reduction strategies have been assessed. Furthermore, the sensitivity of the predictions with respect to resolution in particular relating to the CMC grid has been investigated. Overall, the results suggest that the presented modelling strategy has considerable predictive capability concerning Diesel engine combustion without requiring model constant calibration based on experimental data. This is true particularly for the heat release rates predictions and, to a lesser extent, for NOx emissions where further progress is still necessary. © 2009 SAE International.

Influence of film cooling hole angles and geometries on aerodynamic loss and net heat flux reduction

Turbine design engineers have to ensure that film cooling can provide sufficient protection to turbine blades from the hot mainstream gas, while keeping the losses low. Film cooling hole design parameters include inclination angle (α), compound angle (β ), hole inlet geometry and hole exit geometry. The influence of these parameters on aerodynamic loss and net heat flux reduction is investigated, with loss being the primary focus. Low-speed flat plate experiments have been conducted at momentum flux ratios of IR = 0.16, 0.64 and 1.44. The film cooling aerodynamic mixing loss, generated by the mixing of mainstream and coolant, can be quantified using a three-dimensional analytical model that has been previously reported by the authors. The model suggests that for the same flow conditions, the aerodynamic mixing loss is the same for holes with different α and β but with the same angle between the mainstream and coolant flow directions (angle κ). This relationship is assessed through experiments by testing two sets of cylindrical holes with different α and β : one set with κ = 35°, another set with κ = 60°. The data confirm the stated relationship between α, β, κ and the aerodynamic mixing loss. The results show that the designer should minimise κ to obtain the lowest loss, but maximise β to achieve the best heat transfer performance. A suggestion on improving the loss model is also given. Five different hole geometries (α =35.0°, β =0°) were also tested: cylindrical hole, trenched hole, fan-shaped hole, D-Fan and SD-Fan. The D-Fan and the SD-Fan have similar hole exits to the fan-shaped hole but their hole inlets are laterally expanded. The external mixing loss and the loss generated inside the hole are compared. It was found that the D-Fan and the SD-Fan have the lowest loss. This is attributed to their laterally expanded hole inlets, which lead to significant reduction in the loss generated inside the holes. As a result, the loss of these geometries is ≈ 50 % of the loss of the fan-shaped hole at IR = 0.64 and 1.44. Copyright © 2011 by ASME.

Recycle of transuranium elements in light water reactors for reduction of geological storage requirements

This paper presents results of a feasibility study aimed at developing a zero-transuranic-discharge fuel cycle based on the U-Th-TRU ternary cycle. The design objective is to find a fuel composition (mixture of thorium, enriched uranium, and recycled transuranic components) and fuel management strategy resulting in an equilibrium charge-discharge mass flow. In such a fuel cycle scheme, the quantity and isotopic vector of the transuranium (TRU) component is identical at the charge and discharge time points, thus allowing the whole amount of the TRU at the end of the fuel irradiation period to be separated and reloaded into the following cycle. The TRU reprocessing activity losses are the only waste stream that will require permanent geological storage, virtually eliminating the long-term radiological waste of the commercial nuclear fuel cycle. A detailed three-dimensional full pressurized water reactor (PWR) core model was used to analyze the proposed fuel composition and management strategy. The results demonstrate the neutronic feasibility of the fuel cycle with zero-TRU discharge. The amount of TRU and enriched uranium loaded reach equilibrium after about four TRU recycles. The reactivity coefficients were found to be within a range typical for a reference PWR core. The soluble boron worth is reduced by a factor of ∼2 from a typical PWR value. Nevertheless, the results indicate the feasibility of an 18-month fuel cycle design with an acceptable beginning-of-cycle soluble boron concentration even without application of burnable poisons.

Use of axially graded burnable boron for hot-spot temperature reduction in a pressurized water reactor core

Shortly after the loading of a pressurized water reactor (PWR) core, the axial power distribution in fresh fuel has already attained the characteristic, almost flat shape, typical of burned fuel. At beginning of cycle (BOC), however, the axial distribution is centrally peaked. In assemblies hosting uniform burnable boron rods, this BOC peaking is even more pronounced. A reduction in the axial peaking is today often achieved by shortening the burnable boron rods by some 30 cm at each edge. It is shown that a two-zone grading of the boron rod leads, in a representative PWR cycle, to a reduction of the hot-spot temperature of approximately 70 °C, compared with the nongraded case. However, with a proper three-zone grading of the boron rod, an additional 20 °C may be cut off the hot-spot temperature. Further, with a slightly skewed application of this three-zone grading, an additional 50 °C may be cut off. The representative PWR cycle studied was cycle 11 of the Indian Point 2 station, with a simplification in the number of fuel types and in the burnup distribution. The analysis was based on a complete three-dimensional burnup calculation. The code system was ELCOS, with BOXER as an assembly code for the generation of burnup-dependent cross sections and SILWER as a three-dimensional core code with thermal-hydraulic feedback.

Influence of film cooling hole angles and geometries on aerodynamic loss and net heat flux reduction

Turbine design engineers have to ensure that film cooling can provide sufficient protection to turbine blades from the hot mainstream gas, while keeping the losses low. Film cooling hole design parameters include inclination angle (a), compound angle (b), hole inlet geometry, and hole exit geometry. The influence of these parameters on aerodynamic loss and net heat flux reduction is investigated, with loss being the primary focus. Low-speed flat plate experiments have been conducted at momentum flux ratios of IR=0.16, 0.64, and 1.44. The film cooling aerodynamic mixing loss, generated by the mixing of mainstream and coolant, can be quantified using a three-dimensional analytical model that has been previously reported by the authors. The model suggests that for the same flow conditions, the aerodynamic mixing loss is the same for holes with different a and b but with the same angle between the mainstream and coolant flow directions (angle k). This relationship is assessed through experiments by testing two sets of cylindrical holes with different a and b: one set with k=35 deg, and another set with k=60 deg. The data confirm the stated relationship between α, β, k and the aerodynamic mixing loss. The results show that the designer should minimize k to obtain the lowest loss, but maximize b to achieve the best heat transfer performance. A suggestion on improving the loss model is also given. Five different hole geometries (α=35.0 deg, β=0 deg) were also tested: cylindrical hole, trenched hole, fan-shaped hole, D-Fan, and SD-Fan. The D-Fan and the SD-Fan have similar hole exits to the fan-shaped hole but their hole inlets are laterally expanded. The external mixing loss and the loss generated inside the hole are compared. It was found that the D-Fan and the SD-Fan have the lowest loss. This is attributed to their laterally expanded hole inlets, which lead to significant reduction in the loss generated inside the holes. As a result, the loss of these geometries is≈50% of the loss of the fan-shaped hole at IR=0.64 and 1.44. © 2013 by ASME.

Linear dimensionality reduction: Survey, insights, and generalizations

© 2015 John P. Cunningham and Zoubin Ghahramani. Linear dimensionality reduction methods are a cornerstone of analyzing high dimensional data, due to their simple geometric interpretations and typically attractive computational properties. These methods capture many data features of interest, such as covariance, dynamical structure, correlation between data sets, input-output relationships, and margin between data classes. Methods have been developed with a variety of names and motivations in many fields, and perhaps as a result the connections between all these methods have not been highlighted. Here we survey methods from this disparate literature as optimization programs over matrix manifolds. We discuss principal component analysis, factor analysis, linear multidimensional scaling, Fisher's linear discriminant analysis, canonical correlations analysis, maximum autocorrelation factors, slow feature analysis, sufficient dimensionality reduction, undercomplete independent component analysis, linear regression, distance metric learning, and more. This optimization framework gives insight to some rarely discussed shortcomings of well-known methods, such as the suboptimality of certain eigenvector solutions. Modern techniques for optimization over matrix manifolds enable a generic linear dimensionality reduction solver, which accepts as input data and an objective to be optimized, and returns, as output, an optimal low-dimensional projection of the data. This simple optimization framework further allows straightforward generalizations and novel variants of classical methods, which we demonstrate here by creating an orthogonal-projection canonical correlations analysis. More broadly, this survey and generic solver suggest that linear dimensionality reduction can move toward becoming a blackbox, objective-agnostic numerical technology.