Energy resolution calibration is crucial for gamma-ray spectral analysis,as measured using a scintillation detector.A locally constrained regularization method was proposed to determine the resolution calibration para...Energy resolution calibration is crucial for gamma-ray spectral analysis,as measured using a scintillation detector.A locally constrained regularization method was proposed to determine the resolution calibration parameters.First,a Monte Carlo simulation model consistent with an actual measurement system was constructed to obtain the energy deposition distribution in the scintillation crystal.Subsequently,the regularization objective function is established based on weighted least squares and additional constraints.Additional constraints were designed using a special weighting scheme based on the incident gamma-ray energies.Subsequently,an intelligent algorithm was introduced to search for the optimal resolution calibration parameters by minimizing the objective function.The most appropriate regularization parameter was determined through mathematical experiments.When the regularization parameter was 30,the calibrated results exhibited the minimum RMSE.Simulations and test pit experiments were conducted to verify the performance of the proposed method.The simulation results demonstrate that the proposed algorithm can determine resolution calibration parameters more accurately than the traditional weighted least squares,and the test pit experimental results show that the R-squares between the calibrated and measured spectra are larger than 0.99.The accurate resolution calibration parameters determined by the proposed method lay the foundation for gamma-ray spectral processing and simulation benchmarking.展开更多
For cone beam computed tomography(CBCT),there has long been a desire to modulate the intensity and distribution of the X-rays to accommodate the patient’s anatomy as the gantry rotates from one projection to another....For cone beam computed tomography(CBCT),there has long been a desire to modulate the intensity and distribution of the X-rays to accommodate the patient’s anatomy as the gantry rotates from one projection to another.This would reduce both image artifacts and radiation dose.However,the current beam modulation setups,such as dynamic bowtie filters,may be too complex for practical use in clinical applications.This study aimed to investigate a simplified dynamic beam filtration strategy for CBCT imaging to reduce image artifacts and radiation dose.In this study,the beam filtration was designed to vary dynamically as the CBCT gantry rotates around the object.Specifically,two distinct components were integrated:the sheet filter part and the bowtie filter part.The dynamic beam filtration setup has two working schemes,one is a combination of dynamic sheet filter and dynamic bowtie filter,denoted as dynamic filterdynamic bowtie(DFDB);the other is a combination of dynamic sheet filter and static bowtie filter,denoted as dynamic filter-static bowtie(DFSB).Numerical imaging experiments were performed for three human body parts:the shoulder,chest,and knee.In addition,the Monte Carlo simulation platform MC-GPU was used to generate the dose distribution maps.Results showed that the proposed DFDB and DFSB beam filtration schemes can significantly reduce the image artifacts and thus improve the CBCT image quality.Depending on the scanned object,the total radiation dose could be reduced by 30%.The proposed simple dynamic beam filtration strategy,especially the DFSB approach,could be beneficial in the future to improve the CBCT image quality with reduced image artifacts and radiation dose.展开更多
This paper is a continuation of recent work by Guo-Xiang-Zheng[10].We deduce the sharp Morrey regularity theory for weak solutions to the fourth order nonhomogeneous Lamm-Rivière equation △^{2}u=△(V▽u)+div(w▽...This paper is a continuation of recent work by Guo-Xiang-Zheng[10].We deduce the sharp Morrey regularity theory for weak solutions to the fourth order nonhomogeneous Lamm-Rivière equation △^{2}u=△(V▽u)+div(w▽u)+(▽ω+F)·▽u+f in B^(4),under the smallest regularity assumptions of V,ω,ω,F,where f belongs to some Morrey spaces.This work was motivated by many geometrical problems such as the flow of biharmonic mappings.Our results deepens the Lp type regularity theory of[10],and generalizes the work of Du,Kang and Wang[4]on a second order problem to our fourth order problems.展开更多
Seismic data is commonly acquired sparsely and irregularly, which necessitates the regularization of seismic data with anti-aliasing and anti-leakage methods during seismic data processing. We propose a novel method o...Seismic data is commonly acquired sparsely and irregularly, which necessitates the regularization of seismic data with anti-aliasing and anti-leakage methods during seismic data processing. We propose a novel method of 4D anti-aliasing and anti-leakage Fourier transform using a cube-removal strategy to address the combination of irregular sampling and aliasing in high-dimensional seismic data. We compute a weighting function by stacking the spectrum along the radial lines, apply this function to suppress the aliasing energy, and then iteratively pick the dominant amplitude cube to construct the Fourier spectrum. The proposed method is very efficient due to a cube removal strategy for accelerating the convergence of Fourier reconstruction and a well-designed parallel architecture using CPU/GPU collaborative computing. To better fill the acquisition holes from 5D seismic data and meanwhile considering the GPU memory limitation, we developed the anti-aliasing and anti-leakage Fourier transform method in 4D with the remaining spatial dimension looped. The entire workflow is composed of three steps: data splitting, 4D regularization, and data merging. Numerical tests on both synthetic and field data examples demonstrate the high efficiency and effectiveness of our approach.展开更多
Field emission electric propulsion(FEEP) thrusters possess excellent characteristics, such as high specific impulse, low power requirements, compact size and precise pointing capabilities,making them ideal propulsion ...Field emission electric propulsion(FEEP) thrusters possess excellent characteristics, such as high specific impulse, low power requirements, compact size and precise pointing capabilities,making them ideal propulsion devices for micro-nano satellites. However, the detection of certain aspects, such as the evolution process of the liquid cone and the physical quantities at the cone apex, proves challenging due to the minute size of the needle tip and the vacuum environment in which they operate. Consequently, this paper introduces a computational fluid dynamics(CFD) model to gain insight into the formation process of the liquid cone on the tip apex of indium FEEP. The CFD model is based on electrohydrodynamic(EHD) equations and the volume of fluid(VOF) method. The entire cone formation process can be divided into three stages, and the time-dependent characteristics of the physical quantities at the cone apex are investigated. The influences of film thickness, apex radius size and applied voltage are compared.The results indicate a gradual increase in the values of electrostatic stress and surface tension stress at the cone apex over an initial period, followed by a rapid escalation within a short duration.Apex configurations featuring a small radius, thick film and high voltage exhibit a propensity for liquid cone formation, and the cone growth time decreases as the film thickness increases.Moreover, some unstable behavior is observed during the cone formation process.展开更多
In this paper,we study the three-dimensional regularized MHD equations with fractional Laplacians in the dissipative and diffusive terms.We establish the global existence of mild solutions to this system with small in...In this paper,we study the three-dimensional regularized MHD equations with fractional Laplacians in the dissipative and diffusive terms.We establish the global existence of mild solutions to this system with small initial data.In addition,we also obtain the Gevrey class regularity and the temporal decay rate of the solution.展开更多
Aiming at the time-optimal control problem of hypersonic vehicles(HSV)in ascending stage,a trigonometric regularization method(TRM)is introduced based on the indirect method of optimal control.This method avoids analy...Aiming at the time-optimal control problem of hypersonic vehicles(HSV)in ascending stage,a trigonometric regularization method(TRM)is introduced based on the indirect method of optimal control.This method avoids analyzing the switching function and distinguishing between singular control and bang-bang control,where the singular control problem is more complicated.While in bang-bang control,the costate variables are unsmooth due to the control jumping,resulting in difficulty in solving the two-point boundary value problem(TPBVP)induced by the indirect method.Aiming at the easy divergence when solving the TPBVP,the continuation method is introduced.This method uses the solution of the simplified problem as the initial value of the iteration.Then through solving a series of TPBVP,it approximates to the solution of the original complex problem.The calculation results show that through the above two methods,the time-optimal control problem of HSV in ascending stage under the complex model can be solved conveniently.展开更多
In this manuscript,we consider two kinds of the Fokker-Planck-type systems in the whole space.The first part involves proving the global existence and the algebraic time decay rates of the mild solutions to the Fokker...In this manuscript,we consider two kinds of the Fokker-Planck-type systems in the whole space.The first part involves proving the global existence and the algebraic time decay rates of the mild solutions to the Fokker-Planck-Boltzmann equation near Maxwellians if initial data satisfies some smallness in the function space L_(k)^(1)L_(T)^(∞)L_(v)^(2)∩L_(k)^(p)L_(T)^(∞)L_(v)^(2).The second part proves the global existence of the mild solutions to the Vlasov-Poisson-Fokker-Planck system in the function space L_(k)^(1)L_(T)^(∞)L_(v)^(2),and we also obtain the exponential time decay rates,which are different from the algebraic time decay rates of the Fokker-Planck-Boltzmann equation.Our analysis is based on Lk1LT∞Lv2function space introduced by Duan et al.(Comm Pure Appl Math,2021,74:932-1020),the L_(k)^(1)∩L_(k)^(p) approach developed by Duan et al.(SIAM J Math Anal,2024,56:762-800),and the coercivity property of the Fokker-Planck operator.However,it is worth pointing out that the L_(k)^(1)∩L_(k)^(p)approach is not required for the Vlasov-Poisson-Fokker-Planck system,due to the influence of the electric field term,which is different from the Fokker-Planck-Boltzmann equation in this paper and in the work of Duan et al.(SIAM J Math Anal,2024,56:762-800).展开更多
基金supported by the National Natural Science Foundation of China(No.41804141)。
文摘Energy resolution calibration is crucial for gamma-ray spectral analysis,as measured using a scintillation detector.A locally constrained regularization method was proposed to determine the resolution calibration parameters.First,a Monte Carlo simulation model consistent with an actual measurement system was constructed to obtain the energy deposition distribution in the scintillation crystal.Subsequently,the regularization objective function is established based on weighted least squares and additional constraints.Additional constraints were designed using a special weighting scheme based on the incident gamma-ray energies.Subsequently,an intelligent algorithm was introduced to search for the optimal resolution calibration parameters by minimizing the objective function.The most appropriate regularization parameter was determined through mathematical experiments.When the regularization parameter was 30,the calibrated results exhibited the minimum RMSE.Simulations and test pit experiments were conducted to verify the performance of the proposed method.The simulation results demonstrate that the proposed algorithm can determine resolution calibration parameters more accurately than the traditional weighted least squares,and the test pit experimental results show that the R-squares between the calibrated and measured spectra are larger than 0.99.The accurate resolution calibration parameters determined by the proposed method lay the foundation for gamma-ray spectral processing and simulation benchmarking.
文摘For cone beam computed tomography(CBCT),there has long been a desire to modulate the intensity and distribution of the X-rays to accommodate the patient’s anatomy as the gantry rotates from one projection to another.This would reduce both image artifacts and radiation dose.However,the current beam modulation setups,such as dynamic bowtie filters,may be too complex for practical use in clinical applications.This study aimed to investigate a simplified dynamic beam filtration strategy for CBCT imaging to reduce image artifacts and radiation dose.In this study,the beam filtration was designed to vary dynamically as the CBCT gantry rotates around the object.Specifically,two distinct components were integrated:the sheet filter part and the bowtie filter part.The dynamic beam filtration setup has two working schemes,one is a combination of dynamic sheet filter and dynamic bowtie filter,denoted as dynamic filterdynamic bowtie(DFDB);the other is a combination of dynamic sheet filter and static bowtie filter,denoted as dynamic filter-static bowtie(DFSB).Numerical imaging experiments were performed for three human body parts:the shoulder,chest,and knee.In addition,the Monte Carlo simulation platform MC-GPU was used to generate the dose distribution maps.Results showed that the proposed DFDB and DFSB beam filtration schemes can significantly reduce the image artifacts and thus improve the CBCT image quality.Depending on the scanned object,the total radiation dose could be reduced by 30%.The proposed simple dynamic beam filtration strategy,especially the DFSB approach,could be beneficial in the future to improve the CBCT image quality with reduced image artifacts and radiation dose.
基金supported by the National Natural Science Foundation of China(12271296,12271195).
文摘This paper is a continuation of recent work by Guo-Xiang-Zheng[10].We deduce the sharp Morrey regularity theory for weak solutions to the fourth order nonhomogeneous Lamm-Rivière equation △^{2}u=△(V▽u)+div(w▽u)+(▽ω+F)·▽u+f in B^(4),under the smallest regularity assumptions of V,ω,ω,F,where f belongs to some Morrey spaces.This work was motivated by many geometrical problems such as the flow of biharmonic mappings.Our results deepens the Lp type regularity theory of[10],and generalizes the work of Du,Kang and Wang[4]on a second order problem to our fourth order problems.
文摘Seismic data is commonly acquired sparsely and irregularly, which necessitates the regularization of seismic data with anti-aliasing and anti-leakage methods during seismic data processing. We propose a novel method of 4D anti-aliasing and anti-leakage Fourier transform using a cube-removal strategy to address the combination of irregular sampling and aliasing in high-dimensional seismic data. We compute a weighting function by stacking the spectrum along the radial lines, apply this function to suppress the aliasing energy, and then iteratively pick the dominant amplitude cube to construct the Fourier spectrum. The proposed method is very efficient due to a cube removal strategy for accelerating the convergence of Fourier reconstruction and a well-designed parallel architecture using CPU/GPU collaborative computing. To better fill the acquisition holes from 5D seismic data and meanwhile considering the GPU memory limitation, we developed the anti-aliasing and anti-leakage Fourier transform method in 4D with the remaining spatial dimension looped. The entire workflow is composed of three steps: data splitting, 4D regularization, and data merging. Numerical tests on both synthetic and field data examples demonstrate the high efficiency and effectiveness of our approach.
基金supported by National Natural Science Foundation of China(No.52075334)。
文摘Field emission electric propulsion(FEEP) thrusters possess excellent characteristics, such as high specific impulse, low power requirements, compact size and precise pointing capabilities,making them ideal propulsion devices for micro-nano satellites. However, the detection of certain aspects, such as the evolution process of the liquid cone and the physical quantities at the cone apex, proves challenging due to the minute size of the needle tip and the vacuum environment in which they operate. Consequently, this paper introduces a computational fluid dynamics(CFD) model to gain insight into the formation process of the liquid cone on the tip apex of indium FEEP. The CFD model is based on electrohydrodynamic(EHD) equations and the volume of fluid(VOF) method. The entire cone formation process can be divided into three stages, and the time-dependent characteristics of the physical quantities at the cone apex are investigated. The influences of film thickness, apex radius size and applied voltage are compared.The results indicate a gradual increase in the values of electrostatic stress and surface tension stress at the cone apex over an initial period, followed by a rapid escalation within a short duration.Apex configurations featuring a small radius, thick film and high voltage exhibit a propensity for liquid cone formation, and the cone growth time decreases as the film thickness increases.Moreover, some unstable behavior is observed during the cone formation process.
基金supported by the Opening Project of Guangdong Province Key Laboratory of Cyber-Physical System(20168030301008)supported by the National Natural Science Foundation of China(11126266)+4 种基金the Natural Science Foundation of Guangdong Province(2016A030313390)the Quality Engineering Project of Guangdong Province(SCAU-2021-69)the SCAU Fund for High-level University Buildingsupported by the National Key Research and Development Program of China(2020YFA0712500)the National Natural Science Foundation of China(11971496,12126609)。
文摘In this paper,we study the three-dimensional regularized MHD equations with fractional Laplacians in the dissipative and diffusive terms.We establish the global existence of mild solutions to this system with small initial data.In addition,we also obtain the Gevrey class regularity and the temporal decay rate of the solution.
基金supported by the Na-tional Natural Science Foundation of China(No.52272369).
文摘Aiming at the time-optimal control problem of hypersonic vehicles(HSV)in ascending stage,a trigonometric regularization method(TRM)is introduced based on the indirect method of optimal control.This method avoids analyzing the switching function and distinguishing between singular control and bang-bang control,where the singular control problem is more complicated.While in bang-bang control,the costate variables are unsmooth due to the control jumping,resulting in difficulty in solving the two-point boundary value problem(TPBVP)induced by the indirect method.Aiming at the easy divergence when solving the TPBVP,the continuation method is introduced.This method uses the solution of the simplified problem as the initial value of the iteration.Then through solving a series of TPBVP,it approximates to the solution of the original complex problem.The calculation results show that through the above two methods,the time-optimal control problem of HSV in ascending stage under the complex model can be solved conveniently.
基金supported by the National Natural Science Foundation of China(11801285,12326337)。
文摘In this manuscript,we consider two kinds of the Fokker-Planck-type systems in the whole space.The first part involves proving the global existence and the algebraic time decay rates of the mild solutions to the Fokker-Planck-Boltzmann equation near Maxwellians if initial data satisfies some smallness in the function space L_(k)^(1)L_(T)^(∞)L_(v)^(2)∩L_(k)^(p)L_(T)^(∞)L_(v)^(2).The second part proves the global existence of the mild solutions to the Vlasov-Poisson-Fokker-Planck system in the function space L_(k)^(1)L_(T)^(∞)L_(v)^(2),and we also obtain the exponential time decay rates,which are different from the algebraic time decay rates of the Fokker-Planck-Boltzmann equation.Our analysis is based on Lk1LT∞Lv2function space introduced by Duan et al.(Comm Pure Appl Math,2021,74:932-1020),the L_(k)^(1)∩L_(k)^(p) approach developed by Duan et al.(SIAM J Math Anal,2024,56:762-800),and the coercivity property of the Fokker-Planck operator.However,it is worth pointing out that the L_(k)^(1)∩L_(k)^(p)approach is not required for the Vlasov-Poisson-Fokker-Planck system,due to the influence of the electric field term,which is different from the Fokker-Planck-Boltzmann equation in this paper and in the work of Duan et al.(SIAM J Math Anal,2024,56:762-800).
