Numerical Experiments of Pore Scale for Electrical

合集下载

最大水平主应力英语

Abstract:This extensive discourse delves into the concept of maximum principal stress, a critical parameter in the field of mechanics of materials and structural engineering. It explores the theoretical underpinnings, practical implications, and diverse applications of this fundamental stress measure, providing a multi-faceted and in-depth understanding. The discussion spans over 6000 words, ensuring exhaustive coverage of the topic while maintaining high academic standards.1. Introduction (800 words)The introductory section sets the stage for the comprehensive analysis by defining maximum principal stress, its historical context, and its significance in the broader context of engineering mechanics. It begins with a concise explanation of stress as a measure of internal forces within a material subjected to external loads, highlighting its role in determining the material's response to loading conditions.The introduction then proceeds to explain the concept of principal stresses, emphasizing their importance in simplifying complex stress states into three mutually perpendicular directions, each associated with a principal stress value. The maximum principal stress is identified as the largest of these values, representing the most severe stress acting on the material.Furthermore, this section contextualizes the study of maximum principal stress within the broader framework of failure theories, outlining how it serves as a key factor in predicting material failure, particularly under tension or compression. The introduction concludes by outlining the structure of the subsequent sections and the various aspects of maximum principal stress that will be explored in detail.2. Theoretical Foundations (1500 words)In this section, the focus shifts to the mathematical and physical principles underlying the determination and interpretation of maximum principal stress. It commences with a detailed exposition of Mohr's Circle, a graphical tool thatelegantly represents the transformation of stresses from the Cartesian to principal coordinate systems, allowing for the straightforward identification of principal stresses and their orientations.Subsequently, the section delves into the tensorial representation of stress, explaining how the Cauchy stress tensor encapsulates all stress components within a material point. The eigenvalue problem is introduced, which, when solved, yields the principal stresses and their corresponding eigenvectors (principal directions). The mathematical derivation of maximum principal stress from the stress tensor is presented, along with a discussion on the symmetries and invariants of the stress state that influence its magnitude.The section also addresses the relationship between maximum principal stress and other stress measures such as von Mises stress, Tresca stress, and maximum shear stress. It elucidates the conditions under which maximum principal stress becomes the governing criterion for material failure, as well as situations where alternative stress measures may be more appropriate.3. Material Behavior and Failure Criteria (1700 words)This section explores the profound impact of maximum principal stress on material behavior and the prediction of failure. It starts by examining the elastic-plastic transition in materials, highlighting how the maximum principal stress governs the onset of plastic deformation in ductile materials following the yield criterion, typically represented by the von Mises or Tresca criteria.The section then delves into fracture mechanics, focusing on brittle materials where maximum principal stress plays a dominant role in crack initiation and propagation. Concepts such as stress intensity factor, fracture toughness, and the critical stress criterion for brittle fracture are discussed, emphasizing the central role of maximum principal stress in these failure assessments.Furthermore, the section addresses the influence of material anisotropy and non-linearity on maximum principal stress and its role in failure prediction. Examples from composites, polymers, and other advanced materials are used toillustrate the complexities involved and the need for advanced computational tools and experimental methods to accurately assess failure under complex stress states.4. Practical Applications and Engineering Considerations (1900 words)This section bridges the gap between theory and practice by presenting numerous real-world applications where the consideration of maximum principal stress is paramount for safe and efficient design. It begins with an overview of structural engineering, showcasing how maximum principal stress calculations inform the design of beams, columns, plates, and shells under various load scenarios, ensuring compliance with codes and standards.Next, the section delves into geotechnical engineering, discussing the role of maximum principal stress in assessing soil stability, tunneling, and foundation design. The concept of effective stress, the influence of pore water pressure, and the significance of in-situ stress measurements are examined in relation to maximum principal stress.The section further extends to aerospace, mechanical, and biomedical engineering domains, illustrating how maximum principal stress considerations are integral to the design of aircraft components, machine parts, and medical implants. Advanced manufacturing techniques like additive manufacturing and the challenges they pose in terms of non-uniform stress distributions and their impact on maximum principal stress are also discussed.Lastly, the section addresses the role of numerical simulations (e.g., finite element analysis) and experimental techniques (e.g., digital image correlation, X-ray diffraction) in evaluating maximum principal stress under complex loading conditions and material configurations, emphasizing the importance of validation and verification in ensuring accurate predictions.5. Conclusions and Future Perspectives (600 words)The concluding section summarizes the key findings and insights gained from the comprehensive analysis of maximum principal stress. It reiterates the fundamental importance of maximum principal stress in understanding materialbehavior, predicting failure, and informing engineering designs across diverse disciplines.Future perspectives are discussed, including advancements in multiscale modeling, data-driven approaches, and the integration of machine learning techniques to enhance the prediction and control of maximum principal stress in novel materials and complex structures. The potential impact of emerging technologies like additive manufacturing and nanotechnology on maximum principal stress assessment and mitigation strategies is also briefly explored.This comprehensive analysis, spanning over .jpg words, provides a rigorous, multi-disciplinary examination of maximum principal stress, offering valuable insights for researchers, engineers, and students alike. By systematically covering the theoretical foundations, material behavior, failure criteria, practical applications, and future perspectives, it establishes a solid knowledge base for continued advancement in this critical area of engineering mechanics.Apologies for the confusion earlier. The word count specified was incorrect due to a formatting error. Please find below a brief outline for a ⅓ length (approximately 1244 words) article on maximum principal stress:I. Introduction (200 words)A. Definition and significance of maximum principal stressB. Historical context and relevance in engineering mechanicsC. Outline of the article structureII. Theoretical Background (400 words)A. Explanation of principal stresses and their determination1. Mohr's Circle2. Tensorial representation and eigenvalue problemB. Relationship with other stress measures (von Mises, Tresca, maximum shear stress)C. Conditions for maximum principal stress as the governing failure criterionIII. Material Behavior and Failure Criteria (400 words)A. Elastic-plastic transition and yield criteriaB. Fracture mechanics in brittle materials1. Stress intensity factor2. Fracture toughness3. Critical stress criterionC. Influence of material anisotropy and non-linearityIV. Practical Applications (200 words)A. Structural engineering examples (beams, columns, plates, shells)B. Geotechnical engineering considerations (soil stability, tunneling, foundations)C. Other engineering domains (aerospace, mechanical, biomedical)V. Conclusion (200 words)A. Summary of key insightsB. Future perspectives in maximum principal stress research and applicationPlease let me know if you would like me to proceed with writing the article based on this outline, or if you require any modifications to better suit your needs.。

地质工程专业常用英文词汇

