Exact Solution of a One-Dimensional Multicomponent Lattice Gas with Hyperbolic Interaction

We present the exact solution to a one-dimensional multicomponent quantum lattice model interacting by an exchange operator which falls off as the inverse sinh square of the distance. This interaction contains a variable range as a parameter and can thus interpolate between the known solutions for the nearest-neighbor chain and the inverse-square chain. The energy, susceptibility, charge stiffness, and the dispersion relations for low-lying excitations are explicitly calculated for the absolute ground state, as a function of both the range of the interaction and the number of species of fermions.

Improving the convergence rate of the inverse iteration method

Solution of generalized eigenproblem, K phi = lambda M phi, by the classical inverse iteration method exhibits slow convergence for some eigenproblems. In this paper, a modified inverse iteration algorithm is presented for improving the convergence rate. At every iteration, an optimal linear combination of the latest and the preceding iteration vectors is used as the input vector for the next iteration. The effectiveness of the proposed algorithm is demonstrated for three typical eigenproblems, i.e. eigenproblems with distinct, close and repeated eigenvalues. The algorithm yields 29, 96 and 23% savings in computational time, respectively, for these problems. The algorithm is simple and easy to implement, and this renders the algorithm even more attractive.

Swam Intelligent Based Inverse Model for Laboratory Double Reservoir Diffusion Experiments

Theoretical approaches are of fundamental importance to predict the potential impact of waste disposal facilities on ground water contamination. Appropriate design parameters are, in general, estimated by fitting the theoretical models to a field monitoring or laboratory experimental data. Double-reservoir diffusion (Transient Through-Diffusion) experiments are generally conducted in the laboratory to estimate the mass transport parameters of the proposed barrier material. These design parameters are estimated by manual parameter adjusting techniques (also called eye-fitting) like Pollute. In this work an automated inverse model is developed to estimate the mass transport parameters from transient through-diffusion experimental data. The proposed inverse model uses particle swarm optimization (PSO) algorithm which is based on the social behaviour of animals for finding their food sources. Finite difference numerical solution of the transient through-diffusion mathematical model is integrated with the PSO algorithm to solve the inverse problem of parameter estimation.The working principle of the new solver is demonstrated by estimating mass transport parameters from the published transient through-diffusion experimental data. The estimated values are compared with the values obtained by existing procedure. The present technique is robust and efficient. The mass transport parameters are obtained with a very good precision in less time

Robust solver based on modified particle swarm optimization for improved solution of diffusion transport through containment facilities

Accurate estimation of mass transport parameters is necessary for overall design and evaluation processes of the waste disposal facilities. The mass transport parameters, such as effective diffusion coefficient, retardation factor and diffusion accessible porosity, are estimated from observed diffusion data by inverse analysis. Recently, particle swarm optimization (PSO) algorithm has been used to develop inverse model for estimating these parameters that alleviated existing limitations in the inverse analysis. However, PSO solver yields different solutions in successive runs because of the stochastic nature of the algorithm and also because of the presence of multiple optimum solutions. Thus the estimated mean solution from independent runs is significantly different from the best solution. In this paper, two variants of the PSO algorithms are proposed to improve the performance of the inverse analysis. The proposed algorithms use perturbation equation for the gbest particle to gain information around gbest region on the search space and catfish particles in alternative iterations to improve exploration capabilities. Performance comparison of developed solvers on synthetic test data for two different diffusion problems reveals that one of the proposed solvers, CPPSO, significantly improves overall performance with improved best, worst and mean fitness values. The developed solver is further used to estimate transport parameters from 12 sets of experimentally observed diffusion data obtained from three diffusion problems and compared with published values from the literature. The proposed solver is quick, simple and robust on different diffusion problems. (C) 2012 Elsevier Ltd. All rights reserved.

A note on the inverse mode shape problem for bars, beams, and plates

