Sensitivity of observational data is important in the study of Glacial Isostatic Adjustment(GIA).However,depending on whether sensitivity is used for the Inverse Problem or the Forward Problem,the final formulation an...Sensitivity of observational data is important in the study of Glacial Isostatic Adjustment(GIA).However,depending on whether sensitivity is used for the Inverse Problem or the Forward Problem,the final formulation and display of the sensitivity kernel will be different.Unfortunately,in the past,both perspectives give the same name to their quantity computed/displayed,and that has caused some confusion.To distinguish between the two,their perspective should be added to the names.This paper focuses only on the perspective of the Forward Problem where the input parameters are known.The Perturbation method has been successfully used in the computation of the sensitivity kernels of observations on 1D and 3D viscosity variations from the Forward perspective.One aim of this paper is to review and clarify the physics of the Perturbation method and bring out some important aspects of this method that have been misunderstood or neglected.Another aim is to present sensitivity kernels from the Perturbation method using 3D(both radially and laterally heterogeneous)Earth models with realistic ice history.These new results are now suitable for future comparison with those from new methods using the Forward perspective.Finally,the sensitivity computations for realistic ice histories on a 3D Earth is reviewed and used to search for optimal locations of new GIA observations.展开更多
In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow e...In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow equation.The velocity and pressure are computed simultaneously.The accuracy of velocity is improved one order.The concentration equation is solved by using mixed finite element,multi-step difference and upwind approximation.A multi-step method is used to approximate time derivative for improving the accuracy.The upwind approximation and an expanded mixed finite element are adopted to solve the convection and diffusion,respectively.The composite method could compute the diffusion flux and its gradient.It possibly becomes an eficient tool for solving convection-dominated diffusion problems.Firstly,the conservation of mass holds.Secondly,the multi-step method has high accuracy.Thirdly,the upwind approximation could avoid numerical dispersion.Using numerical analysis of a priori estimates and special techniques of differential equations,we give an error estimates for a positive definite problem.Numerical experiments illustrate its computational efficiency and feasibility of application.展开更多
Two kinds of contact problems, i.e., the frictional contact problem and the adhesive contact problem, in three-dimensional （3D） icosahedral quasicrystals are dis- cussed by a complex variable function method. For th...Two kinds of contact problems, i.e., the frictional contact problem and the adhesive contact problem, in three-dimensional （3D） icosahedral quasicrystals are dis- cussed by a complex variable function method. For the frictional contact problem, the contact stress exhibits power singularities at the edge of the contact zone. For the adhe- sive contact problem, the contact stress exhibits oscillatory singularities at the edge of the contact zone. The numerical examples show that for the two kinds of contact problems, the contact stress exhibits singularities, and reaches the maximum value at the edge of the contact zone. The phonon-phason coupling constant has a significant effect on the contact stress intensity, while has little impact on the contact stress distribution regu- lation. The results are consistent with those of the classical elastic materials when the phonon-phason coupling constant is 0. For the adhesive contact problem, the indentation force has positive correlation with the contact displacement, but the phonon-phason cou- pling constant impact is barely perceptible. The validity of the conclusions is verified.展开更多
A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forwar...A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.展开更多
A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been deve...A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.展开更多
We present an efficient three-dimensional coupled-mode model based on the Fourier synthesis technique. In principle, this model is a one-way model, and hence provides satisfactory accuracy for problems where the forwa...We present an efficient three-dimensional coupled-mode model based on the Fourier synthesis technique. In principle, this model is a one-way model, and hence provides satisfactory accuracy for problems where the forward scattering dominates. At the same time, this model provides an efficiency gain of an order of magnitude or more over two-way coupled-mode models. This model can be applied to three-dimensional range-dependent problems with a slowly varying bathymetry or internal waves. A numerical example of the latter is demonstrated in this work. Comparisons of both accuracy and efficiency between the present model and a benchmark model are also provided.展开更多
In this paper, we use differential game theory to study the three-dimensional two-aircraft air-to-air combat problem. We give the ways to determine the Capture Ranges (CR) and the Dangerous Ranges (DR) for these two a...In this paper, we use differential game theory to study the three-dimensional two-aircraft air-to-air combat problem. We give the ways to determine the Capture Ranges (CR) and the Dangerous Ranges (DR) for these two aircraft according to the target entry directions, barrier and isochronic lines respectively. The simulations are given by referring to two sets of real aircraft parameters. After discussing the simulation results, we have obtained some conclusions that match the real air-to-air combat situation quite well.展开更多
Analytically solving a three-dimensional （3-D） bioheat transfer problem with phase change during a freezing process is extremely difficult but theoretically important. The moving heat source model and the Green func...Analytically solving a three-dimensional （3-D） bioheat transfer problem with phase change during a freezing process is extremely difficult but theoretically important. The moving heat source model and the Green function method are introduced to deal with the cryopreservation process of in vitro biomaterials. Exact solutions for the 3-D temperature transients of tissues under various boundary conditions, such as totally convective cooling, totally fixed temperature cooling and a hybrid between them on tissue surfaces, are obtained. Furthermore, the cryosurgical process in living tissues subject to freezing by a single or multiple cryoprobes is also analytically solved. A closed-form analytical solution to the bioheat phase change process is derived by considering contributions from blood perfusion heat transfer, metabolic heat generation, and heat sink of a cryoprobe. The present method is expected to have significant value for analytically solving complex bioheat transfer problems with phase change.展开更多
The paper discusses the application of three-dimensional textbook in English teaching. The three-dimensional textbook plays an incomparable function than the traditional textbook and thus has an extensive application ...The paper discusses the application of three-dimensional textbook in English teaching. The three-dimensional textbook plays an incomparable function than the traditional textbook and thus has an extensive application in the college English education. It demonstrates the application through the usage of the three-dimensional textbook in the real teaching and shows that excessive usage of sound and pictures or drawings may distract the attention of students. It raises several suggestions at the end.展开更多
The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the tim...The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.展开更多
The interaction of arbitrarily distributed penny-shaped cracks in three-dimensional solids is analyzed in this paper. Using oblate spheroidal coordinates and displacement functions, an analytic method is devel- oped i...The interaction of arbitrarily distributed penny-shaped cracks in three-dimensional solids is analyzed in this paper. Using oblate spheroidal coordinates and displacement functions, an analytic method is devel- oped in which the opening and the sliding displacements on each crack surface are taken as the basic unknown functions. The basic unknown functions can be expanded in series of Legendre polynomials with unknown coefficients. Based on superposition technique, a set of governing equations for the unknown coefficients are formulated from the traction free conditions on each crack surface. The boundary collocation procedure and the average method for crack-surface tractions are used for solving the governing equations. The solution can be obtained for quite closely located cracks. Numerical examples are given for several crack problems. By comparing the present results with other existing results, one can conclude that the present method provides a direct and efficient approach to deal with three-dimensional solids containing multiple cracks.展开更多
Transonic shocks play a pivotal role in designation of supersonic inlets and ramjets.For the three-dimensional steady non-isentropic compressible Euler system with frictions,we constructe a family of transonic shock s...Transonic shocks play a pivotal role in designation of supersonic inlets and ramjets.For the three-dimensional steady non-isentropic compressible Euler system with frictions,we constructe a family of transonic shock solutions in rectilinear ducts with square cross-sections.In this article,we are devoted to proving rigorously that a large class of these transonic shock solutions are stable,under multidimensional small perturbations of the upcoming supersonic flows and back pressures at the exits of ducts in suitable function spaces.This manifests that frictions have a stabilization effect on transonic shocks in ducts,in consideration of previous works which shown that transonic shocks in purely steady Euler flows are not stable in such ducts.Except its implications to applications,because frictions lead to a stronger coupling between the elliptic and hyperbolic parts of the three-dimensional steady subsonic Euler system,we develop the framework established in previous works to study more complex and interesting Venttsel problems of nonlocal elliptic equations.展开更多
Petroleum science has made remarkable progress in organic geochcmistry and in the research into the theories of petroleum origin, its transport and accumulation. In estimating the oil-gas resources of a basin, the kno...Petroleum science has made remarkable progress in organic geochcmistry and in the research into the theories of petroleum origin, its transport and accumulation. In estimating the oil-gas resources of a basin, the knowledge of its evolutionary history and especially the numerical computation of fluid flow and the history of its changes under heat is vital. The mathematical model can be described as a coupled system of nonlinear partial differentical equations with initial-boundary value problems. This thesis, from actual conditions such as the effect of fluid compressibility and the three-dimensional characteristic of large-scale science-engineering computation, we put forward a kind of characteristic finite element alternating-direction schemes and obtain optimal order estimates in L^2 norm for the error in the approximate assumption.展开更多
In this paper. an analytic method is. presented. for the research of nonlinear Three-dimensional problems of composite laminated plates. The perturbation method and the variational principle are used to satisfy the ba...In this paper. an analytic method is. presented. for the research of nonlinear Three-dimensional problems of composite laminated plates. The perturbation method and the variational principle are used to satisfy the basic differential equations and the boundary conditions of the three-dimensional theory of elasticity. The nonlinear three-dimensional problems are studied .for composite anisotropic circular laminas and laminates subjected to transverse loading. The perturbation series solutions of high accuracy are obtained. A large number of results show that transverse normal stress and transverse shear stresses are very important in the nonlinear three-dimensional analysis of laminated plates.展开更多
This paper deals with the forward and backward problems for the nonlinear fractional pseudo-parabolic equation ut+(-Δ)^(s1)ut+β(-Δ)^(s2)u=F(u,x,t)subject o random Gaussian white noise for initial and final data.Und...This paper deals with the forward and backward problems for the nonlinear fractional pseudo-parabolic equation ut+(-Δ)^(s1)ut+β(-Δ)^(s2)u=F(u,x,t)subject o random Gaussian white noise for initial and final data.Under the suitable assumptions s1,s2andβ,we first show the ill-posedness of mild solutions for forward and backward problems in the sense of Hadamard,which are mainly driven by random noise.Moreover,we propose the Fourier truncation method for stabilizing the above ill-posed problems.We derive an error estimate between the exact solution and its regularized solution in an E‖·‖Hs22norm,and give some numerical examples illustrating the effect of above method.展开更多
The Magneto-acoustic Tomography with Current Injection (MAT-CI) is a new biological electrical impedance imaging technique that combines Electrical Impedance Tomography (EIT) with Ultrasonic Imaging (UI), which posses...The Magneto-acoustic Tomography with Current Injection (MAT-CI) is a new biological electrical impedance imaging technique that combines Electrical Impedance Tomography (EIT) with Ultrasonic Imaging (UI), which possesses the non-invasive and high-contrast of the EIT and the high-resolution of the UI. The MAT-CI is expected to acquire high quality image and embraces a wide application. Its principle is to put the conductive sample in the Static Magnetic Field(SMF) and inject a time-varying current, during which the SMF and the current interact and generate the Lorentz Force that inspire ultrasonic signal received by the ultrasonic transducers positioned around the sample. And then according to related reconstruction algorithm and ultrasonic signal, electrical conductivity image is obtained. In this paper, a forward problem mathematical model of the MAT-CI has been set up to deduce the theoretical equation of the electromagnetic field and solve the sound source distribution by Green’s function. Secondly, a sound field restoration by Wiener filtering and reconstruction of current density by time-rotating method have deduced the Laplace’s equation that caters to the current density to further acquire the electrical conductivity distribution image of the sample through iteration method. In the end, double-loop coils experiments have been conducted to verify its feasibility.展开更多
By making use of the direct integration method,an exact analysis of the general three-dimensional thermoelasticity problem is performed for the case of a transversely isotropic homogeneous half-space subject to local ...By making use of the direct integration method,an exact analysis of the general three-dimensional thermoelasticity problem is performed for the case of a transversely isotropic homogeneous half-space subject to local thermal and force loadings.The material plane of isotropy is assumed to be parallel to the limiting surface of the halfspace.By reducing the original thermoelasticity equations to the governing ones for individual stress-tensor components,the effect of material anisotropy in the stress field is analyzed with regard to the feasibility requirement,i.e.,the finiteness of the stress field at a distance from the disturbed area.As a result,the solution is constructed in the form of explicit analytical dependencies on the force and thermal loadings for various kinds of transversely isotropic materials and agrees with the basic principles of the continua mechanics.The solution can be efficiently used as a benchmark one for the direct computation of temperature and thermal stresses in transversely isotropic semi-infinite domains,as well as for the verification of solutions constructed by different means.展开更多
Physics-informed neural networks(PINNs),as a novel artificial intelligence method for solving partial differential equations,are applicable to solve both forward and inverse problems.This study evaluates the performan...Physics-informed neural networks(PINNs),as a novel artificial intelligence method for solving partial differential equations,are applicable to solve both forward and inverse problems.This study evaluates the performance of PINNs in solving the temperature diffusion equation of the seawater across six scenarios,including forward and inverse problems under three different boundary conditions.Results demonstrate that PINNs achieved consistently higher accuracy with the Dirichlet and Neumann boundary conditions compared to the Robin boundary condition for both forward and inverse problems.Inaccurate weighting of terms in the loss function can reduce model accuracy.Additionally,the sensitivity of model performance to the positioning of sampling points varied between different boundary conditions.In particular,the model under the Dirichlet boundary condition exhibited superior robustness to variations in point positions during the solutions of inverse problems.In contrast,for the Neumann and Robin boundary conditions,accuracy declines when points were sampled from identical positions or at the same time.Subsequently,the Argo observations were used to reconstruct the vertical diffusion of seawater temperature in the north-central Pacific for the applicability of PINNs in the real ocean.The PINNs successfully captured the vertical diffusion characteristics of seawater temperature,reflected the seasonal changes of vertical temperature under different topographic conditions,and revealed the influence of topography on the temperature diffusion coefficient.The PINNs were proved effective in solving the temperature diffusion equation of seawater with limited data,providing a promising technique for simulating or predicting ocean phenomena using sparse observations.展开更多
文摘Sensitivity of observational data is important in the study of Glacial Isostatic Adjustment(GIA).However,depending on whether sensitivity is used for the Inverse Problem or the Forward Problem,the final formulation and display of the sensitivity kernel will be different.Unfortunately,in the past,both perspectives give the same name to their quantity computed/displayed,and that has caused some confusion.To distinguish between the two,their perspective should be added to the names.This paper focuses only on the perspective of the Forward Problem where the input parameters are known.The Perturbation method has been successfully used in the computation of the sensitivity kernels of observations on 1D and 3D viscosity variations from the Forward perspective.One aim of this paper is to review and clarify the physics of the Perturbation method and bring out some important aspects of this method that have been misunderstood or neglected.Another aim is to present sensitivity kernels from the Perturbation method using 3D(both radially and laterally heterogeneous)Earth models with realistic ice history.These new results are now suitable for future comparison with those from new methods using the Forward perspective.Finally,the sensitivity computations for realistic ice histories on a 3D Earth is reviewed and used to search for optimal locations of new GIA observations.
基金supported by the Natural Science Foundation of Shandong Province(ZR2021MA019)the National Natural Science Foundation of China(11871312)。
文摘In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow equation.The velocity and pressure are computed simultaneously.The accuracy of velocity is improved one order.The concentration equation is solved by using mixed finite element,multi-step difference and upwind approximation.A multi-step method is used to approximate time derivative for improving the accuracy.The upwind approximation and an expanded mixed finite element are adopted to solve the convection and diffusion,respectively.The composite method could compute the diffusion flux and its gradient.It possibly becomes an eficient tool for solving convection-dominated diffusion problems.Firstly,the conservation of mass holds.Secondly,the multi-step method has high accuracy.Thirdly,the upwind approximation could avoid numerical dispersion.Using numerical analysis of a priori estimates and special techniques of differential equations,we give an error estimates for a positive definite problem.Numerical experiments illustrate its computational efficiency and feasibility of application.
基金supported by the National Natural Science Foundation of China(Nos.11362018,11261045,and 11261401)the Specialized Research Fund for the Doctoral Program of Higher Education of China(No.20116401110002)
文摘Two kinds of contact problems, i.e., the frictional contact problem and the adhesive contact problem, in three-dimensional （3D） icosahedral quasicrystals are dis- cussed by a complex variable function method. For the frictional contact problem, the contact stress exhibits power singularities at the edge of the contact zone. For the adhe- sive contact problem, the contact stress exhibits oscillatory singularities at the edge of the contact zone. The numerical examples show that for the two kinds of contact problems, the contact stress exhibits singularities, and reaches the maximum value at the edge of the contact zone. The phonon-phason coupling constant has a significant effect on the contact stress intensity, while has little impact on the contact stress distribution regu- lation. The results are consistent with those of the classical elastic materials when the phonon-phason coupling constant is 0. For the adhesive contact problem, the indentation force has positive correlation with the contact displacement, but the phonon-phason cou- pling constant impact is barely perceptible. The validity of the conclusions is verified.
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.
基金Projects(41674080,41674079)supported by the National Natural Science Foundation of China
文摘A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.
基金Supported by the National Natural Science Foundation of China under Grant No 11774374the Natural Science Foundation of Shandong Province of China under Grant No ZR2016AL10
文摘We present an efficient three-dimensional coupled-mode model based on the Fourier synthesis technique. In principle, this model is a one-way model, and hence provides satisfactory accuracy for problems where the forward scattering dominates. At the same time, this model provides an efficiency gain of an order of magnitude or more over two-way coupled-mode models. This model can be applied to three-dimensional range-dependent problems with a slowly varying bathymetry or internal waves. A numerical example of the latter is demonstrated in this work. Comparisons of both accuracy and efficiency between the present model and a benchmark model are also provided.
基金research was supported by Aviation Science Fund.
文摘In this paper, we use differential game theory to study the three-dimensional two-aircraft air-to-air combat problem. We give the ways to determine the Capture Ranges (CR) and the Dangerous Ranges (DR) for these two aircraft according to the target entry directions, barrier and isochronic lines respectively. The simulations are given by referring to two sets of real aircraft parameters. After discussing the simulation results, we have obtained some conclusions that match the real air-to-air combat situation quite well.
基金Project supported by the National Natural Science Foundation of China (No. 50776097)
文摘Analytically solving a three-dimensional （3-D） bioheat transfer problem with phase change during a freezing process is extremely difficult but theoretically important. The moving heat source model and the Green function method are introduced to deal with the cryopreservation process of in vitro biomaterials. Exact solutions for the 3-D temperature transients of tissues under various boundary conditions, such as totally convective cooling, totally fixed temperature cooling and a hybrid between them on tissue surfaces, are obtained. Furthermore, the cryosurgical process in living tissues subject to freezing by a single or multiple cryoprobes is also analytically solved. A closed-form analytical solution to the bioheat phase change process is derived by considering contributions from blood perfusion heat transfer, metabolic heat generation, and heat sink of a cryoprobe. The present method is expected to have significant value for analytically solving complex bioheat transfer problems with phase change.
文摘The paper discusses the application of three-dimensional textbook in English teaching. The three-dimensional textbook plays an incomparable function than the traditional textbook and thus has an extensive application in the college English education. It demonstrates the application through the usage of the three-dimensional textbook in the real teaching and shows that excessive usage of sound and pictures or drawings may distract the attention of students. It raises several suggestions at the end.
基金supported by the Joint Fund of Seismological Science(Grant No.U1839206)the National R&D Program on Monitoring,Early Warning and Prevention of Major Natural Disaster(Grant No.2017YFC1500301)+2 种基金supported by IGGCAS Research Start-up Funds(Grant No.E0515402)National Natural Science Foundation of China(Grant No.E1115401)supported by National Natural Science Foundation of China(Grant No.11971258).
文摘The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.
文摘The interaction of arbitrarily distributed penny-shaped cracks in three-dimensional solids is analyzed in this paper. Using oblate spheroidal coordinates and displacement functions, an analytic method is devel- oped in which the opening and the sliding displacements on each crack surface are taken as the basic unknown functions. The basic unknown functions can be expanded in series of Legendre polynomials with unknown coefficients. Based on superposition technique, a set of governing equations for the unknown coefficients are formulated from the traction free conditions on each crack surface. The boundary collocation procedure and the average method for crack-surface tractions are used for solving the governing equations. The solution can be obtained for quite closely located cracks. Numerical examples are given for several crack problems. By comparing the present results with other existing results, one can conclude that the present method provides a direct and efficient approach to deal with three-dimensional solids containing multiple cracks.
基金This work was supported in part by National Nature Science Foundation of China(11371141 and 11871218)by Science and Technology Commission of Shanghai Municipality(18dz2271000).
文摘Transonic shocks play a pivotal role in designation of supersonic inlets and ramjets.For the three-dimensional steady non-isentropic compressible Euler system with frictions,we constructe a family of transonic shock solutions in rectilinear ducts with square cross-sections.In this article,we are devoted to proving rigorously that a large class of these transonic shock solutions are stable,under multidimensional small perturbations of the upcoming supersonic flows and back pressures at the exits of ducts in suitable function spaces.This manifests that frictions have a stabilization effect on transonic shocks in ducts,in consideration of previous works which shown that transonic shocks in purely steady Euler flows are not stable in such ducts.Except its implications to applications,because frictions lead to a stronger coupling between the elliptic and hyperbolic parts of the three-dimensional steady subsonic Euler system,we develop the framework established in previous works to study more complex and interesting Venttsel problems of nonlocal elliptic equations.
基金Project supported by the National Science Foundation,the National Scaling Programthe Doctoral Foundation of the National Education Commission
文摘Petroleum science has made remarkable progress in organic geochcmistry and in the research into the theories of petroleum origin, its transport and accumulation. In estimating the oil-gas resources of a basin, the knowledge of its evolutionary history and especially the numerical computation of fluid flow and the history of its changes under heat is vital. The mathematical model can be described as a coupled system of nonlinear partial differentical equations with initial-boundary value problems. This thesis, from actual conditions such as the effect of fluid compressibility and the three-dimensional characteristic of large-scale science-engineering computation, we put forward a kind of characteristic finite element alternating-direction schemes and obtain optimal order estimates in L^2 norm for the error in the approximate assumption.
文摘In this paper. an analytic method is. presented. for the research of nonlinear Three-dimensional problems of composite laminated plates. The perturbation method and the variational principle are used to satisfy the basic differential equations and the boundary conditions of the three-dimensional theory of elasticity. The nonlinear three-dimensional problems are studied .for composite anisotropic circular laminas and laminates subjected to transverse loading. The perturbation series solutions of high accuracy are obtained. A large number of results show that transverse normal stress and transverse shear stresses are very important in the nonlinear three-dimensional analysis of laminated plates.
基金supported by the Natural Science Foundation of China(11801108)the Natural Science Foundation of Guangdong Province(2021A1515010314)the Science and Technology Planning Project of Guangzhou City(202201010111)。
文摘This paper deals with the forward and backward problems for the nonlinear fractional pseudo-parabolic equation ut+(-Δ)^(s1)ut+β(-Δ)^(s2)u=F(u,x,t)subject o random Gaussian white noise for initial and final data.Under the suitable assumptions s1,s2andβ,we first show the ill-posedness of mild solutions for forward and backward problems in the sense of Hadamard,which are mainly driven by random noise.Moreover,we propose the Fourier truncation method for stabilizing the above ill-posed problems.We derive an error estimate between the exact solution and its regularized solution in an E‖·‖Hs22norm,and give some numerical examples illustrating the effect of above method.
文摘The Magneto-acoustic Tomography with Current Injection (MAT-CI) is a new biological electrical impedance imaging technique that combines Electrical Impedance Tomography (EIT) with Ultrasonic Imaging (UI), which possesses the non-invasive and high-contrast of the EIT and the high-resolution of the UI. The MAT-CI is expected to acquire high quality image and embraces a wide application. Its principle is to put the conductive sample in the Static Magnetic Field(SMF) and inject a time-varying current, during which the SMF and the current interact and generate the Lorentz Force that inspire ultrasonic signal received by the ultrasonic transducers positioned around the sample. And then according to related reconstruction algorithm and ultrasonic signal, electrical conductivity image is obtained. In this paper, a forward problem mathematical model of the MAT-CI has been set up to deduce the theoretical equation of the electromagnetic field and solve the sound source distribution by Green’s function. Secondly, a sound field restoration by Wiener filtering and reconstruction of current density by time-rotating method have deduced the Laplace’s equation that caters to the current density to further acquire the electrical conductivity distribution image of the sample through iteration method. In the end, double-loop coils experiments have been conducted to verify its feasibility.
基金supported by Joint Fund of Advanced Aerospace Manufacturing Technology Research(No. U1937601)the partial financial support of this research by the budget program of Ukraine“Support for the Development of Priority Research Areas”(No.CPCEC 6451230)。
文摘By making use of the direct integration method,an exact analysis of the general three-dimensional thermoelasticity problem is performed for the case of a transversely isotropic homogeneous half-space subject to local thermal and force loadings.The material plane of isotropy is assumed to be parallel to the limiting surface of the halfspace.By reducing the original thermoelasticity equations to the governing ones for individual stress-tensor components,the effect of material anisotropy in the stress field is analyzed with regard to the feasibility requirement,i.e.,the finiteness of the stress field at a distance from the disturbed area.As a result,the solution is constructed in the form of explicit analytical dependencies on the force and thermal loadings for various kinds of transversely isotropic materials and agrees with the basic principles of the continua mechanics.The solution can be efficiently used as a benchmark one for the direct computation of temperature and thermal stresses in transversely isotropic semi-infinite domains,as well as for the verification of solutions constructed by different means.
基金Supported by the National Key Research and Development Program of China(No.2023YFC3008200)the Independent Research Project of Southern Marine Science and Engineering Guangdong Laboratory(Zhuhai)(No.SML2022SP505)。
文摘Physics-informed neural networks(PINNs),as a novel artificial intelligence method for solving partial differential equations,are applicable to solve both forward and inverse problems.This study evaluates the performance of PINNs in solving the temperature diffusion equation of the seawater across six scenarios,including forward and inverse problems under three different boundary conditions.Results demonstrate that PINNs achieved consistently higher accuracy with the Dirichlet and Neumann boundary conditions compared to the Robin boundary condition for both forward and inverse problems.Inaccurate weighting of terms in the loss function can reduce model accuracy.Additionally,the sensitivity of model performance to the positioning of sampling points varied between different boundary conditions.In particular,the model under the Dirichlet boundary condition exhibited superior robustness to variations in point positions during the solutions of inverse problems.In contrast,for the Neumann and Robin boundary conditions,accuracy declines when points were sampled from identical positions or at the same time.Subsequently,the Argo observations were used to reconstruct the vertical diffusion of seawater temperature in the north-central Pacific for the applicability of PINNs in the real ocean.The PINNs successfully captured the vertical diffusion characteristics of seawater temperature,reflected the seasonal changes of vertical temperature under different topographic conditions,and revealed the influence of topography on the temperature diffusion coefficient.The PINNs were proved effective in solving the temperature diffusion equation of seawater with limited data,providing a promising technique for simulating or predicting ocean phenomena using sparse observations.
