We consider the inverse problem to determine the shape of a open cavity embedded in the infinite ground plane from knowledge of the far-field pattern of the scattering of TM polarization. For its approximate solution ...We consider the inverse problem to determine the shape of a open cavity embedded in the infinite ground plane from knowledge of the far-field pattern of the scattering of TM polarization. For its approximate solution we propose a regularized Newton iteration scheme. For a foundation of Newton type methods we establish the Fr^chet differentiability of solution to the scattering problem with respect to the boundary of the cavity. Some numerical examples of the feasibility of the method are presented.展开更多
The electric inversion technique reconstructs the subsurface medium distribution from acquired data.On the basis of electric inversion,objects buried under the earth or seabed,such as pipelines and unexploded ordnance...The electric inversion technique reconstructs the subsurface medium distribution from acquired data.On the basis of electric inversion,objects buried under the earth or seabed,such as pipelines and unexploded ordnance,are detected and located in a contactless manner.However,the process of accurately reconstructing the shape of the target object is challenging because electric inversion is a nonlinear and ill-posed problem.In this work,we present an inverse multiquadric(IMQ)regularization method based on the level set function for reconstructing buried pipelines.In the case of locating underwater objects,the unknown inversion area is split into two parts,the background and the pipeline with known conductivity.The geometry of the pipeline is represented based on the level set function for achieving a noiseless inversion image.To obtain a binary image,the IMQ is used as the regularization term,which‘pushes’the level set function away from 0.We also provide an appropriate method to select the bandwidth and regularization parameters for the IMQ regularization term,resulting in reconstructed images with sharp edges.The simulation results and analysis show that the proposed method performs better than classical inversion methods.展开更多
In this paper,the Cauchy problem of biharmonic equation is considered.This problem is ill-posed,i.e.,the solution(if exists)does not depend on the measurable data.Firstly,we give the conditional stability result under...In this paper,the Cauchy problem of biharmonic equation is considered.This problem is ill-posed,i.e.,the solution(if exists)does not depend on the measurable data.Firstly,we give the conditional stability result under the a priori bound assumption for the exact solution.Secondly,a modified Tikhonov regularization method is used to solve this ill-posed problem.Under the a priori and the a posteriori regularization parameter choice rule,the error estimates between the regularization solutions and the exact solution are obtained.Finally,some numerical examples are presented to verify that our method is effective.展开更多
This article compares the isotropic and anisotropic TV regularizations used in inverse acoustic scattering. It is observed that compared with the traditional Tikhonov regularization, isotropic and anisotropic TV regul...This article compares the isotropic and anisotropic TV regularizations used in inverse acoustic scattering. It is observed that compared with the traditional Tikhonov regularization, isotropic and anisotropic TV regularizations perform better in the sense of edge preserving. While anisotropic TV regularization will cause distortions along axes. To minimize the energy function with isotropic and anisotropic regularization terms, we use split Bregman scheme. We do several 2D numerical experiments to validate the above arguments.展开更多
Source term identification is very important for the contaminant gas emission event. Thus, it is necessary to study the source parameter estimation method with high computation efficiency, high estimation accuracy and...Source term identification is very important for the contaminant gas emission event. Thus, it is necessary to study the source parameter estimation method with high computation efficiency, high estimation accuracy and reasonable confidence interval. Tikhonov regularization method is a potential good tool to identify the source parameters. However, it is invalid for nonlinear inverse problem like gas emission process. 2-step nonlinear and linear PSO （partial swarm optimization）-Tikhonov regularization method proposed previously have estimated the emission source parameters successfully. But there are still some problems in computation efficiency and confidence interval. Hence, a new 1-step nonlinear method combined Tikhonov regularizafion and PSO algorithm with nonlinear forward dispersion model was proposed. First, the method was tested with simulation and experiment cases. The test results showed that 1-step nonlinear hybrid method is able to estimate multiple source parameters with reasonable confidence interval. Then, the estimation performances of different methods were compared with different cases. The estimation values with 1-step nonlinear method were close to that with 2-step nonlinear and linear PSO-Tikhonov regularization method, 1-step nonlinear method even performs better than other two methods in some cases, especially for source strength and downwind distance estimation. Compared with 2-step nonlinear method, 1-step method has higher computation efficiency. On the other hand, the confidence intervals with the method proposed in this paper seem more reasonable than that with other two methods. Finally, single PSO algorithm was compared with 1-step nonlinear PSO-Tikhonov hybrid regularization method. The results showed that the skill scores of 1-step nonlinear hybrid method to estimate source parameters were close to that of single PSO method and even better in some cases. One more important property of 1-step nonlinear PSO-Tikhonov regularization method is its reasonable confidence interval, which is not obtained by single PSO algorithm. Therefore, 1-step nonlinear hybrid regularization method proposed in this paper is a potential good method to estimate contaminant gas emission source term.展开更多
Discrete-type continuation method for solving nonlinear system of equations and Tikhonov's regularization method for solving linear ill-posed problems are combined into a stable and widely convergent one for solvi...Discrete-type continuation method for solving nonlinear system of equations and Tikhonov's regularization method for solving linear ill-posed problems are combined into a stable and widely convergent one for solving nonlinear operator equations with difficultly computed and ill-conditioned derivatives. Some results about their convergence are given The application of this method to solve the inverse problem of one-dimensional diffusion equation is demonstrated.展开更多
The regularized integrodifferential equation for the first kind of Fredholm, integral equation with a complex kernel is derived by generalizing the Tikhonov regularization method and the convergence of approximate reg...The regularized integrodifferential equation for the first kind of Fredholm, integral equation with a complex kernel is derived by generalizing the Tikhonov regularization method and the convergence of approximate regularized solutions is discussed. As an application of the method, an inverse problem in the two-dimensional wave-making problem of a flat plate is solved numerically, and a practical approach of choosing optimal regularization parameter is given.展开更多
In this paper,we mainly study an inverse source problem of time fractional diffusion equation in a bounded domain with an over-specified terminal condition at a fixed time.A novel regularization method,which we call t...In this paper,we mainly study an inverse source problem of time fractional diffusion equation in a bounded domain with an over-specified terminal condition at a fixed time.A novel regularization method,which we call the exponential Tikhonov regularization method with a parameter γ,is proposed to solve the inverse source problem,and the corresponding convergence analysis is given under a-priori and a-posteriori regularization parameter choice rules.Whenγis less than or equal to zero,the optimal convergence rate can be achieved and it is independent of the value of γ.However,when γ is greater than zero,the optimal convergence rate depends on the value of γ which is related to the regularity of the unknown source.Finally,numerical experiments are conducted for showing the effectiveness of the proposed exponential regularization method.展开更多
In this paper, we consider the inverse scattering problem of reconstructing a bounded obstacle in a three-dimensional planar waveguide from the scattered near-field data measured on a finite cylindrical surface contai...In this paper, we consider the inverse scattering problem of reconstructing a bounded obstacle in a three-dimensional planar waveguide from the scattered near-field data measured on a finite cylindrical surface containing the obstacle and corresponding to infinitely many incident point sources also placed on the measurement surface. The obstacle is allowed to be an impenetrable scatterer or a penetrable scatterer. We establish the validity of the factorization method with the nearfield data to characterize the obstacle in the planar waveguide by constructing an outgoing-to-incoming operator which is an integral operator defined on the measurement surface with the kernel given in terms of an infinite series.展开更多
This paper is concerned with the inverse scattering problems for Schrdinger equations with compactly supported potentials.For purpose of reconstructing the support of the potential,we derive a factorization of the sca...This paper is concerned with the inverse scattering problems for Schrdinger equations with compactly supported potentials.For purpose of reconstructing the support of the potential,we derive a factorization of the scattering amplitude operator A and prove that the ranges of （A＊ A） ^1/4 and G which maps more general incident fields than plane waves into the scattering amplitude coincide.As an application we characterize the support of the potential using only the spectral data of the operator A.展开更多
An inverse problem for identification of the coefficient in heat-conduction equation is considered. After reducing the problem to a nonlinear ill-posed operator equation, Newton type iterative methods are considered. ...An inverse problem for identification of the coefficient in heat-conduction equation is considered. After reducing the problem to a nonlinear ill-posed operator equation, Newton type iterative methods are considered. The implicit iterative method is applied to the linearized Newton equation, and the key step in the process is that a new reasonable a posteriori stopping rule for the inner iteration is presented. Numerical experiments for the new method as well as for Tikhonov method and Bakushikskii method are given, and these results show the obvious advantages of the new method over the other ones.展开更多
Ship maneuverability, in the field of ship engineering, is often predicted by maneuvering motion group （MMG） mathematical model. Then it is necessary to determine hydrodynamic coefficients and interaction force coef...Ship maneuverability, in the field of ship engineering, is often predicted by maneuvering motion group （MMG） mathematical model. Then it is necessary to determine hydrodynamic coefficients and interaction force coefficients of the model. Based on the data of free running model test, the problem for obtaining these coefficients is called inverse one. For the inverse problem, ill-posedness is inherent, nonlinearity and great computation happen, and the computation is also insensitive, unstable and time-consuming. In the paper, a regularization method is introduced to solve ill-posed problem and genetic algorithm is used for nonlinear motion of ship maneuvering. In addition, the immunity is applied to solve the prematurity, to promote the global searching ability and to increase the converging speed. The combination of regularization method and immune genetic algorithm（RIGA） applied in MMG mathematical model, showed rapid converging speed and good stability.展开更多
This paper studies to numerical solutions of an inverse heat conduction problem.The effect of algorithms based on the Newton-Tikhonov method and the Newton-implicit iterative method is investigated,and then several mo...This paper studies to numerical solutions of an inverse heat conduction problem.The effect of algorithms based on the Newton-Tikhonov method and the Newton-implicit iterative method is investigated,and then several modifications are presented.Numerical examples show the modified algorithms always work and can greatly reduce the computational costs.展开更多
The elastic plate vibration model is studied under the external force. The size of the source term by the given mode of the source and some observations from the body of the plate is determined over a time interval, w...The elastic plate vibration model is studied under the external force. The size of the source term by the given mode of the source and some observations from the body of the plate is determined over a time interval, which is referred to be an inverse source problem of a plate equation. The uniqueness theorem for this problem is stated, and the fundamental solution to the plate equation is derived. In the case that the plate is driven by the harmonic load, the fundamental solution method （FSM） and the Tikhonov regularization technique axe used to calculate the source term. Numerical experiments of the Euler-Bernoulli beam and the Kirchhoff-Love plate show that the FSM can work well for practical use, no matter the source term is smooth or piecewise.展开更多
In this paper, we present an effective meshless method for solving the inverse heat conduction problems, with the Neumann boundary condition. A PDE-constrained optimization method is developed to get a global approxim...In this paper, we present an effective meshless method for solving the inverse heat conduction problems, with the Neumann boundary condition. A PDE-constrained optimization method is developed to get a global approximation scheme in both spatial and temporal domains, by using the fundamental solution of the governing equation as the basis function.Since the initial measured data contain some noises, and the resulting systems of equations are usually ill-conditioned, the Tikhonov regularization technique with the generalized crossvalidation criterion is applied to obtain more stable numerical solutions. It is shown that the proposed schemes are effective by some numerical tests.展开更多
For the backward diffusion equation,a stable discrete energy regularization algorithm is proposed.Existence and uniqueness of the numerical solution are given.Moreover,the error between the solution of the given backw...For the backward diffusion equation,a stable discrete energy regularization algorithm is proposed.Existence and uniqueness of the numerical solution are given.Moreover,the error between the solution of the given backward diffusion equation and the numerical solution via the regularization method can be estimated.Some numerical experiments illustrate the efficiency of the method,and its application in image deblurring.展开更多
This paper aims to apply a virtual boundary element method(VBEM)to solve the inverse problems of three-dimensional heat conduction in orthotropic media.This method avoids the singular integrations in the conventional ...This paper aims to apply a virtual boundary element method(VBEM)to solve the inverse problems of three-dimensional heat conduction in orthotropic media.This method avoids the singular integrations in the conventional boundary element method,and can be treated as a potential approach for solving the inverse problems of the heat conduction owing to the boundary-only discretization and semi-analytical algorithm.When the VBEM is applied to the inverse problems,the numerical instability may occur if a virtual boundary is not properly chosen.The method encounters a highly illconditioned matrix for the larger distance between the physical boundary and the virtual boundary,and otherwise is hard to avoid the singularity of the source point.Thus,it must adopt an appropriate regularization method to deal with the ill-posed systems of inverse problems.In this study,the VBEM and different regularization techniques are combined to model the inverse problem of three-dimensional heat conduction in orthotropic media.The proper regularization techniques not only make the virtual boundary to be allocated freer,but also solve the ill-conditioned equation of the inverse problem.Numerical examples demonstrate that the proposed method is efficient,accurate and numerically stable for solving the inverse problems of three-dimensional heat conduction in orthotropic media.展开更多
Detecting corrosion by electrical field can be modeled by a Cauchy problem of Laplace equation in annulus domain under the assumption that the thickness of the pipe is relatively small compared with the radius of the ...Detecting corrosion by electrical field can be modeled by a Cauchy problem of Laplace equation in annulus domain under the assumption that the thickness of the pipe is relatively small compared with the radius of the pipe.The interior surface of the pipe is inaccessible and the nondestructive detection is solely based on measurements from the outer layer.The Cauchy problem for an elliptic equation is a typical ill-posed problem whose solution does not depend continuously on the boundary data.In this work,we assume that the measurements are available on the whole outer boundary on an annulus domain.By imposing reasonable assumptions,the theoretical goal here is to derive the stabilities of the Cauchy solutions and an energy regularization method.Relationship between the proposed energy regularization method and the Tikhonov regularization with Morozov principle is also given.A novel numerical algorithm is proposed and numerical examples are given.展开更多
The accurate material physical properties, initial and boundary conditions are indispensable to the numerical simulation in the casting process, and they are related to the simulation accuracy directly. The inverse he...The accurate material physical properties, initial and boundary conditions are indispensable to the numerical simulation in the casting process, and they are related to the simulation accuracy directly. The inverse heat conduction method can be used to identify the mentioned above parameters based on the temperature measurement data. This paper presented a new inverse method according to Tikhonov regularization theory. A regularization functional was established and the regularization parameter was deduced, the Newton-Raphson iteration method was used to solve the equations. One detailed case was solved to identify the thermal conductivity and specific heat of sand mold and interfacial heat transfer coefficient (IHTC) at the meantime. This indicates that the regularization method is very efficient in decreasing the sensitivity to the temperature measurement data, overcoming the ill-posedness of the inverse heat conduction problem (IHCP) and improving the stability and accuracy of the results. As a general inverse method, it can be used to identify not only the material physical properties but also the initial and boundary conditions' parameters.展开更多
基金The NNSF(10626017)of Chinathe Science Foundation(11511276)of the Education Committee of Heilongjiang Provincethe Foundation(LBH-Q05114)of Heilongjiang Province
文摘We consider the inverse problem to determine the shape of a open cavity embedded in the infinite ground plane from knowledge of the far-field pattern of the scattering of TM polarization. For its approximate solution we propose a regularized Newton iteration scheme. For a foundation of Newton type methods we establish the Fr^chet differentiability of solution to the scattering problem with respect to the boundary of the cavity. Some numerical examples of the feasibility of the method are presented.
基金supported by the National Natural Sci-ence Foundation of China(No.52101383)the Fundamen-tal Research Funds for the Central Universities(No.3072021CF0802)+3 种基金the Key Laboratory of Advanced Marine Communication and Information Technology,Ministry of Industry and Information Technology(No.AMCIT2101-02)the Sino-Russian Cooperation Fund of Harbin Engi-neering University(No.2021HEUCRF006)the Ministry of Science and Higher Education of the Russian Federation(No.075-15-2020-934)the International Science&Technology Cooperation Program of China(No.2014DF R10240).
文摘The electric inversion technique reconstructs the subsurface medium distribution from acquired data.On the basis of electric inversion,objects buried under the earth or seabed,such as pipelines and unexploded ordnance,are detected and located in a contactless manner.However,the process of accurately reconstructing the shape of the target object is challenging because electric inversion is a nonlinear and ill-posed problem.In this work,we present an inverse multiquadric(IMQ)regularization method based on the level set function for reconstructing buried pipelines.In the case of locating underwater objects,the unknown inversion area is split into two parts,the background and the pipeline with known conductivity.The geometry of the pipeline is represented based on the level set function for achieving a noiseless inversion image.To obtain a binary image,the IMQ is used as the regularization term,which‘pushes’the level set function away from 0.We also provide an appropriate method to select the bandwidth and regularization parameters for the IMQ regularization term,resulting in reconstructed images with sharp edges.The simulation results and analysis show that the proposed method performs better than classical inversion methods.
基金Supported by the National Natural Science Foundation of China(Grant No.11961044)。
文摘In this paper,the Cauchy problem of biharmonic equation is considered.This problem is ill-posed,i.e.,the solution(if exists)does not depend on the measurable data.Firstly,we give the conditional stability result under the a priori bound assumption for the exact solution.Secondly,a modified Tikhonov regularization method is used to solve this ill-posed problem.Under the a priori and the a posteriori regularization parameter choice rule,the error estimates between the regularization solutions and the exact solution are obtained.Finally,some numerical examples are presented to verify that our method is effective.
文摘This article compares the isotropic and anisotropic TV regularizations used in inverse acoustic scattering. It is observed that compared with the traditional Tikhonov regularization, isotropic and anisotropic TV regularizations perform better in the sense of edge preserving. While anisotropic TV regularization will cause distortions along axes. To minimize the energy function with isotropic and anisotropic regularization terms, we use split Bregman scheme. We do several 2D numerical experiments to validate the above arguments.
基金Supported by the National Natural Science Foundation of China(21676216)China Postdoctoral Science Foundation(2015M582667)+2 种基金Natural Science Basic Research Plan in Shaanxi Province of China(2016JQ5079)Key Research Project of Shaanxi Province(2015ZDXM-GY-115)the Fundamental Research Funds for the Central Universities(xjj2017124)
文摘Source term identification is very important for the contaminant gas emission event. Thus, it is necessary to study the source parameter estimation method with high computation efficiency, high estimation accuracy and reasonable confidence interval. Tikhonov regularization method is a potential good tool to identify the source parameters. However, it is invalid for nonlinear inverse problem like gas emission process. 2-step nonlinear and linear PSO （partial swarm optimization）-Tikhonov regularization method proposed previously have estimated the emission source parameters successfully. But there are still some problems in computation efficiency and confidence interval. Hence, a new 1-step nonlinear method combined Tikhonov regularizafion and PSO algorithm with nonlinear forward dispersion model was proposed. First, the method was tested with simulation and experiment cases. The test results showed that 1-step nonlinear hybrid method is able to estimate multiple source parameters with reasonable confidence interval. Then, the estimation performances of different methods were compared with different cases. The estimation values with 1-step nonlinear method were close to that with 2-step nonlinear and linear PSO-Tikhonov regularization method, 1-step nonlinear method even performs better than other two methods in some cases, especially for source strength and downwind distance estimation. Compared with 2-step nonlinear method, 1-step method has higher computation efficiency. On the other hand, the confidence intervals with the method proposed in this paper seem more reasonable than that with other two methods. Finally, single PSO algorithm was compared with 1-step nonlinear PSO-Tikhonov hybrid regularization method. The results showed that the skill scores of 1-step nonlinear hybrid method to estimate source parameters were close to that of single PSO method and even better in some cases. One more important property of 1-step nonlinear PSO-Tikhonov regularization method is its reasonable confidence interval, which is not obtained by single PSO algorithm. Therefore, 1-step nonlinear hybrid regularization method proposed in this paper is a potential good method to estimate contaminant gas emission source term.
文摘Discrete-type continuation method for solving nonlinear system of equations and Tikhonov's regularization method for solving linear ill-posed problems are combined into a stable and widely convergent one for solving nonlinear operator equations with difficultly computed and ill-conditioned derivatives. Some results about their convergence are given The application of this method to solve the inverse problem of one-dimensional diffusion equation is demonstrated.
文摘The regularized integrodifferential equation for the first kind of Fredholm, integral equation with a complex kernel is derived by generalizing the Tikhonov regularization method and the convergence of approximate regularized solutions is discussed. As an application of the method, an inverse problem in the two-dimensional wave-making problem of a flat plate is solved numerically, and a practical approach of choosing optimal regularization parameter is given.
基金supported by National Natural Science Foundation of China(11961002,11761007,11861007)Key Project of the Natural Science Foundation of Jiangxi Province(20212ACB201001).
文摘In this paper,we mainly study an inverse source problem of time fractional diffusion equation in a bounded domain with an over-specified terminal condition at a fixed time.A novel regularization method,which we call the exponential Tikhonov regularization method with a parameter γ,is proposed to solve the inverse source problem,and the corresponding convergence analysis is given under a-priori and a-posteriori regularization parameter choice rules.Whenγis less than or equal to zero,the optimal convergence rate can be achieved and it is independent of the value of γ.However,when γ is greater than zero,the optimal convergence rate depends on the value of γ which is related to the regularity of the unknown source.Finally,numerical experiments are conducted for showing the effectiveness of the proposed exponential regularization method.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.61421062 and 61520106004)the Microsoft Research Fund of Asia
文摘In this paper, we consider the inverse scattering problem of reconstructing a bounded obstacle in a three-dimensional planar waveguide from the scattered near-field data measured on a finite cylindrical surface containing the obstacle and corresponding to infinitely many incident point sources also placed on the measurement surface. The obstacle is allowed to be an impenetrable scatterer or a penetrable scatterer. We establish the validity of the factorization method with the nearfield data to characterize the obstacle in the planar waveguide by constructing an outgoing-to-incoming operator which is an integral operator defined on the measurement surface with the kernel given in terms of an infinite series.
基金The Major State Basic Research Development Program Grant (2005CB321701)the Heilongjiang Education Committee Grant (11551364) of China
文摘This paper is concerned with the inverse scattering problems for Schrdinger equations with compactly supported potentials.For purpose of reconstructing the support of the potential,we derive a factorization of the scattering amplitude operator A and prove that the ranges of （A＊ A） ^1/4 and G which maps more general incident fields than plane waves into the scattering amplitude coincide.As an application we characterize the support of the potential using only the spectral data of the operator A.
文摘An inverse problem for identification of the coefficient in heat-conduction equation is considered. After reducing the problem to a nonlinear ill-posed operator equation, Newton type iterative methods are considered. The implicit iterative method is applied to the linearized Newton equation, and the key step in the process is that a new reasonable a posteriori stopping rule for the inner iteration is presented. Numerical experiments for the new method as well as for Tikhonov method and Bakushikskii method are given, and these results show the obvious advantages of the new method over the other ones.
文摘Ship maneuverability, in the field of ship engineering, is often predicted by maneuvering motion group （MMG） mathematical model. Then it is necessary to determine hydrodynamic coefficients and interaction force coefficients of the model. Based on the data of free running model test, the problem for obtaining these coefficients is called inverse one. For the inverse problem, ill-posedness is inherent, nonlinearity and great computation happen, and the computation is also insensitive, unstable and time-consuming. In the paper, a regularization method is introduced to solve ill-posed problem and genetic algorithm is used for nonlinear motion of ship maneuvering. In addition, the immunity is applied to solve the prematurity, to promote the global searching ability and to increase the converging speed. The combination of regularization method and immune genetic algorithm（RIGA） applied in MMG mathematical model, showed rapid converging speed and good stability.
基金Project supported by the Key Disciplines of Shanghai Municipality (Grant No.S30104)the Shanghai Leading Academic Discipline Project (Grant No.J50101)
文摘This paper studies to numerical solutions of an inverse heat conduction problem.The effect of algorithms based on the Newton-Tikhonov method and the Newton-implicit iterative method is investigated,and then several modifications are presented.Numerical examples show the modified algorithms always work and can greatly reduce the computational costs.
文摘The elastic plate vibration model is studied under the external force. The size of the source term by the given mode of the source and some observations from the body of the plate is determined over a time interval, which is referred to be an inverse source problem of a plate equation. The uniqueness theorem for this problem is stated, and the fundamental solution to the plate equation is derived. In the case that the plate is driven by the harmonic load, the fundamental solution method （FSM） and the Tikhonov regularization technique axe used to calculate the source term. Numerical experiments of the Euler-Bernoulli beam and the Kirchhoff-Love plate show that the FSM can work well for practical use, no matter the source term is smooth or piecewise.
基金Supported by the National Natural Science Foundation of China(Grant Nos.1129014311471066+3 种基金11572081)the Fundamental Research of Civil Aircraft(Grant No.MJ-F-2012-04)the Fundamental Research Funds for the Central Universities(Grant No.DUT15LK44)the Scientific Research Funds of Inner Mongolia University for the Nationalities(Grant No.NMD1304)
文摘In this paper, we present an effective meshless method for solving the inverse heat conduction problems, with the Neumann boundary condition. A PDE-constrained optimization method is developed to get a global approximation scheme in both spatial and temporal domains, by using the fundamental solution of the governing equation as the basis function.Since the initial measured data contain some noises, and the resulting systems of equations are usually ill-conditioned, the Tikhonov regularization technique with the generalized crossvalidation criterion is applied to obtain more stable numerical solutions. It is shown that the proposed schemes are effective by some numerical tests.
基金National Natural Science Foundation of China(No.10471073)。
文摘For the backward diffusion equation,a stable discrete energy regularization algorithm is proposed.Existence and uniqueness of the numerical solution are given.Moreover,the error between the solution of the given backward diffusion equation and the numerical solution via the regularization method can be estimated.Some numerical experiments illustrate the efficiency of the method,and its application in image deblurring.
基金This study was supported by“the Fundamental Research Funds for the Central Universities”(Grant No.2015B37814)the Postgraduate Research and Practice Innovation Program of Jiangsu Province(Grant No.KYLX15_0489)+1 种基金the National Natural Science Foundation of China(Grant No.51679081)“the Fundamental Research Funds for the Central Universities”(Grant No.2018B48514).
文摘This paper aims to apply a virtual boundary element method(VBEM)to solve the inverse problems of three-dimensional heat conduction in orthotropic media.This method avoids the singular integrations in the conventional boundary element method,and can be treated as a potential approach for solving the inverse problems of the heat conduction owing to the boundary-only discretization and semi-analytical algorithm.When the VBEM is applied to the inverse problems,the numerical instability may occur if a virtual boundary is not properly chosen.The method encounters a highly illconditioned matrix for the larger distance between the physical boundary and the virtual boundary,and otherwise is hard to avoid the singularity of the source point.Thus,it must adopt an appropriate regularization method to deal with the ill-posed systems of inverse problems.In this study,the VBEM and different regularization techniques are combined to model the inverse problem of three-dimensional heat conduction in orthotropic media.The proper regularization techniques not only make the virtual boundary to be allocated freer,but also solve the ill-conditioned equation of the inverse problem.Numerical examples demonstrate that the proposed method is efficient,accurate and numerically stable for solving the inverse problems of three-dimensional heat conduction in orthotropic media.
基金supported by a CERG Grant of Hong Kong Research Grant Council,a FRG grant of Hong Kong Baptist University,and was partially supported by the NSFC Project No.19971116.
文摘Detecting corrosion by electrical field can be modeled by a Cauchy problem of Laplace equation in annulus domain under the assumption that the thickness of the pipe is relatively small compared with the radius of the pipe.The interior surface of the pipe is inaccessible and the nondestructive detection is solely based on measurements from the outer layer.The Cauchy problem for an elliptic equation is a typical ill-posed problem whose solution does not depend continuously on the boundary data.In this work,we assume that the measurements are available on the whole outer boundary on an annulus domain.By imposing reasonable assumptions,the theoretical goal here is to derive the stabilities of the Cauchy solutions and an energy regularization method.Relationship between the proposed energy regularization method and the Tikhonov regularization with Morozov principle is also given.A novel numerical algorithm is proposed and numerical examples are given.
文摘The accurate material physical properties, initial and boundary conditions are indispensable to the numerical simulation in the casting process, and they are related to the simulation accuracy directly. The inverse heat conduction method can be used to identify the mentioned above parameters based on the temperature measurement data. This paper presented a new inverse method according to Tikhonov regularization theory. A regularization functional was established and the regularization parameter was deduced, the Newton-Raphson iteration method was used to solve the equations. One detailed case was solved to identify the thermal conductivity and specific heat of sand mold and interfacial heat transfer coefficient (IHTC) at the meantime. This indicates that the regularization method is very efficient in decreasing the sensitivity to the temperature measurement data, overcoming the ill-posedness of the inverse heat conduction problem (IHCP) and improving the stability and accuracy of the results. As a general inverse method, it can be used to identify not only the material physical properties but also the initial and boundary conditions' parameters.
