Reactive transport modeling(RTM)is an emerging method used to address geological issues in diagenesis research.However,the extrapolation of RTM results to practical reservoir prediction is not sufficiently understood....Reactive transport modeling(RTM)is an emerging method used to address geological issues in diagenesis research.However,the extrapolation of RTM results to practical reservoir prediction is not sufficiently understood.This paper presents a case study of the Eocene Qaidam Basin that combines RTM results with petrological and mineralogical evidence.The results show that the Eocene Xiaganchaigou Formation is characterized by mixed siliciclastic-carbonate-evaporite sedimentation in a semiclosed saline lacustrine environment.Periodic evaporation and salinization processes during the syngeneticpenecontemporaneous stage gave rise to the replacive genesis of dolomites and the cyclic enrichment of dolomite in the middle-upper parts of the meter-scale depositional sequences.The successive change in mineral paragenesis from terrigenous clastics to carbonates to evaporites was reconstructed using RTM simulations.Parametric uncertainty analyses further suggest that the evaporation intensity(brine salinity)and particle size of sediments(reactive surface area)were important rate-determining factors in the dolomitization,as shown by the relatively higher reaction rates under conditions of higher brine salinity and fine-grained sediments.Combining the simulation results with measured mineralogical and reservoir physical property data indicates that the preservation of original intergranular pores and the generation of porosity via replacive dolomitization were the major formation mechanisms of the distinctive lacustrine dolomite reservoirs(widespread submicron intercrystalline micropores)in the Eocene Qaidam Basin.The results confirm that RTM can be effectively used in geological studies,can provide a better general understanding of the dolomitizing fluid-rock interactions,and can shed light on the spatiotemporal evolution of mineralogy and porosity during dolomitization and the formation of lacustrine dolomite reservoirs.展开更多
The elliptic azimuthal anisotropy coefficient(v_(2))of the identified particles at midrapidity(|η|<0.8)was investigated in p-Pb collisions at√s_(NN)=5.02 TeV using a multi-phase transport model(AMPT).The calculat...The elliptic azimuthal anisotropy coefficient(v_(2))of the identified particles at midrapidity(|η|<0.8)was investigated in p-Pb collisions at√s_(NN)=5.02 TeV using a multi-phase transport model(AMPT).The calculations of differential v_(2)based on the advanced flow extraction method of light flavor hadrons(pions,kaons,protons,andΛ)in small collision systems were extended to a wider transverse momentum(p_(T))range of up to 8 GeV/c for the first time.The string-melting version of the AMPT model provides a good description of the measured p_(T)-differential v_(2)of the mesons but exhibits a slight deviation from the baryon v_(2).In addition,we observed the features of mass ordering at low p_(T)and the approximate number-of-constituentquark(NCQ)scaling at intermediate p_(T).Moreover,we demonstrate that hadronic rescattering does not have a significant impact on v_(2)in p-Pb collisions for different centrality selections,whereas partonic scattering dominates in generating the elliptic anisotropy of the final particles.This study provides further insight into the origin of collective-like behavior in small collision systems and has referential value for future measurements of azimuthal anisotropy.展开更多
Transport analysis and impact evaluations are important input for decisions about infrastructure projects.The impacts on transport from fjord crossing tunnels or bridges are the foundation for the cost benefit analysi...Transport analysis and impact evaluations are important input for decisions about infrastructure projects.The impacts on transport from fjord crossing tunnels or bridges are the foundation for the cost benefit analysis,and also the basis for estimating the income from toll collection.Based on experiences from concept evaluations of several fixed link projects on E39,and an ongoing overall analysis,we question the results from transport analysis made by the official tools for such analysis:the RTM(regional transport model)which estimates the demand for trips below 10 km,the NTM(national transport model),for trips of 10 km or more,and the freight transport model.Both the NTM and the freight transport model are integrated in the RTM in the net assignment stage.We will demonstrate strengths and weaknesses in the transport models by showing contra intuitive or questionable results using the model as it is.The following questions arose as the initial results from the transport model were presented:Are the transport models able to capture immediate as well as long-term impacts?How would different assumptions about the monetary costs on these projects affect the forecasted demand and the cost benefit analysis?Are there other and wider ranges of impacts,if the analysis covers the total coastal highway as a whole,compared to evaluating impacts of each fixed link project individually?Do we have enough data to include transport effects of wider impacts of the fixed link projects?We had to deal with these questions in the concept evaluations carried out for the various fixed links project and in the current overall evaluation.We would like to suggest improvements in the analysis tools and emphasize requirements for knowledge about impacts of fixed links projects.展开更多
A multi-phase transport(AMPT)model was constructed as a self-contained kinetic theory-based description of relativistic nuclear collisions as it contains four main components:the fluctuating initial condition,a parton...A multi-phase transport(AMPT)model was constructed as a self-contained kinetic theory-based description of relativistic nuclear collisions as it contains four main components:the fluctuating initial condition,a parton cascade,hadronization,and a hadron cascade.Here,we review the main developments after the first public release of the AMPT source code in 2004 and the corre-sponding publication that described the physics details of the model at that time.We also discuss possible directions for future developments of the AMPT model to better study the properties of the dense matter created in relativistic collisions of small or large systems.展开更多
A regional chemical transport model,RAMS-CMAQ,was employed to assess the impacts of biosphere-atmosphere CO2 exchange on seasonal variations in atmospheric CO2 concentrations over East Asia.Simulated CO2 concentration...A regional chemical transport model,RAMS-CMAQ,was employed to assess the impacts of biosphere-atmosphere CO2 exchange on seasonal variations in atmospheric CO2 concentrations over East Asia.Simulated CO2 concentrations were compared with observations at 12 surface stations and the comparison showed they were generally in good agreement.Both observations and simulations suggested that surface CO2 over East Asia features a summertime trough due to biospheric absorption,while in some urban areas surface CO2 has a distinct summer peak,which could be attributed to the strong impact from anthropogenic emissions.Analysis of the model results indicated that biospheric fluxes and fossil-fuel emissions are comparably important in shaping spatial distributions of CO2 near the surface over East Asia.Biospheric flux plays an important role in the prevailing spatial pattern of CO2 enhancement and reduction on the synoptic scale due to the strong seasonality of biospheric CO2 flux.The elevation of CO2 levels by the biosphere during winter was found to be larger than 5 ppm in North China and Southeast China,and during summertime a significant depletion （≥ 7 ppm） occurred in most areas,except for the Indo-China Peninsula where positive bioflux values were found.展开更多
One of the largest known megafloods on earth resulted from a glacier dam-break,which occurred during the Late Quaternary in the Altai Mountains in Southern Siberia.Computational modeling is one of the viable approache...One of the largest known megafloods on earth resulted from a glacier dam-break,which occurred during the Late Quaternary in the Altai Mountains in Southern Siberia.Computational modeling is one of the viable approaches to enhancing the understanding of the flood events.The computational domain of this flood is over 9460 km2 and about 3.784 × 106 cells are involved as a 50 m × 50 m mesh is used,which necessitates a computationally efficient model.Here the Open MP(Open Multiprocessing) technique is adopted to parallelize the code of a coupled 2D hydrodynamic and sediment transport model.It is shown that the computational efficiency is enhanced by over 80% due to the parallelization.The floods over both fixed and mobile beds are well reproduced with specified discharge hydrographs at the dam site.Qualitatively,backwater effects during the flood are resolved at the bifurcation between the Chuja and Katun rivers.Quantitatively,the computed maximum stage and thalweg are physically consistent with the field data of the bars and deposits.The effects of sediment transport and morphological evolution on the flood are considerable.Sensitivity analyses indicate that the impact of the peak discharge is significant,whilst those of the Manningroughness,medium sediment size and shape of the inlet discharge hydrograph are marginal.展开更多
In this article, the global existence and the large time behavior of smooth solutions to the initial boundary value problem for a degenerate compressible energy transport model are established.
A global transport model is proposed in which a multimoment constrained finite volume （MCV） scheme is applied to a Yin-Yang overset grid. The MCV scheme defines 16 degrees of freedom （DOFs） within each element to ...A global transport model is proposed in which a multimoment constrained finite volume （MCV） scheme is applied to a Yin-Yang overset grid. The MCV scheme defines 16 degrees of freedom （DOFs） within each element to build a 2D cubic reconstruction polynomial. The time evolution equations for DOFs are derived from constraint conditions on moments of line-integrated averages （LIA）, point values （PV）, and values of first-order derivatives （DV）. The Yin-Yang grid eliminates polar singularities and results in a quasi-uniform mesh. A limiting projection is designed to remove nonphysical oscillations around discontinuities. Our model was tested against widely used benchmarks; the competitive results reveal that the model is accurate and promising for developing general circulation models.展开更多
A positivity-preserving conservative semi-Lagrangian transport model by multi-moment finite volume method has been developed on the cubed-sphere grid.Two kinds of moments(i.e.,point values(PV moment) at cell interface...A positivity-preserving conservative semi-Lagrangian transport model by multi-moment finite volume method has been developed on the cubed-sphere grid.Two kinds of moments(i.e.,point values(PV moment) at cell interfaces and volume integrated average(VIA moment) value) are defined within a single cell.The PV moment is updated by a conventional semi-Lagrangian method,while the VIA moment is cast by the flux form formulation to assure the exact numerical conservation.Different from the spatial approximation used in the CSL2(conservative semi-Lagrangian scheme with second order polynomial function) scheme,a monotonic rational function which can effectively remove non-physical oscillations is reconstructed within a single cell by the PV moments and VIA moment.To achieve exactly positive-definite preserving,two kinds of corrections are made on the original conservative semi-Lagrangian with rational function(CSLR)scheme.The resulting scheme is inherently conservative,non-negative,and allows a Courant number larger than one.Moreover,the spatial reconstruction can be performed within a single cell,which is very efficient and economical for practical implementation.In addition,a dimension-splitting approach coupled with multi-moment finite volume scheme is adopted on cubed-sphere geometry,which benefitsthe implementation of the 1 D CSLR solver with large Courant number.The proposed model is evaluated by several widely used benchmark tests on cubed-sphere geometry.Numerical results show that the proposed transport model can effectively remove nonphysical oscillations and preserve the numerical nonnegativity,and it has the potential to transport the tracers accurately in a real atmospheric model.展开更多
Carbon dioxide injection into deep saline aquifers results in a variety of strongly coupled physical and chemical processes. In this study, reactive transport simulations using a 2-D radial model were performed to inv...Carbon dioxide injection into deep saline aquifers results in a variety of strongly coupled physical and chemical processes. In this study, reactive transport simulations using a 2-D radial model were performed to investigate the fate of the injected CO2, the effect of CO2-water-rock interactions on mineral alteration, and the long-term CO2 sequestration mechanisms of the Liujiagou Formation sandstone at the Shenhua CCS（carbon capture and storage） pilot site of China. Carbon dioxide was injected at a constant rate of 0.1 Mt/year for 30 years, and the fluid flow and geochemical transport simulation was run for a period of 10 000 years by the TOUGHREACT code according to the underground conditions of the Liujiagou Formation. The results show that different trapping phases of CO2 vary with time. Sensitivity analyses indicate that plagioclase composition and chlorite presence are the most significant determinants of stable carbonate minerals and CO2 mineral trapping capacity. For arkosic arenite in the Liujiagou Formation, CO2 can be immobilized by precipitation of ankerite, magnesite, siderite, dawsonite, and calcite for different mineral compositions, with Ca（2＋）, Mg（2＋）, Fe（2＋） and Na＋ provided by dissolution of calcite, albite（or oligoclase） and chlorite. This study can provide useful insights into the geochemistry of CO2 storage in other arkosic arenite（feldspar rich sandstone） formations at other pilots or target sites.展开更多
A 2-dimensional atmospheric transport model is deployed in a simplified CO 2 inverse study. Calculated carbon flux distribution for the interval from 1981 to 1997 confirms the existence of a terrestrial carbon sink i...A 2-dimensional atmospheric transport model is deployed in a simplified CO 2 inverse study. Calculated carbon flux distribution for the interval from 1981 to 1997 confirms the existence of a terrestrial carbon sink in mid-high latitude area of North Hemisphere. Strong interannual variability exists in carbon flux patterns, implying a possible link with ENSO and other natural episodes such as Pinatubo volcano eruption in 1991. Mechanism of this possible link was investigated with statistic method. Correlation analysis indicated that in North Hemisphere, climatic factors such as temperature and precipitation, to some extend, could influence the carbon cycle process of land and ocean, thus cause considerable change in carbon flux distribution. In addition, correlation study also demonstrated the possible important role of Asian terrestrial ecosystems in carbon cycle.展开更多
The recent work of Nation et al., in which the Hawking radiation energy and entropy flow from a black hole is considered to be produced in a one-dimensional Landauer transport process, is extended to the case of a Rei...The recent work of Nation et al., in which the Hawking radiation energy and entropy flow from a black hole is considered to be produced in a one-dimensional Landauer transport process, is extended to the case of a Reissner- Nordstrom black hole. The energy flow contains not only the contribution of the thermal flux but also that of the particle flux. It is found that the charge can also be transported via the one-dimensional quantum tunnel. Because of the existence of the electrostatic potential, the entropy production rate is shown to be smaller than that of the Schwarzschild black hole.展开更多
Within the framework of a semiclassical Boltzmann-Uehling-Uhlenbeck （BUU） transport model, the high mo- mentum tail （HMT） effects of nucleon momentum distribution in the nucleus on the nucleon collective flows are...Within the framework of a semiclassical Boltzmann-Uehling-Uhlenbeck （BUU） transport model, the high mo- mentum tail （HMT） effects of nucleon momentum distribution in the nucleus on the nucleon collective flows are studied in semieentral Au＋Au collisions. The HMT due to the isospin-dependent short-range correlations causes a smaller value of the collective flows. We find that the HMT effects on the nucleon collective flows are remarkable at beam energy of 300 MeV/nucleon and become weak as the incident beam energy increases. The results indicate that for the collective flow studies at intermediate energies, the HMT of nucleon momentum distribution in nucleus should be taken into account in transport models.展开更多
We investigate the Hawking radiation of a Kaluza-Klein black hole by using one-dimensional(1D),nonequilibrium,Landauer transport model.The derived Hawking radiation temperature is in consistence with that obtained by ...We investigate the Hawking radiation of a Kaluza-Klein black hole by using one-dimensional(1D),nonequilibrium,Landauer transport model.The derived Hawking radiation temperature is in consistence with that obtained by using the usual anomaly method.With the Landauer transport model,we calculate the entropy flow out of the Kaluza-Klein black hole and the relevant entropy production rate.How these quantities depending on the physical parameters of the black hole is also discussed.展开更多
In the high-energy environment of coastal seas and estuaries,strong sediment resuspension/ deposition events are driven by surface waves,tides,winds and buoyancy driven currents.In recent years,A POM based three-dimen...In the high-energy environment of coastal seas and estuaries,strong sediment resuspension/ deposition events are driven by surface waves,tides,winds and buoyancy driven currents.In recent years,A POM based three-dimensional,wave-current coupled,sediment transport model has been developed by the University of New South Wales.This paper presents several examples of the model applications to study sediment dynamics in the environments where forcings such as waves,tides,and winds are equally important to affect sediment fluxes and distributions.Firstly,the sediment transport model coupled to the Yellow Sea general circulation model and a third generation wave model SWAN was implemented in the Yellow Sea to study the dynamics of the sediment transport and resuspension in the northern Jiangsu shoal-wate(rNJSW).The sediment distributions and fluxes and their inter-annual variability were studied by realistic numerical simulations.The study found that the surface waves played a dominant role over the tides to form the turbidity maxima along the muddy coast of NJSW. Secondly,the sediment transport model was used to explore the effect of suspended sediment-induced stratification in the bottom boundary laye(rBBL).The model uses a re-parameterized bottom drag coefficient Cd that incorporates a linear stability function of flux Richardson number Rf.The study has shown that the sediment induced stratification in the BBL reduces the vertical eddy viscosity and bottom shear stress in comparison with the model prediction in a neutrally stratified BBL.In response to these apparent reductions,the tidal current shear is increased and sediments are abnormally concentrated within a thin wall layer that is overlain by a thicker layer with much smaller concentration.The formation of this fluid-mud layer near the seabed has led to a significant reduction in the total sediment transport.This study contributes to the understanding of formations of tidal flats along the coasts of turbid seas and estuaries.展开更多
Acid mine drainage（AMD） is commonly treated by neutralization with alkaline substances.This treatment is supported by titration experiments that illustrate the buffering mechanisms and estimate the base neutralizati...Acid mine drainage（AMD） is commonly treated by neutralization with alkaline substances.This treatment is supported by titration experiments that illustrate the buffering mechanisms and estimate the base neutralization capacity（BNC） of the AMD.Detailed explanation of titration curves requires modeling with a hydro-chemical model.In this study the titration curves of water samples from the drainage of the As Pontes mine and the corresponding dumps have been investigated and six buffers are selerted by analyzing those curves.Titration curves have been simulated by a reactive transport model to discover the detailed buffering mechanisms.These simulations show seven regions involving different buffering mechanism.The BNC is primarily from buffers of dissolved Fe,Al and hydrogen sulfate.The BNC can be approximated by：BNC = 3（C_（Fe） ＋ C_（Al）） ＋ 0.05C_（sulfate）,where the units are mol/L.The BNC of the sample from the mine is 9.25×10^（-3） mol/L and that of the dumps sample is 1.28×10^（-2） mol/L展开更多
The elliptic flow v2, for π±, K±, p and p in Au＋Au collisions at center-of-mass energies √sNN=7.7, 11.5, 14.5 and 19.6 GeV, is analyzed using a multiphase transport model. A significant difference in the ...The elliptic flow v2, for π±, K±, p and p in Au＋Au collisions at center-of-mass energies √sNN=7.7, 11.5, 14.5 and 19.6 GeV, is analyzed using a multiphase transport model. A significant difference in the v2 values for p and p is observed, and the values of v2 splitting are larger compared with π＋ and π-, K＋ and K-. The difference increases with decreasing the center-of-mass energy. The effect of the quark coalescence mechanism in a multi-phase transport model to the value of elliptic difference △v2 between p and p- has been discussed. The simulation of Au＋Au collisions at 14.5 GeV shows that the effect of hadron cascade to △v2 is not obvious, and a larger patton-scattering cross section can lead to a larger △v2.展开更多
Identifying the unknown geometric and material information of a multi-shield object by analyzing the radiation signature measurements is always an important problem in national and global security. In order to identif...Identifying the unknown geometric and material information of a multi-shield object by analyzing the radiation signature measurements is always an important problem in national and global security. In order to identify the unknown shielding layer thicknesses of a source/shield system with gamma-ray spectra, we have developed a derivative-free inverse radiation transport model based on a differential evolution algorithm with global and local neighbourhoods(IRT-DEGL). In the present paper, the IRT-DEGL model is further extended for estimating the unknown thicknesses with random initial guesses and material mass densities of multi-shielding layers as well as their combinations. Using the detected gamma-ray spectra,the illustration of inverse studies is implemented and the main influence factors for inverse results are also analyzed.展开更多
文摘Reactive transport modeling(RTM)is an emerging method used to address geological issues in diagenesis research.However,the extrapolation of RTM results to practical reservoir prediction is not sufficiently understood.This paper presents a case study of the Eocene Qaidam Basin that combines RTM results with petrological and mineralogical evidence.The results show that the Eocene Xiaganchaigou Formation is characterized by mixed siliciclastic-carbonate-evaporite sedimentation in a semiclosed saline lacustrine environment.Periodic evaporation and salinization processes during the syngeneticpenecontemporaneous stage gave rise to the replacive genesis of dolomites and the cyclic enrichment of dolomite in the middle-upper parts of the meter-scale depositional sequences.The successive change in mineral paragenesis from terrigenous clastics to carbonates to evaporites was reconstructed using RTM simulations.Parametric uncertainty analyses further suggest that the evaporation intensity(brine salinity)and particle size of sediments(reactive surface area)were important rate-determining factors in the dolomitization,as shown by the relatively higher reaction rates under conditions of higher brine salinity and fine-grained sediments.Combining the simulation results with measured mineralogical and reservoir physical property data indicates that the preservation of original intergranular pores and the generation of porosity via replacive dolomitization were the major formation mechanisms of the distinctive lacustrine dolomite reservoirs(widespread submicron intercrystalline micropores)in the Eocene Qaidam Basin.The results confirm that RTM can be effectively used in geological studies,can provide a better general understanding of the dolomitizing fluid-rock interactions,and can shed light on the spatiotemporal evolution of mineralogy and porosity during dolomitization and the formation of lacustrine dolomite reservoirs.
基金This work was supported by the Key Laboratory of Quark and Lepton Physics(MOE)in Central China Normal University(Nos.QLPL2022P01,QLPL202106)Natural Science Foundation of Hubei Provincial Education Department(No.Q20131603)+2 种基金National key research,development program of China(No.2018YFE0104700)National Natural Science Foundation of China(No.12175085)Fundamental research funds for the Central Universities(No.CCNU220N003).
文摘The elliptic azimuthal anisotropy coefficient(v_(2))of the identified particles at midrapidity(|η|<0.8)was investigated in p-Pb collisions at√s_(NN)=5.02 TeV using a multi-phase transport model(AMPT).The calculations of differential v_(2)based on the advanced flow extraction method of light flavor hadrons(pions,kaons,protons,andΛ)in small collision systems were extended to a wider transverse momentum(p_(T))range of up to 8 GeV/c for the first time.The string-melting version of the AMPT model provides a good description of the measured p_(T)-differential v_(2)of the mesons but exhibits a slight deviation from the baryon v_(2).In addition,we observed the features of mass ordering at low p_(T)and the approximate number-of-constituentquark(NCQ)scaling at intermediate p_(T).Moreover,we demonstrate that hadronic rescattering does not have a significant impact on v_(2)in p-Pb collisions for different centrality selections,whereas partonic scattering dominates in generating the elliptic anisotropy of the final particles.This study provides further insight into the origin of collective-like behavior in small collision systems and has referential value for future measurements of azimuthal anisotropy.
文摘Transport analysis and impact evaluations are important input for decisions about infrastructure projects.The impacts on transport from fjord crossing tunnels or bridges are the foundation for the cost benefit analysis,and also the basis for estimating the income from toll collection.Based on experiences from concept evaluations of several fixed link projects on E39,and an ongoing overall analysis,we question the results from transport analysis made by the official tools for such analysis:the RTM(regional transport model)which estimates the demand for trips below 10 km,the NTM(national transport model),for trips of 10 km or more,and the freight transport model.Both the NTM and the freight transport model are integrated in the RTM in the net assignment stage.We will demonstrate strengths and weaknesses in the transport models by showing contra intuitive or questionable results using the model as it is.The following questions arose as the initial results from the transport model were presented:Are the transport models able to capture immediate as well as long-term impacts?How would different assumptions about the monetary costs on these projects affect the forecasted demand and the cost benefit analysis?Are there other and wider ranges of impacts,if the analysis covers the total coastal highway as a whole,compared to evaluating impacts of each fixed link project individually?Do we have enough data to include transport effects of wider impacts of the fixed link projects?We had to deal with these questions in the concept evaluations carried out for the various fixed links project and in the current overall evaluation.We would like to suggest improvements in the analysis tools and emphasize requirements for knowledge about impacts of fixed links projects.
基金Z.-W.L.is supported in part by the National Science Foundation under Grant No.PHY-2012947L.Z.is supported in part by the National Natural Science Foundation of China under Grant No.11905188.
文摘A multi-phase transport(AMPT)model was constructed as a self-contained kinetic theory-based description of relativistic nuclear collisions as it contains four main components:the fluctuating initial condition,a parton cascade,hadronization,and a hadron cascade.Here,we review the main developments after the first public release of the AMPT source code in 2004 and the corre-sponding publication that described the physics details of the model at that time.We also discuss possible directions for future developments of the AMPT model to better study the properties of the dense matter created in relativistic collisions of small or large systems.
基金supported by the Strategic Priority Research Program–Climate Change: Carbon Budget and Relevant Issues (Grant No. XDA05040404)the National High Technology Research and Development Program of China (Grant No. 2013AA122002)+1 种基金the National Natural Science Foundation of China (Grant No. 41130528)support staff of the stations for providing their CO2 observation records on the WDCGG website
文摘A regional chemical transport model,RAMS-CMAQ,was employed to assess the impacts of biosphere-atmosphere CO2 exchange on seasonal variations in atmospheric CO2 concentrations over East Asia.Simulated CO2 concentrations were compared with observations at 12 surface stations and the comparison showed they were generally in good agreement.Both observations and simulations suggested that surface CO2 over East Asia features a summertime trough due to biospheric absorption,while in some urban areas surface CO2 has a distinct summer peak,which could be attributed to the strong impact from anthropogenic emissions.Analysis of the model results indicated that biospheric fluxes and fossil-fuel emissions are comparably important in shaping spatial distributions of CO2 near the surface over East Asia.Biospheric flux plays an important role in the prevailing spatial pattern of CO2 enhancement and reduction on the synoptic scale due to the strong seasonality of biospheric CO2 flux.The elevation of CO2 levels by the biosphere during winter was found to be larger than 5 ppm in North China and Southeast China,and during summertime a significant depletion （≥ 7 ppm） occurred in most areas,except for the Indo-China Peninsula where positive bioflux values were found.
基金funded by Natural Science Foundation of China (Grants No. 11172217 and 11432015)National Key Basic Research and Development Program (i.e., 973 Program) of China (Grant No. 2007CB714106)
文摘One of the largest known megafloods on earth resulted from a glacier dam-break,which occurred during the Late Quaternary in the Altai Mountains in Southern Siberia.Computational modeling is one of the viable approaches to enhancing the understanding of the flood events.The computational domain of this flood is over 9460 km2 and about 3.784 × 106 cells are involved as a 50 m × 50 m mesh is used,which necessitates a computationally efficient model.Here the Open MP(Open Multiprocessing) technique is adopted to parallelize the code of a coupled 2D hydrodynamic and sediment transport model.It is shown that the computational efficiency is enhanced by over 80% due to the parallelization.The floods over both fixed and mobile beds are well reproduced with specified discharge hydrographs at the dam site.Qualitatively,backwater effects during the flood are resolved at the bifurcation between the Chuja and Katun rivers.Quantitatively,the computed maximum stage and thalweg are physically consistent with the field data of the bars and deposits.The effects of sediment transport and morphological evolution on the flood are considerable.Sensitivity analyses indicate that the impact of the peak discharge is significant,whilst those of the Manningroughness,medium sediment size and shape of the inlet discharge hydrograph are marginal.
基金Supported by the Foundation for Talents of Beijing (20081D0501500171)the Funds of Beijing University of Technology
文摘In this article, the global existence and the large time behavior of smooth solutions to the initial boundary value problem for a degenerate compressible energy transport model are established.
基金supported by National Key Technology R&D Program of China (Grant No. 2012BAC22B01)Natural Science Foundation of China (Grant Nos. 10902116, 40805045, and 41175095)Grants-in-Aid for Scientific Research, Japan Society for the Promotion of Science (Grant No. 24560187)
文摘A global transport model is proposed in which a multimoment constrained finite volume （MCV） scheme is applied to a Yin-Yang overset grid. The MCV scheme defines 16 degrees of freedom （DOFs） within each element to build a 2D cubic reconstruction polynomial. The time evolution equations for DOFs are derived from constraint conditions on moments of line-integrated averages （LIA）, point values （PV）, and values of first-order derivatives （DV）. The Yin-Yang grid eliminates polar singularities and results in a quasi-uniform mesh. A limiting projection is designed to remove nonphysical oscillations around discontinuities. Our model was tested against widely used benchmarks; the competitive results reveal that the model is accurate and promising for developing general circulation models.
基金supported by the National Key Research and Development Program of China (Grant Nos.2017YFC1501901 and 2017YFA0603901)the Beijing Natural Science Foundation (Grant No.JQ18001)。
文摘A positivity-preserving conservative semi-Lagrangian transport model by multi-moment finite volume method has been developed on the cubed-sphere grid.Two kinds of moments(i.e.,point values(PV moment) at cell interfaces and volume integrated average(VIA moment) value) are defined within a single cell.The PV moment is updated by a conventional semi-Lagrangian method,while the VIA moment is cast by the flux form formulation to assure the exact numerical conservation.Different from the spatial approximation used in the CSL2(conservative semi-Lagrangian scheme with second order polynomial function) scheme,a monotonic rational function which can effectively remove non-physical oscillations is reconstructed within a single cell by the PV moments and VIA moment.To achieve exactly positive-definite preserving,two kinds of corrections are made on the original conservative semi-Lagrangian with rational function(CSLR)scheme.The resulting scheme is inherently conservative,non-negative,and allows a Courant number larger than one.Moreover,the spatial reconstruction can be performed within a single cell,which is very efficient and economical for practical implementation.In addition,a dimension-splitting approach coupled with multi-moment finite volume scheme is adopted on cubed-sphere geometry,which benefitsthe implementation of the 1 D CSLR solver with large Courant number.The proposed model is evaluated by several widely used benchmark tests on cubed-sphere geometry.Numerical results show that the proposed transport model can effectively remove nonphysical oscillations and preserve the numerical nonnegativity,and it has the potential to transport the tracers accurately in a real atmospheric model.
基金supported by the Global Climate and Energy Project(No.2384638-43106-A)the National Natural Science Foundation of China(No.41072180)+1 种基金the Special Scientific Research Fund of Public Welfare Profession of the Ministry of Land and Resources of China(No.201211063)a bilateral project of China Australia Geological Storage of CO2 Project Phase 2(CAGS2)
文摘Carbon dioxide injection into deep saline aquifers results in a variety of strongly coupled physical and chemical processes. In this study, reactive transport simulations using a 2-D radial model were performed to investigate the fate of the injected CO2, the effect of CO2-water-rock interactions on mineral alteration, and the long-term CO2 sequestration mechanisms of the Liujiagou Formation sandstone at the Shenhua CCS（carbon capture and storage） pilot site of China. Carbon dioxide was injected at a constant rate of 0.1 Mt/year for 30 years, and the fluid flow and geochemical transport simulation was run for a period of 10 000 years by the TOUGHREACT code according to the underground conditions of the Liujiagou Formation. The results show that different trapping phases of CO2 vary with time. Sensitivity analyses indicate that plagioclase composition and chlorite presence are the most significant determinants of stable carbonate minerals and CO2 mineral trapping capacity. For arkosic arenite in the Liujiagou Formation, CO2 can be immobilized by precipitation of ankerite, magnesite, siderite, dawsonite, and calcite for different mineral compositions, with Ca（2＋）, Mg（2＋）, Fe（2＋） and Na＋ provided by dissolution of calcite, albite（or oligoclase） and chlorite. This study can provide useful insights into the geochemistry of CO2 storage in other arkosic arenite（feldspar rich sandstone） formations at other pilots or target sites.
文摘A 2-dimensional atmospheric transport model is deployed in a simplified CO 2 inverse study. Calculated carbon flux distribution for the interval from 1981 to 1997 confirms the existence of a terrestrial carbon sink in mid-high latitude area of North Hemisphere. Strong interannual variability exists in carbon flux patterns, implying a possible link with ENSO and other natural episodes such as Pinatubo volcano eruption in 1991. Mechanism of this possible link was investigated with statistic method. Correlation analysis indicated that in North Hemisphere, climatic factors such as temperature and precipitation, to some extend, could influence the carbon cycle process of land and ocean, thus cause considerable change in carbon flux distribution. In addition, correlation study also demonstrated the possible important role of Asian terrestrial ecosystems in carbon cycle.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 10773002, 10875012, and 11175019)the Fundamental Research Funds for the Central Universities, China (Grant No. 105116)
文摘The recent work of Nation et al., in which the Hawking radiation energy and entropy flow from a black hole is considered to be produced in a one-dimensional Landauer transport process, is extended to the case of a Reissner- Nordstrom black hole. The energy flow contains not only the contribution of the thermal flux but also that of the particle flux. It is found that the charge can also be transported via the one-dimensional quantum tunnel. Because of the existence of the electrostatic potential, the entropy production rate is shown to be smaller than that of the Schwarzschild black hole.
基金Supported by the Fundamental Research Funds for the Central Universities under Grant No lzujbky-2014-170the Research Fund for the Doctoral Program of Higher Education of China under Grant No 20120211120002the National Natural Science Foundation of China under Grant Nos 11205075 and 11375076
文摘Within the framework of a semiclassical Boltzmann-Uehling-Uhlenbeck （BUU） transport model, the high mo- mentum tail （HMT） effects of nucleon momentum distribution in the nucleus on the nucleon collective flows are studied in semieentral Au＋Au collisions. The HMT due to the isospin-dependent short-range correlations causes a smaller value of the collective flows. We find that the HMT effects on the nucleon collective flows are remarkable at beam energy of 300 MeV/nucleon and become weak as the incident beam energy increases. The results indicate that for the collective flow studies at intermediate energies, the HMT of nucleon momentum distribution in nucleus should be taken into account in transport models.
基金Supported by the National Natural Science Foundation of China under Grant Nos.90921010,11174373,11005086,11205125, 2010CB923104the Youth Foundation of Education Department of Sichuan Province under Grant No.11ZB032the Sichuan Youth Science and Technology Foundation under Grant No.2011JQ0019
文摘We investigate the Hawking radiation of a Kaluza-Klein black hole by using one-dimensional(1D),nonequilibrium,Landauer transport model.The derived Hawking radiation temperature is in consistence with that obtained by using the usual anomaly method.With the Landauer transport model,we calculate the entropy flow out of the Kaluza-Klein black hole and the relevant entropy production rate.How these quantities depending on the physical parameters of the black hole is also discussed.
文摘In the high-energy environment of coastal seas and estuaries,strong sediment resuspension/ deposition events are driven by surface waves,tides,winds and buoyancy driven currents.In recent years,A POM based three-dimensional,wave-current coupled,sediment transport model has been developed by the University of New South Wales.This paper presents several examples of the model applications to study sediment dynamics in the environments where forcings such as waves,tides,and winds are equally important to affect sediment fluxes and distributions.Firstly,the sediment transport model coupled to the Yellow Sea general circulation model and a third generation wave model SWAN was implemented in the Yellow Sea to study the dynamics of the sediment transport and resuspension in the northern Jiangsu shoal-wate(rNJSW).The sediment distributions and fluxes and their inter-annual variability were studied by realistic numerical simulations.The study found that the surface waves played a dominant role over the tides to form the turbidity maxima along the muddy coast of NJSW. Secondly,the sediment transport model was used to explore the effect of suspended sediment-induced stratification in the bottom boundary laye(rBBL).The model uses a re-parameterized bottom drag coefficient Cd that incorporates a linear stability function of flux Richardson number Rf.The study has shown that the sediment induced stratification in the BBL reduces the vertical eddy viscosity and bottom shear stress in comparison with the model prediction in a neutrally stratified BBL.In response to these apparent reductions,the tidal current shear is increased and sediments are abnormally concentrated within a thin wall layer that is overlain by a thicker layer with much smaller concentration.The formation of this fluid-mud layer near the seabed has led to a significant reduction in the total sediment transport.This study contributes to the understanding of formations of tidal flats along the coasts of turbid seas and estuaries.
文摘Acid mine drainage（AMD） is commonly treated by neutralization with alkaline substances.This treatment is supported by titration experiments that illustrate the buffering mechanisms and estimate the base neutralization capacity（BNC） of the AMD.Detailed explanation of titration curves requires modeling with a hydro-chemical model.In this study the titration curves of water samples from the drainage of the As Pontes mine and the corresponding dumps have been investigated and six buffers are selerted by analyzing those curves.Titration curves have been simulated by a reactive transport model to discover the detailed buffering mechanisms.These simulations show seven regions involving different buffering mechanism.The BNC is primarily from buffers of dissolved Fe,Al and hydrogen sulfate.The BNC can be approximated by：BNC = 3（C_（Fe） ＋ C_（Al）） ＋ 0.05C_（sulfate）,where the units are mol/L.The BNC of the sample from the mine is 9.25×10^（-3） mol/L and that of the dumps sample is 1.28×10^（-2） mol/L
基金Supported by the National Natural Science Foundation of China under Grant No U1332125the Program for Innovation Research of Science in Harbin Institute of Technology under Grant No B201408
文摘The elliptic flow v2, for π±, K±, p and p in Au＋Au collisions at center-of-mass energies √sNN=7.7, 11.5, 14.5 and 19.6 GeV, is analyzed using a multiphase transport model. A significant difference in the v2 values for p and p is observed, and the values of v2 splitting are larger compared with π＋ and π-, K＋ and K-. The difference increases with decreasing the center-of-mass energy. The effect of the quark coalescence mechanism in a multi-phase transport model to the value of elliptic difference △v2 between p and p- has been discussed. The simulation of Au＋Au collisions at 14.5 GeV shows that the effect of hadron cascade to △v2 is not obvious, and a larger patton-scattering cross section can lead to a larger △v2.
基金supported by the National Natural Science Foundation of China(Nos.11605163 and 21504085)the China Academy of Engineering Physics Foundation for Development of Science and Technology(No.201580103014 and No.2015B0301063)+1 种基金the Foundation for Special Talents in China Academy of Engineering Physics(No.TP201502-3)the Sichuan Science and Technology Development Foundation for Young Scientists(No.2017Q0050)
文摘Identifying the unknown geometric and material information of a multi-shield object by analyzing the radiation signature measurements is always an important problem in national and global security. In order to identify the unknown shielding layer thicknesses of a source/shield system with gamma-ray spectra, we have developed a derivative-free inverse radiation transport model based on a differential evolution algorithm with global and local neighbourhoods(IRT-DEGL). In the present paper, the IRT-DEGL model is further extended for estimating the unknown thicknesses with random initial guesses and material mass densities of multi-shielding layers as well as their combinations. Using the detected gamma-ray spectra,the illustration of inverse studies is implemented and the main influence factors for inverse results are also analyzed.