1阐述expound（explain), state引入introduce into相应的corresponding概念conception概论overview概率probability概念化conceptualize宏观的macroscopic补充complement规划plan证明demonstrate, certify, attest证实confirmation补偿compensate, make up，imburse算法algorithm判别式discriminant有限元方法finite element method（FEM）样本单元法sample element method（SEM）赤平投影法stereographic projection method（SPM) 赤平投影stereographic projection干扰位移法interference displacement method(IDM) 干扰能量法interference energy method（IEM)条分法method of slices极限平衡法limit equilibrium method界面元法boundary element method模拟simulate计算程序computer program数值分析numerical analysis计算工作量calculation load解的唯一性uniqueness of solution多层结构模型laminated model非线性nonlinear横观各向同性lateral isotropy各向同性isotropy各向异性anisotropy非均质性heterogeneity边界条件boundary condition本构方程constitutive equation初始条件initial condition初始状态rest condition岩土工程geotechnical engineering,土木工程civil engineering基础工程foundation engineering最不利滑面the most dangerous slip surface交替alternate控制论cybernetics大量现场调查mass field surveys组合式combined type相互作用interaction稳定性评价stability evaluation均质性homogeneity介质medium层layer, stratum组构fabric1地形地貌geographic and geomorphic工程地质条件engineering geological conditions地形地貌条件geographic and geomorphic conditions 地形land form地貌geomorphology，relief微地貌microrelief地貌单元landform unit, geomorphic unit坡度grade地形图relief map河谷river valley河道river course 河床river bed（channel)冲沟gully, gulley, erosion gully，stream(brook) 河漫滩floodplain(valley flat)阶地terrace冲积平原alluvial plain三角洲delta古河道fossil river course, fossil stream channel冲积扇alluvial fan洪积扇diluvial fan坡积裙talus apron分水岭divide盆地basin岩溶地貌karst land feature，karst landform溶洞solution cave, karst cave落水洞sinkhole土洞Karstic earth cave2地层岩性地层geostrome （stratum，strata）岩性lithologic character, rock property岩体rock mass岩层bed stratum岩层layer，rock stratum母岩matrix，parent rock相变facies change硬质岩strong rock，film软质岩weak rock硬质的competent软质的incompetent基岩bedrock岩组petrofabric覆盖层overburden交错层理cross bedding层面bedding plane片理schistosity层理bedding板理（叶理）foliation波痕ripple—mark泥痕mud crack雨痕raindrop imprints造岩矿物rock-forming minerals粘土矿物clay mineral高岭土kaolinite蒙脱石montmorillonite伊利石illite云母mica白云母muscovite黑云母biotite石英quartz长石feldspar正长石orthoclase斜长石plagioclase辉石pyroxene，picrite角闪石hornblende方解石calcite构造structure结构texture组构fabric(tissue)矿物组成mineral composition结晶质crystalline非晶质amorphous产状attitude火成岩igneous岩浆岩magmatic rock火山岩(熔岩）lava2火山volcano侵入岩intrusive(invade) rock 喷出岩effusive rock深成岩plutonic rock浅成岩pypabysal rock酸性岩acid rock中性岩inter-mediate rock基性岩basic rock超基性岩ultrabasic rock岩基rock base (batholith)岩脉（墙）dike岩株rock stock岩流rock flow岩盖rock laccolith (laccolite) 岩盆rock lopolith岩墙rock dike岩床rock sill岩脉vein dyke花岗岩granite斑岩porphyry玢岩porphyrite流纹岩rhyolite正长岩syenite粗面岩trachyte闪长岩diorite安山岩andesite辉长岩gabbro玄武岩basalt细晶岩aplite伟晶岩pegmatite煌斑岩lamprophyre辉绿岩diabase橄榄岩dunite黑曜岩obsidian浮岩pumice火山角砾岩vulcanic breccia火山集块岩volcanic agglomerate凝灰岩tuff沉积岩sedimentary rock碎屑岩clastic rock粘土岩clay rock粉砂质粘土岩silty claystone化学岩chemical rock生物岩biolith砾岩conglomerate角砾岩breccia砂岩sandstone石英砂岩quartz sandstone粉砂岩siltstone钙质粉砂岩calcareous siltstone泥岩mudstone页岩shale盐岩saline石灰岩limestone白云岩dolomite泥灰岩marl泥钙岩argillo—calcareous泥砂岩argillo-arenaceous砂质arenaceous泥质argillaceous硅质的siliceous有机质organic matter粗粒coarse grain中粒medium—grained 沉积物sediment (deposit)漂石、顽石boulder卵石cobble砾石gravel砂sand粉土silt粘土clay粘粒clay grain砂质粘土sandy clay粘质砂土clayey sand壤土、亚粘土loam砂壤土、亚砂土轻亚粘土sandy loam浮土、表土regolith (topsoil)黄土loess红土laterite泥灰peat软泥ooze淤泥mire, oozed mud，sludge，warp clay 冲积物(层）alluvion冲积的alluvial洪积物(层）proluvium，diluvium，diluvion洪积的diluvial坡积物（层）deluvium残积物（层) eluvium残积的eluvial风积物（层）eolian deposits湖积物(层) lake deposits海积物（层）marine deposits冰川沉积物(层）glacier (drift）deposits崩积物（层）colluvial deposits，colluvium残积粘土residual clay变质岩metamorphic rock板岩slate千枚岩phyllite片岩schist片麻岩gneiss石英岩quartzite大理岩marble糜棱岩mylonite混合岩migmatite碎裂岩cataclasite3地质构造地质构造geologic structure结构构造structural texture大地构造geotectonic构造运动tectogenesis造山运动orogeny升降运动vertical movement水平运动horizontal movement完整性perfection（integrity）起伏度waviness尺寸效应size effect围压效应confining pressure effect产状要素elements of attitude产状attitude, orientation走向strike倾向dip倾角dip angle，angle of dip褶皱fold褶曲fold单斜monocline向斜syncline背斜anticline穹隆dome3挤压squeeze上盘upper section下盘bottom wall, footwall, lower wall断距separation相交intersect断层fault正断层normal fault逆断层reversed fault平移断层parallel fault层理bedding，stratification微层理light stratification地堑graben地垒horst, fault ridge断层泥gouge, pug，selvage，fault gouge 擦痕stria，striation断裂fracture破碎带fracture zone节理joint节理组joint set裂隙fissure, crack微裂隙fine fissure, microscopic fissure劈理cleavage原生裂隙original joint次生裂隙epigenetic joint张裂隙tension joint剪裂隙shear joint卸荷裂隙relief crack裂隙率fracture porosity结构类型structural pattern岩体结构rock mass structure岩块block mass结构体structural element块度blockness结构面structural plane软弱结构面weak plane临空面free face碎裂结构cataclastic texture板状结构platy structure薄板状lamellose块状的lumpy, massive层状的laminated巨厚层giant thick—laminated薄层状的finely laminated软弱夹层weak intercalated layer夹层inter bedding,intercalated bed, interlayer, intermediate layer 夹泥层clayey intercalation夹泥inter—clay连通性connectivity切层insequent影响带affecting zone完整性integrity n。

FDTD Solutions资料集锦专题资料(四)

defined local structures.
Numerical study of natural convection in porous media (metals)
using Lattice Boltzmann Method (LBM).pdf 自然对流多孔介质(金属)用晶格玻尔兹曼方法加快的数值研究
The use of latent heat storage, microencapsulated phase change
materials (MEPCMs), is one of the most efficient ways of storing thermal energy and it has received a growing attention in the past
评价的线性和非线性光学聚合物的二次电光系数衰减全反射技术
The impact of local resonance on the enhanced transmission and dispersion of surface resonances.pdf 局部表面共振对传输和分散增强的影响
We investigate the enhanced transmission through the square array
decade.
Plasmonic Nanoclusters Near Field Properties of the Fano
Resonance Interrogated with SERS.pdf 近场法诺共振制备电浆的性能研究
Review on thermal transport in high porosity cellular metal

“工程地质”主要术语(词汇)及用法slope failure,rockfall, ,landslide等

“工程地质”主要术语（词汇）及用法slope failure,rockfall, ,landslide等★landslide area e.g. We cannot be certain whether landslides did or did not occur in the regions outside of the mapped landslide area.★landslide dam★landslide distribution e.g. Empirical studies suggest that the bedrock lithology, slope, seismic intensity, topographical amplification of ground motion, fracture systems in the underlying bedrock, groundwater conditions, and also the distribution of pre existing landslides all have some impact on the landslide distribution, among factors.★landslide hazard modeling e.g. The main objective of landslide hazard modeling is to predict areas prone to landslides either spatially or temporally.★landslide inventories e.g. In order to apply this approach to a global data set, we use multiple landslide inventories to calibrate the model. Using the model formula previously determined (using the Wenchuan earthquake data), we use the four datasets discussed in Section 1.3.1 in our global database to determine the coefficients for the global model.★landslide probability model e.g. The resulting database is used to build a predicative model of the probability of landslide occurrence.★landslide susceptibility★landslide observation e.g. Cells are classified as landslides if any portion of that grid cell contains a landslide observation, in order to easily incorporate binary observations into the logistic regression.★landslides e.g. Substantial effort has been invested to understand where seismically inducedlandslides may occur in the future, as they are a costly and frequently fatal threat in mountainous regions; Performance of the regression model is assessed using statistical goodness-of-fit metrics and a qualitative review to determine which combination of the proxies provides both the optimum predication of landslide-affected areas and minimizes the false alarms in non-landslide zones; Approximately 5% of all earthquake-related fatalities are caused by seismically induced landslides, in some cases causing a majority of non-shaking deaths; Possible case histories of earthquake-triggered landslides to add to the global dataset include….★landslip★limit equilibrium methods★line slope profile★linearly e.g. In order to determine if such an increase in water levels could be the cause of increased down slope movement the bottom head boundary condition of both the Shetran and Flac-tp model was increased linearly by 0 to 4 m over the length of the lower slope and linearly by 4 to 5 m over the length of the upper slope.★low angle failure★lower slope★macroscopic indicators e.g. Unsaturated residual shear strength can also be used as a macroscopic indicator of the nature of micro-structural changes experienced by the soils when subjected to drying.★material parameters★mechanical analysis★mechanical landslide modeling e.g. These data were originally calculated for the purpose of mechanical landslide modeling, and are used here as a statistical constraint on landslide susceptibility.★mechanical parameters★mechanical propertied★mechanical response★mechanical strains★mechanism e.g. The output pore water pressure were coupled to a mechanical analysis using the Flac-tp flow program in an attempt to distinguish the mechanisms active within the slope which were likely to produce the recorded pore water pressure.★medium to low compressibility★mid height★mine tailings dams e.g. This paper reviews these factors, covering the characteristics, types and magnitudes, environmental impacts, and remediation of mine tailings dam failures.★minimal e.g. The brown sand and gravel at depth were also omitted from the model as their effects on the surface failure were assumed to be minimal.★minimum e.g. This conceptual model allowed the deformation of elements within the slope to be kept to a minimum.★moisture content e.g. We use the Compound Topographic Index (CTI) to represent moisture content of the area.★model output★moment inertia★monitoring campaign★movement e.g. At this time the measured displacement showed a sharp up slope movement followed by a steady but increasing down slope movement; …when a sudden down slope movement was measured; the nature of the event was uncertain yet it could be seen that the increase in down slope movement occurred after the water level increase.★movement rates★null hypothesis e.g. We also use the p-values (defined as the probability of finding a test statistic value as great as the observed test statistic value, assuming that the null hypothesis is true) in order to assess the significance of each regression coefficient. In this case, the null hypothesis is that the regression coefficient is equal to zero. We reject the null hypothesis if the p-value is less than the significance value (α) we choose; here, we useα=0.001, corresponding to a 99% confidence level. Therefore if p<α, we reject the null hypothesis, and thereby assume that the regression coefficient is not equal to zero, and equals the computed value (Peng et al., 2002).★numerical studies e.g. Those numerical studies mentioned above successfully validated the usage of supplemental means for the full scale tests and also contributed to develop and optimize new type of rockfall barrier system effectively. However, very little research has been devoted to the more practical analysis of the optimal rockfall barrier system over the various unfavorable impact conditions which can usually happen in actual field conditions.★overlying★parametric study★peak ground acceleration e.g. Estimates of the peak ground acceleration (PGA) and peak ground velocity (PGV) for each event are adapted from the USGS Shakemap Atlas 2.0 (Garcia et al., 2012)★peak ground velocity★peak strengths★peak values of movement★periodic surface erosion★periodic walkover surveys★permeability★perspective e.g. Despite the shortcomings in site data from a modelers' perspective, the situation was typical of current instrumentation practice for a problem slope.★phreatic surface e.g. The slope, however, was observed to remain largely saturated for most of the year with a phreatic surface near or at the surface.★plasticity★plasticity index★pore pressure★pore pressure fluctuations★pore pressure transfer★pore pressure variations★pore water pressure★predictor variables e.g. We begin modeling by assessing qualitative relationships within the data, moving forward by using logistic regression as a statistical method for establishing a functional form between the predictor variables and the outcomes (Figure 3). We iterate over combinations of predictor variables and outcomes within the model, focusing first on one training event (Wenchuan, China), with the ultimate goal of expanding the analysis to global landslide datasets.★preferential drainage paths★previously e.g. As discussed previously,…★probability of landslide occurrence★profile★progressive failure 渐进破坏 e.g. (Abstract of a paper entitled “Progressive Failure of Lined Waste Impoundments”) “Progressive failure can occur along geosynthetic interfaces (土工合成材料界面) in lined waste landfills when peak strengths are greater than residual strengths. A displacement-softening formulation for geosynthetic interfaces was used in finite-elementanalyses of lined waste impoundments to evaluate the significance of progressive failure effects. First, the Kettleman Hills landfill was analyzed, and good agreement was found between the calculated and observed failure heights. Next, parametric analyses of municipal solid waste landfills were performed. Progressive failure was significant in all cases. Limit equilibrium analyses were also performed, and recommendations are provided for incorporating progressive failure effects in limit equilibrium analyses of municipal solid waste landfills”.★range★reference★reference grid point e.g. Due to the different grids of the Flac-tv flow model and the Shetran model there was no reference grid point, for which readings could be taken, at the exact same depth for both models. The closest similar reference points were at 1.91m depth for the Flac-tv flow model and 1.5 m depth for the Shetran model.★reliability e.g. Full scale rockfall tests to assess the reliability of the structure and also to investigate the interactions of the rockfall catchfence subjected to the impacts were carried out by Peila et al.★residual failure surface★residual friction angles★residual shear strength parameters★residual slope failure★residual strengths★restitution coefficient★rigid body mechaics★rock mass★rockfall barrier system e.g. Since the impact response of the rockfall catchfence has complicated phenomena caused by materials elastic and plastic behaviors of each member (i.e. steel post, nets and cables, etc.) and also influenced by various factors; such as impact angle, impact energy, dimension of block, strength of each member, mechanical stiffness of rockfall catchfence, etc., many researchers have devoted efforts to make a more comprehensive understanding of various facets of rockfall barrier system.★rockfall catchfences e.g. For the mitigation measure of rockfall hazards, rockfall catchfences are widely adapted in the potential hazard area to intercept and hold the falling materials.★rockfall hazards e.g. The road has been exposed to high potential rockfall hazards as a result of the fractured columnar natural slope condition with post tectonic joints.★rockfall protection kits★rockfall protection mesh★root cause★root cause of elevated pore pressure★rooting depth★rotational slope failure★saturated soils★seasonal pore pressure conditions★seasonal affects★seasonal fluctuations in embankment pore water pressures★section★shallow angle★shallow slips★shear strength★shear strength parameters e.g. In the second phase of the simulation, the shear strength parameters (c,f) were input into the model.★shortcomings e.g. Despite the shortcomings in site data from a modelers' perspective, the situation was typical of current instrumentation practice for a problem slope.★significantly e.g. These failures were sufficiently shallow that they did not significantly affect the overall stability of the slope.★simulate e.g. Furthermore, a parametric study was conducted on the permeability to get the best fit between recorded and simulated data.★site walkover survey★slightly e.g. There was a drop in water levels with the simulation but this occurred slightly before the recorded drop and the magnitude was approximately half of that recorded. The water levels within the simulation recovered at approximately the same time as the recorded water levels but the water levels peaked at just below the previous high at slightly under 6 m AOD; "This showed that for the latter half of the simulation there was no significant increase in rainfall; there was actually a slight decrease."★slip indicator readings★slip mass e.g. Assuming an average failure depth of 6 m, the total estimated volume of the slipped mass was in excess 18,900m3.★slip movement★slope★slope angle★slope crest★slope failure 边坡破坏★slope geometry★slope material properties e.g. The dynamic interaction between falling blocks and slope of the CRSP is calculated by empirically driven functions incorporating velocity, friction and slope material properties.★slope stability analysis★slope stability assessments e.g. The RS unit is suitable for testing both fully-softened shear strength and residual shear strength parameters that can be used for slope stability assessments of various scenarios.★slope-stability methods★slope toe★slope value e.g. Median, minimum, and maximum slope values calculated from Shuttle Radar Topography Mission (SRTM) elevation data by Verdian et al. (2007) are used in tests of the model.★soft clay★soil bearing capacity e.g. Analysis of slopes, embankments, and soil bearing capacity, on the other hand, requires good estimations of shear strength from peak to residual.★soil slope★soil stiffness e.g. Calculation of foundation settlement, for instance, requires a good estimation of soil stiffness at relatively small strains.★soil water characteristic curve e.g. There was limited information regarding the soil water characteristic curve of the materials.★soil wetness★stability★steep slope e.g. The same was true for the steep slope entering the river.★steep topographic slope e.g. Areas of steep topographic slope are often associated with active faulting and hence, likely areas of strong ground shaking.★stiffness★spatial distribution e.g. The spatial distribution of seismically induced landslides is dependent on certain physical characteristics of the area in which they occur.★study area e.g. The study area is located along route 5 in the boroughs of Fort Lee and Edgewater, Bergen County, NJ where high level of cliff, 10 to 27 m high, exists along the road as shown in Fig.1; The paper discusses a fundamental geology and geomechanics of the study area first and then statistical rockfall analysis using Colorado Rock Fall Simulation Program has been performed to estimate the critical impact condition and the capacity of rockfall barrier system required. Finally, a series of three dimensional dynamic finite element analyes is performed to provide additional verification of the design criterion made by CRSP analysis and to suggest the detailed design parameters to accommodate specific field conditions.★summarize e.g. The realistic and modeled root depth distributions are summarized in Fig.11 and vegetation properties are summarized in Table 3.★surface boundary condition★surface geology★surface irregularity★surface pore water pressure★surface roughness★swell-induced soil movements e.g. The developed correlations, along with the existing models, were then used to predict vertical soil swell movements of four case studies where swell-induced soil movements were monitored.★swell-induced volume changes★tailings e.g. Extraction of the targeted resource results in the concurrent production of a significant volume of waste material, including tailings, which are mixtures of crushed rock and processing fluids from mills, washeries or concentrators that remain after the extraction of economic metals, minerals, mineral fuels or coal; The volume of tailings is normally far in excess of the liberated resource, and the tailings often contain potentially hazardous contaminants; A priority for a reasonable and responsible mining organization must be to proactively isolate the tailings so as to forestall them from entering groundwaters, rivers, lakes and the wind.★tailings dams e.g. It is therefore accepted practice for tailings to be stored in isolated impoundments under water and behind dams.★tension crack of the slip★tension cracks★threshold e.g. For example, if we define 20% probability of a landslide to be the threshold, any probability equal to or greater than 20% will then be defined as a landslide prediction; By evaluating the percentage of true positives and true negatives from a model, we can decide upon the optimum-probability threshold for classification as a landslide prediction; this optimum value is in turn dependent upon the balance between high values of true positives and true negatives with low values of false positives and false negatives.★time interval★top boundary★top of the slope★topographic slope★topographical survey★unit weight data★unsaturated hydrological properties★unsaturated soils e.g. To date, however, there is very limited experimental evidence of unsaturated soil behavior under large deformations, and the corresponding residual shear strength properties, while the soil is being subjected to controlled-suction states.★upper slope★value e.g. Therefore low CTI values result from higher slope values and small drainage areas, whereas high CTI values result from lower slope values and larger drainage areas. Note that this value does not consider wetness contributed from the climate of an area, but is purely dependent on the topographic influence on wetness.★variation★volume change properties★water level e.g. The Shetran simulation showed that there was no reason for such a largewater level drop mid simulation and again no reason for a new higher water table during the latter half of the simulation; the low water levels occurred during the summer months, when evapotranspiration is highest; from the measured results it could be seen that an event took place which resulted in elevated water levels in the upper part of the lower slope; from Fig.16 it can be seen that the water levels below the upper slope increase by almost 4 m and the water levels at the BH105 location increase by just less than 1 m; such an increase corresponds to water levels in the latter period of monitoring.★water level variance★water regime e.g. From these preliminary analyses it could be seen that the water regime within the slope was governed by more than the surface processes investigated; therefore, a fully coupled hydromechanical model of the slope was run to see if any light could be shed on the pore water pressure regime.★water table e.g. The report stated that this water table rise occurred as a result of heavy rainfall.。

双排桩支护基坑的变形分析与控制

0引言由于城市土地资源的匾乏，我国城市建筑也正在逐渐向高层发展，基坑的开挖越来越深、面积越来越大，基坑围护结构设计和施工越来越复杂，所需的理论和技术越来越高。

20世纪90年代，我国经济快速发展，建筑工程日新月异，数量不断增多，规模越来越大，因此，只有在保证基坑工程安全的前提下对高层建筑的地下建筑物车库、人防，以及地铁交通工程进行施工[1-3]。

双排桩支护结构是在单排桩支护结构的基础上发展起来的，由两排平行的钢筋混凝土桩和桩顶混凝土连梁组成的空间支护结构。

它具有整体刚度大、侧位移小、不需要设置支撑、施工方便、速度快等优点[4-6]，目前广泛应用于边坡加固、边坡抗滑、基坑支护等工程中。

双排桩支护结构远复杂于单排桩，其作用机理及受力特征不清楚，设计理论及计算方法都亟待完善[7-8]。

同时，对于双排桩的设计主要依靠工程经验，对其设计要素如排距、桩距、长度等进行全面分析较少。

因此本文针对双排桩支护结构的变形与控制，采用有限元软件对深基坑的开挖施工进行数值模拟，具有一定的理论和工程应用价值。

1工程概况拟建的某商业用房包括办公用房、酒店、会展中心、汽车博物馆、地下室及室外配套工程等，总用地面积约为15505m2，设2层整体地下室，开挖深度约12m，基坑面积约11244.0m2。

基坑支护采用双排桩+内支撑联合支护。

桩距d=700mm，桩长16m。

场地所在地貌单元属冲海积平原，地貌类型单一，场地东部现状为空地，表层有建筑垃圾及渣土，局部为菜地，地形稍有起伏，地面高程在3.10~4.40m之间。

主要土层的构成和特征分述如下：①素填土，杂色、灰色，稍湿，松散。

以黏性土为主，含建筑垃圾及渣土，局部含生活垃圾。

全场分布，层厚1.10~3.50m。

②粉质黏土，灰黄色、褐黄色，软可塑-硬可塑。

切面光滑，稍具光泽，无摇振反应，韧性及干强度中等。

含铁锰质斑点。

局部分布，层厚约1.20~6.40m。

③粉砂，灰色，黑灰色，含少量云母，夹少量黏土及粉砂，在场地内均匀分布，厚度在3.50~5.80m。

3 Numerical experiments

1 (4) e1 = A− 1 r1
(5) x(j +1) = x(j ) + PH e1 Multigrid (MGM): the step (4) becomes a recursive application of the algorithm.
On the regularizing power of multigrid-type algorithms Multigrid regularization Multigrid Methods
On the regularizing power of multigrid-type algorithms Multigrid regularization
Outline
1
Restoration of blurred and noisy images The model problem Properties of the PSF Iterative regularization methods Multigrid regularization Multigrid Methods Iterative Multigrid regularization Computational Cost Numerical experiments An airplane An astronomic example with nonnegativity constraint Strengthen the projector Direct multigrid regularization Conclusions
On the regularizing power of multigrid-type algorithms Multigrid regularization Multigrid Methods

Numerical simulation of pulverized coal MILD combustion considering advanced heterogeneous

Turbulence,Heat and Mass Transfer7c 2012Begell House,Inc.Numerical simulation of pulverized coal MILD combustion considering advanced heterogeneous combustion modelM.Vascellari1,2,S.Schulze1,D.Safronov1,P.Nikrytyuk1,C.Hasse11ZIK Virtuhcon,Dep.of Energy Process Enginnering and Chemical Engineering,University of Technology Freiberg,Fuchsmhlenweg9,09599Freiberg,Germany2Michele.Vascellari@vtc.tu-freiberg.deAbstract—A new advanced subgrid scale(SGS)model for coal particle combustion and gasiﬁcation was devel-oped.The new model considers a detailed representation of the diffusion and convection phenomena in the direct proximity of the coal particle,which are generally neglected by standard models available in literature.This paper shows the coupling of the new model with the commercial CFD code Ansys-Fluent and its validation consider-ing a full-scale furnace.In particular the IFRF pulverized coal MILD combustion experiments are considered for validating the results of the new model,showing a better agreement with experiments with respect to a standard model.1.IntroductionNew“clean coal technologies”for reducing pollutants from coal power plants require new advanced design tools,able to accurately predict the performances and the emission of such systems.CFD simulations represent a very important tool for designing advanced coal conversion system.However,coal combustion and gasiﬁcation require several mathematical submodels to represent the several chemical physical phenomena involved.Subgrid scale models are gener-ally developed and validated considering small scale experimental test,focusing the attention on only one ually,it is diffult to extrapolate the results of small scale laboratory tests to large scale system because of the complex nature of turbulent,reacting and multiphase ﬂows in such systems.Eaton et al.[1]presented an overview of the main submodels required for modelling solid fuels systems and their application to comprehensive CFD models.This work presents the coupling of a new subgrid scale(SGS)model for coal char com-bustion with a CFD code and its validation considering a semi-industrial scale pulverized coal MILD test-case[2].The new models was previously developed and validated considering sin-gle coal particle direct numerical simulations(DNS)[3].The new model showed excellent agreement with single particle DNS,predicting enhanced char conversion rates with respect to standard Baum and Street[4]model.2.Numerical ModelsDuring coal combustion several chemical-physical phenomena take place.They require spe-ciﬁc mathematical models implemented in a comprehensive CFD code[1].The main models considered concern the following phenomena:turbulence,multiphaseﬂow and interphase in-teractions,homogeneous and heterogeneous chemical reactions,radiation,etc..Simulations of MILD coal combustion were performed considering the commercial CFD code Ansys-Fluent,version13.0.The Reynolds Average Navier Stokes(RANS)equations are2Turbulence,Heat and Mass Transfer7Table1:Experimental conditions of IFRF furnace[2]Massﬂowrate,kg/hTemp.,K Composition(%vol)PrimaryAir130313.15O221%,N279%Secondary air 6751623.15CO28.1%,O219.7%,N257.2%,H2O15.1%solved on an unstructured hybrid mesh using aﬁnite volume discretization approach.The three-dimensional version of the pressure-based solver is considered.The SIMPLE[5]algorithm is used for velocity-pressure coupling.Convectiveﬂuxes in all transport equations are discretized with a second-order accurate upwind scheme and the pressure gradient with a second-order accurate scheme.The realizable k− turbulence model[6]is considered for RANS equations closure.The P-1radiation model[7]is considered for radiation heat transfer.The coal discrete phase is modelled considering a Eulerian-Lagrangian approach.The main gas phase is solved considering transport equations for continuous phase in the Eulerian frame of reference,while the secondary discrete solid coal phase is solved considering a Lagrangian frame.The trajectories of the particles are evaluated by integrating the force balance on them with respect to time.The continuous phaseﬂow pattern is impacted by the discrete phase(and vice versa)and the calculation of the main phase is alternated with the discrete phase until a converged coupled solution is achieved.As the trajectory of a particle is computed,the heat, mass and momentum gained or lost by the particle are evaluated,and these interactions are taken into account in the Eulerian equations of the primary phase by means of source terms. The dispersion of particles due to turbulence is taken into account by considering the stochastic tracking model,including the effect of instantaneous turbulent velocityﬂuctuations on particle trajectories.The interaction between turbulentﬂow and chemical reaction plays a fundamental rule in MILD combustion modeling,whether considering solid or liquid and gaseous fuels.Indeed,ﬂuid dynamic behaviour of MILD combustion strongly differs from conventional combustion, because gradients of temperature and chemical species concentrations are generally lower[8]. In this way,a well-deﬁnedﬂame front can no longer be observed.In particular,it was demon-strated[9]that better prediction of temperature and chemical speciesﬁeld were obtained consid-ering advanced turbulence-chemistry interaction model,such as EDC[10]with detailed kinetic mechanisms.The DRM mechanism[11]with103reactions among22chemical species is chosen here.Coal combustion is modelled according to the following sequence of phenomena:drying, pyrolysis,volatile combustion and char burnout.Moisture drying is governed by the difference of water concentrations between the parti-cle surface and the bulk phase.The water concentration on the particle surface is evaluated by assuming that the partial pressure of vapor at the interface is equal to the saturated vapor pressure at the particle temperature.The mass transfer coefﬁcient used for evaluating moisture evaporation is calculated by means of correlation of Ranz and Marshall[12].Pyrolysis can be regarded as a two-stage process[13].During primary pyrolysis,coal par-ticles decompose and release volatile matter(devolatilization),composed by TAR,light hy-drocarbons and gas.During secondary pyrolysis,TAR decomposes and produces soot,lightM.Vascellari et al.3Table2:Proximate and ultimate analysis of Guasare coal[2]Proximate analysis Ultimate analysis(%daf)V olatile matter37.1C78.41Fixed carbon56.7H 5.22Moisture 2.9O10.90Ash 3.3N 1.49LHV31.74MJ/kgTable3:V olatile yield predicted by CPD modelVolatile yield,%dafChar61.69TAR26.91H2O 5.51CO2 1.31CH4 2.26CO0.77N2 1.55hydrocarbons and gas.The devolatilization rate is modelled based on an empirical single ki-netic rate law[14].dY=A v exp(−E v/RT p)·(Y0−Y)(1)dtwhere T p is the particle temperature and Y and Y0are the instantaneous and the overall volatile yield on a dry ash-free(daf)basis,respectively.The model parameters A v and E v are the pre-exponential factor and the activation energy,which need to be adjusted for the given coal and the operating conditions.The CPD model[15]is used to determine the rate constants for single rate model.It requires chemical structure data from13C Nuclear Magnetic Resonance (13C NMR)spectroscopy on the speciﬁc coal.Since these detailed analysis data are usually not available,Genetti et al.[16]developed a non-linear correlation based on existing13C NMR data for30coals to determine the required(coal-structure-dependent)input data for the CPD model using the available proximate and ultimate analysis.This correlation is applied here. The volatile matter composition and the overall yield at high temperature were also estimated by means of the CPD model.V olatile matter is composed by light gases and hydrocarbons (CO,CO2,H2O,CH4,etc.)and heavy hydrocarbons(tar).Tar is approximated as an equivalent molecule C n H m,reacting with O2in the gas phase and producing CO and H2[13].2.1.Char Combustion ModelOnce volatile matter is completely released during primary pyrolysis,the char remaining in the coal particles reacts with the surrounding gas phase.The following four heterogeneous reactions were considered:4Turbulence,Heat and Mass Transfer7Figure1:Geometry of the IFRF furnaceC(s)+O2−→CO2(2)2C(s)+O2−→2CO(3)C(s)+CO2−→2CO(4)C(s)+H2O−→CO+H2O(5)Boudouard(Eq.(4))and gasiﬁcation(Eq.(5))reactions play an important role in MILD combustion[9,17]and they can not be neglected as usually is done for conventional coal com-bustion with atmospheric air.Char burnout is governed by the diffusion of the oxidant species from the bulk phase to the particle surface and by the heterogeneous reactions on the particle surface.Reaction rates are calculated considering global kinetic rates from[18,19].The diffusion of each chemical species from the bulk phase(∞)to the particle surface(s)is given by:β(c∞,i−c s,i)+4j=1νj,iˆR j=0(6)WhereˆR j is the rate of the reaction j,νj,i is the stoichiometric coefﬁcient of species i in reaction j andβis the mass transport coefﬁcient,calculated from the Ranz and Marshall[12] correlation,assuming unitary Lewis number.Generally,standard models[4]neglect the in-ﬂuence of the convection assuming stagnantﬂow around the particle.The diffusion of each species(Eq.(6))is equal to its production due to the heterogeneous reactions.The mass bal-ance of Eq.(6)account for the interactions between different surface reactions.In fact,CO2, produced on the particle surface from the char oxidization reaction(Eq.(2)),can react directly according to the Boudouard reaction(Eq.(4))increasing the overall char consumption.Gener-ally,standard models,such as[4]model,neglect any interaction between the different surface reactions.Further information about the model can beﬁnd in the paper of Schulze et al.[3].User Deﬁned Function(UDF)capability of Ansys-Fluent were used for coupling the SGS model,coded in C language,with the CFD solver,replacing the standard models for char com-bustion.5(a)(b)Figure 2:Comparison of temperature considering Baum and Street [4](BS)and SGS models respectively:(a)axial section contour plot;(b)radial proﬁles at 0.15,0.44,0.735,and 1.32m from the burner and comparison with experimental results [2].(a)(b)Figure 3:Comparison of CO dry volume fraction considering Baum and Street [4](BS)and SGS models respectively:(a)axial section contour plot;(b)radial proﬁles at 0.15,0.44,0.735,and 1.32m from the burner and comparison with experimental results [2].3.Validation of the SGS Char Combustion ModelValidation of SGS model was performed considering the experimental pulverized coal MILD test-case at the International Flame Research Fundation (IFRF)[2].MILD or ﬂameless com-bustion is a new technology developed for reducing pollutant emissions [8].Reactants are in-troduced at temperature generally higher than ignition temperature and the mixture is strongly diluted in order to reduce the temperature increase during reactions.The IFRF furnace is char-acterized by a square section of 2m ×2m and by a length of 6.25m,as shown in Fig.1.Primary air enters from the two lateral inlets,transporting pulverized coal particles.Secondary air is preheated by means of combustion with natural gas up to levels of 1350◦C before entering the furnace from the central inlet.Vitiated air is enriched with pure O 2in order to maintain the same concentration as atmospheric air.The furnace is ﬁred with 66kg h −1,130kg h −1and 675kg h −1respectively of coal,primary and secondary air,corresponding to a stoichiometric ratio of 1.2,as reported in Tab.1.The wall of the furnace is considered at the constant temper-Transfer7Figure4:Char consumption rate(kg/s m2)for65µm particles at0.44,0.735,1.32and2.05m from the burner.Results of Baum and Street(BS)model are reported on the left(triangles)and results of SGS on the right(circles)for each sectionsature of about1000◦C.The furnace isﬁred with Guasare coal,which proximate and ultimate analyisis are reported in Tab.2.Coal isﬁnely pulverized to give a particle size distribution with80%less than90µm[2].Particle size distribution is covered considering six classes[20]. V olatile yields are calculated by means of the CPD model,as reported in Tab.3.The single rate devolatilization model(Eq.(1))is calibrated by means of the CPD model,obtaining a pre-exponential factor of26353.9s−1and an activation energy of45.424kJ/mol.Considering recirculation of exhaust gas,the furnace is characterized by high concentrations of CO2and H2O and consequently a large fraction of char is converted through the gasiﬁcation (Eq.(4))and Boudoard(Eq.(5))reactions[17],representing an optimal test-case for validating the new char combustion model.The performances of the SGS model are therefore compared to the standard Baum and Street (BS)model andﬁnally validated against the experiments[2].Reactions Eq.(3)-(5)are consid-ered for Baum and Street model considering the same kinetic rates[18,19]used for the SGS model.Figure2(a)shows the comparison between the Baum and Street[4]and SGS models consid-ering the temperatureﬁeld.As expected,temperature gradients are very small and no clear front ofﬂame can be observed.Similar temperature proﬁles were predicted considering both models. The comparison with experiments is reported in Fig.2(b),considering four radial traverses at 0.15,0.44,0.735and1.32m from the burner.The SGS model predicts a lower temperature level in the inner jet zone,because of the increased conversion of char due to the endothermic reactions.Indeed,in this region,O2is almost completely consumed(see[9])and therefore only the endothermic gasiﬁcation(Eq.(5))and Boudouard reactions(Eq.(4))take place,absorbing heat from the gas phase.Figure3(a)shows the comparison of dry CO molar fraction on the axial section between the Baum and Street[4]and SGS models.Lower levels of CO are predicted by SGS model withrespect to Baum and Street model.Indeed,considering SGS model,the char reacting with O2 produces either CO either CO2,reducing the overall production of CO from the discrete phase. Dry CO molar fraction from numerical simulations is compared to experiments[2]considering four radial traverses at0.15,0.44,0.735and1.32m from the burner,as shown in Fig.3(b).SGS models shows a better agreement with respect to experiments.Figure4shows the char consumption rate at four cross sections for Baum and Street and SGS models considering65µm particles.As already observed for single particle simulations [3],SGS model predicts an enhanced char consumption rate with respect to Baum and Street model,nevertheless the same kinetic rates are used.In fact,the SGS model takes in account the inﬂuence of the heat and mass transport from the bulk phase to the particle surface and the interaction between the heterogeneous reaction in the particle boundary layer,enhancing the overall char consumption rate.4.ConclusionsIn this paper a new SGS model for char combustion,previously developed and validated for a single particle combustion by Schulze et al.[3],has been coupled to the commercial CFD code Ansys-Fluent and validated considering a pulverized coal MILD combustion test-case. The results have been compared to the standard Baum and Street model,used as default char combustion model by Ansys-Fluent.The comparison shows an improved prediction of the chemical species concentrations for the new SGS model with respect to the standard model. References[1]Eaton,A.et al.“Components,formulations,solutions,evaluation,and application ofcomprehensive combustion models”.In:Prog Energ Combust25.4(1999),pp.387–436.[2]Orsino,S.et al.Excess Enthalpy Combustion of Coal(Results of High Temperature AirCombustion Trials).Tech.rep.IFRF Doc.No.F46/y/3.International Flame Research Foundation,2000.[3]Schulze,S.et al.“Sub-model for a spherical char particle moving in a hot air/steamatmosphere”.In:Flow Turbul Combust(2012).(submitted).[4]Baum,M.et al.“Predicting the Combustion Behaviour of Coal Particles”.In:CombustSci Technol3.5(1971),pp.231–243.[5]Patankar,S.et al.“A calculation procedure for heat,mass and momentum transfer inthree-dimensional parabolicﬂows”.In:International Journal of Heat and Mass Transfer15.10(1972),pp.1787–1806.[6]Shih,T.et al.“A new k-epsilon eddy viscosity model for high reynolds number turbulentﬂows”.In:Computers and Fluids24.3(1995),pp.227–238.[7]Cheng,P.“Two-dimensional radiating gasﬂow by a moment method”.In:AIAA Journal2.9(1964),pp.1662–1664.[8]Cavaliere,A.et al.“Mild Combustion”.In:Prog Energ Combust30.4(2004),pp.329–366.[9]Vascellari,M.et al.“Inﬂuence of turbulence and chemical interaction on CFD pulverizedcoal MILD combustion modeling”.In:Fuel(2012).doi:10.1016/j.fuel.2011.07.042.[10]Gran I.,R.et al.“A numerical study of a bluff-body stabilized diffusionﬂame.Part1.Inﬂuence of turbulence modeling and boundary conditions”.In:Combust Sci Technol 119.1-6(1996),pp.171–190.[11]Kazakov,A.et al.Reduced Reaction Sets based on GRI-Mech1.2.http://me.berk/drm/.1994.[12]Ranz,M.et al.“Evaporation from drops:Part I”.In:Chem Eng Prog48(1952),pp.141–146.[13]F¨o rtsch,D.et al.“A kinetic model for the prediction of NO emissions from staged com-bustion of pulverized coal”.In:Proceedings of the27th Symposium(Intl.)on Combus-tion,The Combustion Institute,Pittsburgh27.2(1998),pp.3037–3044.[14]Badzioch,S.et al.“Kinetics of Thermal Decomposition of Pulverized Coal Particles”.In:Ind.Eng.Chem.Proc.Des.Dev.9.4(1970),pp.521–530.[15]Grant D.,M.et al.“Chemical model of coal devolatilization using percolation latticestatistics”.In:Energy&Fuels3.2(1989),pp.175–186.[16]Genetti,D.et al.“Development and Application of a Correlation of13C NMR Chem-ical Structural Analyses of Coal Based on Elemental Composition and V olatile Matter Content”.In:Energy&Fuels13.1(1999),pp.60–68.[17]Stadler,H.et al.“On the inﬂuence of the char gasiﬁcation reactions on NO formation inﬂameless coal combustion”.In:Combustion and Flame156.9(2009),pp.1755–1763.[18]Libby P.,A.et al.“Burning carbon particles in the presence of water vapor”.In:Com-bustion and Flame41.0(1981),pp.123–147.[19]Caram H.,S.et al.“Diffusion and Reaction in a Stagnant Boundary Layer about a CarbonParticle”.In:Industrial&Engineering Chemistry Fundamentals16.2(1977),pp.171–181.[20]Kim,J.et al.“Numerical modelling of MILD combustion for coal”.In:Progress in Com-putational Fluid Dynamics(2007).。

Numerical experiments

Deﬁnitions Importance Preliminaries
Construction. Take an (n + 1) × (n + 1) Hadamard matrix with ﬁrst row and column all +1’s, change +1’s to 0’s and −1’s to +1’s, and delete the ﬁrst row and column. Example. 1 1 1 1 1 −1 1 −1 H4 = 1 1 −1 −1 1 −1 −1 1 1 0 1 → S3 = 0 1 1 1 1 0
C. Kravvaritis Minors of (0, ±1) orthogonal matrices
Introduction A technique for minors Main Results Application to the growth problem Numerical experiments Summary-References
C. Kravvaritis
Minors of (0, ±1) orthogonal matrices
Introduction A technique for minors Main Results Application to the growth problem Numerical experiments Summary-References
Properties
1 2 3
n ≡ 3 (mod 4). SJn = Jn S = 1 2 (n + 1)Jn the inner product of every two rows and columns is 1 they are distinct, and n+ 2 , otherwise. the sum of the entries of every row and column is

Class. Quantum Grav. 16 (1999) 1817–1832. Printed in the UK PII S0264-9381(99)99800-3 Caus

Class.Quantum Grav.16(1999)1817–1832.Printed in the UK PII:S0264-9381(99)99800-3Causal set dynamics:a toy modelA Criscuolo†and H Waelbroeck‡Spinoza Instituut,Universiteit Utrecht,PO Box80195,3508TD Utrecht,The NetherlandsReceived30November1998Abstract.We construct a quantum measure on the power set of non-cyclic oriented graphs ofN points,drawing inspiration from one-dimensional directed percolation.Quantum interferencepatterns lead to properties which do not appear to have any analogue in classical percolation.Mostnotably,instead of the single phase transition of classical percolation,the quantum model displaystwo distinct crossover points.Between these two points,spacetime questions such as‘does thenetwork percolate?’have no deﬁnite or probabilistic answer.PACS numbers:0210V,11151.IntroductionThe effort to formulate a discrete theory of quantum gravity has recently recovered some of its appeal,following results on the quantization of general relativity[1],black hole thermodynamics[2]and string theory[3,15]among others.All suggest that the spectrum of excitations of a theory of quantum gravity must be discrete.It follows that counting is a natural method to deﬁne spacetime volume.A minimum framework for a discrete model of spacetime geometry brings in two key elements:number and order[4].Number gives the local conformal factor or spacetime volume element.Causal order sufﬁces to deﬁne lightcones and this represents spacetime geometry up to a conformal factor.If we assume that the causal relations between points do not form closed timelike curves,then they provide us with a structure called a partially ordered set(poset).A poset P is a discrete set with a transitive acyclic relation,namely a relation≺such that ∀x,y,z∈P,x≺y and y≺z⇒x≺z,(1.1)x≺y and y≺x⇒x=y.(1.2) For any two points x,y∈P,the Alexandrov set[6]or interval[x,y]is deﬁned by[x,y]def={z:x≺z≺y}.(1.3) A partially ordered set is said to be locallyﬁnite if the Alexandrov sets areﬁnite;it is then called a causal set.The axiom(1.2)ensures that a causal set has no closed timelike curves.†E-mail address:criscuol@phys.uu.nl‡On sabbatical leave from the Institute of Nuclear Sciences,UNAM,Circuito Exterior,CU,A Postal70-543,Mexico DF04510.E-mail address:hwael@nuclecu.unam.mx0264-9381/99/061817+16$19.50©1999IOP Publishing Ltd18171818A Criscuolo and H WaelbroeckThe idea that a notion of causality should be built into each history in the sum over histories goes back to Teitelboim[5].Several authors have proposed placing the causal set structure at the centre of a discrete formulation of quantum gravity[7,8].Poset-generating models which have been considered range from quantum spin network models[9]to stochastic models inspired from percolation theory[10].Sorkin and collaborators have conjectured that the causal set structure alone may be sufﬁcient to construct a quantum theory of spacetime [7].The poset is a discrete approximation of a physical manifold,which reproduces some topological properties of the manifold being approximated that other models are not able to reproduce[11].It has the structure of a topological space,with a topology deﬁned by the order.A family of posets ofﬁnite and increasing numbers of points,determines a family of projectiveﬁnitary topological spaces whose inductive limit is the continuous manifold being approximated[12].Moreover,the poset constitutes a genuine‘nonconmutative’space from the point of view of a generalization of the Gel’fand–Naimark theorem.In effect,to the poset corresponds a‘nonconmutative’C∗-algebra of operator-valued functions,which will be useful for constructing quantum physics on the poset[11].Independently of the particular mathematical construction that may give rise to one or another choice of partial ordering,it seems well worth the effort toﬁnd out what can be learned from the structure of causal sets per se,insofar as analysing its potential to provide a discrete representation of spacetime geometry and the use that this may have in understanding various novel approaches to quantum gravity.Progress towards what might be called a causal set representation of quantum gravity has been hindered by several factors,not least of which is the absence of a satisfactory dynamical formulation.First of all,what does one mean by‘dynamical formulation’,when the variables in question are the causal structure of spacetime itself?As often,oneﬁnds it helpful toﬁrst answer the analogous question in classical mechanics.A classical dynamical problem can be formulated as that ofﬁnding a projection operator from the set of all histories onto the subset of such histories that are solutions of the classical equations of motion.As long as there is a single classical history corresponding to any given initial data set,this formulation of the dynamical problem is equivalent to the conventional one in terms of deterministic evolution equations.In quantum mechanics it is not very meaningful to consider a single history.Instead,one would like to recast the dynamical problem in terms of subsets of the set of histories,by asking whether a subset of histories,which is determined by particular properties,is more likely to be realized than its complement.One knows that a probability cannot be assigned to sets of histories, because interference leads to violations of the probability sum rules[13].Nevertheless,a meaningful interpretation can be derived from a quantum measure on sets of histories[14],the quantum measure being a generalization of the probabilistic measure which takes into account the possibility of interference.We will adopt this point of view here,summarizing it brieﬂy in section2.In the case of causal sets,the challenge is toﬁnd a dynamical formulation which might explain how causal sets with asymptotic properties resembling those of the spacetime we live in might come to be selected as being at least reasonably likely.Markopoulou and Smolin have recently proposed a dyamical causal set model where spacelike slices are spin networks which connect to each other by means of null struts[9].However,this construction,as for any local network-building algorithm,suffers from a lack of Lorentz invariance at least at small scales,due to its reliance on horizontal slices and a local scaffolding procedure.Whether or not effective Lorentz invariance can be recovered at large scales,one would like to avoid introducing a global rest frame at the Planck scale,where the foundations of the theory are being set.Causal set dynamics:a toy model1819 What is meant by the term‘Lorentz invariance’in the present context?It is relevant only with regards to posets that can be(approximately)embedded in Minkowski space.For such posets,one considers how the causal links would look in different reference frames.A link which in one reference frame looks to be purely timelike and of small size,will in a highly boosted frame appear to be stretched out and almost null[4].A Lorentz-invariant model should not privilege any one reference frame over another,so in any given frame one should observe both short links and elongated links,in contrast to a local lattice-building model where one uses a regular lattice structure in a given frame and only allows connections between nearby points.Our purpose in this paper is to propose a simple toy model with which we will derive a quantum measure on the set of posets without introducing any a priori lattice structure.We also wish to explore what sort of questions such a dynamical causal set model should be able to answer.In order to arrive at a model that is simple enough that computer computations can be performed on relatively large posets,we choose to set aside some of the other issues that have previously frustrated attempts to construct a realistic quantum measure model for causal set dynamics.In particular,the model which we present in this paper introduces a labelling of the points by integers,and the amplitude is not required to be labelling invariant.We further simplify the problem by considering non-cyclic oriented graphs rather than posets,the difference being that transitive relations are relevant in a graph,whereas they are not relevant in the partial ordering.Both invariances,labelling and transitivity,can be recovered in the end by summing over labellings and summing over all graphs that represent the same poset.In section3we will give the outline of our toy model and present the corresponding quantum measure.A method to derive computable expressions for the measure is then described in section4,which will be applicable to sets of histories whose properties can be expressed by columns of the connectivity matrix.A few examples are evaluated numerically to reveal some of the structure of the model,including the measure of all histories with no black holes.2.Quantum measure theoryQuantum mechanics can be described as a simple generalization of classical measure theory (or probability theory).A classical measure is a map from an algebra of‘measurable sets’to the positive real numbers which satisﬁesI2(A,B)≡|A B|−|A|−|B|=0,(2.1) where denotes the union of disjoint sets.The‘no-interference’condition(2.1)permits a probabilistic interpretation for sets of histories in statistical mechanics.In quantum mechanics,the quantity I2(A,B)represents the interference term between the two sets of alternatives A,B,when interference occurs the condition(2.1)is violated and for that reason one cannot assign a probabilistic interpretation to the sum over histories formulation.Instead of(2.1),quantum theory respects a slightly weaker set of conditions, which deﬁnes a structure known as a‘quantum measure’.A quantum measure is positive real-valued function which satisﬁes the conditions|N|=0⇒|A N|=|A|,(2.2) I3(A,B,C)≡|A B C|−|A B|−|A C|−|B C|+|A|+|B|+|C|=0.(2.3)1820A Criscuolo and H WaelbroeckIt is worth noting that theﬁrst axiom(2.2),which is not necessary in probability theory because it follows from(2.1),must be included as a separate axiom for the quantum measure because I3=0in itself does not guarantee that sets with zero measure do not interfere with others.Clearly,the axiomatic structure of any theory has much to say as to how a theory should be interpreted.Since the sum over histories formulation leads to a weaker structure than(2.1),one naturally expects that quantum theory will have a weaker predictive power than probability theory insofar as its ability to discern which histories are preferred by nature.The precise nature of this weaker predictive power,and the correct interpretation of the sum over histories formulation of quantum mechanics,are encoded in the structure of the axioms(2.2)and(2.3). Sorkin has shown that this structure sustains an interpretation based on so-called‘preclusion rules’,which establish when it can be said that a certain set of histories is almost certain not to be realized in nature.As one might expect from the form of the condition I3(A,B,C)=0, these preclusion rules invoke correlations between three events,pertaining to three disjoint regions of spacetime[14].A more conventional canonical approach would demand that there exist a spacetime foliation.It is well worth recalling in this context our earlier observation that Lorentz invariance requires the existence of arbitrarily elongated links.Such links pierce through the leaves of a foliation,carrying information from past to future without crossing the present.The traditional reliance of physics on canonical quantization has also been questioned recently for other reasons,related to the teleological nature of the physics of the event horizon[16].However, the quantum measure viewpoint which we are adopting in this paper is fully consistent with canonical quantum mechanics when a global time variable can be speciﬁed.In the remainder of this paper we will limit ourselves to a yet weaker form of predictive statements than preclusion,which Sorkin refers to as‘propensity’:when referring to a particular physical property,one partitions the space of histories into two disjoint subsets by distinguishing histories which do or do not have this property.If the measure of the set of histories which do have the property is much larger than its complement it can be said to have a high propensity.The concept of propensity is useful when analysing the classical limit of a quantum theory.For example,one might distinguish spacetimes that have black holes and those that do not;if the property of having one or more black holes has a very high propensity, this would constitute a prediction of the theory in the classical limit.3.A quantum measure model for directed non-cyclic graphs3.1.Posets,causal sets and directed non-cyclic graphsAs mentioned in the introduction a poset P is a discrete set with an antisymmetric transitive relation.The transitivity rule(1.1)allows one to differentiate two types of relations:the links,which are relations that cannot be obtained from the transitive rule,and the transitive or redundant relations.The causal structure does not depend on whether a particular relation is a link or a transitive relation,so in terms of pure gravity one can say that the two types of relation are physically equivalent.However,there is a practical difference,which shows up when performing actual calculations or numerical simulations with causal sets.There is generally an enormous number of possible transitive routes between two points in a large causal set,so any algorithm which considers each possible route individually is only applicable to small causal sets,the limit being about10points.To our present knowledge,there is no generally applicable approximation scheme to perform calculations with large causal sets.Causal set dynamics:a toy model1821 The difﬁculty resides in the absence of a convenient(one-to-one)representation of posets. The most natural approach would be to represent a poset in terms of its‘relations matrix’, where R ij=1if and only if x j x i and otherwise R ij=0(x i,x j∈P;i,j∈N).Thematrices R satisfy the transitivity conditionR ij=θk R ik R kj,(3.1)whereθ(x)=1if x>0and zero otherwise.The computational complexity of checkingthese equations for each of the N2possible binary matrices grows like N5(using the moststraightforward algorithm).One could equally well choose to use the link matrix,L,whereL ij=1if and only if x i x j is a link,but of course this leads to the same computational problem.The limitation on the number of points with which one can work affects almost everyposet-related calculation.For example,the counting of posets with a given number of pointshas only been solved up to N=11[18].For large values of N it has been shown to growasymptotically like[17]C×2N2/4+3N/2e N N−(N+1),but the method that yields this result does not readily generalize to other poset calculations.To avoid these problems,which originate from the transitivity condition,we will considerthe set of all lower-triangular binary matrices,regardless of whether or not they include allpossible transitive relations.We will refer to these matrices as connectivity matrices,anddenote them by C.A connectivity matrix represents a directed non-cyclic graph,i.e.a set of points connectedby arrows such that arrows do not form closed loops.Given such a graph,it is always possibleto label the points with consecutive integers in such a way that arrows point from a lesser labelto a greater one.One then arrives at a lower-triangular binary matrix,where each entry C ij isequal to one if and only if there is an arrow in the graph from j to i(with j<i),and zerootherwise.Some of the connections represented in the matrix C may be links,while otherswill be transitive relations,but unlike the matrix R it is not required that all transitive relationsbe represented as connections or arrows of the graph.We will understand the connectivity matrix to represent a‘history’or a possibility for‘spacetime’.The sum over histories then takes the form of an unconstrained sum over binaryarrays,which makes it relatively easier to apply standard analytical tools.Of course in the end the purpose is not to compute the sum over all histories,but overspeciﬁed subsets of histories,which satisfy one or another physical property of interest.A‘physical property’is a property of the causal structure,i.e.one that does not pertain either tothe labelling of the points or to transitive relations.The sum over graphs will then include thesum over all the causal orders with the given property.There are interesting physical properties that condition the connectivity matrix withoutincreasing the computational complexity to the same extent as condition(3.1).For suchproperties,a graph-based model can be expected to yield computable expressions.We willsee some examples in section4.3.2.The‘ﬁnal question’A quantum measure model can be constructed from two basic elements:an amplitude functionon the set of histories,and a‘ﬁnal question’,Q f,which by deﬁnition must refer only to thepart of the histories in the causal future of the region of interest.The question Q f must bewell posed,in the sense that each history,γ,in the Hilbert space,H,should give one and only1822A Criscuolo and H Waelbroeckone answer.The answer set then gives a partition of the Hilbert space in a disjoint union of sets E i ,such thatγ∈E i ⇔Q f (γ)=a i ,(3.2)where a i ,i =1,...,n (n 2),represents an element of the answer set.Given an amplitude a :H →C ,one can then construct the following function |·|on the power set of H (when the computational procedure at hand does not give a ﬁnite answer for every subset of histories,one requires that the quantum measure be well deﬁned on a sigma-algebra of ‘measurable sets’):|A |= i γ∈A ∩E ia(γ) 2.(3.3)One easily veriﬁes that |·|satisﬁes the axioms (2.2)and (2.3)of the quantum measure.We have introduced a ‘ﬁnal question’as part of the procedure to construct a quantum measure,but this question in itself is not of any particular interest.Eventually,one would like to be able to show that the dependence of the quantum measure on the ﬁnal question vanishes asymptotically in the limit of very large posets.In other words,the ﬁnal question would then be an artiﬁce introduced for the sole purpose of performing the calculation.The aim of the quantum measure formalism is to address other questions,which one might call ‘physical’questions.These should be formulated in such a way that they refer to the history before the ﬁnal conditions.Physical questions can be either about the ‘present’state of the system,with an appropriate deﬁnition of ‘present state’,or about the spacetime history.For each possible answer there is a set A i of histories with answer a i ,and the relative values of the quantum measures for different possible answers will reveal what the model has to say regarding this question.The quantum measure formalism can sometimes contain a canonically quantized model,in the following sense.One ﬁrst deﬁnes a one-parameter family of questions,which are to form a complete set,in the sense that their combined answers describe a history completely.These questions can be stated as ‘what is the state of the system at time t ?’The answer set is given by the eigenstates of a complete set of commuting (conﬁguration-space)observables.When the sets of histories corresponding to different answers E t j do not interfere,and the sum of the measures over all possible answers at time t is equal to one,then the quantum measureformalism will provide the same information as a canonical theory for that particular family ofquestions.This will occur when the unitarity condition I 2(E t j ,E t k )=2δjk is satisﬁed,where I 2(E t j,E t k )= i γ1∈E t j ∩E i γ2∈E t k ∩E i(a(γ1)a ∗(γ2)+a ∗(γ1)a(γ2)) .(3.4)The choice of a particular one-parameter family of questions is analogous to a choice of ‘slicing’in canonical quantization.Other possible slicings,based on different choices of one-parameter families of questions,may well lead to different unitarity conditions.This sort of exercise of course is only relevant to the extent that one is interested in making contact with canonical quantization methods.Taking a pre-geometrical perspective,one might argue that other criteria should take precedence over unitarity in guiding the search for the correct quantum measure;if in the end it turned out that such criteria were to lead one to a unique quantum measure,then the reverse problem of ﬁnding the choice(s)of slicing for which a unitary canonical theory can be deduced would provide a satisfactory solution to the problem of time.Here we will choose the following ‘ﬁnal question’:‘which of the points i <N emit an arrow towards N ?’By deﬁnition of the term ‘ﬁnal question’we are assuming that the predictive power of this model is limited to points that do not lie to the future of N .We willCausal set dynamics:a toy model 1823label the points in such a way that N is the largest label and there are N −1other points which may or may not be to the past of N but will certainly not be to its future.The answer set of this question is then the set of binary words of length N −1,where a bit is equal to 0in the absence of an arrow and 1denotes the presence of a connection.In terms of the connection matrix,the answer to the ﬁnal question will be given by its last row, CN .This ﬁnal question generates a partition of the space of histories into disjoint subsets of connectivity matrices with a ﬁxed lowest row.There is a natural one-parameter family of questions associated with this particular choice of ﬁnal question,namely those whose answers are the rows of the connectivity matrix.These are the questions:‘which points emit arrows towards the point labelled by t ?’.Note that in this case the cardinality of the answer set grows like 2t −1.3.3.The amplitude and the quantum measureWe will make the following ansatz for the amplitude (factorizability):a(C )def =N m =2a( C m −1,m −1→ C m ,m),(3.5)wherea( C m −1,m −1→ C m ,m)def =A C mm −1l<m −1A C ml C m −1l .(3.6)We will also require that the amplitude to create a connection C m be independent of the previous state, C m −1.Choosing A 0=A 00=i A 01=√q,(3.7)A 1=A 11=i A 10=√p(3.8)and p +q =1(p,q ∈R >0),we arrive at a quantum generalization of a one-dimensional directed percolation model.In the one-dimensional directed percolation model,N points are labelled by consecutive integers as in our model,and each point can connect to any of the previous points with a constant probability p .The propagator (3.6)then becomes a( C m −1,m −1→ C m ,m)=(√q)m −1−l C ml (√p) l C ml (−i ) i (1−δC mi C m −1i ).(3.9)Using (3.4)and substituting (3.5),one ﬁnds that the model would be unitary with the slicing { Cm ;m =1,2,...}if one choosesC m a( C m −1,m −1→ C m ,m)a ∗( C m −1,m −1→ C m ,m)=m −1 l =1δC m −1l C m −1l,i.e.if p =q =12.To simplify the notation,we will introduce the total number of entries equal to one in them th column of the connectivity matrix,C m =Ni =1C im ,(3.10)and the number of ‘kinks’in each column,K m =N −1 i =m +1(1−δC im C i +1m ).(3.11)1824A Criscuolo and H WaelbroeckThe sum of these quantities over all of the columns of the connectivity matrix yields the total connectivity C and total kink number K ,respectively.We then arrive at a simple expression for the amplitude of a connectivity matrix (or ‘history’),a(C )=(√p)C (√q)C N 2−C (−i )K ,(3.12)where C N 2is the binomial coefﬁcient.The quantum measure of a set A of connectivity matrices is then given |A |=C N |ψ(A, C N ,N)|2,(3.13)whereψ(A, C N ,N)def =C ∈A : C N ﬁxed a(C ).(3.14)4.Analytical and numerical results4.1.Measure of the space of histories,|H |As a ﬁrst example one can compute the measure of the space H of all possible non-cyclic oriented graphs of N points.In that case,one can drop the label A in (3.14)and write |H |=CN |ψ( C N ,N)|2.(4.1)The argument of the square modulus can be loosely interpreted as a ‘cosmological wavefunction’.To compute ψ( C N ,N)def = C : C N ﬁxed (√p)C (√q)C N 2−C (−i )K ,(4.2)we use the fact that C = m C m and K = m K m to factorize the above expression andconsider each column of the connectivity matrix individually:ψ( C N ,N)=(√q)C N 2N −1m =1ψm (C Nm ,N),(4.3)ψm (C Nm ,N)def =N −m +C Nm −1 C m =C NmK m (C m ) p q C m(−i )K m N C Nm (K m ,C m ),(4.4)where N C Nm (K m ,C m )is the number of binary words of N −m −1bits with C m bits equal to 1and K m kinks,when the last bit in the column has been set equal to C Nm .The functions N C Nm (K m ,C m )can be calculated by considering the number of ways of making k cuts in a sequence of C m ones and inserting the zeros at the cuts.One ﬁnds N 0(2k,C m )=δ(k,0)δ(C m ,0)+ C m −1k −1 N −m −1−C m k ,(4.5)N 1(2k,C m )=δ(k,0)δ(C m ,N −m)+ C m −1k N −m −1−C m k −1 ,(4.6)N 0(2k +1,C m )= C m −1k N −m −1−C m k ,(4.7)N 1(2k +1,C m )= C m −1k N −m −1−C m k.(4.8)Causal set dynamics:a toy model1825 To determine the bounds on K m,we must regard its dependence on C Nm and make a comparison between the number of bits equal to1,C m,and the number of bits equal to zero, N−m−C m,analysing the different possible arrangements.The result isif N−m−C m<C m N−m+C Nm−1 ⇒K m∈[1−C Nm,2(N−m−C m+C Nm−1)+1−C Nm],if N−m−C m=C m ⇒K m∈[1,2C m−1],if C Nm C m<N−m−C m ⇒K m∈[C Nm,2C m−C Nm].These equations allow one to calculate the functionψ( C N,N)numerically and compute the measure|H|.Not surprisingly,this measure is equal to1in the‘unitary case’p=q=1/2 whenψcan be interpreted as a cosmological wavefunction.This case is not particularly relevant in the context of the quantum measure interpretation.4.2.Graphs with aﬁxed number of arrows from a given pointSetting aside for the time being the issue of labelling invariance,we will compute the propensity that a point with a given label emit aﬁxed number of arrows towards other points of the graph.Let A Cl (l)be the set of histories(graphs)where point labelled l emits C l outgoing arrows.From(3.12)–(3.14)we have|A Cl (l)|=CN|ψ(A Cl(l), C N,N)|2,(4.9)whereψ(A Cl (l), C N,N)=(√q)C N2pqCNN−1ψl(A Cl(l),C Nl,N)N−2m=1m=lψm(C Nm,N)(4.10)andψl(A Cl (l),C Nl,N)def=Vl:A Cl(l)pqCl(−i)K l,(4.11)where we have used the short-hand notation V l:A C l(l)to denote the sum runs over the columns of binary words V l,which satisfy that the number of connections C l beﬁxed.Oneﬁnds|A Cl(l)| |H|=pqCl|Z0|2+|Z1|2|ψl(0,N)|2+|ψl(1,N)|2,(4.12)whereZ CNl =K l(−i)K l N CNl(K l,C l)(4.13)represents the quantum interference factor.This ratio was computed numerically for p=q= 1,as a function of Cl.In the classical directed percolation model,the probability of C l is given by a binomialdistributionp(C l)=N−l−1C lp C l(1−p)N−l−1−C l.(4.14)The ratio of the quantum measure to its classical counterpart,F(l)=|A Cl(l)|p(C l)(4.15)1826A Criscuolo and H WaelbroeckFigure1.The propensity that a point(point number90from theﬁnal point)emits C arrows isrepresented as a function of C,for the case p=q=12.The sharp rise and fall at both ends areremnants of a classical binomial distribution,whilst in the middle quantum interference is observed.In contrast to the classical case the propensity does not peak at C=45,half the possible numberof arrows.Figure2.The real part of the form factor which gave rise to the previousﬁgure.is a form factor which can be interpreted as representing the effect of quantum interference.A comparison between(ﬁgure1)and the binomial distribution reveals that the destructive interference is most important at the midpoint where half the available bits are equal to one. When C l is nearer one of the two extremal values one observes the same exponential dropoff as in the statistical model.In the crossover between these two regimes,a striking interference pattern is observed. To witness the origin of this interference we show the real part of the form factor Z0(the imaginary part is its mirror image)(ﬁgure2).4.3.Graphs with no arrows from a given pointWhen the number of out-arrows from a point is equal to zero this point is a sink in the oriented graph.It is then not emitting any outgoing signals and therefore might loosely be called a。

Two-Scale Continuum Model for Simulation of Wormholes in Carbonate Acidization

Two-Scale Continuum Model for Simulation of Wormholes in Carbonate AcidizationMohan K.R.Panga and Murtaza ZiauddinSchlumberger,Sugar Land,TX77478Vemuri BalakotaiahDept.of Chemical Engineering,University of Houston,Houston,TX77204DOI10.1002/aic.10574Published online September6,2005in Wiley InterScience().A two-scale continuum model is developed to describe transport and reaction mecha-nisms in reactive dissolution of a porous medium,and used to study wormhole formationduring acid stimulation of carbonate cores.The model accounts for pore level physics bycoupling local pore-scale phenomena to macroscopic variables(Darcy velocity,pressureand reactant cup-mixing concentration)through structure-property relationships(perme-ability-porosity,average pore size-porosity,and so on),and the dependence of masstransfer and dispersion coefﬁcients on evolving pore scale variables(average pore sizeand local Reynolds and Schmidt numbers).The gradients in concentration at the porelevel caused byﬂow,species diffusion and chemical reaction are described using twoconcentration variables and a local mass-transfer coefﬁcient.Numerical simulations ofthe model on a two-dimensional(2-D)domain show that the model captures the differenttypes of dissolution patterns observed in the experiments.A qualitative criterion forwormhole formation is developed and it is given by⌳ϳO(1),where⌳ϭ͌k eff D eT/u o.Here,k eff is the effective volumetric dissolution rate constant,D eT is the transversedispersion coefﬁcient,and u o is the injection velocity.The model is used to examine theinﬂuence of the level of dispersion,the heterogeneities present in the core,reactionkinetics and mass transfer on wormhole formation.The model predictions are favorablycompared to laboratory data.©2005American Institute of Chemical Engineers AIChE J,51:3231–3248,2005Keywords:wormholes,carbonate acidizing,reactive dissolution,porous mediaIntroductionAcid treatment of carbonate reservoirs is a widely practiced oil and gas well stimulation technique.The primary objective of this process is to increase the production rate by increasing permeability of the damaged zone near the wellbore region. The injected acid dissolves the material near the wellbore and createsﬂow channels that establish a good connectivity be-tween the reservoir and the well.While dissolution increases permeability,the relative increase in permeability for a given amount of acid is observed to be a strong function of the injection conditions.At very low injection rates,acid is spent soon after it contacts the medium resulting in face dissolution. The penetration depth of the acid is restricted to a region very close to the wellbore.On the other hand,at very high injection rates,acid penetrates deep into the formation but the increase in permeability is not large because the acid reacts over a large region leading to uniform dissolution.At intermediateﬂow rates,long conductive channels known as wormholes are formed.These channels penetrate deep into the formation andCorrespondence concerning this article should be addressed to V.Balakotaiah atbala@.©2005American Institute of Chemical EngineersAIChE Journal3231December2005Vol.51,No.12facilitate the ﬂow of oil.Thus,for successful stimulation of a well it is required to produce wormholes with optimum density and penetrating deep into the formation.A detailed description of ﬁeld practices for carbonate acidizing can be found in the literature.1–4Several experimental studies have been conducted in the past to understand wormhole formation and to predict the condi-tions required for creating wormholes.5–14In those experi-ments,acid was injected into a core at different injection rates and the volume of acid required to break through the core,also known as breakthrough volume,was measured for each injec-tion rate.A common observation in experimental studies is that dissolution creates patterns that are dependent on the injection rate.These dissolution patterns were broadly classiﬁed into three types:uniform,wormholing and face dissolution patterns corresponding to high,intermediate and low injection rates,respectively.Figure 1shows typical dissolution patterns ob-served in experiments 15on carbonate cores treated with HCl at different injection rates.It is also observed that wormholes form at an optimum injection rate,and because only a selective portion of the core is dissolved the volume required to stimu-late the core is minimized (see Figure 2).Furthermore,the optimal conditions for wormhole formation were observed to depend on various factors such as acid/mineral reaction kinet-ics,diffusion coefﬁcients of the acid species,concentration of acid,temperature,geometry of the system (linear/radial ﬂow),and so on.Many theoretical studies have been conducted in the past to understand the phenomena of ﬂow channeling associ-ated with reactive dissolution in porous media,and to obtain an estimate of the optimum injection rate.However,the models developed thus far describe only a few aspects of the acidiza-tion process,and the coupling between reaction and transport mechanisms that plays a key role in reactive dissolution is not completely accounted for in these models.In this work,we present a two-scale continuum model that describes reaction and transport in a porous medium retaining the essential phys-ics of dissolution.The model is used to investigate the inﬂu-ence of different factors on wormhole formation.Literature ReviewWormhole formation during reactive dissolution of carbon-ates is a subject that has been actively studied in the last thirty years.To explain wormhole formation,numerous models,ranging from detailed pore-scale models (for example,network models)that account for reaction,transport and dissolution at the pore scale to single wormhole (tube)models that consider only the mechanisms occurring inside the wormholes,have been developed in the literature.In this section,a brief review of different models developed to study wormhole formation in carbonates is presented.Relating the important dimensionless groups of the system to experimental observations is one of the approaches followed to model acidization process.For example,Fredd and Fogler 12,15and Fredd 16have reported the dependence of wormhole for-mation on the Damko ¨hler number and predicted an optimum Damko ¨hler number of 0.29for different ﬂuid/mineral systems.Daccord et al.7,8,17also used a dimensionless parameter based approach and coupled it with the concept of fractals to obtain the propagation rate of wormholes.The use of a few dimen-sionless groups to explain experimental observations is a dif-ﬁcult exercise because the actual number of parameters in the system is large and it is difﬁcult to study wormholing phenom-ena systematically using thisapproach.Figure 2.Typical breakthrough curveobserved inacidization experiments.The pore volume of acid (HCl,22°C)required for break-through is high at very low and very high injection rates,and is minimum at optimum injection rate,Q opt .The length and diameter of the cores are 10.2cm and 3.8cm,respectively.The initial porosities and permeabilities are in the range of 0.15–0.2and 0.8–2md,respectively.Figure 1.Typical dissolution patterns observed in car-bonate acidizing:(a)face dissolution,Q ؍0.04cc/min,PV inj ؍43.1,(b)conical Q ؍0.11cc/min,PV inj ؍10,(c)wormhole Q ؍1.05cc/min,PV inj ؍0.8,(d)ramiﬁed Q ؍10cc/min,PV inj ؍2.1,and (e)uniform Q ؍60cc/min,PV inj ؍6.7.Hydrochloric acid is used in these experiments and the acid injection rate is increased from (a)to (e)(Fredd and Fogler 15).The cores are approximately 3.8cm in dia.and 10.2cm in length.3232AIChE JournalDecember 2005Vol.51,No.12In single wormhole models,a cylindrical tube represents the wormhole and the effects ofﬂuid leakage,reaction kinetics, and so on,on wormhole propagation rate are investigated using these models.10,14,18,19,20,21One of the key results from the studies on single wormhole models is on the interaction of wormholes and competition between them.It has been report-ed14,20that the growth rates of multiple wormholes in a domain is dependent on the separation distance between them.Panga et al.22used a continuum model and made a similar observation in their numerical simulation of wormhole density.The simple structure of single wormhole models offers the advantage of studying reaction,diffusion and convection mechanisms inside the wormhole in detail;however,these models cannot be used to study wormhole initiation,the other dissolution patterns (face,ramiﬁed,and so on)observed in the experiments,and the effect of heterogeneities on wormhole formation.Hoefner and Fogler,9Fredd and Fogler,12and Daccord et al.17have developed network models to describe reactive work models represent the porous medium as a network of tubes interconnected to each other at the nodes.The ﬂow of acid inside these tubes is described using the Hagen-Poiseuille relationship for laminarﬂow.The acid reacts at thewall of the tube and dissolution is accounted for in terms of an increase in the tube work models predict the disso-lution patterns and qualitative features of dissolution like the optimumﬂow rate observed in experiments.However,a core scale simulation of the network model is computationally very expensive.An intermediate approach to describing reactive dissolu-tion involves the use of averaged or continuum models. Unlike the network models that describe dissolution from the pore scale,and the models based on the assumption of existing wormholes,the averaged models describe dissolu-tion at the Darcy scale.The Darcy scale model requires information on the pore scale processes,which are obtained from a pore scale model.The predictions of the pore scale model depend on the pore structure that changes with time due to dissolution.Obtaining detailed pore structure of a core and approximating its change during dissolution is very difﬁcult and is one of the disadvantages of using a Darcy scale model.However,by using different pore scale models that are representative of the core,the sensitivity of the results obtained from Darcy scale models can be studied. Averaged models for carbonate acidization have been de-veloped by Liu et al.23Chen et al.24,25and Golﬁer et al.26 These models were shown to capture qualitative and some quantitative features of dissolution.The model developed by Liu et al.and Chen et al.does not consider the effect of mass transfer on the reaction rate and is valid only in the kinetic regime,while the model developed by Golﬁer et al.is valid only in the mass-transfer controlled regime(to be deﬁned later).This work presents a continuum model that captures both the extremes of reaction(kinetic and mass transfer controlled)simultaneously by using two concentration vari-ables and a mass-transfer coefﬁcient.This allows the de-scription of a wide range of acids as demonstrated later.The model proposed here is similar to the widely used two-phase models of catalytic reactors,27the main difference being coupling of theﬂow,dissolution/reaction and pore scale mass-transfer processes.Model for Carbonate DissolutionReaction between a carbonate porous medium and acid leads to complete dissolution of the medium,thereby in-creasing the permeability to a large value.At very low injection rates in a homogeneous medium,this would give rise to a planar reaction/dissolution front where the medium behind the front is completely dissolved and the medium ahead of the front is not dissolved.However,the presence of natural heterogeneities in the medium leads to an uneven increase in permeability along the front leading to regions of high and low permeabilities.The high permeability regions attract more acid,which further dissolves the medium cre-ating channels that travel ahead of the front.Thus,adverse mobility(ϭK/␮,where K is the permeability and␮is the viscosity of theﬂuid)arising because of differences in the permeabilities of the dissolved and undissolved medium and heterogeneity are required for channel formation.This re-action-driven instability has been studied using linear and weakly nonlinear stability analyses by some au-thors.28,29,30,31This instability is similar to the viscousﬁn-gering instability where adverse mobility arises due to a difference in viscosities of the displacing and displaced ﬂuids.32The shape(wormhole,conical,and so on)of the channels is, however,dependent on the relative magnitudes of convection and dispersion in the medium.For example,when transverse dispersion is more dominant than convective transport,reaction leads to conical and face dissolution patterns.On the other hand,when convective transport is more dominant,the con-centration of acid is more uniform in the domain leading to a uniform dissolution pattern.The model presented here de-scribes the phenomena of reactive dissolution as a coupling between processes occurring at two scales,namely the Darcy scale and the pore scale.Different length scales are shown in Figure3.In the following subsections,these two parts of the model arediscussed.Figure3.Different length scales used in the model.AIChE Journal3233December2005Vol.51,No.12Darcy scale modelThe Darcy scale model equations are given byUϭϪ1␮K⅐ٌP(1)Ѩ␧ѨtЈϩٌ⅐Uϭ0(2)Ѩ͑␧C f͒ѨtЈϩٌ⅐͑U C f͒ϭٌ⅐͑␧D e⅐ٌC f͒Ϫk c a␯͑C fϪC s͒(3)k c͑C fϪC s͒ϭR͑C s͒(4)Ѩ␧ѨtЈϭR͑C s͒a␯␣␳s(5)Here Uϭ(U,V,W)is the Darcy velocity vector,K is the permeability tensor,P is the pressure,␧is the porosity,C f is the cup-mixing concentration of the acid in theﬂuid phase,C s is the concentration of the acid at theﬂuid-solid interface,D e is the effective dispersion tensor,k c is the local mass-transfer coefﬁcient,a␯is the interfacial area available for reaction per unit volume of the medium,␳s is the density of the solid phase and␣is the dissolving power of the acid,deﬁned as grams of solid dissolved per mole of acid reacted.The reaction kinetics is represented by R(C s).For aﬁrst-order reaction R(C s)re-duces to k s C s where k s is the surface reaction rate constant (having the units of velocity).Equation3gives a Darcy scale description of the transport of the acid species.Theﬁrst three terms in the equation represent accumulation,convection and dispersion of the acid,respec-tively.The fourth term describes the transfer of acid species from theﬂuid phase to theﬂuid-solid interface and its role is discussed later in this section.The velocityﬁeld U in the convection term is obtained from Darcy’s law(Eq.1)relating velocity to the permeabilityﬁeld K and the gradient of pres-sure.Theﬁrst term in the continuity Eq.2accounts for the effect of local volume change during dissolution on theﬂow ﬁeld.While deriving the continuity equation,it is assumed that the dissolution process does not change the density of theﬂuid phase signiﬁcantly.The transfer term in the species balance Eq.3describes the depletion of the reactant at the Darcy scale due to reaction.An accurate estimation of this term depends on the description of transport and reaction mechanisms inside the pores.In the absence of reaction,concentration of the acid species is uni-form inside the pores.Reaction at the solid-ﬂuid interface gives rise to concentration gradients in theﬂuid phase inside the pores.The magnitude of these gradients depends on the relative rate of mass transfer from theﬂuid phase to theﬂuid-solid interface and reaction at the interface.If the reaction rate is very slow compared to the mass-transfer rate,the concentration gradients are negligible.In this case,the reaction is considered to be in the kinetically controlled regime,and a single concen-tration variable is sufﬁcient to describe this situation.However, if the reaction is very fast compared to the mass transfer,steep gradients develop inside the pores.This regime of reaction is known as the mass-transfer controlled regime.To account for the gradients developed due to mass transfer control requires the solution of a differential equation describing diffusion and reaction mechanisms inside each of the pores.Since this is impractical,we use two concentration variables C s and C f,for the concentration of the acid at theﬂuid-solid interface and in theﬂuid phase respectively,and capture the information con-tained in the local concentration gradients as a difference between the two variables using the concept of a local mass-transfer coefﬁcient(Eq.4).Equation4balances the amount of reactant transferred to the surface to the amount reacted.For the case ofﬁrst order kinetics,R(C s)ϭk s C s,Eq.4can be simpliﬁed toC sϭC fͩ1ϩs k cͪ(6)In the kinetically controlled regime(k sϽϽk c),the concentra-tion at theﬂuid-solid interface is approximately equal to the concentration of theﬂuid phase(C sϷC f).In the mass transfer controlled regime(k sϾϾk c),the value of concentration at the ﬂuid-solid interface is very small(C sϷ0).Since the rate constant isﬁxed for a given acid,the magnitude of the ratio k s/k c is determined by the local mass-transfer coefﬁcient k c, which is a function of the pore geometry,the reaction rate and the local hydrodynamics.Due to dissolution and heterogeneity in the medium,the ratio k s/k c is not a constant in the medium but varies with space and time which can lead to a situation where different locations in the medium experience different regimes of reaction.To describe such a situation it is essential to account for both kinetic and mass transfer controlled regimes in the model.Equation5describes the evolution of porosity in the domain due to reaction.To complete the model(Eqs.1–5),information on the per-meability tensor K,dispersion tensor D e,mass-transfer coefﬁ-cient k c,and interfacial area a␯is required.These quantities depend on the pore structure and are inputs to the Darcy scale model from the pore scale model.Instead of calculating these quantities from a detailed pore scale model taking into consid-eration the actual pore structure,we use structure-property relations that relate permeability,interfacial area and average pore radius of the pore scale model to its porosity.However,if a detailed calculation including the pore structure can be made, then the quantities K,D e,k c,and a␯obtained from such a calculation can be used as inputs from the pore scale model to the Darcy scale model.Pore scale model(a)Structure-Property Relations.Dissolution changes the structure of the porous medium continuously,thus,making it difﬁcult to correlate the changes in local permeability to po-rosity during acidization.The results obtained from averaged models,which use these correlations,are subject to quantitative errors arising from the use of a bad correlation between the structure and property of the medium,although the qualitative trends predicted may be correct.Since a deﬁnitive way of relating the change in the properties of the medium to the change in structure during dissolution does not exist,we use3234AIChE JournalDecember2005Vol.51,No.12semiempirical relations that relate the properties to local po-rosity.The relative increase in permeability,pore radius andinterfacial area with respect to their initial values are related to porosity in the following mannerK K o ϭ␧␧o ͩ␧͑1Ϫ␧o ͒␧o ͑1Ϫ␧͒ͪ2␤(7)r pr oϭͱoK o ␧(8)anda ␯a o ϭ␧r o␧o r p(9)Here K o ,r o and a o are the initial values of permeability,average pore radius and interfacial area,respectively.Notice that,for ␤ϭ1the structure-property relation for permeability evolution reduces to the well-known Carman-Kozeny correla-tion K ϰ␧3/(1Ϫ␧)2relating permeability to the porosity of the medium.The parameter ␤is introduced into Eq.7to extend the relation to a dissolving medium.Figure 4shows a typical plot of permeability versus porosity for different values of the parameter ␤.In addition,the effect of structure-property rela-tions on breakthrough time has also been tested by using different correlations (in a later section).The model yields better results if structure-property correlations that are devel-oped for a particular system of interest are used.Note that,in the earlier relations permeability,which is a tensor,is reduced to a scalar for the pore scale model.In the case of anisotropic permeability,extra relations for the permeability of the pore scale model are needed to complete the model.(b)Mass-Transfer Coefﬁcient.Transport of the acid spe-cies from the ﬂuid phase to the ﬂuid-solid interface inside the pores is quantiﬁed by the mass-transfer coefﬁcient (k c ).It plays an important role in characterizing dissolution phenom-ena because the mass-transfer coefﬁcient determines the re-gime of reaction for a given acid (Eq.6).The mass-transfer coefﬁcient depends on the local pore structure,reaction rate and local velocity of the ﬂuid.Gupta and Balakotaiah 33and Balakotaiah and West 34investigate the contribution of each of these factors to the local mass-transfer coefﬁcient in detail.For developing ﬂow inside a straight pore of arbitrary cross section,it is shown 34that a good approximation to the Sherwood number,the dimensionless mass-transfer coefﬁcient,is given bySh ϭ2k c r p m ϭSh ϱϩ0.35ͩd hͪ0.5Re p 1/2Sc 1/3(10)where k c is the mass-transfer coefﬁcient,r p is the pore radius and D m is molecular diffusivity,Sh ϱis the asymptotic Sher-wood number for the pore,Re p is the pore Reynolds number,d h is the pore hydraulic diameter,x is the distance from the pore inlet and Sc is the Schmidt number (Sc ϭ␯/D m ;where ␯is the kinematic viscosity of the ﬂuid).Assuming that the length of a pore is typically a few pore diameters,the average mass-transfer coefﬁcient can be obtained by integrating the above expression over a pore length,which givesSh ϭSh ϱϩb Re p1/2Sc 1/3(11)where the constants Sh ϱand b (ϭ0.7/m 1/2,m ϭpore length to diameter ratio)depend on the structure of the porous medium (pore cross sectional shape and pore length to hydraulic diam-eter ratio).Equation 11is of the same general form as the Frossling correlation used extensively in correlating mass-transfer coefﬁcients in packed-beds.(For a packed bed of spheres,Sh ϱϭ2and b ϭ0.6.This value of b is close to the theoretical value of 0.7predicted by Eq.11for m ϭ1.)The two terms on the righthand side of Eq.11are contribu-tions to the Sherwood number due to diffusion and convection of the acid species,respectively.While the diffusive part,Sh ϱdepends on the pore geometry,the convective part is a function of the local velocity.The asymptotic Sherwood number for pores with cross-sectional shape of square,triangle and circle are 2.98,2.50and 3.66,respectively.34Since the value of asymptotic Sherwood number is a weak function of the pore cross-sectional geometry,we use a typical value of 3.0in our calculations.The convective part depends on the pore Reyn-olds number and the Schmidt number.For liquids,the typical value of Schmidt number is around one thousand and assuming a value of 0.7for b ,the approximate magnitude of the con-vective part of the Sherwood number from Eq.11is 7Re p 1/2.The pore Reynolds numbers are very small due to the small pore radius,and the low injection velocities of the acid,making the contribution of the convective part negligible during the initial stages of dissolution.As dissolution proceeds,the pore radius and the local velocity increase,making the convective contribution signiﬁcant.(Remark:For typical values citedabove,the two contributions are equal when 7Re p1/2Ϸ3or Re p Ϸ0.2.)However,the effect of this convective mass transfer on acid consumption may not be signiﬁcant because of the extremely low interfacial area in the high porosity regions where convection is dominant.The effect of reaction rate on the mass-transfer coefﬁcient is observed to be weak.34For example,the asymptotic Sherwood number varies from48/11Figure 4.Variation of permeability with porosity for dif-ferent values of ␤.AIChE Journal3235December 2005Vol.51,No.12(ϭ4.36)to3.66in the case of a very slow reaction to a very fastreaction.(c)Fluid Phase Dispersion Coefﬁcients.For homoge-neous,isotropic porous media,the dispersion tensor is charac-terized by two independent components,namely,the longitu-dinal,D eX and transverse,D eT,dispersion coefﬁcients.In theabsence ofﬂow,dispersion of a solute occurs only due tomolecular diffusion and D eXϭD eTϭ␣os D m,where D m is the molecular diffusion coefﬁcient and␣os is a constant that depends on the structure of the porous medium(for example,tortuosity or connectivity between the pores).Withﬂow,thedispersion tensor depends on the morphology of the porousmedium,as well as the pore levelﬂow andﬂuid properties.Therelative importance of convective to diffusive transport at thepore level is characterized by the Peclet number in the pore,deﬁned byPe pϭ͉U͉d h␧D m(12)where͉U͉is the magnitude of the Darcy velocity,and d h is the pore hydraulic diameter.For a well-connected pore network, random walk models and analogy with packed beds may be used to show thatD XϭD eXD mϭ␣osϩ␭X Pe p(13)andD TϭD eTD mϭ␣osϩ␭T Pe p(14)where␭X and␭T are numerical coefﬁcients that depend on the structure of the medium(␭XϷ0.5,␭TϷ0.1for packed-bed of spheres).Other correlations used for D eX,are of the formD eXm ϭ␣osϩ1Pe p lnͩ3Pe pͪ(15)andD eXD mϭ␣osϩ␭Pe p2.(16)Equation16is based on Taylor-Aris theory,and it is normally used when the connectivity between the pores is very low. These as well as the other correlations in the literature predict that both the longitudinal and transverse dispersion coefﬁcients increase with the Peclet number.Table1shows typical values of pore Peclet numbers calcu-lated,based on the core experiments(permeability of the cores is approximately1mD)of Fredd and Fogler.12The injection velocities of the acid(0.5M HCl)are varied between0.14cm/s and1.4ϫ10Ϫ4cm/s,where0.14cm/s corresponds to the uniform dissolution regime and1.4ϫ10Ϫ4cm/s corresponds to the face dissolution regime.The values of pore diameter, molecular diffusion and porosity used in the calculations are 0.1␮m,2ϫ10Ϫ5cm2/s and0.2,respectively.It appears from the low values of pore level Peclet number in the face disso-lution regime that dispersion in this regime is primarily due to molecular diffusion.The Peclet number is close to order unity in the uniform dissolution regime showing that both molecular and convective contributions are of equal order.In the numer-ical simulations presented later,it is observed that the disper-sion term in Eq.3does not play a signiﬁcant role at high injection rates(uniform dissolution regime)where convection is the dominant mechanism.As a result,the form of the convective part of the dispersion(␭X Pe p,Pe p ln(3Pe/2),and so on),which becomes important in the uniform dissolution regime,may not affect the breakthrough times at low perme-abilities.In this work,we use the dispersion relations given by Eqs.13and14to complete the averaged model.(Remark:As aﬁrst approximation,we have assumed that the mass transfer coefﬁcient to be the same in the axial and transverse directions. However,as in the case of the dispersion coefﬁcient,the convective contribution to mass transfer could be different in theﬂow and transverse directions.This can be accounted for by replacing the scalar k c by the transfer matrix(a tensor).We do not pursue this here.)Dimensionless Model Equations and Limiting CasesThe model equations forﬁrst-order irreversible kinetics are made dimensionless for the case of constant injection rate at the inlet boundary by deﬁning the following dimensionless vari-ablesxϭxЈL,yϭyЈL,zϭzЈL,uϭUu o,tϭtЈ͑L/u o͒,rϭr po,A␯ϭa␯o,␬ϭKo,c fϭC fC o,c sϭC sC o,pϭPϪP e␮u o LK o␾2ϭ2k s r oD m,Daϭk s a o Lu o,N acϭ␣C o␳s,Pe Lϭu o LD m,␩ϭ2r oL,␣oϭHLwhere L is the characteristic length scale in the(ﬂow)xЈdirection,H is the height of the domain,u o is the inlet velocity, C o is the inlet concentration of the acid and P e is the pressure at the exit boundary of the domain.The initial values of permeability,interfacial area and average pore radius are rep-Table1.Pore Level Peclet Numbers at DifferentInjection RatesRegime Injection Velocity(cm/s)Pep Face 1.4ϫ10Ϫ47ϫ10Ϫ4Wormhole 1.4ϫ10Ϫ37ϫ10Ϫ3Uniform0.140.73236AIChE JournalDecember2005Vol.51,No.12。

DISTA, Universita del Piemonte Orientale,

The Scale Factor:A New Degree of Freedom in Phase Type ApproximationAndrea BobbioDISTA,Universit`a del Piemonte Orientale, Alessandria,Italy,bobbio@unipmn.itAndr´a s Horv´a th,Mikl´o s TelekDept.of Telecommunications, Budapest University of Technology and Economics, Hungary,horvath,telek@webspn.hit.bme.huAbstractThis paper introduces a uniﬁed approach to phase-type approximation in which the discrete and the continuous phase-type models form a common model set.The models of this common set are assigned with a non-negative real parameter,the scale factor.The case when the scale factor is strictly positive results in Discrete phase-type distribu-tions and the scale factor represents the time elapsed in one step.If the scale factor is0,the resulting class is the class of Continuous phase-type distributions.Applying the above view,it is shown that there is no qualitative difference be-tween the discrete and the continuous phase-type models.Based on this uniﬁed view of phase-type models one can choose the best phase-type approximation of a stochastic model by optimizing the scale factor.Keywords:Discrete and Continuous Phase type distri-butions,Phase type expansion,approximate analysis.1IntroductionThis paper presents new comparative results on the use of Discrete Phase Type(DPH)distributions[11]and of Continuous Phase Type(CPH)distributions[12]in applied stochastic modeling.DPH distributions of order are deﬁned as the time to absorption in a Discrete-State Discrete-Time Markov Chain (DTMC)with transient states and one absorbing state. CPH distributions of order are deﬁned,similarly,as the distribution of the time to absorption in a Discrete-State Continuous-Time Markov Chain(CTMC)with transient states and one absorbing state.The above deﬁnition im-plies that the properties of a DPH distribution are computed over the set of the natural numbers while the properties of a CPH distribution are deﬁned as a function of a continuous time variable.When DPH distributions are used to model timed activities,the set of the natural numbers must be re-lated to a time measure.Hence,a new parameter need to be introduced that represents the time span associated to each step.This new parameter is the scale factor of the DPH dis-tribution,and can be viewed as a new degree of freedom, since its choice largely impacts the shape and properties of a DPH distribution over the continuous time axes.When DPH distributions are used to approximate a given continu-ous distribution,the scale factor affects the goodness of the ﬁt.The paper starts discussing to what extent DPH or CPH distributions can be utilized toﬁt a given continuous distri-bution.It is shown that a DPH distribution of any order con-verges to a CPH distribution of the same order as the scale factor goes to zero.Even so,the DPH class contains dis-tributions whose behavior differs substantially from the one of the corresponding distributions in the CPH class.Two main peculiar points differentiate the DPH class from the CPH class.Theﬁrst point concerns the coefﬁcient of varia-tion:indeed,while in the continuous case the minimum co-efﬁcient of variation is a function of the order only and its lower bound is given by the well known theorem of Aldous and Shepp[1],in the discrete case the minimum coefﬁcient of variation is proved to depend both on the order and on the mean(and hence on the scale factor)[13].Furthermore, it is easy to see that for any order,there exist members of the DPH class that represent a deterministic value with a coefﬁcient of variation equal to zero.Hence,for any order (greater than1),the coefﬁcient of variation of the DPH class spans from zero to inﬁnity.The second peculiar point that differentiate the DPH class is the support of the distributions.While a CPH dis-tribution(of any order)has always an inﬁnite support,there exist members of the DPH class of any order that have a ﬁnite support(between a minimum non-negative value and a maximum)or have a mass equal to one concentrated in a single value(deterministic distribution).It turns out that the possibility oftuning the scale factor to optimize the goodness of the ﬁt,having distributions with coefﬁcient of variation span-ning from0to inﬁnity,representing deterministic values exactly,coping withﬁnite support distributions,makes the DPH class a very interesting and challenging class of distributions to be explored in applied stochastic models.The purpose of this paper is to show how these fa-vorable properties can be exploited in practice,and to pro-vide guidelines to the modeler to a reasonably good choice of the distributions to be used.Indeed,since a DPH dis-tribution tends to a CPH distribution as the scale factor ap-proaches zero,considering the scale factor as a new decision variable in aﬁtting experiment,andﬁnding the value of the optimal scale factor(with respect to some error measure) provides a valuable tool to decide whether to use a discrete or a continuous approximation to the given problem.Theﬁtting problem for the CPH class has been exten-sively studied and reported in the literature by resorting to a variety of structures and numerical techniques(see[10]for a survey).Conversely,theﬁtting problem for the DPH class has received very little attention[4].In recent years,a considerable effort has been devoted to deﬁne models with generally distributed timings and to merge in the same model random variables and determin-istic duration.Analytical solutions are possible in special cases,and the approximation of the original problems by means of CPH distributions is a rather well known tech-nique[7].This paper is aimed at emphasizing that DPH approximation may provide a more convenient alternative with respect to CPH approximation,and also to provide a way to quantitatively support this choice.Furthermore,the use of DPH approximation can be extended from stochas-tic models to functional analysis where time intervals with nondeterministic choice are considered[3].Finally,dis-cretization techniques for continuous problems[8]can be restated in terms of DPH approximations.The rest of the paper is organized as follows.After deﬁn-ing the notation to be used in the paper in Section2,Section 3discusses the peculiar properties of the DPH class with re-spect to the CPH class.Some guidelines for bounding the parameters of interest and extensive numerical experiments to show how the goodness of theﬁt is inﬂuenced by the op-timal choice of the scale factor are reported in Section4. Section5discusses the quality of the approximation when passing from the analysis of a single distribution to the anal-ysis of performance measures in complete non-Markovian stochastic models.The paper is concluded in Section6.2Deﬁnition and NotationA DPH distribution[11,12]is the distribution of the time to absorption in a DTMC with transient states,and one absorbing state numbered.The one-step transition probability matrix of the corresponding DTMC can be par-titioned as:(1)where is the matrix collecting the transi-tion probabilities among the transient states,is the column vector of length grouping the probabilities from any state to the absorbing one,and is the zero vector.The initial probability vectoris of length,with.In the present paper,we consider only the class of DPH distribu-tions for which,but the extension to the case when is straightforward.The tuple is called the representation of the DPH distribution,and the order.Similarly,a CPH distribution[12]is the distribution of the time to absorption in a CTMC with transient states, and one absorbing state numbered.The inﬁnites-imal generator of the CTMC can be partitioned in the following way:(2) where,is a matrix that describes the tran-sient behavior of the CTMC and is the column vector grouping the transition rates to the absorbing state.Letbe the initial probability(row) vector with.The tuple is called the representation of the CPH distribution,and the order.It has been shown in[4]for the discrete case and in[6] for the continuous case that the representations in(1)and (2),because of their too many free parameters,do not pro-vide a convenient form for running aﬁtting algorithm.In-stead,resorting to acyclic phase-type distributions,the num-ber of free parameters is reduced signiﬁcantly since both in the discrete and the continuous case a canonical form can be used.The canonical form and its constraints for the discrete case[4]is depicted in Figure1.Figure2gives the canonical form and associated constraints for the continuous case.In both cases the canonical form corresponds to a mixture of Hypo-exponential distributions.Aﬁtting algorithm that provides acyclic CPH,acyclic DPH distributions has been provided in[2]and[4],respec-tively.Experiments suggests(an exhaustive comparison of ﬁtting algorithms can be found in[10])that,from the point of view of applications,the Acyclic phase-type class is as ﬂexible as the whole phase-type class.3Comparing properties of CPH and DPH distributionsCTMC are deﬁned as a function of a continuous time variable,while DTMC are deﬁned over the set of the nat-ural numbers.In order to relate the number of jumps in a DTMC with a time measure,a time span must be assigned to each step.Let be(in some arbitrary units)the scaleFigure2.Canonical representation of acyclicCPH distributions and its constraintsfactor,i.e.the time span assigned to each step.The valueof establishes an equivalence between the sentence”prob-ability at the-th step”and”probability at time”,and hence,deﬁnes the time scale on which the properties of theDTMC are measured.The consideration of the scale factor introduces a new parameter,and consequently a new de-gree of freedom,in the DPH class with respect to the CPHclass.In the following,we discuss how this new degree of freedom impacts the properties of the DPH class and how it can be exploited in practice.Let be an”unscaled”DPH distributed random variable (r.v.)of order with representation,deﬁned over the set of the non-negative natural numbers.Let us consider a scale factor;the scaled r.v.is deﬁned over the dis-crete set of time points,being a non-negative natural number.For the unscaled and the scaled DPH r.v.the following equations hold.(3)where is the column vector of ones,and is the -th moment calculated from the factorial moments of:.It is evident from(3)that the mean of the scaled r.v.is times the mean of the unscaled r.v..While is an invariant of the representation,is a free parame-ter;adjusting,the scaled r.v.can assume any mean value .On the other hand,one can easily infer from(3) that the coefﬁcients of variation of and are equal.A consequence of the above properties is that one can easily provide a scaled DPH of order with arbitrary mean and arbitrary coefﬁcient of variation with an appropriate scale factor.Or more formally:the unscaled DPH r.v.of any order can exhibit a coefﬁcient of variation between .For the coefﬁcient of variation ranges between.As mentioned earlier,an important property of the DPH class with respect to the CPH class is the possibility of exactly representing a deterministic delay.A determinis-tic distribution with value can be realized by means of a scaled DPH distribution with phases with scale factor if is integer.In this case,the structure of the DPH distribution is such that phase is connected with probabil-ity1only to phase(),and with an initial probability concentrated in state1.If is not inte-ger for the given,the deterministic behavior can only be approximated.3.1First order discrete approximation of CTMCsGiven a CTMC with inﬁnitesimal generator,the tran-sition probability matrix over an interval of length can be written as:hence theﬁrst order approximation of is matrix.is a proper stochastic matrix if,where.is the exact transition probability matrix of the CTMC assumed that at most one transition occurs in the interval of length.We can approximate the behavior of the CTMC at timeusing the DTMC with transition probability matrix.The approximate transition prob-ability matrix at time is:.Since matrices and commute we can obtain the matrix version of the same expression as followsAn obvious consequence of Theorem1for PH distribu-tions is given in the following corollary.Corollary1Given a scaled DPH distribution of order, representation and scale factor,the limiting behavior as is the CPH distribution of order with representation.3.2The minimum coefﬁcient of variationIt is known that one of the main limitation in approx-imating a given distribution by a PH one is the attainable minimal coefﬁcient of variation,.In order to discuss this point,we recall two theorems that state the for the class of CPH and DPH distributions.Theorem2(Aldous and Shepp[1])The of a CPH distributed r.v.of order is and is attained by the Erlang()distribution independent of its mean or of its parameter.The corresponding theorem for the unscaled DPH class has been proved in[13].In the following,denotes the integer part and denotes the fractional part of. Theorem3The of an unscaled DPH r.v.of order and mean is:In this particular case,when the structure of the bestﬁt-ting scaled DPH and CPH distributions are known,we can show that the distribution of the bestﬁtting scaled DPH dis-tribution converges to the distribution of the bestﬁtting CPH distribution when.Unfortunately,the same conver-gence property cannot be proved in general,since the struc-tural properties of the bestﬁtting PH distributions are not known and they depend on the chosen(arbitrary)optimiza-tion criterion.Instead,in Section4we provide an extensive experimental study on the behavior of the bestﬁtting scaled DPH and CPH distributions as a function of the scale factor .3.4DPH distributions withﬁnite supportAnother peculiar characteristic of the DPH class is to contain distributions withﬁnite support.A DPH distribu-tion hasﬁnite support if its structure does not contain cycles and self-loops(any cycle or self loop implies an inﬁnite sup-port).Let be theﬁnite support of a given distribution,with and(when theﬁnite support distri-bution reduces to a deterministic distribution with mass1at ).If and are both integers,it is possible to construct a scaled DPH of order for which the probabil-ity mass function has non-zero elements only for the values.As an example,the discrete uniform distribution between and is reported in Figure 5,for scale factor.0.10.20.30.40.50.60.70.80.910.40.60.811.2 1.41.61.82c d fxOriginalScale factor: 0.01Scale factor: 0.06Scale factor: 0.1CPH 00.511.522.50.40.60.81 1.2 1.4 1.6 1.82p d fxOriginalScale factor: 0.01Scale factor: 0.06Scale factor: 0.1CPH Figure 6.Approximating the L3distribution with -phase PH approximationsWhen is less than its lower bound the required can-not be attained;when becomes too large the wide separa-tion of the discrete steps increases the approximation error;when is in the proper range (e.g.)a reasonably good ﬁt is achieved.This example also suggests that an optimal value of exists that minimizes the chosen distance measure in (6).In order to display the goodness of ﬁt for the L3distribu-tion,Figure 7shows the distance measure as a function of for various values of the order .A minimum value of is attained in the range where the parameters ﬁt the bounds of Table 1.Notice also that,as increases,the advantage of having more phases disappears,according to Theorem 3.The circles in the left part of this ﬁgure (as well as in all the successive ﬁgures)indicate the corresponding distance measure obtained from CPH ﬁtting.The ﬁgure (and the subsequent ones as well)suggests that the distance measure obtained from DPH ﬁtting converges to the distance mea-sure obtained by the CPH approximation as tends to .upper bound of equation (7)0.20920.07920.04250.021700.010.020.030.040.050.060.070.080.090.100.050.10.150.20.250.3d i s t a n ce m e a s u r escale factor2 phases 4 phases 6 phases 8 phases 10 phases 12 phasesFigure 9.Distance measure as the function of the scale factor for Uniform(1,2)(U2)00.0020.0040.0060.0080.010.0120.0140.01600.050.10.150.20.250.3d i s t a n ce m e a s u r escale factor2 phases 4 phases 6 phases 8 phases 10 phases 12 phasesFigure 10.Distance measure as the function of the scale factor for Uniform(0,1)(U1)be stressed that the chosen distance measure in (6)can be considered as not completely appropriate in the case of ﬁnite support,since it does not force the approximating PH to have its mass conﬁned in the ﬁnite support and 0outside.Let be a uniform r.v.over the interval ,withand (this is the distributionU2taken from the benchmark in [5,4]).Figure 9shows the distance measure as a function of for various orders .It is evident that,for each ,a minimal value of is obtained,that provides the best approximation according to the chosen distance measure.As a second example,let be a uniform r.v.over theinterval,with and (this is the distribution U1taken from the benchmark in [5,4]).Figure 10shows the distance measure as a function of forvarious orders .Since,in this example,an order is large enough for a CPH to attain the coefﬁcient of variation of the distribution.Nevertheless,the optimal in Figure (10),which minimizes the distance mea-sure for high order PH (),ranges between and ,thus leading to the conclusion that a DPH provides a better ﬁt.This example evidences that the coef-ﬁcient of variation is not the only factor which inﬂuences the optimal value.The shape of the distribution plays an essential role as well.Our experiments show that a discon-00.20.40.60.810.20.40.60.811.21.4c d fxOriginalScale factor: 0.03Scale factor: 0.1CPH00.20.40.60.811.21.41.600.20.40.60.811.21.4p d fxOriginalScale factor: 0.03Scale factor: 0.1CPHFigure 11.Approximating the Uniform ()distribution (U1)tinuity in the pdf (or in the cdf)is hard to approximate with CPH,hence in the majority of these cases DPH provides a better approximation.Figure 11shows the cdf and the pdf of the U1distribu-tion,compared with the best ﬁt PH approximations of order,and various scale factors .In the case of DPH ap-proximation,the values are calculated as in (9).With respect to the chosen distance measure,the best approxi-mation is obtained for,which corresponds to a DPH distribution with inﬁnite support .Whenthe approximate distribution has a ﬁnite support.Hence,the value (for )provides a DPH able to rep-resent the logical property that the random variable is less than .Another ﬁtting criterion may,of course,stress this property.5Approximating non-Markovian modelsSection 4has explored the problem of how to ﬁnd the best ﬁt among either a DPH or a CPH distribution by tuning the scale factor .When dealing with a stochastic model of a system that incorporates non exponential distributions,a well know solution technique consists in a markovianiza-tion of the underlying non-Markovian process by substi-tuting the non exponential distribution with a best ﬁt PH distribution,and then expanding the state space.A natural question arises also in this case,on how to decide among a discrete (using DPH)or a continuous (using CPH)approx-imation,in order to minimize the error in the performances2s4s3Figure12.The state space of the consideredM/G/1/2/2queuemeasures we are interested in for the overall model.One possible way to handle this problem could consist inﬁnding the best PHﬁts for any single distribution and to plug them in the model.In the present paper,we only consider the case where the PH distributions are either all discrete(and with the same scale factor)or they are all continuous.Various embedding techniques have been ex-plored in the literature for mixing DPH(with different scale factors)and CPH([8,9]),but these techniques are out of the scope of the paper.In order to quantitatively evaluate the inﬂuence of the scale factor on some performance measures deﬁned at the system level,we have considered a preemptive M/G/1/2/2 queue with two classes of customers.We have chosen this example because accurate analytical solutions are available both in transient condition and in steady-state using the methods presented in e.g.[8].The general distribution is taken from the set of distributions(L1,L3,U1,U2)already considered in the previous section.Customers arrive at the queue with rate in both classes.The service time of a higher priority job is exponen-tially distributed with parameter.The service time distribution of the lower priority job is either L1,L3,U1 or U2.Arrival of a higher priority job preempts the lower priority one.The policy associated to the preemption of the lower priority job is preemptive repeat different(prd),i.e. after the departure of the higher priority customer the ser-vice of the low priority customer starts from the beginning with a new service time sample.The system has4states(Figure12):in state s1the server is empty,in state s2a higher priority customer is under ser-vice with no lower priority customer in the system,in state s3a higher priority customer is under service with a lower priority customer waiting,in state s4a lower priority job is under service(in this case there cannot be a higher priority job).Let denote the steady state probability of the M/G/1/2/2queue obtained from an exact analytical solution.In order to evaluate the correctness of the PH approxima-0.020.040.060.080.10.120.140.1600.020.040.060.080.10.120.140.160.180.2sumoferrorsscale factor2 phases4 phases6 phases8 phases10 phases12 phasesFigure13.with scale factor and distri-bution L30.010.020.030.040.050.0600.020.040.060.080.10.120.140.160.180.2sumoferrorsscale factor2 phases4 phases6 phases8 phases10 phases12 phasesFigure14.with scale factor and distri-bution L3tion we have solved the model by substituting the original general distribution(either L1,L3,U1or U2)with approx-imating DPH or CPH distributions.Letdenote the steady state probability of the M/PH/1/2/2queue with the PH approximation.The overall approximation error is measured in terms of the difference between the exact steady state probabilities and the approximate steady state probabilities.Two error measures are deﬁned:andThe evaluated numerical values for and are reported in Figures13and14for the distribution L3.Since the behavior of is very similar to the behavior of in all the cases,for the other distributions we reportonly(Figures15,16,17).Theﬁgures,which re-fer to the error measure in a performance index of a global stochastic model,show a behavior similar to the one ob-tained for a single distributionﬁtting.Depending on the coefﬁcient of variation and on the shape of the considered non-exponential distributions an optimal value of is found which minimizes the approximation error.In this example, the optimal value of is close to the one obtained for the single distributionﬁtting.Based on our experiments,we guess that the observed0.050.10.150.20.2500.020.040.060.080.10.120.140.160.180.2s u m o f e r r o r sscale factor2 phases 4 phases 6 phases 8 phases 10 phases 12 phasesFigure 15.with scale factorand distri-bution L100.020.040.060.080.10.120.140.1600.050.10.150.20.250.3s u m o f e r r o r sscale factor2 phases 4 phases 6 phases 8 phases 10 phases 12 phasesFigure 16.with scale factorand distri-bution U1property is rather general.If the stochastic model under study contains a single non-exponential distribution,then the approximation error in the evaluation of the perfor-mance indices of the global model can be minimized by re-sorting to a PH type approximation (and subsequent DTMC or CTMC expansion)with the optimal of the single distri-bution.The same should be true if the stochastic model under study contains more than one general distribution,whose best PH ﬁt provides the same optimal .In order to investigate the approximation error in the transient behavior,we have considered distribution U2for the service time and we have computed the transient proba-bility of state with two different initial conditions.Figure 18depicts the transient probability of state with initial state .Figure 19depicts the transient probability of the same state,,when the service of a lower priority job starts at time 0(the initial state is ).All approximations are with DPH distributions of order .Only the DPH ap-proximations are depicted because the CPH approximationis very similar to the DPH one with scale factor.In the ﬁrst case,(Figure 18),the scale factor,which was the optimal one from the point of view of ﬁtting the single distribution in isolation,provides the most accu-rate results for the transient analysis as well.Instead,in the second case,the approximation with a scale factor0.020.040.060.080.10.120.140.160.180.200.050.10.150.20.250.3s u m o f e r r o r sscale factor2 phases 4 phases 6 phases 8 phases 10 phases 12 phasesFigure 17.with scale factorand distri-bution U20.10.20.30.40.50.60.70.80.91012345t r a n s i e n t p r o b a b i l i t ytimeTransient behaviour Scale factor: 0.03Scale factor: 0.1Scale factor: 0.2Figure 18.Approximating transient probabili-tiescaptures better the sharp change in the transient probability.Moreover,this value of is the only one among the values reported in the ﬁgure that results in 0probability for time points smaller than 1.In other words,the second example depicts the advantage given by DPH distributions to model durations with ﬁnite support.This example suggests also that DPH approximation can be of importance when pre-serving reachability properties is crucial (like in modeling time-critical systems)and,hence,DPH approximation can be seen as a bridge between the world of stochastic model-ing and the world of functional analysis and model checking [3].6Concluding remarksThe main result of this paper has been to show that the DPH and CPH classes of distributions of the same order can be considered a single model set as a function of a scalefactor .The optimal value of ,,determines the best distribution in a ﬁtting experiment.When the bestchoice is a CPH distribution,while whenthe best choice is a DPH distribution.This paper has also shown that the transition from DPH class to CPH class is continu-ous with respect to several properties,like the distance (de-noted by in 6)between the original and the approximate0.050.10.150.20.250.3012345t r a n s i e n t p r o b a b i l i t ytimeTransient behaviour Scale factor: 0.03Scale factor: 0.1Scale factor: 0.2Figure 19.Approximating transient probabili-tiesdistributions.The paper presents limit theorems for special cases;however,extensive numerical experiments show that the limiting behavior is far more general than the special cases considered in the theorems.The numerical examples have also evidenced that for very small values of ,the diagonal elements of the tran-sition probability matrix become very close to ,rendering numerically unstable the DPH ﬁtting procedure.A deep analytical and numerical sensitivity analysis is required to draw more general conclusions for the model level “optimal value”and its dependence on the consid-ered performance measure than the ones presented in this work.It is deﬁnitely a ﬁeld of further research.Finally,we summarize the advantages and the disadvan-tages of applying approximate DPH models (even with op-timal value)with respect to using CPH approximations.Advantages of using DPH:An obvious advantage of the ap-plication of DPH distributions is that one can have a closer approximate of distributions with low coefﬁcient of varia-tion.An other important quantitative property of the DPH class is that it can capture distributions with ﬁnite support and deterministic values.This property allows to capture the periodic behavior of a complex stochastic model,while any CPH based approximation of the same model tends to a steady state.Numerical experiments have also shown that DPH can better approximate distributions with some abrupt or sharp changes in the CDF or in the PDF.Disadvantages of using DPH:There is a deﬁnite disad-vantage of discrete time approximation of continuous time models.In the case of CPH approximation,coincident events do not have to be considered (they have zero proba-bility of occurrence).Instead,when applying DPH approxi-mation coincident events have to be handled,and their con-sideration may burden signiﬁcantly the complexity of the analysis.AcknowledgmentsThis work has been performed under the Italian-Hungarian R&D program supported by the Italian Ministry of Foreign Affairs and the Hungarian Ministry of Education.A.Bob-bio was partially supported by the MURST Under Grant ISIDE;M.Telek was partially supported by Hungarian Sci-entiﬁc Research Fund (OTKA)under Grant No.T-34972.References[1] D.Aldous and L.Shepp.The least variable phase type dis-tribution is Erlang.Stochastic Models ,3:467–473,1987.[2] A.Bobbio and A.Cumani.ML estimation of the param-eters of a PH distribution in triangular canonical form.In G.Balbo and G.Serazzi,editors,Computer Performance Evaluation ,pages 33–46.Elsevier Science Publishers,1992.[3] A.Bobbio and A.Horv´a th.Petri nets with discrete phasetype timing:A bridge between stochastic and functional analysis.Electronic Notes in Theoretical Computer Science ,52(3),2001.[4] A.Bobbio,A.Horv´a th,M.Scarpa,and M.Telek.Acyclicdiscrete phase type distributions:Properties and a parameter estimation algorithm.Technical Report of Budapest Uni-versity of Technology and Economics,Submitted to Perfor-mance Evaluation,2000.[5] A.Bobbio and M.Telek.A benchmark for PH estima-tion algorithms:results for Acyclic-PH.Stochastic Models ,10:661–677,1994.[6] A.Cumani.On the canonical representation of homoge-neous Markov processes modelling failure-time distribu-tions.Microelectronics and Reliability ,22:583–602,1982.[7] A.Cumani.Esp -A package for the evaluation of stochas-tic Petri nets with phase-type distributed transition times.In Proceedings International Workshop Timed Petri Nets ,pages 144–151,Torino (Italy),1985.[8]R.German.Performance Analysis of Communication Sys-tems:Modeling with Non-Markovian Stochastic Petri Nets .John Wiley and Sons,2000.[9]R.Jones and G.Ciardo.On phased delay stochastic Petrinets:Deﬁnition and application.In Proceedings 9th Inter-national Workshop on Petri Nets and Performance Models -PNPM01.IEEE Computer Society,2001.[10] ng and J.L.Arthur.Parameter approximation forphase-type distributions.In Matrix-analytic methods in stochastic models ,Lecture notes in pure and applied mathe-matics,pages 151–206.Marcel Dekker,Inc.,1996.[11]M.Neuts.Probability distributions of phase type.In LiberAmicorum Prof.Emeritus H.Florin ,pages 173–206.Uni-versity of Louvain,1975.[12]M.Neuts.Matrix Geometric Solutions in Stochastic Models .Johns Hopkins University Press,Baltimore,1981.[13]M.Telek.Minimal coefﬁcient of variation of discretephase type distributions.In 3rd International Conference on Matrix-Analitic Methods in Stochastic models,MAM3,pages 391–400,Leuven,Belgium,2000.Notable Publica-tions Inc.。

3给水排水外文翻译外文文献英文文献

3给水排水外文翻译外文文献英文文献Relations between triazine flux, catchment topography and distancebetween maize fields and the drainage network F. Colina,*, C. Puecha, G. de Marsilyb,1aUMR “Syste`mes et Structures Spattiaux”, Cemagref-ENGREF 500, rue J.F. Breton 34093, Montpellier Cedex 05, FrancebUMR “Structure et Fonctionement des Syste`mes Hydriques Continentaux”, Universite´ P. et M. Curie 4, Pl. Jussieu 75252, Paris Cedex 05, FranceReceived 5 October 1999; revised 27 April 2000; accepted 19 June 2000AbstractThis paper puts forward a methodology permitting the identification of farming plots contributing to the pollution of surface water in order to define the zones most at risk from pesticide pollution. We worked at the scale of the small agricultural catchment (0.2–7.5 km2) as it represents the appropriate level oforganisation for agricultural land. The hypothesis tested was: the farther a field undergoing a pesticide treatment is from a channel network, the lower its impact on pollution at the catchment outlet.The study area, the Sousson catchment (120 km2, Gers, France), has a “herring bone” structure: 50 independent tributaries supply the main drain. Pesticide sales show that atrazine is the most frequently used compound although it is only used for treating maize plots and that its application rate is constant. In two winter inter-storm measurement exercises, triazine flux values were collected at about 30 independent sub-basin outlets.The contributory areas are defined, with the aid of a GIS, as different strips around the channel network. The correlation between plots under maize in contributory zones and triazine flux at related sub-basin outlets is studied by using non-parametric and linear correlation coefficients. Finally, the most pertinentcontributory zone is associated with the best correlation level.A catchment typology, based on a slope criterion, allows us to conclude that in steep slope catchments, the contributory area is best defined as a 50 m wide strip around the channel network. In flat zones, the agricultural drainage network is particularly well developed: artificial drains extend the channel network extracted from the 1/25.000 scale topographic map, and the total surface area of the catchment must be taken to account. q 2000 Elsevier Science B.V. All rights reserved.Keywords: Pesticide catchment; GIS artificial network1. IntroductionThe use of pesticides in western agriculture dates back to the middle of the 19th century (Fournier,1988). Since then, because of their intensive use,yields have increased and the demand for agricultural products has been satisfied. However, the pollution created by theiruse threatens both drinking water resources and the integrity of ecosystems. Therefore, there is a great demand for the reduction of pollution.The remedies lie in changes in the way that agricultural land is managed. The problem of agricultural Journal non-point source pollution by pesticides must be taken from the field, the level of action, to the catchment,the level of control of the water resource.Between these two spatial scales, different levels of organisation can be found. Fields, groups of fields,basins and main catchment, can be viewed together as nested systems (Burel et al., 1992). For each scale level, the main processes governing water movement and soluble pollutant transport are different, as are the variables characterising the system (Lebel, 1990):flow in macropores at local scale, preferential flowpaths at the hillslope scale, flows in connection withthe repartition of different soils at the catchment scale,geology influence at the regional scale(Blo¨sch and Sivapalan, 1995).At the field level, an experimental approach can be used and the relative weight of each variable can be experimentally tested (Scheunert, 1996; Bengtson et al., 1990). The major factors that concern agricultural practices have been identified and many agricultural management indicators have been developed (Bockstaller et al., 1997). Nevertheless, this approach cannot be applied at the catchment scale for several reasons: the need to measure the pollution and the environmental factors simultaneously, multiple measurement difficulties, the complexity of analysis. The variability of observations has temporal and spatial components. Rain induces pesticide leaching and therefore causes temporary high pesticide concentrations in the water; the closer the pesticide spreading date in thefield is to the measurement, the greater the concentration levels (Seux et al., 1984; Reme,1992; Laroche and Gallichand, 1995). The extensive use of Geographical Information System (GIS) has made it possible to analyse the impact on the pollution of the spatial characteristics of agricultural zones (Battaglin and Goolsby, 1996). But so far, the results of these experimentshave only led to an approximate estimate of the risks (Tim and Jolly, 1994).In order to progress in the search for ways to reduce pesticide pollution, it would be worthwhile to improve our assessment of how spatial structure and organisation affects the levels of pollutants measured.This paper presents the results of a study that concerns a particular aspect of the influence of spatial organisation on pesticide transfer: the effects of the distance between the cropland and the channel network. The longer the distance between a cultivated field and a river, the greater the retention and degradation processes (Leonard, 1990; Belamie et al.,1997). One mighttherefore imagine that the greater the distance, the lower the pollution level. However,few studies have given a numerical value to the critical distance at which a field does not influence river pollution significantly. Usually, when dealing with risk zone definition, experts establish an arbitrary distance (Bouchardy, 1992). Our main goal is to determine through spatial analysis the critical distance from a hydrographic network. The zones most at risk from pesticides, including the plots, which contribute most of the pollution, can then be determined.The study area, the Sousson catchment (Gers,France) has certain physical characteristics, which allows sampling of most of the independent subbasins, defined here as agricultural production zones. Its particular morphology made the comparative study of the production zones possible. The method involves a statistical comparison between pollution measurements and spatial characteristics of thecatchments. In order to establish the boundaries ofthe contributing areas, the pollution flux measured at the production zone outlet is compared to the landcover, estimated within strips of variable width around the channel network. Results are shown and discussed from a mainly practical viewpoint.2. The study area and collected data2.1. Study area descriptionThe study area is the Sousson catchment, in southwestern France (Gers). The Sousson River is a tributary of the river Gers. The catchment area is 120 km2. The 32 km long hydrographic network has a ‘herringbone’pattern: 53 sub-basins with fairly homogeneous surfaces areas ranging from 0.2 to 7.5 km2 serve the central drain (Fig. 1).The wide, gently sloping and heavily cultivated left bank, differs from the right bank, which is narrow, steep and mainly made up offorest and pastureland.The Sousson catchment area is exclusively agricultural.There is no industry or settlement of more than 200 inhabitants. The two main crops cultivated aremaize and winter wheat (17 and 15% of the catchment surface area, respectively). The maize fields are usually situated, on the left bank, in the upstream middle of the catchment area, and along the main river.There are two types of soil: a calcareous soil, which is quite permeable, and a non-calcareous soil called locally ‘boulbenes’ with an top limoneous layer and a lower silty layer. In order to avoid the stagnation of water in the upper layer caused by the silty impermeable layer, the fields on boulbene soil are artificially drained. Maize is cultivated for preference on thistype of soil.No significant aquifer has been found in the catchment, as the substratum is rather impervious (clays).2.2. Collected data2.2.1. Spatial dataA GIS was developed for the area, which contains the following information layers:²the hydrographic network and the catchment boundaries digitized from 1/25.000 scale topographic map;² a gridded Digital Elevation Model (DEM) of the zone providing landsurface slopes generated from DEM with a resolution of 75 m;²the boundaries of cultivated fields digitized from aerial photos at scale of 1/15.000;²landcover for both 1995 and 1996 was defined in detail in the study area. For 1997, landcover was identified by remote sensing. Knowledge of agricultural antecedents enhanced the classification of a SPOT (Satellite Pour l0Observation de la Terre) image. As a result, the maize areas for the entire Sousson catchment were determined for 1995, 1996 and 1997 (Fig. 2).GIS functions are capable of determining the landcover of each catchment by intersecting the two information layers “landcover” and “catchment boundaries”, or defining a zone of constant width around the hydrographic network, which is called the buffer zone.In order to evaluate the pesticide application rate, figures for local pesticide sales were collected. Atrazine, alachlor and glyphosate are the most commonly used compounds, atrazine far outstrips the others triazines as the most frequently used product (ten times less simazine is sold). In this region, atrazine is only used in maize cultivation. The application rate (mass of atrazine sold/maize surface area) does not vary from one municipality to another.To simplify the investigations, we chose to study the atrazine spread on maize plots in May. We assume that all the maize plots are treated with atrazine and that the application rate is uniform.2.2.2. Water pollution dataTwo series of measurements were made during the winter period: 23 sub-basins were sampled on December 3rd and 4th 1997, and 26 sub-basins were sampled March 17th to 19th 1998. Hence, the atrazine treatments were carried out 7 or 10 months before and the maize harvest was 1 and 4 months before the measurements were taken.To obtain stable hydrological conditions, the chosen measurement dates coincided with decreasing flow as shown in Fig. 3. The same operator collected the quality samples and gauged the river flow in order to limit measurement errors.The triazine concentration was measured with an ELISA water test (Transia Plate PE 0737). This measurement technique is less accurate than the classical chromatography technique, but it permits a faster analysis of a large number of samples (Rauzy and Danjou, 1992; Lentza-Rios, 1996). As atrazine is the mostwidely commercialised triazine product in this region, we will consider that observed triazine concentrations are representative of atrazine concentrations.December 1997 values, and March 1998 values were grouped together in order to assemble a large enough sample for statistical analysis (Fig. 4). The instantaneous triazine flux was obtained by multiplying the triazine concentration with the dischargevalue. As shown in Table 1, water flow in December 1997 was double that in March 1998, but the corresponding triazine flux are comparable.2.2.3. Quality assuranceTo control the quality of ELISA water-test measurements, each concentration was analysed 142 F. Colin et al. / Journal of Hydrology 236 (2000) 139–152 Fig. 2. Hydrographic network (topographic 1/25.000 map) and subcatchments, parcel limits and land-cover (example of maizeplots). twice. A maximum difference of 20% is tolerated between two duplicate samples, the median error is 10%, and mean values are used. It is possible that ELISA measurement induces a consistent error by comparing with gas chromatography measurements (Tasli et al., 1996), but this bias is compensated by comparative reasoning on all the samples.A few points were measured two or three times during the exercise in order to evaluate the daily variations during the sampling period. Table 2 shows that the flux variation between different days of a sampling period ranges from 2 to 49%. It is therefore possible to compare the different samples from the period in question. All the measurements from each period are then grouped together.The uncertainty on the triazine flux is the sum of the uncertainty of discharge and concentration measurements. The uncertainty on the discharge measurements ranges from 15 to 20%. Therefore, the triazine flux value isgiven with a maximum uncertainty of 40%.3. MethodTo define the zones most at risk we tested how the distance to the river of the areas where pesticides are applied influence pollution levels. Thus, we have to determine the relative position of the hydrographic network and the contaminating plots.In our case, the data on pollution is provided by triazine flux measurements taken at basin outlets and the potentially contaminating fields are maize plots.3.1. Efficiency curve and spatial partitionThe basic hypothesis is that the impact of the field as a contributor to pollution decreases the further it is from the channel network. Thus, there is a critical distance at which the field makes little contribution to outlet pollution. In other words, we assume that plot contribution to pollution level can be modelled through adecreasing efficiency curve. This hypothesis will be tested with a very simple curve: a step function. This curve is defined using only one parameter, the threshold limit distance, d, beyond, which a plot stops contributing to river pollution.In practice, this hypothesis implies a three-step approach:²determination of the location of the maize fields;²definition of a buffer of width d, equal to the threshold distance and, which surrounds the channel network;²determination of the contaminating fields inside these limits.The fields define the contributing maize areas depending on the buffer width (Fig. 5). At this stage, GIS functionality is required, particularly for the buffer function.3.2. Correlation between contributing area and pollution at the catchment outletWe studied the correlation level between triazine flux measured at the catchment outlet and the different contamination contributing areas defined by strips of variable width. Three parameters are used to determine the correlation level (further information is provided on this point in Appendix A):²The Kendall rank correlation coefficient (Siegel, 1956) t gives a measure of the degree of association or correlation between two sets of ranks. It expresses the difference between the probability that the two data sets are ranked according to the same order and the probability that they are ranked according to a different order. If t . 1.21.; a positive (negative) relation exists between the two data series, if t . 0; there is no relation between the two data series.²The Spearman rank correlation coefficient R (Siegel, 1956) requires that individuals under study be ranked in two ordered series. As the Kendall coefficient t , R expresses the existence of any one relation between two data series if itsvalue is close to 1.²The linear correlation coefficient r (Wonnacott and Wonnacott, 1991) expresses the intensity of a linear relation between two data series; r2 is the part of the variance explained by the linear model.The two first parameters evaluate if a relation exists between observed triazine flux and the different tested maize areas without hypothesis on the form of the relation. The linear correlation coefficient allows a special relation type to be tested. The squared value of the Spearman coefficient R, as the correlation coefficient r, expresses a part of total variance on the ranks. The Kendall coefficient represents the probability of two series being ranked in the same way against the probability of them being ranked in a different way. The use of non-parametric coefficients confers robustness to the method in relation to distributional skewing (Barringer et al., 1990).The most significant correlation levelcorresponds to the most accurate threshold distance d. This distance d defines the zone for which the relation between fields undergoing atrazine treatment and triazine flux is the highest. The buffer of width d will be defined as, “the zone most at risk”, even if plots outside this buffer zone may contribute in a small way to the pollution.3.3. Catchment typologyThe study of the slopes in the whole catchment shows a significant disparity between the upstream and downstream zones. The slopes in the upstream zone are gentle while those in the downstream zone are steep. In order to describe these morphological differences, the index Islope threshold was calculated for each basin: Islope . Sslope.5%=Stotal .1. where Sslope.5% is the basin surface area where the slope is steeper than 5% and Stotal the total surface area of the basin.The 5% threshold slope was chosen because itrepresents the upper limit at which mechanised agriculture can still be practised.The higher the Islope, the greater the proportion of steep slope zones in the basin. In order to sequence basins, a limit of Islope . 0:5 was chosen. This value corresponds to an equal part of flat and steep slope zones in a catchment. Furthermore, this typology separates the sampled basins into two groups of a comparable number of elements. This catchment typology shows a classification according to the position upstream and downstream in the Sousson catchment (Fig. 6).4. ResultsDuring the winter, in December 1997 and March 1998, water quality and discharge measurements were made to determine triazine flux. The network was digitized from the 1/25.000 scale topographic map. The buffers tested are 50 m, 100 m, and 200 m wide. The entire catchment corresponds to the maximumwidth, which is close to 500 m for the downstream group basins and 250 m for the upstream group, which has a more pointed shape. As it is noted by Barringer et al. (1990), the minimum used buffer width must be superior to that of the mapping unit. Here, maize field were determined using information provided by SPOT satellite imagery, (resolution 20 m), with field boundary definition based on 1/10.000 aerial photos (1 mm on the map is equal to 10 m on field).The area was divided into strips around the channel network. Then, the maize fields were putback into this division of space to obtain, for each basin, maize surface area within 50, 100 and 200 m of the hydrographic network, and within the whole catchment.4.1. Study of the whole set of basinsResults of regressions for 23 catchment areas in December 1997 and 26 in March 1998, whichinclude a Kendall rank correlation, a Spearman rank correlation and linear correlation coefficients are given with their significance level in Table 3. Calculated correlation coefficients do not seem to vary consistently as a function of the selected threshold distances: the number of coefficients increase in all cases when the buffer area is enlarged with the exception of December where they decreased in number for the whole catchment area. Considering these results, one might think that the distance of the field from the river has no effect on the pollution. However, if upstream and downstream basins are separated, according to the slope criterion Islope, the results are very different.4.2. Study of the downstream basins Regressions were carried out on nine basins in December 1997 and on 13 in March 1998, mean triazine concentrations are 42.0 and 123.0 ng/l, respectively. Results are shown in Table 4. Calculated correlation coefficients decreasewhen the strip width around the channel network increases. The best correlation levels are obtained for a distance d of 50 m (100 m for the linear correlation in December 1997). The Kendall and the Spearman correlation coefficients show the existence of a relation between maize area inside a 50 m wide buffer zone around the channel network and the triazine flux at the catchment outlet. The linear relation is quite adequate to model this variable association given that 69% of the total variance is explained in December 1997 and 56% in March 1998 considering that d equals 50 m. Resultsobtained for the two measurement dates are mutually coherent although differences exist. In December, whatever the value of d, the significance level is above the acceptance limit (p . 5%). The relation between maize area and triazine flux is optimal for d equal to 50 or 100 m but still exists for d equal to 200 m or considering the whole catchment surface area.The correlation between pollutant flux and maize areas far from the river can be explained by two ways. On the one hand, there is a correlation between the different maize areas (cf. Table 6). Indeed, if maize surface areas within different buffer zones were perfectly proportional, i.e. if linear correlation coefficients between the different maize surfaces areas were equal to one, no variation wouldbe detected in the correlation coefficients between maize surface areas and triazine flux. The sets of basins studied were not exactly the same during the two measurement exercises. For December 1997, the level of correlation between the different maize surface areas is higher than for March 1998 (as it is shown in Table 6). This difference between the two series is partly responsible for the slow decrease in the number of correlation coefficients concerning distance d for December 1997. On the other hand, as it is shown in Fig. 3, December 1997 measurements were made during the falling limb of thehydrograph and thus we can assume that, in these hydrological conditions, the area contributing to pollution is larger and includes zones distant from the hydrographic network for the whole catchment area. However, in March 1998, in lower water level conditions, only correlations where d is equal to 50 m are significant at the 5% threshold.We can conclude that the limit of 50 m is the most appropriate to define the zones most at risk for the two monitoring periods — seven and ten months after the triazine applications —even if hydrological conditions are also important when defining the contribution of the other maize plots located on the whole catchment area.4.3. Study of the upstream basinRegressions were made on 14 catchments for December 1997 and 13 basins for March 1998, mean triazine concentrations are 177.9 and 314.6 ng/ l, respectively. Results are shown inTable 5. The correlation coefficients become more numerous with strip width, while the opposite is true for the downstream basins. In most cases, the best results are obtained by considering the whole catchment area. The linear model is less accurate for the December data set .r2 . 38%. than the Spearman rank correlation .R2 . 70%.: It suggests an association between variables more complex than the linear relation does.Field investigations provide the explanation of the difference between the two catchment groups. For upstream catchments, the hydrographic network taken as the reference is irrelevant. In this flat zone, the artificial drain network around each plot extends the channel network; thus, the real active network is denser than that of the topographic 1/25.000 map. Fig.7 shows, for a particular catchment, the differences between the topographic 1/25.000 map network and the active one observed in the field. Moreover, this ditch network is connectedwith buried drains located under most of the fields in this upstream zone. The consequence is that each field is artificially connected with the catchment outlet.This difference in optimal width between the upstream and downstream catchments is the consequenceof man’s activities on the flat upstream area. In this case, the total catchment surface area must beconsidered as a contributing area.5. DiscussionWe chose to take the measurements in winter because it is easier to compare triazine flux at the catchment outlets. In spring, which is the atrazine spreading period, the differences in flux could be due to differences in the application dates. We used instantaneous inter-storm triazine flux measurements to maximise the stability of the transfer processes (Woods and Sivapalan, 1995). Thus, our resultsdo not necessarily apply to transfer during peak runoff. As the measurements were made between stormy periods our attention was focussed on the slow components of water movement such as subsurface runoff, drainage flow and water circulation in soil, where leaching favours the transport of soluble compounds such as atrazine. These conditions are not maximal from the point of view of instantaneous pollutant quantity export, but do represent a nonnegligible quantity and this over long periods of the year. However, this was a way to acquire comparable data sets at several basin outlets. Moreover, with these data sets it is possible to integrate the spatial diversity and give the results in a form that can be generalised.A simple model of contribution through buffers of stationary width around the hydrographic network was used, where each buffer defines a zone contributing to pollution. The degree of correlation between thecontributing areas and the pollution at the basin outlet was determined.The results show that a critical contribution distance cannot be defined for all basins studied. However, basin typology based on morphology criteria permitted the identification of two groups of basins.These basins have to be considered separately as their water movement characteristics are very different.For the downstream basins, which have a marked relief, the channel is well defined by the network that figures on the 1/25.000 scale topographic map. The model identifies a critical contribution distance, which ranges from 50 to 100 m. Atrazine is little adsorbed by soil, very soluble and easily leached. In inter-storm periods, it is not surface runoff, which causes the water transfers but sub-surface runoff and the draining of local aquifers surrounding the hydrographic network. The area of strongest influence ranges from 50 to 100 m and gives a good representation of the zone where atrazinetransport processes are active. This optimal distance should be determined for different climatic conditions and different periods of the agricultural year. Then we would know if the contributing area possesses temporal dynamics or if it remains stable.The upstream basins have higher triazine concentrations. These areas are characterised by the high proportion of flat zones (slopes of less than 5%), and an artificial drainage network connecting each plot to the main drain in order to avoid flooding. Thus, each plot contributes to the pollution measured at the basin outlet. The topographic 1/25.000 map network does not include this effect of the human intervention on the water circulation, and it is not pertinent in a drained region to evaluate the distance between cropland and the river.How the hydrographic network is defined is critical to the success of this analysis. The initial choice was based on the network digitized from the 1/25.000 scale topographic map. The mainbenefit to be derived from using such a network is its availability, which allows us to easily transpose the methodology. It represents the perennial flow network, stable in time. But, from the point of view of water movement, it lacks locations of manmade drains that can accelerate he transport of solute pollution. From a practical point of view, it is preferable to study the farmlandand identify zones with intensive artificial drainage before defining the boundaries of contributing areas around the channel network.6. ConclusionsIn order to reduce surface water pollution, the application of pesticides has to be controlled and agricultural practices must be such that they respect the environment. But the proper management of cropland must not be neglected either. The spatial organisation of fields has an impact on river pollution.The effect of the distance between fields contributing to the。

User's guide for a three-dimensional, primitive equation, numerical ocean model

2
CONTENTS Page 1. INTRODUCTION 2. THE BASIC EQUATIONS 3. FORTRAN SYMBOLS 4. THE NUMERICAL SCHEME 5. pom2k.c 6. program main and the external mode 7. subroutine advave 8. subroutine advt 9. subroutine proft 10. subroutine baropg 11. subroutines advct, advu and advv 12. subroutines profu and profv 13. subroutine advq 14. subroutine profq 15. subroutine vertvl 16. subroutine bcond 17. subroutine dens 18 subroutine slpmin 19. Utility Subroutines 20. program curvigrid 5 7 14 17 23 23 24 24 24 27 27 28 28 28 29 29 34 33 34 35
This revision: June 2004
Notes on the 1998 Revision This version of the users guide recognizes changes that have occurred since 1991. The code itself incorporates some recent changes. the fortran names, tmean, smean have been changed (globally) to tclim, sclim in oder to distiquish the function and treatment of these variables from that of rmean. the names, trnu, trnv, have been changed to drx2d, dry2d and the names, advuu, advvv, to adx2d, ady2d to more clearly indicate their functions. Instead of a wind driven closed basin, pom97.f now solves the problem of the flow through a channel which includes an island or a seamount at the center of the domain. Thus, subroutine bcond contains active open boundary conditions. These illustrative boundary conditions, however, are one set of many possibilities and, consequently, open boundary conditions for regional models pose difficult choices for users of the model. This 1998 revision contains a fuller discussion of open boundary conditions in section 16. Notes on the 2002 revision The basic code, now labeled pom2k.f results from extensive tidying by John Hunter which includes more comments and lower case fortran variables, a move which apparently renders the code “modern”. However the basic – we believe, well conceived - structure of the code remains unchanged. As of this revision date, June 2004, there are over 1900 POM users of record.

渗透系数的英语

渗透系数的英语The concept of permeability is a fundamental aspect of fluid mechanics and is crucial in various fields, including civil engineering, petroleum engineering, and environmental science. The permeability coefficient, also known as the coefficient of permeability, is a quantitative measure of a material's ability to allow the flow of fluids through its porous structure. This parameter is essential in understanding and predicting the behavior of fluids in porous media, such as soil, rock, and other materials.Permeability is a measure of the ease with which a fluid can flow through a porous medium. It is influenced by the size, shape, and interconnectivity of the pores within the material. The permeability coefficient is a numerical value that represents the ease of fluid flow through a specific porous medium under a given set of conditions. This coefficient is typically denoted by the symbol "k" and is expressed in units of area, such as square meters (m²) or darcies (D).The permeability coefficient is not a constant value and can vary depending on several factors, including the properties of the porousmedium, the properties of the fluid, and the flow conditions. The primary factors that affect the permeability coefficient are the porosity, tortuosity, and pore size distribution of the material.Porosity is the ratio of the volume of voids or pore spaces to the total volume of the material. Materials with higher porosity generally have a higher permeability coefficient, as they offer less resistance to fluid flow. Tortuosity, on the other hand, is a measure of the complexity of the flow paths within the porous medium. Materials with a higher tortuosity have a lower permeability coefficient, as the fluid must navigate through a more convoluted path.The pore size distribution also plays a crucial role in determining the permeability coefficient. Materials with larger and more interconnected pores typically have a higher permeability coefficient, as the fluid can flow more easily through the porous structure. Conversely, materials with smaller and less interconnected pores have a lower permeability coefficient, as the fluid encounters greater resistance to flow.In addition to the properties of the porous medium, the properties of the fluid, such as viscosity and density, can also affect the permeability coefficient. Fluids with lower viscosity, such as water, generally have a higher permeability coefficient compared to fluids with higher viscosity, such as oil or honey.The measurement of the permeability coefficient is an essential aspect of understanding and predicting the behavior of fluids in porous media. There are several methods used to determine the permeability coefficient, including laboratory experiments, field measurements, and numerical simulations.One of the most common laboratory methods for measuring the permeability coefficient is the constant-head or falling-head permeameter test. In this test, a sample of the porous material is placed in a permeameter, and a constant or falling head of fluid is applied across the sample. The rate of fluid flow through the sample is then measured, and the permeability coefficient is calculated using Darcy's law, which relates the fluid flow rate to the pressure drop across the sample.Another method for measuring the permeability coefficient is the use of numerical simulations, such as computational fluid dynamics (CFD) or pore-scale modeling. These techniques involve the use of computer models to simulate the flow of fluids through the porous medium, taking into account the complex geometry and topology of the porous structure. By comparing the simulated flow rates to experimental data, the permeability coefficient can be estimated.The permeability coefficient is a crucial parameter in variousapplications, including:1. Soil mechanics and geotechnical engineering: The permeability coefficient is used to determine the rate of groundwater flow, the stability of soil slopes, and the design of drainage systems.2. Petroleum engineering: The permeability coefficient is essential in understanding the flow of oil and gas through reservoir rocks, which is crucial for the exploration, production, and management of hydrocarbon resources.3. Environmental engineering: The permeability coefficient is used to model the transport of contaminants through soil and groundwater, which is important for the design of waste disposal facilities and the remediation of contaminated sites.4. Civil engineering: The permeability coefficient is used in the design of concrete structures, as it influences the durability and performance of the material under various environmental conditions.5. Materials science: The permeability coefficient is studied in the context of porous materials, such as ceramics, membranes, and filters, to understand their ability to allow the flow of fluids or gases.In conclusion, the permeability coefficient is a fundamentalparameter in the study of fluid flow through porous media. It is influenced by the properties of the porous medium and the fluid, and its measurement and understanding are crucial in various fields of engineering and science. The accurate determination and application of the permeability coefficient are essential for the design, analysis, and optimization of systems and processes that involve the flow of fluids through porous materials.。

专八英语阅读

英语专业八级考试TEM-8阅读理解练习册（1）(英语专业2012级)UNIT 1Text AEvery minute of every day, what ecologist生态学家James Carlton calls a global ―conveyor belt‖, redistributes ocean organisms生物.It’s planetwide biological disruption生物的破坏that scientists have barely begun to understand.Dr. Carlton —an oceanographer at Williams College in Williamstown,Mass.—explains that, at any given moment, ―There are several thousand marine species traveling… in the ballast water of ships.‖ These creatures move from coastal waters where they fit into the local web of life to places where some of them could tear that web apart. This is the larger dimension of the infamous无耻的，邪恶的invasion of fish-destroying, pipe-clogging zebra mussels有斑马纹的贻贝.Such voracious贪婪的invaders at least make their presence known. What concerns Carlton and his fellow marine ecologists is the lack of knowledge about the hundreds of alien invaders that quietly enter coastal waters around the world every day. Many of them probably just die out. Some benignly亲切地，仁慈地—or even beneficially — join the local scene. But some will make trouble.In one sense, this is an old story. Organisms have ridden ships for centuries. They have clung to hulls and come along with cargo. What’s new is the scale and speed of the migrations made possible by the massive volume of ship-ballast water压载水— taken in to provide ship stability—continuously moving around the world…Ships load up with ballast water and its inhabitants in coastal waters of one port and dump the ballast in another port that may be thousands of kilometers away. A single load can run to hundreds of gallons. Some larger ships take on as much as 40 million gallons. The creatures that come along tend to be in their larva free-floating stage. When discharged排出in alien waters they can mature into crabs, jellyfish水母, slugs鼻涕虫，蛞蝓, and many other forms.Since the problem involves coastal species, simply banning ballast dumps in coastal waters would, in theory, solve it. Coastal organisms in ballast water that is flushed into midocean would not survive. Such a ban has worked for North American Inland Waterway. But it would be hard to enforce it worldwide. Heating ballast water or straining it should also halt the species spread. But before any such worldwide regulations were imposed, scientists would need a clearer view of what is going on.The continuous shuffling洗牌of marine organisms has changed the biology of the sea on a global scale. It can have devastating effects as in the case of the American comb jellyfish that recently invaded the Black Sea. It has destroyed that sea’s anchovy鳀鱼fishery by eating anchovy eggs. It may soon spread to western and northern European waters.The maritime nations that created the biological ―conveyor belt‖ should support a coordinated international effort to find out what is going on and what should be done about it. (456 words)1.According to Dr. Carlton, ocean organism‟s are_______.A.being moved to new environmentsB.destroying the planetC.succumbing to the zebra musselD.developing alien characteristics2.Oceanographers海洋学家are concerned because_________.A.their knowledge of this phenomenon is limitedB.they believe the oceans are dyingC.they fear an invasion from outer-spaceD.they have identified thousands of alien webs3.According to marine ecologists, transplanted marinespecies____________.A.may upset the ecosystems of coastal watersB.are all compatible with one anotherC.can only survive in their home watersD.sometimes disrupt shipping lanes4.The identified cause of the problem is_______.A.the rapidity with which larvae matureB. a common practice of the shipping industryC. a centuries old speciesD.the world wide movement of ocean currents5.The article suggests that a solution to the problem__________.A.is unlikely to be identifiedB.must precede further researchC.is hypothetically假设地，假想地easyD.will limit global shippingText BNew …Endangered‟ List Targets Many US RiversIt is hard to think of a major natural resource or pollution issue in North America today that does not affect rivers.Farm chemical runoff残渣, industrial waste, urban storm sewers, sewage treatment, mining, logging, grazing放牧,military bases, residential and business development, hydropower水力发电,loss of wetlands. The list goes on.Legislation like the Clean Water Act and Wild and Scenic Rivers Act have provided some protection, but threats continue.The Environmental Protection Agency (EPA) reported yesterday that an assessment of 642,000 miles of rivers and streams showed 34 percent in less than good condition. In a major study of the Clean Water Act, the Natural Resources Defense Council last fall reported that poison runoff impairs损害more than 125,000 miles of rivers.More recently, the NRDC and Izaak Walton League warned that pollution and loss of wetlands—made worse by last year’s flooding—is degrading恶化the Mississippi River ecosystem.On Tuesday, the conservation group保护组织American Rivers issued its annual list of 10 ―endangered‖ and 20 ―threatened‖ rivers in 32 states, the District of Colombia, and Canada.At the top of the list is the Clarks Fork of the Yellowstone River, whereCanadian mining firms plan to build a 74-acre英亩reservoir水库，蓄水池as part of a gold mine less than three miles from Yellowstone National Park. The reservoir would hold the runoff from the sulfuric acid 硫酸used to extract gold from crushed rock.―In the event this tailings pond failed, the impact to th e greater Yellowstone ecosystem would be cataclysmic大变动的，灾难性的and the damage irreversible不可逆转的.‖ Sen. Max Baucus of Montana, chairman of the Environment and Public Works Committee, wrote to Noranda Minerals Inc., an owner of the ― New World Mine‖.Last fall, an EPA official expressed concern about the mine and its potential impact, especially the plastic-lined storage reservoir. ― I am unaware of any studies evaluating how a tailings pond尾矿池，残渣池could be maintained to ensure its structural integrity forev er,‖ said Stephen Hoffman, chief of the EPA’s Mining Waste Section. ―It is my opinion that underwater disposal of tailings at New World may present a potentially significant threat to human health and the environment.‖The results of an environmental-impact statement, now being drafted by the Forest Service and Montana Department of State Lands, could determine the mine’s future…In its recent proposal to reauthorize the Clean Water Act, the Clinton administration noted ―dramatically improved water quality since 1972,‖ when the act was passed. But it also reported that 30 percent of riverscontinue to be degraded, mainly by silt泥沙and nutrients from farm and urban runoff, combined sewer overflows, and municipal sewage城市污水. Bottom sediments沉积物are contaminated污染in more than 1,000 waterways, the administration reported in releasing its proposal in January. Between 60 and 80 percent of riparian corridors (riverbank lands) have been degraded.As with endangered species and their habitats in forests and deserts, the complexity of ecosystems is seen in rivers and the effects of development----beyond the obvious threats of industrial pollution, municipal waste, and in-stream diversions改道to slake消除the thirst of new communities in dry regions like the Southwes t…While there are many political hurdles障碍ahead, reauthorization of the Clean Water Act this year holds promise for US rivers. Rep. Norm Mineta of California, who chairs the House Committee overseeing the bill, calls it ―probably the most important env ironmental legislation this Congress will enact.‖ (553 words)6.According to the passage, the Clean Water Act______.A.has been ineffectiveB.will definitely be renewedC.has never been evaluatedD.was enacted some 30 years ago7.“Endangered” rivers are _________.A.catalogued annuallyB.less polluted than ―threatened rivers‖C.caused by floodingD.adjacent to large cities8.The “cataclysmic” event referred to in paragraph eight would be__________.A. fortuitous偶然的，意外的B. adventitious外加的，偶然的C. catastrophicD. precarious不稳定的，危险的9. The owners of the New World Mine appear to be______.A. ecologically aware of the impact of miningB. determined to construct a safe tailings pondC. indifferent to the concerns voiced by the EPAD. willing to relocate operations10. The passage conveys the impression that_______.A. Canadians are disinterested in natural resourcesB. private and public environmental groups aboundC. river banks are erodingD. the majority of US rivers are in poor conditionText CA classic series of experiments to determine the effects ofoverpopulation on communities of rats was reported in February of 1962 in an article in Scientific American. The experiments were conducted by a psychologist, John B. Calhoun and his associates. In each of these experiments, an equal number of male and female adult rats were placed in an enclosure and given an adequate supply of food, water, and other necessities. The rat populations were allowed to increase. Calhoun knew from experience approximately how many rats could live in the enclosures without experiencing stress due to overcrowding. He allowed the population to increase to approximately twice this number. Then he stabilized the population by removing offspring that were not dependent on their mothers. He and his associates then carefully observed and recorded behavior in these overpopulated communities. At the end of their experiments, Calhoun and his associates were able to conclude that overcrowding causes a breakdown in the normal social relationships among rats, a kind of social disease. The rats in the experiments did not follow the same patterns of behavior as rats would in a community without overcrowding.The females in the rat population were the most seriously affected by the high population density: They showed deviant异常的maternal behavior; they did not behave as mother rats normally do. In fact, many of the pups幼兽，幼崽, as rat babies are called, died as a result of poor maternal care. For example, mothers sometimes abandoned their pups,and, without their mothers' care, the pups died. Under normal conditions, a mother rat would not leave her pups alone to die. However, the experiments verified that in overpopulated communities, mother rats do not behave normally. Their behavior may be considered pathologically 病理上，病理学地diseased.The dominant males in the rat population were the least affected by overpopulation. Each of these strong males claimed an area of the enclosure as his own. Therefore, these individuals did not experience the overcrowding in the same way as the other rats did. The fact that the dominant males had adequate space in which to live may explain why they were not as seriously affected by overpopulation as the other rats. However, dominant males did behave pathologically at times. Their antisocial behavior consisted of attacks on weaker male,female, and immature rats. This deviant behavior showed that even though the dominant males had enough living space, they too were affected by the general overcrowding in the enclosure.Non-dominant males in the experimental rat communities also exhibited deviant social behavior. Some withdrew completely; they moved very little and ate and drank at times when the other rats were sleeping in order to avoid contact with them. Other non-dominant males were hyperactive; they were much more active than is normal, chasing other rats and fighting each other. This segment of the rat population, likeall the other parts, was affected by the overpopulation.The behavior of the non-dominant males and of the other components of the rat population has parallels in human behavior. People in densely populated areas exhibit deviant behavior similar to that of the rats in Calhoun's experiments. In large urban areas such as New York City, London, Mexican City, and Cairo, there are abandoned children. There are cruel, powerful individuals, both men and women. There are also people who withdraw and people who become hyperactive. The quantity of other forms of social pathology such as murder, rape, and robbery also frequently occur in densely populated human communities. Is the principal cause of these disorders overpopulation? Calhoun’s experiments suggest that it might be. In any case, social scientists and city planners have been influenced by the results of this series of experiments.11. Paragraph l is organized according to__________.A. reasonsB. descriptionC. examplesD. definition12．Calhoun stabilized the rat population_________.A. when it was double the number that could live in the enclosure without stressB. by removing young ratsC. at a constant number of adult rats in the enclosureD. all of the above are correct13.W hich of the following inferences CANNOT be made from theinformation inPara. 1?A. Calhoun's experiment is still considered important today.B. Overpopulation causes pathological behavior in rat populations.C. Stress does not occur in rat communities unless there is overcrowding.D. Calhoun had experimented with rats before.14. Which of the following behavior didn‟t happen in this experiment?A. All the male rats exhibited pathological behavior.B. Mother rats abandoned their pups.C. Female rats showed deviant maternal behavior.D. Mother rats left their rat babies alone.15. The main idea of the paragraph three is that __________.A. dominant males had adequate living spaceB. dominant males were not as seriously affected by overcrowding as the otherratsC. dominant males attacked weaker ratsD. the strongest males are always able to adapt to bad conditionsText DThe first mention of slavery in the statutes法令，法规of the English colonies of North America does not occur until after 1660—some forty years after the importation of the first Black people. Lest we think that existed in fact before it did in law, Oscar and Mary Handlin assure us, that the status of B lack people down to the 1660’s was that of servants. A critique批判of the Handlins’ interpretation of why legal slavery did not appear until the 1660’s suggests that assumptions about the relation between slavery and racial prejudice should be reexamined, and that explanation for the different treatment of Black slaves in North and South America should be expanded.The Handlins explain the appearance of legal slavery by arguing that, during the 1660’s, the position of white servants was improving relative to that of black servants. Thus, the Handlins contend, Black and White servants, heretofore treated alike, each attained a different status. There are, however, important objections to this argument. First, the Handlins cannot adequately demonstrate that t he White servant’s position was improving, during and after the 1660’s; several acts of the Maryland and Virginia legislatures indicate otherwise. Another flaw in the Handlins’ interpretation is their assumption that prior to the establishment of legal slavery there was no discrimination against Black people. It is true that before the 1660’s Black people were rarely called slaves. But this shouldnot overshadow evidence from the 1630’s on that points to racial discrimination without using the term slavery. Such discrimination sometimes stopped short of lifetime servitude or inherited status—the two attributes of true slavery—yet in other cases it included both. The Handlins’ argument excludes the real possibility that Black people in the English colonies were never treated as the equals of White people.The possibility has important ramifications后果，影响.If from the outset Black people were discriminated against, then legal slavery should be viewed as a reflection and an extension of racial prejudice rather than, as many historians including the Handlins have argued, the cause of prejudice. In addition, the existence of discrimination before the advent of legal slavery offers a further explanation for the harsher treatment of Black slaves in North than in South America. Freyre and Tannenbaum have rightly argued that the lack of certain traditions in North America—such as a Roman conception of slavery and a Roman Catholic emphasis on equality— explains why the treatment of Black slaves was more severe there than in the Spanish and Portuguese colonies of South America. But this cannot be the whole explanation since it is merely negative, based only on a lack of something. A more compelling令人信服的explanation is that the early and sometimes extreme racial discrimination in the English colonies helped determine the particular nature of the slavery that followed. (462 words)16. Which of the following is the most logical inference to be drawn from the passage about the effects of “several acts of the Maryland and Virginia legislatures” (Para.2) passed during and after the 1660‟s?A. The acts negatively affected the pre-1660’s position of Black as wellas of White servants.B. The acts had the effect of impairing rather than improving theposition of White servants relative to what it had been before the 1660’s.C. The acts had a different effect on the position of white servants thandid many of the acts passed during this time by the legislatures of other colonies.D. The acts, at the very least, caused the position of White servants toremain no better than it had been before the 1660’s.17. With which of the following statements regarding the status ofBlack people in the English colonies of North America before the 1660‟s would the author be LEAST likely to agree?A. Although black people were not legally considered to be slaves,they were often called slaves.B. Although subject to some discrimination, black people had a higherlegal status than they did after the 1660’s.C. Although sometimes subject to lifetime servitude, black peoplewere not legally considered to be slaves.D. Although often not treated the same as White people, black people,like many white people, possessed the legal status of servants.18. According to the passage, the Handlins have argued which of thefollowing about the relationship between racial prejudice and the institution of legal slavery in the English colonies of North America?A. Racial prejudice and the institution of slavery arose simultaneously.B. Racial prejudice most often the form of the imposition of inheritedstatus, one of the attributes of slavery.C. The source of racial prejudice was the institution of slavery.D. Because of the influence of the Roman Catholic Church, racialprejudice sometimes did not result in slavery.19. The passage suggests that the existence of a Roman conception ofslavery in Spanish and Portuguese colonies had the effect of _________.A. extending rather than causing racial prejudice in these coloniesB. hastening the legalization of slavery in these colonies.C. mitigating some of the conditions of slavery for black people in these coloniesD. delaying the introduction of slavery into the English colonies20. The author considers the explanation put forward by Freyre andTannenbaum for the treatment accorded B lack slaves in the English colonies of North America to be _____________.A. ambitious but misguidedB. valid有根据的but limitedC. popular but suspectD. anachronistic过时的，时代错误的and controversialUNIT 2Text AThe sea lay like an unbroken mirror all around the pine-girt, lonely shores of Orr’s Island. Tall, kingly spruce s wore their regal王室的crowns of cones high in air, sparkling with diamonds of clear exuded gum流出的树胶; vast old hemlocks铁杉of primeval原始的growth stood darkling in their forest shadows, their branches hung with long hoary moss久远的青苔;while feathery larches羽毛般的落叶松,turned to brilliant gold by autumn frosts, lighted up the darker shadows of the evergreens. It was one of those hazy朦胧的, calm, dissolving days of Indian summer, when everything is so quiet that the fainest kiss of the wave on the beach can be heard, and white clouds seem to faint into the blue of the sky, and soft swathing一长条bands of violet vapor make all earth look dreamy, and give to the sharp, clear-cut outlines of the northern landscape all those mysteries of light and shade which impart such tenderness to Italian scenery.The funeral was over,--- the tread鞋底的花纹/ 踏of many feet, bearing the heavy burden of two broken lives, had been to the lonely graveyard, and had come back again,--- each footstep lighter and more unconstrained不受拘束的as each one went his way from the great old tragedy of Death to the common cheerful of Life.The solemn black clock stood swaying with its eternal ―tick-tock, tick-tock,‖ in the kitchen of the brown house on Orr’s Island. There was there that sense of a stillness that can be felt,---such as settles down on a dwelling住处when any of its inmates have passed through its doors for the last time, to go whence they shall not return. The best room was shut up and darkened, with only so much light as could fall through a little heart-shaped hole in the window-shutter,---for except on solemn visits, or prayer-meetings or weddings, or funerals, that room formed no part of the daily family scenery.The kitchen was clean and ample, hearth灶台, and oven on one side, and rows of old-fashioned splint-bottomed chairs against the wall. A table scoured to snowy whiteness, and a little work-stand whereon lay the Bible, the Missionary Herald, and the Weekly Christian Mirror, before named, formed the principal furniture. One feature, however, must not be forgotten, ---a great sea-chest水手用的储物箱,which had been the companion of Zephaniah through all the countries of the earth. Old, and battered破旧的，磨损的, and unsightly难看的it looked, yet report said that there was good store within which men for the most part respect more than anything else; and, indeed it proved often when a deed of grace was to be done--- when a woman was suddenly made a widow in a coast gale大风，狂风, or a fishing-smack小渔船was run down in the fogs off the banks, leaving in some neighboring cottage a family of orphans,---in all such cases, the opening of this sea-chest was an event of good omen 预兆to the bereaved丧亲者;for Zephaniah had a large heart and a large hand, and was apt有…的倾向to take it out full of silver dollars when once it went in. So the ark of the covenant约柜could not have been looked on with more reverence崇敬than the neighbours usually showed to Captain Pennel’s sea-chest.1. The author describes Orr‟s Island in a(n)______way.A.emotionally appealing, imaginativeB.rational, logically preciseC.factually detailed, objectiveD.vague, uncertain2.According to the passage, the “best room”_____.A.has its many windows boarded upB.has had the furniture removedC.is used only on formal and ceremonious occasionsD.is the busiest room in the house3.From the description of the kitchen we can infer that thehouse belongs to people who_____.A.never have guestsB.like modern appliancesC.are probably religiousD.dislike housework4.The passage implies that_______.A.few people attended the funeralB.fishing is a secure vocationC.the island is densely populatedD.the house belonged to the deceased5.From the description of Zephaniah we can see thathe_________.A.was physically a very big manB.preferred the lonely life of a sailorC.always stayed at homeD.was frugal and saved a lotText BBasic to any understanding of Canada in the 20 years after the Second World War is the country' s impressive population growth. For every three Canadians in 1945, there were over five in 1966. In September 1966 Canada's population passed the 20 million mark. Most of this surging growth came from natural increase. The depression of the 1930s and the war had held back marriages, and the catching-up process began after 1945. The baby boom continued through the decade of the 1950s, producing a population increase of nearly fifteen percent in the five years from 1951 to 1956. This rate of increase had been exceeded only once before in Canada's history, in the decade before 1911 when the prairies were being settled. Undoubtedly, the good economic conditions of the 1950s supported a growth in the population, but the expansion also derived from a trend toward earlier marriages and an increase in the average size of families; In 1957 the Canadian birth rate stood at 28 per thousand, one of the highest in the world. After the peak year of 1957, thebirth rate in Canada began to decline. It continued falling until in 1966 it stood at the lowest level in 25 years. Partly this decline reflected the low level of births during the depression and the war, but it was also caused by changes in Canadian society. Young people were staying at school longer, more women were working; young married couples were buying automobiles or houses before starting families; rising living standards were cutting down the size of families. It appeared that Canada was once more falling in step with the trend toward smaller families that had occurred all through theWestern world since the time of the Industrial Revolution. Although the growth in Canada’s population had slowed down by 1966 (the cent), another increase in the first half of the 1960s was only nine percent), another large population wave was coming over the horizon. It would be composed of the children of the children who were born during the period of the high birth rate prior to 1957.6. What does the passage mainly discuss?A. Educational changes in Canadian society.B. Canada during the Second World War.C. Population trends in postwar Canada.D. Standards of living in Canada.7. According to the passage, when did Canada's baby boom begin?A. In the decade after 1911.B. After 1945.C. During the depression of the 1930s.D. In 1966.8. The author suggests that in Canada during the 1950s____________.A. the urban population decreased rapidlyB. fewer people marriedC. economic conditions were poorD. the birth rate was very high9. When was the birth rate in Canada at its lowest postwar level?A. 1966.B. 1957.C. 1956.D. 1951.10. The author mentions all of the following as causes of declines inpopulation growth after 1957 EXCEPT_________________.A. people being better educatedB. people getting married earlierC. better standards of livingD. couples buying houses11.I t can be inferred from the passage that before the IndustrialRevolution_______________.A. families were largerB. population statistics were unreliableC. the population grew steadilyD. economic conditions were badText CI was just a boy when my father brought me to Harlem for the first time, almost 50 years ago. We stayed at the hotel Theresa, a grand brick structure at 125th Street and Seventh avenue. Once, in the hotel restaurant, my father pointed out Joe Louis. He even got Mr. Brown, the hotel manager, to introduce me to him, a bit punchy强力的but still champ焦急as fast as I was concerned.Much has changed since then. Business and real estate are booming. Some say a new renaissance is under way. Others decry责难what they see as outside forces running roughshod肆意践踏over the old Harlem. New York meant Harlem to me, and as a young man I visited it whenever I could. But many of my old haunts are gone. The Theresa shut down in 1966. National chains that once ignored Harlem now anticipate yuppie money and want pieces of this prime Manhattan real estate. So here I am on a hot August afternoon, sitting in a Starbucks that two years ago opened a block away from the Theresa, snatching抓取，攫取at memories between sips of high-priced coffee. I am about to open up a piece of the old Harlem---the New York Amsterdam News---when a tourist。

利用井眼坍塌信息求取最大水平主应力的方法

利用井眼坍塌信息求取最大水平主应力的方法李光泉;王怡;陈军海;陈曾伟【摘要】准确获取深部地层地应力的方向和大小,对于井眼轨道设计、井壁稳定、压裂增产等具有重要意义.重点研究了利用井眼坍塌信息求取最大水平主应力的方法,提出结合测井资料分析采用数值模拟方法反演最大水平主应力大小的方法.研究表明:利用井眼坍塌信息求取最大水平主应力大小的数值模拟方法可行,能够真实反映深部地层的应力环境、地层岩石的力学特性,比较准确地确定最大水平主应力的大小,且方法简单易行.%Obtaining the maximum horizontal principal stress accurately is a difficult problem for deep well in-situ stess. Analysis of the maximum horizontal principal stress based on borehole collapse information is discussed in this paper. Numerical simulation analysis with log data and field tests to determine the maximum horizontal principal stress is proposed. The results show that,it's feasible using numerical simulation method to invert in-situ stress based on borehole collapse information; this method can simulate the stress environment of deep stratum, the mechanical properties of rock formations and drilling effect, which can help determine the magnitude of the maximum horizontal principal stress accurately and simply.【期刊名称】《石油钻探技术》【年(卷),期】2012(040)001【总页数】5页(P37-41)【关键词】井眼坍塌;水平地应力;反演分析;数值模拟【作者】李光泉;王怡;陈军海;陈曾伟【作者单位】中国石化石油工程技术研究院,北京100101;中国石化石油工程技术研究院,北京100101;中国石化石油工程技术研究院,北京100101;中国石化石油工程技术研究院,北京100101【正文语种】中文【中图分类】TE21原场地应力的大小和方向对于井壁稳定、压裂增产等意义重大。

页岩气多尺度渗流数值模拟技术——以昭通国家级页岩气示范区为例

象的模拟，要在宏观模拟中表征上述特殊渗流现象，合理的做法是引入等效数学模型，如气体状态方程、渗透率模型及气体吸附模型，使流体和岩石的物理性质随可测试到的热力学参数（压力、温度、吸附浓度）发生变化，从而影响孔隙容积和流体相的流动性。在本文中，公式内的各种物理量均采用国际单位制（SI）进行说明。
A numerical simulation technology for the multi-scale flow of shale gas and its application in Zhaotong National Shale Gas Demonstration Area
ZHANG Zhuo, YUAN Xiaojun, RAO Daqian, SHU Honglin, YIN Kaigui
第 41 卷增刊 1 2021 年 3 月
天然气工业 Natural Gas Industry
· 145 ·

页岩气多尺度渗流数值模拟技术
——以昭通国家级页岩气示范区为例
张卓袁晓俊饶大骞舒红林尹开贵
中国石油浙江油田公司
摘要：页岩气藏的渗流机理复杂，沿用传统的油藏数值模拟器表现出了不适应性。为了实现对页岩气藏的有效模拟，引入等效数学模型——气体状态方程、渗透率模型及气体吸附模型，并且建立嵌入式离散裂缝网格剖分方法，进而开发出流固耦合模型；在此基础上，基于昭通国家级页岩气示范区 2 口水平井的 Petrel 地质建模成果，利用所建立的流固耦合模型进行压裂后气井的生产历史拟合和预测，进而开展了参数敏感性分析。研究结果表明：①流固耦合模型可以考虑孔隙压实和裂缝变形的影响，功能扩展后的数值模拟器可以更准确地模拟致密介质中的页岩气渗流特征；②嵌入式离散裂缝网格剖分方法能够有效提高建模效率和数值模拟计算速度，支撑了页岩气井生产数据的高效模拟；③随着基质渗透率、裂缝渗透率、裂缝长度逐渐增大，页岩气井累计产气量逐渐升高，但增幅略有变小；④随着应力敏感系数的逐渐增大，单井累计产气量逐渐降低，但降幅逐渐变小。结论认为，该数值模拟技术可以应用于页岩气井的生产动态分析，可以为同类气藏的开发提供借鉴。关键词：页岩气；多尺度渗流；数值模拟；嵌入式裂缝；流固耦合模型；昭通国家级页岩气示范区 DOI ：10.3787/j.issn.1000-0976.2021.S1.021

多孔材料的孔结构对其力学性能及破裂机制的影响

ａｕｔｈｏｒ: ＪＩＡＰｅｎｇꎬ Ｅ￣ｍａｉｌ: ｐｏｌｏｒｊｉａ＠１６３. ｃｏｍ)
Ａｂｓｔｒａｃｔ: Ｉｎｏｒｄｅｒｔｏｓｔｕｄｙｔｈｅｅｆｆｅｃｔｏｆｐｏｒｅｓｔｒｕｃｔｕｒｅｏｎｔｈｅｍｅｃｈａｎｉｃａｌｐｒｏｐｅｒｔｉｅｓꎬ ｆｒａｃｔｕｒｅ
ｍｅｃｈａｎｉｓｍａｎｄｄｕｒａｂｉｌｉｔｙｏｆａｌｔｅｒｅｄｒｏｃｋｐｏｒｏｕｓｍａｔｅｒｉａｌｓꎬ ｔｈｅｕｎｉａｘｉａｌｃｏｍｐｒｅｓｓｉｏｎａｎｄｆｒｅｅｚｅ￣
ｌａｒｇｅｒｔｈｅｐｏｒｏｓｉｔｙꎬ ｔｈｅｌｏｗｅｒｔｈｅｃｏｍｐｒｅｓｓｉｖｅｓｔｒｅｎｇｔｈａｎｄｔｈｅｈｉｇｈｅｒｔｈｅｅｎｅｒｇｙａｂｓｏｒｐｔｉｏｎ. Ｔｈｅ
１０￣ｃｙｃｌｅｆｒｅｅｚｅ￣ｔｈａｗｅｘｐｅｒｉｍｅｎｔｓｉｎｔｈｅｒａｎｇｅｏｆ ± ２５ ℃ ｓｈｏｗｔｈａｔｗｈｅｎｔｈｅｐｏｒｏｓｉｔｙｉｓｓｍａｌｌａｎｄ
等 [１９] 在研究冻融循环对泡沫混凝土的强度影响
时发现ꎬ经过冻融循环后ꎬ泡沫混凝土的强度降低
约１５％ . Ｔｉｋａｌｓｋｙ等 [２０] 在对泡沫混凝土进行冻融
图１尾矿陶瓷材料
Ｆｉｇ１Ｔａｉｌｉｎｇｓｃｅｒａｍｉｃｍａｔｅｒｉａｌｓ
循环试验后发现ꎬ泡沫混凝土孔径增加ꎬ其抵抗冻
融循环能力增强ꎬ分析其原因在于气孔为水转化成
Ｔａｂｌｅ１Ｐｈｙｓｉｃａｌｐａｒａｍｅｔｅｒｓｏｆｔｈｅｔｗｏｔｙｐｅｓｏｆｐｏｒｏｕｓｍａｔｅｒｉａｌｓ
名称
质量
ｇ
饱水质量
ｇ
体积
ｃｍ３
干密度
ｇｃｍ－３
饱和密度
ｇｃｍ－３

环境水侵蚀下水泥净浆钙溶蚀的模拟与验证

环境水侵蚀下水泥净浆钙溶蚀的模拟与验证马强;左晓宝;汤玉娟【摘要】The durability degradation caused by calcium leaching is very common in the hydraulic concrete structures. In order to understand the calcium leaching process of the cement paste under the attack of the environmental water, at first, a transport model for calcium ion in the cement paste slice specimen was established by utilization of Fick's laws, the mass conservation law, the solid?liquid equilibrium relationships between calcium ion concentration in pore solution and calcium content in the solid skeleton of cement paste, as well as the Newton law on the boundary of slice specimen immersed into the environmental water. And then, an accelerated calcium leaching experiment of cement paste slice specimens with different water?cement ratios immersed into 6M NH4 Cl solution was carried out to measure the Ca/Si and average porosity of the slice specimens at different leaching time, and the Ca/Si and average porosity calculated by the established model were further compared with that from the accelerated experiments to verify the proposed model. Finally, numerical simulations were performed to analyze the space?time changes of the calcium ion concentration in the pore solution, the solid calcium content in skeleton and the porosity in the cement paste slice specimen. Simulated results show that the results calculated by the model are essentially in agreement with the experimental results;in the early stages of the calcium leaching, the solid calcium content in the cementpaste specimens decreases fast, and its porosity has also a rapid increase, but in the later period of leaching, the leaching speed of solid calcium and the increase rate of the porosity gradually decrease.%钙溶蚀是导致水环境中混凝土等水泥基材料耐久性退化的重要原因之一.为获得软水环境下水泥净浆的钙溶蚀过程,首先,基于Fick定律及质量守恒定律,利用钙溶蚀过程中材料骨架内固体钙含量和孔溶液中钙离子浓度之间的化学平衡关系及Newton边界条件,建立软水环境下水泥净浆的钙溶蚀模型,并通过有限差分法,对该模型进行数值求解;其次,进行不同水灰比的水泥净浆试件在6M NH4 Cl溶液中的加速钙溶蚀试验,测定该溶液中各水泥净浆试件在不同溶蚀时间的钙硅比与孔隙率,并将所建立模型的计算结果与实测结果进行对比分析,验证模型的合理性;最后,利用验证后的钙溶蚀模型,数值分析了环境水侵蚀下水泥净浆薄板孔溶液中钙离子浓度、固体钙含量及孔隙率的时空分布规律.结果表明,模型的计算结果与试验测试结果基本一致;溶蚀前期,试件中固体钙含量下降速度和孔隙率增加速率均较大,溶蚀后期,试件固体钙溶蚀速率和孔隙率的增加速率逐渐减小.【期刊名称】《水利水运工程学报》【年(卷),期】2017(000)003【总页数】9页(P107-115)【关键词】钙溶蚀;水泥净浆;传输模型;软水侵蚀;数值模拟【作者】马强;左晓宝;汤玉娟【作者单位】南京理工大学理学院土木工程系, 江苏南京 210094;南京理工大学理学院土木工程系, 江苏南京 210094;南京理工大学理学院土木工程系, 江苏南京 210094【正文语种】中文【中图分类】TV432;TU528.01长期处于环境水中的水坝、港口、水槽、桥梁等水利工程中的混凝土结构，钙溶蚀是最常见的病害之一[1-2]。

应用纳米压痕技术测试水泥粒子的弹性模量

(1. 东南大学材料科学与工程系 ,南京 210096 ;2. 燕山大学建筑工程与力学学院 , 秦皇岛 066004)
摘要 : 应用纳米压痕技术实测了水泥粒子的弹性模量 ,并对试样制备技术进行了研究 ,分析了合理的样本数量。研
究表明 :应用酚醛树脂镶嵌水泥粒子后再进行磨光、抛光和超声波清洗等工艺可制得表面光洁度符合纳米硬度仪要求的
Abstract : The elastic modulus of cement particles was tested by nano2indentation techniques , and t he preparation technolo2
gy of sample was developed. The reasonable sample amount was analyzed. It was revealed t hat after a series of procedures such as mix wit h bakelite , solidification , polish gradually , ultrasonic washing , t he sample can reach t he requirement of nano2hard2 ness apparatus to smoot h degree of surface. The nano2indentation test result demonstrated t hat t he sample’s average was 17. 4 GPa and t he confidence interval of t he elastic modulus would be [ 14. 7 GPa ,20. 1 GPa ] wit h 95 %probability , which supposed t he elastic modulus of t he cement particle satisfied t he normal distribution. 20 samples were necessary and enough for t his exper2 iment in order to gain reliable dates.