Numerical Models and Methods of Atmospheric Parameters Originating in the Formation of the Earth’s Climatic Cycle

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 in2/s for humidity, 2.993 in2/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. .

2-D Modeling and Calculations of Stratospheric Ozone and Influences of Convection, Diffusion, and Time

An engineering system approach of 2-D cylindrical model of transient mass balance calculations of ozone and other concerned chemicals along with fourteen photolysis, ozone-generating and ozone-depleting chemical reaction equations was developed, validated, and used for studying the ozone concentrations, distribution and peak of the layer, ozone depletion and total ozone abundance in the stratosphere. The calculated ozone concentrations and profile at both the Equator and a 60˚N location were found to follow closely with the measured data. The calculated average ozone concentration was within 1% of the measured average, and the deviation of ozone profiles was within 14%. The monthly evolution of stratospheric ozone concentrations and distribution above the Equator was studied with results discussed in details. The influences of slow air movement in both altitudinal and radial directions on ozone concentrations and profile in the stratosphere were explored and discussed. Parametric studies of the influences of gas diffusivities of ozone DO3 and active atomic oxygen DO on ozone concentrations and distributions were also studied and delineated. Having both influences through physical diffusion and chemical reactions, the diffusivity (and diffusion) of atomic oxygen DO was found to be more sensitive and important than that of ozone DO3 on ozone concentrations and distribution. The 2-D ozone model present in this paper for stratospheric ozone and its layer and depletion is shown to be robust, convenient, efficient, and executable for analyzing the complex ozone phenomena in the stratosphere. .

Numerical Analysis of Steady Smoldering of Biomass Rods

Understanding the steady mechanism of biomass smoldering plays a great role in the utilization of smoldering technology.In this study numerical analysis of steady smoldering of biomass rods was performed.A two-dimensional(2D)steady model taking into account both char oxidation and pyrolysis was developed on the basis of a calculated propagation velocity according to empirical correlation.The model was validated against the smoldering experiment of biomass rods under natural conditions,and the maximum error was smaller than 31%.Parameter sensitivity analysis found that propagation velocity decreases significantly while oxidation area and pyrolysis zone increase significantly with the increasing diameter of rod fuel.

Numerical study on residual current and its impact on mass transport in the Hangzhou Bay and the Changjiang Estuary I. A3-D joint model of the Hangzhou Bay and the Changjiang Estuary

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.

Geomechanics model test and numerical simulation of 2G-NPR bolt support effect in an active fault tunnel

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.

Two Improvements on Numerical Simulation of 2-DOF Vortex-Induced Vibration with Low Mass Ratio

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 flow simulation 2D(DFS 2D):Numerical modelling of debris flows and calibration of friction parameters

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.

AN INTEGRATION METHOD WITH FITTING CUBIC SPLINE FUNCTIONS TO A NUMERICAL MODEL OF 2ND-ORDER SPACE-TIME DIFFERENTIAL REMAINDER——FOR AN IDEAL GLOBAL SIMULATION CASE WITH PRIMITIVE ATMOSPHERIC EQUATIONS

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.

Investigation on Numerical Modeling of Water Vapour Condensation from a Flue Gas with High COR2R Content

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.

2D Coupled 3D:A New Numerical Model for Dual - Structured - Aquifer System

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.

Numerical Modeling of Sediment Transport and Its Effect on Algal Biomass Distribution in Lake Pontchartrain Due to Flood Release from Bonnet CarréSpillway

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.

A 2-D Non-local Closure Model for Atmospheric Boundary Layer Simulations

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.

Numerical modeling of SiC by low-pressure chemical vapor deposition from methyltrichlorosilane

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.

Using TOUGH2 numerical simulation to analyse the geothermal formation in Guide basin,China

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.

Workflow-Based Platform Design and Implementation for Numerical Weather Prediction Models and Meteorological Data Service

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.

A numerical study to assess the effect of heterogeneity on CO_(2) storage potential of saline aquifers

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.

鄂尔多斯盆地铁边城地区长2油层组剩余油模拟

为了实现油田增储上产及老井挖潜的目标,提升剩余油资源开发潜力,通过岩心、薄片、电镜、物性、模拟方法系统的研究了鄂尔多斯盆地西北部姬塬油田铁边城地区长2段主力油层的剩余油分布特征。研究显示,长2_(1)^(3-2)和长2_(1)^(3-3)小层为主力产油层,平均厚度分别为12.2 m及11.1 m。目的层砂岩储层的孔隙度通常介于12%~18%,渗透率则通常低于8 mD。根据相控原则,建立了主力产层的物理属性模型及流体模型;进而,开展了剩余油分布数值模拟研究。通过基于实测产水及产油量约束的历史生产曲线拟合,获得了剩余油的纵横向分布结果。研究显示,剩余油分布受砂体边界、物性差异及注采关系的影响。剩余油整体上成片分布;但是,受油藏物性差异的影响,注水多呈不均匀受效,因而造成局部剩余油富集。剩余油主要分布在砂体边界处,这些区域无井网控制,因而水驱程度低,剩余油分布较多。受长期注采的影响,井间会形成高渗区,这些其余类似高速通道,会发生快速水淹,其周边则形成剩余油滞留区。可见,剩余油富集明显受微相、物性、注水及井网因素的控制。

Numerical Modeling and Simulation of CIGS-Based Solar Cells with ZnS Buffer Layer

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%.

Calculations and Sensitivity Analysis of Chlorine-,NO_(x)-,and Bromine-Depleting Cycles of Stratospheric Ozone

This paper presents an engineering system approach using a 2D model of conservation of mass to study the dynamics of ozone and concerned chemical species in the stratosphere.By considering all fourteen photolysis,ozone-generating,and-depleting chemical reactions,the model calculated the transient,spatial changes of ozone under different physical-chemical-radiative conditions.Validation against the measured data demonstrated good accuracy,close match of our model with the observed ozone concentrations at both 20°S and 90°N locations.The deviation in the average concentration was less than 1% and in ozone profiles less than 17%.The impacts of various chlorine-(Cl),nitrogen oxides-(NO_(x)),and bromine-(Br)depleting cycles on ozone concentrations and distribution were investigated.The chlorine catalytic depleting cycle was found to exhibit the most significant impact on ozone dynamics,confirming the key role of chlorine in the problem of ozone depletion.Sensitivity analysis was conducted with levels of 25%,50%,100%,200%,and 400% of the baseline value.The combined cycles(Cl+NO_(x)+Br)showed the most significant influence on ozone behavior.The total ozone abundance above the South Pole could decrease by a small 3%,from 281 DU(Dubson Units)to 273 DU for the 25% level,or by a huge thinning of 60%to 114 DU for the 400% concentration level.When the level of chlorine gases increased beyond 200%,it would cause ozone depletion to a level of ozone hole(below 220 DU).The 2D Ozone Model presented in this paper demonstrates robustness,convenience,efficiency,and executability for analyzing complex ozone phenomena in the stratosphere.

Interannual variability of carbon cycle implied by a 2-D atmospheric transport model

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.