Modelling of submarine landslides and generated water waves

The slide of unstable sedimentary bodies and their hydraulic effects are studied by numerical means. A two-dimensional fluid mechanics model based on Navier-Stokes equations has been developed considering the sediments and water as a mixture. Viscoplastic and diffusion laws for the sediments have been introduced into the model. The numerical model is validated with an analytical solution for a Bingham flow. Laboratory experiments consisting in the slide of gravel mass have been carried out. The results of these experiments have shown the importance of the sediment rheology and the diffusion. The model parameters are adjusted by trial and error to match the observed “sandflow”.

岩石中地震波速度和衰减的应力变化响应研究

Stress change is one of key factors in seismic nucleating and triggering; therefore for understanding and forecasting earthquakes, it is necessary to research on stress status and its changes in rocks. Propagating in underground structures, wave velocity and attenuation contain information on stress changes of the Earth’s interior. For a better understanding of relationship between seismic data and stress changes, modeling and ultrasonic test supply significant references. In this article, acoustoelastic theory is introduced to explain nonlinear elastic characteristics of rocks. Based on the acoustoelastic theory, a solid-fluid coupled model is given to calculate velocity under different stress for porous and liquid fulfilled rocks. Except for the stress-velocity relationship, effects of pore pressure induced stress changes on ultrasonic coda attenuation are also studied. Intrinsic attenuation quality factors are calculated for a comparison purpose. Finally, the relationship between elastic constants and stress changes is thoroughly investigated, a mixture model from two phases of Hooke media is introduced to explain the differences between dynamic and static moduli, a relation among wave length, wave velocities and elastic moduli considering dimension of microstructure, dimension and state of surface between phases is presented. The most important aspect of this work is exploring and establishing relationships between the seismic properties of rocks and changes of their stress conditions, which will have its application in earthquake forecast and seismic hazard.

班公湖带多不杂超大型富金斑岩铜矿床的成岩成矿年代学、岩石学及高氧化岩浆¬流体-成矿作用

