It is of crucial importance to investigate the spatial structures of ancient landslides in the eastern Tibetan Plateau’s alpine canyons as they could provide valuable insights into the evolutionary history of the lan...It is of crucial importance to investigate the spatial structures of ancient landslides in the eastern Tibetan Plateau’s alpine canyons as they could provide valuable insights into the evolutionary history of the landslides and indicate the potential for future reactivation.This study examines the Deda ancient landslide,situated in the Chalong-ranbu fault zone,where creep deformation suggests a complex underground structure.By integrating remote sensing,field surveys,Audio-frequency Magnetotellurics(AMT),and Microtremor Survey Method(MSM)techniques,along with engineering geological drilling for validation,to uncover the landslide’s spatial feature s.The research indicates that a fault is developed in the upper part of the Deda ancient landslide,and the gully divides it into Deda landslide accumulation zoneⅠand Deda landslide accumulation zoneⅡin space.The distinctive geological characteristics detectable by MSM in the shallow subsurface and by AMT in deeper layers.The findings include the identification of two sliding zones in the Deda I landslide,the shallow sliding zone(DD-I-S1)depth is approximately 20 m,and the deep sliding zone(DD-I-S2)depth is 36.2-49.9 m.The sliding zone(DD-Ⅱ-S1)depth of the DedaⅡlandslide is 37.6-43.1 m.A novel MSM-based method for sliding zone identification is proposed,achieving less than 5%discrepancy in depth determination when compared with drilling data.These results provide a valuable reference for the spatial structural analysis of large-deepseated landslides in geologically complex regions like the eastern Tibetan Plateau.展开更多
Cross-sections for electronic excitation and de-excitation among the ground state and lowest-lying seven electronic excited states of carbon monoxide(CO)by low-energy electron impact are computed using the R-matrix me...Cross-sections for electronic excitation and de-excitation among the ground state and lowest-lying seven electronic excited states of carbon monoxide(CO)by low-energy electron impact are computed using the R-matrix method.The excitation cross-sections from the ground state to the electronic states a^(3)Π,a'^(3)Σ^(+)+and A^(1)Πagree with previous experimental and theoretical results.In addition,the cross-sections for the I^(1)Σ^(+)-and D^(1)Δstates of CO,which will cascade to CO a'^(3)Σ^(+)+and A^(1)Πstates,are calculated.Furthermore,in contrast to the typical increase in electronic excitation cross-sections with collision energy,the de-excitation cross-sections show a negative trend with increasing energy.展开更多
A numerical method is proposed to calculate the eigenvalues of the Zakharov–Shabat system based on Chebyshev polynomials. A mapping in the form of tanh(ax) is constructed according to the asymptotic of the potential ...A numerical method is proposed to calculate the eigenvalues of the Zakharov–Shabat system based on Chebyshev polynomials. A mapping in the form of tanh(ax) is constructed according to the asymptotic of the potential function for the Zakharov–Shabat eigenvalue problem. The mapping can distribute Chebyshev nodes very well considering the gradient for the potential function. Using Chebyshev polynomials, tanh(ax) mapping, and Chebyshev nodes, the Zakharov–Shabat eigenvalue problem is transformed into a matrix eigenvalue problem. This method has good convergence for the Satsuma–Yajima potential and the convergence rate is faster than the Fourier collocation method. This method is not only suitable for simple potential functions but also converges quickly for a complex Y-shape potential. It can also be further extended to other linear eigenvalue problems.展开更多
An anisotropic diffusion filter can be used to model a flow-dependent background error covariance matrix,which can be achieved by solving the advection-diffusion equation.Because of the directionality of the advection...An anisotropic diffusion filter can be used to model a flow-dependent background error covariance matrix,which can be achieved by solving the advection-diffusion equation.Because of the directionality of the advection term,the discrete method needs to be chosen very carefully.The finite analytic method is an alternative scheme to solve the advection-diffusion equation.As a combination of analytical and numerical methods,it not only has high calculation accuracy but also holds the characteristic of the auto upwind.To demonstrate its ability,the one-dimensional steady and unsteady advection-diffusion equation numerical examples are respectively solved by the finite analytic method.The more widely used upwind difference method is used as a control approach.The result indicates that the finite analytic method has higher accuracy than the upwind difference method.For the two-dimensional case,the finite analytic method still has a better performance.In the three-dimensional variational assimilation experiment,the finite analytic method can effectively improve analysis field accuracy,and its effect is significantly better than the upwind difference and the central difference method.Moreover,it is still a more effective solution method in the strong flow region where the advective-diffusion filter performs most prominently.展开更多
With the development of industrial activities,global warming has accelerated due to excessive emission of CO_(2).Enhanced Geothermal System(EGS)utilizes deep geothermal heat for power generation.Although porous medium...With the development of industrial activities,global warming has accelerated due to excessive emission of CO_(2).Enhanced Geothermal System(EGS)utilizes deep geothermal heat for power generation.Although porous medium theory is commonly employed to model geothermal reservoirs in EGS,Hot Dry Rock(HDR)presents a challenge as it consists of impermeable granite with zero porosity,potentially distorting the physical interpretation.To address this,the Lattice Boltzmann Method(LBM)is employed to simulate CO_(2)flow within geothermal reservoirs and the Finite Volume Method(FVM)to solve the energy conservation equation for temperature distribution.This combined method of LBM and FVM is imple-mented using MATLAB.The results showed that the Reynolds numbers(Re)of 3,000 and 8,000 lead to higher heat extraction rates from geothermal reservoirs.However,higher Re values may accelerate thermal breakthrough,posing challenges to EGS operation.Meanwhile,non-equilibrium of density in fractures becomes more pronounced during the system's life cycle,with non-Darcy's law becoming significant at Re values of 3,000 and 8,000.Density stratification due to buoyancy effects significantly impacts temperature distribution within geothermal reservoirs,with buoyancy effects at Re=100 under gravitational influence being noteworthy.Larger Re values(3,000 and 8,000)induce stronger forced convection,leading to more uniform density distribution.The addition of proppant negatively affects heat transfer performance in geothermal reservoirs,especially in single fractures.Practical engineering considerations should determine the quantity of proppant through detailed numerical simulations.展开更多
Deriving atmospheric parameters of a large sample of stars is of vital importance to understand the formation and evolution of the Milky Way.Photometric surveys,especially those with near-ultraviolet filters,can offer...Deriving atmospheric parameters of a large sample of stars is of vital importance to understand the formation and evolution of the Milky Way.Photometric surveys,especially those with near-ultraviolet filters,can offer accurate measurements of stellar parameters,with the precision comparable to that from low/medium resolution spectroscopy.In this study,we explore the capability of measuring stellar atmospheric parameters from Chinese Space Station Telescope(CSST)broad-band photometry(particularly in the near-ultraviolet bands),based on synthetic colors derived from model spectra.We find that colors from the optical and near-ultraviolet filter systems adopted by CSST show significant sensitivities to the stellar atmospheric parameters,especially the metallicity.According to our mock data tests,the precision of the photometric metallicity is quite high,with typical values of0.17 and 0.20 dex for dwarf and giant stars,respectively.The precision of the effective temperature estimated from broad-band colors are within 50 K.展开更多
A coupled numerical calculation method combining smooth particle hydrodynamics(SPH)and the finite element method(FEM)was implemented to investigate the seismic response of horizontal storage tanks.Anumericalmodel of a...A coupled numerical calculation method combining smooth particle hydrodynamics(SPH)and the finite element method(FEM)was implemented to investigate the seismic response of horizontal storage tanks.Anumericalmodel of a horizontal storage tank featuring a free liquid surface under seismic action was constructed using the SPH–FEM coupling method.The stored liquid was discretized using SPH particles,while the tank and supports were discretized using the FEM.The interaction between the stored liquid and the tank was simulated by using the meshless particle contact method.Then,the numerical simulation results were compared and analyzed against seismic simulation shaking table test data to validate the method.Subsequently,a series of numerical models,considering different liquid storage volumes and seismic effects,were constructed to obtain time history data of base shear and top center displacement,which revealed the seismic performance of horizontal storage tanks.Numerical simulation results and experimental data showed good agreement,with an error rate of less than 18.85%.And this conformity signifies the rationality of the SPH-FEM coupling method.The base shear and top center displacement values obtained by the coupled SPH-FEM method were only 53.3% to 69.1% of those calculated by the equivalent mass method employed in the current code.As the stored liquid volume increased,the seismic response of the horizontal storage tank exhibited a gradual upward trend,with the seismic response increasing from 73% to 388% for every 35% increase in stored liquid volume.The maximum von Mises stress of the tank and the supports remained below the steel yield strength during the earthquake.The coupled SPH-FEM method holds certain advantages in studying the seismic problems of tanks with complex structural forms,particularly due to the representation of the flow field distribution during earthquakes by involving reservoir fluid participation.展开更多
A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization...A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.展开更多
The two-component cold atom systems with anisotropic hopping amplitudes can be phenomenologically described by a two-dimensional Ising-XY coupled model with spatial anisotropy.At low temperatures,theoretical predictio...The two-component cold atom systems with anisotropic hopping amplitudes can be phenomenologically described by a two-dimensional Ising-XY coupled model with spatial anisotropy.At low temperatures,theoretical predictions[Phys.Rev.A 72053604(2005)]and[arXiv:0706.1609]indicate the existence of a topological ordered phase characterized by Ising and XY disorder but with 2XY ordering.However,due to ergodic difficulties faced by Monte Carlo methods at low temperatures,this topological phase has not been numerically explored.We propose a linear cluster updating Monte Carlo method,which flips spins without rejection in the anisotropy limit but does not change the energy.Using this scheme and conventional Monte Carlo methods,we succeed in revealing the nature of topological phases with half-vortices and domain walls.In the constructed global phase diagram,Ising and XY-type transitions are very close to each other and differ significantly from the schematic phase diagram reported earlier.We also propose and explore a wide range of quantities,including magnetism,superfluidity,specific heat,susceptibility,and even percolation susceptibility,and obtain consistent and reliable results.Furthermore,we observed first-order transitions characterized by common intersection points in magnetizations for different system sizes,as opposed to the conventional phase transition where Binder cumulants of various sizes share common intersections.The critical exponents of different types of phase transitions are reasonably fitted.The results are useful to help cold atom experiments explore the half-vortex topological phase.展开更多
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.展开更多
The initial shape of the secondary arc considerably influences its subsequent shape.To establish the model for the arcing time of the secondary arc and modify the single-phase reclosing sequence,theoretical and experi...The initial shape of the secondary arc considerably influences its subsequent shape.To establish the model for the arcing time of the secondary arc and modify the single-phase reclosing sequence,theoretical and experimental analysis of the evolution process of the short-circuit arc to the secondary arc is critical.In this study,an improved charge simulation method was used to develop the internal-space electric-field model of the short-circuit arc.The intensity of the electric field was used as an independent variable to describe the initial shape of the secondary arc.A secondary arc evolution model was developed based on this model.Moreover,the accuracy of the model was evaluated by comparison with physical experimental results.When the secondary arc current increased,the arcing time and dispersion increased.There is an overall trend of increasing arc length with increasing arcing time.Nevertheless,there is a reduction in arc length during arc ignition due to short circuits between the arc columns.Furthermore,the arcing time decreased in the range of 0°-90°as the angle between the wind direction and the x-axis increased.This work investigated the method by which short-circuit arcs evolve into secondary arcs.The results can be used to develop the secondary arc evolution model and to provide both a technical and theoretical basis for secondary arc suppression.展开更多
The optical potential ambiguity is a long-standing problem in the analysis of elastic scattering data.For a specific collid-ing system,ambiguous potential families can lead to different behaviors in the nearside and f...The optical potential ambiguity is a long-standing problem in the analysis of elastic scattering data.For a specific collid-ing system,ambiguous potential families can lead to different behaviors in the nearside and farside scattering components.By contrast,the envelope method can decompose the experimental data into two components with negative and positive deflection angles,respectively.Hence,a question arises as to whether the comparison between the calculated nearside(or farside)component and the derived positive-deflection-angle(or negative-deflection-angle)component can help analyze the potential ambiguity problem.In this study,we conducted a trial application of the envelope method to the potential ambiguity problem.The envelope method was improved by including uncertainties in the experimental data.The colliding systems of 16O+28Si at 215.2 MeV and 12C+12C at 1016 MeV were considered in the analyses.For each colliding system,the angular distribution experimental data were described nearly equally well by two potential sets,one of which is“surface transpar-ent”and the other is refractive.The calculated angular distributions were decomposed into nearside and farside scattering components.Using the improved envelope method,the experimental data were decomposed into the positive-deflection-angle and negative-deflection-angle components,which were then compared with the calculated nearside and farside components.The capability of the envelope method to analyze the potential ambiguities was also discussed.展开更多
Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridyna...Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.展开更多
The synthesis of α-calcium sulfate hemihydrate (α-CSH) from flue gas desulfurization (FGD)gypsum is a good way to realize the comprehensive utilization of FGD gypsum. To obtainα-CSH with the satisfactory performanc...The synthesis of α-calcium sulfate hemihydrate (α-CSH) from flue gas desulfurization (FGD)gypsum is a good way to realize the comprehensive utilization of FGD gypsum. To obtainα-CSH with the satisfactory performances, a facile hydrothermal-aging pretreatment process for FGD gypsum raw materials was proposed, where FGD gypsum was firstly hydrothermally converted to α-CSH whiskers, and α-CSH whiskers were further hydrated to synthesize CaSO4·2H2O (CSD) by aging under the regulation of N,N'-methylenebisacrylamide (MBA). The effects of aging time, MBA addition, aging temperature, and pH on the morphology of the synthesized CSD were investigated. The synthesized CSD crystals exhibit highly uniform prismatic morphology with the length of ca 100μm and the whiteness of 91.56%. The regulation mechanism of MBA was also illustrated. The synthesized CSD crystals with prismatic morphology were further used as raw materials to synthesize the short columnar α-CSH. The absolute dry compressive strength of paste prepared from the short columnar α-CSH is 40.85 MPa, which reaches α40 strength grade.展开更多
Mathematical modeling of the interaction between solar radiation and the Earth's atmosphere is formalized by the radiative transfer equation(RTE), whose resolution calls for two-stream approximations among other m...Mathematical modeling of the interaction between solar radiation and the Earth's atmosphere is formalized by the radiative transfer equation(RTE), whose resolution calls for two-stream approximations among other methods. This paper proposes a new two-stream approximation of the RTE with the development of the phase function and the intensity into a third-order series of Legendre polynomials. This new approach, which adds one more term in the expression of the intensity and the phase function, allows in the conditions of a plane parallel atmosphere a new mathematical formulation of γparameters. It is then compared to the Eddington, Hemispheric Constant, Quadrature, Combined Delta Function and Modified Eddington, and second-order approximation methods with reference to the Discrete Ordinate(Disort) method(δ –128 streams), considered as the most precise. This work also determines the conversion function of the proposed New Method using the fundamental definition of two-stream approximation(F-TSA) developed in a previous work. Notably,New Method has generally better precision compared to the second-order approximation and Hemispheric Constant methods. Compared to the Quadrature and Eddington methods, New Method shows very good precision for wide domains of the zenith angle μ 0, but tends to deviate from the Disort method with the zenith angle, especially for high values of optical thickness. In spite of this divergence in reflectance for high values of optical thickness, very strong correlation with the Disort method(R ≈ 1) was obtained for most cases of optical thickness in this study. An analysis of the Legendre polynomial series for simple functions shows that the high precision is due to the fact that the approximated functions ameliorate the accuracy when the order of approximation increases, although it has been proven that there is a limit order depending on the function from which the precision is lost. This observation indicates that increasing the order of approximation of the phase function of the RTE leads to a better precision in flux calculations. However, this approach may be limited to a certain order that has not been studied in this paper.展开更多
Food safety problems caused by excessive nitrite addition have been frequently reported and the detection of nitrite in food is particularly important. The standing time during the pretreatment of primary sample has a...Food safety problems caused by excessive nitrite addition have been frequently reported and the detection of nitrite in food is particularly important. The standing time during the pretreatment of primary sample has a great influence on the concentration of nitrite tested by spectrophotometric method. In this context, three kinds of food samples are prepared, including canned mustard, canned fish and home-made pickled water. A series of standing times are placed during the sample pretreatments and the corresponding nitrite contents in these samples are detected by spectrophotometric method based on N-ethylenediamine dihydrochloride. This study aims to find out a reasonable standing time during the pretreatment of food sample, providing influence factor for precise detection of nitrite.展开更多
This paper presents the results of a research project on the intercultural climate in schools and the academic success of immigrant students in Quebec.The intercultural climate in eight secondary public schools is doc...This paper presents the results of a research project on the intercultural climate in schools and the academic success of immigrant students in Quebec.The intercultural climate in eight secondary public schools is documented from the standpoint of various school and community stakeholders,students,and parents.More specifically,the paper focuses on the commitment to promoting cultural equity,which is identified through interviews with the targeted stakeholders on exclusionary practices and processes and the practices implemented to counter them.The results identify five areas of vulnerability among certain minoritized groups of students:the transition from welcoming to regular classes for those students learning French,evaluation and support for those experiencing learning or behavioral difficulties,the referral of these students to adult education programs or work-oriented training programs,their participation in extracurricular activities,and their access to special school projects.Nonetheless,some schools do report the introduction of many equity practices in response to the exclusionary ones.However,these initiatives are mostly occasional and voluntary and pertain to interactions between individuals.These findings call for a discussion on the importance of systemic and concerted initiatives towards achieving equity.展开更多
In connection with the historical context of social segregation and the assimilative aim of formal education imposed on them,the Indigenous Peoples of Québec have long been excluded from higher education.Today,ev...In connection with the historical context of social segregation and the assimilative aim of formal education imposed on them,the Indigenous Peoples of Québec have long been excluded from higher education.Today,even if their graduation rates are increasing,a persistent gap with the non-Indigenous population is maintained.The data in this article are drawn from 23 interviews with students and university graduates of the First Nations of Québec,as part of a thesis in educational sciences.Our paper analyzes how these students manage to combine the contributions of Indigenous education with those of Western education by developing paths that are part of the decolonization movement of education.Their relationship to Indigenous identity and cultures clearly influences their educational background and is analyzed taking into account power relations with the dominant culture.Our paper emphasizes more particularly the relationship to studies,the adaptation to the student profession,and the balance between studies,work,and family among Indigenous students,in a context of transition from university to multiversity which contributed to the increase of Indigenous university attendance.展开更多
The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent ...The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent singularities on the one-dimensional heat equation. The method is based on a Fourier decomposition of the solution and an extraction formula of the coefficients of the singularities coupled with a predictor-corrector algorithm. The method recovers the optimal convergence rate of the finite element method on a quasi-uniform mesh refinement. Numerical results are carried out to show the efficiency of the method.展开更多
The inversion of a non-singular square matrix applying a Computer Algebra System (CAS) is straightforward. The CASs make the numeric computation efficient but mock the mathematical characteristics. The algorithms cond...The inversion of a non-singular square matrix applying a Computer Algebra System (CAS) is straightforward. The CASs make the numeric computation efficient but mock the mathematical characteristics. The algorithms conducive to the output are sealed and inaccessible. In practice, other than the CPU timing, the applied inversion method is irrelevant. This research-oriented article discusses one such process, the Cayley-Hamilton (C.H.) [1]. Pursuing the process symbolically reveals its unpublished hidden mathematical characteristics even in the original article [1]. This article expands the general vision of the original named method without altering its practical applications. We have used the famous CAS Mathematica [2]. We have briefed the theory behind the method and applied it to different-sized symbolic and numeric matrices. The results are compared to the named CAS’s sealed, packaged library commands. The codes are given, and the algorithms are unsealed.展开更多
基金supported by the National Natural Science Foundation of China(42372339)the China Geological Survey Project(DD20221816,DD20190319)。
文摘It is of crucial importance to investigate the spatial structures of ancient landslides in the eastern Tibetan Plateau’s alpine canyons as they could provide valuable insights into the evolutionary history of the landslides and indicate the potential for future reactivation.This study examines the Deda ancient landslide,situated in the Chalong-ranbu fault zone,where creep deformation suggests a complex underground structure.By integrating remote sensing,field surveys,Audio-frequency Magnetotellurics(AMT),and Microtremor Survey Method(MSM)techniques,along with engineering geological drilling for validation,to uncover the landslide’s spatial feature s.The research indicates that a fault is developed in the upper part of the Deda ancient landslide,and the gully divides it into Deda landslide accumulation zoneⅠand Deda landslide accumulation zoneⅡin space.The distinctive geological characteristics detectable by MSM in the shallow subsurface and by AMT in deeper layers.The findings include the identification of two sliding zones in the Deda I landslide,the shallow sliding zone(DD-I-S1)depth is approximately 20 m,and the deep sliding zone(DD-I-S2)depth is 36.2-49.9 m.The sliding zone(DD-Ⅱ-S1)depth of the DedaⅡlandslide is 37.6-43.1 m.A novel MSM-based method for sliding zone identification is proposed,achieving less than 5%discrepancy in depth determination when compared with drilling data.These results provide a valuable reference for the spatial structural analysis of large-deepseated landslides in geologically complex regions like the eastern Tibetan Plateau.
基金Project supported by the National Natural Science Foundation of China (Grant No.11974253)。
文摘Cross-sections for electronic excitation and de-excitation among the ground state and lowest-lying seven electronic excited states of carbon monoxide(CO)by low-energy electron impact are computed using the R-matrix method.The excitation cross-sections from the ground state to the electronic states a^(3)Π,a'^(3)Σ^(+)+and A^(1)Πagree with previous experimental and theoretical results.In addition,the cross-sections for the I^(1)Σ^(+)-and D^(1)Δstates of CO,which will cascade to CO a'^(3)Σ^(+)+and A^(1)Πstates,are calculated.Furthermore,in contrast to the typical increase in electronic excitation cross-sections with collision energy,the de-excitation cross-sections show a negative trend with increasing energy.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.52171251,U2106225,and 52231011)Dalian Science and Technology Innovation Fund (Grant No.2022JJ12GX036)。
文摘A numerical method is proposed to calculate the eigenvalues of the Zakharov–Shabat system based on Chebyshev polynomials. A mapping in the form of tanh(ax) is constructed according to the asymptotic of the potential function for the Zakharov–Shabat eigenvalue problem. The mapping can distribute Chebyshev nodes very well considering the gradient for the potential function. Using Chebyshev polynomials, tanh(ax) mapping, and Chebyshev nodes, the Zakharov–Shabat eigenvalue problem is transformed into a matrix eigenvalue problem. This method has good convergence for the Satsuma–Yajima potential and the convergence rate is faster than the Fourier collocation method. This method is not only suitable for simple potential functions but also converges quickly for a complex Y-shape potential. It can also be further extended to other linear eigenvalue problems.
基金The National Key Research and Development Program of China under contract Nos 2022YFC3104804,2021YFC3101501,and 2017YFC1404103the National Programme on Global Change and Air-Sea Interaction of China under contract No.GASI-IPOVAI-04the National Natural Science Foundation of China under contract Nos 41876014,41606039,and 11801402.
文摘An anisotropic diffusion filter can be used to model a flow-dependent background error covariance matrix,which can be achieved by solving the advection-diffusion equation.Because of the directionality of the advection term,the discrete method needs to be chosen very carefully.The finite analytic method is an alternative scheme to solve the advection-diffusion equation.As a combination of analytical and numerical methods,it not only has high calculation accuracy but also holds the characteristic of the auto upwind.To demonstrate its ability,the one-dimensional steady and unsteady advection-diffusion equation numerical examples are respectively solved by the finite analytic method.The more widely used upwind difference method is used as a control approach.The result indicates that the finite analytic method has higher accuracy than the upwind difference method.For the two-dimensional case,the finite analytic method still has a better performance.In the three-dimensional variational assimilation experiment,the finite analytic method can effectively improve analysis field accuracy,and its effect is significantly better than the upwind difference and the central difference method.Moreover,it is still a more effective solution method in the strong flow region where the advective-diffusion filter performs most prominently.
基金supported by the Hebei Province Graduate Innovation Funding Project(CXZZBS2022029).
文摘With the development of industrial activities,global warming has accelerated due to excessive emission of CO_(2).Enhanced Geothermal System(EGS)utilizes deep geothermal heat for power generation.Although porous medium theory is commonly employed to model geothermal reservoirs in EGS,Hot Dry Rock(HDR)presents a challenge as it consists of impermeable granite with zero porosity,potentially distorting the physical interpretation.To address this,the Lattice Boltzmann Method(LBM)is employed to simulate CO_(2)flow within geothermal reservoirs and the Finite Volume Method(FVM)to solve the energy conservation equation for temperature distribution.This combined method of LBM and FVM is imple-mented using MATLAB.The results showed that the Reynolds numbers(Re)of 3,000 and 8,000 lead to higher heat extraction rates from geothermal reservoirs.However,higher Re values may accelerate thermal breakthrough,posing challenges to EGS operation.Meanwhile,non-equilibrium of density in fractures becomes more pronounced during the system's life cycle,with non-Darcy's law becoming significant at Re values of 3,000 and 8,000.Density stratification due to buoyancy effects significantly impacts temperature distribution within geothermal reservoirs,with buoyancy effects at Re=100 under gravitational influence being noteworthy.Larger Re values(3,000 and 8,000)induce stronger forced convection,leading to more uniform density distribution.The addition of proppant negatively affects heat transfer performance in geothermal reservoirs,especially in single fractures.Practical engineering considerations should determine the quantity of proppant through detailed numerical simulations.
基金the science research grants from the China Manned Space Project with No.CMS-CSST-2021-A08.Y.Hthe science research grants from the China Manned Space Project with No.CMS-CSST-2021-B05+3 种基金the NSFC for grant Nos.11903027 and11833006the NSFC for grant Nos.11973001,12090040,and 12090044the National Key R&D Program of China for grant No.2019YFA0405503.H.W.Zthe National Key R&D Program of China for grant No.2019YFA0405504。
文摘Deriving atmospheric parameters of a large sample of stars is of vital importance to understand the formation and evolution of the Milky Way.Photometric surveys,especially those with near-ultraviolet filters,can offer accurate measurements of stellar parameters,with the precision comparable to that from low/medium resolution spectroscopy.In this study,we explore the capability of measuring stellar atmospheric parameters from Chinese Space Station Telescope(CSST)broad-band photometry(particularly in the near-ultraviolet bands),based on synthetic colors derived from model spectra.We find that colors from the optical and near-ultraviolet filter systems adopted by CSST show significant sensitivities to the stellar atmospheric parameters,especially the metallicity.According to our mock data tests,the precision of the photometric metallicity is quite high,with typical values of0.17 and 0.20 dex for dwarf and giant stars,respectively.The precision of the effective temperature estimated from broad-band colors are within 50 K.
基金supported by Scientific Research Fund of Institute of Engineering Mechanics,China Earthquake Administration(Grant Nos.2021B06,2021C05)Heilongjiang Natural Science Foundation Joint Guidance Project(Grant No.LH2021E122).
文摘A coupled numerical calculation method combining smooth particle hydrodynamics(SPH)and the finite element method(FEM)was implemented to investigate the seismic response of horizontal storage tanks.Anumericalmodel of a horizontal storage tank featuring a free liquid surface under seismic action was constructed using the SPH–FEM coupling method.The stored liquid was discretized using SPH particles,while the tank and supports were discretized using the FEM.The interaction between the stored liquid and the tank was simulated by using the meshless particle contact method.Then,the numerical simulation results were compared and analyzed against seismic simulation shaking table test data to validate the method.Subsequently,a series of numerical models,considering different liquid storage volumes and seismic effects,were constructed to obtain time history data of base shear and top center displacement,which revealed the seismic performance of horizontal storage tanks.Numerical simulation results and experimental data showed good agreement,with an error rate of less than 18.85%.And this conformity signifies the rationality of the SPH-FEM coupling method.The base shear and top center displacement values obtained by the coupled SPH-FEM method were only 53.3% to 69.1% of those calculated by the equivalent mass method employed in the current code.As the stored liquid volume increased,the seismic response of the horizontal storage tank exhibited a gradual upward trend,with the seismic response increasing from 73% to 388% for every 35% increase in stored liquid volume.The maximum von Mises stress of the tank and the supports remained below the steel yield strength during the earthquake.The coupled SPH-FEM method holds certain advantages in studying the seismic problems of tanks with complex structural forms,particularly due to the representation of the flow field distribution during earthquakes by involving reservoir fluid participation.
基金supported by a Major Research Project in Higher Education Institutions in Henan Province,with Project Number 23A560015.
文摘A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.
基金Project supported by the Hefei National Research Center for Physical Sciences at the Microscale (Grant No.KF2021002)the Natural Science Foundation of Shanxi Province,China (Grant Nos.202303021221029 and 202103021224051)+2 种基金the National Natural Science Foundation of China (Grant Nos.11975024,12047503,and 12275263)the Anhui Provincial Supporting Program for Excellent Young Talents in Colleges and Universities (Grant No.gxyq ZD2019023)the National Key Research and Development Program of China (Grant No.2018YFA0306501)。
文摘The two-component cold atom systems with anisotropic hopping amplitudes can be phenomenologically described by a two-dimensional Ising-XY coupled model with spatial anisotropy.At low temperatures,theoretical predictions[Phys.Rev.A 72053604(2005)]and[arXiv:0706.1609]indicate the existence of a topological ordered phase characterized by Ising and XY disorder but with 2XY ordering.However,due to ergodic difficulties faced by Monte Carlo methods at low temperatures,this topological phase has not been numerically explored.We propose a linear cluster updating Monte Carlo method,which flips spins without rejection in the anisotropy limit but does not change the energy.Using this scheme and conventional Monte Carlo methods,we succeed in revealing the nature of topological phases with half-vortices and domain walls.In the constructed global phase diagram,Ising and XY-type transitions are very close to each other and differ significantly from the schematic phase diagram reported earlier.We also propose and explore a wide range of quantities,including magnetism,superfluidity,specific heat,susceptibility,and even percolation susceptibility,and obtain consistent and reliable results.Furthermore,we observed first-order transitions characterized by common intersection points in magnetizations for different system sizes,as opposed to the conventional phase transition where Binder cumulants of various sizes share common intersections.The critical exponents of different types of phase transitions are reasonably fitted.The results are useful to help cold atom experiments explore the half-vortex topological phase.
基金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.
基金supported by National Natural Science Foundation of China(Nos.92066108 and 51277061)。
文摘The initial shape of the secondary arc considerably influences its subsequent shape.To establish the model for the arcing time of the secondary arc and modify the single-phase reclosing sequence,theoretical and experimental analysis of the evolution process of the short-circuit arc to the secondary arc is critical.In this study,an improved charge simulation method was used to develop the internal-space electric-field model of the short-circuit arc.The intensity of the electric field was used as an independent variable to describe the initial shape of the secondary arc.A secondary arc evolution model was developed based on this model.Moreover,the accuracy of the model was evaluated by comparison with physical experimental results.When the secondary arc current increased,the arcing time and dispersion increased.There is an overall trend of increasing arc length with increasing arcing time.Nevertheless,there is a reduction in arc length during arc ignition due to short circuits between the arc columns.Furthermore,the arcing time decreased in the range of 0°-90°as the angle between the wind direction and the x-axis increased.This work investigated the method by which short-circuit arcs evolve into secondary arcs.The results can be used to develop the secondary arc evolution model and to provide both a technical and theoretical basis for secondary arc suppression.
基金This work was supported by the National Natural Science Foundation of China(Nos.12005047 and U1832105).
文摘The optical potential ambiguity is a long-standing problem in the analysis of elastic scattering data.For a specific collid-ing system,ambiguous potential families can lead to different behaviors in the nearside and farside scattering components.By contrast,the envelope method can decompose the experimental data into two components with negative and positive deflection angles,respectively.Hence,a question arises as to whether the comparison between the calculated nearside(or farside)component and the derived positive-deflection-angle(or negative-deflection-angle)component can help analyze the potential ambiguity problem.In this study,we conducted a trial application of the envelope method to the potential ambiguity problem.The envelope method was improved by including uncertainties in the experimental data.The colliding systems of 16O+28Si at 215.2 MeV and 12C+12C at 1016 MeV were considered in the analyses.For each colliding system,the angular distribution experimental data were described nearly equally well by two potential sets,one of which is“surface transpar-ent”and the other is refractive.The calculated angular distributions were decomposed into nearside and farside scattering components.Using the improved envelope method,the experimental data were decomposed into the positive-deflection-angle and negative-deflection-angle components,which were then compared with the calculated nearside and farside components.The capability of the envelope method to analyze the potential ambiguities was also discussed.
文摘Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.
基金Funded by National Natural Science Foundation of China(No.22008049)Natural Science Foundation of Hebei Province,China (Nos.B2020202081 and B2018202330)+1 种基金Key Laboratory of Gas Hydrate,Guangzhou Institute of Energy Conversion,Chinese Academy of Sciences,China (No.E029kf1601)Research Fund Program of Science and Technology of Colleges and Universities of Hebei Province,China (No.QN2019012)。
文摘The synthesis of α-calcium sulfate hemihydrate (α-CSH) from flue gas desulfurization (FGD)gypsum is a good way to realize the comprehensive utilization of FGD gypsum. To obtainα-CSH with the satisfactory performances, a facile hydrothermal-aging pretreatment process for FGD gypsum raw materials was proposed, where FGD gypsum was firstly hydrothermally converted to α-CSH whiskers, and α-CSH whiskers were further hydrated to synthesize CaSO4·2H2O (CSD) by aging under the regulation of N,N'-methylenebisacrylamide (MBA). The effects of aging time, MBA addition, aging temperature, and pH on the morphology of the synthesized CSD were investigated. The synthesized CSD crystals exhibit highly uniform prismatic morphology with the length of ca 100μm and the whiteness of 91.56%. The regulation mechanism of MBA was also illustrated. The synthesized CSD crystals with prismatic morphology were further used as raw materials to synthesize the short columnar α-CSH. The absolute dry compressive strength of paste prepared from the short columnar α-CSH is 40.85 MPa, which reaches α40 strength grade.
文摘Mathematical modeling of the interaction between solar radiation and the Earth's atmosphere is formalized by the radiative transfer equation(RTE), whose resolution calls for two-stream approximations among other methods. This paper proposes a new two-stream approximation of the RTE with the development of the phase function and the intensity into a third-order series of Legendre polynomials. This new approach, which adds one more term in the expression of the intensity and the phase function, allows in the conditions of a plane parallel atmosphere a new mathematical formulation of γparameters. It is then compared to the Eddington, Hemispheric Constant, Quadrature, Combined Delta Function and Modified Eddington, and second-order approximation methods with reference to the Discrete Ordinate(Disort) method(δ –128 streams), considered as the most precise. This work also determines the conversion function of the proposed New Method using the fundamental definition of two-stream approximation(F-TSA) developed in a previous work. Notably,New Method has generally better precision compared to the second-order approximation and Hemispheric Constant methods. Compared to the Quadrature and Eddington methods, New Method shows very good precision for wide domains of the zenith angle μ 0, but tends to deviate from the Disort method with the zenith angle, especially for high values of optical thickness. In spite of this divergence in reflectance for high values of optical thickness, very strong correlation with the Disort method(R ≈ 1) was obtained for most cases of optical thickness in this study. An analysis of the Legendre polynomial series for simple functions shows that the high precision is due to the fact that the approximated functions ameliorate the accuracy when the order of approximation increases, although it has been proven that there is a limit order depending on the function from which the precision is lost. This observation indicates that increasing the order of approximation of the phase function of the RTE leads to a better precision in flux calculations. However, this approach may be limited to a certain order that has not been studied in this paper.
文摘Food safety problems caused by excessive nitrite addition have been frequently reported and the detection of nitrite in food is particularly important. The standing time during the pretreatment of primary sample has a great influence on the concentration of nitrite tested by spectrophotometric method. In this context, three kinds of food samples are prepared, including canned mustard, canned fish and home-made pickled water. A series of standing times are placed during the sample pretreatments and the corresponding nitrite contents in these samples are detected by spectrophotometric method based on N-ethylenediamine dihydrochloride. This study aims to find out a reasonable standing time during the pretreatment of food sample, providing influence factor for precise detection of nitrite.
文摘This paper presents the results of a research project on the intercultural climate in schools and the academic success of immigrant students in Quebec.The intercultural climate in eight secondary public schools is documented from the standpoint of various school and community stakeholders,students,and parents.More specifically,the paper focuses on the commitment to promoting cultural equity,which is identified through interviews with the targeted stakeholders on exclusionary practices and processes and the practices implemented to counter them.The results identify five areas of vulnerability among certain minoritized groups of students:the transition from welcoming to regular classes for those students learning French,evaluation and support for those experiencing learning or behavioral difficulties,the referral of these students to adult education programs or work-oriented training programs,their participation in extracurricular activities,and their access to special school projects.Nonetheless,some schools do report the introduction of many equity practices in response to the exclusionary ones.However,these initiatives are mostly occasional and voluntary and pertain to interactions between individuals.These findings call for a discussion on the importance of systemic and concerted initiatives towards achieving equity.
文摘In connection with the historical context of social segregation and the assimilative aim of formal education imposed on them,the Indigenous Peoples of Québec have long been excluded from higher education.Today,even if their graduation rates are increasing,a persistent gap with the non-Indigenous population is maintained.The data in this article are drawn from 23 interviews with students and university graduates of the First Nations of Québec,as part of a thesis in educational sciences.Our paper analyzes how these students manage to combine the contributions of Indigenous education with those of Western education by developing paths that are part of the decolonization movement of education.Their relationship to Indigenous identity and cultures clearly influences their educational background and is analyzed taking into account power relations with the dominant culture.Our paper emphasizes more particularly the relationship to studies,the adaptation to the student profession,and the balance between studies,work,and family among Indigenous students,in a context of transition from university to multiversity which contributed to the increase of Indigenous university attendance.
文摘The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent singularities on the one-dimensional heat equation. The method is based on a Fourier decomposition of the solution and an extraction formula of the coefficients of the singularities coupled with a predictor-corrector algorithm. The method recovers the optimal convergence rate of the finite element method on a quasi-uniform mesh refinement. Numerical results are carried out to show the efficiency of the method.
文摘The inversion of a non-singular square matrix applying a Computer Algebra System (CAS) is straightforward. The CASs make the numeric computation efficient but mock the mathematical characteristics. The algorithms conducive to the output are sealed and inaccessible. In practice, other than the CPU timing, the applied inversion method is irrelevant. This research-oriented article discusses one such process, the Cayley-Hamilton (C.H.) [1]. Pursuing the process symbolically reveals its unpublished hidden mathematical characteristics even in the original article [1]. This article expands the general vision of the original named method without altering its practical applications. We have used the famous CAS Mathematica [2]. We have briefed the theory behind the method and applied it to different-sized symbolic and numeric matrices. The results are compared to the named CAS’s sealed, packaged library commands. The codes are given, and the algorithms are unsealed.