基于窗口精细结构描述的岩体质量评价及力学参数综合优选

This thesis bases on horizontal research project “The research about the fine structure and mechanical parameters of abutment jointed rock mass of high arch dam on Jinping Ⅰ Hydropower Station, Yalong River” and “The research about the fine structure and mechanical parameters of the columnar basalt rock mass on Baihetan Hydropower Station, Jinsha River”. A rounded system about the fine structure description and rock mass classification is established. This research mainly contains six aspects as follow: (1) Methods about fine structure description of the window rock mass; (2) The window rock mass classification about the fine structure; (3) Model test study of intermittent joints; (4) Window rock mass strength theory; (5) Numerical experimentations about window rock mass; (6) The multi-source fusion of mechanical parameters based on Bayes principle. Variation of intact rock strength and joint conditions with the weathering and relaxation degree is studied through the description of window rock mass. And four principal parameters: intact rock point load strength, integration degree of window rock mass, joint conditions, and groundwater condition is selected to assess the window rock mass. Window rock mass is classified into three types using the results of window rock mass fine structure description combined with joints develop model. Scores about intact rock strength, integrality condition, divisional plane condition and groundwater conditions are given based on window rock mass fine structure description. Then quality evaluation about two different types of rock mass: general joint structure and columnar jointing structure are carried out to use this window rock mass classification system. Application results show that the window rock mass classification system is effective and applicable. Aimed at structural features of window structure of “the rock mass damaged by recessive fracture”, model tests and numerical models are designed about intermittent joints. By conducting model tests we get shear strength under different normal stress in integrated samples, through samples and intermittent joints samples. Also, the changing trends of shear strength in various connectivity rates are analyzed. We numerically simulate the entire process of direct shear tests by using PFC2D. In order to tally the stress-strain curve of numerical simulation with experimental tests about both integrated samples and through samples, we adjust mechanical factors between particles. Through adopting the same particle geometric parameter, the numerical sample of intermittent joints in different connective condition is re-built. At the same time, we endow the rock bridges and joints in testing samples with the fixed particle contacting parameters, and conduct a series of direct shear tests. Then the destructive process and mechanical parameters in both micro-prospective and macro-prospective are obtained. By synthesizing the results of numerical and sample tests and analyzing the evolutionary changes of stress and strain on intermittent joints plane, we conclude that the centralization of compressive stress on rock bridges increase the shear strength of it. We discuss the destructive mechanics of intermittent joints rock under direct shear condition, meanwhile, divide the whole shear process into five phases, which are elasticity phase, fracture initiation phase, peak value phase, after-peak phase and residual phase. In development of strength theory, the shear strength mechanisms of joint and rock bridge are analyzed respectively. In order to apply the deducted formulation conveniently in the real projects, a relationship between these formulations and Mohr-Coulomb hypothesis is built up. Some sets of numerical simulation methods, i.e. the distinct element method (UDEC) based on in-situ geology mapping are developed and introduced. The working methods about determining mechanical parameters of intact rock and joints in numerical model are studied. The operation process and analysis results are demonstrated detailed from the research on parameters of rock mass based on numerical test in the Jinping Ⅰ Hydropower Station and Baihetan Hydropower Station. By comparison，the advantages and disadvantages are discussed. Results about numerical simulation study show that we can get the shear strength mechanical parameters by changing the load conditions. The multi-source rock mass mechanical parameters can be fused by the Bayes theory, which are test value, empirical value and theoretical value. Then the value range and its confidence probability of different rock mass grade are induced and these data supports the reliability design.

喜马拉雅西构造结区岩石圈演化三维有限元动力学数值模拟