Motivated by the idea of designing a structure for a desired mode shape, intended towards applications such as resonant sensors, actuators and vibration confinement, we present the inverse mode shape problem for bars, beams and plates in this work. The objective is to determine the cross-sectional profile of these structures, given a mode shape, boundary condition and the mass. The contribution of this article is twofold: (i) A numerical method to solve this problem when a valid mode shape is provided in the finite element framework for both linear and nonlinear versions of the problem. (ii) An analytical result to prove the uniqueness and existence of the solution in the case of bars. This article also highlights a very important question of the validity of a mode shape for any structure of given boundary conditions.

A BIEM optimization method for fracture dynamics inverse problem

In the present paper, based on the theory of dynamic boundary integral equation, an optimization method for crack identification is set up in the Laplace frequency space, where the direct problem is solved by the author's new type boundary integral equations and a method for choosing the high sensitive frequency region is proposed. The results show that the method proposed is successful in using the information of boundary elastic wave and overcoming the ill-posed difficulties on solution, and helpful to improve the identification precision.

I. Exact solution of some nonlinear evolution equations. II. The similarity solution for the Korteweg-De Vries equation and the related Painleve Transcendent

In Part I, a method for finding solutions of certain diffusive dispersive nonlinear evolution equations is introduced. The method consists of a straightforward iteration procedure, applied to the equation as it stands (in most cases), which can be carried out to all terms, followed by a summation of the resulting infinite series, sometimes directly and other times in terms of traces of inverses of operators in an appropriate space.

We first illustrate our method with Burgers' and Thomas' equations, and show how it quickly leads to the Cole-Hopft transformation, which is known to linearize these equations.

We also apply this method to the Korteweg and de Vries, nonlinear (cubic) Schrödinger, Sine-Gordon, modified KdV and Boussinesq equations. In all these cases the multisoliton solutions are easily obtained and new expressions for some of them follow. More generally we show that the Marcenko integral equations, together with the inverse problem that originates them, follow naturally from our expressions.

Only solutions that are small in some sense (i.e., they tend to zero as the independent variable goes to ∞) are covered by our methods. However, by the study of the effect of writing the initial iterate u_1 = u_(1)(x,t) as a sum u_1 = ^∼/u_1 + ^≈/u_1 when we know the solution which results if u_1 = ^∼/u_1, we are led to expressions that describe the interaction of two arbitrary solutions, only one of which is small. This should not be confused with Backlund transformations and is more in the direction of performing the inverse scattering over an arbitrary “base” solution. Thus we are able to write expressions for the interaction of a cnoidal wave with a multisoliton in the case of the KdV equation; these expressions are somewhat different from the ones obtained by Wahlquist (1976). Similarly, we find multi-dark-pulse solutions and solutions describing the interaction of envelope-solitons with a uniform wave train in the case of the Schrodinger equation.

Other equations tractable by our method are presented. These include the following equations: Self-induced transparency, reduced Maxwell-Bloch, and a two-dimensional nonlinear Schrodinger. Higher order and matrix-valued equations with nonscalar dispersion functions are also presented.

In Part II, the second Painleve transcendent is treated in conjunction with the similarity solutions of the Korteweg-de Vries equat ion and the modified Korteweg-de Vries equation.

Plasma inverse scattering theory

The object of this report is to calculate the electron density profile of plane stratified inhomogeneous plasmas. The electron density profile is obtained through a numerical solution of the inverse scattering algorithm.

The inverse scattering algorithm connects the time dependent reflected field resulting from a δ-function field incident normally on the plasma to the inhomogeneous plasma density.

Examples show that the method produces uniquely the electron density on or behind maxima of the plasma frequency.

It is shown that the δ-function incident field used in the inverse scattering algorithm can be replaced by a thin square pulse.

Formation equilibria of ternary metal complexes with citric acid and glutamine (alanine) in aqueous solution

The species and their formation constants in the ternary, systems were obtained by the Scogs2 software from potentiometric titration data. The Comics software was used to calculate the distribution of species in the ternary systems. MLXH, MLXH2 and MLXH3 are the common species in these systems. The coordination behaviors of the rare earths are very similar and their stability is closely matched. The ternary rare earth complexes are more stable than the corresponding ternary complexes of calcium. The ternary zinc complex with glutamine as the secondary ligand is more stable than the corresponding complexes of rare earths, but the ternary complex with alanine as the secondary ligand shows an inverse trend. The distributions of species in the ternary systems vary with pH changing. A prediction can be made that exogenous rare earths can affect the species of Ca and Zn in human body.

