A CG-Type Method for Inverse Quadratic Eigenvalue Problems in Model Updating of Structural Dynamics

In this paper we first present a CG-type method for inverse eigenvalue problem of constructing real and symmetric matrices $M$, $D$ and $K$ for the quadratic pencil $Q(\lambda)=\lambda^2M+\lambda D+K$, so that $Q(\lambda)$ has a prescribed subset of eigenvalues and eigenvectors. This method can determine the solvability of the inverse eigenvalue problem automatically. We then consider the least squares model for updating a quadratic pencil $Q(\lambda)$. More precisely, we update the model coefficient matrices $M$, $C$ and $K$ so that (i) the updated model reproduces the measured data, (ii) the symmetry of the original model is preserved, and (iii) the difference between the analytical triplet $(M, D, K)$ and the updated triplet $(M_{\text{new}}, D_{\text{new}}, K_{\text{new}})$ is minimized. In this paper a computationally efficient method is provided for such model updating and numerical examples are given to illustrate the effectiveness of the proposed method.     

A Novel Perturbative Iteration Algorithm for Effective and Efficient Solution of Frequency-Dependent Eigenvalue Problems

Many engineering structures exhibit frequency dependent characteristics and analyses of these structures lead to frequency dependent eigenvalue problems. This paper presents a novel perturbative iteration (PI) algorithm which can be used to effectively and efficiently solve frequency dependent eigenvalue problems of general frequency dependent systems. Mathematical formulations of the proposed method are developed and based on these formulations, a computer algorithm is devised. Extensive numerical case examples are given to demonstrate the practicality of the proposed method. When all modes are included, the method is exact and when only a subset of modes are used, very accurate results are obtained.     

Tailored Finite Point Method for Numerical Solutions of Singular Perturbed Eigenvalue Problems

We propose two variants of tailored finite point (TFP) methods for discretizing two dimensional singular perturbed eigenvalue (SPE) problems. A continuation method and an iterative method are exploited for solving discretized systems of equations to obtain the eigen-pairs of the SPE. We study the analytical solutions of two special cases of the SPE, and provide an asymptotic analysis for the solutions. The theoretical results are verified in the numerical experiments. The numerical results demonstrate that the proposed schemes effectively resolve the delta function like of the eigenfunctions on relatively coarse grid.     

Eigenvalue Problems of non-Hermitian Systems via Improved Asymptotic Iteration Method

We simply use the relation between the asymptotic iteration method and the Nikiforov-Uvarov method for the analytical solution of the second order linear ordinary differential equations. We apply this relation to study the Schroedinger equation with potentials admitting quasinormal modes. Non-Hermitian PT symmetric potentials have also been studied. Energy eigenvalues in all the cases by the relation are found to be consistent with exact results.     

Inverse Spectral Results for AKNS Systems with Partial Information on the Potentials