Duobuza copper deposit, newly discovered typical gold-rich porphyry copper deposit with superlarge potential, is located in the Tiegelong Mesozoic tectonic -magmatic arc of the southern edge of Qiangtang block and the northern margin of Bangonghu-Nujiang suture. Quartz diorite porphyrite and grandiorite porphyry, occurred in stock, are the main ore-bearing porphyries. As the emplacement of porphyry stock, a wide range of hydrothermal alteration has developed. Within the framework of the ore district, abundant hydrothermal magnetite developed, and the relationship between precipitation of copper and gold and hydrothermal magnetite seems much close. Correspondingly, a series of veinlets and network veinlets occurred in all alteration zones. Therefore, systematic research on such a superlarge high-grade Duobuza gold-rich porphyry copper deposit can fully revealed the metallogenic characteristics of gold-rich porphyry copper deposits in this region, establish metallogenetic model and prospecting criteria, and has important practical significance on the promotion of regional exploration. In addition, this research on it can enrich metallogenic theory of strong oxidation magma-fluid to gold-rich porphyry copper deposit, and will be helpful to understand the metallogenic characteristics in early of subduction of Gangdese arc stages and its entire evolution history of the Qinghai-Tibet Plateau, the temporal and spatial distribution of ore deposits and their geodynamics settings. Northern ore body of Duobuza copper deposit have been controlled with width (north-south) about 100 ~ 400 m, length (east-west) about 1400 m, dip of 200 °, angle of dip 65 °~ 80 °. And controlled resource amount is of 2.7 million tons Cu with grade 0.94% and 13 tons Au with 0.21g/tAu. Overall features of ore body are large scale, higher grade copper, gold-rich. Ore occurred in the body of granodiotite porphyry and quartz diorite porphyrite and its contact zone with wall rock. Through the detailed mapping and field work studies, some typies of alteration are identificated as follows: albitization, biotititation, sericitization, silication, epidotization, chloritization, carbonatization, illitization, kaolinization and so on. The range of alteration is more than 10km2. Wall alteration zone can be divided into potassic alteration, moderate argillization alteration, argillization, illite-hydromuscovite or propylitization from ore-bearing porphyry center outwards, but phyllic alteration has not well developed and only sericite-quartz veins occurred in local area. Moreover, micro-fracture is development in ore district , and correspondingly a series of veinlets are development as follows: biotite vein (EB type), K-feldspar-biotite-chalcopyrite-quartz vein, magnetite-antinolite-K-feldspar vein, quartz-chalcopyrite-magnetite veins (A-type), quartz-magnetite-biotite-K-feldspar vein, chalcopyrite veinlets in potassic alteration zone; (2) chalcopyrite occurring in the center vein–quartz vein (B type), chalcopyrite veinlets, chalcopyrite-gypsum vein in intermediate argillization alteration; (3) chalcopyrite- pyrite-quartz vein, pyrite-quartz vein, chalcopyrite-gypsum veins, quartz-gypsum- molybdenite-chalcopyrite vein in argillization alteration; (4) gypsum veins, quartz-(molybdenite)-chalcopyrite vein, quartz-pyrite vein, gypsum- chalcopyrite vein, potassium feldspar veinlets, Carbonate veins, quartz-magnetite veins in the wall rock. In short, various veins are very abundant within the framework of the ore district. The results of electronic probe microscopy analysis (EMPA) indicate that Albite (Ab 91.5~99.7%) occurred along the rim of plagioclase phenocryst and fracture, and respresents the earliest stages of alteration. K-feldspar (Or 75.1~96.9%) altered plagioclase phenocryst and matrix or formed secondary potassium feldspar veinlets. Secondary biotite occurred mainly in phenocryst, matrix and veinlets, belong to magnesium-rich biotite formed under the conditions of high-oxidation magma- hydrothermal. Chloritization developed in all alteration zones and alterd iron- magnesium minerals such as biotite and hornblende and then formed chlorite veinlets. As the temperature rises, Si in the tetrahedral site of chlorite decreased, and chlorite component evolved from diabantite to ripiolite. The consistent 280℃~360℃ of formation temperature hinted that chlorite formed on the same temperature range in all alteration zones. However, formation temperature range of chlorite from the gypsum-carbonate-chlorite vein was 190℃~220℃, and it may be the product of the latest stage of hydrothermal activity. The closely relationship between biotite and rutile indicate that most of rutiles are precipitated in the process of biotite alteration and recrystallization. In addition, the V2O3 concentration of rutile from ore body in Duobuza gold-rich porphyry copper deposit is >0.4％, indicate that V concentration in rutile has important significance on marking main ore body of porphyry copper deposit. Apatites from Duobuza deposit all are F-rich. And apatite in the wall rock contained low MnO content and relatively high FeO content, which may due to the basaltic composition of the wall rocks. The MnO in apatite from altered porphyry show a strong positive correlation with FeO. In addition, Cl/F ratio of apatite from wall rock was highest, followed by the potassic alteration zone and potassic alteration zone overprinted by moderate argillization alteration was the lowest. SO2 in Apatite are in the scope of 0 to 0.66%, biotite in the apatite has the highest SO2, followed by the potassic alteration zone, potassic alteration zone overprinted by moderate argillization alteration, and the lowest in the surrounding rocks, which may be caused by the decrease of oxygen fugacity of hydrothermal fluid and S exhaust by sulfide precipitation in potassic alteration. Magnetite in the wall rock have higher Cr2O3 and lower Al2O3 features compared with altered porphyry, this may be due to basalt wall rock generally has high Cr content. And magnetites have higher TiO2 content in potassic alteration than moderate argillization alteration overprinted by potassic alteration, argillization and wall rock, suggested that its formation temperature in potassic alteration was the highest among them. The ore minerals mainly are chalcopyrite and bornite, and Au contents of chalcopyrite, bornite, and pyrite are similar with chalcopyrite slightly higher. The Eu* negative anomaly of disseminated chalcopyrite was relatively lower than chalcopyrite in veinlets. Within a drill hole, the Eu* negative anomaly of disseminated chalcopyrite was gradually larger from bottom to top. Magnetite has the same distribution model, with obvious negative Eu* abnormal, and ΣREE in great changes. The gypsum has the highest ΣREE content and the obvious negative anomaly, and biotite obviously has the Eu* abnormal. Based on the petrographic and geochemical characteristics, five series of magmatic rocks can be broadly classified; they are volcanic rocks of the normal island arc, high-Nb basaltic rocks, adakites, altered porphyry and diorite. The Sr, Nd, Hf isotopes and geochemistry of various series of magmatic rock show that they may be the result of mixing between basic magma and various degrees of acid magma coming from lower crust melted by high temperature basic underplating from partial melting of the subduction sediment melt metasomatic mantle wedge. Furthermore S isotope and Pb isotope of the sulfide, ore-bearing porphyries and volcanic rocks indicated ore-forming source is the mantle wedge metasomatied by subduction sediment melt. Oxygen fugacity of magma estimated by Fe2O3/FeO of whole rock and zircon Ce4+/Ce3+ indicated that the oxidation of basalt-andesitic rocks is higher than ore-forming porphyry, and might imply high-oxidation characteristics of underplated basic magma. Its high oxidative mechanism is likely mantle sources metasomatied by subduction sediment magma, including water and Fe3+. And such high oxidation of basaltic magma is conducive to the mantle of sulfides in the effective access to melt. And the An component of dark part within plagioclase phenocryst zoning belong to bytownite (An 74%), and its may be a result of magma composition changes refreshment by basaltic magma injection. SHRIMP zircon U-Pb and LA-ICP-MS zircon U-Pb geochronology study showed that the intrusions and volcanic rocks from Duobuza porphyry copper deposit belong to early Cretaceous magma series (126~105Ma). The magma evolution series are as follows: the earliest diorite and diorite porphyrite → ore-bearing porphyry and barren grandiorite porphyry →basaltic andesite → diorite porphyrite → andesite → basaltic andesite, and magma component shows a evolution trend from intermediate to intermediate-acid to basic. Based on the field evidences, the formation age of high-Nb basalt may be the latest. The Ar-Ar geochronology of altered secondary biotite, K-feldspar and sericite shows that the main mineralization lasting a interval of about 4 Ma, the duration limit of whole magma-hydrothermal evolution of about 6 Ma, and possibly such a long duration limit may result in the formation of Duobuza super-large copper deposit. Moreover, tectonic diagram and trace element geochemistry of volcanic rocks and diorite from Duobuza porphyry copper deposit confirm that it formed in a continental margin arc environment. Zircon U-Pb age of volcanic rocks and porphyry fall in the range of 105~121Ma, and Duobuza porphyry copper deposit locating in the north of the Bangonghu- Nujiang suture zone, suggested that Neo-Tethys ocean still subducted northward at least early Cretaceous, and its closure time should be later than 105 Ma. Three major inclusion types and ten subtypes are distinguished from quartz phenocrysts and various quartz veins. Vapor generally coexisting with brine inclusions, suggest that fluid boiling may be the main ore-forming mechanism. Raman spectrums of fluid inclusions display that the content of vapor and liquid inclusion mainly contain water, and vapor occasionally contain a little CO2. In addition, the component of liquid inclusions mainly include Cl-, SO42-, Na+, K+, a small amount of Ca2+, F-; and Cl- and Na+ show good correlation. Vapor mainly contains water, a small amount of CO2, CH4 and C2H6 and so on. The daughter minerals identified by Laman spectroscopy and SEM include gypsum, chalcopyrite, halite, sylvite, rutile, potassium feldspar, Fe-Mn-chloride and other minerals, and ore-forming fluid belong to a complex hydrothermal system containing H2O-NaCl-KClFeCl2CaCl2. H and O isotopic analysis of quartz phenocryst, vein quartz, magnetite, chlorite and gypsum from all alteration zones show that the ore-forming fluid of Duobuza gold-rich porphyry copper deposit consisted mainly of magmatic water, without addition of meteric water. Duobuza gold-rich porphyry copper deposit formed by the primary magmatic fluid (600-950C), which has high oxidation, ultra-high salinity and metallogenic element-rich, exsolution direct from the magma, and it is representative of the typical orthomagmatic end member of the porphyry continuum. Moreover, the fluid evolution model of Duobuza gold-rich porphyry copper deposit has been established. Furthermore, two key factors for formation of large Au-rich porphyry copper deposit have been summed up, which are ore-forming fluids earlier separated from magma and high oxidation magma-mineralization fluid system.

