An integral equation method for a mixed initial boundary value problem for unsteady Stokes system in a doubly-connected domain

We present a novel numerical method for a mixed initial boundary value problem for the unsteady Stokes system in a planar doubly-connected domain. Using a Laguerre transformation the unsteady problem is reduced to a system of boundary value problems for the Stokes resolvent equations. Employing a modied potential approach we obtain a system of boundary integral equations with various singularities and we use a trigonometric quadrature method for their numerical solution. Numerical examples are presented showing that accurate approximations can be obtained with low computational cost.

Numerical solution of the Dirichlet initial boundary value problem for the heat equation in exterior 3-dimensional domains using integral equations

A numerical method for the Dirichlet initial boundary value problem for the heat equation in the exterior and unbounded region of a smooth closed simply connected 3-dimensional domain is proposed and investigated. This method is based on a combination of a Laguerre transformation with respect to the time variable and an integral equation approach in the spatial variables. Using the Laguerre transformation in time reduces the parabolic problem to a sequence of stationary elliptic problems which are solved by a boundary layer approach giving a sequence of boundary integral equations of the first kind to solve. Under the assumption that the boundary surface of the solution domain has a one-to-one mapping onto the unit sphere, these integral equations are transformed and rewritten over this sphere. The numerical discretisation and solution are obtained by a discrete projection method involving spherical harmonic functions. Numerical results are included.

An iterative method for a Cauchy problem for the heat equation

An iterative method for reconstruction of the solution to a parabolic initial boundary value problem of second order from Cauchy data is presented. The data are given on a part of the boundary. At each iteration step, a series of well-posed mixed boundary value problems are solved for the parabolic operator and its adjoint. The convergence proof of this method in a weighted L2-space is included.

Optimization of chemical plant simulation using double collocation

A method has been constructed for the solution of a wide range of chemical plant simulation models including differential equations and optimization. Double orthogonal collocation on finite elements is applied to convert the model into an NLP problem that is solved either by the VF 13AD package based on successive quadratic programming, or by the GRG2 package, based on the generalized reduced gradient method. This approach is termed simultaneous optimization and solution strategy. The objective functional can contain integral terms. The state and control variables can have time delays. Equalities and inequalities containing state and control variables can be included into the model as well as algebraic equations and inequalities. The maximum number of independent variables is 2. Problems containing 3 independent variables can be transformed into problems having 2 independent variables using finite differencing. The maximum number of NLP variables and constraints is 1500. The method is also suitable for solving ordinary and partial differential equations. The state functions are approximated by a linear combination of Lagrange interpolation polynomials. The control function can either be approximated by a linear combination of Lagrange interpolation polynomials or by a piecewise constant function over finite elements. The number of internal collocation points can vary by finite elements. The residual error is evaluated at arbitrarily chosen equidistant grid-points, thus enabling the user to check the accuracy of the solution between collocation points, where the solution is exact. The solution functions can be tabulated. There is an option to use control vector parameterization to solve optimization problems containing initial value ordinary differential equations. When there are many differential equations or the upper integration limit should be selected optimally then this approach should be used. The portability of the package has been addressed converting the package from V AX FORTRAN 77 into IBM PC FORTRAN 77 and into SUN SPARC 2000 FORTRAN 77. Computer runs have shown that the method can reproduce optimization problems published in the literature. The GRG2 and the VF I 3AD packages, integrated into the optimization package, proved to be robust and reliable. The package contains an executive module, a module performing control vector parameterization and 2 nonlinear problem solver modules, GRG2 and VF I 3AD. There is a stand-alone module that converts the differential-algebraic optimization problem into a nonlinear programming problem.

Determining the temperature from incomplete boundary data

An iterative procedure for determining temperature fields from Cauchy data given on a part of the boundary is presented. At each iteration step, a series of mixed well-posed boundary value problems are solved for the heat operator and its adjoint. A convergence proof of this method in a weighted L2-space is included, as well as a stopping criteria for the case of noisy data. Moreover, a solvability result in a weighted Sobolev space for a parabolic initial boundary value problem of second order with mixed boundary conditions is presented. Regularity of the solution is proved. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

Time complexity and convergence analysis of domain theoretic Picard method

We present an implementation of the domain-theoretic Picard method for solving initial value problems (IVPs) introduced by Edalat and Pattinson [1]. Compared to Edalat and Pattinson's implementation, our algorithm uses a more efficient arithmetic based on an arbitrary precision floating-point library. Despite the additional overestimations due to floating-point rounding, we obtain a similar bound on the convergence rate of the produced approximations. Moreover, our convergence analysis is detailed enough to allow a static optimisation in the growth of the precision used in successive Picard iterations. Such optimisation greatly improves the efficiency of the solving process. Although a similar optimisation could be performed dynamically without our analysis, a static one gives us a significant advantage: we are able to predict the time it will take the solver to obtain an approximation of a certain (arbitrarily high) quality.