For the AKNS operator on L ²([0,1],C ²) it is well known that the data of two spectra uniquely determine the corresponding potential a.e. on [0,1] (Borg's type Theorem). We prove that, in the case where is a-priori known on [a,1], then only a part (depending on a) of two spectra determine on [0,1]. Our results include generalizations for Dirac systems of classical results obtained by Hochstadt and Lieberman for the Sturm–Liouville case, where they showed that half of the potential and one spectrum determine all the potential functions. An important ingredient in our strategy is the link between the rate of growth of an entire function and the distribution of its zeros.     

Dynamical r-Matrices for Hitchin's Systems on Schottky Curves

We express Hitchin's systems on curves in Schottky parametrization, and construct dynamical r-matrices attached to them.     

Boundary Value Problems for Boussinesq Type Systems

We characterise the boundary conditions that yield a linearly well posed problem for the so-called KdV–KdV system and for the classical Boussinesq system. Each of them is a system of two evolution PDEs modelling two-way propagation of water waves. We study these problems with the spatial variable in either the half-line or in a finite interval. The results are obtained by extending a spectral transform approach, recently developed for the analysis of scalar evolution PDEs, to the case of systems of PDEs. The knowledge of the boundary conditions that should be imposed in order for the problem to be linearly well posed can be used to obtain an integral representation of the solution. This knowledge is also necessary in order to conduct numerical simulations for the fully nonlinear systems. Mathematics Subject Classifications (2000) 34A30, 34A34, 35F10.     

Deterministic Models of the Simplest Chemical Reactions

We present a general mathematical framework for constructing deterministic models of simple chemical reactions. In such a model, an underlying dynamical system drives a process in which a particle undergoes a reaction (changes color) when it enters a certain subset (the catalytic site) of the phase space and (possibly) some other conditions are satisfied. The framework we suggest allows us to define the entropy of reaction precisely and does not rely, as was the case in previous studies, on a stochastic mechanism to generate additional entropy. Thus our approach provides a natural setting in which to derive macroscopic chemical reaction laws from microscopic deterministic dynamics without invoking any random mechanisms.     

Inverse statistical problems: from the inverse Ising problem to data science

Inverse problems in statistical physics are motivated by the challenges of ‘big data’ in different fields, in particular high-throughput experiments in biology. In inverse problems, the usual procedure of statistical physics needs to be reversed: Instead of calculating observables on the basis of model parameters, we seek to infer parameters of a model based on observations. In this review, we focus on the inverse Ising problem and closely related problems, namely how to infer the coupling strengths between spins given observed spin correlations, magnetizations, or other data. We review applications of the inverse Ising problem, including the reconstruction of neural connections, protein structure determination, and the inference of gene regulatory networks. For the inverse Ising problem in equilibrium, a number of controlled and uncontrolled approximate solutions have been developed in the statistical mechanics community. A particularly strong method, pseudolikelihood, stems from statistics. We also review the inverse Ising problem in the non-equilibrium case, where the model parameters must be reconstructed based on non-equilibrium statistics.     

An Improvement of the Asymptotic Iteration Method for Exactly Solvable Eigenvalue Problems

We derive a formula that simplifies the original asymptotic iteration method formulation to find the energy eigenvalues for the analytically solvable cases. We then show that there is a connection between the asymptotic iteration and the Nikiforo-Uvarov methods, which both solve the second order linear ordinary differential equations analytically.     

Inverse problem of the variational calculus for higher KdV equations

It is shown that the usual Hamilton's variational principle supplemented by the methodology of the integer-programming problem can be used to construct expressions for the Lagrangian densities of higher KdV fields. This is demonstrated with special emphasis on the second and third members of the hierarchy. However, the method is general enough for applications to equations of any order. The expressions for Lagrangian densities are used to calculate results for Hamiltonian densities that characterize Zakharov-Faddeev-Gardner equation. Received 27 January 2002 / Received in final form 6 May 2002 Published online 24 September 2002     

Dynamics of twists on antiferromagnetic spin chains: Theory

The equation of motion of twists on classical antiferromagnetic Heisenberg spin chains are derived. It is shown that twists interact via position- and velocity-dependent long-range two-body forces. A quiescent regime is identified wherein the system conserves momentum. With increasing kinetic energy the system exits this regime and momentum conservation is violated due to walls annihilation. A bitwist system is shown to be integrable and its exact solution is analysed. Many-twist systems are discussed and novel periodic twist lattice solutions are found on closed chains. The stability of these solutions is discussed. Received 12 June 2002 Published online 2 October 2002 RID="a" ID="a"e-mail: rbbll@phy.cam.ac.uk     

An Alternative Approach to Obtain the Exact Wave Functions of Time-Dependent Hamiltonian Systems Involving Quadratic, Inverse Quadratic, and (1/x)p+p(1/x) Terms

By using the Lewis-Riesenfeld theory and algebraic method, we present an alternative approach to obtain the exact solution of time-dependent Hamiltonian systems involving quadratic, inverse quadratic and (1/x)p+p(1/x) terms. This solution is discussed and compared with that obtained by Choi, J. R. (2003). International Journal of Theoretical Physics 42, 853]. PACS: 03.65Ge; 03.65Fd; 03.65Bz     

Dynamics of a single peak of the Rosensweig instability in a magnetic fluid

