Existence and uniqueness for a nonlinear inverse reaction‐diffusion problem with a nonlinear source in higher dimensions

This paper analyzes the existence and the uniqueness problem for an n‐dimensional nonlinear inverse reaction‐diffusion problem with a nonlinear source. A transformation is used to obtain a new inverse coefficient problem. Then, a parabolic differential operator L_λ is defined to establish the relation between the solution of L_λ = 0 and the new inverse problem. Following this, it is shown that the inverse problem has at least one solution in the class of admissible coefficients. Furthermore, it is proved that this solution is the unique solution of the undertaken inverse problem. A numerical example is given to illustrate ill‐posedness of the inverse problem. Copyright © 2013 John Wiley & Sons, Ltd.     

A coupled complex boundary method for an inverse conductivity problem with one measurement

We recently proposed in [Cheng, XL et al. A novel coupled complex boundary method for inverse source problems Inverse Problem 2014 30 055002] a coupled complex boundary method (CCBM) for inverse source problems. In this paper, we apply the CCBM to inverse conductivity problems (ICPs) with one measurement. In the ICP, the diffusion coefficient q is to be determined from both Dirichlet and Neumann boundary data. With the CCBM, q is sought such that the imaginary part of the solution of a forward Robin boundary value problem vanishes in the problem domain. This brings in advantages on robustness and computation in reconstruction. Based on the complex forward problem, the Tikhonov regularization is used for a stable reconstruction. Some theoretical analysis is given on the optimization models. Several numerical examples are provided to show the feasibility and usefulness of the CCBM for the ICP. It is illustrated that as long as all the subdomains share some portion of the boundary, our CCBM-based Tikhonov regularization method can reconstruct the diffusion parameters stably and effectively.     

A conformal mapping algorithm for the Bernoulli free boundary value problem

We propose a new numerical method for the solution of the Bernoulli free boundary value problem for harmonic functions in a doubly connected domain D in where an unknown free boundary Γ₀ is determined by prescribed Cauchy data on Γ₀ in addition to a Dirichlet condition on the known boundary Γ₁. Our main idea is to involve the conformal mapping method as proposed and analyzed by Akduman, Haddar, and Kress for the solution of a related inverse boundary value problem. For this, we interpret the free boundary Γ₀ as the unknown boundary in the inverse problem to construct Γ₀ from the Dirichlet condition on Γ₀ and Cauchy data on the known boundary Γ₁. Our method for the Bernoulli problem iterates on the missing normal derivative on Γ₁ by alternating between the application of the conformal mapping method for the inverse problem and solving a mixed Dirichlet–Neumann boundary value problem in D. We present the mathematical foundations of our algorithm and prove a convergence result. Some numerical examples will serve as proof of concept of our approach. Copyright © 2015 John Wiley & Sons, Ltd.     

A numerical method for solving an inverse thermoacoustic problem

In this paper, we consider an inverse problem of determining the initial condition of an initial boundary value problem for the wave equation with some additional information about solving a direct initial boundary value problem. The information is obtained from measurements at the boundary of the solution domain. The purpose of our paper is to construct a numerical algorithm for solving the inverse problem by an iterative method called a method of simple iteration (MSI) and to study the resolution quality of the inverse problem as a function of the number and location of measurement points. Three two-dimensional inverse problem formulations are considered. The results of our numerical calculations are presented. It is shown that the MSI decreases the objective functional at each iteration step. However, due to the ill-posedness of the inverse problem the difference between the exact and approximate solutions decreases up to some fixed number k _min, and then monotonically increases. This shows the regularizing properties of the MSI, and the iteration number can be considered a regularization parameter.     

A Study of the Direct Spectral Transform for the Defocusing Davey-Stewartson II Equation the Semiclassical Limit