Coronal loops heated by turbulence-driven Alfven waves

Li, Xing, Habbal, S. R., 'Coronal loops heated by turbulence-driven Alfven waves', The Astrophysical Journal, (2003) 598(2) pp.L125-L128 RAE2008

Nonlinear MHD stability of aluminium reduction cells

An industrial electrolysis cell used to produce primary aluminium is sensitive to waves at the interface of liquid aluminium and electrolyte. The interface waves are similar to stratified sea layers [1], but the penetrating electric current and the associated magnetic field are intricately involved in the oscillation process, and the observed wave frequencies are shifted from the purely hydrodynamic ones [2]. The interface stability problem is of great practical importance because the electrolytic aluminium production is a major electrical energy consumer, and it is related to environmental pollution rate. The stability analysis was started in [3] and a short summary of the main developments is given in [2]. Important aspects of the multiple mode interaction have been introduced in [4], and a widely used linear friction law first applied in [5]. In [6] a systematic perturbation expansion is developed for the fluid dynamics and electric current problems permitting reduction of the three-dimensional problem to a two dimensional one. The procedure is more generally known as “shallow water approximation” which can be extended for the case of weakly non-linear and dispersive waves. The Boussinesq formulation permits to generalise the problem for non-unidirectionally propagating waves accounting for side walls and for a two fluid layer interface [1]. Attempts to extend the electrolytic cell wave modelling to the weakly nonlinear case have started in [7] where the basic equations are derived, including the nonlinearity and linear dispersion terms. An alternative approach for the nonlinear numerical simulation for an electrolysis cell wave evolution is attempted in [8 and references there], yet, omitting the dispersion terms and without a proper account for the dissipation, the model can predict unstable waves growth only. The present paper contains a generalisation of the previous non linear wave equations [7] by accounting for the turbulent horizontal circulation flows in the two fluid layers. The inclusion of the turbulence model is essential in order to explain the small amplitude self-sustained oscillations of the liquid metal surface observed in real cells, known as “MHD noise”. The fluid dynamic model is coupled to the extended electromagnetic simulation including not only the fluid layers, but the whole bus bar circuit and the ferromagnetic effects [9].