To describe the dynamics of a single peak of the Rosensweig instability a model is proposed which approximates the peak by a half-ellipsoid atop a layer of magnetic fluid. The resulting nonlinear equation for the height of the peak leads to the correct subcritical character of the bifurcation for static induction. For a time-dependent induction the effects of inertia and damping are incorporated. The results of the model show qualitative agreement with the experimental findings, as in the appearance of period doubling, trebling, and higher multiples of the driving period. Furthermore, a quantitative agreement is also found for the parameter ranges of frequency and induction in which these phenomena occur.     

Numerical Approximation of Hopf Bifurcation for Tumor-Immune System Competition Model with Two Delays

This paper is concerned with the Hopf bifurcation analysis of tumor-immune system competition model with two delays. First, we discuss the stability of state points with different kinds of delays. Then, a sufficient condition to the existence of the Hopf bifurcation is derived with parameters at different points. Furthermore, under this condition, the stability and direction of bifurcation are determined by applying the normal form method and the center manifold theory. Finally, a kind of Runge-Kutta methods is given out to simulate the periodic solutions numerically. At last, some numerical experiments are given to match well with the main conclusion of this paper.     

Dynamics and stability of a new class of periodic solutions of the optical parametric oscillator

The stability and dynamics of a new class of periodic solutions is investigated when a degenerate optical parametric oscillator system is forced by an external pumping field with a periodic spatial profile modeled by Jacobi elliptic functions. Both sinusoidal behavior as well as localized hyperbolic (front and pulse) behavior can be considered in this model. The stability and bifurcation behaviors of these transverse electromagnetic structures are studied numerically. The periodic solutions are shown to be stabilized by the nonlinear parametric interaction between the pump and signal fields interacting with the cavity diffraction, attenuation, and periodic external pumping. Specifically, sinusoidal solutions result in robust and stable configurations while well-separated and more localized field structures often undergo bifurcation to new steady-state solutions having the same period as the external forcing. Extensive numerical simulations and studies of the solutions are provided.     

Limiting Dynamics for Spherical Models of Spin Glasses at High Temperature

We analyze the coupled non-linear integro-differential equations whose solution is the thermodynamical limit of the empirical correlation and response functions in the Langevin dynamics for spherical p-spin disordered mean-field models. We provide a mathematically rigorous derivation of their FDT solution (for the high temperature regime) and of certain key properties of this solution, which are in agreement with earlier derivations based on physical grounds. AMS (2000) Subject Classification: Primary: 82C44 Secondary: 82C31, 60H10, 60F15, 60K35     

Limiting Dynamics for Spherical Models of Spin Glasses at High Temperature

We analyze the coupled non-linear integro-differential equations whose solution is the thermodynamical limit of the empirical correlation and response functions in the Langevin dynamics for spherical p−spin disordered mean-field models. We provide a mathematically rigorous derivation of their FDT solution (for the high temperature regime) and of certain key properties of this solution, which are in agreement with earlier derivations based on physical grounds. AMS (2000) Subject Classification: Primary: 82C44, Secondery: 82C31, 60H10, 60F15, 60K35     

On the use of the resting potential and level set methods for identifying ischemic heart disease: An inverse problem

The electrical activity in the heart is modeled by a complex, nonlinear, fully coupled system of differential equations. Several scientists have studied how this model, referred to as the bidomain model, can be modified to incorporate the effect of heart infarctions on simulated ECG (electrocardiogram) recordings.We are concerned with the associated inverse problem; how can we use ECG recordings and mathematical models to identify the position, size and shape of heart infarctions? Due to the extreme CPU efforts needed to solve the bidomain equations, this model, in its full complexity, is not well-suited for this kind of problems. In this paper we show how biological knowledge about the resting potential in the heart and level set techniques can be combined to derive a suitable stationary model, expressed in terms of an elliptic PDE, for such applications. This approach leads to a nonlinear ill-posed minimization problem, which we propose to regularize and solve with a simple iterative scheme.Finally, our theoretical findings are illuminated through a series of computer simulations for an experimental setup involving a realistic heart in torso geometry. More specifically, experiments with synthetic ECG recordings, produced by solving the bidomain model, indicate that our method manages to identify the physical characteristics of the ischemic region(s) in the heart. Furthermore, the ill-posed nature of this inverse problem is explored, i.e. several quantitative issues of our scheme are explored.