The past three decades have seen numerous attempts to numerically model stress and strain patterns in the lithosphere of the Earth on both global and regional scales. This efforts have been indispensable in identifying the features we need to include in our endeavour to develop better models of our planet’s lithosphere and they have also raised our awareness for the many unresolved issue in the deep geodynamical issues that need to be addressed in the future. Nonetheless, in most models, the lithosphere is treated as a single layer with depth-averaged properties, and as the same distribution in the stress and strain fields, and as deforming under plane strain. All these above make a great hander for its reality and degree of recognition. As the beginning in this paper, some principal numerical models and results on the evolution of Tibetan plateau are reviewed and analyzed. Then, the geological and geophysical expedition on the Western Himalayan Syntaxis is briefly reviewed. Furthermore, we analysis the feature in deep geophysical field studies in this area and adjacent regions. Because, for most continents, stress models driven by plate boundary forces have successfully reproduced the main characteristics of the stress and strain field, we present a set of three-dimensional models of lithosphere system for a simplified geometry of the Western Himalayan Syntaxis area and its adjacent regions, where we try to match the first-order characteristics of the stress and strain fields of lithosphere since 10 Ma, and deformation and geodynamical evolution process in former 2Ma. Of course, the kinematic boundary conditions of the stress models driven by plate boundary forces were applied. The rheology plays a significant role in the lithospheric tectonics, which lead to different rheological parameters were used in different works although the have the same constitutive equations in models. So, in this paper we do not aim to produce all characteristics of the Western Himalayan Syntaxis areas’ stress and strain fields by the choices of various parameters, but rather the dynamic response between various rheological parameters and stress and strain fields. We have chosen to concentrate on the importance of rheology and lateral strength variations for lithospheric stress and strain patterns and use our findings to build a model of the Western Himalayan Syntaxis areas. In doing so, we want to go beyond purely elastic models or purely viscoelastic models. Compared the results of the crust viscosity in the Western Himalayan Syntaxis areas, we believed that, when various viscoelastic models are adopted, the selection of the coefficient of viscosity in the Western Syntaxis area has important influence on the its uplifts and evolutions. A wider uplift ranges and gently elevation was observed at the same time when a lower viscosity was used in our models, and vice versa. Data of stress magnitudes are not available, but it is clear that the stress levels must be at or below the failure threshold of rock under compression. Under these criteria, the calculation results show that the viscosity in the Western Syntaxis area should be smaller than 1023Pa.s When elastic model is adopted in relatively rigid Tarim basin, obvious changes are induced to the stress and strain fields of the whole Western Syntaxis area. We found that rigid block of lithosphere reduced stress levels within its interior and that, at the edges of such regions, stress orientation can change. Furthermore there is no evidence that such rigid regions act as stress barriers in that they shield areas in opposite sides of the structure from the influence of one another. In our models, the upper crustal material of the Western Syntaxis area does not turns to move westward. Whereas, because of the stress and strain fields have been decoupling at the interior of the lithosphere, we can get the results that the deep material must not move westward.

The operation of multiquantum-well barriers

The multiquantum barrier (MQB), proposed by Iga et al in 1986, has been shown by several researchers to be an effective structure for improving the operating characteristics of laser diodes. These improvements include a reduction in the laser threshold current and increased characteristic temperatures. The operation of the MQB has been described as providing an increased barrier to electron overflow by reflecting high energy electrons trying to escape from the active region of the laser.This is achieved in a manner analogous to a Bragg reflector in optics. This thesis presents an investigation of the effectiveness of the MQB as an electron reflector. Numerical models have been developed for calculating the electron reflection due to MQB. Novel optical and electrical characterisation techniques have been used to try to measure an increase in barrier height due to the MQB in AlGaInP.It has been shown that the inclusion of MQB structures in bulk double heterostructure visible laser diodes can halve the threshold current above room temperature and the characteristic temperature of these lasers can be increased by up to 20K.These improvements are shown to occur in visible laser diodes even with the inclusion of theoretically ineffective MQB structures, hence the observed improvement in the characteristics of the laser diodes described above cannot be uniquely attributed to an increased barrier height due to enhance electron reflection. It is proposed here that the MQB improves the performance of laser diodes by proventing the diffusion of zinc into the active region of the laser. It is also proposed that the trapped zinc in the MQB region of the laser diode locally increases the p-type doping bringing the quasi-Fermi level for holes closer to the valence band edge thus increasing the barrier to electron overflow in the conduction band.

High speed optical communication systems: From modulation formats to radically new fibres

High volumes of data traffic along with bandwidth hungry applications, such as cloud computing and video on demand, is driving the core optical communication links closer and closer to their maximum capacity. The research community has clearly identifying the coming approach of the nonlinear Shannon limit for standard single mode fibre [1,2]. It is in this context that the work on modulation formats, contained in Chapter 3 of this thesis, was undertaken. The work investigates the proposed energy-efficient four-dimensional modulation formats. The work begins by studying a new visualisation technique for four dimensional modulation formats, akin to constellation diagrams. The work then carries out one of the first implementations of one such modulation format, polarisation-switched quadrature phase-shift keying (PS-QPSK). This thesis also studies two potential next-generation fibres, few-mode and hollow-core photonic band-gap fibre. Chapter 4 studies ways to experimentally quantify the nonlinearities in few-mode fibre and assess the potential benefits and limitations of such fibres. It carries out detailed experiments to measure the effects of stimulated Brillouin scattering, self-phase modulation and four-wave mixing and compares the results to numerical models, along with capacity limit calculations. Chapter 5 investigates hollow-core photonic band-gap fibre, where such fibres are predicted to have a low-loss minima at a wavelength of 2μm. To benefit from this potential low loss window requires the development of telecoms grade subsystems and components. The chapter will outline some of the development and characterisation of these components. The world's first wavelength division multiplexed (WDM) subsystem directly implemented at 2μm is presented along with WDM transmission over hollow-core photonic band-gap fibre at 2μm. References: [1]P. P. Mitra, J. B. Stark, Nature, 411, 1027-1030, 2001 [2] A. D. Ellis et al., JLT, 28, 423-433, 2010.

