A critical assessment of the finite element method for calculating electric and magnetic fields

The present dissertation is concerned with the determination of the magnetic field distribution in ma[.rnetic electron lenses by means of the finite element method. In the differential form of this method a Poisson type equation is solved by numerical methods over a finite boundary. Previous methods of adapting this procedure to the requirements of digital computers have restricted its use to computers of extremely large core size. It is shown that by reformulating the boundary conditions, a considerable reduction in core store can be achieved for a given accuracy of field distribution. The magnetic field distribution of a lens may also be calculated by the integral form of the finite element rnethod. This eliminates boundary problems mentioned but introduces other difficulties. After a careful analysis of both methods it has proved possible to combine the advantages of both in a .new approach to the problem which may be called the 'differential-integral' finite element method. The application of this method to the determination of the magnetic field distribution of some new types of magnetic lenses is described. In the course of the work considerable re-programming of standard programs was necessary in order to reduce the core store requirements to a minimum.

Dynamic simulation of chemical plant

A CSSL- type modular FORTRAN package, called ACES, has been developed to assist in the simulation of the dynamic behaviour of chemical plant. ACES can be harnessed, for instance, to simulate the transients in startups or after a throughput change. ACES has benefited from two existing simulators. The structure was adapted from ICL SLAM and most plant models originate in DYFLO. The latter employs sequential modularisation which is not always applicable to chemical engineering problems. A novel device of twice- round execution enables ACES to achieve general simultaneous modularisation. During the FIRST ROUND, STATE-VARIABLES are retrieved from the integrator and local calculations performed. During the SECOND ROUND, fresh derivatives are estimated and stored for simultaneous integration. ACES further includes a version of DIFSUB, a variable-step integrator capable of handling stiff differential systems. ACES is highly formalised . It does not use pseudo steady- state approximations and excludes inconsistent and arbitrary features of DYFLO. Built- in debug traps make ACES robust. ACES shows generality, flexibility, versatility and portability, and is very convenient to use. It undertakes substantial housekeeping behind the scenes and thus minimises the detailed involvement of the user. ACES provides a working set of defaults for simulation to proceed as far as possible. Built- in interfaces allow for reactions and user supplied algorithms to be incorporated . New plant models can be easily appended. Boundary- value problems and optimisation may be tackled using the RERUN feature. ACES is file oriented; a STATE can be saved in a readable form and reactivated later. Thus piecewise simulation is possible. ACES has been illustrated and verified to a large extent using some literature-based examples. Actual plant tests are desirable however to complete the verification of the library. Interaction and graphics are recommended for future work.

Developments and extensions of Roth's method of double Fourier series with applications to the solution of electromagnetic field problems

The first part of the thesis compares Roth's method with other methods, in particular the method of separation of variables and the finite cosine transform method, for solving certain elliptic partial differential equations arising in practice. In particular we consider the solution of steady state problems associated with insulated conductors in rectangular slots. Roth's method has two main disadvantages namely the slow rate of con­vergence of the double Fourier series and the restrictive form of the allowable boundary conditions. A combined Roth-separation of variables method is derived to remove the restrictions on the form of the boundary conditions and various Chebyshev approximations are used to try to improve the rate of convergence of the series. All the techniques are then applied to the Neumann problem arising from balanced rectangular windings in a transformer window. Roth's method is then extended to deal with problems other than those resulting from static fields. First we consider a rectangular insulated conductor in a rectangular slot when the current is varying sinusoidally with time. An approximate method is also developed and compared with the exact method.The approximation is then used to consider the problem of an insulated conductor in a slot facing an air gap. We also consider the exact method applied to the determination of the eddy-current loss produced in an isolated rectangular conductor by a transverse magnetic field varying sinusoidally with time. The results obtained using Roth's method are critically compared with those obtained by other authors using different methods. The final part of the thesis investigates further the application of Chebyshdev methods to the solution of elliptic partial differential equations; an area where Chebyshev approximations have rarely been used. A poisson equation with a polynomial term is treated first followed by a slot problem in cylindrical geometry.

Eye movement problems and differential diagnosis of the parkinsonian syndromes

Differential clinical diagnosis of the parkinsonian syndromes, viz., Parkinson’s disease (PD), progressive supranuclear palsy (PSP), dementia with Lewy bodies (DLB), and multiple system atrophy (MSA) can be difficult. Eye movement problems, however, are a chronic complication of many of these disorders and may be a useful aid to diagnosis. Hence, the presence in PSP of vertical supranuclear gaze palsy, fixation instability, lid retraction, blepharospasm, and apraxia of eyelid opening and closing is useful in separating PD from PSP. Moreover, atypical features of PSP include slowing of upward saccades, moderate slowing of downward saccades, the presence of a full range of voluntary vertical eye movements, a curved trajectory of oblique saccades, and absence of square-wave jerks. Downgaze palsy is probably the most useful diagnostic clinical symptom of PSP. By contrast, DLB patients are specifically impaired in both reflexive and saccadic execution and in the performance of more complex saccadic eye movement tasks. Problems in convergence in DLB are also followed by akinesia and rigidity. Abnormal ocular fixation may occur in a significant proportion of MSA patients along with excessive square-wave jerks, a mild supranuclear gaze palsy, a gaze-evoked nystagmus, a positioning down-beat nystagmus, mild-moderate saccadic hypometria, impaired smooth pursuit movements, and reduced vestibulo-ocular reflex (VOR) suppression. There may be considerable overlap between the eye movement problems characteristic of the various parkinsonian disorders, but taken together with other signs and symptoms, can be a useful aid in differential diagnosis, especially in the separation of PD and PSP.

Quasi-separation of the biharmonic partial differential equation

In this paper, we consider analytical and numerical solutions to the Dirichlet boundary-value problem for the biharmonic partial differential equation on a disc of finite radius in the plane. The physical interpretation of these solutions is that of the harmonic oscillations of a thin, clamped plate. For the linear, fourth-order, biharmonic partial differential equation in the plane, it is well known that the solution method of separation in polar coordinates is not possible, in general. However, in this paper, for circular domains in the plane, it is shown that a method, here called quasi-separation of variables, does lead to solutions of the partial differential equation. These solutions are products of solutions of two ordinary linear differential equations: a fourth-order radial equation and a second-order angular differential equation. To be expected, without complete separation of the polar variables, there is some restriction on the range of these solutions in comparison with the corresponding separated solutions of the second-order harmonic differential equation in the plane. Notwithstanding these restrictions, the quasi-separation method leads to solutions of the Dirichlet boundary-value problem on a disc with centre at the origin, with boundary conditions determined by the solution and its inward drawn normal taking the value 0 on the edge of the disc. One significant feature for these biharmonic boundary-value problems, in general, follows from the form of the biharmonic differential expression when represented in polar coordinates. In this form, the differential expression has a singularity at the origin, in the radial variable. This singularity translates to a singularity at the origin of the fourth-order radial separated equation; this singularity necessitates the application of a third boundary condition in order to determine a self-adjoint solution to the Dirichlet boundary-value problem. The penultimate section of the paper reports on numerical solutions to the Dirichlet boundary-value problem; these results are also presented graphically. Two specific cases are studied in detail and numerical values of the eigenvalues are compared with the results obtained in earlier studies.

An iterative method based on boundary integrals for elliptic Cauchy problems in semi-infinite domains

In this study, we investigate the problem of reconstruction of a stationary temperature field from given temperature and heat flux on a part of the boundary of a semi-infinite region containing an inclusion. This situation can be modelled as a Cauchy problem for the Laplace operator and it is an ill-posed problem in the sense of Hadamard. We propose and investigate a Landweber-Fridman type iterative method, which preserve the (stationary) heat operator, for the stable reconstruction of the temperature field on the boundary of the inclusion. In each iteration step, mixed boundary value problems for the Laplace operator are solved in the semi-infinite region. Well-posedness of these problems is investigated and convergence of the procedures is discussed. For the numerical implementation of these mixed problems an efficient boundary integral method is proposed which is based on the indirect variant of the boundary integral approach. Using this approach the mixed problems are reduced to integral equations over the (bounded) boundary of the inclusion. Numerical examples are included showing that stable and accurate reconstructions of the temperature field on the boundary of the inclusion can be obtained also in the case of noisy data. These results are compared with those obtained with the alternating iterative method.

Numerical solution of parabolic Cauchy problems in planar corner domains

An iterative method for the parabolic Cauchy problem in planar domains having a finite number of corners is implemented based on boundary integral equations. At each iteration, mixed well-posed problems are solved for the same parabolic operator. The presence of corner points renders singularities of the solutions to these mixed problems, and this is handled with the use of weight functions together with, in the numerical implementation, mesh grading near the corners. The mixed problems are reformulated in terms of boundary integrals obtained via discretization of the time-derivative to obtain an elliptic system of partial differential equations. To numerically solve these integral equations a Nyström method with super-algebraic convergence order is employed. Numerical results are presented showing the feasibility of the proposed approach. © 2014 IMACS.

Trade-offs in customer–employee focus:implications on boundary-spanning capabilities

In pursuing their particular strategic goals, organisations place differential emphasis in their activities on employees and customers. Those who focus predominantly on employees see them as key resources contributing to successful achievement of goals, while those who focus on customers believe that customer commitment is the prime source of financial prosperity and competitive advantage. In all firms, there is potential for conflict and subsequent trade-offs in the emphasis placed on employee resources and customers. This paper investigates the implications of these potential trade-offs on firms’ boundary-spanning capabilities (those capabilities that enable firms to interact effectively with their customers, their suppliers, their strategic allies and others in their business networks) in European and Australasian markets. The results suggest that those companies which achieve a strong but balanced focus on both employees and customers have heightened boundary-spanning capabilities. The comparison between UK and Australian businesses indicates that the balance is more highly developed in Australasian than European (UK) markets.

Trade-offs in customer-employee focus:Implications on boundary-spanning capabilities

In pursuing their particular strategic goals, organisations place differential emphasis in their activities on employees and customers. Those who focus predominantly on employees see them as key resources contributing to successful achievement of goals, while those who focus on customers believe that customer commitment is the prime source of financial prosperity and competitive advantage. In all firms, there is potential for conflict and subsequent trade-offs in the emphasis placed on employee resources and customers. This paper investigates the implications of these potential trade-offs on firms' boundary-spanning capabilities (those capabilities that enable firms to interact effectively with their customers, their suppliers, their strategic allies and others in their business networks) in European and Australasian markets. The results suggest that those companies which achieve a strong but balanced focus on both employees and customers have heightened boundary-spanning capabilities. The comparison between UK and Australian businesses indicates that the balance is more highly developed in Australasian than European (UK) markets. © 2004 by The Haworth Press, Inc. All rights reserved.

Gaussian process approximations of stochastic differential equations

Stochastic differential equations arise naturally in a range of contexts, from financial to environmental modeling. Current solution methods are limited in their representation of the posterior process in the presence of data. In this work, we present a novel Gaussian process approximation to the posterior measure over paths for a general class of stochastic differential equations in the presence of observations. The method is applied to two simple problems: the Ornstein-Uhlenbeck process, of which the exact solution is known and can be compared to, and the double-well system, for which standard approaches such as the ensemble Kalman smoother fail to provide a satisfactory result. Experiments show that our variational approximation is viable and that the results are very promising as the variational approximate solution outperforms standard Gaussian process regression for non-Gaussian Markov processes.

Differential perception of innovation with particular reference to the organisational context

Differential perception of innovation is a research area which has been advocated as a suitable topic for study in recent years. It developed from the problems encountered within earlier perception of innovation studies which sought to establish what characteristics of an innovation affected the ease of its adoption. While some success was achieved In relating perception of innovation to adoption behaviour, variability encountered Within groups expected - to fercelve innovation similarly suggested that the needs and experiences of the potential adopter were significantly affecting the research findings. Such analysis being supported by both sociological and psychological perceptual research. The present study sought to identify the presence of differential perception of innovation and explore the nature of the process. It was decided to base the research in an organisational context and to concentrate upon manufacturing innovation. It has been recognised that such adoption of technological innovation is commonly the product of a collective decision-making process, involving individuals from a variety of occupational backgrounds, both in terms of occupational speciality and level within the hierarchy. Such roles appeared likely to significantly influence perception of technological innovation, as gathered through an appropriate measure and were readily identifiable. Data vas collected by means of a face-to-face card presentation technique, a questionnaire and through case study material. Differential perception of innovation effects were apparent In the results, many similarities and differences of perception being related to the needs and experiences of the individuals studied. Phenomenological analysis, which recognises the total nature of experience in infiuencing behaviour, offered the best means of explaining the findings. It was also clear that the bureaucratic model of role definition was not applicable to the area studied, it seeming likely that such definitions are weaker under conditions of uncertainty, such as encountered in innovative decision-making.

Stock selection with full-scale optimization and differential evolution

When composing stock portfolios, managers frequently choose among hundreds of stocks. The stocks' risk properties are analyzed with statistical tools, and managers try to combine these to meet the investors' risk profiles. A recently developed tool for performing such optimization is called full-scale optimization (FSO). This methodology is very flexible for investor preferences, but because of computational limitations it has until now been infeasible to use when many stocks are considered. We apply the artificial intelligence technique of differential evolution to solve FSO-type stock selection problems of 97 assets. Differential evolution finds the optimal solutions by self-learning from randomly drawn candidate solutions. We show that this search technique makes large scale problem computationally feasible and that the solutions retrieved are stable. The study also gives further merit to the FSO technique, as it shows that the solutions suit investor risk profiles better than portfolios retrieved from traditional methods.

Ingroup and outgroup minorities:differential impact upon public and private responses

The experiment which is presented in this paper was designed to overcome some of the problems associated with previous research investigating the effects of social categorization and minority influence. Sixty-eight fourteen-year-old British Secondary School pupils indicated their attitudes towards a 'grant for pupils' before and after reading a text which advocated a minority position. The text was attributed as being the work of either pupils from their own school (ingroup minority) or from a school they discriminated against (outgroup minority). Responses were either made in ‘public’ (by telling subjects that other pupils would see their responses) or in ‘private’ (by subjects putting their responses into a ‘ballot box’). The results showed that on public responses ingroup minorities had more influence than outgroup minorities while there was no difference on private responses. Also, greater change occurred when responses were made in private than in public. These results are compatible with the intergroup analysis of minority influence.

An alternating boundary integral based method for a Cauchy problem for the Laplace equation in a quadrant

We study the Cauchy problem for the Laplace equation in a quadrant (quarter-plane) containing a bounded inclusion. Given the values of the solution and its derivative on the edges of the quadrant the solution is reconstructed on the boundary of the inclusion. This is achieved using an alternating iterative method where at each iteration step mixed boundary value problems are being solved. A numerical method is also proposed and investigated for the direct mixed problems reducing these to integral equations over the inclusion. Numerical examples verify the efficiency of the proposed scheme.

A variational method and approximations of a Cauchy problem for elliptic equations

A Cauchy problem for general elliptic second-order linear partial differential equations in which the Dirichlet data in H½(?1 ? ?3) is assumed available on a larger part of the boundary ? of the bounded domain O than the boundary portion ?1 on which the Neumann data is prescribed, is investigated using a conjugate gradient method. We obtain an approximation to the solution of the Cauchy problem by minimizing a certain discrete functional and interpolating using the finite diference or boundary element method. The minimization involves solving equations obtained by discretising mixed boundary value problems for the same operator and its adjoint. It is proved that the solution of the discretised optimization problem converges to the continuous one, as the mesh size tends to zero. Numerical results are presented and discussed.