Block co-polymerization by transformation reactions

The aim of this study was to use the transformation of anionic to metathesis polymerization to produce block co-polymers of styrene-b-pentenylene using WC16 /PStLi and WC16/PStLi/ AlEtC12 catalyst systems. Analysis of the products using SEC and 1H and 13C NMR spectroscopy enabled mechanisms for metathesis initiation reactions to be proposed. The initial work involved preparation of the constituent homo-polymers. Solutions of polystyryllithium in cyclohexane were prepared and diluted so that the [PStLi]o<2x10-3M. The dilution produced initial rapid decay of the active species, followed by slower spontaneous decay within a period of days. This was investigated using UV / visible spectrophotometry and the wavelength of maximum absorbance of the PStLi was found to change with the decay from an initial value of 328mn. to λmax of approximately 340nm. after 4-7 days. SEC analysis of solutions of polystyrene, using RI and UV / visible (set at 254nm.) detectors; showed the UV:RI peak area was constant for a range of polystyrene samples of different moleculor weight. Samples of polypentenylene were prepared and analysed using SEC. Unexpectedly the solutions showed an absorbance at 254nm. which had to be considered when this technique was used subsequently to analyse polymer samples to determine their styrene/ pentenylene co-polymer composition. Cyclohexane was found to be a poor solvent for these ring-opening metathesis polymerizations of cyclopentene. Attempts to produce styrene-b-pentenylene block co-polymers, using a range of co-catalyst systems, were generally unsuccessful as the products were shown to be mainly homopolymers. The character of the polymers did suggest that several catalytic species are present in these systems and mechanisms have been suggested for the formation of initiating carbenes. Evidence of some low molecular weight product with co-polymer character has been obtained. Further investigation indicated that this is most likely to be ABA block copolymer, which led to a mechanism being proposed for the termination of the polymerization.

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.

Stereoacuity as an indicator of prism adaptation

Purpose: The purpose of this study was to determine whether stereoacuity can be used as an indicator of prism adaptation. In particular, we wanted to know whether the time required for stereoacuity to return to the initial level after viewing through a prism can be used to determine the degree of adaptation. Materials and Methods: Eighteen subjects participated in this study. Stereoacuity and dissociated phoria were determined using the TNO stereotest and the Maddox rod, respectively. Prism vergences were measured using a prism bar. For each participant, prism power equivalent to the blur point of base-in (BI) and base-out (BO) fusional vergence at 40 cm was divided and placed in front of both eyes. At 0, 3, 6, 9 and 12 min after prism introduction, the stereoacuity was measured, and at 0 and 12 min, the heterophoria was measured. Results: The repeated measures ANOVA showed a significant difference between the mean stereoacuity for BI and BO prisms at the different measurement times (p < 0.05). For BO prism, the initial value was different between 0 and 3 min after the prism introduction, whereas for BI prism, a difference in stereoacuity was found between the pre-prism value and the value at 0, 3 and 6 min. The size of the heterophoria with BO and BI prisms was different from 0 to 12 min (p < 0.05). Conclusion: The time required for stereoacuity to return to baseline level was more than 3 min for BO, and more than 6 min for BI prism. In addition, the time required to return to baseline values was not similar for the stereoacuity and heterophoria. The recovery of stereoacuity is slower when adapting to divergence, as when looking from near to far. This implies that stereopsis responds faster to near targets than to distant one, and may precede complete phoria adaptation. © 2014 Informa Healthcare USA, Inc. All rights reserved: reproduction in whole or part not permitted.

An aspect oriented performance analysis environment

AOSD'03 Practitioner Report Performance analysis is motivated as an ideal domain for benefiting from the application of Aspect Oriented (AO) technology. The experience of a ten week project to apply AO to the performance analysis domain is described. We show how all phases of a performance analysts’ activities – initial profiling, problem identification, problem analysis and solution exploration – were candidates for AO technology assistance – some being addressed with more success than others. A Profiling Workbench is described that leverages the capabilities of AspectJ, and delivers unique capabilities into the hands of developers exploring caching opportunities.

New insights in the deactivation of sulfonic modified SBA-15 catalysts for biodiesel production from low-grade oleaginous feedstock