CBE-conveyor: a case-based reasoning system to assist engineers in designing conveyor systems

In this paper, we address the use of CBR in collaboration with numerical engineering models. This collaborative combination has a particular application in engineering domains where numerical models are used. We term this domain “Case Based Engineering” (CBE), and present the general architecture of a CBE system. We define and discuss the general characteristics of CBE and the special problems which arise. These are: the handling of engineering constraints of both continuous and nominal kind; interpolation over both continuous and nominal variables, and conformability for interpolation. In order to illustrate the utility of the method proposed, and to provide practical examples of the general theory, the paper describes a practical application of the CBE architecture, known as CBE-CONVEYOR, which has been implemented by the authors.Pneumatic conveying is an important transportation technology in the solid bulks conveying industry. One of the major industry concerns is the attrition of powders and granules during pneumatic conveying. To minimize the fraction of particles during pneumatic conveying, engineers want to know what design parameters they should use in building a conveyor system. To do this, engineers often run simulations in a repetitive manner to find appropriate input parameters. CBE-Conveyor is shown to speed up conventional methods for searching for solutions, and to solve problems directly that would otherwise require considerable intervention from the engineer.

Levitated liquid droplets in AC and DC magnetic field

Electromagnetic levitation of liquid metal droplets can be used to measure the properties of highly reactive liquid materials. Two independent numerical models, the commercial COMSOL and the spectral-collocation based free surface code SPHINX, have been applied to solve the transient electromagnetic, fluid flow and thermodynamic equations, which describe the levitated liquid motion and heating processes. The SPHINX model incorporates free surface deformation to accurately model the oscillations that result from the interaction between the electromagnetic and gravity forces, temperature dependent surface tension, magnetically controlled turbulent momentum transport. The models are adapted to incorporate periodic laser heating at the top of the droplet, which is used to measure the thermal conductivity of the material. Novel effects in the levitated droplet of magnetically damped turbulence and nonlinear growth of velocities in high DC magnetic field are analysed.

Determining surface tension using DC magnetic levitation

The values of material physical properties are vital for the successful use of numerical simulations for electromagnetic processing of materials. The surface tension of materials can be determined from the experimental measurement of the surface oscillation frequency of liquid droplets. In order for this technique to be used, a positioning field is required that results in a modification to the oscillation frequency. A number of previous analytical models have been developed that mainly focus on electrically conducting droplets positioned using an A.C. electromagnetic field, but due to the turbulent flow resulting from the high electromagnetic fields required to balance gravity, reliable measurements have largely been limited to microgravity. In this work axisymmetric analytical and numerical models are developed, which allow the surface tension of a diamagnetic droplet positioned in a high DC magnetic field to be determined from the surface oscillations. In the case of D.C. levitation there is no internal electric currents with resulting Joule heating, Marangoni flow and other effects that introduce additional physics that complicates the measurement process. The analytical solution uses the linearised Navier-Stokes equations in the inviscid case. The body force from a DC field is potential, in contrast to the AC case, and it can be derived from Maxwell equations giving a solution for the magnetic field in the form of a series expansion of Legendre polynomials. The first few terms in this expansion represent a constant and gradient magnetic field valid close to the origin, which can be used to position the droplet. Initially the mathematical model is verified in microgravity conditions using a numerical model developed to solve the transient electromagnetics, fluid flow and thermodynamic equations. In the numerical model (as in experiment) the magnetic field is obtained using electrical current carrying coils, which provides the confinement force for a liquid droplet. The model incorporates free surface deformation to accurately model the oscillations that result from the interaction between the droplet and the non-uniform external magnetic field. A comparison is made between the analytical perturbation theory and the numerical pseudo spectral approximation solutions for small amplitude oscillations.