Modulated electrostatic modes in pair plasmas: Modulational stability profile and envelope excitations

A pair plasma consisting of two types of ions, possessing equal masses and opposite charges, is considered. The nonlinear propagation of modulated electrostatic wave packets is studied by employing a two-fluid plasma model. Considering propagation parallel to the external magnetic field, two distinct electrostatic modes are obtained, namely a quasiacoustic lower moddfe and a Langmuir-like, as optic-type upper one, in agreement with experimental observations and theoretical predictions. Considering small yet weakly nonlinear deviations from equilibrium, and adopting a multiple-scale technique, the basic set of model equations is reduced to a nonlinear Schrodinger equation for the slowly varying electric field perturbation amplitude. The analysis reveals that the lower (acoustic) mode is stable and may propagate in the form of a dark-type envelope soliton (a void) modulating a carrier wave packet, while the upper linear mode is intrinsically unstable, and may favor the formation of bright-type envelope soliton (pulse) modulated wave packets. These results are relevant to recent observations of electrostatic waves in pair-ion (fullerene) plasmas, and also with respect to electron-positron plasma emission in pulsar magnetospheres. (c) 2006 American Institute of Physics.

Electrostatic mode envelope excitations in e-p-i plasmas - application in warm pair ion plasmas with a small fraction of stationary ions

The nonlinear propagation of amplitude-modulated electrostatic wavepackets in an electron-positron-ion (e-p-i) plasma is considered, by employing a two-fluid plasma model. Considering propagation parallel to the external magnetic field, two distinct electrostatic modes are obtained, namely a quasi-thermal acoustic-like lower mode and a Langmuir-like optic-type upper one. These results equally apply in warm pair ion ( e. g. fullerene) plasmas contaminated by a small fraction of stationary ions ( or dust), in agreement with experimental observations and theoretical predictions in pair plasmas. Considering small yet weakly nonlinear deviations from equilibrium, and adopting a multiple-scales perturbation technique, the basic set of model equations is reduced to a nonlinear Schrodinger (NLS) equation for the slowly varying electric field perturbation amplitude. The analysis reveals that the lower ( acoustic) mode is mostly stable for large wavelengths, and may propagate in the form of a dark-type envelope soliton ( a void) modulating a carrier wavepacket, while the upper linear mode is intrinsically unstable, and thus favours the formation of bright-type envelope soliton ( pulse) modulated wavepackets. The stability ( instability) range for the acoustic ( Langmuir-like optic) mode shifts to larger wavenumbers as the positive-to-negative ion temperature ( density) ratio increases. These results may be of relevance in astrophysical contexts, where e-p-i plasmas are encountered, and may also serve as prediction of the behaviour of doped ( or dust-contaminated) fullerene plasmas, in the laboratory.

