The equivalent source(ES)method in the spherical coordinate system has been widely applied to processing,reduction,field modeling,and geophysical and geological interpretation of satellite magnetic anomaly data.Howeve...The equivalent source(ES)method in the spherical coordinate system has been widely applied to processing,reduction,field modeling,and geophysical and geological interpretation of satellite magnetic anomaly data.However,the inversion for the ES model suffers from nonuniqueness and instability,which remain unresolved.To mitigate these issues,we introduce both the minimum and flattest models into the model objective function as an alternative regularization approach in the spherical ES method.We first present the methods,then analyze the accuracy of forward calculation and test the proposed ES method in this study by using synthetic data.The experimental results from simulation data indicate that our proposed regularization effectively suppresses the Backus effect and mitigates inversion instability in the low-latitude region.Finally,we apply the proposed method to magnetic anomaly data from China Seismo-Electromagnetic Satellite-1(CSES-1)and Macao Science Satellite-1(MSS-1)magnetic measurements over Africa by constructing an ES model of the large-scale lithospheric magnetic field.Compared with existing global lithospheric magnetic field models,our ES model demonstrates good consistency at high altitudes and predicts more stable fields at low altitudes.Furthermore,we derive the reduction to the pole(RTP)magnetic anomaly fields and the apparent susceptibility contrast distribution based on the ES model.The latter correlates well with the regional tectonic framework in Africa and surroundings.展开更多
In order to overcome the non-unique solution problem of the conventional equivalent source method at the corresponding eigenfrequency,the basic guidelines of adding virtual damping into the equivalent source method is...In order to overcome the non-unique solution problem of the conventional equivalent source method at the corresponding eigenfrequency,the basic guidelines of adding virtual damping into the equivalent source method is proposed according to the damping characteristics of the imaginary part of the coordinates of the equivalent source point.The complex source point constructed by complex retraction is used to add virtual damping to the equivalent source system,and a complex retractive equivalent source method with unique solution in full wavenumber domain is proposed.The mechanism of the complex retraction equivalent source method in overcoming the non-unique solution problem is described in accordance with the zero-point characteristic of the characteristic indicator function,involving the theory,modeling method and selection of damping parameter.The simulation and analysis of typical acoustic radiation and acoustic scattering problems indicate that the present method can obtain the unique solution in the entire frequency domain simply and effectively for arbitrarily shaped acoustic source.Its computational accuracy and efficiency are equivalent to the conventional equivalent source method,which solves the problem of poor adaptability of the complex radius vector method.展开更多
The Local Monte Carlo(LMC)method is used to solve the problems of deep penetration and long time in the neutronics calculation of the radial neutron camera(RNC)diagnostic system on the experimental advanced supercondu...The Local Monte Carlo(LMC)method is used to solve the problems of deep penetration and long time in the neutronics calculation of the radial neutron camera(RNC)diagnostic system on the experimental advanced superconducting tokamak(EAST),and the radiation distribution of the RNC and the neutron flux at the detector positions of each channel are obtained.Compared with the results calculated by the global variance reduction method,it is shown that the LMC calculation is reliable within the reasonable error range.The calculation process of LMC is analyzed in detail,and the transport process of radiation particles is simulated in two steps.In the first step,an integrated neutronics model considering the complex window environment and a neutron source model based on EAST plasma shape are used to support the calculation.The particle information on the equivalent surface is analyzed to evaluate the rationality of settings of equivalent surface source and boundary.Based on the characteristic that only a local geometric model is needed in the second step,it is shown that the LMC is an advantageous calculation method for the nuclear shielding design of tokamak diagnostic systems.展开更多
The time-domain inverse technique based on the time-domain rotating equivalent source method has been proposed to localize and quantify rotating sound sources. However, this technique encounters two problems to be add...The time-domain inverse technique based on the time-domain rotating equivalent source method has been proposed to localize and quantify rotating sound sources. However, this technique encounters two problems to be addressed: one is the time-consuming process of solving the transcendental equation at each time step, and the other is the difculty of controlling the instability problem due to the time-varying transfer matrix. In view of that, an improved technique is proposed in this paper to resolve these two problems. In the improved technique, a de-Dopplerization method in the time-domain rotating reference frame is frst applied to eliminate the Doppler efect caused by the source rotation in the measured pressure signals, and then the restored pressure signals without the Doppler efect are used as the inputs of the time-domain stationary equivalent source method to locate and quantify sound sources. Compared with the original technique, the improved technique can avoid solving the transcendental equation at each time step, and facilitate the treatment of the instability problem because the transfer matrix does not change with time. Numerical simulation and experimental results show that the improved technique can eliminate the Doppler efect efectively, and then localize and quantify the rotating nonstationary or broadband sources accurately. The results also demonstrate that the improved technique can guarantee a more stable reconstruction and compute more efciently than the original one.展开更多
针对声学边界元法中解的非唯一性和奇异积分问题,基于组合亥姆霍兹积分方程公式(combined helmholtz integral equation formulation,CHIEF)法思想,将常规边界元方程和等效源方程进行联立,并利用两者方程系数矩阵间的耦合等价关系,间接...针对声学边界元法中解的非唯一性和奇异积分问题,基于组合亥姆霍兹积分方程公式(combined helmholtz integral equation formulation,CHIEF)法思想,将常规边界元方程和等效源方程进行联立,并利用两者方程系数矩阵间的耦合等价关系,间接替换计算常规边界元法中的奇异系数矩阵,进而提出一种具有全频域唯一解、高计算精度和高稳定性的耦合CHIEF法。该方法将等效源方程作为补充方程,不仅解决了传统CHIEF法内点补充方程失效的问题,而且矩阵的间接替换计算避免了直接计算奇异积分,显著提高了计算效率和精度。通过声辐射和声散射的典型算例对比了所提方法、常规边界元法、常规Burton-Miller法和等效源法的计算效果。结果表明,所提方法不仅在全波数域内均能获得唯一解,且其计算精度和效率均优于常规边界元法和常规Burton-Miller方法,其系数矩阵条件数远低于等效源法。展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.42250103 and 42174090)the Opening Fund of Key Laboratory of Geological Survey and Evaluation of Ministry of Education(Grant No.GLAB2023ZR02)the MOST Special Fund from the State Key Laboratory of Geological Processes and Mineral Resources(Grant No.MSFGPMR2022-4).
文摘The equivalent source(ES)method in the spherical coordinate system has been widely applied to processing,reduction,field modeling,and geophysical and geological interpretation of satellite magnetic anomaly data.However,the inversion for the ES model suffers from nonuniqueness and instability,which remain unresolved.To mitigate these issues,we introduce both the minimum and flattest models into the model objective function as an alternative regularization approach in the spherical ES method.We first present the methods,then analyze the accuracy of forward calculation and test the proposed ES method in this study by using synthetic data.The experimental results from simulation data indicate that our proposed regularization effectively suppresses the Backus effect and mitigates inversion instability in the low-latitude region.Finally,we apply the proposed method to magnetic anomaly data from China Seismo-Electromagnetic Satellite-1(CSES-1)and Macao Science Satellite-1(MSS-1)magnetic measurements over Africa by constructing an ES model of the large-scale lithospheric magnetic field.Compared with existing global lithospheric magnetic field models,our ES model demonstrates good consistency at high altitudes and predicts more stable fields at low altitudes.Furthermore,we derive the reduction to the pole(RTP)magnetic anomaly fields and the apparent susceptibility contrast distribution based on the ES model.The latter correlates well with the regional tectonic framework in Africa and surroundings.
基金supported by the National Natural Science Foundation of China(12364055,51775121).
文摘In order to overcome the non-unique solution problem of the conventional equivalent source method at the corresponding eigenfrequency,the basic guidelines of adding virtual damping into the equivalent source method is proposed according to the damping characteristics of the imaginary part of the coordinates of the equivalent source point.The complex source point constructed by complex retraction is used to add virtual damping to the equivalent source system,and a complex retractive equivalent source method with unique solution in full wavenumber domain is proposed.The mechanism of the complex retraction equivalent source method in overcoming the non-unique solution problem is described in accordance with the zero-point characteristic of the characteristic indicator function,involving the theory,modeling method and selection of damping parameter.The simulation and analysis of typical acoustic radiation and acoustic scattering problems indicate that the present method can obtain the unique solution in the entire frequency domain simply and effectively for arbitrarily shaped acoustic source.Its computational accuracy and efficiency are equivalent to the conventional equivalent source method,which solves the problem of poor adaptability of the complex radius vector method.
基金support and help in this research.This work was supported by Users with Excellence Program of Hefei Science Center CAS(No.2020HSC-UE012)Comprehensive Research Facility for Fusion Technology Program of China(No.2018-000052-73-01-001228)National Natural Science Foundation of China(No.11605241)。
文摘The Local Monte Carlo(LMC)method is used to solve the problems of deep penetration and long time in the neutronics calculation of the radial neutron camera(RNC)diagnostic system on the experimental advanced superconducting tokamak(EAST),and the radiation distribution of the RNC and the neutron flux at the detector positions of each channel are obtained.Compared with the results calculated by the global variance reduction method,it is shown that the LMC calculation is reliable within the reasonable error range.The calculation process of LMC is analyzed in detail,and the transport process of radiation particles is simulated in two steps.In the first step,an integrated neutronics model considering the complex window environment and a neutron source model based on EAST plasma shape are used to support the calculation.The particle information on the equivalent surface is analyzed to evaluate the rationality of settings of equivalent surface source and boundary.Based on the characteristic that only a local geometric model is needed in the second step,it is shown that the LMC is an advantageous calculation method for the nuclear shielding design of tokamak diagnostic systems.
基金Supported by National Natural Science Foundation of China(Grant Nos.51875147,12174082,51675149)。
文摘The time-domain inverse technique based on the time-domain rotating equivalent source method has been proposed to localize and quantify rotating sound sources. However, this technique encounters two problems to be addressed: one is the time-consuming process of solving the transcendental equation at each time step, and the other is the difculty of controlling the instability problem due to the time-varying transfer matrix. In view of that, an improved technique is proposed in this paper to resolve these two problems. In the improved technique, a de-Dopplerization method in the time-domain rotating reference frame is frst applied to eliminate the Doppler efect caused by the source rotation in the measured pressure signals, and then the restored pressure signals without the Doppler efect are used as the inputs of the time-domain stationary equivalent source method to locate and quantify sound sources. Compared with the original technique, the improved technique can avoid solving the transcendental equation at each time step, and facilitate the treatment of the instability problem because the transfer matrix does not change with time. Numerical simulation and experimental results show that the improved technique can eliminate the Doppler efect efectively, and then localize and quantify the rotating nonstationary or broadband sources accurately. The results also demonstrate that the improved technique can guarantee a more stable reconstruction and compute more efciently than the original one.
文摘针对声学边界元法中解的非唯一性和奇异积分问题,基于组合亥姆霍兹积分方程公式(combined helmholtz integral equation formulation,CHIEF)法思想,将常规边界元方程和等效源方程进行联立,并利用两者方程系数矩阵间的耦合等价关系,间接替换计算常规边界元法中的奇异系数矩阵,进而提出一种具有全频域唯一解、高计算精度和高稳定性的耦合CHIEF法。该方法将等效源方程作为补充方程,不仅解决了传统CHIEF法内点补充方程失效的问题,而且矩阵的间接替换计算避免了直接计算奇异积分,显著提高了计算效率和精度。通过声辐射和声散射的典型算例对比了所提方法、常规边界元法、常规Burton-Miller法和等效源法的计算效果。结果表明,所提方法不仅在全波数域内均能获得唯一解,且其计算精度和效率均优于常规边界元法和常规Burton-Miller方法,其系数矩阵条件数远低于等效源法。
