A Wentzel-Kramers-Brillouin(WKB)method is introduced for obtaining a uniform asymptotic solution for underwater sound propagation at very low frequencies in deep ocean.The method utilizes a mode sum and employs the re...A Wentzel-Kramers-Brillouin(WKB)method is introduced for obtaining a uniform asymptotic solution for underwater sound propagation at very low frequencies in deep ocean.The method utilizes a mode sum and employs the reference functions method to describe the solution to the depth-separated wave equation approximately using parabolic cylinder functions.The conditions for the validity of this approximation are also discussed.Furthermore,a formula that incorporates waveguide effects for the modal group velocity is derived,revealing that boundary effects at very low frequencies can have a significant impact on the propagation characteristics of even low-order normal modes.The present method not only offers improved accuracy compared to the classical WKB approximation and the uniform asymptotic approximation based on Airy functions,but also provides a wider range of depth applicability.Additionally,this method exhibits strong agreement with numerical methods and offers valuable physical insights.Finally,the method is applied to the study of very-low-frequency sound propagation in the South China Sea,leading to sound transmission loss predictions that closely align with experimental observations.展开更多
We report a linear-scaling random Green's function(rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projec...We report a linear-scaling random Green's function(rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projecting onto Krylov subspace. With the rGF method, the Fermi–Dirac operator can be obtained directly, avoiding the polynomial expansion to Fermi–Dirac function. To demonstrate the applicability, we implement the rGF method with the density-functional tight-binding method. It is shown that the Krylov subspace can maintain at small size for materials with different gaps at zero temperature, including H_(2)O and Si clusters. We find with a simple deflation technique that the rGF self-consistent calculation of H_(2)O clusters at T = 0 K can reach an error of~ 1 me V per H_(2)O molecule in total energy, compared to deterministic calculations. The rGF method provides an effective stochastic method for large-scale electronic structure simulation.展开更多
Recent advancements in additive manufacturing(AM)have revolutionized the design and production of complex engineering microstructures.Despite these advancements,their mathematical modeling and computational analysis r...Recent advancements in additive manufacturing(AM)have revolutionized the design and production of complex engineering microstructures.Despite these advancements,their mathematical modeling and computational analysis remain significant challenges.This research aims to develop an effective computational method for analyzing the free vibration of functionally graded(FG)microplates under high temperatures while resting on a Pasternak foundation(PF).This formulation leverages a new thirdorder shear deformation theory(new TSDT)for improved accuracy without requiring shear correction factors.Additionally,the modified couple stress theory(MCST)is incorporated to account for sizedependent effects in microplates.The PF is characterized by two parameters including spring stiffness(k_(w))and shear layer stiffness(k_(s)).To validate the proposed method,the results obtained are compared with those of the existing literature.Furthermore,numerical examples explore the influence of various factors on the high-temperature free vibration of FG microplates.These factors include the length scale parameter(l),geometric dimensions,material properties,and the presence of the elastic foundation.The findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the results of this research will have great potential in military and defense applications such as components of submarines,fighter aircraft,and missiles.展开更多
Some doubly-periodic solutions of the Zakharov-Kuznetsov equation are presented. Our approach is to introduce an auxiliary ordinary differential equation and use its Jacobi elliptic function solutions to construct dou...Some doubly-periodic solutions of the Zakharov-Kuznetsov equation are presented. Our approach is to introduce an auxiliary ordinary differential equation and use its Jacobi elliptic function solutions to construct doubly-periodic solutions of the Zakharov-Kuznetsov equation, which has been derived by Gottwald as a two-dimensional model for nonlinear Rossby waves. When the modulus k →1, these solutions reduce to the solitary wave solutions of the equation.展开更多
The distribution of the final surface subsidence basin induced by longwall operations in inclined coal seam could be significantly different from that in flat coal seam and demands special prediction methods. Though m...The distribution of the final surface subsidence basin induced by longwall operations in inclined coal seam could be significantly different from that in flat coal seam and demands special prediction methods. Though many empirical prediction methods have been developed, these methods are inflexible for varying geological and mining conditions. An influence function method has been developed to take the advantage of its fundamentally sound nature and flexibility. In developing this method, significant modifications have been made to the original Knothe function to produce an asymmetrical influence function. The empirical equations for final subsidence parameters derived from US subsidence data and Chinese empirical values have been incorpo- rated into the mathematical models to improve the prediction accuracy. A corresponding computer program is developed. A number of subsidence cases for longwall mining operations in coal seams with varying inclination angles have been used to demonstrate the applicability of the developed subsidence prediction model.展开更多
In this paper, the trial function method is extended to study the generalized nonlinear Schrodinger equation with time- dependent coefficients. On the basis of a generalized traveling wave transformation and a trial f...In this paper, the trial function method is extended to study the generalized nonlinear Schrodinger equation with time- dependent coefficients. On the basis of a generalized traveling wave transformation and a trial function, we investigate the exact envelope traveling wave solutions of the generalized nonlinear Schrodinger equation with time-dependent coefficients. Taking advantage of solutions to trial function, we successfully obtain exact solutions for the generalized nonlinear Schrodinger equation with time-dependent coefficients under constraint conditions.展开更多
In this paper, a new algorithm-approximate penalty function method is designed, which can be used to solve a bilevel optimization problem with linear constrained function. In this kind of bilevel optimization problem....In this paper, a new algorithm-approximate penalty function method is designed, which can be used to solve a bilevel optimization problem with linear constrained function. In this kind of bilevel optimization problem. the evaluation of the objective function is very difficult, so that only their approximate values can be obtained. This algorithm is obtained by combining penalty function method and approximation in bilevel programming. The presented algorithm is completely different from existing methods. That convergence for this algorithm is proved.展开更多
The possible exotic nuclear properties in the neutron-rich Ca,Ni,Zr,and Sn isotopes are examined with the continuum Skyrme Hartree-Fock-Bogoliubov theory in the framework of the Green’s function method.The pairing co...The possible exotic nuclear properties in the neutron-rich Ca,Ni,Zr,and Sn isotopes are examined with the continuum Skyrme Hartree-Fock-Bogoliubov theory in the framework of the Green’s function method.The pairing correlation,the couplings with the continuum,and the blocking effects for the unpaired nucleon in odd-A nuclei are properly treated.The Skyrme interaction SLy4 is adopted for the ph channel and the density-dependentinteraction is adopted for the pp chan-nel,which well reproduce the experimental two-neutron separation energies S_(2n)and one-neutron separation energies Sn.It is found that the criterion S_(n)>0 predicts a neutron drip line with neutron numbers much smaller than those for S_(2n)>0.Owing to the unpaired odd neutron,the neutron pairing energies−E_(pair)in odd-A nuclei are much lower than those in the neighbor-ing even-even nuclei.By investigating the single-particle structures,the possible halo structures in the neutron-rich Ca,Ni,and Sn isotopes are predicted,where sharp increases in the root-mean-square(rms)radii with significant deviations from the traditional rA^(1∕3)rule and diffuse spatial density distributions are observed.Analyzing the contributions of various partial waves to the total neutron densityρlj(r)∕ρ(r)reveals that the orbitals located around the Fermi surface-particularly those with small angular momenta-significantly affect the extended nuclear density and large rms radii.The number of neutrons Nλ(N_(0))occupying above the Fermi surfacen(continuum threshold)is discussed,whose evolution as a function of the mass number A in each isotope is consistent with that of the pairing energy,supporting the key role of the pairing correlation in halo phenomena.展开更多
Water function classification is to divide the waters in a basin or a region into different water function regions according to the natural features such as water resources condition, physical geographical location, e...Water function classification is to divide the waters in a basin or a region into different water function regions according to the natural features such as water resources condition, physical geographical location, environmental condition, etc. and social features such as the status quo of development and utilization, the requirement of social and economic development on water quantity and quality etc. Water function division has not ever been carried out in China and no ready-made theory and method can be complied with, in the meantime, it is a fundamental work with strong practical function. Therefore, the basic concept and theory foundation and research method are put forward on the base of summary of water resources management and developed and will be perfected in practice process of water function division. The Yangtze function zone covers the Yangtze river basin, the Lancang river basin and the rivers in the western area of the Lancang river in Southwest China. According to the technical outline of national water function division and combined with division practice of the Yangtze function zone, this paper presents the scope determination, procedure and method of water function division.展开更多
The finite volume time domain(FVTD) algorithm and Green function algorithm are extended to Schwarzschild spacetime for numerical simulation of electromagnetic scattering. The FVTD method in Schwarzschild spacetime is ...The finite volume time domain(FVTD) algorithm and Green function algorithm are extended to Schwarzschild spacetime for numerical simulation of electromagnetic scattering. The FVTD method in Schwarzschild spacetime is developed by filling the flat spacetime with an equivalent medium. The Green function in Schwarzschild spacetime is acquired by solving initial value problems. Both the FVTD code and the Green function code are validated by numerical results. Scattering in Schwarzschild spacetime is simulated with these methods.展开更多
The main purpose of this paper is to present numerical results of static bending and free vibration of functionally graded porous(FGP) variable-thickness plates by using an edge-based smoothed finite element method(ES...The main purpose of this paper is to present numerical results of static bending and free vibration of functionally graded porous(FGP) variable-thickness plates by using an edge-based smoothed finite element method(ES-FEM) associate with the mixed interpolation of tensorial components technique for the three-node triangular element(MITC3), so-called ES-MITC3. This ES-MITC3 element is performed to eliminate the shear locking problem and to enhance the accuracy of the existing MITC3 element. In the ES-MITC3 element, the stiffness matrices are obtained by using the strain smoothing technique over the smoothing domains formed by two adjacent MITC3 triangular elements sharing an edge. Materials of the plate are FGP with a power-law index(k) and maximum porosity distributions(U) in the forms of cosine functions. The influences of some geometric parameters, material properties on static bending, and natural frequency of the FGP variable-thickness plates are examined in detail.展开更多
The atmospheric scattering optical transfer function （OTF） is solved by applying the multi-coupled single scattering （MCSS） method to the three-dimensional radiative transfer equation （RTE） under the periodic gr...The atmospheric scattering optical transfer function （OTF） is solved by applying the multi-coupled single scattering （MCSS） method to the three-dimensional radiative transfer equation （RTE） under the periodic ground condition. This approach is a direct hit to the atmospheric scattering OTF using the same original context of modulation transfer function （MTF） measurement, i.e., images of sinusoidal grating at different spatial frequencies. Both the amplitude and phase shift of the OTF at various zenith and azimuth angles can be obtained at an arbitrary spatial frequency.展开更多
The multilevel characteristic basis function method(MLCBFM)with the adaptive cross approximation(ACA)algorithm for accelerated solution of electrically large scattering problems is studied in this paper.In the convent...The multilevel characteristic basis function method(MLCBFM)with the adaptive cross approximation(ACA)algorithm for accelerated solution of electrically large scattering problems is studied in this paper.In the conventional MLCBFM based on Foldy-Lax multiple scattering equations,the improvement is only made in the generation of characteristic basis functions(CBFs).However,it does not provide a change in impedance matrix filling and reducing matrix calculation procedure,which is time-consuming.In reality,all the impedance and reduced matrix of each level of the MLCBFM have low-rank property and can be calculated efficiently.Therefore,ACA is used for the efficient generation of two-level CBFs and the fast calculation of reduced matrix in this study.Numerical results are given to demonstrate the accuracy and efficiency of the method.展开更多
This paper presents a self-contained description on the configuration of propagator method(PM)to calculate the electron velocity distribution function(EVDF) of electron swarms in gases under DC electric and magnetic f...This paper presents a self-contained description on the configuration of propagator method(PM)to calculate the electron velocity distribution function(EVDF) of electron swarms in gases under DC electric and magnetic fields crossed at a right angle. Velocity space is divided into cells with respect to three polar coordinates v,θ and f. The number of electrons in each cell is stored in three-dimensional arrays. The changes of electron velocity due to acceleration by the electric and magnetic fields and scattering by gas molecules are treated as intercellular electron transfers on the basis of the Boltzmann equation and are represented using operators called the propagators or Green’s functions. The collision propagator, assuming isotropic scattering, is basically unchanged from conventional PMs performed under electric fields without magnetic fields. On the other hand, the acceleration propagator is customized for rotational acceleration under the action of the Lorentz force. The acceleration propagator specific to the present cell configuration is analytically derived. The mean electron energy and average electron velocity vector in a model gas and SF6 were derived from the EVDF as a demonstration of the PM under the Hall deflection and they were in a fine agreement with those obtained by Monte Carlo simulations. A strategy for fast relaxation is discussed, and extension of the PM for the EVDF under AC electric and DC/AC magnetic fields is outlined as well.展开更多
The Zakharov equation to describe the laser plasma interaction process has very important sense, this paper gives the solitary wave solutions for Zakharov equation by using Jacobi elliptic function method.
The original temporal clustering analysis (OTCA) is an effective technique for obtaining brain activation maps when the timing and location of the activation are completely unknown, but its deficiency of sensitivity i...The original temporal clustering analysis (OTCA) is an effective technique for obtaining brain activation maps when the timing and location of the activation are completely unknown, but its deficiency of sensitivity is exposed in processing brain activation signal which is relatively weak. The time slice analysis method based on OTCA is proposed considering the weakness of the functional magnetic resonance imaging (fMRI) signal of the rat model. By dividing the stimulation period into several time slices and analyzing each slice to detect the activated pixels respectively after the background removal, the sensitivity is significantly improved. The inhibitory response in the hypothalamus after glucose loading is detected successfully with this method in the experiment on rat. Combined with the OTCA method, the time slice analysis method based on OTCA is effective on detecting when, where and which type of response will happen after stimulation, even if the fMRI signal is weak.展开更多
In this paper,we propose an iterative algorithm to find the optimal incentive mechanism for the principal-agent problem under moral hazard where the number of agent action profiles is infinite,and where there are an i...In this paper,we propose an iterative algorithm to find the optimal incentive mechanism for the principal-agent problem under moral hazard where the number of agent action profiles is infinite,and where there are an infinite number of results that can be observed by the principal.This principal-agent problem has an infinite number of incentive-compatibility constraints,and we transform it into an optimization problem with an infinite number of constraints called a semi-infinite programming problem.We then propose an exterior penalty function method to find the optimal solution to this semi-infinite programming and illustrate the convergence of this algorithm.By analyzing the optimal solution obtained by the proposed penalty function method,we can obtain the optimal incentive mechanism for the principal-agent problem with an infinite number of incentive-compatibility constraints under moral hazard.展开更多
The realization of protein functional movement is usually accompanied by specific conformational changes,and there exist some key residues that mediate and control the functional motions of proteins in the allosteric ...The realization of protein functional movement is usually accompanied by specific conformational changes,and there exist some key residues that mediate and control the functional motions of proteins in the allosteric process.In the present work,the perturbation-response scanning method developed by our group was combined with the molecular dynamics(MD)simulation to identify the key residues controlling the functional movement of proteins.In our method,a physical quantity that is directly related to protein specific function was introduced,and then based on the MD simulation trajectories,the perturbation-response scanning method was used to identify the key residues for functional motions,in which the residues that highly correlated with the fluctuation of the function-related quantity were identified as the key residues controlling the specific functional motions of the protein.Two protein systems,i.e.,the heat shock protein 70 and glutamine binding protein,were selected as case studies to validate the effectiveness of our method.Our calculated results are in good agreement with the experimental results.The location of the key residues in the two proteins are similar,indicating the similar mechanisms behind the performance of their biological functions.展开更多
The aim of this study is to evaluate the uncertainty of 2πα and 2πβ surface emission rates using the windowless multiwire proportional counter method.This study used the Monte Carlo method (MCM) to validate the co...The aim of this study is to evaluate the uncertainty of 2πα and 2πβ surface emission rates using the windowless multiwire proportional counter method.This study used the Monte Carlo method (MCM) to validate the conventional Guide to the Expression of Uncertainty in Measurement (GUM) method.A dead time measurement model for the two-source method was established based on the characteristics of a single-channel measurement system,and the voltage threshold correction factor measurement function was indirectly obtained by fitting the threshold correction curve.The uncertainty in the surface emission rate was calculated using the GUM method and the law of propagation of uncertainty.The MCM provided clear definitions for each input quantity and its uncertainty distribution,and the simulation training was realized with a complete and complex mathematical model.The results of the surface emission rate uncertainty evaluation for four radioactive plane sources using both methods showed the uncertainty’s consistency E_(n)<0.070 for the comparison of each source,and the uncertainty results of the GUM were all lower than those of the MCM.However,the MCM has a more objective evaluation process and can serve as a validation tool for GUM results.展开更多
For an arbitrary solution to the Volterra lattice hierarchy,the logarithmic derivatives of the tau-function of the solution can be computed by the matrix-resolvent method.In this paper,we define a pair of wave functio...For an arbitrary solution to the Volterra lattice hierarchy,the logarithmic derivatives of the tau-function of the solution can be computed by the matrix-resolvent method.In this paper,we define a pair of wave functions of the solution and use them to give an expression of the matrix resolvent;based on this we obtain a new formula for the k-point functions for the Volterra lattice hierarchy in terms of wave functions.As an application,we give an explicit formula of k-point functions for the even GUE(Gaussian Unitary Ensemble)correlators.展开更多
基金Project supported by the National Natural Science Foundation of China(Grant Nos.12174048 and 12204128)。
文摘A Wentzel-Kramers-Brillouin(WKB)method is introduced for obtaining a uniform asymptotic solution for underwater sound propagation at very low frequencies in deep ocean.The method utilizes a mode sum and employs the reference functions method to describe the solution to the depth-separated wave equation approximately using parabolic cylinder functions.The conditions for the validity of this approximation are also discussed.Furthermore,a formula that incorporates waveguide effects for the modal group velocity is derived,revealing that boundary effects at very low frequencies can have a significant impact on the propagation characteristics of even low-order normal modes.The present method not only offers improved accuracy compared to the classical WKB approximation and the uniform asymptotic approximation based on Airy functions,but also provides a wider range of depth applicability.Additionally,this method exhibits strong agreement with numerical methods and offers valuable physical insights.Finally,the method is applied to the study of very-low-frequency sound propagation in the South China Sea,leading to sound transmission loss predictions that closely align with experimental observations.
基金financial support from the National Natural Science Foundation of China (Grant No. 12227901)the financial support from the National Natural Science Foundation of China (Grant Nos. 11974263 and 12174291)。
文摘We report a linear-scaling random Green's function(rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projecting onto Krylov subspace. With the rGF method, the Fermi–Dirac operator can be obtained directly, avoiding the polynomial expansion to Fermi–Dirac function. To demonstrate the applicability, we implement the rGF method with the density-functional tight-binding method. It is shown that the Krylov subspace can maintain at small size for materials with different gaps at zero temperature, including H_(2)O and Si clusters. We find with a simple deflation technique that the rGF self-consistent calculation of H_(2)O clusters at T = 0 K can reach an error of~ 1 me V per H_(2)O molecule in total energy, compared to deterministic calculations. The rGF method provides an effective stochastic method for large-scale electronic structure simulation.
文摘Recent advancements in additive manufacturing(AM)have revolutionized the design and production of complex engineering microstructures.Despite these advancements,their mathematical modeling and computational analysis remain significant challenges.This research aims to develop an effective computational method for analyzing the free vibration of functionally graded(FG)microplates under high temperatures while resting on a Pasternak foundation(PF).This formulation leverages a new thirdorder shear deformation theory(new TSDT)for improved accuracy without requiring shear correction factors.Additionally,the modified couple stress theory(MCST)is incorporated to account for sizedependent effects in microplates.The PF is characterized by two parameters including spring stiffness(k_(w))and shear layer stiffness(k_(s)).To validate the proposed method,the results obtained are compared with those of the existing literature.Furthermore,numerical examples explore the influence of various factors on the high-temperature free vibration of FG microplates.These factors include the length scale parameter(l),geometric dimensions,material properties,and the presence of the elastic foundation.The findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the results of this research will have great potential in military and defense applications such as components of submarines,fighter aircraft,and missiles.
文摘Some doubly-periodic solutions of the Zakharov-Kuznetsov equation are presented. Our approach is to introduce an auxiliary ordinary differential equation and use its Jacobi elliptic function solutions to construct doubly-periodic solutions of the Zakharov-Kuznetsov equation, which has been derived by Gottwald as a two-dimensional model for nonlinear Rossby waves. When the modulus k →1, these solutions reduce to the solitary wave solutions of the equation.
文摘The distribution of the final surface subsidence basin induced by longwall operations in inclined coal seam could be significantly different from that in flat coal seam and demands special prediction methods. Though many empirical prediction methods have been developed, these methods are inflexible for varying geological and mining conditions. An influence function method has been developed to take the advantage of its fundamentally sound nature and flexibility. In developing this method, significant modifications have been made to the original Knothe function to produce an asymmetrical influence function. The empirical equations for final subsidence parameters derived from US subsidence data and Chinese empirical values have been incorpo- rated into the mathematical models to improve the prediction accuracy. A corresponding computer program is developed. A number of subsidence cases for longwall mining operations in coal seams with varying inclination angles have been used to demonstrate the applicability of the developed subsidence prediction model.
基金Project supported in part by the National Natural Science Foundation of China(Grant No.11071177)
文摘In this paper, the trial function method is extended to study the generalized nonlinear Schrodinger equation with time- dependent coefficients. On the basis of a generalized traveling wave transformation and a trial function, we investigate the exact envelope traveling wave solutions of the generalized nonlinear Schrodinger equation with time-dependent coefficients. Taking advantage of solutions to trial function, we successfully obtain exact solutions for the generalized nonlinear Schrodinger equation with time-dependent coefficients under constraint conditions.
文摘In this paper, a new algorithm-approximate penalty function method is designed, which can be used to solve a bilevel optimization problem with linear constrained function. In this kind of bilevel optimization problem. the evaluation of the objective function is very difficult, so that only their approximate values can be obtained. This algorithm is obtained by combining penalty function method and approximation in bilevel programming. The presented algorithm is completely different from existing methods. That convergence for this algorithm is proved.
基金the National Natural Science Foundation of China(No.U2032141)the Open Project of Guangxi Key Laboratory of Nuclear Physics and Nuclear Technology(No.NLK2022-02)+4 种基金the Central Government Guidance Funds for Local Scientific and Technological Development,China(Guike ZY22096024)the Natural Science Foundation of Henan Province(No.202300410479)the Guizhou Provincial Science and Technology Projects(No.ZK[2022]203)the Foundation of Fundamental Research for Young Teachers of Zhengzhou University(No.JC202041041)the Physics Research and Development Program of Zhengzhou University(No.32410217).
文摘The possible exotic nuclear properties in the neutron-rich Ca,Ni,Zr,and Sn isotopes are examined with the continuum Skyrme Hartree-Fock-Bogoliubov theory in the framework of the Green’s function method.The pairing correlation,the couplings with the continuum,and the blocking effects for the unpaired nucleon in odd-A nuclei are properly treated.The Skyrme interaction SLy4 is adopted for the ph channel and the density-dependentinteraction is adopted for the pp chan-nel,which well reproduce the experimental two-neutron separation energies S_(2n)and one-neutron separation energies Sn.It is found that the criterion S_(n)>0 predicts a neutron drip line with neutron numbers much smaller than those for S_(2n)>0.Owing to the unpaired odd neutron,the neutron pairing energies−E_(pair)in odd-A nuclei are much lower than those in the neighbor-ing even-even nuclei.By investigating the single-particle structures,the possible halo structures in the neutron-rich Ca,Ni,and Sn isotopes are predicted,where sharp increases in the root-mean-square(rms)radii with significant deviations from the traditional rA^(1∕3)rule and diffuse spatial density distributions are observed.Analyzing the contributions of various partial waves to the total neutron densityρlj(r)∕ρ(r)reveals that the orbitals located around the Fermi surface-particularly those with small angular momenta-significantly affect the extended nuclear density and large rms radii.The number of neutrons Nλ(N_(0))occupying above the Fermi surfacen(continuum threshold)is discussed,whose evolution as a function of the mass number A in each isotope is consistent with that of the pairing energy,supporting the key role of the pairing correlation in halo phenomena.
文摘Water function classification is to divide the waters in a basin or a region into different water function regions according to the natural features such as water resources condition, physical geographical location, environmental condition, etc. and social features such as the status quo of development and utilization, the requirement of social and economic development on water quantity and quality etc. Water function division has not ever been carried out in China and no ready-made theory and method can be complied with, in the meantime, it is a fundamental work with strong practical function. Therefore, the basic concept and theory foundation and research method are put forward on the base of summary of water resources management and developed and will be perfected in practice process of water function division. The Yangtze function zone covers the Yangtze river basin, the Lancang river basin and the rivers in the western area of the Lancang river in Southwest China. According to the technical outline of national water function division and combined with division practice of the Yangtze function zone, this paper presents the scope determination, procedure and method of water function division.
基金Supported by the National Natural Science Foundation of China under Grant No 61601105
文摘The finite volume time domain(FVTD) algorithm and Green function algorithm are extended to Schwarzschild spacetime for numerical simulation of electromagnetic scattering. The FVTD method in Schwarzschild spacetime is developed by filling the flat spacetime with an equivalent medium. The Green function in Schwarzschild spacetime is acquired by solving initial value problems. Both the FVTD code and the Green function code are validated by numerical results. Scattering in Schwarzschild spacetime is simulated with these methods.
基金funded by Vietnam National Foundation for Science and Technology Development (NAFOSTED) under Grant number 107.02-2019.330。
文摘The main purpose of this paper is to present numerical results of static bending and free vibration of functionally graded porous(FGP) variable-thickness plates by using an edge-based smoothed finite element method(ES-FEM) associate with the mixed interpolation of tensorial components technique for the three-node triangular element(MITC3), so-called ES-MITC3. This ES-MITC3 element is performed to eliminate the shear locking problem and to enhance the accuracy of the existing MITC3 element. In the ES-MITC3 element, the stiffness matrices are obtained by using the strain smoothing technique over the smoothing domains formed by two adjacent MITC3 triangular elements sharing an edge. Materials of the plate are FGP with a power-law index(k) and maximum porosity distributions(U) in the forms of cosine functions. The influences of some geometric parameters, material properties on static bending, and natural frequency of the FGP variable-thickness plates are examined in detail.
基金supported by the Main Direction Program of Knowledge Innovation of Chinese Academy of Sciences(Grant No.KGFZD-125-13-006)
文摘The atmospheric scattering optical transfer function （OTF） is solved by applying the multi-coupled single scattering （MCSS） method to the three-dimensional radiative transfer equation （RTE） under the periodic ground condition. This approach is a direct hit to the atmospheric scattering OTF using the same original context of modulation transfer function （MTF） measurement, i.e., images of sinusoidal grating at different spatial frequencies. Both the amplitude and phase shift of the OTF at various zenith and azimuth angles can be obtained at an arbitrary spatial frequency.
基金supported by the National Natural Science Foundation of China (No.61401003)the Specialized Research Fund for the Doctoral Program of Higher Education of China (No.20123401110006)the Natural Science Research Project of Anhui Education ( No. KJ2015A436)
文摘The multilevel characteristic basis function method(MLCBFM)with the adaptive cross approximation(ACA)algorithm for accelerated solution of electrically large scattering problems is studied in this paper.In the conventional MLCBFM based on Foldy-Lax multiple scattering equations,the improvement is only made in the generation of characteristic basis functions(CBFs).However,it does not provide a change in impedance matrix filling and reducing matrix calculation procedure,which is time-consuming.In reality,all the impedance and reduced matrix of each level of the MLCBFM have low-rank property and can be calculated efficiently.Therefore,ACA is used for the efficient generation of two-level CBFs and the fast calculation of reduced matrix in this study.Numerical results are given to demonstrate the accuracy and efficiency of the method.
文摘This paper presents a self-contained description on the configuration of propagator method(PM)to calculate the electron velocity distribution function(EVDF) of electron swarms in gases under DC electric and magnetic fields crossed at a right angle. Velocity space is divided into cells with respect to three polar coordinates v,θ and f. The number of electrons in each cell is stored in three-dimensional arrays. The changes of electron velocity due to acceleration by the electric and magnetic fields and scattering by gas molecules are treated as intercellular electron transfers on the basis of the Boltzmann equation and are represented using operators called the propagators or Green’s functions. The collision propagator, assuming isotropic scattering, is basically unchanged from conventional PMs performed under electric fields without magnetic fields. On the other hand, the acceleration propagator is customized for rotational acceleration under the action of the Lorentz force. The acceleration propagator specific to the present cell configuration is analytically derived. The mean electron energy and average electron velocity vector in a model gas and SF6 were derived from the EVDF as a demonstration of the PM under the Hall deflection and they were in a fine agreement with those obtained by Monte Carlo simulations. A strategy for fast relaxation is discussed, and extension of the PM for the EVDF under AC electric and DC/AC magnetic fields is outlined as well.
文摘The Zakharov equation to describe the laser plasma interaction process has very important sense, this paper gives the solitary wave solutions for Zakharov equation by using Jacobi elliptic function method.
基金the National Natural Science Foundation of China (30370432)
文摘The original temporal clustering analysis (OTCA) is an effective technique for obtaining brain activation maps when the timing and location of the activation are completely unknown, but its deficiency of sensitivity is exposed in processing brain activation signal which is relatively weak. The time slice analysis method based on OTCA is proposed considering the weakness of the functional magnetic resonance imaging (fMRI) signal of the rat model. By dividing the stimulation period into several time slices and analyzing each slice to detect the activated pixels respectively after the background removal, the sensitivity is significantly improved. The inhibitory response in the hypothalamus after glucose loading is detected successfully with this method in the experiment on rat. Combined with the OTCA method, the time slice analysis method based on OTCA is effective on detecting when, where and which type of response will happen after stimulation, even if the fMRI signal is weak.
基金supported by National Natural Science Foundation of China(72031009 and 71871171)the National Social Science Foundation of China(20&ZD058).
文摘In this paper,we propose an iterative algorithm to find the optimal incentive mechanism for the principal-agent problem under moral hazard where the number of agent action profiles is infinite,and where there are an infinite number of results that can be observed by the principal.This principal-agent problem has an infinite number of incentive-compatibility constraints,and we transform it into an optimization problem with an infinite number of constraints called a semi-infinite programming problem.We then propose an exterior penalty function method to find the optimal solution to this semi-infinite programming and illustrate the convergence of this algorithm.By analyzing the optimal solution obtained by the proposed penalty function method,we can obtain the optimal incentive mechanism for the principal-agent problem with an infinite number of incentive-compatibility constraints under moral hazard.
文摘The realization of protein functional movement is usually accompanied by specific conformational changes,and there exist some key residues that mediate and control the functional motions of proteins in the allosteric process.In the present work,the perturbation-response scanning method developed by our group was combined with the molecular dynamics(MD)simulation to identify the key residues controlling the functional movement of proteins.In our method,a physical quantity that is directly related to protein specific function was introduced,and then based on the MD simulation trajectories,the perturbation-response scanning method was used to identify the key residues for functional motions,in which the residues that highly correlated with the fluctuation of the function-related quantity were identified as the key residues controlling the specific functional motions of the protein.Two protein systems,i.e.,the heat shock protein 70 and glutamine binding protein,were selected as case studies to validate the effectiveness of our method.Our calculated results are in good agreement with the experimental results.The location of the key residues in the two proteins are similar,indicating the similar mechanisms behind the performance of their biological functions.
文摘The aim of this study is to evaluate the uncertainty of 2πα and 2πβ surface emission rates using the windowless multiwire proportional counter method.This study used the Monte Carlo method (MCM) to validate the conventional Guide to the Expression of Uncertainty in Measurement (GUM) method.A dead time measurement model for the two-source method was established based on the characteristics of a single-channel measurement system,and the voltage threshold correction factor measurement function was indirectly obtained by fitting the threshold correction curve.The uncertainty in the surface emission rate was calculated using the GUM method and the law of propagation of uncertainty.The MCM provided clear definitions for each input quantity and its uncertainty distribution,and the simulation training was realized with a complete and complex mathematical model.The results of the surface emission rate uncertainty evaluation for four radioactive plane sources using both methods showed the uncertainty’s consistency E_(n)<0.070 for the comparison of each source,and the uncertainty results of the GUM were all lower than those of the MCM.However,the MCM has a more objective evaluation process and can serve as a validation tool for GUM results.
基金supported by the National Key R and D Program of China(2020YFA0713100).
文摘For an arbitrary solution to the Volterra lattice hierarchy,the logarithmic derivatives of the tau-function of the solution can be computed by the matrix-resolvent method.In this paper,we define a pair of wave functions of the solution and use them to give an expression of the matrix resolvent;based on this we obtain a new formula for the k-point functions for the Volterra lattice hierarchy in terms of wave functions.As an application,we give an explicit formula of k-point functions for the even GUE(Gaussian Unitary Ensemble)correlators.