Low frequency localized wavepackets in dusty plasmas with opposite charge polarity dust components

The reductive perturbation technique is employed to investigate the modulational instability of dust-acoustic (DA) waves propagating in a four-component dusty plasma. The dusty plasma consists of both positive- and negative-charge dust grains, characterized by a different mass, temperature and density, in addition to a background of Maxwellian electrons and ions. Relying on a multi-fluid plasma model and employing a multiple scales technique, a nonlinear Schrodinger type equation (NLSE) is obtained for the electric potential amplitude perturbation. The occurrence of localized electrostatic wavepackets is shown, in the form of oscillating structures whose modulated envelope is modelled as a soliton (or multi-soliton) solution of the NLSE. The DA wave characteristics, as well as the associated stability thresholds, are studied analytically and numerically. The relevance of these theoretical results with dusty plasmas observed in cosmic and laboratory environments is analysed in detail, by considering realistic multi-component plasma configurations observed in the polar mesosphere, as well as in laboratory experiments.

Nonlinear Dynamics of Rotating Multi-Component Pair Plasmas and e-p-i Plasmas

The propagation of small amplitude stationary profile nonlinear electrostatic excitations in a pair plasma is investigated, mainly drawing inspiration from experiments on fullerene pair-ion plasmas. Two distinct pair ion species are considered of opposite polarity and same mass, in addition to a massive charged background species, which is assumed to be stationary, given the frequency scale of interest. In the pair-ion context, the third species is thought of as a background defect (e.g. charged dust) component. On the other hand, the model also applies formally to electron-positron-ion (e-p-i) plasmas, if one neglects electron-positron annihilation. A two-fluid plasma model is employed, incorporating both Lorentz and Coriolis forces, thus taking into account the interplay between the gyroscopic (Larmor) frequency ?c and the (intrinsic) plasma rotation frequency O0. By employing a multi-dimensional reductive perturbation technique, a Zakharov-Kuznetsov (ZK) type equation is derived for the evolution of the electric potential perturbation. Assuming an arbitrary direction of propagation, with respect to the magnetic field, we derive the exact form of nonlinear solutions, and study their characteristics. A parametric analysis is carried out, as regards the effect of the dusty plasma composition (background number density), species temperature(s) and the relative strength of rotation to Larmor frequencies. It is shown that the Larmor and mechanical rotation affect the pulse dynamics via a parallel-to-transverse mode coupling diffusion term, which in fact diverges at ?c ? ±2O0. Pulses collapse at this limit, as nonlinearity fails to balance dispersion. The analysis is complemented by investigating critical plasma compositions, in fact near-symmetric (T- ˜ T+) “pure” (n- ˜ n+) pair plasmas, i.e. when the concentration of the 3rd background species is negligible, case in which the (quadratic) nonlinearity vanishes, so one needs to resort to higher order nonlinear theory. A modified ZK equation is derived and analyzed. Our results are of relevance in pair-ion (fullerene) experiments and also potentially in astrophysical environments, e.g. in pulsars.

Modeling relativistic soliton interactions in overdense plasmas: A perturbed nonlinear Schrodinger equation framework