Transfer Across the Air-Sea Interface

The efficiency of transfer of gases and particles across the air-sea interface is controlled by several physical, biological and chemical processes in the atmosphere and water which are described here (including waves, large- and small-scale turbulence, bubbles, sea spray, rain and surface films). For a deeper understanding of relevant transport mechanisms, several models have been developed, ranging from conceptual models to numerical models. Most frequently the transfer is described by various functional dependencies of the wind speed, but more detailed descriptions need additional information. The study of gas transfer mechanisms uses a variety of experimental methods ranging from laboratory studies to carbon budgets, mass balance methods, micrometeorological techniques and thermographic techniques. Different methods resolve the transfer at different scales of time and space; this is important to take into account when comparing different results. Air-sea transfer is relevant in a wide range of applications, for example, local and regional fluxes, global models, remote sensing and computations of global inventories. The sensitivity of global models to the description of transfer velocity is limited; it is however likely that the formulations are more important when the resolution increases and other processes in models are improved. For global flux estimates using inventories or remote sensing products the accuracy of the transfer formulation as well as the accuracy of the wind field is crucial.

Challenges in integrative approaches to modelling the marine ecosystems of the North Atlantic: Physics to fish and coasts to ocean

It has long been recognised that there are strong interactions and feedbacks between climate, upper ocean biogeochemistry and marine food webs, and also that food web structure and phytoplankton community distribution are important determinants of variability in carbon production and export from the euphotic zone. Numerical models provide a vital tool to explore these interactions, given their capability to investigate multiple connected components of the system and the sensitivity to multiple drivers, including potential future conditions. A major driver for ecosystem model development is the demand for quantitative tools to support ecosystem-based management initiatives. The purpose of this paper is to review approaches to the modelling of marine ecosystems with a focus on the North Atlantic Ocean and its adjacent shelf seas, and to highlight the challenges they face and suggest ways forward. We consider the state of the art in simulating oceans and shelf sea physics, planktonic and higher trophic level ecosystems, and look towards building an integrative approach with these existing tools. We note how the different approaches have evolved historically and that many of the previous obstacles to harmonisation may no longer be present. We illustrate this with examples from the on-going and planned modelling effort in the Integrative Modelling Work Package of the EURO-BASIN programme.

Radiotherapy in the presence of contrast agents: a general figure of merit and its application to gold nanoparticles

Delivering sufficient dose to tumours while sparing surrounding tissue is one of the primary challenges of radiotherapy, and in common practice this is typically achieved by using highly penetrating MV photon beams and spatially shaping dose. However, there has been a recent increase in interest in the possibility of using contrast agents with high atomic number to enhance the dose deposited in tumours when used in conjunction with kV x-rays, which see a significant increase in absorption due to the heavy element's high-photoelectric cross-section at such energies. Unfortunately, the introduction of such contrast agents significantly complicates the comparison of different source types for treatment efficacy, as the dose deposited now depends very strongly on the exact composition of the spectrum, making traditional metrics such as beam quality less valuable. To address this, a 'figure of merit' is proposed, which yields a value which enables the direct comparison of different source types for tumours at different depths inside a patient. This figure of merit is evaluated for a 15 MV LINAC source and two 150 kVp sources (both of which make use of a tungsten target, one with conventional aluminium filtration, while the other uses a more aggressive thorium filter) through analytical methods as well as numerical models, considering tissue treated with a realistic concentration and uptake ratio of gold nanoparticle contrast agents (10 mg ml(-1) concentration in 'tumour' volume, 10: 1 uptake ratio). Finally, a test case of human neck phantom is considered with a similar contrast agent to compare the abstract figure to a more realistic treatment situation. Good agreement was found both between the different approaches to calculate the figure of merit, and between the figure of merit and the effectiveness in a more realistic patient scenario. Together, these observations suggest that there is the potential for contrast-enhanced kilovoltage radiation to be a useful therapeutic tool for a number of classes of tumour on dosimetric considerations alone, and they point to the need for further research in this area.

Observation and characterization of laser-driven Phase Space Electron Holes

The direct observation and full characterization of a phase space electron hole (EH) generated during laser-matter interaction is presented. This structure, propagating in a tenuous, nonmagnetized plasma, has been detected via proton radiography during the irradiation with a ns laser pulse (I?2 ˜ 1014 W/cm2) of a gold hohlraum. This technique has allowed the simultaneous detection of propagation velocity, potential, and electron density spatial profile across the EH with fine spatial and temporal resolution allowing a detailed comparison with theoretical and numerical models.

