二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.

软土亚塑性模型与一维非线性固结理论研究

With the great development of Tianjing New Coastal District economy, people need more land to build and live. Land subsidence, which is caused by its special engineering geological conditions, has restricted the further development in the district. Soft soil consolidation is main factor of land subsidence ;thus , on the basis of consolidation theory, the paper make further study on soft soils one-dimension nonlinear consolidation which contains two parts:(1) the nonlinear consolidation of permeability coefficient and compressibility coefficient changing with time and depth, which means real one-dimension nonlinear consolidation;(2) the non-homogeneous consolidation of permeability coefficient and compressibility coefficient only changing with depth. Firstly, nonlinear characteristics of soft soils are elaborated. Hypoplastic theory is introduced to establish a modified soft soils nonlinear constitutive model; the nonlinear governing equation of compressibility coefficient is built, and the nonlinear characteristics of compressibility coefficient are analyzed. Secondly, Considering Load Fluctuation and soil thickness changing ,the consolidation characteristics of single layer is discussed in the paper; meanwhile, on the basis of the Davis and Raymond’s hypothesis and single layer nonlinear consolidation equation, the doubled-layer one-dimension nonlinear consolidation equation is also derived. The solution of the equation is obtained by analytical method, and the consolidation characteristics of doubled-layer soft soil nonlinear theory is also analyzed. Finally, based on assumption that permeability coefficient and compressibility coefficient is varying along depth, single layer soil one-dimension non-homogeneous consolidation differential equation is derived; and the approximate solution is obtained. Furthermore, the single layer non-homogeneous consolidation is extended to double layer non-homogeneous consolidation theory. By using parabolic differential scheme, the matrix equation is established; and the solution of the matrix equation is obtained by chase method. Consolidation characteristics of soil soft single (double) layer non-homogeneous consolidation theory and Terzaghi’s theory are also discussed.

济阳-昌潍坳陷火成岩油气藏体系形成机理和分布规律研究