Formation and characterisation of ordered porous vanadium oxide inverse opal materials for Li-ion batteries

This thesis presents several routes towards achieving artificial opal templates by colloidal self-assembly of polystyrene (PS) or poly(methyl methacrylate) (PMMA) spheres and the use of these template for the fabrication of V2O5 inverse opals as cathode materials for lithium ion battery applications. First, through the manipulation of different experimental factors, several methods of affecting or directing opal growth towards realizing different structures, improving order and/or achieving faster formation on a variety of substrates are presented. The addition of the surfactant sodium dodecyl sulphate (SDS) at a concentration above the critical micelle concentration for SDS to a 5 wt% solution of PMMA spheres before dip-coating is presented as a method of achieving ordered 2D PhC monolayers on hydrophobic Au-coated silicon substrates at fast and slow rates of withdrawal. The effect that the degree of hydrophilicity of glass substrates has on the ordering of PMMA spheres is next investigated for a slow rate of withdrawal under noise agitation. Heating of the colloidal solution is also presented as a means of affecting order and thickness of opal deposits formed using fast rate dip coating. E-beam patterned substrates are shown as a means of altering the thermodynamically favoured FCC ordering of polystyrene spheres (PS) when dip coated at slow rate. Facile routes toward the synthesis of ordered V2O5 inverse opals are presented with direct infiltration of polymer sphere templates using liquid precursor. The use of different opal templates, both 2D and 3D partially ordered templates, is compared and the composition and arrangement of the subsequent IO structures post infiltration and calcination for various procedures is characterised. V2O5 IOs are also synthesised by electrodeposition from an aqueous VOSO4 solution at constant voltage. Electrochemical characterisation of these structures as cathode material for Li-ion batteries is assessed in a half cell arrangement for samples deposited on stainless steel foil substrates. Improved rate capabilities are demonstrated for these materials over bulk V2O5, with the improvement attributed to the shorter Li ion diffusion distances and increased electrolyte infiltration provided by the IO structure.

3D vanadium oxide inverse opal growth by electrodeposition

Three-dimensional vanadium pentoxide (V2O5) material architectures in the form of inverse opals (IOs) were fabricated using a simple electrodeposition process into artificial opal templates on stainless steel foil using an aqueous solution of VOSO4.χH2O with added ethanol. The direct deposition of V2O5 IOs was compared with V2O5 planar electrodeposition and confirms a similar progressive nucleation and growth mechanism. An in-depth examination of the chemical and morphological nature of the IO material was performed using X-ray crystallography, X-ray photoelectron spectroscopy, Raman scattering and scanning/transmission electron microscopy. Electrodeposition is demonstrated to be a function of the interstitial void fraction of the artificial opal and ionic diffusivity that leads to high quality, phase pure V2O5 inverse opals is not adversely affected by diffusion pathway tortuosity. Methods to alleviate electrodeposited overlayer formation on the artificial opal templates for the fabrication of the porous 3D structures are also demonstrated. Such a 3D material is ideally suited as a cathode for lithium ion batteries, electrochromic devices, sensors and for applications requiring high surface area electrochemically active metal oxides.

Positive regularised solutions in electromagnetic inverse scattering

Fredholm integral equations of the first kind are the mathematical model common to several electromagnetic, optical and acoustical inverse scattering problems. In most of these problems the solution must be positive in order to satisfy physical plausibility. We consider ill-posed deconvolution problems and investigate several linear regularization algorithms which provide positive approximate solutions at least in the absence of errors on the data.

Linear inverse problems with discrete data: II. Stability and regularization