The defocusing Davey-Stewartson II equation has been shown in numerical experiments to exhibit behavior in the semiclassical limit that qualitatively resembles that of its one-dimensional reduction, the defocusing nonlinear Schrödinger equation, namely the generation from smooth initial data of regular rapid oscillations occupying domains of space-time that become well-defined in the limit. As a first step to studying this problem analytically using the inverse scattering transform, we consider the direct spectral transform for the defocusing Davey-Stewartson II equation for smooth initial data in the semiclassical limit. The direct spectral transform involves a singularly perturbed elliptic Dirac system in two dimensions. We introduce a WKB-type method for this problem, proving that it makes sense formally for sufficiently large values of the spectral parameter k by controlling the solution of an associated nonlinear eikonal problem, and we give numerical evidence that the method is accurate for such k in the semiclassical limit. Producing this evidence requires both the numerical solution of the singularly perturbed Dirac system and the numerical solution of the eikonal problem. The former is carried out using a method previously developed by two of the authors, and we give in this paper a new method for the numerical solution of the eikonal problem valid for sufficiently large k. For a particular potential we are able to solve the eikonal problem in closed form for all k, a calculation that yields some insight into the failure of the WKB method for smaller values of k. Informed by numerical calculations of the direct spectral transform, we then begin a study of the singularly perturbed Dirac system for values of k so small that there is no global solution of the eikonal problem. We provide a rigorous semiclassical analysis of the solution for real radial potentials at k=0, which yields an asymptotic formula for the reflection coefficient at k=0 and suggests an annular structure for the solution that may be exploited when k ≠ 0 is small. The numerics also suggest that for some potentials the reflection coefficient converges pointwise as ɛ↓ 0 to a limiting function that is supported in the domain of k-values on which the eikonal problem does not have a global solution. It is expected that singularities of the eikonal function play a role similar to that of turning points in the one-dimensional theory. © 2019 Wiley Periodicals, Inc.     

Isospectral flow method for nonnegative inverse eigenvalue problem with prescribed structure

The nonnegative inverse eigenvalue problem is that given a family of complex numbers λ={λ₁,…,λ_n}, find a nonnegative matrix of order n with spectrum λ. This problem is difficult and remains unsolved partially. In this paper, we focus on its generalization that the reconstructed nonnegative matrices should have some prescribed entries. It is easy to see that this new problem will come back to the common nonnegative inverse eigenvalue problem if there is no constraint of the locations of entries. A numerical isospectral flow method which is developed by hybridizing the optimization theory and steepest descent method is used to study the reconstruction. Moreover, an error estimate of the numerical iteration for ordinary differential equations on the matrix manifold is presented. After that, a numerical method for the nonnegative symmetric inverse eigenvalue problem with prescribed entries and its error estimate are considered. Finally, the approaches are verified by the numerical test results.     

A Richardson-type iterative approach for identification of delamination boundaries

A direct problem of applied mathematical modelling is to determine the response of a system given the governing partial differential equations, the geometry of interest, the complete boundary and initial conditions, and material properties. When one or more of the conditions for the solution of the direct problem are unknown, an inverse problem can be formulated. One of the methods frequently used for the solution of inverse problems involves finding the values of the unknowns in a mathematical formulation such that the behavior calculated with the model matches the measured response to a degree evaluated in terms of the classical L ₂ norm. Considered in this sense, the inverse problem is equivalent to an ill-posed optimization problem for the estimation of parameters whose solution in the majority of cases is a real mathematical challenge. In this contribution, we report a novel approach that avoids the mathematical difficulties inspired by the ill-posed character of the model. Our method is devoted to the computation of inverse problems furnished by second-order elliptical systems of partial differential equations and falls in the same conceptual line with the method initiated by Kozlov et al. and further extended and algorithmized by Weikl et al. We construct and employ a weak version of the algorithm found by Weikl et al. Proofs for the convergence and regularity of this version are given for the case of a single layer. The computational realization of the algorithm (called briefly AICRA) is applied and numerical results are obtained. Comparison with experiments demonstrates a good significance and representativeness. __________ Translated from Fundamentalnaya i Prikladnaya Matematika, Vol. 12, No. 4, pp. 209–230, 2006.     

Numerical solution of the Rosenau–KdV–RLW equation by operator splitting techniques based on B‐spline collocation method

In the present study, the operator splitting techniques based on the quintic B‐spline collocation finite element method are presented for calculating the numerical solutions of the Rosenau–KdV–RLW equation. Two test problems having exact solutions have been considered. To demonstrate the efficiency and accuracy of the present methods, the error norms L₂ and L_∞ with the discrete mass Q and energy E conservative properties have been calculated. The results obtained by the method have been compared with the exact solution of each problem and other numerical results in the literature, and also found to be in good agreement with each other. A Fourier stability analysis of each presented method is also investigated.     

Inverse problems and solution methods for a class of nonlinear complementarity problems