Jiyang & Changwei depressions are two neighboring depressions in Bahai Bay Basin, the famous oil rich basin in East China. The exploration activities in the past 40 years has proved that, within the basins, there exists not only plentiful sandstone hydrocarbon reservoirs (conventional), but also abundant special reservoirs as igneous rock, mudstone and conglomerate ones which have been knowing as the unconventional in the past, and with the prospecting activity is getting more and more detailed, the unconventional reservoirs are also getting more and more important for further resources, among which, the igneous lithological reservoir be of significance as a new research and exploration area. The purpose of this paper is, with the historical researches and data as base, the System Theory, Practice Theory and Modern Comprehensive Petroleum Geology Theory as guide, the theoretical and practice break through as the goal, and the existing problems in the past as the break through direction, to explore and establish a valid reservoir formation and distribution models for igneous strata in the profile of the eastern faulted basins. After investigating the distribution of the igneous rocks and review the history of the igneous rocks reservoirs in basins, the author focused on the following issues and correspondingly the following progresses have been made: 1.Come to a new basin evolution and structure model named "Combined-Basin-bodies Model" for Jiyang even Eastern faulted basins based on the study on the origin and evolution of Jiyang & Changwei basins, depending on this model, every faulted basin in the Bo-hai Bay Basin is consisted of three Basin-Bodies including the Lower (Mesozoic), Middle (Early Tertiary) and the Upper (Late Tertiary) Bodies, each evolved in different geo-stress setting and with different basin trend, shape and igneous-sedimentary buildings system, and from this one to next one, the basin experienced a kind of process named "shape changing" and "Style changing". 2. Supposed a serious of new realizations as follows (1) There were "multi-level magma sources" including Upper mantel and the Lower, Middle and even the Upper Shell magma Chambers in the historical Magma Processes in the basins; (2) There were "multi-magma accessing or pass" from the first level (Mantel faults) to the second, third and fourth levels (that is the different levels of fault in the basin sediment strata) worked in the geo-historical and magma processes; (3) Three tectonic magma cycles and more periods have been recognized those are matched with the "Basin -body-Model" and (4)The geo-historical magma processes were non-homogeneous in time and space scale and so the magma rocks distributed in "zones" or "belts". 3. The study of magma process's effect on basin petroleum conditions have been made and the following new conclusions were reached: (1) the eruptive rocks were tend to be matched with the "caped source rock", and the magma process were favorable to the maturing of the source rocks. (2) The magma process were fruitful to the accumulation of the non-hydrocarbon reservoirs however a over magma process may damage the grade of resource rock; (3) Eruptive activity provided a fruitful environment for the formation of such new reservoir rocks as "co-eruptive turbidity sandstones" and "thermal water carbonate rocks" and the intrusive process can lead to the origin of "metamorphism rock reservoir"; (4) even if the intrusive process may cause the cap rock broken, the late Tertiary intrusive rocks may indeed provide the lateral seal and act as the cap rock locally even regionally. All above progresses are valuable for reconstructing the magma-sedimentary process history and enriching the theory system of modem petroleum geology. 4. A systematic classification system has been provided and the dominating factors for the origin and distribution of igneous rock reservoirs have been worked out based on the systematic case studies, which are as follows: (1) The classification is given based on multi-factors as the origin type, litho-phase, type of reservoir pore, reservoir ability etc., (2) Each type of reservoir was characterized in a detailed way; (3) There are 7 factors dominated the intrusive reservoir's characteristics including depth of intrusion, litho-facies of surrounding rocks, thickness of intrusive rock, intrusive facies, frequency and size of the working faults, shape and tectonic deformation of rock, erosion strength of the rock and the time of the intrusion ect., in the contrast, 4 factors are for eruptive rocks as volcanic facies, frequency and size of the working faults, strength of erosion and the thermal water processing. 5. Several new concept including "reservoir litho-facies", "composite-volcanic facies" and "reservoir system" ect. Were suggested, based on which the following models were established: (1) A seven reservoir belts model for a intrusive unit profile and further more, (2) a three layers cubic model consisted of three layer as "metamorphic roe layer", "marginal layer" and "the core"; (3) A five zones vertical reservoir sequence model consisted of five litho-facies named A, B, C, D and E for a original lava unit and furthermore three models respectively for a erosion, subsidence and faulted lava unit; (4) A composite volcanic face model for a lava cone or a composite cone that is consisted of three facies as "crater and nearby face", "middle slope" and "far slope", among which, the middle slope face is the most potential reservoir area and producible for oil & gas. 6. The concept of "igneous reservoir" was redefined as the igneous, and then a new concept of "igneous reservoir system" was supposed which means the reservoir system consisted of igneous and associated non-igneous reservoirs, with non-hydrocarbon reservoir included. 7. The origin and distribution of igneous reservoir system were probed and generalized for the exploration applications, and origin models of the main reservoir sub-systems have been established including those of igneous, related non-igneous and non-hydrocarbon. For intrusive rocks, two reservoir formation models have been suggested, one is called "Original or Primary Model", and the another one is "Secondary Model"; Similarly, the eruptive rock reservoirs were divided in three types including "Highly Produced", "Moderately Produced" and "Lowly Produced" and accordingly their formation models were given off; the related non-igneous reservoir system was considered combination of eight reservoirs, among which some ones like the Above Anticline Trap are highly produced; Also, the non-hydrocarbon. Trap system including five kinds of traps was discussed. 8. The concept models for four reservoir systems were suggested, which include the intrusive system consisted of 7 kinds of traps, the land eruptive system with 6 traps, the under water eruptive system including 6 kinds of traps and the non-hydrocarbon system combined by 5 kinds of traps. In this part, the techniques for exploration of igneous reservoir system were also generalized and probed, and based on which and the geological progresses of this paper, the potential resources and distributions of every reservoir system was evaluated and about 186 millions of reserves and eight most potential non-hydrocarbon areas were predicted and outlined. The author believe that the igneous reservoir system is a very important exploration area and its study is only in its early stage, the framework of this paper should be filled with more detailed studies, and only along way, the exploration of igneous reservoir system can go into it's really effective stage.

三维油藏数值模拟正反演算法研究