A Comparison of the Flow Structures and Losses Within Vaned and Vaneless Stators for Radial Turbines

This paper details the numerical analysis of different vaned and vaneless radial inflow turbine stators. Selected results are presented from a test program carried out to determine performance differences between the radial turbines with vaned stators and vaneless volutes under the same operating conditions. A commercial computational fluid dynamics code was used to develop numerical models of each of the turbine configurations, which were validated using the experimental results. From the numerical models, areas of loss generation in the different stators were identified and compared, and the stator losses were quantified. Predictions showed the vaneless turbine stators to incur lower losses than the corresponding vaned stator at matching operating conditions, in line with the trends in measured performance.. Flow conditions at rotor inlet were studied and validated with internal static pressure measurements so as to judge the levels of circumferential nonuniformity for each stator design. In each case, the vaneless volutes were found to deliver a higher level of uniformity in the rotor inlet pressure field. [DOI: 10.1115/1.2988493]

Transonic aeroelastic simulation for instability searches and uncertainty analysis

In this paper the use of eigenvalue stability analysis of very large dimension aeroelastic numerical models arising from the exploitation of computational fluid dynamics is reviewed. A formulation based on a block reduction of the system Jacobian proves powerful to allow various numerical algorithms to be exploited, including frequency domain solvers, reconstruction of a term describing the fluid–structure interaction from the sparse data which incurs the main computational cost, and sampling to place the expensive samples where they are most needed. The stability formulation also allows non-deterministic analysis to be carried out very efficiently through the use of an approximate Newton solver. Finally, the system eigenvectors are exploited to produce nonlinear and parameterised reduced order models for computing limit cycle responses. The performance of the methods is illustrated with results from a number of academic and large dimension aircraft test cases.

Development and Validation of a Compact Thermal Model for an Aircraft Compartment

The development of accurate structural/thermal numerical models of complex systems, such as aircraft fuselage barrels, is often limited and determined by the smallest scales that need to be modelled. The development of reduced order models of the smallest scales and consequently their integration with higher level models can be a way to minimise the bottle neck present, while still having efficient, robust and accurate numerical models. In this paper a methodology on how to develop compact thermal fluid models (CTFMs) for compartments where mixed convection regimes are present is demonstrated. Detailed numerical simulations (CFD) have been developed for an aircraft crown compartment and validated against experimental data obtained from a 1:1 scale compartment rig. The crown compartment is defined as the confined area between the upper fuselage and the passenger cabin in a single aisle commercial aircraft. CFD results were utilised to extract average quantities (temperature and heat fluxes) and characteristic parameters (heat transfer coefficients) to generate CTFMs. The CTFMs have then been compared with the results obtained from the detailed models showing average errors for temperature predictions lower than 5%. This error can be deemed acceptable when compared to the nominal experimental error associated with the thermocouple measurements.

The CTFMs methodology developed allows to generate accurate reduced order models where accuracy is restricted to the region of Boundary Conditions applied. This limitation arises from the sensitivity of the internal flow structures to the applied boundary condition set. CTFMs thus generated can be then integrated in complex numerical modelling of whole fuselage sections.

Further steps in the development of an exhaustive methodology would be the implementation of a logic ruled based approach to extract directly from the CFD simulations numbers and positions of the nodes for the CTFM.

Optimising the locations of thermally sensitive equipment in an aircraft crown compartment

A Design of Experiments (DoE) analysis was undertaken to generate a list of configurations for CFD numerical simulation of an aircraft crown compartment. Fitted regression models were built to predict the convective heat transfer coefficients of thermally sensitive dissipating elements located inside this compartment. These are namely the SEPDC and the Route G. Currently they are positioned close to the fuselage and it is of interest to optimise the heat transfer for reliability and performance purposes. Their locations and the external fuselage surface temperature were selected as input variables for the DoE. The models fit the CFD data with values ranging from 0.878 to 0.978, and predict that the optimum locations in terms of heat transfer are when the elements are positioned as close to the crown floor as possible ( and ?min. limits), where they come in direct contact with the air flow from the cabin ventilation system, and when they are positioned close to the centreline ( and ?CL). The methodology employed allows aircraft thermal designers to optimise equipment placement in confined areas of an aircraft during the design phase. The determined models should be incorporated into global aircraft numerical models to improve accuracy and reduce model size and computational time. © 2012 Elsevier Masson SAS. All rights reserved.