H∞ sampled-data synthesis and related numerical issues

In this paper we present a new, compact derivation of state-space formulae for the so-called discretisation-based solution of the H∞ sampled-data control problem. Our approach is based on the established technique of continuous time-lifting, which is used to isometrically map the continuous-time, linear, periodically time-varying, sampled-data problem to a discretetime, linear, time-invariant problem. State-space formulae are derived for the equivalent, discrete-time problem by solving a set of two-point, boundary-value problems. The formulae accommodate a direct feed-through term from the disturbance inputs to the controlled outputs of the original plant and are simple, requiring the computation of only a single matrix exponential. It is also shown that the resultant formulae can be easily re-structured to give a numerically robust algorithm for computing the state-space matrices. © 1997 Elsevier Science Ltd. All rights reserved.

Numerical back-analysis of energy pile test at Lambeth College, London

Energy Piles present an efficient solution for long-term carbon emission reduction and sustainable construction. However, they have received only partial acceptance by the industry, because of concerns regarding the impact of cyclic thermal changes on the serviceability of energy pile foundations. This paper investigates the applicability of the hybrid load transfer approach to load-settlement analysis of single piles behavior during thermal energy exchange processes. Back-analysis results in terms of the thermal and mechanical response of energy piles show good agreement with field test results from Lambeth College in London. © ASCE 2011.

Application of the fixed point method for solution in time stepping finite element analysis using the inverse vector Jiles-Atherton model

An implementation of the inverse vector Jiles-Atherton model for the solution of non-linear hysteretic finite element problems is presented. The implementation applies the fixed point method with differential reluctivity values obtained from the Jiles-Atherton model. Differential reluctivities are usually computed using numerical differentiation, which is ill-posed and amplifies small perturbations causing large sudden increases or decreases of differential reluctivity values, which may cause numerical problems. A rule based algorithm for conditioning differential reluctivity values is presented. Unwanted perturbations on the computed differential reluctivity values are eliminated or reduced with the aim to guarantee convergence. Details of the algorithm are presented together with an evaluation of the algorithm by a numerical example. The algorithm is shown to guarantee convergence, although the rate of convergence depends on the choice of algorithm parameters. © 2011 IEEE.

Finite element solution of Navier Stokes equations adapted to a priori error estimates

The paper is based on qualitative properties of the solution of the Navier-Stokes equations for incompressible fluid, and on properties of their finite element solution. In problems with corner-like singularities (e.g. on the well-known L-shaped domain) usually some adaptive strategy is used. In this paper we present an alternative approach. For flow problems on domains with corner singularities we use the a priori error estimates and asymptotic expansion of the solution to derive an algorithm for refining the mesh near the corner. It gives very precise solution in a cheap way. We present some numerical results.

Boundary integral technique for the numerical modelling of the air flow for the Aaberg exhaust system

A boundary integral technique has been developed for the numerical simulation of the air flow for the Aaberg exhaust system. For the steady, ideal, irrotational air flow induced by a jet, the air velocity is an analytical function. The solution of the problem is formulated in the form of a boundary integral equation by seeking the solution of a mixed boundary-value problem of an analytical function based on the Riemann-Hilbert technique. The boundary integral equation is numerically solved by converting it into a system of linear algebraic equations, which are solved by the process of the Gaussian elimination. The air velocity vector at any point in the solution domain is then computed from the air velocity on the boundary of the solution domains.

Efficient numerical models for investigating the statistics of random dynamic systems

The study of random dynamic systems usually requires the definition of an ensemble of structures and the solution of the eigenproblem for each member of the ensemble. If the process is carried out using a conventional numerical approach, the computational cost becomes prohibitive for complex systems. In this work, an alternative numerical method is proposed. The results for the response statistics are compared with values obtained from a detailed stochastic FE analysis of plates. The proposed method seems to capture the statistical behaviour of the response with a reduced computational cost.

A numerical study into the effects of Bio-synthetic Paraffinic Kerosine blends with Jet-A fuel for civil aircraft engine

Growing concerns regarding fluctuating fuel costs and pollution targets for gas emissions, have led the aviation industry to seek alternative technologies to reduce its dependency on crude oil, and its net emissions. Recently blends of bio-fuel with kerosine, have become an alternative solution as they offer "greener" aircraft and reduce demand on crude oil. Interestingly, this technique is able to be implemented in current aircraft as it does not require any modification to the engine. Therefore, the present study investigates the effect of blends of bio-synthetic paraffinic kerosine with Jet-A in a civil aircraft engine, focusing on its performance and exhaust emissions. Two bio-fuels are considered: Jatropha Bio-synthetic Paraffinic Kerosine (JSPK) and Camelina Bio-synthetic Paraffinic Kerosine (CSPK); there are evaluated as pure fuels, and as 10% and 50% blend with Jet-A. Results obtained show improvement in thrust, fuel flow and SFC as composition of bio-fuel in the blend increases. At design point condition, results on engine emissions show reduction in NO x, and CO, but increases of CO is observed at fixed fuel condition, as the composition of bio-fuel in the mixture increases. Copyright © 2012 by ASME.

Numerical analysis of shallow foundations under combined horizontal and torsional loading

Foundations of subsea infrastructure in deep water subjected to asymmetric environmental loads have underscored the importance of combined torsional and horizontal loading effects on the bearing capacity of rectangular shallow foundations. The purpose of this study is to investigate the undrained sliding and torsional bearing capacity of rectangular and square shallow foundations together with the interaction response under combined loading using three-dimensional finite element (3D-FE) analysis. Upper bound plastic limit analysis is employed to establish a reference value for horizontal and torsional bearing capacity, and an interaction relationship for the combined loading condition. Satisfactory agreement of plastic limit analysis (PLA) and 3D-FE results for ultimate capacity and interaction curves ensures that simple PLA solution could be used to evaluate the bearing capacity problem of foundation under combined sliding and torsion.