The dynamic prediction of complex reservoir development is one of the important research contents of dynamic analysis of oil and gas development. With the increase development of time, the permeabilities and porosities of reservoirs and the permeability of block reservoir at its boundaries are dynamically changing. How to track the dynamic change of permeability and porosity and make certain the permeability of block reservoir at its boundary is an important practical problem. To study developing dynamic prediction of complex reservoir, the key problem of research of dynamic prediction of complex reservoir development is realizing inversion of permeability and porosity. To realize the inversion, first of all, the fast forward and inverse method of 3-dimension reservoir simulation must be studied. Although the inversion has been widely applied to exploration and logging, it has not been applied to3-dimension reservoir simulation. Therefore, the study of fast forward and inverse method of 3-dimension reservoir simulation is a cutting-edge problem, takes on important realistic signification and application value. In this dissertation, 2-dimension and 3-dimension fluid equations in porous media are discretized by finite difference, obtaining finite difference equations to meet the inner boundary conditions by Peaceman's equations, giving successive over relaxation iteration of 3-dimension fluid equations in porous media and the dimensional analysis. Several equation-solving methods are compared in common use, analyzing its convergence and convergence rate. The alternating direction implicit procedure of 2-dimension has been turned into successive over relaxation iteration of alternating direction implicit procedure of 3-dimension fluid equations in porous media, which possesses the virtues of fast computing speed, needing small memory of computer, good adaptability for heterogeneous media and fast convergence rate. The geological model of channel-sandy reservoir has been generated with the help of stochastic simulation technique, whose cross sections of channel-sandy reservoir are parabolic shapes. This method makes the hard data commendably meet, very suit for geological modeling of containing complex boundary surface reservoir. To verify reliability of the method, theoretical solution and numerical solution are compared by simplifying model of 3-dimension fluid equations in porous media, whose results show that the only difference of the two pressure curves is that the numerical solution is lower than theoretical at the wellbore in the same space. It proves that using finite difference to solve fluid equations in porous media is reliable. As numerical examples of 3-dimension heterogeneous reservoir of the single-well and multi-well, the pressure distributions have been computed respectively, which show the pressure distributions there are clearly difference as difference of the permeabilities is greater than one order of magnitude, otherwise there are no clearly difference. As application, the pressure distribution of the channel-sandy reservoir have been computed, which indicates that the space distribution of pressure strongly relies on the direction of permeability, and is sensitive for space distributions of permeability. In this dissertation, the Peaceman's equations have been modified into solving vertical well problem and horizontal well problem simultaneously. In porous media, a 3D layer reservoir in which contain vertical wells and horizontal wells has been calculated with iteration. For channel-sandy reservoir in which there are also vertical wells and horizontal wells, a 3D transient heterogeneous fluid equation has been discretized. As an example, the space distribution of pressure has been calculated with iteration. The results of examples are accord with the fact, which shows the modification of Peaceman's equation is correct. The problem has been solved in the space where there are vertical and horizontal wells. In the dissertation, the nonuniform grid permeability integration equation upscaling method, the nonuniform grid 2D flow rate upscaling method and the nonuniform grid 3D flow rate upscaling method have been studied respectively. In those methods, they enhance computing speed greatly, but the computing speed of 3D flow rate upscaling method is faster than that of 2D flow rate upscaling method, and the precision of 3D flow rate upscaling method is better than that of 2D flow rate upscaling method. The results also show that the solutions of upscaling method are very approximating to that of fine grid blocks. In this paper, 4 methods of fast adaptive nonuniform grid upscaling method of 3D fluid equations in porous media have been put forward, and applied to calculate 3D heterogeneous reservoir and channel-sandy reservoir, whose computing results show that the solutions of nonuniform adaptive upscaling method of 3D heterogeneous fluid equations in porous media are very approximating to that of fine grid blocks in the regions the permeability or porosity being abnormity and very approximating to that of coarsen grid blocks in the other region, however, the computing speed of adaptive upscaling method is 100 times faster than that of fine grid block method. The formula of sensitivity coefficients are derived from initial boundary value problems of fluid equations in porous media by Green's reciprocity principle. The sensitivity coefficients of wellbore pressure to permeability parameters are given by Peaceman's equation and calculated by means of numerical calculation method of 3D transient anisotropic fluid equation in porous media and verified by direct method. The computing results are in excellent agreement with those obtained by the direct method, which shows feasibility of the method. In the dissertation, the calculating examples are also given for 3D reservoir, channel-sandy reservoir and 3D multi-well reservoir, whose numerical results indicate: around the well hole, the value of the sensitivity coefficients of permeability is very large, the value of the sensitivity coefficients of porosity is very large too, but the sensitivity coefficients of porosity is much less than the sensitivity coefficients of permeability, so that the effect of the sensitivity coefficients of permeability for inversion of reservoir parameters is much greater than that of the sensitivity coefficients of porosity. Because computing the sensitivity coefficients needs to call twice the program of reservoir simulation in one iteration, realizing inversion of reservoir parameters must be sustained by the fast forward method. Using the sensitivity coefficients of permeability and porosity, conditioned on observed valley erosion thickness in wells (hard data), the inversion of the permeabilities and porosities in the homogeneous reservoir, homogeneous reservoir only along the certain direction and block reservoir are implemented by Gauss-Newton method or conjugate gradient method respectively. The results of our examples are very approximating to the real data of permeability and porosity, but the convergence rate of conjugate gradient method is much faster than that of Gauss-Newton method.