In this paper, motivated by the KKT optimality conditions for a sort of quadratic programs, we first introduce a class of nonlinear complementarity problems (NCPs). Then we present and discuss a kind of inverse problems of the NCPs, i.e., for a given feasible decision [`(x)]\bar{x} , we aim to characterize the set of parameter values for which there exists a point [`(y)]\bar{y} such that ([`(x)],[`(y)])(\bar{x},\bar{y}) forms a solution of the NCP and require the parameter values to be adjusted as little as possible. This leads to an inverse optimization problem. In particular, under ℓ _∞, ℓ ₁ and Frobenius norms as well as affine maps, this paper presents three simple and efficient solution methods for the inverse NCPs. Finally, some preliminary numerical results show that the proposed methods are very promising.     

Application of Rothe’s Method to a Parabolic Inverse Problem with Nonlocal Boundary Condition

We consider an inverse problem for the determination of a purely time-dependent source in a semilinear parabolic equation with a nonlocal boundary condition. An approximation scheme for the solution together with the well-posedness of the problem with the initial value u₀ ∈ H¹(Ω) is presented by means of the Rothe time-discretization method. Further approximation scheme via Rothe’s method is constructed for the problem when u₀ ∈ L²(Ω) and the integral kernel in the nonlocal boundary condition is symmetric.

     

An improved Serre model: Efficient simulation and comparative evaluation

The so-called Serre or Green and Naghdi equations are a well-known set of fully nonlinear and weakly dispersive equations that describe the propagation of long surface waves in shallow water. In order to extend its range of application to intermediate water depths, some modifications have been proposed in the literature. In this work, we analyze a new Serre model with improved linear dispersion characteristics. This new Serre system, herein denoted by Serre_{α, β}, presents additional terms of dispersive origin, thus extending its applicability to more general depth to wavelength ratios.A careful development of the Serre_{α, β} model allows a straightforward and efficient numerical implementation. This model is suitable for numerical integration by a splitting strategy which requires the solution of a hyperbolic problem and a dispersive problem. The hyperbolic part is discretized using a high-order finite volume method. For the dispersive part standard finite differences are used. A set of numerical experiments are conducted to validate the Serre_{α, β} model and to test the robustness of our numerical scheme. Theoretical solutions and benchmark experimental data are used. Moreover, comparisons against the classical Serre equations and against another well established Serre model with improved dispersion characteristics are also made.     

Splitting-based conjugate gradient method for a multi-dimensional linear inverse heat conduction problem

In this paper we consider a multi-dimensional inverse heat conduction problem with time-dependent coefficients in a box, which is well-known to be severely ill-posed, by a variational method. The gradient of the functional to be minimized is obtained by the aid of an adjoint problem, and the conjugate gradient method with a stopping rule is then applied to this ill-posed optimization problem. To enhance the stability and the accuracy of the numerical solution to the problem, we apply this scheme to the discretized inverse problem rather than to the continuous one. The difficulties with large dimensions of discretized problems are overcome by a splitting method which only requires the solution of easy-to-solve one-dimensional problems. The numerical results provided by our method are very good and the techniques seem to be very promising.     

Conservation laws with a random source

We study the scalar conservation law with a noisy nonlinear source, namely,u _l + f(u)_x = h(u, x, t) + g(u)W(t), whereW(t) is the white noise in the time variable, and we analyse the Cauchy problem for this equation where the initial data are assumed to be deterministic. A method is proposed to construct approximate weak solutions, and we then show that this yields a convergent sequence. This sequence converges to a (pathwise) solution of the Cauchy problem. The equation can be considered as a model of deterministic driven phase transitions with a random perturbation in a system of two constituents. Finally we show some numerical results motivated by two-phase flow in porous media. This research has been supported by VISTA (a research cooperation between the Norwegian Academy of Science and Letters and Den norske stats oljeselskap, Statoil) and NAVF (the Norwegian Research Council for Science and the Humanities).     

An inverse problem of identifying the time-dependent potential in a fourth-order pseudo-parabolic equation from additional condition

The aim of this work is to identify numerically, for the first time, the time-dependent potential coefficient in a fourth-order pseudo-parabolic equation with nonlocal initial data, nonlocal boundary conditions, and the boundary data as overdetermination condition. This problem emerges significantly in the modeling of various phenomena in physics and engineering. From literature we already know that this inverse problem has a unique solution. However, the problem is still ill-posed by being unstable to noise in the input data. For the numerical realization, we apply the quintic B-spline (QB-spline) collocation method for discretizing the pseudo-parabolic problem and the Tikhonov regularization for finding a stable and accurate solution. The resulting nonlinear minimization problem is solved using the MATLAB subroutine lsqnonlin. Moreover, the von Neumann stability analysis is also discussed.     

A reconstruction method for a two-dimensional inverse eigenvalue problem