Experimental and numerical study on AlGaInAs/AlGaAs distributed feedback lasers with GaInP gratings

Ridge-waveguide AlGaInAs/AlGaAs distributed feedback lasers with lattice-matched GaInP gratings were fabricated and their light-current characteristics, spectrum and far-field characteristics were measured. On the basis of our experimental results we analyze the effect of the electron stopper layer on light-current performance using the commercial laser simulation software PICS3D. The simulator is based on the self-consistent solution of drift diffusion equations, the Schrodinger equation, and the photon rate equation. The simulation results suggest that, with the use of a 80 nm-width p-doped Al0.6GaAs electron stopper layer, the slope efficiency can be increased and the threshold current can be reduced by more than 10 mA.

Experimental and numerical investigations on dissolution and recrystallization processes of GaSb/InSb/GaSb under microgravity and terrestrial conditions

The effects of gravity and crystal orientation on the dissolution of GaSb into InSb melt and the recrystallization of InGaSb were investigated under microgravity condition using a Chinese recoverable satellite and under normal gravity condition on earth. To investigate the effect of gravity on the solid/liquid interface and compositional profiles. a numerical simulation was carried out. The InSb crystal melted at 525 degrees C and then a part of GaSb dissolved into the InSb melt during heating to 706 degrees C and this process led to the formation of InGaSb solution. InGaSb solidified during the cooling process. The experimental and calculation results clearly show that the shape of the solid/liquid interface and compositional profiles in the solution were significantly affected by gravity. Under microgravity, as the Ga compositional profiles were uniform in the radial direction. the interfaces were almost parallel. On the contrary, for normal gravity condition, as large amounts of Ga moved up in the upper region due to buoyancy, the dissolved zone broadened towards gravitational direction. Also. during the cooling process, needle crystals of InGaSb started appearing and the value of x of InxGa1-xSb crystals increased with the decrease of temperature. The GaSb with the (111)B plane dissolved into the InSb melt much more than that of the (111)A plane. (C) 2000 Elsevier Science B.V. All rights reserved.

Two-dimensional numerical simulation of the channel electron in an In0.52Al0.48As/In0.53Ga0.48As HEMT

A two-dimensional quantum model based on the solution of Schrodinger and Poisson equations is first presented for In0.52Al0.48As/In0.53Ga0.47As/InP HEMT. According to the model, the two-dimensional distributions of electron density and transverse electric field in the channel of InAlAs/InGaAs HEMT are discussed.

Study on the origin of inverted phase in drying solution-cast block copolymer films

In a previous study, we reported observation of the novel inverted phase (the minority blocks comprising the continuum phase) in kinetically controlled phase separating solution-cast poly(styrene-b-butadiene-b-styrene) (SBS) triblock copolymer films [Zhang et al. Macromolecules 2000, 33, 9561-7]. In this study, we adopt the same approach to investigate the formation of inverted phase in a series of solution-cast poly(styrene-b-butadiene) (SB) asymmetric diblock copolymers having nearly equal polystyrene (PS) weight fraction (about 30 wt %) but different molecular weights. The microstructure of the solution-cast block copolymer films resulting from different solvent evaporation rates, R, was inspected, from which the kinetically frozen-in phase structures at qualitatively different block copolymer concentrations and correspondingly different effective interaction parameter, chieff, can be deduced. Our result shows that there is a threshold molecular weight or range of molecular weight below which the unusual inverted phase is accessible by controlling the solvent evaporation rate. In comparing the present result with that of our previous study on the SBS triblock copolymer, we find that the formation of the inverted phase has little bearing on the chain architecture. We performed numerical calculations for the free energy of block copolymer cylinders and found that the normal phase is always preferred irrespective of the interaction parameter and molecular weight, which suggests the formation of the inverted phase to have a kinetic origin.

The Image Irradiance Equation: Its Solution and Application

How much information about the shape of an object can be inferred from its image? In particular, can the shape of an object be reconstructed by measuring the light it reflects from points on its surface? These questions were raised by Horn [HO70] who formulated a set of conditions such that the image formation can be described in terms of a first order partial differential equation, the image irradiance equation. In general, an image irradiance equation has infinitely many solutions. Thus constraints necessary to find a unique solution need to be identified. First we study the continuous image irradiance equation. It is demonstrated when and how the knowledge of the position of edges on a surface can be used to reconstruct the surface. Furthermore we show how much about the shape of a surface can be deduced from so called singular points. At these points the surface orientation is uniquely determined by the measured brightness. Then we investigate images in which certain types of silhouettes, which we call b-silhouettes, can be detected. In particular we answer the following question in the affirmative: Is there a set of constraints which assure that if an image irradiance equation has a solution, it is unique? To this end we postulate three constraints upon the image irradiance equation and prove that they are sufficient to uniquely reconstruct the surface from its image. Furthermore it is shown that any two of these constraints are insufficient to assure a unique solution to an image irradiance equation. Examples are given which illustrate the different issues. Finally, an overview of known numerical methods for computing solutions to an image irradiance equation are presented.

Asymptotic behaviour of the elastic solution near the tip of a crack situated at a nonideal interface

Mishuris, G; Kuhn, G., (2001) 'Asymptotic behaviour of the elastic solution near the tip of a crack situated at a nonideal interface', Zeitschrift f?r Angewandte Mathematik und Mechanik 81(12) pp.811-826 RAE2008