In the past decade,notable progress has been achieved in the development of the generalized finite difference method(GFDM).The underlying principle of GFDM involves dividing the domain into multiple sub-domains.Within...In the past decade,notable progress has been achieved in the development of the generalized finite difference method(GFDM).The underlying principle of GFDM involves dividing the domain into multiple sub-domains.Within each sub-domain,explicit formulas for the necessary partial derivatives of the partial differential equations(PDEs)can be obtained through the application of Taylor series expansion and moving-least square approximation methods.Consequently,the method generates a sparse coefficient matrix,exhibiting a banded structure,making it highly advantageous for large-scale engineering computations.In this study,we present the application of the GFDM to numerically solve inverse Cauchy problems in two-and three-dimensional piezoelectric structures.Through our preliminary numerical experiments,we demonstrate that the proposed GFDMapproach shows great promise for accurately simulating coupled electroelastic equations in inverse problems,even with 3%errors added to the input data.展开更多
To solve the Laplacian problems,we adopt a meshless method with the multiquadric radial basis function(MQRBF)as a basis whose center is distributed inside a circle with a fictitious radius.A maximal projection techniq...To solve the Laplacian problems,we adopt a meshless method with the multiquadric radial basis function(MQRBF)as a basis whose center is distributed inside a circle with a fictitious radius.A maximal projection technique is developed to identify the optimal shape factor and fictitious radius by minimizing a merit function.A sample function is interpolated by theMQ-RBF to provide a trial coefficient vector to compute the merit function.We can quickly determine the optimal values of the parameters within a preferred rage using the golden section search algorithm.The novel method provides the optimal values of parameters and,hence,an optimal MQ-RBF;the performance of the method is validated in numerical examples.Moreover,nonharmonic problems are transformed to the Poisson equation endowed with a homogeneous boundary condition;this can overcome the problem of these problems being ill-posed.The optimal MQ-RBF is extremely accurate.We further propose a novel optimal polynomial method to solve the nonharmonic problems,which achieves high precision up to an order of 10^(−11).展开更多
In this paper, the inverse spectral problem of Sturm-Liouville operator with boundary conditions and jump conditions dependent on the spectral parameter is investigated. Firstly, the self-adjointness of the problem an...In this paper, the inverse spectral problem of Sturm-Liouville operator with boundary conditions and jump conditions dependent on the spectral parameter is investigated. Firstly, the self-adjointness of the problem and the eigenvalue properties are given, then the asymptotic formulas of eigenvalues and eigenfunctions are presented. Finally, the uniqueness theorems of the corresponding inverse problems are given by Weyl function theory and inverse spectral data approach.展开更多
This paper considers the estimation of an unknown function h that can be characterized as a solution to a nonlinear operator equation mapping between two infinite dimensional Hilbert spaces. The nonlinear operator is ...This paper considers the estimation of an unknown function h that can be characterized as a solution to a nonlinear operator equation mapping between two infinite dimensional Hilbert spaces. The nonlinear operator is unknown but can be consistently estimated, and its inverse is discontinuous, rendering the problem ill-posed. We establish the consistency for the class of estimators that are regularized using general lower semicompact penalty functions. We derive the optimal convergence rates of the estimators under the Hilbert scale norms. We apply our results to two important problems in economics and finance: (1) estimating the parameters of the pricing kernel of defaultable bonds; (2) recovering the volatility surface implied by option prices allowing for measurement error in the option prices and numerical error in the computation of the operator.展开更多
This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeabi...This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.展开更多
A novel method based on the relevance vector machine(RVM) for the inverse scattering problem is presented in this paper.The nonlinearity and the ill-posedness inherent in this problem are simultaneously considered.T...A novel method based on the relevance vector machine(RVM) for the inverse scattering problem is presented in this paper.The nonlinearity and the ill-posedness inherent in this problem are simultaneously considered.The nonlinearity is embodied in the relation between the scattered field and the target property,which can be obtained through the RVM training process.Besides,rather than utilizing regularization,the ill-posed nature of the inversion is naturally accounted for because the RVM can produce a probabilistic output.Simulation results reveal that the proposed RVM-based approach can provide comparative performances in terms of accuracy,convergence,robustness,generalization,and improved performance in terms of sparse property in comparison with the support vector machine(SVM) based approach.展开更多
Model error is one of the key factors restricting the accuracy of numerical weather prediction (NWP). Considering the continuous evolution of the atmosphere, the observed data (ignoring the measurement error) can ...Model error is one of the key factors restricting the accuracy of numerical weather prediction (NWP). Considering the continuous evolution of the atmosphere, the observed data (ignoring the measurement error) can be viewed as a series of solutions of an accurate model governing the actual atmosphere. Model error is represented as an unknown term in the accurate model, thus NWP can be considered as an inverse problem to uncover the unknown error term. The inverse problem models can absorb long periods of observed data to generate model error correction procedures. They thus resolve the deficiency and faultiness of the NWP schemes employing only the initial-time data. In this study we construct two inverse problem models to estimate and extrapolate the time-varying and spatial-varying model errors in both the historical and forecast periods by using recent observations and analogue phenomena of the atmosphere. Numerical experiment on Burgers' equation has illustrated the substantial forecast improvement using inverse problem algorithms. The proposed inverse problem methods of suppressing NWP errors will be useful in future high accuracy applications of NWP.展开更多
This paper presents an inverse problem in analytical dynamics. The inverse problem is to construct the Lagrangian when the integrals of a system are given. Firstly, the differential equations are obtained by using the...This paper presents an inverse problem in analytical dynamics. The inverse problem is to construct the Lagrangian when the integrals of a system are given. Firstly, the differential equations are obtained by using the time derivative of the integrals. Secondly, the differential equations can be written in the Lagrange equations under certain conditions and the Lagrangian can be obtained. Finally, two examples are given to illustrate the application of the result.展开更多
Errors inevitably exist in numerical weather prediction (NWP) due to imperfect numeric and physical parameterizations. To eliminate these errors, by considering NWP as an inverse problem, an unknown term in the pred...Errors inevitably exist in numerical weather prediction (NWP) due to imperfect numeric and physical parameterizations. To eliminate these errors, by considering NWP as an inverse problem, an unknown term in the prediction equations can be estimated inversely by using the past data, which are presumed to represent the imperfection of the NWP model (model error, denoted as ME). In this first paper of a two-part series, an iteration method for obtaining the MEs in past intervals is presented, and the results from testing its convergence in idealized experiments are reported. Moreover, two batches of iteration tests were applied in the global forecast system of the Global and Regional Assimilation and Prediction System (GRAPES-GFS) for July-August 2009 and January-February 2010. The datasets associated with the initial conditions and sea surface temperature (SST) were both based on NCEP (National Centers for Environmental Prediction) FNL (final) data. The results showed that 6th h forecast errors were reduced to 10% of their original value after a 20-step iteration. Then, off-line forecast error corrections were estimated linearly based on the 2-month mean MEs and compared with forecast errors. The estimated error corrections agreed well with the forecast errors, but the linear growth rate of the estimation was steeper than the forecast error. The advantage of this iteration method is that the MEs can provide the foundation for online correction. A larger proportion of the forecast errors can be expected to be canceled out by properly introducing the model error correction into GRAPES-GFS.展开更多
The spectral distribution exp( ), where {} are the eigenvalues of the negative Laplacian -△=- in the (x^1,x^2)-plane, is studied for a variety of domains, where -∞< t <∞ and i=(1/2)(-1) . The dependence of (...The spectral distribution exp( ), where {} are the eigenvalues of the negative Laplacian -△=- in the (x^1,x^2)-plane, is studied for a variety of domains, where -∞< t <∞ and i=(1/2)(-1) . The dependence of (t)on the connectivity of a domain and the boundary conditions are analyzed. Particular attention is given to a general bounded domain Ω in R^2 with a smooth boundary Ω, where a finite number of piecewise smooth Dirichlet, Neumann and Robin boundary conditions on the piecewise smooth parts Γj(j = 1,……,n) of Ω are considered such that Some geometrical properties of Ω(e.g., the area of Ω, the total lengths of the boundary, the curvature of its boundary, etc.) are determined, from the asymptotic expansions of (t) for |t| → 0.展开更多
This paper researches the following inverse eigenvalue problem for arrow-like matrices. Give two characteristic pairs, get a generalized arrow-like matrix, let the two characteristic pairs are the characteristic pairs...This paper researches the following inverse eigenvalue problem for arrow-like matrices. Give two characteristic pairs, get a generalized arrow-like matrix, let the two characteristic pairs are the characteristic pairs of this generalized arrow-like matrix. The expression and an algorithm of the solution of the problem is given, and a numerical example is provided.展开更多
Newton type methods are one kind of the efficient methods to solve nonlinear ill-posed problems, which have attracted extensive attention. However, computational cost of Newton type methods is high because practical p...Newton type methods are one kind of the efficient methods to solve nonlinear ill-posed problems, which have attracted extensive attention. However, computational cost of Newton type methods is high because practical problems are complicated. We propose a mixed Newton-Tikhonov method, i.e., one step Newton-Tikhonov method with several other steps of simplified Newton-Tikhonov method. Convergence and stability of this method are proved under some conditions. Numerical experiments show that the proposed method has obvious advantages over the classical Newton method in terms of computational costs.展开更多
In an indentation test,the effective Young's modulus of a film/substrate bilayer heterostructure varies with the indentation depth,a phenomenon known as the substrate effect.In previous studies investigating this,...In an indentation test,the effective Young's modulus of a film/substrate bilayer heterostructure varies with the indentation depth,a phenomenon known as the substrate effect.In previous studies investigating this,only the Young's modulus of the film was unknown.Once the effective Young's modulus of a film/substrate structure is determined at a given contact depth,the Young's modulus of the film can be uniquely determined,i.e.,there is a one-to-one relation between the Young's modulus of the film and the film/substrate effective Young's modulus.However,at times it is extremely challenging or even impossible to measure the film thickness.Furthermore,the precise definition of the layer/film thickness for a two-dimensional material can be problematic.In the current study,therefore,the thickness of the film and its Young's modulus are treated as two unknowns that must be determined.Unlike the case with one unknown,there are infinite combinations of film thickness and Young's modulus which can yield the same effective Young's modulus for the film/substrate.An inverse problem is formulated and solved to extract the Young's modulus and thickness of the film from the indentation depth-load curve.The accuracy and robustness of the inverse problem-solving method are also demonstrated.展开更多
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.展开更多
This article is devoted to the regularization of nonlinear ill-posed problems with accretive operators in Banach spaces. The data involved are assumed to be known approximately. The authors concentrate their discussio...This article is devoted to the regularization of nonlinear ill-posed problems with accretive operators in Banach spaces. The data involved are assumed to be known approximately. The authors concentrate their discussion on the convergence rates of regular solutions.展开更多
An inverse problem of elastica of a variable-arclength beam subjected to a concentrated load is investigated. The beam is fixed at one end, and can slide freely over a hinge support at the other end. The inverse probl...An inverse problem of elastica of a variable-arclength beam subjected to a concentrated load is investigated. The beam is fixed at one end, and can slide freely over a hinge support at the other end. The inverse problem is to determine the value of the load when the deflection of the action point of the load is given. Based on the elasitca equations and the elliptic integrals, a set of nonlinear equations for the inverse problem are derived, and an analytical solution by means of iterations and Quasi-Newton method is presented. From the results, the relationship between the loads and deflections of the loading point is obtained.展开更多
A completely new method is presented to solve the inverse black body radiation problem.The solution makes use of the widely known and easily implemented Fourier transform,and the expression is very simple and clear.Th...A completely new method is presented to solve the inverse black body radiation problem.The solution makes use of the widely known and easily implemented Fourier transform,and the expression is very simple and clear.This method is of major significance with respect not only to the numerical computation,but also to uniqueness or existence of this problem.At the same time,this method is also very useful for many other fields.展开更多
Aim The purpose of this study was to develop a mathe-matical model to quantitatively describe the passive trans-port of macromolecules within dental biofilms. Methodology Fluorescently labeled dextrans with different ...Aim The purpose of this study was to develop a mathe-matical model to quantitatively describe the passive trans-port of macromolecules within dental biofilms. Methodology Fluorescently labeled dextrans with different molecular mass (3 kD,10 kD,40 kD,70 kD,2 000 kD) were used as a series of diffusion probes. Streptococcus mutans,Streptococcus sanguinis,Actinomyces naeslundii and Fusobacterium nucleatum were used as inocula for biofilm formation. The diffusion processes of different probes through the in vitro biofilm were recorded with a confocal laser microscope. Results Mathematical function of biofilm penetration was constructed on the basis of the inverse problem method. Based on this function,not only the relationship between average concentration of steady-state and molecule weights can be analyzed,but also that between penetrative time and molecule weights. Conclusion This can be used to predict the effective concentration and the penetrative time of anti-biofilm medicines that can diffuse through oral biofilm. Further-more,an improved model for large molecule is proposed by considering the exchange time at the upper boundary of the dental biofilm.展开更多
The computational complexity of inverse mimimum capacity path problem with lower bound on capacity of maximum capacity path is examined, and it is proved that solution of this problem is NP-complete. A strong polynomi...The computational complexity of inverse mimimum capacity path problem with lower bound on capacity of maximum capacity path is examined, and it is proved that solution of this problem is NP-complete. A strong polynomial algorithm for a local optimal solution is provided.展开更多
This paper studies the influence of a finite container on an ideal gas.The trace of the heat kernel (t) =exp, where are the eigenvalues of the negative Laplacian -in Rn(n = 2 or 3), is studied for a general multi-conn...This paper studies the influence of a finite container on an ideal gas.The trace of the heat kernel (t) =exp, where are the eigenvalues of the negative Laplacian -in Rn(n = 2 or 3), is studied for a general multi-connected bounded drum ft which is surrounded by simply connected bounded domains Ωi with smooth boundaries Ωi(i = 1,… ,m) where the Dirichlet, Neumann and Robin boundary conditions on Ωi(i = 1,…,m) are considered. Some geometrical properties of Ω are determined. The thermodynamic quantities for an ideal gas enclosed in Ω are examined by using the asymptotic expansions of (t) for short-time t. It is shown that the ideal gas can not feel the shape of its container Ω, although it can feel some geometrical properties of it.展开更多
基金the Natural Science Foundation of Shandong Province of China(Grant No.ZR2022YQ06)the Development Plan of Youth Innovation Team in Colleges and Universities of Shandong Province(Grant No.2022KJ140)the Key Laboratory ofRoad Construction Technology and Equipment(Chang’an University,No.300102253502).
文摘In the past decade,notable progress has been achieved in the development of the generalized finite difference method(GFDM).The underlying principle of GFDM involves dividing the domain into multiple sub-domains.Within each sub-domain,explicit formulas for the necessary partial derivatives of the partial differential equations(PDEs)can be obtained through the application of Taylor series expansion and moving-least square approximation methods.Consequently,the method generates a sparse coefficient matrix,exhibiting a banded structure,making it highly advantageous for large-scale engineering computations.In this study,we present the application of the GFDM to numerically solve inverse Cauchy problems in two-and three-dimensional piezoelectric structures.Through our preliminary numerical experiments,we demonstrate that the proposed GFDMapproach shows great promise for accurately simulating coupled electroelastic equations in inverse problems,even with 3%errors added to the input data.
基金supported by the the National Science and Technology Council(Grant Number:NSTC 112-2221-E239-022).
文摘To solve the Laplacian problems,we adopt a meshless method with the multiquadric radial basis function(MQRBF)as a basis whose center is distributed inside a circle with a fictitious radius.A maximal projection technique is developed to identify the optimal shape factor and fictitious radius by minimizing a merit function.A sample function is interpolated by theMQ-RBF to provide a trial coefficient vector to compute the merit function.We can quickly determine the optimal values of the parameters within a preferred rage using the golden section search algorithm.The novel method provides the optimal values of parameters and,hence,an optimal MQ-RBF;the performance of the method is validated in numerical examples.Moreover,nonharmonic problems are transformed to the Poisson equation endowed with a homogeneous boundary condition;this can overcome the problem of these problems being ill-posed.The optimal MQ-RBF is extremely accurate.We further propose a novel optimal polynomial method to solve the nonharmonic problems,which achieves high precision up to an order of 10^(−11).
文摘In this paper, the inverse spectral problem of Sturm-Liouville operator with boundary conditions and jump conditions dependent on the spectral parameter is investigated. Firstly, the self-adjointness of the problem and the eigenvalue properties are given, then the asymptotic formulas of eigenvalues and eigenfunctions are presented. Finally, the uniqueness theorems of the corresponding inverse problems are given by Weyl function theory and inverse spectral data approach.
基金supported by US National Science Foundation (Grant No. SES-0631613)the Cowles Foundation for Research in Economics
文摘This paper considers the estimation of an unknown function h that can be characterized as a solution to a nonlinear operator equation mapping between two infinite dimensional Hilbert spaces. The nonlinear operator is unknown but can be consistently estimated, and its inverse is discontinuous, rendering the problem ill-posed. We establish the consistency for the class of estimators that are regularized using general lower semicompact penalty functions. We derive the optimal convergence rates of the estimators under the Hilbert scale norms. We apply our results to two important problems in economics and finance: (1) estimating the parameters of the pricing kernel of defaultable bonds; (2) recovering the volatility surface implied by option prices allowing for measurement error in the option prices and numerical error in the computation of the operator.
文摘This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.
基金Project supported by the National Natural Science Foundation of China (Grant No. 61071022)the Graduate Student Research and Innovation Program of Jiangsu Province,China (Grant No. CXZZ11-0381)
文摘A novel method based on the relevance vector machine(RVM) for the inverse scattering problem is presented in this paper.The nonlinearity and the ill-posedness inherent in this problem are simultaneously considered.The nonlinearity is embodied in the relation between the scattered field and the target property,which can be obtained through the RVM training process.Besides,rather than utilizing regularization,the ill-posed nature of the inversion is naturally accounted for because the RVM can produce a probabilistic output.Simulation results reveal that the proposed RVM-based approach can provide comparative performances in terms of accuracy,convergence,robustness,generalization,and improved performance in terms of sparse property in comparison with the support vector machine(SVM) based approach.
基金Project supported by the Special Scientific Research Project for Public Interest(Grant No.GYHY201206009)the Fundamental Research Funds for the Central Universities,China(Grant Nos.lzujbky-2012-13 and lzujbky-2013-11)the National Basic Research Program of China(Grant Nos.2012CB955902 and 2013CB430204)
文摘Model error is one of the key factors restricting the accuracy of numerical weather prediction (NWP). Considering the continuous evolution of the atmosphere, the observed data (ignoring the measurement error) can be viewed as a series of solutions of an accurate model governing the actual atmosphere. Model error is represented as an unknown term in the accurate model, thus NWP can be considered as an inverse problem to uncover the unknown error term. The inverse problem models can absorb long periods of observed data to generate model error correction procedures. They thus resolve the deficiency and faultiness of the NWP schemes employing only the initial-time data. In this study we construct two inverse problem models to estimate and extrapolate the time-varying and spatial-varying model errors in both the historical and forecast periods by using recent observations and analogue phenomena of the atmosphere. Numerical experiment on Burgers' equation has illustrated the substantial forecast improvement using inverse problem algorithms. The proposed inverse problem methods of suppressing NWP errors will be useful in future high accuracy applications of NWP.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 10272021, 10572021) and the Doctoral Programme Foundation of Institution of Higher Education of China (Grant No 20040007022).
文摘This paper presents an inverse problem in analytical dynamics. The inverse problem is to construct the Lagrangian when the integrals of a system are given. Firstly, the differential equations are obtained by using the time derivative of the integrals. Secondly, the differential equations can be written in the Lagrange equations under certain conditions and the Lagrangian can be obtained. Finally, two examples are given to illustrate the application of the result.
基金funded by the National Natural Science Foundation Science Fund for Youth (Grant No.41405095)the Key Projects in the National Science and Technology Pillar Program during the Twelfth Fiveyear Plan Period (Grant No.2012BAC22B02)the National Natural Science Foundation Science Fund for Creative Research Groups (Grant No.41221064)
文摘Errors inevitably exist in numerical weather prediction (NWP) due to imperfect numeric and physical parameterizations. To eliminate these errors, by considering NWP as an inverse problem, an unknown term in the prediction equations can be estimated inversely by using the past data, which are presumed to represent the imperfection of the NWP model (model error, denoted as ME). In this first paper of a two-part series, an iteration method for obtaining the MEs in past intervals is presented, and the results from testing its convergence in idealized experiments are reported. Moreover, two batches of iteration tests were applied in the global forecast system of the Global and Regional Assimilation and Prediction System (GRAPES-GFS) for July-August 2009 and January-February 2010. The datasets associated with the initial conditions and sea surface temperature (SST) were both based on NCEP (National Centers for Environmental Prediction) FNL (final) data. The results showed that 6th h forecast errors were reduced to 10% of their original value after a 20-step iteration. Then, off-line forecast error corrections were estimated linearly based on the 2-month mean MEs and compared with forecast errors. The estimated error corrections agreed well with the forecast errors, but the linear growth rate of the estimation was steeper than the forecast error. The advantage of this iteration method is that the MEs can provide the foundation for online correction. A larger proportion of the forecast errors can be expected to be canceled out by properly introducing the model error correction into GRAPES-GFS.
文摘The spectral distribution exp( ), where {} are the eigenvalues of the negative Laplacian -△=- in the (x^1,x^2)-plane, is studied for a variety of domains, where -∞< t <∞ and i=(1/2)(-1) . The dependence of (t)on the connectivity of a domain and the boundary conditions are analyzed. Particular attention is given to a general bounded domain Ω in R^2 with a smooth boundary Ω, where a finite number of piecewise smooth Dirichlet, Neumann and Robin boundary conditions on the piecewise smooth parts Γj(j = 1,……,n) of Ω are considered such that Some geometrical properties of Ω(e.g., the area of Ω, the total lengths of the boundary, the curvature of its boundary, etc.) are determined, from the asymptotic expansions of (t) for |t| → 0.
文摘This paper researches the following inverse eigenvalue problem for arrow-like matrices. Give two characteristic pairs, get a generalized arrow-like matrix, let the two characteristic pairs are the characteristic pairs of this generalized arrow-like matrix. The expression and an algorithm of the solution of the problem is given, and a numerical example is provided.
基金supported by the Key Disciplines of Shanghai Municipality (Operations Research & Cybernetics, No. S30104)Shanghai Leading Academic Discipline Project (No. J50101)
文摘Newton type methods are one kind of the efficient methods to solve nonlinear ill-posed problems, which have attracted extensive attention. However, computational cost of Newton type methods is high because practical problems are complicated. We propose a mixed Newton-Tikhonov method, i.e., one step Newton-Tikhonov method with several other steps of simplified Newton-Tikhonov method. Convergence and stability of this method are proved under some conditions. Numerical experiments show that the proposed method has obvious advantages over the classical Newton method in terms of computational costs.
基金This research was supported by the National Natural Science Foundation of China(11772335,21622304,61674045,and 21203038)by the Ministry of Science and Technology(MOST)of China(2016YFA0200700).Z.H.Cheng was supported by the Distinguished Technical Talents Project and the Youth Innovation Promotion Association of Chinese Academy of Sciences.
文摘In an indentation test,the effective Young's modulus of a film/substrate bilayer heterostructure varies with the indentation depth,a phenomenon known as the substrate effect.In previous studies investigating this,only the Young's modulus of the film was unknown.Once the effective Young's modulus of a film/substrate structure is determined at a given contact depth,the Young's modulus of the film can be uniquely determined,i.e.,there is a one-to-one relation between the Young's modulus of the film and the film/substrate effective Young's modulus.However,at times it is extremely challenging or even impossible to measure the film thickness.Furthermore,the precise definition of the layer/film thickness for a two-dimensional material can be problematic.In the current study,therefore,the thickness of the film and its Young's modulus are treated as two unknowns that must be determined.Unlike the case with one unknown,there are infinite combinations of film thickness and Young's modulus which can yield the same effective Young's modulus for the film/substrate.An inverse problem is formulated and solved to extract the Young's modulus and thickness of the film from the indentation depth-load curve.The accuracy and robustness of the inverse problem-solving method are also demonstrated.
基金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.
文摘This article is devoted to the regularization of nonlinear ill-posed problems with accretive operators in Banach spaces. The data involved are assumed to be known approximately. The authors concentrate their discussion on the convergence rates of regular solutions.
基金The project supported by the National Natural Science Foundation of China(10272011)
文摘An inverse problem of elastica of a variable-arclength beam subjected to a concentrated load is investigated. The beam is fixed at one end, and can slide freely over a hinge support at the other end. The inverse problem is to determine the value of the load when the deflection of the action point of the load is given. Based on the elasitca equations and the elliptic integrals, a set of nonlinear equations for the inverse problem are derived, and an analytical solution by means of iterations and Quasi-Newton method is presented. From the results, the relationship between the loads and deflections of the loading point is obtained.
文摘A completely new method is presented to solve the inverse black body radiation problem.The solution makes use of the widely known and easily implemented Fourier transform,and the expression is very simple and clear.This method is of major significance with respect not only to the numerical computation,but also to uniqueness or existence of this problem.At the same time,this method is also very useful for many other fields.
基金supported by a grant from the National Natural Science Foundation of China (NSFC) No. 81070826/30872886/30400497Sponsored by Shanghai Rising-Star Program No. 09QA1403700+1 种基金funded by Shanghai Leading Academic Discipline Project (Project Number: S30206)the Science and Technology Commission of Shanghai (08DZ2271100)
文摘Aim The purpose of this study was to develop a mathe-matical model to quantitatively describe the passive trans-port of macromolecules within dental biofilms. Methodology Fluorescently labeled dextrans with different molecular mass (3 kD,10 kD,40 kD,70 kD,2 000 kD) were used as a series of diffusion probes. Streptococcus mutans,Streptococcus sanguinis,Actinomyces naeslundii and Fusobacterium nucleatum were used as inocula for biofilm formation. The diffusion processes of different probes through the in vitro biofilm were recorded with a confocal laser microscope. Results Mathematical function of biofilm penetration was constructed on the basis of the inverse problem method. Based on this function,not only the relationship between average concentration of steady-state and molecule weights can be analyzed,but also that between penetrative time and molecule weights. Conclusion This can be used to predict the effective concentration and the penetrative time of anti-biofilm medicines that can diffuse through oral biofilm. Further-more,an improved model for large molecule is proposed by considering the exchange time at the upper boundary of the dental biofilm.
基金The authors gratefully acknowledge the partial support of national natural Founda-tion (Grant 70071011)
文摘The computational complexity of inverse mimimum capacity path problem with lower bound on capacity of maximum capacity path is examined, and it is proved that solution of this problem is NP-complete. A strong polynomial algorithm for a local optimal solution is provided.
文摘This paper studies the influence of a finite container on an ideal gas.The trace of the heat kernel (t) =exp, where are the eigenvalues of the negative Laplacian -in Rn(n = 2 or 3), is studied for a general multi-connected bounded drum ft which is surrounded by simply connected bounded domains Ωi with smooth boundaries Ωi(i = 1,… ,m) where the Dirichlet, Neumann and Robin boundary conditions on Ωi(i = 1,…,m) are considered. Some geometrical properties of Ω are determined. The thermodynamic quantities for an ideal gas enclosed in Ω are examined by using the asymptotic expansions of (t) for short-time t. It is shown that the ideal gas can not feel the shape of its container Ω, although it can feel some geometrical properties of it.