We investigate the dynamics of localized solutions of the relativistic cold-fluid plasma model in the small but finite amplitude limit, for slightly overcritical plasma density. Adopting a multiple scale analysis, we derive a perturbed nonlinear Schrodinger equation that describes the evolution of the envelope of circularly polarized electromagnetic field. Retaining terms up to fifth order in the small perturbation parameter, we derive a self-consistent framework for the description of the plasma response in the presence of localized electromagnetic field. The formalism is applied to standing electromagnetic soliton interactions and the results are validated by simulations of the full cold-fluid model. To lowest order, a cubic nonlinear Schrodinger equation with a focusing nonlinearity is recovered. Classical quasiparticle theory is used to obtain analytical estimates for the collision time and minimum distance of approach between solitons. For larger soliton amplitudes the inclusion of the fifth-order terms is essential for a qualitatively correct description of soliton interactions. The defocusing quintic nonlinearity leads to inelastic soliton collisions, while bound states of solitons do not persist under perturbations in the initial phase or amplitude

Community strategies in England: reshaping spaces of governance?

The preparation of Community Strategies (CS) has been required of LSPs and Local Authorities in England since the passing of the Local Government Act 2000. This paper examines the process and content of two Community Strategies in southern England as part of an ongoing project to understand their impact and explore ways that CSs may be carried through in a meaningful and effective manner. The paper concludes that the two CSs studied illustrate the challenge faced by LSPs in producing Strategies that are meaningful, inclusive and which follow the spirit of the government CS guidance. LAs and LSPs are also posed with a difficult challenge of seeing through an implicitly required transition from a traditional representative democratic structure/process with a more fluid participatory model. Thus we detect that at least two forms of conflict may arise – firstly with elected councillors threatened by a loss of power and secondly between communities and the LAs who are encouraged to problematise local policy and service delivery in the context of limited resource availability.

Effect of operating parameters and cleaning on the performance of ceramic membranes treating partially clarified sugar cane juice

The performance of ceramic membranes with pore sizes of 0.05 and 0.10 mm in purifying limed and partially clarified sugar cane juice was investigated under different operating conditions. From various operating conditions and strategies, switching off the permeate for 5 seconds for every 5 minutes (S5sT5 m) by an automated control valve provided higher flux. From the three pH experiments conducted on the 0.05 mm membrane, the best performance was observed at a pH of 7.5. Amongst the four fouling models tested, the cake filtration model fitted the performance of both membranes with higher accuracy at a transmembrane pressure of 0.5 bar. Filtering the cane juice through the membrane reduced the turbidity by 99.7%, color by 15%, and the starch concentration by 80% as well as increased the purity by 1.4%. The effective cleaning chemical composition from experimental results showed that 1% NaOH and 3000 ppm NaOCl solution performed the best but only for the experiments that were treating limed and partially clarified juice at pH 7.5.

Cocoon of the silkworm Antheraea pernyi as an example of a thermally insulating biological interface.

Biological materials are hierarchically organized complex composites, which embrace multiple practical functionalities. As an example, the wild silkworm cocoon provides multiple protective functions against environmental and physical hazards, promoting the survival chance of moth pupae that resides inside. In the present investigation, the microstructure and thermal property of the Chinese tussah silkworm (Antheraea pernyi) cocoon in both warm and cold environments under windy conditions have been studied by experimental and numerical methods. A new computational fluid dynamics model has been developed according to the original fibrous structure of the Antheraea pernyi cocoon to simulate the unique heat transfer process through the cocoon wall. The structure of the Antheraea pernyi cocoon wall can promote the disorderness of the interior air, which increases the wind resistance by stopping most of the air flowing into the cocoon. The Antheraea pernyi cocoon is wind-proof due to the mineral crystals deposited on the outer layer surface and its hierarchical structure with low porosity and high tortuosity. The research findings have important implications to enhancing the thermal function of biomimetic protective textiles and clothing.

Modelagem analítica e experimental da filtração em meios porosos

Deep bed filtration occurs in several industrial and environmental processes like water filtration and soil contamination. In petroleum industry, deep bed filtration occurs near to injection wells during water injection, causing injectivity reduction. It also takes place during well drilling, sand production control, produced water disposal in aquifers, etc. The particle capture in porous media can be caused by different physical mechanisms (size exclusion, electrical forces, bridging, gravity, etc). A statistical model for filtration in porous media is proposed and analytical solutions for suspended and retained particles are derived. The model, which incorporates particle retention probability, is compared with the classical deep bed filtration model allowing a physical interpretation of the filtration coefficients. Comparison of the obtained analytical solutions for the proposed model with the classical model solutions allows concluding that the larger the particle capture probability, the larger the discrepancy between the proposed and the classical models