Au nanorods-semiconductor nanowire hybrid nanostructures : nanofabrication techniques and optoelectronic properties

The objective of this thesis is the exploration and characterization of novel Au nanorod-semiconductor nanowire hybrid nanostructures. I provide a comprehensive bottom-up approach in which, starting from the synthesis and theoretical investigation of the optical properties of Au nanorods, I design, nanofabricate and characterize Au nanorods-semiconductor nanowire hybrid nanodevices with novel optoelectronic capabilities compared to the non-hybrid counterpart. In this regards, I first discuss the seed-mediated protocols to synthesize Au nanorods with different sizes and the influence of nanorod geometries and non-homogeneous surrounding medium on the optical properties investigated by theoretical simulation. Novel methodologies for assembling Au nanorods on (i) a Si/SiO2 substrate with highly-ordered architecture and (ii) on semiconductor nanowires with spatial precision are developed and optimized. By exploiting these approaches, I demonstrate that Raman active modes of an individual ZnO nanowire can be detected in non-resonant conditions by exploring the longitudinal plasmonic resonance mediation of chemical-synthesized Au nanorods deposited on the nanowire surface otherwise not observable on bare ZnO nanowire. Finally, nanofabrication and detailed electrical characterization of ZnO nanowire field-effect transistor (FET) and optoelectronic properties of Au nanorods - ZnO nanowire FET tunable near-infrared photodetector are investigated. In particular we demonstrated orders of magnitude enhancement in the photocurrent intensity in the explored range of wavelengths and 40 times faster time response compared to the bare ZnO FET detector. The improved performance, attributed to the plasmonicmediated hot-electron generation and injection mechanism underlying the photoresponse is investigated both experimentally and theoretically. The miniaturized, tunable and integrated capabilities offered by metal nanorodssemicondictor nanowire device architectures presented in this thesis work could have an important impact in many application fields such as opto-electronic sensors, photodetectors and photovoltaic devices and open new avenues for designing of novel nanoscale optoelectronic devices.

Waveguide QED: Many-body bound-state effects in coherent and Fock-state scattering from a two-level system

Strong coupling between a two-level system (TLS) and bosonic modes produces dramatic quantum optics effects. We consider a one-dimensional continuum of bosons coupled to a single localized TLS, a system which may be realized in a variety of plasmonic, photonic, or electronic contexts. We present the exact many-body scattering eigenstate obtained by imposing open boundary conditions. Multiphoton bound states appear in the scattering of two or more photons due to the coupling between the photons and the TLS. Such bound states are shown to have a large effect on scattering of both Fock- and coherent-state wave packets, especially in the intermediate coupling-strength regime. We compare the statistics of the transmitted light with a coherent state having the same mean photon number: as the interaction strength increases, the one-photon probability is suppressed rapidly, and the two- and three-photon probabilities are greatly enhanced due to the many-body bound states. This results in non-Poissonian light. © 2010 The American Physical Society.

A computational model of coupled heat and moisture transfer with phase change in granular sugar during varying environmental conditions