For pt.I. see ibid. vol.1, p.301 (1985). In the first part of this work a general definition of an inverse problem with discrete data has been given and an analysis in terms of singular systems has been performed. The problem of the numerical stability of the solution, which in that paper was only briefly discussed, is the main topic of this second part. When the condition number of the problem is too large, a small error on the data can produce an extremely large error on the generalised solution, which therefore has no physical meaning. The authors review most of the methods which have been developed for overcoming this difficulty, including numerical filtering, Tikhonov regularisation, iterative methods, the Backus-Gilbert method and so on. Regularisation methods for the stable approximation of generalised solutions obtained through minimisation of suitable seminorms (C-generalised solutions), such as the method of Phillips (1962), are also considered.

Inverse planned constant dose rate volumetric modulated arc therapy (VMAT) as an efficient alternative to five-field intensity modulated radiation therapy (IMRT) for prostate

Purpose: The aim of this work was to determine if volumetric modulated arc therapy (VMAT) plans, created for constant dose-rate (cdrVMAT) delivery are a viable alternative to step and shoot five-field intensity modulated radiation therapy (IMRT). Materials and methods: The cdrVMAT plans, inverse planned on a treatment planning system with no solution to account for couch top or rails, were created for delivery on a linear accelerator with no variable dose rate control system. A series of five-field IMRT and cdrVMAT plans were created using dual partial arcs (gantry rotating between 260° and 100°) with 4° control points for ten prostate patients with the average rectal constraint incrementally increased. Pareto fronts were compared for the planning target volume homogeneity and average rectal dose between the two techniques for each patient. Also investigated were tumour control probability and normal tissue complication probability values for each technique. The delivery parameters [monitor units (MU) and time] and delivery accuracy of the IMRT and VMAT plans were also compared. Results: Pareto fronts showed that the dual partial arc plans were superior to the five-field IMRT plans, particularly for the clinically acceptable plans where average rectal doses were less for rotational plans (p = 0·009) with no statistical difference in target homogeneity. The cdrVMAT plans had significantly more MU (p = 0·005) but the average delivery time was significantly less than the IMRT plans by 42%. All clinically acceptable cdrVMAT plans were accurate in their delivery (gamma 99·2 ± 1·1%, 3%3 mm criteria). Conclusions Accurate delivery of dual partial arc cdrVMAT avoiding the couch top and rails has been demonstrated.

Strategies for Inverse Profiling of Two Dimensional Dielectric Scatterers Using Electromagnetic Illumination

Electromagnetic tomography has been applied to problems in nondestructive evolution, ground-penetrating radar, synthetic aperture radar, target identification, electrical well logging, medical imaging etc. The problem of electromagnetic tomography involves the estimation of cross sectional distribution dielectric permittivity, conductivity etc based on measurement of the scattered fields. The inverse scattering problem of electromagnetic imaging is highly non linear and ill posed, and is liable to get trapped in local minima. The iterative solution techniques employed for computing the inverse scattering problem of electromagnetic imaging are highly computation intensive. Thus the solution to electromagnetic imaging problem is beset with convergence and computational issues. The attempt of this thesis is to develop methods suitable for improving the convergence and reduce the total computations for tomographic imaging of two dimensional dielectric cylinders illuminated by TM polarized waves, where the scattering problem is defmed using scalar equations. A multi resolution frequency hopping approach was proposed as opposed to the conventional frequency hopping approach employed to image large inhomogeneous scatterers. The strategy was tested on both synthetic and experimental data and gave results that were better localized and also accelerated the iterative procedure employed for the imaging. A Degree of Symmetry formulation was introduced to locate the scatterer in the investigation domain when the scatterer cross section was circular. The investigation domain could thus be reduced which reduced the degrees of freedom of the inverse scattering process. Thus the entire measured scattered data was available for the optimization of fewer numbers of pixels. This resulted in better and more robust reconstructions of the scatterer cross sectional profile. The Degree of Symmetry formulation could also be applied to the practical problem of limited angle tomography, as in the case of a buried pipeline, where the ill posedness is much larger. The formulation was also tested using experimental data generated from an experimental setup that was designed. The experimental results confirmed the practical applicability of the formulation.