NUMERICAL SOLUTION OF THE SCATTERING PROBLEM FOR ACOUSTIC WAVES BY A TWO-SIDED CRACK IN 2-DIMENSIONAL SPACE

The wave scattering problem by a crack F in R2 with impedance type boundary is considered. This problem models the diffraction of waves by thin two-sided cylindrical screens. A numerical method for solving the problem is developed. The solution of the problem is represented in the form of the combined angular potential and single-layer potential. The linear integral equations satisfied by the density functions are derived for general parameterized arcs. The weakly singular integrals and the Cauchy singular integral arising in these equations are computed using a highly accurate scheme with a truncation error analysis. The advantage of the scheme proposed in this paper is, in one hand, the fact that we do not need the analyticity property of the crack and we allow different complex valued surface impedances in both sides of the crack. In the other hand, we avoid the hyper-singular integrals. Numerical implementations showing the validity of the scheme are presented.

Identification of the Exchange Coefficient from Indirect Data for a Coupled Continuum Pipe-Flow Model

Calibration and identification of the exchange effect between the karst aquifers and the underlying conduit network are important issues in order to gain a better understanding of these hydraulic systems. Based on a coupled continuum pipe-flow(CCPF for short) model describing flows in karst aquifers, this paper is devoted to the identification of an exchange rate function, which models the hydraulic interaction between the fissured volume(matrix) and the conduit, from the Neumann boundary data, i.e., matrix/conduit seepage velocity. The authors formulate this parameter identification problem as a nonlinear operator equation and prove the compactness of the forward mapping. The stable approximate solution is obtained by two classic iterative regularization methods, namely,the Landweber iteration and Levenberg-Marquardt method. Numerical examples on noisefree and noisy data shed light on the appropriateness of the proposed approaches.

Inverse Eigenvalue Problems for Exploring the Dynamics of Systems Biology Models

This paper describes inverse eigenvalue problems that arise in studying qualitative dynamics in systems biology models.An algorithm based on lift-andproject iterations is proposed,where the lifting step entails solving a constrained matrix inverse eigenvalue problem.In particular,prior to carrying out the iterative steps,a-priori bounds on the entries of the Jacobian matrix are computed by relying on the reaction network structure as well as the form of the rate law expressions for the model under consideration.Numerical results on a number of models show that the proposed algorithm can be used to computationally explore the possible dynamical scenarios while identifying the important mechanisms via the use of sparsity-promoting regularization.

Halton-type sequences from global function fields

For any prime power q and any dimension s≥1, a new construction of (t, s)-sequences in base q using global function fields is presented. The construction yields an analog of Halton sequences for global function fields. It is the first general construction of (t, s)-sequences that is not directly based on the digital method. The construction can also be put into the framework of the theory of (u, e, s)-sequences that was recently introduced by Tezuka and leads in this way to better discrepancy bounds for the constructed sequences.

Rational Solutions of High-Order Algebraic Ordinary Differential Equations

This paper considers algebraic ordinary differential equations(AODEs)and study their polynomial and rational solutions.The authors first prove a sufficient condition for the existence of a bound on the degree of the possible polynomial solutions to an AODE.An AODE satisfying this condition is called noncritical.Then the authors prove that some common classes of low-order AODEs are noncritical.For rational solutions,the authors determine a class of AODEs,which are called maximally comparable,such that the possible poles of any rational solutions are recognizable from their coefficients.This generalizes the well-known fact that any pole of rational solutions to a linear ODE is contained in the set of zeros of its leading coefficient.Finally,the authors develop an algorithm to compute all rational solutions of certain maximally comparable AODEs,which is applicable to 78.54%of the AODEs in Kamke's collection of standard differential equations.