As part of a comprehensive effort to predict the development of caking in granular materials, a mathematical model is introduced to model simultaneous heat and moisture transfer with phase change in porous media when undergoing temperature oscillations/cycling. The resulting model partial differential equations were solved using finite-volume procedures in the context of the PHYSICA framework and then applied to the analysis of sugar in storage. The influence of temperature on absorption/desorption and diffusion coefficients is coupled into the transport equations. The temperature profile, the depth of penetration of the temperature oscillation into the bulk solid, and the solids moisture content distribution were first calculated, and these proved to be in good agreement with experimental data. Then, the influence of temperature oscillation on absolute humidity, moisture concentration, and moisture migration for different parameters and boundary conditions was examined. As expected, the results show that moisture near boundary regions responds faster than farther away from them with surface temperature changes. The moisture absorption and desorption in materials occurs mainly near boundary regions (where interactions with the environment are more pronounced). Small amounts of solids moisture content, driven by both temperature and vapour concentration gradients, migrate between boundary and center with oscillating temperature.

Periodic transonic flow and control

The current understanding of periodic transonic flow is reviewed briefly. The effects of boundary-layer transition, non-adiabatic wall conditions and modifications to the aerofoil surface geometry at the shock interactions on periodic transonic flow are discussed. Through the methods presented, it is proposed that the frequency of periodic motion can be predicted with reasonable accuracy, but there are limitations on the prediction of buffet boundaries associated with periodic transonic flows. Several methods have been proposed by which the periodic motion may be virtually eliminated, most relevantly by altering the position of transition fix, contouring the aerofoils surface or adding a porous surface and a cavity in the region of shock interaction. In addition, it has been shown that heat transfer can have a significant effect on buffet.

Analysis of compressive membrane action in concrete slabs

This paper summarises the results obtained from non-linear finite-element analysis (NLFEA) of a series of reinforced-concrete one-way slabs with various boundary conditions representative of a bridge deck slab strip in which compressive membrane action governs the structural behaviour. The application of NLFEA for the optimum analysis and design of in-plane restrained concrete slabs is explored. An accurate material model and various equation solution methods were assessed to find a suitable finite-element method for the analysis of concrete slabs in which arching action occurs. Finally, the results from the NLFEA are compared and validated with those from various experimental test data. Significantly, the numerical analysis was able to model the arching action that occurred as a result of external in-plane restraint at the supports and which enhanced the ultimate strength of the slab. The NLFEA gave excellent predictions for the ultimate load-carrying capacity and far more accurate predictions than those obtained using standard flexural or elastic theory.

Modelling the mechanical response of the Reed-mouthpiece-lip system of a clarinet. Part I. A one-dimensional distributed model

The motion of a clarinet reed that is clamped to a mouthpiece and supported by a lip is simulated in the time-domain using finite difference methods. The reed is modelled as a bar with non-uniform cross section, and is described using a one-dimensional, fourth-order partial differential equation. The interactions with the mouthpiece Jay and the player's lip are taken into account by incorporating conditional contact forces in the bar equation. The model is completed by clamped-free boundary conditions for the reed. An implicit finite difference method is used for discretising the system, and values for the physical parameters are chosen both from laboratory measurements and by accurate tuning of the numerical simulations. The accuracy of the numerical system is assessed through analysis of frequency warping effects and of resonance estimation. Finally, the mechanical properties of the system are studied by analysing its response to external driving forces. In particular, the effects of reed curling are investigated.

Asymmetric Carbon-Carbon Bond Forming Reactions Catalysed by Metal(II) Bis(oxazoline) Complexes Immobilized using Supported Ionic Liquids

A series of bis(oxazoline) metal(II) complexes has been supported on silica and carbon supports by non-covalent immobilisation using an ionic liquid. The catalytic performance of these solids was compared for the enantioselective Diels-Alder reaction between N-acryloyloxazolidinone and cyclopentadiene and the Mukaiyama-aldol reaction between methyl pyruvate and 1-methoxy-1-trimethylsilyloxy-propene. In both reactions the enantioselectivity was strongly influenced by the choice of support displaying enantioselectivies (ee values) up to 40% higher than those conducted under homogeneous reaction conditions.

Epitaxial strain and electric boundary condition effects on the structural and ferroelectric properties of BiFeO3 films

