In this study,the structural characters,antioxidant activities and bile acid-binding ability of sea buckthorn polysaccharides(HRPs)obtained by the commonly used hot water(HRP-W),pressurized hot water(HRP-H),ultrasonic...In this study,the structural characters,antioxidant activities and bile acid-binding ability of sea buckthorn polysaccharides(HRPs)obtained by the commonly used hot water(HRP-W),pressurized hot water(HRP-H),ultrasonic(HRP-U),acid(HRP-C)and alkali(HRP-A)assisted extraction methods were investigated.The results demonstrated that extraction methods had significant effects on extraction yield,monosaccharide composition,molecular weight,particle size,triple-helical structure,and surface morphology of HRPs except for the major linkage bands.Thermogravimetric analysis showed that HRP-U with filamentous reticular microstructure exhibited better thermal stability.The HRP-A with the lowest molecular weight and highest arabinose content possessed the best antioxidant activities.Moreover,the rheological analysis indicated that HRPs with higher galacturonic acid content and molecular weight showed higher viscosity and stronger crosslinking network(HRP-C,HRP-W and HRP-U),which exhibited stronger bile acid binding capacity.The present findings provide scientific evidence in the preparation technology of sea buckthorn polysaccharides with good antioxidant and bile acid binding capacity which are related to the structure affected by the extraction methods.展开更多
The nomenclature for compounds that are modified with isotopes is growing every day. Compounds can be modified with isotopes either individually, in a functional group or groups, or completely with all atomic centers ...The nomenclature for compounds that are modified with isotopes is growing every day. Compounds can be modified with isotopes either individually, in a functional group or groups, or completely with all atomic centers of the element. This diversity of isotope-modified compounds increases the range of researches that can be studied using them. Compounds modified with isotopes of carbon-13 or nitrogen-15 can be converted into carbon monoxide, carbon dioxide and molecular nitrogen. Currently, only the average value of carbon-13 or nitrogen-15 isotopes can be determined. However, by directly determining the atomic share of these isotopes in organic compounds modified with isotopes, information about the isotopic centers of the element can be obtained. The atomic fraction of an element is defined as a single carbon or nitrogen isotope-modified center or centers, or all centers that are isotope-modified with that element at the same time. Carbon-13 or nitrogen-15 isotopes’ atomic fraction can be determined molecularly or with fragment ions of different elemental content, or both. This makes the method self-verifying, increasing the accuracy and reliability of the results obtained. Amino acids, such as asparagine, aspartic acid, methionine, and threonine, are essential for the human body. This proposed method of isotopic analysis will increase the possibilities for scientific research using these compounds.展开更多
The inverse and direct piezoelectric and circuit coupling are widely observed in advanced electro-mechanical systems such as piezoelectric energy harvesters.Existing strongly coupled analysis methods based on direct n...The inverse and direct piezoelectric and circuit coupling are widely observed in advanced electro-mechanical systems such as piezoelectric energy harvesters.Existing strongly coupled analysis methods based on direct numerical modeling for this phenomenon can be classified into partitioned or monolithic formulations.Each formulation has its advantages and disadvantages,and the choice depends on the characteristics of each coupled problem.This study proposes a new option:a coupled analysis strategy that combines the best features of the existing formulations,namely,the hybrid partitioned-monolithic method.The analysis of inverse piezoelectricity and the monolithic analysis of direct piezoelectric and circuit interaction are strongly coupled using a partitioned iterative hierarchical algorithm.In a typical benchmark problem of a piezoelectric energy harvester,this research compares the results from the proposed method to those from the conventional strongly coupled partitioned iterative method,discussing the accuracy,stability,and computational cost.The proposed hybrid concept is effective for coupled multi-physics problems,including various coupling conditions.展开更多
Analyses of stable isotopes(C,O,H)in tree rings are increasingly important cross-disciplinary programs.The rapid development in this field documented in an increasing number of publications requires a comprehensive re...Analyses of stable isotopes(C,O,H)in tree rings are increasingly important cross-disciplinary programs.The rapid development in this field documented in an increasing number of publications requires a comprehensive review.This study includes a bibliometric analysis-based review to better understand research trends in tree ring stable isotope research.Overall,1475 publications were selected from the Web of Science Core Collection for 1974-2023.The findings are that:(1)numbers of annual publications and citations increased since 1974.From 1974 to 1980,there were around two relevant publications per year.However,from 2020 to 2022,this rose sharply to 109 publications per year.Likewise,average article citations were less than four per year before 1990,but were around four per article per year after 2000;(2)the major subjects using tree ring stable isotopes include forestry,geosciences,and environmental sciences,contributing to 42.5%of the total during 1974-2023;(3)the top three most productive institutions are the Chinese Academy of Sciences(423),the Swiss Federal Institute for Forest,Snow and Landscape Research(227),and the University of Arizona(204).These achievements result from strong collaborations;(4)review papers,for example,(Dawson et al.,Annu Rev Ecol Syst 33:507-559,2002)and(McCarroll and Loader,Quat Sci Rev 23:771-801,2004),are among the most cited,with more than 1000 citations;(5)tree ring stable isotope studies mainly focus on climatology and ecology,with atmospheric CO_(2) one of the most popular topics.Since 2010,precipitation and drought have received increasing attention.Based on this analysis,the research stages,key findings,debated issues,limitations and direc-tions for future research are summarized.This study serves as an important attempt to understand the progress on the use of stable isotopes in tree rings,providing scientific guid-ance for young researchers in this field.展开更多
Geochemistry, zircon U–Pb geochronology, and Hf isotope data for the Early Paleozoic granites in the Baoshan Block reveal the Early Paleozoic tectonic evolution of the Proto-Tethys. The samples are high-K, calcalkali...Geochemistry, zircon U–Pb geochronology, and Hf isotope data for the Early Paleozoic granites in the Baoshan Block reveal the Early Paleozoic tectonic evolution of the Proto-Tethys. The samples are high-K, calcalkaline, strongly peraluminous rocks with A/CNK values of 1.37–1.46, are enriched in SiO2, K2O, and Rb, and are depleted in Nb, P, Ti, Eu, and heavy rare earth elements,which indicates the crystallization fractionation of the granitic magma. Zircon U–Pb dating indicates that they formed in ca. 480 Ma. The Nansa granites have εHf(t) values ranging from-16.04 to 4.36 with corresponding TC DMages of 2.10–0.81 Ga, which suggests the magmas derived from the partial melting of ancient metasedimentary with minor involvement of mantle-derived components. A synthesis of data for the Early Paleozoic igneous rocks in the Baoshan block and adjacent(Tengchong,Qiangtang, Sibumasu, Himalaya, etc.) blocks indicates that these blocks were all aligned along the proto-Tethyan margin of East Gondwana in the Early Paleozoic. The Early Paleozoic S-type granites from Nansa were generated in a high-temperature and low-pressure(HTLP) extensional tectonic setting, which resulted from Andean-type orogeny instead of the final assembly of Gondwana or crustal extension in a non-arc environment. In certain places, an expanding environment may exist in opposition to the tectonic backdrop of the lithosphere’s thickening and shortening, leading the crust to melt and decompress,mantle-derived materials to mix, and a small quantity of peraluminous granite to emerge.展开更多
Additive Runge-Kutta methods designed for preserving highly accurate solutions in mixed-precision computation were previously proposed and analyzed.These specially designed methods use reduced precision for the implic...Additive Runge-Kutta methods designed for preserving highly accurate solutions in mixed-precision computation were previously proposed and analyzed.These specially designed methods use reduced precision for the implicit computations and full precision for the explicit computations.In this work,we analyze the stability properties of these methods and their sensitivity to the low-precision rounding errors,and demonstrate their performance in terms of accuracy and efficiency.We develop codes in FORTRAN and Julia to solve nonlinear systems of ODEs and PDEs using the mixed-precision additive Runge-Kutta(MP-ARK)methods.The convergence,accuracy,and runtime of these methods are explored.We show that for a given level of accuracy,suitably chosen MP-ARK methods may provide significant reductions in runtime.展开更多
Condensed and hydrolysable tannins are non-toxic natural polyphenols that are a commercial commodity industrialized for tanning hides to obtain leather and for a growing number of other industrial applications mainly ...Condensed and hydrolysable tannins are non-toxic natural polyphenols that are a commercial commodity industrialized for tanning hides to obtain leather and for a growing number of other industrial applications mainly to substitute petroleum-based products.They are a definite class of sustainable materials of the forestry industry.They have been in operation for hundreds of years to manufacture leather and now for a growing number of applications in a variety of other industries,such as wood adhesives,metal coating,pharmaceutical/medical applications and several others.This review presents the main sources,either already or potentially commercial of this forestry by-materials,their industrial and laboratory extraction systems,their systems of analysis with their advantages and drawbacks,be these methods so simple to even appear primitive but nonetheless of proven effectiveness,or very modern and instrumental.It constitutes a basic but essential summary of what is necessary to know of these sustainable materials.In doing so,the review highlights some of the main challenges that remain to be addressed to deliver the quality and economics of tannin supply necessary to fulfill the industrial production requirements for some materials-based uses.展开更多
It is a challenge to make thorough but efficient experimental designs for the coupled mineral dissolution and precipitation studies in a multi-mineral system, because it is difficult to speculate the best experimental...It is a challenge to make thorough but efficient experimental designs for the coupled mineral dissolution and precipitation studies in a multi-mineral system, because it is difficult to speculate the best experimental duration, optimal sampling schedule, effects of different experimental conditions, and how to maximize the experimental outputs prior to the actual experiments. Geochemical modeling is an efficient and effective tool to assist the experimental design by virtually running all scenarios of interest for the studied system and predicting the experimental outcomes. Here we demonstrated an example of geochemical modeling assisted experimental design of coupled labradorite dissolution and calcite and clayey mineral precipitation using multiple isotope tracers. In this study, labradorite(plagioclase) was chosen as the reactant because it is both a major component and one of the most reactive minerals in basalt. Following our isotope doping studies of single minerals in the last ten years, initial solutions in the simulations were doped withmultiple isotopes(e.g., Ca and Si). Geochemical modeling results show that the use of isotope tracers gives us orders of magnitude more sensitivity than the conventional method based on concentrations and allows us to decouple dissolution and precipitation reactions at near-equilibrium condition. The simulations suggest that the precise unidirectional dissolution rates can inform us which rate laws plagioclase dissolution has followed. Calcite precipitation occurred at near-equilibrium and the multiple isotope tracer experiments would provide near-equilibrium precipitation rates, which was a challenge for the conventional concentration-based experiments. In addition, whether the precipitation of clayey phases is the rate-limiting step in some multi-mineral systems will be revealed. Overall, the modeling results of multimineral reaction kinetics will improve the understanding of the coupled dissolution–precipitation in the multi-mineral systems and the quality of geochemical modeling prediction of CO_(2) removal and storage efficacy in the basalt systems.展开更多
In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-depe...In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-dependent problems.We use the convex splitting method,the variant energy quadratization method,and the scalar auxiliary variable method coupled with the LDG method to construct first-order temporal accurate schemes based on the gradient flow structure of the models.These semi-implicit schemes are decoupled,energy stable,and can be extended to high accuracy schemes using the semi-implicit spectral deferred correction method.Many bound preserving DG discretizations are only worked on explicit time integration methods and are difficult to get high-order accuracy.To overcome these difficulties,we use the Lagrange multipliers to enforce the implicit or semi-implicit LDG schemes to satisfy the bound constraints at each time step.This bound preserving limiter results in the Karush-Kuhn-Tucker condition,which can be solved by an efficient active set semi-smooth Newton method.Various numerical experiments illustrate the high-order accuracy and the effect of bound preserving.展开更多
The collocation method is a widely used numerical method for science and engineering problems governed by partial differential equations.This paper provides a comprehensive review of collocation methods and their appl...The collocation method is a widely used numerical method for science and engineering problems governed by partial differential equations.This paper provides a comprehensive review of collocation methods and their applications,focused on elasticity,heat conduction,electromagnetic field analysis,and fluid dynamics.The merits of the collocation method can be attributed to the need for element mesh,simple implementation,high computational efficiency,and ease in handling irregular domain problems since the collocation method is a type of node-based numerical method.Beginning with the fundamental principles of the collocation method,the discretization process in the continuous domain is elucidated,and how the collocation method approximation solutions for solving differential equations are explained.Delving into the historical development of the collocation methods,their earliest applications and key milestones are traced,thereby demonstrating their evolution within the realm of numerical computation.The mathematical foundations of collocation methods,encompassing the selection of interpolation functions,definition of weighting functions,and derivation of integration rules,are examined in detail,emphasizing their significance in comprehending the method’s effectiveness and stability.At last,the practical application of the collocation methods in engineering contexts is emphasized,including heat conduction simulations,electromagnetic coupled field analysis,and fluid dynamics simulations.These specific case studies can underscore collocation method’s broad applicability and effectiveness in addressing complex engineering challenges.In conclusion,this paper puts forward the future development trend of the collocation method through rigorous analysis and discussion,thereby facilitating further advancements in research and practical applications within these fields.展开更多
It is well-known that the equilibrium isotope fractionation of mercury(Hg)includes classical massdependent fractionations(MDFs)and nuclear volume effect(NVE)induced mass-independent fractionations(MIFs).However,the ef...It is well-known that the equilibrium isotope fractionation of mercury(Hg)includes classical massdependent fractionations(MDFs)and nuclear volume effect(NVE)induced mass-independent fractionations(MIFs).However,the effect of the NVE on these kinetic processes is not known.The total fractionations(MDFs+NVEinduced MIFs)of several representative Hg-incorporated substances were selected and calculated with ab initio calculations in this work for both equilibrium and kinetic processes.NVE-induced MIFs were calculated with scaled contact electron densities at the nucleus through systematic evaluations of their accuracy and errors using the Gaussian09 and DIRAC19 packages(named the electron density scaling method).Additionally,the NVE-induced kinetic isotope effect(KIE)of Hg isotopes are also calculated with this method for several representative Hg oxidation reactions by chlorine species.Total KIEs for 202 Hg/^(198)Hg ranging from−2.27‰to 0.96‰are obtained.Three anomalous^(202)Hg-enriched KIEs(δ^(202)Hg/^(198)Hg=0.83‰,0.94‰,and 0.96‰,)caused by the NVE are observed,which are quite different from the classical view(i.e.,light isotopes react faster than the heavy ones).The electron density scaling method we developed in this study can provide an easier way to calculate the NVE-induced KIEs for heavy isotopes and serve to better understand the fractionation mechanisms of mercury isotope systems.展开更多
Understanding the relationship between forest management and water use efficiency(WUE)is important for evaluating forest adaptability to climate change.However,the effects of thinning and understory removal on WUE and...Understanding the relationship between forest management and water use efficiency(WUE)is important for evaluating forest adaptability to climate change.However,the effects of thinning and understory removal on WUE and its key controlling processes are not well understood,which limits our comprehension of the physiological mechanisms of various management practices.In this study,four forest management measures(no thinning:NT;understory removal:UR;light thinning:LT;and heavy thinning:HT)were carried out in Pinus massoniana plantations in a subtropical region of China.Photosynthetic capacity and needle stable carbon isotope composition(δ^(13)C)were measured to assess instantaneous water use efficiency(WUE_(inst))and long-term water use efficiency(WUE_(i)).Multiple regression models and structural equation modelling(SEM)identified the effects of soil properties and physiological performances on WUE_(inst)and WUE_(i).The results show that WUE_(inst)values among the four treatments were insignificant.However,compared with the NT stand(35.8μmol·mol^(-1)),WUE_(i)values significantly increased to 41.7μmol·mol^(-1)in the UR,50.1μmol·mol^(-1)in the LT and 46.6μmol·mol^(-1)in HT treatments,largely explained by photosynthetic capacity and soil water content.Understory removal did not change physiological performance(needle water potential and photosynthetic capacity).Thinning increased the net photosynthetic rate(A_n)but not stomatal conductance(g_s)or predawn needle water potential(ψ_(pd)),implying that the improvement in water use efficiency for thinned stands was largely driven by radiation interception than by soil water availability.In general,thinning may be an appropriate management measure to promote P.massoniana WUE to cope with seasonal droughts under future extreme climates.展开更多
Hydrogen is the new age alternative energy source to combat energy demand and climate change.Storage of hydrogen is vital for a nation’s growth.Works of literature provide different methods for storing the produced h...Hydrogen is the new age alternative energy source to combat energy demand and climate change.Storage of hydrogen is vital for a nation’s growth.Works of literature provide different methods for storing the produced hydrogen,and the rational selection of a viable method is crucial for promoting sustainability and green practices.Typically,hydrogen storage is associated with diverse sustainable and circular economy(SCE)criteria.As a result,the authors consider the situation a multi-criteria decision-making(MCDM)problem.Studies infer that previous models for hydrogen storage method(HSM)selection(i)do not consider preferences in the natural language form;(ii)weights of experts are not methodically determined;(iii)hesitation of experts during criteria weight assessment is not effectively explored;and(iv)three-stage solution of a suitable selection of HSM is unexplored.Driven by these gaps,in this paper,authors put forward a new integrated framework,which considers double hierarchy linguistic information for rating,criteria importance through inter-criteria correlation(CRITIC)for expert weight calculation,evidence-based Bayesian method for criteria weight estimation,and combined compromise solution(CoCoSo)for ranking HSMs.The applicability of the developed framework is testified by using a case example of HSM selection in India.Sensitivity and comparative analysis reveal the merits and limitations of the developed framework.展开更多
By studying the light isotopic compositions of carbon,oxygen,and hydrogen,combined with previous research results on the ore-forming source of the deposit,the authors try to uncover its metallogenic origin.The δ^(18)...By studying the light isotopic compositions of carbon,oxygen,and hydrogen,combined with previous research results on the ore-forming source of the deposit,the authors try to uncover its metallogenic origin.The δ^(18)O and δ^(13)C isotope signatures of dolomite samples vary between 10.2 and 13.0‰,and between−7.2 and−5.2‰,respectively,implying that the carbon derives from the upper mantle.δD and δ^(18) O of quartz,biotite,and muscovite from diff erent ore veins of the deposit vary between−82 and−59‰,and between 11.6 and 12.4‰,respectively,implying that the metallogenic solutions are mainly magmatic.According to the relevant research results of many isotope geologists,the fractionation degree of hydrogen isotopes increases as the depth to the Earth’s core increases,and the more diff erentiated the hydrogen isotopes are,the lower their values will be.In other words,mantle-derived solutions can have extremely low hydrogen isotope values.This means that the δD‰ value−134 of the pyrrhotite sample numbered SD-34 in this article may indicate mantle-derived oreforming fl uid of the deposit.The formation of the Dashuigou tellurium deposit occurred between 91.71 and 80.19 Ma.展开更多
The Nianzi granite unit,which includes the Nianzi,Xiaolianghou and Xiawopu granitic intrusions,is a significant component of the northern part of the North China Craton(NCC)and is situated in the Yanshan fold and thru...The Nianzi granite unit,which includes the Nianzi,Xiaolianghou and Xiawopu granitic intrusions,is a significant component of the northern part of the North China Craton(NCC)and is situated in the Yanshan fold and thrust belt(YFTB).However,there is still debate regarding the tectonic evolutionary history of the YFTB during the late Permian to Triassic period,specifically regarding the timing of subduction and collision between the NCC and the Paleo-Asian Ocean.The Nianzi granite unit exhibits unique petrological,geochronological and geochemical signatures that shed light on the tectonic evolutionary history of the YFTB.This study presents detailed petrology,whole-rock geochemistry,together with Sr-Nd isotopic,zircon U-Pb dating and Lu-Hf isotopic data of the granites within the Nianzi granite unit.Our findings demonstrate that the granites primarily consist of subhedral K-feldspar,plagioclase,quartz,minor biotite and hornblende,with accessory titanite,apatite,magnetite and zircon.Zircon U-Pb dating indicates that the Xiaolianghou granite was emplaced at 247.5±0.62 Ma.Additionally,the adakitic characteristics of the Nianzi,Xiawopu and Xiaolianghou granitic intrusions,such as high Sr and Ba contents and high ratios of Sr/Y and(La/Yb)N,combined with negative Sr-Nd and Lu-Hf isotopes(87Sr/86Sr)i=0.705681–0.7057433,εNd(t)=−21.98 to−20.97,zirconεHf(t)=−20.26 to−9.92,as well as the I-type granite features of high SiO_(2),Na_(2)O and K_(2)O/Na_(2)O ratios,enriched Rb,K,Sr and Ba,along with depleted Th,U,Nb,Ta,P and Ti,suggest that the Nianzi granitic unit was mainly derived from the partial melting of a thickened lower crust containing hydrous,calc-alkaline to high-K calc-alkaline,mafic to intermediate metamorphic rocks.In light of these parameters,we further integrate our data with previous studies and conclude that the Nianzi granitic unit was generated in a post-collisional extensional environment during the Early Triassic.展开更多
In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal e...In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal energy are positive,and the mass fraction of each species is between 0 and 1.Second,due to the rapid reaction rate,the system may contain stiff sources,and the strong-stability-preserving explicit Runge-Kutta method may result in limited time-step sizes.To obtain physically relevant numerical approximations,we apply the bound-preserving technique to the DG methods.Though traditional positivity-preserving techniques can successfully yield positive density,internal energy,and mass fractions,they may not enforce the upper bound 1 of the mass fractions.To solve this problem,we need to(i)make sure the numerical fluxes in the equations of the mass fractions are consistent with that in the equation of the density;(ii)choose conservative time integrations,such that the summation of the mass fractions is preserved.With the above two conditions,the positive mass fractions have summation 1,and then,they are all between 0 and 1.For time discretization,we apply the modified Runge-Kutta/multi-step Patankar methods,which are explicit for the flux while implicit for the source.Such methods can handle stiff sources with relatively large time steps,preserve the positivity of the target variables,and keep the summation of the mass fractions to be 1.Finally,it is not straightforward to combine the bound-preserving DG methods and the Patankar time integrations.The positivity-preserving technique for DG methods requires positive numerical approximations at the cell interfaces,while Patankar methods can keep the positivity of the pre-selected point values of the target variables.To match the degree of freedom,we use polynomials on rectangular meshes for problems in two space dimensions.To evolve in time,we first read the polynomials at the Gaussian points.Then,suitable slope limiters can be applied to enforce the positivity of the solutions at those points,which can be preserved by the Patankar methods,leading to positive updated numerical cell averages.In addition,we use another slope limiter to get positive solutions used for the bound-preserving technique for the flux.Numerical examples are given to demonstrate the good performance of the proposed schemes.展开更多
Hydrogen storage alloys(HSAs)are attracting widespread interest in the nuclear industry because of the generation of stable metal hydrides after tritium absorption,thus effectively preventing the leakage of radioactiv...Hydrogen storage alloys(HSAs)are attracting widespread interest in the nuclear industry because of the generation of stable metal hydrides after tritium absorption,thus effectively preventing the leakage of radioactive tritium.Commonly used HSAs in the hydrogen isotopes field are Zr2M(M=Co,Ni,Fe)alloys,metallic Pd,depleted U,and ZrCo alloy.Specifically,Zr2M(M=Co,Ni,Fe)alloys are considered promising tritium-getter materials,and metallic Pd is utilized to separate and purify hydrogen isotopes.Furthermore,depleted U and ZrCo alloy are well suited for storing and delivering hydrogen isotopes.Notably,all the aforementioned HSAs need to modulate their hydrogen storage properties for complex operating conditions.In this review,we present a comprehensive overview of the reported modification methods applied to the above alloys.Alloying is an effective amelioration method that mainly modulates the properties of HSAs by altering their local geometrical/electronic structures.Besides,microstructural modifications such as nano-sizing and nanopores have been used to increase the specific surface area and active sites of metallic Pd and ZrCo alloys for enhancing de-/hydrogenation kinetics.The combination of metallic Pd with support materials can significantly reduce the cost and enhance the pulverization resistance.Moreover,the poisoning resistance of ZrCo alloy is improved by constructing active surfaces with selective permeability.Overall,the review is constructive for better understanding the properties and mechanisms of hydrogen isotope storage alloys and provides effective guidance for future modification research.展开更多
In this paper,numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin(DG)methods(Dong et al...In this paper,numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin(DG)methods(Dong et al.in J Sci Comput 66:321–345,2016;Dong and Wang in J Comput Appl Math 380:1–11,2020)for a one-dimensional stationary Schrödinger equation.Previous work showed that penalty parameters were required to be positive in error analysis,but the methods with zero penalty parameters worked fine in numerical simulations on coarse meshes.In this work,by performing extensive numerical experiments,we discover that zero penalty parameters lead to resonance errors in the multiscale DG methods,and taking positive penalty parameters can effectively reduce resonance errors and make the matrix in the global linear system have better condition numbers.展开更多
This paper investigates superconvergence properties of the direct discontinuous Galerkin(DDG)method with interface corrections and the symmetric DDG method for diffusion equations.We apply the Fourier analysis techniq...This paper investigates superconvergence properties of the direct discontinuous Galerkin(DDG)method with interface corrections and the symmetric DDG method for diffusion equations.We apply the Fourier analysis technique to symbolically compute eigenvalues and eigenvectors of the amplification matrices for both DDG methods with different coefficient settings in the numerical fluxes.Based on the eigen-structure analysis,we carry out error estimates of the DDG solutions,which can be decomposed into three parts:(i)dissipation errors of the physically relevant eigenvalue,which grow linearly with the time and are of order 2k for P^(k)(k=2,3)approximations;(ii)projection error from a special projection of the exact solution,which is decreasing over the time and is related to the eigenvector corresponding to the physically relevant eigenvalue;(iii)dissipative errors of non-physically relevant eigenvalues,which decay exponentially with respect to the spatial mesh sizeΔx.We observe that the errors are sensitive to the choice of the numerical flux coefficient for even degree P^(2)approximations,but are not for odd degree P^(3)approximations.Numerical experiments are provided to verify the theoretical results.展开更多
The isotope composition in precipitation has been widely considered as a tracer of monsoon activity.Compared with the coastal region,the monsoon margin usually has limited precipitation with large fluctuation and is u...The isotope composition in precipitation has been widely considered as a tracer of monsoon activity.Compared with the coastal region,the monsoon margin usually has limited precipitation with large fluctuation and is usually sensitive to climate change.The water resource management in the monsoon margin should be better planned by understanding the composition of precipitation isotope and its influencing factors.In this study,the precipitation samples were collected at five sampling sites(Baiyin City,Kongtong District,Maqu County,Wudu District,and Yinchuan City)of the monsoon margin in the northwest of China in 2022 to analyze the characteristics of stable hydrogen(δD)and oxygen(δ18O)isotopes.We analyzed the impact of meteorological factors(temperature,precipitation,and relative humidity)on the composition of precipitation isotope at daily level by regression analysis,utilized the Hybrid Single-Particle Lagrangian Integrated Trajectory(HYSPLIT)-based backward trajectory model to simulate the air mass trajectory of precipitation events,and adopted the potential source contribution function(PSCF)and concentration weighted trajectory(CWT)to analyze the water vapor sources.The results showed that compared with the global meteoric water line(GMWL),the slope of the local meteoric water line(LMWL;δD=7.34δ^(18)O-1.16)was lower,indicating the existence of strong regional evaporation in the study area.Temperature significantly contributed toδ18O value,while relative humidity had a significant negative effect onδ18O value.Through the backward trajectory analysis,we found eight primary locations that were responsible for the water vapor sources of precipitation in the study area,of which moisture from the Indian Ocean to South China Sea(ITSC)and the western continental(CW)had the greatest influence on precipitation in the study area.The hydrogen and oxygen isotopes in precipitation are significantly influenced by the sources and transportation paths of air mass.In addition,the results of PSCF and CWT analysis showed that the water vapor source areas were primarily distributed in the south and northwest direction of the study area.展开更多
基金The Guangdong Basic and Applied Basic Research Foundation(2022A1515010730)National Natural Science Foundation of China(32001647)+2 种基金National Natural Science Foundation of China(31972022)Financial and moral assistance supported by the Guangdong Basic and Applied Basic Research Foundation(2019A1515011996)111 Project(B17018)。
文摘In this study,the structural characters,antioxidant activities and bile acid-binding ability of sea buckthorn polysaccharides(HRPs)obtained by the commonly used hot water(HRP-W),pressurized hot water(HRP-H),ultrasonic(HRP-U),acid(HRP-C)and alkali(HRP-A)assisted extraction methods were investigated.The results demonstrated that extraction methods had significant effects on extraction yield,monosaccharide composition,molecular weight,particle size,triple-helical structure,and surface morphology of HRPs except for the major linkage bands.Thermogravimetric analysis showed that HRP-U with filamentous reticular microstructure exhibited better thermal stability.The HRP-A with the lowest molecular weight and highest arabinose content possessed the best antioxidant activities.Moreover,the rheological analysis indicated that HRPs with higher galacturonic acid content and molecular weight showed higher viscosity and stronger crosslinking network(HRP-C,HRP-W and HRP-U),which exhibited stronger bile acid binding capacity.The present findings provide scientific evidence in the preparation technology of sea buckthorn polysaccharides with good antioxidant and bile acid binding capacity which are related to the structure affected by the extraction methods.
文摘The nomenclature for compounds that are modified with isotopes is growing every day. Compounds can be modified with isotopes either individually, in a functional group or groups, or completely with all atomic centers of the element. This diversity of isotope-modified compounds increases the range of researches that can be studied using them. Compounds modified with isotopes of carbon-13 or nitrogen-15 can be converted into carbon monoxide, carbon dioxide and molecular nitrogen. Currently, only the average value of carbon-13 or nitrogen-15 isotopes can be determined. However, by directly determining the atomic share of these isotopes in organic compounds modified with isotopes, information about the isotopic centers of the element can be obtained. The atomic fraction of an element is defined as a single carbon or nitrogen isotope-modified center or centers, or all centers that are isotope-modified with that element at the same time. Carbon-13 or nitrogen-15 isotopes’ atomic fraction can be determined molecularly or with fragment ions of different elemental content, or both. This makes the method self-verifying, increasing the accuracy and reliability of the results obtained. Amino acids, such as asparagine, aspartic acid, methionine, and threonine, are essential for the human body. This proposed method of isotopic analysis will increase the possibilities for scientific research using these compounds.
基金supported by the Japan Society for the Promotion of Science,KAKENHI Grant No.23H00475.
文摘The inverse and direct piezoelectric and circuit coupling are widely observed in advanced electro-mechanical systems such as piezoelectric energy harvesters.Existing strongly coupled analysis methods based on direct numerical modeling for this phenomenon can be classified into partitioned or monolithic formulations.Each formulation has its advantages and disadvantages,and the choice depends on the characteristics of each coupled problem.This study proposes a new option:a coupled analysis strategy that combines the best features of the existing formulations,namely,the hybrid partitioned-monolithic method.The analysis of inverse piezoelectricity and the monolithic analysis of direct piezoelectric and circuit interaction are strongly coupled using a partitioned iterative hierarchical algorithm.In a typical benchmark problem of a piezoelectric energy harvester,this research compares the results from the proposed method to those from the conventional strongly coupled partitioned iterative method,discussing the accuracy,stability,and computational cost.The proposed hybrid concept is effective for coupled multi-physics problems,including various coupling conditions.
基金This study was supported by the National Natural Science Foundation of China(Grant Number:42007407,42022059)the Sino-German mobility program(M-0393)+1 种基金the Key Research Program of the Institute of Geology and Geophysics(CAS Grant IGGCAS-201905)the CAS Youth Interdisciplinary Team(JCTD-2021-05).
文摘Analyses of stable isotopes(C,O,H)in tree rings are increasingly important cross-disciplinary programs.The rapid development in this field documented in an increasing number of publications requires a comprehensive review.This study includes a bibliometric analysis-based review to better understand research trends in tree ring stable isotope research.Overall,1475 publications were selected from the Web of Science Core Collection for 1974-2023.The findings are that:(1)numbers of annual publications and citations increased since 1974.From 1974 to 1980,there were around two relevant publications per year.However,from 2020 to 2022,this rose sharply to 109 publications per year.Likewise,average article citations were less than four per year before 1990,but were around four per article per year after 2000;(2)the major subjects using tree ring stable isotopes include forestry,geosciences,and environmental sciences,contributing to 42.5%of the total during 1974-2023;(3)the top three most productive institutions are the Chinese Academy of Sciences(423),the Swiss Federal Institute for Forest,Snow and Landscape Research(227),and the University of Arizona(204).These achievements result from strong collaborations;(4)review papers,for example,(Dawson et al.,Annu Rev Ecol Syst 33:507-559,2002)and(McCarroll and Loader,Quat Sci Rev 23:771-801,2004),are among the most cited,with more than 1000 citations;(5)tree ring stable isotope studies mainly focus on climatology and ecology,with atmospheric CO_(2) one of the most popular topics.Since 2010,precipitation and drought have received increasing attention.Based on this analysis,the research stages,key findings,debated issues,limitations and direc-tions for future research are summarized.This study serves as an important attempt to understand the progress on the use of stable isotopes in tree rings,providing scientific guid-ance for young researchers in this field.
基金funded by the National Natural Science Foundation of China (2019M653840XB)the National Natural Science Foundation of China (41972043 and 42062006)。
文摘Geochemistry, zircon U–Pb geochronology, and Hf isotope data for the Early Paleozoic granites in the Baoshan Block reveal the Early Paleozoic tectonic evolution of the Proto-Tethys. The samples are high-K, calcalkaline, strongly peraluminous rocks with A/CNK values of 1.37–1.46, are enriched in SiO2, K2O, and Rb, and are depleted in Nb, P, Ti, Eu, and heavy rare earth elements,which indicates the crystallization fractionation of the granitic magma. Zircon U–Pb dating indicates that they formed in ca. 480 Ma. The Nansa granites have εHf(t) values ranging from-16.04 to 4.36 with corresponding TC DMages of 2.10–0.81 Ga, which suggests the magmas derived from the partial melting of ancient metasedimentary with minor involvement of mantle-derived components. A synthesis of data for the Early Paleozoic igneous rocks in the Baoshan block and adjacent(Tengchong,Qiangtang, Sibumasu, Himalaya, etc.) blocks indicates that these blocks were all aligned along the proto-Tethyan margin of East Gondwana in the Early Paleozoic. The Early Paleozoic S-type granites from Nansa were generated in a high-temperature and low-pressure(HTLP) extensional tectonic setting, which resulted from Andean-type orogeny instead of the final assembly of Gondwana or crustal extension in a non-arc environment. In certain places, an expanding environment may exist in opposition to the tectonic backdrop of the lithosphere’s thickening and shortening, leading the crust to melt and decompress,mantle-derived materials to mix, and a small quantity of peraluminous granite to emerge.
基金supported by ONR UMass Dartmouth Marine and UnderSea Technology(MUST)grant N00014-20-1-2849 under the project S31320000049160by DOE grant DE-SC0023164 sub-award RC114586-UMD+2 种基金by AFOSR grants FA9550-18-1-0383 and FA9550-23-1-0037supported by Michigan State University,by AFOSR grants FA9550-19-1-0281 and FA9550-18-1-0383by DOE grant DE-SC0023164.
文摘Additive Runge-Kutta methods designed for preserving highly accurate solutions in mixed-precision computation were previously proposed and analyzed.These specially designed methods use reduced precision for the implicit computations and full precision for the explicit computations.In this work,we analyze the stability properties of these methods and their sensitivity to the low-precision rounding errors,and demonstrate their performance in terms of accuracy and efficiency.We develop codes in FORTRAN and Julia to solve nonlinear systems of ODEs and PDEs using the mixed-precision additive Runge-Kutta(MP-ARK)methods.The convergence,accuracy,and runtime of these methods are explored.We show that for a given level of accuracy,suitably chosen MP-ARK methods may provide significant reductions in runtime.
文摘Condensed and hydrolysable tannins are non-toxic natural polyphenols that are a commercial commodity industrialized for tanning hides to obtain leather and for a growing number of other industrial applications mainly to substitute petroleum-based products.They are a definite class of sustainable materials of the forestry industry.They have been in operation for hundreds of years to manufacture leather and now for a growing number of applications in a variety of other industries,such as wood adhesives,metal coating,pharmaceutical/medical applications and several others.This review presents the main sources,either already or potentially commercial of this forestry by-materials,their industrial and laboratory extraction systems,their systems of analysis with their advantages and drawbacks,be these methods so simple to even appear primitive but nonetheless of proven effectiveness,or very modern and instrumental.It constitutes a basic but essential summary of what is necessary to know of these sustainable materials.In doing so,the review highlights some of the main challenges that remain to be addressed to deliver the quality and economics of tannin supply necessary to fulfill the industrial production requirements for some materials-based uses.
基金partially supported by U.S. National Science Foundation grants EAR-2221907partly sponsored by agencies of the United States Government。
文摘It is a challenge to make thorough but efficient experimental designs for the coupled mineral dissolution and precipitation studies in a multi-mineral system, because it is difficult to speculate the best experimental duration, optimal sampling schedule, effects of different experimental conditions, and how to maximize the experimental outputs prior to the actual experiments. Geochemical modeling is an efficient and effective tool to assist the experimental design by virtually running all scenarios of interest for the studied system and predicting the experimental outcomes. Here we demonstrated an example of geochemical modeling assisted experimental design of coupled labradorite dissolution and calcite and clayey mineral precipitation using multiple isotope tracers. In this study, labradorite(plagioclase) was chosen as the reactant because it is both a major component and one of the most reactive minerals in basalt. Following our isotope doping studies of single minerals in the last ten years, initial solutions in the simulations were doped withmultiple isotopes(e.g., Ca and Si). Geochemical modeling results show that the use of isotope tracers gives us orders of magnitude more sensitivity than the conventional method based on concentrations and allows us to decouple dissolution and precipitation reactions at near-equilibrium condition. The simulations suggest that the precise unidirectional dissolution rates can inform us which rate laws plagioclase dissolution has followed. Calcite precipitation occurred at near-equilibrium and the multiple isotope tracer experiments would provide near-equilibrium precipitation rates, which was a challenge for the conventional concentration-based experiments. In addition, whether the precipitation of clayey phases is the rate-limiting step in some multi-mineral systems will be revealed. Overall, the modeling results of multimineral reaction kinetics will improve the understanding of the coupled dissolution–precipitation in the multi-mineral systems and the quality of geochemical modeling prediction of CO_(2) removal and storage efficacy in the basalt systems.
文摘In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-dependent problems.We use the convex splitting method,the variant energy quadratization method,and the scalar auxiliary variable method coupled with the LDG method to construct first-order temporal accurate schemes based on the gradient flow structure of the models.These semi-implicit schemes are decoupled,energy stable,and can be extended to high accuracy schemes using the semi-implicit spectral deferred correction method.Many bound preserving DG discretizations are only worked on explicit time integration methods and are difficult to get high-order accuracy.To overcome these difficulties,we use the Lagrange multipliers to enforce the implicit or semi-implicit LDG schemes to satisfy the bound constraints at each time step.This bound preserving limiter results in the Karush-Kuhn-Tucker condition,which can be solved by an efficient active set semi-smooth Newton method.Various numerical experiments illustrate the high-order accuracy and the effect of bound preserving.
基金the National Natural Science Foundation of China for financial support to this work under Grant NSFC No.12072064.
文摘The collocation method is a widely used numerical method for science and engineering problems governed by partial differential equations.This paper provides a comprehensive review of collocation methods and their applications,focused on elasticity,heat conduction,electromagnetic field analysis,and fluid dynamics.The merits of the collocation method can be attributed to the need for element mesh,simple implementation,high computational efficiency,and ease in handling irregular domain problems since the collocation method is a type of node-based numerical method.Beginning with the fundamental principles of the collocation method,the discretization process in the continuous domain is elucidated,and how the collocation method approximation solutions for solving differential equations are explained.Delving into the historical development of the collocation methods,their earliest applications and key milestones are traced,thereby demonstrating their evolution within the realm of numerical computation.The mathematical foundations of collocation methods,encompassing the selection of interpolation functions,definition of weighting functions,and derivation of integration rules,are examined in detail,emphasizing their significance in comprehending the method’s effectiveness and stability.At last,the practical application of the collocation methods in engineering contexts is emphasized,including heat conduction simulations,electromagnetic coupled field analysis,and fluid dynamics simulations.These specific case studies can underscore collocation method’s broad applicability and effectiveness in addressing complex engineering challenges.In conclusion,this paper puts forward the future development trend of the collocation method through rigorous analysis and discussion,thereby facilitating further advancements in research and practical applications within these fields.
基金This paper is supported by Chinese NSF project(42130114)the strategic priority research program(B)of CAS(XDB41000000)the pre-research Project on Civil Aerospace Technologies No.D020202 funded by Chinese National Space Administration(CNSA)。
文摘It is well-known that the equilibrium isotope fractionation of mercury(Hg)includes classical massdependent fractionations(MDFs)and nuclear volume effect(NVE)induced mass-independent fractionations(MIFs).However,the effect of the NVE on these kinetic processes is not known.The total fractionations(MDFs+NVEinduced MIFs)of several representative Hg-incorporated substances were selected and calculated with ab initio calculations in this work for both equilibrium and kinetic processes.NVE-induced MIFs were calculated with scaled contact electron densities at the nucleus through systematic evaluations of their accuracy and errors using the Gaussian09 and DIRAC19 packages(named the electron density scaling method).Additionally,the NVE-induced kinetic isotope effect(KIE)of Hg isotopes are also calculated with this method for several representative Hg oxidation reactions by chlorine species.Total KIEs for 202 Hg/^(198)Hg ranging from−2.27‰to 0.96‰are obtained.Three anomalous^(202)Hg-enriched KIEs(δ^(202)Hg/^(198)Hg=0.83‰,0.94‰,and 0.96‰,)caused by the NVE are observed,which are quite different from the classical view(i.e.,light isotopes react faster than the heavy ones).The electron density scaling method we developed in this study can provide an easier way to calculate the NVE-induced KIEs for heavy isotopes and serve to better understand the fractionation mechanisms of mercury isotope systems.
基金supported by the National Key Research and Development Program of China(2016YFD0600201)the National Nonprofit Institute Research Grant of CAF(CAFYBB2017ZB003)+1 种基金the National Natural Science Foundation of China(3187071631670720)。
文摘Understanding the relationship between forest management and water use efficiency(WUE)is important for evaluating forest adaptability to climate change.However,the effects of thinning and understory removal on WUE and its key controlling processes are not well understood,which limits our comprehension of the physiological mechanisms of various management practices.In this study,four forest management measures(no thinning:NT;understory removal:UR;light thinning:LT;and heavy thinning:HT)were carried out in Pinus massoniana plantations in a subtropical region of China.Photosynthetic capacity and needle stable carbon isotope composition(δ^(13)C)were measured to assess instantaneous water use efficiency(WUE_(inst))and long-term water use efficiency(WUE_(i)).Multiple regression models and structural equation modelling(SEM)identified the effects of soil properties and physiological performances on WUE_(inst)and WUE_(i).The results show that WUE_(inst)values among the four treatments were insignificant.However,compared with the NT stand(35.8μmol·mol^(-1)),WUE_(i)values significantly increased to 41.7μmol·mol^(-1)in the UR,50.1μmol·mol^(-1)in the LT and 46.6μmol·mol^(-1)in HT treatments,largely explained by photosynthetic capacity and soil water content.Understory removal did not change physiological performance(needle water potential and photosynthetic capacity).Thinning increased the net photosynthetic rate(A_n)but not stomatal conductance(g_s)or predawn needle water potential(ψ_(pd)),implying that the improvement in water use efficiency for thinned stands was largely driven by radiation interception than by soil water availability.In general,thinning may be an appropriate management measure to promote P.massoniana WUE to cope with seasonal droughts under future extreme climates.
文摘Hydrogen is the new age alternative energy source to combat energy demand and climate change.Storage of hydrogen is vital for a nation’s growth.Works of literature provide different methods for storing the produced hydrogen,and the rational selection of a viable method is crucial for promoting sustainability and green practices.Typically,hydrogen storage is associated with diverse sustainable and circular economy(SCE)criteria.As a result,the authors consider the situation a multi-criteria decision-making(MCDM)problem.Studies infer that previous models for hydrogen storage method(HSM)selection(i)do not consider preferences in the natural language form;(ii)weights of experts are not methodically determined;(iii)hesitation of experts during criteria weight assessment is not effectively explored;and(iv)three-stage solution of a suitable selection of HSM is unexplored.Driven by these gaps,in this paper,authors put forward a new integrated framework,which considers double hierarchy linguistic information for rating,criteria importance through inter-criteria correlation(CRITIC)for expert weight calculation,evidence-based Bayesian method for criteria weight estimation,and combined compromise solution(CoCoSo)for ranking HSMs.The applicability of the developed framework is testified by using a case example of HSM selection in India.Sensitivity and comparative analysis reveal the merits and limitations of the developed framework.
基金Support for this study was received from Orient Resources Ltd.in Canada,Wuhan Institute of Technology,China,and College of Earth Sciences,Jilin University,China.
文摘By studying the light isotopic compositions of carbon,oxygen,and hydrogen,combined with previous research results on the ore-forming source of the deposit,the authors try to uncover its metallogenic origin.The δ^(18)O and δ^(13)C isotope signatures of dolomite samples vary between 10.2 and 13.0‰,and between−7.2 and−5.2‰,respectively,implying that the carbon derives from the upper mantle.δD and δ^(18) O of quartz,biotite,and muscovite from diff erent ore veins of the deposit vary between−82 and−59‰,and between 11.6 and 12.4‰,respectively,implying that the metallogenic solutions are mainly magmatic.According to the relevant research results of many isotope geologists,the fractionation degree of hydrogen isotopes increases as the depth to the Earth’s core increases,and the more diff erentiated the hydrogen isotopes are,the lower their values will be.In other words,mantle-derived solutions can have extremely low hydrogen isotope values.This means that the δD‰ value−134 of the pyrrhotite sample numbered SD-34 in this article may indicate mantle-derived oreforming fl uid of the deposit.The formation of the Dashuigou tellurium deposit occurred between 91.71 and 80.19 Ma.
基金funded by the National Natural Science Foundation of China(41872232)the Beijing Geological Survey Project(PXM 2016-158203-000008,PXM 2018-158203-000014)the Beijing Innovation Studio(Urban Geology,Active Structure,and Monitoring).
文摘The Nianzi granite unit,which includes the Nianzi,Xiaolianghou and Xiawopu granitic intrusions,is a significant component of the northern part of the North China Craton(NCC)and is situated in the Yanshan fold and thrust belt(YFTB).However,there is still debate regarding the tectonic evolutionary history of the YFTB during the late Permian to Triassic period,specifically regarding the timing of subduction and collision between the NCC and the Paleo-Asian Ocean.The Nianzi granite unit exhibits unique petrological,geochronological and geochemical signatures that shed light on the tectonic evolutionary history of the YFTB.This study presents detailed petrology,whole-rock geochemistry,together with Sr-Nd isotopic,zircon U-Pb dating and Lu-Hf isotopic data of the granites within the Nianzi granite unit.Our findings demonstrate that the granites primarily consist of subhedral K-feldspar,plagioclase,quartz,minor biotite and hornblende,with accessory titanite,apatite,magnetite and zircon.Zircon U-Pb dating indicates that the Xiaolianghou granite was emplaced at 247.5±0.62 Ma.Additionally,the adakitic characteristics of the Nianzi,Xiawopu and Xiaolianghou granitic intrusions,such as high Sr and Ba contents and high ratios of Sr/Y and(La/Yb)N,combined with negative Sr-Nd and Lu-Hf isotopes(87Sr/86Sr)i=0.705681–0.7057433,εNd(t)=−21.98 to−20.97,zirconεHf(t)=−20.26 to−9.92,as well as the I-type granite features of high SiO_(2),Na_(2)O and K_(2)O/Na_(2)O ratios,enriched Rb,K,Sr and Ba,along with depleted Th,U,Nb,Ta,P and Ti,suggest that the Nianzi granitic unit was mainly derived from the partial melting of a thickened lower crust containing hydrous,calc-alkaline to high-K calc-alkaline,mafic to intermediate metamorphic rocks.In light of these parameters,we further integrate our data with previous studies and conclude that the Nianzi granitic unit was generated in a post-collisional extensional environment during the Early Triassic.
基金supported by the NSF under Grant DMS-1818467Simons Foundation under Grant 961585.
文摘In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal energy are positive,and the mass fraction of each species is between 0 and 1.Second,due to the rapid reaction rate,the system may contain stiff sources,and the strong-stability-preserving explicit Runge-Kutta method may result in limited time-step sizes.To obtain physically relevant numerical approximations,we apply the bound-preserving technique to the DG methods.Though traditional positivity-preserving techniques can successfully yield positive density,internal energy,and mass fractions,they may not enforce the upper bound 1 of the mass fractions.To solve this problem,we need to(i)make sure the numerical fluxes in the equations of the mass fractions are consistent with that in the equation of the density;(ii)choose conservative time integrations,such that the summation of the mass fractions is preserved.With the above two conditions,the positive mass fractions have summation 1,and then,they are all between 0 and 1.For time discretization,we apply the modified Runge-Kutta/multi-step Patankar methods,which are explicit for the flux while implicit for the source.Such methods can handle stiff sources with relatively large time steps,preserve the positivity of the target variables,and keep the summation of the mass fractions to be 1.Finally,it is not straightforward to combine the bound-preserving DG methods and the Patankar time integrations.The positivity-preserving technique for DG methods requires positive numerical approximations at the cell interfaces,while Patankar methods can keep the positivity of the pre-selected point values of the target variables.To match the degree of freedom,we use polynomials on rectangular meshes for problems in two space dimensions.To evolve in time,we first read the polynomials at the Gaussian points.Then,suitable slope limiters can be applied to enforce the positivity of the solutions at those points,which can be preserved by the Patankar methods,leading to positive updated numerical cell averages.In addition,we use another slope limiter to get positive solutions used for the bound-preserving technique for the flux.Numerical examples are given to demonstrate the good performance of the proposed schemes.
基金supported by the National Key Research and Development Program of China(2022YFE03170002)the National Natural Science Foundation of China(52071286 and U2030208)the Scientific Research Fund of Zhejiang Provincial Education Department(Y202353551).
文摘Hydrogen storage alloys(HSAs)are attracting widespread interest in the nuclear industry because of the generation of stable metal hydrides after tritium absorption,thus effectively preventing the leakage of radioactive tritium.Commonly used HSAs in the hydrogen isotopes field are Zr2M(M=Co,Ni,Fe)alloys,metallic Pd,depleted U,and ZrCo alloy.Specifically,Zr2M(M=Co,Ni,Fe)alloys are considered promising tritium-getter materials,and metallic Pd is utilized to separate and purify hydrogen isotopes.Furthermore,depleted U and ZrCo alloy are well suited for storing and delivering hydrogen isotopes.Notably,all the aforementioned HSAs need to modulate their hydrogen storage properties for complex operating conditions.In this review,we present a comprehensive overview of the reported modification methods applied to the above alloys.Alloying is an effective amelioration method that mainly modulates the properties of HSAs by altering their local geometrical/electronic structures.Besides,microstructural modifications such as nano-sizing and nanopores have been used to increase the specific surface area and active sites of metallic Pd and ZrCo alloys for enhancing de-/hydrogenation kinetics.The combination of metallic Pd with support materials can significantly reduce the cost and enhance the pulverization resistance.Moreover,the poisoning resistance of ZrCo alloy is improved by constructing active surfaces with selective permeability.Overall,the review is constructive for better understanding the properties and mechanisms of hydrogen isotope storage alloys and provides effective guidance for future modification research.
基金supported by the National Science Foundation grant DMS-1818998.
文摘In this paper,numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin(DG)methods(Dong et al.in J Sci Comput 66:321–345,2016;Dong and Wang in J Comput Appl Math 380:1–11,2020)for a one-dimensional stationary Schrödinger equation.Previous work showed that penalty parameters were required to be positive in error analysis,but the methods with zero penalty parameters worked fine in numerical simulations on coarse meshes.In this work,by performing extensive numerical experiments,we discover that zero penalty parameters lead to resonance errors in the multiscale DG methods,and taking positive penalty parameters can effectively reduce resonance errors and make the matrix in the global linear system have better condition numbers.
基金supported by the National Natural Science Foundation of China(Grant Nos.11871428 and 12071214)the Natural Science Foundation for Colleges and Universities of Jiangsu Province of China(Grant No.20KJB110011)+1 种基金supported by the National Science Foundation(Grant No.DMS-1620335)and the Simons Foundation(Grant No.637716)supported by the National Natural Science Foundation of China(Grant Nos.11871428 and 12272347).
文摘This paper investigates superconvergence properties of the direct discontinuous Galerkin(DDG)method with interface corrections and the symmetric DDG method for diffusion equations.We apply the Fourier analysis technique to symbolically compute eigenvalues and eigenvectors of the amplification matrices for both DDG methods with different coefficient settings in the numerical fluxes.Based on the eigen-structure analysis,we carry out error estimates of the DDG solutions,which can be decomposed into three parts:(i)dissipation errors of the physically relevant eigenvalue,which grow linearly with the time and are of order 2k for P^(k)(k=2,3)approximations;(ii)projection error from a special projection of the exact solution,which is decreasing over the time and is related to the eigenvector corresponding to the physically relevant eigenvalue;(iii)dissipative errors of non-physically relevant eigenvalues,which decay exponentially with respect to the spatial mesh sizeΔx.We observe that the errors are sensitive to the choice of the numerical flux coefficient for even degree P^(2)approximations,but are not for odd degree P^(3)approximations.Numerical experiments are provided to verify the theoretical results.
基金supported by the National Natural Science Foundation of China(42161007)the Scientific Research Program for Higher Education Institutions of Gansu Province(2021B-081)the Natural Science Foundation of Gansu Province(22JR5RA074).
文摘The isotope composition in precipitation has been widely considered as a tracer of monsoon activity.Compared with the coastal region,the monsoon margin usually has limited precipitation with large fluctuation and is usually sensitive to climate change.The water resource management in the monsoon margin should be better planned by understanding the composition of precipitation isotope and its influencing factors.In this study,the precipitation samples were collected at five sampling sites(Baiyin City,Kongtong District,Maqu County,Wudu District,and Yinchuan City)of the monsoon margin in the northwest of China in 2022 to analyze the characteristics of stable hydrogen(δD)and oxygen(δ18O)isotopes.We analyzed the impact of meteorological factors(temperature,precipitation,and relative humidity)on the composition of precipitation isotope at daily level by regression analysis,utilized the Hybrid Single-Particle Lagrangian Integrated Trajectory(HYSPLIT)-based backward trajectory model to simulate the air mass trajectory of precipitation events,and adopted the potential source contribution function(PSCF)and concentration weighted trajectory(CWT)to analyze the water vapor sources.The results showed that compared with the global meteoric water line(GMWL),the slope of the local meteoric water line(LMWL;δD=7.34δ^(18)O-1.16)was lower,indicating the existence of strong regional evaporation in the study area.Temperature significantly contributed toδ18O value,while relative humidity had a significant negative effect onδ18O value.Through the backward trajectory analysis,we found eight primary locations that were responsible for the water vapor sources of precipitation in the study area,of which moisture from the Indian Ocean to South China Sea(ITSC)and the western continental(CW)had the greatest influence on precipitation in the study area.The hydrogen and oxygen isotopes in precipitation are significantly influenced by the sources and transportation paths of air mass.In addition,the results of PSCF and CWT analysis showed that the water vapor source areas were primarily distributed in the south and northwest direction of the study area.