A 3-D numerical model is set up in a large domain covering the Hangzhou Bay and the Changjiang Estuary based on the ECOM model in orthogonal curvilinear coordinates.The numerical schemes for baroclinic pressure gradie...A 3-D numerical model is set up in a large domain covering the Hangzhou Bay and the Changjiang Estuary based on the ECOM model in orthogonal curvilinear coordinates.The numerical schemes for baroclinic pressure gradient (BPG)terms and convective terms are improved in the paper according to the characteristics of velocity field and mass transport in the area.The model is validated by the simulations of residual current and salinity transport in the Hangzhou Bay and the Changjiang Estuary.展开更多
Active faults are a common adverse geological phenomenon that can occur during tunnel excavation and has a very negative impact on the construction and operation of the tunnel.In this paper,the grade IV rock surroundi...Active faults are a common adverse geological phenomenon that can occur during tunnel excavation and has a very negative impact on the construction and operation of the tunnel.In this paper,the grade IV rock surrounding the cross-fault tunnel with poor geological conditions has been chosen for the study.The support capacity of 2^(nd) Generation-Negative Poisson’s Ratio(2G-NPR)bolt in an active fault tunnel has been carried out on the basis of relevant results obtained from the geomechanical model test and numerical investigations of failure model for existing unsupported fault tunnel.The investigation shows that surrounding rock of the tunnel is prone to shear deformation and crack formation along the fault,as a result,the rock mass on the upper part of the fault slips as a whole.Furthermore,small-scale deformation and loss of blocks are observed around the tunnel;however,the 2G-NPR bolt support is found to be helpful in keeping the overall tunnel intact without any damage and instability.Due to the blocking effect of fault,the stress of the surrounding rock on the upper and lower parts of the fault is significantly different,and the stress at the left shoulder of the tunnel is greater than that at the right shoulder.The asymmetrical arrangement of 2G-NPR bolts can effectively control the asymmetric deformation and instability of the surrounding rock.The present numerical scheme is in good agreement with the model test results,and can reasonably reflect the stress and displacement characteristics of the surrounding rock of the tunnel.In comparison to unsupported and ordinary PR(Poisson’s Ratio)bolt support,2G-NPR bolt can effectively limit the fault slip and control the stability of the surrounding rock of the fault tunnel.The research findings may serve as a guideline for the use of 2G-NPR bolts in fault tunnel support engineering.展开更多
Till now, there have been lots of researches on numerical simulation of vortex-induced vibration. Acceptable results have been obtained for fixed cylinders with low Reynolds number. However, for responses of 2-DOF vor...Till now, there have been lots of researches on numerical simulation of vortex-induced vibration. Acceptable results have been obtained for fixed cylinders with low Reynolds number. However, for responses of 2-DOF vortex-induced vibration with low mass ratio, the accuracy is not satisfactory, especially for the maximum amplitudes. In Jauvtis and Williamson＇s work, the maximum amplitude of the cylinder with low mass ratio m＊=2.6 can reach as large as 1.5D to be called as the ＂super-upper branch＂, but from current literatures, few simulation results can achieve such value, even fail to capture the upper branch. Besides, it is found that the amplitude decays too fast in the lower branch with the RANS-based turbulence model. The reason is likely to be the defects of the turbulence model itself in the prediction of unsteady separated flows as well as the unreasonable setting of the numerical simulation parameters. Aiming at above issues, a modified turbulence model is proposed in this paper, and the effect of the acceleration of flow field on the response of vortex-induced vibration is studied based on OpenFOAM. By analyzing the responses of amplitude, phase and trajectory, frequency and vortex mode, it is proved that the vortex-induced vibration can be predicted accurately with the modified turbulence model under appropriate flow field acceleration.展开更多
Debris flows are rapid mass movements with a mixture of rock,soil and water.High-intensity rainfall events have triggered multiple debris flows around the globe,making it an important concern from the disaster managem...Debris flows are rapid mass movements with a mixture of rock,soil and water.High-intensity rainfall events have triggered multiple debris flows around the globe,making it an important concern from the disaster management perspective.This study presents a numerical model called debris flow simulation 2D(DFS 2D)and applicability of the proposed model is investigated through the values of the model parameters used for the reproduction of an occurred debris flow at Yindongzi gully in China on 13 August 2010.The model can be used to simulate debris flows using three different rheologies and has a userfriendly interface for providing the inputs.Using DFS 2D,flow parameters can be estimated with respect to space and time.The values of the flow resistance parameters of model,dry-Coulomb and turbulent friction,were calibrated through the back analysis and the values obtained are 0.1 and 1000 m/s^(2),respectively.Two new methods of calibration are proposed in this study,considering the crosssectional area of flow and topographical changes induced by the debris flow.The proposed methods of calibration provide an effective solution to the cumulative errors induced by coarse-resolution digital elevation models(DEMs)in numerical modelling of debris flows.The statistical indices such as Willmott's index of agreement,mean-absolute-error,and normalized-root-mean-square-error of the calibrated model are 0.5,1.02 and 1.44,respectively.The comparison between simulated and observed values of topographic changes indicates that DFS 2D provides satisfactory results and can be used for dynamic modelling of debris flows.展开更多
In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubi...In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubic spline numerical model(Spline Model for short),which is with a quasi-Lagrangian time-split integration scheme of fitting cubic spline/bicubic surface to all physical variable fields in the atmospheric equations on spherical discrete latitude-longitude mesh.A new algorithm of"fitting cubic spline—time step integration—fitting cubic spline—……"is developed to determine their first-and2nd-order derivatives and their upstream points for time discrete integral to the governing equations in Spline Model.And the cubic spline function and its mathematical polarities are also discussed to understand the Spline Model’s mathematical foundation of numerical analysis.It is pointed out that the Spline Model has mathematical laws of"convergence"of the cubic spline functions contracting to the original functions as well as its 1st-order and 2nd-order derivatives.The"optimality"of the 2nd-order derivative of the cubic spline functions is optimal approximation to that of the original functions.In addition,a Hermite bicubic patch is equivalent to operate on a grid for a 2nd-order derivative variable field.Besides,the slopes and curvatures of a central difference are identified respectively,with a smoothing coefficient of 1/3,three-point smoothing of that of a cubic spline.Then the slopes and curvatures of a central difference are calculated from the smoothing coefficient 1/3 and three-point smoothing of that of a cubic spline,respectively.Furthermore,a global simulation case of adiabatic,non-frictional and"incompressible"model atmosphere is shown with the quasi-Lagrangian time integration by using a global Spline Model,whose initial condition comes from the NCEP reanalysis data,along with quasi-uniform latitude-longitude grids and the so-called"shallow atmosphere"Navier-Stokes primitive equations in the spherical coordinates.The Spline Model,which adopted the Navier-Stokes primitive equations and quasi-Lagrangian time-split integration scheme,provides an initial ideal case of global atmospheric circulation.In addition,considering the essentially non-linear atmospheric motions,the Spline Model could judge reasonably well simple points of any smoothed variable field according to its fitting spline curvatures that must conform to its physical interpretation.展开更多
In this paper, condensation of water vapor from a mixture of COR2R/HR2RO is studied numerically. To simplify the study and focus on the physical model, a simple vertical plate was chosen. Two condensation models are d...In this paper, condensation of water vapor from a mixture of COR2R/HR2RO is studied numerically. To simplify the study and focus on the physical model, a simple vertical plate was chosen. Two condensation models are developed and numerical approach is considered to implement these models. The main objective in the cur-rent paper was to study the capability of numerical modeling in prediction of complex process. Results showed that developed condensation models in combination with numerical approach can predict the trends in condensation behavior of binary mixture very well. Results from this study can be developed further to be used in design of condensers which are suitable for oxy-fuel power plants.展开更多
This paper points out that a successful numerical simulation is to construct a correct conceptional model which is very dose to the natural condition. A new model, two dimensional coupled three dimensional model (2D -...This paper points out that a successful numerical simulation is to construct a correct conceptional model which is very dose to the natural condition. A new model, two dimensional coupled three dimensional model (2D -3D ) is presented in the Present paper,which is the most suitable one for the dual - structured - aquifer system. An example of Wenyinghu area is shown.By using the 2D-3D model, a satisfied result of the simulated area is achieved.展开更多
In order to protect the city of New Orleans from the Mississippi River flooding, the Bonnet Carré Spillway (BCS) was constructed from 1929 to 1936 to divert flood water from the river into Lake Pontchartrain and ...In order to protect the city of New Orleans from the Mississippi River flooding, the Bonnet Carré Spillway (BCS) was constructed from 1929 to 1936 to divert flood water from the river into Lake Pontchartrain and then into the Gulf of Mexico. During the BCS opening for flood release, large amounts of freshwater, nutrients, sediment, etc. were discharged into Lake Pontchartrain, and caused a lot of environmental problems. To evaluate the environmental impacts of the flood water on lake ecosystems, a two-dimensional numerical model was developed based on CCHE2D and applied to simulate the flow circulation, sediment transport and algal biomass distribution in Lake Pontchartrain. The effect of sediment concentration on the growth of algae was considered in the model. The numerical model was calibrated using field measured data provided by USGS, and then it was validated by the BCS Opening Event in 1997. The simulated results were generally in good agreement with filed data and satellite imagery. The field observation and numerical model show that during the spillway opening for flood release, the sediment concentration is very high, which greatly restricts the growth of algae, so there is no algal bloom observed in the lake. After the closure of BCS, the sediment concentration in the lake reduces gradually, and the nutrient concentration of the lake is still high. Under these conditions, numerical results and satellite imagery showed that the chlorophyll concentration was high and algal bloom might occur.展开更多
Atmospheric models are physical equations based on the ideal gas law. Applied to the atmosphere, this law yields equations for water, vapor (gas), ice, air, humidity, dryness, fire, and heat, thus defining the model o...Atmospheric models are physical equations based on the ideal gas law. Applied to the atmosphere, this law yields equations for water, vapor (gas), ice, air, humidity, dryness, fire, and heat, thus defining the model of key atmospheric parameters. The distribution of these parameters across the entire planet Earth is the origin of the formation of the climatic cycle, which is a normal climatic variation. To do this, the Earth is divided into eight (8) parts according to the number of key parameters to be defined in a physical representation of the model. Following this distribution, numerical models calculate the constants for the formation of water, vapor, ice, dryness, thermal energy (fire), heat, air, and humidity. These models vary in complexity depending on the indirect trigonometric direction and simplicity in the sum of neighboring models. Note that the constants obtained from the equations yield 275.156˚K (2.006˚C) for water, 273.1596˚K (0.00963˚C) for vapor, 273.1633˚K (0.0133˚C) for ice, 0.00365 in/s for atmospheric dryness, 1.996 in<sup>2</sup>/s for humidity, 2.993 in<sup>2</sup>/s for air, 1 J for thermal energy of fire, and 0.9963 J for heat. In summary, this study aims to define the main parameters and natural phenomena contributing to the modification of planetary climate. .展开更多
In this paper a new approach for PBL simulation, the non-local closure scheme based on the transient turbulence theory has been used. It was set up as an alternative to local closure schemes which physical concept is ...In this paper a new approach for PBL simulation, the non-local closure scheme based on the transient turbulence theory has been used. It was set up as an alternative to local closure schemes which physical concept is reasonable and distinct. A 2-D non-local closure model was developed in order to study the PBL structure and simulatesome interesting atmospheric processes over non-ulliform underlying surface, especially under the convective and unique weather conditions, such as sea-land circulation and the TIBL structure. The modelled results show good agreement with field measurement.展开更多
The development of functional relationships between the observed deposition rate and the experimental conditions is an important step toward understanding and optimizing low-pressure chemical vapor deposition(LPCVD)or...The development of functional relationships between the observed deposition rate and the experimental conditions is an important step toward understanding and optimizing low-pressure chemical vapor deposition(LPCVD)or low-pressure chemical vapor infiltration(LPCVI).In the field of ceramic matrix composites(CMCs),methyltrichlorosilane(CH3 SiCl3,MTS)is the most widely used source gas system for SiC,because stoichiometric SiC deposit can be facilitated at 900°C–1300°C.However,the reliability and accuracy of existing numerical models for these processing conditions are rarely reported.In this study,a comprehensive transport model was coupled with gas-phase and surface kinetics.The resulting gas-phase kinetics was confirmed via the measured concentration of gaseous species.The relationship between deposition rate and 24 gaseous species has been effectively evaluated by combining the special superiority of the novel extreme machine learning method and the conventional sticking coefficient method.Surface kinetics were then proposed and shown to reproduce the experimental results.The proposed simulation strategy can be used for different material systems.展开更多
The Guide sedimentary basin is located in the northeastern part of Qinghai-Xizang Plateau,which is rich in geothermal resources.However,exploitation of the geothermal resources has so far been limited,because of limit...The Guide sedimentary basin is located in the northeastern part of Qinghai-Xizang Plateau,which is rich in geothermal resources.However,exploitation of the geothermal resources has so far been limited,because of limited understanding of the resources quantity and storage gained from scientific researches.In this study,using a typical cross section across the basin and taking into account its geothermal and geological conditions,a new waterheat coupled model was built and associated modelling was done by the software TOUGH2.During modelling process,the accuracy and applicability of the model was confirmed through the calibration of relevant parameters for modelling the heat and water transport and the formation of geothermal reservoir across the basin,with particular focus on the Neogene geothermal field.Results show that the groundwater that flows from the basin margins to the center is heated by the Neogene and Paleogene sedimentary rocks with high geothermal gradients.Since the east-west extending fault F1 is conductive,it acts as preferential flow paths which on one hand provide additional and rapid flows to the thermal reservoir;and on the other hand,cool down the thermal water to a certain extent due to the infiltration of shallower water sources in the vicinity of the fault.Furthermore,the estimated geothermal resources quantity is close to that of previous studies.In comparison with the Paleogene rock formations,the Neogene geothermal reservoir shows a better nature in terms of water content,aquifer permeability and resources exploitability,although the resource quantity of the Paleogene reservoir is considerable.展开更多
We highlighted the flexibility of using unstructured mesh together with the local refinement by a resistivity model with complicated topography. The effect of topography is emphasized. Based on this, we calculated a s...We highlighted the flexibility of using unstructured mesh together with the local refinement by a resistivity model with complicated topography. The effect of topography is emphasized. Based on this, we calculated a specific class of layered models and found that the accuracy is not always satisfactory by utilizing the standard approach. As an improvement, we employed the layered earth as the reference model to calculate the wavenumbers. The comparison demonstrates that the accuracy is considerably improved by using this enhanced approach.展开更多
In this paper, we present a set of best practices for workflow design and implementation for numerical weather prediction models and meteorological data service, which have been in operation in China Meteorological Ad...In this paper, we present a set of best practices for workflow design and implementation for numerical weather prediction models and meteorological data service, which have been in operation in China Meteorological Administration (CMA) for years and have been proven effective in reliably managing the complexities of large-scale meteorological related workflows. Based on the previous work on the platforms, we argue that a minimum set of guidelines including workflow scheme, module design, implementation standards and maintenance consideration during the whole establishment of the platform are highly recommended, serving to reduce the need for future maintenance and adjustment. A significant gain in performance can be achieved through the workflow-based projects. We believe that a good workflow system plays an important role in the weather forecast service, providing a useful tool for monitoring the whole process, fixing the errors, repairing a workflow, or redesigning an equivalent workflow pattern with new components.展开更多
Many parameters have been indicated crucial for the selection of a saline aquifer as a carbon dioxide(CO_(2))storage site.However,less attention has been given to the impact of heterogeneity on the performance of thes...Many parameters have been indicated crucial for the selection of a saline aquifer as a carbon dioxide(CO_(2))storage site.However,less attention has been given to the impact of heterogeneity on the performance of these storage media.Thus,the heterogeneity effect was evaluated in this paper by adopting a numerical modeling approach and the existing screening criterion developed for the aquifers was updated.The updated criterion for CO_(2)storage purpose would enhance the confidence level during the selection of deep saline aquifer and thus,help to address the climate change issue.The numerical modeling was carried out via CO_(2)STORE module of Eclipse300 Simulator to evaluate the effect of different levels of heterogeneity on CO_(2)storage potential.Different degrees of heterogeneity from homogenous systems to highly heterogeneous systems in the model were incorporated through the Lorenz coefficient.In this way,simulation of nine cases was carried out for three different aquifers with different porosity values.A comparison of these results showed that heterogeneity causes the aquifer to have lower storage capacity.On the trapping potential,dissolution trapping was significant and the amount of free gas in all cases was minimum.In addition,the aquifer with the highest level of heterogeneity(HLH)had a minimum fraction of residual trapping regardless of porosity.It was also found that final pressure at the end of 30 years is the same and high for low-level heterogeneity(LLH)and medium level heterogeneity(MLH)cases and low for HLH,while the injection rate stability duration is least for HLH and maximum for LLH.Based on the results obtained,it can be concluded that low to medium level heterogeneous aquifers with a good porosity can be a suitable choice for CO_(2)storage.展开更多
Usually a buffer layer of cadmium sulphide is used in high efficiency solar cells based on Cu(In,Ga)Se2(CIGS). Because of cadmium toxicity, many in-vestigations have been conducted to use Cd-free buffer layers. Our wo...Usually a buffer layer of cadmium sulphide is used in high efficiency solar cells based on Cu(In,Ga)Se2(CIGS). Because of cadmium toxicity, many in-vestigations have been conducted to use Cd-free buffer layers. Our work focuses on this type of CIGS-based solar cells where CdS is replaced by a ZnS buffer layer. In this contribution, AFORS-HET software is used to simulate n-ZnO: Al/i-ZnO/n-ZnS/p-CIGS/Mo polycrystalline thin-film solar cell where the key parts are p-CIGS absorber layer and n-ZnS buffer layer. The characteristics of these key parts: thickness and Ga-content of the absorber layer, thickness of the buffer layer and doping concentrations of absorber and buffer layers have been investigated to optimize the conversion efficiency. We find a maximum conversion efficiency of 26% with a short-circuit current of 36.9 mA/cm2, an open circuit voltage of 824 mV, and a fill factor of 85.5%.展开更多
Interest in understanding the structural behavior of marine floating photovoltaic(FPV)systems has grown significantly over the last decade.Numerical models are the preferred approach for understanding FPV responses un...Interest in understanding the structural behavior of marine floating photovoltaic(FPV)systems has grown significantly over the last decade.Numerical models are the preferred approach for understanding FPV responses under environmental loads,but they require validation.Several methods are commonly used to validate numerical results,such as comparison with analytical,field data,and experimental data.The use of analytical approaches to validate numerical results can sometimes be inaccurate due to the complexity of the problems;nevertheless,field data is commonly restricted and frequently unavailable for numerical model validation.Thus,physical models play a crucial role in validating numerical results.This study focuses on the two-dimensional(2-D)modeling process and sensors development for an FPV system with taut mooring,aiming to investigate wave-structure interaction while considering hydroelastic effects.The model is developed in accordance with the Froude-Cauchy similitude law and is made from composite materials to capture structural stiffness.Structural motions,specifically heave and pitch,are measured using an Inertial Measurement Unit(IMU),while strain gauges measure structural stress and mooring tension.The sensors provide precise measurements for strain and pitch;however,heave,as a result of time-domain integration from acceleration,requires further validation.The motion responses of the model align with reference results.展开更多
Geological storage and utilization of CO_(2)involve complex interactions among Thermo-hydromechanical-chemical(THMC)coupling processes,which significantly affect storage integrity and efficiency.To address the challen...Geological storage and utilization of CO_(2)involve complex interactions among Thermo-hydromechanical-chemical(THMC)coupling processes,which significantly affect storage integrity and efficiency.To address the challenges in accurately simulating these coupled phenomena,this paper systematically reviews recent advances in the mathematical modeling and numerical solution of THMC coupling in CO_(2)geological storage.The study focuses on the derivation and structure of governing and constitutive equations,the classification and comparative performance of fully coupled,iteratively coupled,and explicitly coupled solution methods,and the modeling of dynamic changes in porosity,permeability,and fracture evolution induced by multi-field interactions.Furthermore,the paper evaluates the capabilities,application scenarios,and limitations of major simulation platforms,including TOUGH,CMG-GEM,and COMSOL.By establishing a comparative framework integrating model formulations and solver strategies,this work clarifies the strengths and gaps of current approaches and contributes to the development of robust,scalable,and mechanism-oriented numerical models for long-term prediction of CO_(2)behavior in geological formations.展开更多
文摘A 3-D numerical model is set up in a large domain covering the Hangzhou Bay and the Changjiang Estuary based on the ECOM model in orthogonal curvilinear coordinates.The numerical schemes for baroclinic pressure gradient (BPG)terms and convective terms are improved in the paper according to the characteristics of velocity field and mass transport in the area.The model is validated by the simulations of residual current and salinity transport in the Hangzhou Bay and the Changjiang Estuary.
基金supported by the National Natural Science Foundation of China(NSFC)(41941018)the Program of China Scholarship Council(202106430031)。
文摘Active faults are a common adverse geological phenomenon that can occur during tunnel excavation and has a very negative impact on the construction and operation of the tunnel.In this paper,the grade IV rock surrounding the cross-fault tunnel with poor geological conditions has been chosen for the study.The support capacity of 2^(nd) Generation-Negative Poisson’s Ratio(2G-NPR)bolt in an active fault tunnel has been carried out on the basis of relevant results obtained from the geomechanical model test and numerical investigations of failure model for existing unsupported fault tunnel.The investigation shows that surrounding rock of the tunnel is prone to shear deformation and crack formation along the fault,as a result,the rock mass on the upper part of the fault slips as a whole.Furthermore,small-scale deformation and loss of blocks are observed around the tunnel;however,the 2G-NPR bolt support is found to be helpful in keeping the overall tunnel intact without any damage and instability.Due to the blocking effect of fault,the stress of the surrounding rock on the upper and lower parts of the fault is significantly different,and the stress at the left shoulder of the tunnel is greater than that at the right shoulder.The asymmetrical arrangement of 2G-NPR bolts can effectively control the asymmetric deformation and instability of the surrounding rock.The present numerical scheme is in good agreement with the model test results,and can reasonably reflect the stress and displacement characteristics of the surrounding rock of the tunnel.In comparison to unsupported and ordinary PR(Poisson’s Ratio)bolt support,2G-NPR bolt can effectively limit the fault slip and control the stability of the surrounding rock of the fault tunnel.The research findings may serve as a guideline for the use of 2G-NPR bolts in fault tunnel support engineering.
基金financially supported by the National Natural Science Foundation of China(Grant No.51509045)
文摘Till now, there have been lots of researches on numerical simulation of vortex-induced vibration. Acceptable results have been obtained for fixed cylinders with low Reynolds number. However, for responses of 2-DOF vortex-induced vibration with low mass ratio, the accuracy is not satisfactory, especially for the maximum amplitudes. In Jauvtis and Williamson＇s work, the maximum amplitude of the cylinder with low mass ratio m＊=2.6 can reach as large as 1.5D to be called as the ＂super-upper branch＂, but from current literatures, few simulation results can achieve such value, even fail to capture the upper branch. Besides, it is found that the amplitude decays too fast in the lower branch with the RANS-based turbulence model. The reason is likely to be the defects of the turbulence model itself in the prediction of unsteady separated flows as well as the unreasonable setting of the numerical simulation parameters. Aiming at above issues, a modified turbulence model is proposed in this paper, and the effect of the acceleration of flow field on the response of vortex-induced vibration is studied based on OpenFOAM. By analyzing the responses of amplitude, phase and trajectory, frequency and vortex mode, it is proved that the vortex-induced vibration can be predicted accurately with the modified turbulence model under appropriate flow field acceleration.
基金financially supported by Department of Space,India(Grant No.ISRO/RES/4/663/18-19)。
文摘Debris flows are rapid mass movements with a mixture of rock,soil and water.High-intensity rainfall events have triggered multiple debris flows around the globe,making it an important concern from the disaster management perspective.This study presents a numerical model called debris flow simulation 2D(DFS 2D)and applicability of the proposed model is investigated through the values of the model parameters used for the reproduction of an occurred debris flow at Yindongzi gully in China on 13 August 2010.The model can be used to simulate debris flows using three different rheologies and has a userfriendly interface for providing the inputs.Using DFS 2D,flow parameters can be estimated with respect to space and time.The values of the flow resistance parameters of model,dry-Coulomb and turbulent friction,were calibrated through the back analysis and the values obtained are 0.1 and 1000 m/s^(2),respectively.Two new methods of calibration are proposed in this study,considering the crosssectional area of flow and topographical changes induced by the debris flow.The proposed methods of calibration provide an effective solution to the cumulative errors induced by coarse-resolution digital elevation models(DEMs)in numerical modelling of debris flows.The statistical indices such as Willmott's index of agreement,mean-absolute-error,and normalized-root-mean-square-error of the calibrated model are 0.5,1.02 and 1.44,respectively.The comparison between simulated and observed values of topographic changes indicates that DFS 2D provides satisfactory results and can be used for dynamic modelling of debris flows.
文摘In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubic spline numerical model(Spline Model for short),which is with a quasi-Lagrangian time-split integration scheme of fitting cubic spline/bicubic surface to all physical variable fields in the atmospheric equations on spherical discrete latitude-longitude mesh.A new algorithm of"fitting cubic spline—time step integration—fitting cubic spline—……"is developed to determine their first-and2nd-order derivatives and their upstream points for time discrete integral to the governing equations in Spline Model.And the cubic spline function and its mathematical polarities are also discussed to understand the Spline Model’s mathematical foundation of numerical analysis.It is pointed out that the Spline Model has mathematical laws of"convergence"of the cubic spline functions contracting to the original functions as well as its 1st-order and 2nd-order derivatives.The"optimality"of the 2nd-order derivative of the cubic spline functions is optimal approximation to that of the original functions.In addition,a Hermite bicubic patch is equivalent to operate on a grid for a 2nd-order derivative variable field.Besides,the slopes and curvatures of a central difference are identified respectively,with a smoothing coefficient of 1/3,three-point smoothing of that of a cubic spline.Then the slopes and curvatures of a central difference are calculated from the smoothing coefficient 1/3 and three-point smoothing of that of a cubic spline,respectively.Furthermore,a global simulation case of adiabatic,non-frictional and"incompressible"model atmosphere is shown with the quasi-Lagrangian time integration by using a global Spline Model,whose initial condition comes from the NCEP reanalysis data,along with quasi-uniform latitude-longitude grids and the so-called"shallow atmosphere"Navier-Stokes primitive equations in the spherical coordinates.The Spline Model,which adopted the Navier-Stokes primitive equations and quasi-Lagrangian time-split integration scheme,provides an initial ideal case of global atmospheric circulation.In addition,considering the essentially non-linear atmospheric motions,the Spline Model could judge reasonably well simple points of any smoothed variable field according to its fitting spline curvatures that must conform to its physical interpretation.
文摘In this paper, condensation of water vapor from a mixture of COR2R/HR2RO is studied numerically. To simplify the study and focus on the physical model, a simple vertical plate was chosen. Two condensation models are developed and numerical approach is considered to implement these models. The main objective in the cur-rent paper was to study the capability of numerical modeling in prediction of complex process. Results showed that developed condensation models in combination with numerical approach can predict the trends in condensation behavior of binary mixture very well. Results from this study can be developed further to be used in design of condensers which are suitable for oxy-fuel power plants.
文摘This paper points out that a successful numerical simulation is to construct a correct conceptional model which is very dose to the natural condition. A new model, two dimensional coupled three dimensional model (2D -3D ) is presented in the Present paper,which is the most suitable one for the dual - structured - aquifer system. An example of Wenyinghu area is shown.By using the 2D-3D model, a satisfied result of the simulated area is achieved.
文摘In order to protect the city of New Orleans from the Mississippi River flooding, the Bonnet Carré Spillway (BCS) was constructed from 1929 to 1936 to divert flood water from the river into Lake Pontchartrain and then into the Gulf of Mexico. During the BCS opening for flood release, large amounts of freshwater, nutrients, sediment, etc. were discharged into Lake Pontchartrain, and caused a lot of environmental problems. To evaluate the environmental impacts of the flood water on lake ecosystems, a two-dimensional numerical model was developed based on CCHE2D and applied to simulate the flow circulation, sediment transport and algal biomass distribution in Lake Pontchartrain. The effect of sediment concentration on the growth of algae was considered in the model. The numerical model was calibrated using field measured data provided by USGS, and then it was validated by the BCS Opening Event in 1997. The simulated results were generally in good agreement with filed data and satellite imagery. The field observation and numerical model show that during the spillway opening for flood release, the sediment concentration is very high, which greatly restricts the growth of algae, so there is no algal bloom observed in the lake. After the closure of BCS, the sediment concentration in the lake reduces gradually, and the nutrient concentration of the lake is still high. Under these conditions, numerical results and satellite imagery showed that the chlorophyll concentration was high and algal bloom might occur.
文摘Atmospheric models are physical equations based on the ideal gas law. Applied to the atmosphere, this law yields equations for water, vapor (gas), ice, air, humidity, dryness, fire, and heat, thus defining the model of key atmospheric parameters. The distribution of these parameters across the entire planet Earth is the origin of the formation of the climatic cycle, which is a normal climatic variation. To do this, the Earth is divided into eight (8) parts according to the number of key parameters to be defined in a physical representation of the model. Following this distribution, numerical models calculate the constants for the formation of water, vapor, ice, dryness, thermal energy (fire), heat, air, and humidity. These models vary in complexity depending on the indirect trigonometric direction and simplicity in the sum of neighboring models. Note that the constants obtained from the equations yield 275.156˚K (2.006˚C) for water, 273.1596˚K (0.00963˚C) for vapor, 273.1633˚K (0.0133˚C) for ice, 0.00365 in/s for atmospheric dryness, 1.996 in<sup>2</sup>/s for humidity, 2.993 in<sup>2</sup>/s for air, 1 J for thermal energy of fire, and 0.9963 J for heat. In summary, this study aims to define the main parameters and natural phenomena contributing to the modification of planetary climate. .
文摘In this paper a new approach for PBL simulation, the non-local closure scheme based on the transient turbulence theory has been used. It was set up as an alternative to local closure schemes which physical concept is reasonable and distinct. A 2-D non-local closure model was developed in order to study the PBL structure and simulatesome interesting atmospheric processes over non-ulliform underlying surface, especially under the convective and unique weather conditions, such as sea-land circulation and the TIBL structure. The modelled results show good agreement with field measurement.
基金the National Key R&D Program of China(Grants No.2017YFB0703200)National Natural Science Foundation of China(Grants Nos.51702100,51972268)China Postdoctoral Science Foundation(Grants No.2018M643075)for financial support。
文摘The development of functional relationships between the observed deposition rate and the experimental conditions is an important step toward understanding and optimizing low-pressure chemical vapor deposition(LPCVD)or low-pressure chemical vapor infiltration(LPCVI).In the field of ceramic matrix composites(CMCs),methyltrichlorosilane(CH3 SiCl3,MTS)is the most widely used source gas system for SiC,because stoichiometric SiC deposit can be facilitated at 900°C–1300°C.However,the reliability and accuracy of existing numerical models for these processing conditions are rarely reported.In this study,a comprehensive transport model was coupled with gas-phase and surface kinetics.The resulting gas-phase kinetics was confirmed via the measured concentration of gaseous species.The relationship between deposition rate and 24 gaseous species has been effectively evaluated by combining the special superiority of the novel extreme machine learning method and the conventional sticking coefficient method.Surface kinetics were then proposed and shown to reproduce the experimental results.The proposed simulation strategy can be used for different material systems.
基金This work was supported by grants from National Natural Science Foundation of China(No.41402231)geothermal survey project of China Geological Survey(No.DD20190128)+1 种基金Natural Science Foundation of Hebei Province China(No.D2019330003)S&T Program of Hebei China(No.20374201D).
文摘The Guide sedimentary basin is located in the northeastern part of Qinghai-Xizang Plateau,which is rich in geothermal resources.However,exploitation of the geothermal resources has so far been limited,because of limited understanding of the resources quantity and storage gained from scientific researches.In this study,using a typical cross section across the basin and taking into account its geothermal and geological conditions,a new waterheat coupled model was built and associated modelling was done by the software TOUGH2.During modelling process,the accuracy and applicability of the model was confirmed through the calibration of relevant parameters for modelling the heat and water transport and the formation of geothermal reservoir across the basin,with particular focus on the Neogene geothermal field.Results show that the groundwater that flows from the basin margins to the center is heated by the Neogene and Paleogene sedimentary rocks with high geothermal gradients.Since the east-west extending fault F1 is conductive,it acts as preferential flow paths which on one hand provide additional and rapid flows to the thermal reservoir;and on the other hand,cool down the thermal water to a certain extent due to the infiltration of shallower water sources in the vicinity of the fault.Furthermore,the estimated geothermal resources quantity is close to that of previous studies.In comparison with the Paleogene rock formations,the Neogene geothermal reservoir shows a better nature in terms of water content,aquifer permeability and resources exploitability,although the resource quantity of the Paleogene reservoir is considerable.
基金supported by the National High Technology Research and Development Program of China (863 Program) (No. 2007AA06Z134)the National Natural Science Foundation of China (No. 40874072)
文摘We highlighted the flexibility of using unstructured mesh together with the local refinement by a resistivity model with complicated topography. The effect of topography is emphasized. Based on this, we calculated a specific class of layered models and found that the accuracy is not always satisfactory by utilizing the standard approach. As an improvement, we employed the layered earth as the reference model to calculate the wavenumbers. The comparison demonstrates that the accuracy is considerably improved by using this enhanced approach.
文摘In this paper, we present a set of best practices for workflow design and implementation for numerical weather prediction models and meteorological data service, which have been in operation in China Meteorological Administration (CMA) for years and have been proven effective in reliably managing the complexities of large-scale meteorological related workflows. Based on the previous work on the platforms, we argue that a minimum set of guidelines including workflow scheme, module design, implementation standards and maintenance consideration during the whole establishment of the platform are highly recommended, serving to reduce the need for future maintenance and adjustment. A significant gain in performance can be achieved through the workflow-based projects. We believe that a good workflow system plays an important role in the weather forecast service, providing a useful tool for monitoring the whole process, fixing the errors, repairing a workflow, or redesigning an equivalent workflow pattern with new components.
文摘Many parameters have been indicated crucial for the selection of a saline aquifer as a carbon dioxide(CO_(2))storage site.However,less attention has been given to the impact of heterogeneity on the performance of these storage media.Thus,the heterogeneity effect was evaluated in this paper by adopting a numerical modeling approach and the existing screening criterion developed for the aquifers was updated.The updated criterion for CO_(2)storage purpose would enhance the confidence level during the selection of deep saline aquifer and thus,help to address the climate change issue.The numerical modeling was carried out via CO_(2)STORE module of Eclipse300 Simulator to evaluate the effect of different levels of heterogeneity on CO_(2)storage potential.Different degrees of heterogeneity from homogenous systems to highly heterogeneous systems in the model were incorporated through the Lorenz coefficient.In this way,simulation of nine cases was carried out for three different aquifers with different porosity values.A comparison of these results showed that heterogeneity causes the aquifer to have lower storage capacity.On the trapping potential,dissolution trapping was significant and the amount of free gas in all cases was minimum.In addition,the aquifer with the highest level of heterogeneity(HLH)had a minimum fraction of residual trapping regardless of porosity.It was also found that final pressure at the end of 30 years is the same and high for low-level heterogeneity(LLH)and medium level heterogeneity(MLH)cases and low for HLH,while the injection rate stability duration is least for HLH and maximum for LLH.Based on the results obtained,it can be concluded that low to medium level heterogeneous aquifers with a good porosity can be a suitable choice for CO_(2)storage.
文摘Usually a buffer layer of cadmium sulphide is used in high efficiency solar cells based on Cu(In,Ga)Se2(CIGS). Because of cadmium toxicity, many in-vestigations have been conducted to use Cd-free buffer layers. Our work focuses on this type of CIGS-based solar cells where CdS is replaced by a ZnS buffer layer. In this contribution, AFORS-HET software is used to simulate n-ZnO: Al/i-ZnO/n-ZnS/p-CIGS/Mo polycrystalline thin-film solar cell where the key parts are p-CIGS absorber layer and n-ZnS buffer layer. The characteristics of these key parts: thickness and Ga-content of the absorber layer, thickness of the buffer layer and doping concentrations of absorber and buffer layers have been investigated to optimize the conversion efficiency. We find a maximum conversion efficiency of 26% with a short-circuit current of 36.9 mA/cm2, an open circuit voltage of 824 mV, and a fill factor of 85.5%.
基金funded by the Center for Higher Education Funding and Assessment(Pusat Pelayanan Pembiayaan dan Asesmen Pendidkan Tinggi)-Ministry of Higher Education,Science,and Technology of Republic Indonesia.
文摘Interest in understanding the structural behavior of marine floating photovoltaic(FPV)systems has grown significantly over the last decade.Numerical models are the preferred approach for understanding FPV responses under environmental loads,but they require validation.Several methods are commonly used to validate numerical results,such as comparison with analytical,field data,and experimental data.The use of analytical approaches to validate numerical results can sometimes be inaccurate due to the complexity of the problems;nevertheless,field data is commonly restricted and frequently unavailable for numerical model validation.Thus,physical models play a crucial role in validating numerical results.This study focuses on the two-dimensional(2-D)modeling process and sensors development for an FPV system with taut mooring,aiming to investigate wave-structure interaction while considering hydroelastic effects.The model is developed in accordance with the Froude-Cauchy similitude law and is made from composite materials to capture structural stiffness.Structural motions,specifically heave and pitch,are measured using an Inertial Measurement Unit(IMU),while strain gauges measure structural stress and mooring tension.The sensors provide precise measurements for strain and pitch;however,heave,as a result of time-domain integration from acceleration,requires further validation.The motion responses of the model align with reference results.
基金supported by the China Postdoctoral Science Foundation(No.2024M752803)the National Natural Science Foundation of China(No.52179112)the Open Fund of National Key Laboratory of Oil and Gas Reservoir Geology and Exploitation(Southwest Petroleum University)(No.PLN2023-02)。
文摘Geological storage and utilization of CO_(2)involve complex interactions among Thermo-hydromechanical-chemical(THMC)coupling processes,which significantly affect storage integrity and efficiency.To address the challenges in accurately simulating these coupled phenomena,this paper systematically reviews recent advances in the mathematical modeling and numerical solution of THMC coupling in CO_(2)geological storage.The study focuses on the derivation and structure of governing and constitutive equations,the classification and comparative performance of fully coupled,iteratively coupled,and explicitly coupled solution methods,and the modeling of dynamic changes in porosity,permeability,and fracture evolution induced by multi-field interactions.Furthermore,the paper evaluates the capabilities,application scenarios,and limitations of major simulation platforms,including TOUGH,CMG-GEM,and COMSOL.By establishing a comparative framework integrating model formulations and solver strategies,this work clarifies the strengths and gaps of current approaches and contributes to the development of robust,scalable,and mechanism-oriented numerical models for long-term prediction of CO_(2)behavior in geological formations.