In this paper we describe a method for constructing approximate solutions of a two-dimensional inverse eigenvalue problem. Here we consider the problem of recovering a functionq(x, y) from the eigenvalues of — +q(x, y) on a rectangle with Dirichlet boundary conditions. The potentialq(x, y) is assumed to be symmetric with respect to the midlines of the rectangle. Our method is a generalization of an algorithm Hald presented for the construction of symmetric potentials in the one-dimensional inverse Sturm-Liouville problem. Using a projection method, the inverse spectral problem is reduced to an inverse eigenvalue problem for a matrix. We show that if the given eigenvalues are small perturbations of simple eigenvalues ofq=0, then the matrix problem has a solution. This solution is used to construct a functionq which has the same lowest eigenvalues as the unknownq, and several numerical examples are given to illustrate the methods.     

Flow and shape reconstructions from remote measurements

We develop a point source method (PSM) to obtain flow field reconstructions from remote measurements. The PSM belongs to the class of decomposition methods in inverse scattering because it solves the nonlinear and ill‐posed inverse shape reconstruction problem by a decomposition into a linear ill‐posed problem and a nonlinear well‐posed problem. As a model problem, we investigate the reconstruction of the flow field of two‐dimensional stationary Oseen equation, which is obtained by linearizing the Navier–Stokes equation with kinematic viscosity μ > 0 around the constant velocity u₀. In contrast to acoustics or electromagnetics, the use of the PSM in fluid dynamics leads to a number of challenges in terms of the analysis and the proper setup of the scheme, in particular, because the null‐spaces of the integral operators under consideration are no longer trivial and the fundamental solution is not symmetric in its spatial coordinate. We provide a suitable formulation of the method and prove convergence of flow reconstructions by the PSM. For the realization of the reconstruction when the inclusions are not known, we employ domain sampling. We will demonstrate the feasibility of the method for reconstructing one or several objects by numerical examples. Copyright © 2012 John Wiley & Sons, Ltd.     

Marching schemes for inverse acoustic scattering problems

Summary For the numerical solution of inverse Helmholtz problems the boundary value problem for a Helmholtz equation with spatially variable wave number has to be solved repeatedly. For large wave numbers this is a challenge. In the paper we reformulate the inverse problem as an initial value problem, and describe a marching scheme for the numerical computation that needs only n² log n operations on an n × n grid. We derive stability and error estimates for the marching scheme. We show that the marching solution is close to the low-pass filtered true solution. We present numerical examples that demonstrate the efficacy of the marching scheme.     

A hybrid radial basis functions collocation technique to numerically solve fractional advection–diffusion models

In this work, we propose a hybrid radial basis functions (RBFs) collocation technique for the numerical solution of fractional advection–diffusion models. In the formulation of hybrid RBFs (HRBFs), there exist shape parameter (c* ) and weight parameter (ϵ) that control numerical accuracy and stability. For these parameters, an adaptive algorithm is developed and validated. The proposed HRBFs method is tested for numerical solutions of some fractional Black–Sholes and diffusion models. Numerical simulations performed for several benchmark problems verified the proposed method accuracy and efficiency. The quantitative analysis is made in terms of L_∞, L₂, L_rms , and L_rel error norms as well as number of nodes N over space domain and time-step δt. Numerical convergence in space and time is also studied for the proposed method. The unconditional stability of the proposed HRBFs scheme is obtained using the von Neumann methodology. It is observed that the HRBFs method circumvented the ill-conditioning problem greatly, a major issue in the Kansa method.     

On a backward problem for fractional diffusion equation with Riemann-Liouville derivative

In the present paper, we study the initial inverse problem (backward problem) for a two-dimensional fractional differential equation with Riemann-Liouville derivative. Our model is considered in the random noise of the given data. We show that our problem is not well-posed in the sense of Hadamard. A truncated method is used to construct an approximate function for the solution (called the regularized solution). Furthermore, the error estimate of the regularized solution in L² and H^τ norms is considered and illustrated by numerical example.     

On approximation in the fictitious domain method for a bidimensional transport equation

In this article a numerical method for solving a two‐dimensional transport equation in the stationary case is presented. Using the techniques of the variational calculus, we find the approximate solution for a homogeneous boundary‐value problem that corresponds to a square domain D₂. Then, using the method of the fictitious domain, we extend our algorithm to a boundary value problem for a set D that has an arbitrary shape. In this approach, the initial computation domain D (called physical domain) is immersed in a square domain D₂. We prove that the solution obtained by this method is a good approximation of the exact solution. The theoretical results are verified with the help of a numerical example. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010