The influence of both compressive and tensile epitaxial strain along with the electrical boundary conditions on the ferroelastic and ferroelectric domain patterns of bismuth ferrite films was studied. BiFeO3 films were grown on SrTiO3(001), DyScO3(110), GdScO3(110), and SmScO3(110) substrates to investigate the effect of room temperature in-plane strain ranging from -1.4% to +0.75%. Piezoresponse force microscopy, transmission electron microscopy, x-ray diffraction measurements, and ferroelectric polarization measurements were performed to study the properties of the films. We show that BiFeO3 films with and without SrRuO3 bottom electrode have different growth mechanisms and that in both cases reduction of the domain variants is possible. Without SrRuO3, stripe domains with reduced variants are formed on all rare earth scandate substrates because of their monoclinic symmetry. In addition, tensile strained films exhibit a rotation of the unit cell with increasing film thickness. On the other side, the presence of SrRuO3 promotes step flow growth of BiFeO3. In case of vicinal SrTiO3 and DyScO3 substrates with high quality SrRuO3 bottom electrode and a low miscut angle of approximate to 0.15 degrees we observed suppression of the formation of certain domain variants. The quite large in-plane misfit of SrRuO3 with GdScO3 and SmScO3 prevents the growth of high quality SrRuO3 films and subsequent domain variants reduction in BiFeO3 on these substrates, when SrRuO3 is used as a bottom electrode.

Finite Element Analysis of the Enhanced Strength of Laterally Restrained RC Slabs

This paper presents the numerical simulation of the ultimate behaviour of 85 one-way and two-way spanning laterally restrained concrete slabs of variable thickness, span, reinforcement ratio, strength and boundary conditions reported in literature by different authors. The developed numerical model was described and all the assumptions were illustrated. ABAQUS, a Finite Element Analysis suite of software, was employed. Non-linear implicit static general analysis method offered by ABAQUS was used. Other analysis methods were also discussed in general in terms of application such as Explicit Dynamic Analysis and Riks method. The aim is to demonstrate the ability and efficacy of FEA to simulate the ultimate load behaviour of slabs considering different material properties and boundary conditions. The authors intended to present a numerical model that provides consistent predictions of the ultimate behaviour of laterally restrained slabs that could be used as an alternative for expensive real life testing as well as for the design and assessment of new and existing structures respectively. The enhanced strength of laterally-restrained slabs compared with conventional design methods predictions is believed to be due to compressive membrane action (CMA). CMA is an inherent phenomenon of laterally restrained concrete beams/slabs. The numerical predictions obtained from the developed model were in good correlation with the experimental results and with those obtained from the CMA method developed at the Queen’s University Belfast, UK.

The Reliability of the Arc-Length Method in the Analysis of Mode-Jumping Problems

The Arc-Length Method is a solution procedure that enables a generic non-linear problem to pass limit points. Some examples are provided of mode-jumping problems solutions using a commercial nite element package, and other investigations are carried out on a simple structure of which the numerical solution can be compared with an analytical one. It is shown that Arc-Length Method is not reliable when bifurcations are present in the primary equilibrium path; also the presence of very sharp snap-backs or special boundary conditions may cause convergence diÆculty at limit points. An improvement to the predictor used in the incremental procedure is suggested, together with a reliable criteria for selecting either solution of the quadratic arc-length constraint. The gap that is sometimes observed between the experimantal load level of mode-jumping and its arc-length prediction is explained through an example.

Finite-element investigation of cold-formed steel portal frames in fire

This paper presents the results of a full-scale site fire test performed on a cold-formed steel portal frame building with semi-rigid joints. The purpose of the study is to establish a performance-based approach for the design of such structures in fire boundary conditions. In the full-scale site fire test, the building collapsed asymmetrically at a temperature of 714°C. A non-linear elasto-plastic finite-element shell model is described and is validated against the results of the full-scale test. A parametric study is presented that highlights the importance of in-plane restraint from the side rails in preventing an outwards sway failure for both a single portal and full building geometry model. The study also demonstrates that the semi-rigidity of the joints should be taken into account in the design. The single portal and full building geometry models display a close match to site test results with failure at 682°C and 704°C, respectively. A design case is described in accordance with Steel Construction Institute design recommendations. The validated single portal model is tested with pinned bases, columns protected, realistic loading and rafters subject to symmetric uniform heating in accordance with the ISO 834 standard fire curve; failure occurs at 703°C.