Arenesulfonic-acid functionalized SBA-15 materials have been used in the production of biodiesel from low grade oleaginous feedstock. These materials display an outstanding catalytic activity, being able to promote the transformation of crude palm oil with methanol into fatty acid methyl esters with high yield (85%) under mild reaction conditions. However, high sensitivity of the catalyst against poisoning by different substances has also been detected. Thus, alkaline metal cations, such as sodium or potassium exert a negative influence on the catalytic activity of these materials, being necessary amounts around 500 ppm of sodium in the reaction media to decrease the catalytic activity of these materials to a half of its initial value in just two reaction runs. The deactivation of arenesulfonic acid functionalized SBA-15 materials seems to occur in this case by ion exchange of the acid protons at the sulfonic groups. Organic unsaponifiable compounds like lecithin or retinol also induce a negative influence in the catalytic activity of these sulfonic acid-based materials, though not so intense as in the case of alkaline metals. The deactivating mechanism associated to the influence of the organic compounds seems to be linked to the adsorption of such substances onto the catalytic acid sites as well as on the silica surface. The accumulation of lecithin in the surface of catalyst, observed by means of thermogravimetric analysis, suggest the creation of a strong interaction, probably by ion pair, between this compound and the sulfonic acid group.

Impact of the carrier relaxation paths on two-state operation in quantum dot lasers

We study InGaAs QD laser operating simultaneously at ground (GS) and excited (ES) states under 30ns pulsed-pumping and distinguish three regimes of operation depending on the pump current and the carrier relaxation pathways. An increased current leads to an increase in ES intensity and to a decrease in GS intensity (or saturation) for low pump range, as typical for the cascade-like pathway. Both the GS and ES intensities are steadily increased for high current ranges, which prove the dominance of the direct capture pathway. The relaxation oscillations are not pronounced for these ranges. For the mediate currents, the interplay between the both pathways leads to the damped large amplitude relaxation oscillations with significant deviation of the relaxation oscillation frequency from the initial value during the pulse.

Quid pro quo:reflections on the value of problem structuring workshops

Attracting clients who are willing to invest in using a problem structuring method (PSM) can be particularly difficult for the emerging generation of modellers. There are many reasons for this, not least that the benefits of a problem structuring intervention are vague and evidence of benefits are often anecdotal for example, claims of constructing a deeper understanding of the problem or building the commitment of a group to implementing an outcome. This paper contributes to the evaluation of problem structuring methods by reflecting on the quid pro quo that a client and problem structuring modeller can enjoy from collaboration. The paper reflects on 21 cases, where Journey Making (a problem structuring method) was used with 16 organizations to help managers agree a suite of actions to tackle a complex strategic issue. The reflections are clustered around those benefits that pertain to: PSMs in general; PSMs that use computer-supported workshops; the Journey Making methodology.

Premium e-grocery:exploring value in logistics integrated service solutions

E-grocery is gradually becoming viable or a necessity for many families. Yet, most e-supermarkets are seen as providers of low value "staple" and bulky goods mainly. While each store has a large number of SKU available, these products are mainly necessity goods with low marginal value for hedonistic consumption. A need to acquire diverse products (e.g., organic), premium priced products (e.g., wine) for special occasions (e.g., anniversary, birthday), or products just for health related reasons (e.g., allergies, diabetes) are yet to be served via one-stop e-tailers. In this paper, we design a mathematical model that takes into account consumers' geo-demographics and multi-product sourcing capacity for creating critical mass and profit. Our mathematical model is a variant of Capacitated Vehicle Routing Problem with Time Windows (CVRPTW), which we extend by adding intermediate locations for trucks to meet and exchange goods. We illustrate our model for the city of Istanbul using GIS maps, and discuss its various extensions as well as managerial implications.

A hybrid genetic algorithm for the multi-depot vehicle routing problem

The distribution of finished products from depots to customers is a practical and challenging problem in logistics management. Better routing and scheduling decisions can result in higher level of customer satisfaction because more customers can be served in a shorter time. The distribution problem is generally formulated as the vehicle routing problem (VRP). Nevertheless, there is a rigid assumption that there is only one depot. In cases, for instance, where a logistics company has more than one depot, the VRP is not suitable. To resolve this limitation, this paper focuses on the VRP with multiple depots, or multi-depot VRP (MDVRP). The MDVRP is NP-hard, which means that an efficient algorithm for solving the problem to optimality is unavailable. To deal with the problem efficiently, two hybrid genetic algorithms (HGAs) are developed in this paper. The major difference between the HGAs is that the initial solutions are generated randomly in HGA1. The Clarke and Wright saving method and the nearest neighbor heuristic are incorporated into HGA2 for the initialization procedure. A computational study is carried out to compare the algorithms with different problem sizes. It is proved that the performance of HGA2 is superior to that of HGA1 in terms of the total delivery time.