The Equation of State and Petrogenesis of Komatiite

(1) Equation of State of Komatiite

The equation of state (EOS) of a molten komatiite (27 wt% MgO) was detennined in the 5 to 36 GPa pressure range via shock wave compression from 1550°C and 0 bar. Shock wave velocity, U_S, and particle velocity, U_P, in km/s follow the linear relationship U_S = 3.13(±0.03) + 1.47(±0.03) U_P. Based on a calculated density at 1550°C, 0 bar of 2.745±0.005 glee, this U_S-U_P relationship gives the isentropic bulk modulus K_S = 27.0 ± 0.6 GPa, and its first and second isentropic pressure derivatives, K'_S = 4.9 ± 0.1 and K"_S = -0.109 ± 0.003 GPa^-1.

The calculated liquidus compression curve agrees within error with the static compression results of Agee and Walker [1988a] to 6 GPa. We detennine that olivine (FO₉₄) will be neutrally buoyant in komatiitic melt of the composition we studied near 8.2 GPa. Clinopyroxene would also be neutrally buoyant near this pressure. Liquidus garnet-majorite may be less dense than this komatiitic liquid in the 20-24 GPa interval, however pyropic-garnet and perovskite phases are denser than this komatiitic liquid in their respective liquidus pressure intervals to 36 GPa. Liquidus perovskite may be neutrally buoyant near 70 GPa.

At 40 GPa, the density of shock-compressed molten komatiite would be approximately equal to the calculated density of an equivalent mixture of dense solid oxide components. This observation supports the model of Rigden et al. [1989] for compressibilities of liquid oxide components. Using their theoretical EOS for liquid forsterite and fayalite, we calculate the densities of a spectrum of melts from basaltic through peridotitic that are related to the experimentally studied komatiitic liquid by addition or subtraction of olivine. At low pressure, olivine fractionation lowers the density of basic magmas, but above 14 GPa this trend is reversed. All of these basic to ultrabasic liquids are predicted to have similar densities at 14 GPa, and this density is approximately equal to the bulk (PREM) mantle. This suggests that melts derived from a peridotitic mantle may be inhibited from ascending from depths greater than 400 km.

The EOS of ultrabasic magmas was used to model adiabatic melting in a peridotitic mantle. If komatiites are formed by >15% partial melting of a peridotitic mantle, then komatiites generated by adiabatic melting come from source regions in the lower transition zone (≈500-670 km) or the lower mantle (>670 km). The great depth of incipient melting implied by this model, and the melt density constraint mentioned above, suggest that komatiitic volcanism may be gravitationally hindered. Although komatiitic magmas are thought to separate from their coexisting crystals at a temperature =200°C greater than that for modern MORBs, their ultimate sources are predicted to be diapirs that, if adiabatically decompressed from initially solid mantle, were more than 700°C hotter than the sources of MORBs and derived from great depth.

We considered the evolution of an initially molten mantle, i.e., a magma ocean. Our model considers the thermal structure of the magma ocean, density constraints on crystal segregation, and approximate phase relationships for a nominally chondritic mantle. Crystallization will begin at the core-mantle boundary. Perovskite buoyancy at > 70 GPa may lead to a compositionally stratified lower mantle with iron-enriched mangesiowiistite content increasing with depth. The upper mantle may be depleted in perovskite components. Olivine neutral buoyancy may lead to the formation of a dunite septum in the upper mantle, partitioning the ocean into upper and lower reservoirs, but this septum must be permeable.

(2) Viscosity Measurement with Shock Waves

We have examined in detail the analytical method for measuring shear viscosity from the decay of perturbations on a corrugated shock front The relevance of initial conditions, finite shock amplitude, bulk viscosity, and the sensitivity of the measurements to the shock boundary conditions are discussed. The validity of the viscous perturbation approach is examined by numerically solving the second-order Navier-Stokes equations. These numerical experiments indicate that shock instabilities may occur even when the Kontorovich-D'yakov stability criteria are satisfied. The experimental results for water at 15 GPa are discussed, and it is suggested that the large effective viscosity determined by this method may reflect the existence of ice VII on the Rayleigh path of the Hugoniot This interpretation reconciles the experimental results with estimates and measurements obtained by other means, and is consistent with the relationship of the Hugoniot with the phase diagram for water. Sound waves are generated at 4.8 MHz at in the water experiments at 15 GPa. The existence of anelastic absorption modes near this frequency would also lead to large effective viscosity estimates.

(3) Equation of State of Molybdenum at 1400°C

Shock compression data to 96 GPa for pure molybdenum, initially heated to 1400°C, are presented. Finite strain analysis of the data gives a bulk modulus at 1400°C, K'_S. of 244±2 GPa and its pressure derivative, K'_OS of 4. A fit of shock velocity to particle velocity gives the coefficients of U_S = C_O+S U_P to be C_O = 4.77±0.06 km/s and S = 1.43±0.05. From the zero pressure sound speed, C_O, a bulk modulus of 232±6 GPa is calculated that is consistent with extrapolation of ultrasonic elasticity measurements. The temperature derivative of the bulk modulus at zero pressure, θK_OSθT|_P, is approximately -0.012 GPa/K. A thermodynamic model is used to show that the thermodynamic Grüneisen parameter is proportional to the density and independent of temperature. The Mie-Grüneisen equation of state adequately describes the high temperature behavior of molybdenum under the present range of shock loading conditions.

溶胶-凝胶条波导Sagnac环的乙醇蒸气传感特性

采用有机/无机混合溶胶-凝胶法制作条形光波导,并将条波导接入光纤Sagnac 环中,测量了输出光功率随环境气氛中乙醇蒸气体积分数变化的特性,表明在实验研究的范围内,输出信号与乙醇蒸气体积分数呈正弦变化。根据Sagnac环结构输出特性的基本关系,反映了溶胶-凝胶条波导在乙醇蒸气气氛下产生了双折射效应。观察到双折射相移与乙醇体积分数的亚线性关系。对实验数据拟合,计算了偏振相移的线性项和二次项系数,得到所制备的条波导的双折射对乙醇体积分数的响应为Δn≈4.4×10-2。测量了信号变化的时间演变特性,典型的上升和

Biossorção de vanádio pela macroalga marinha Sargassum filipendula

Muitos biossorventes naturais têm sido pesquisados por possuírem baixo custo e apresentarem propriedades ligantes, como é o caso da macroalga marinha Sargassum filipendula (S. filipendula) que vem sendo utilizada como material biossorvente no processo de biossorção de metais pesados. No presente trabalho a alga marrom foi utilizada para estudos realizados em batelada, onde se determinou o pH ideal de biossorção de vanádio, a relação sólido/líquido ideal e a importância da velocidade de agitação. O estudo cinético e de equilíbrio dos íons metálicos também foram realizados em bateladas nas seguintes condições de ensaio: (1) 0,10 g de biomassa, 25,0 mL de solução de vanádio igual a 18,0 mg L-1, temperatura 25,0 C e 150 rpm de agitação; (2) 0,10 g de biomassa, 25,0 mL de solução de vanádio igual a 36,0 mg L-1, temperatura 25,0 C e 150 rpm de agitação. A melhor condição para biossorção de vanádio foi encontrada para 36,0 mg L-1 e pH= 2,0. O estudo cinético de biossorção de vanádio mostrou que o modelo de segunda ordem descreve melhor os dados experimentais em 36,0 mg L-1 (R2= 0,9825). O estudo de equilíbrio mostrou um perfil crescente de remoção de vanádio. A melhor eficiência de captação dos íons de vanádio foi de 61,0 % para Co= 40,0 mg L-1 em pH= 2,0. Os dados experimentais da isoterma de vanádio mostraram-se mais adequados ao modelo de Langmuir para pH= 2,0, Os parâmetros de equilíbrio calculados a partir do modelo de Langmuir (b, qmax ) 0,009 e 43,3 mg/g, respectivamente, corroboram melhor para a interpretação dos resultados quando comparados com o modelo de Freundlich (kF, n) 1,56 e 2,41, visto que o coeficiente de correlação é maior para Langmuir

Kinetics and mechanisms of adsorption of Escherichia coli bacteriophage T_4 to activated carbon

A study was conducted on the adsorption of Escherichia coli bacteriophage T₄ to activated carbon. Preliminary adsorption experiments were also made with poliovirus Type III. The effectiveness of such adsorbents as diatomaceous earth, Ottawa sand, and coconut charcoal was also tested for virus adsorption.

The kinetics of adsorption were studied in an agitated solution containing virus and carbon. The mechanism of attachment and site characteristics were investigated by varying pH and ionic strength and using site-blocking reagents.

Plaque assay procedures were developed for bacteriophage T₄ on Escherichia coli cells and poliovirus Type III on monkey kidney cells. Factors influencing the efficiency of plaque formation were investigated.

The kinetics of bacteriophage T₄ adsorption to activated carbon can be described by a reversible second-order equation. The reaction order was first order with respect to both virus and carbon concentration. This kinetic representation, however, is probably incorrect at optimum adsorption conditions, which occurred at a pH of 7.0 and ionic strength of 0.08. At optimum conditions the adsorption rate was satisfactorily described by a diffusion-limited process. Interpretation of adsorption data by a development of the diffusion equation for Langmuir adsorption yielded a diffusion coefficient of 12 X 10^-8 cm²/sec for bacteriophage T₄. This diffusion coefficient is in excellent agreement with the accepted value of 8 X 10^-8 cm²/sec. A diffusion-limited theory may also represent adsorption at conditions other than the maximal. A clear conclusion on the limiting process cannot be made.

Adsorption of bacteriophage T₄ to activated carbon obeys the Langmuir isotherm and is thermodynamically reversible. Thus virus is not inactivated by adsorption. Adsorption is unimolecular with very inefficient use of the available carbon surface area. The virus is probably completely excluded from pores due to its size.

Adsorption is of a physical nature and independent of temperature. Attraction is due to electrostatic forces between the virus and carbon. Effects of pH and ionic strength indicated that carboxyl groups, amino groups, and the virus's tail fibers are involved in the attachment of virus to carbon. The active sites on activated carbon for adsorption of bacteriophage T₄ are carboxyl groups. Adsorption can be completely blocked by esterifying these groups.

I. The rate and mechanism of the vapor-phase reaction between water and parts-per-million concentrations of nitrogen dioxide. II. The rate and mechanism of the air oxidation of parts-per-million concentrations of nitric oxide in the presence of water vapor.

Part I

A study of the thermal reaction of water vapor and parts-per-million concentrations of nitrogen dioxide was carried out at ambient temperature and at atmospheric pressure. Nitric oxide and nitric acid vapor were the principal products. The initial rate of disappearance of nitrogen dioxide was first order with respect to water vapor and second order with respect to nitrogen dioxide. An initial third-order rate constant of 5.5 (± 0.29) x 10⁴ liter² mole^-2 sec^-1 was found at 25˚C. The rate of reaction decreased with increasing temperature. In the temperature range of 25˚C to 50˚C, an activation energy of -978 (± 20) calories was found.

The reaction did not go to completion. From measurements as the reaction approached equilibrium, the free energy of nitric acid vapor was calculated. This value was -18.58 (± 0.04) kilocalories at 25˚C.

The initial rate of reaction was unaffected by the presence of oxygen and was retarded by the presence of nitric oxide. There were no appreciable effects due to the surface of the reactor. Nitric oxide and nitrogen dioxide were monitored by gas chromatography during the reaction.

Part II

The air oxidation of nitric oxide, and the oxidation of nitric oxide in the presence of water vapor, were studied in a glass reactor at ambient temperatures and at atmospheric pressure. The concentration of nitric oxide was less than 100 parts-per-million. The concentration of nitrogen dioxide was monitored by gas chromatography during the reaction.

For the dry oxidation, the third-order rate constant was 1.46 (± 0.03) x 10⁴ liter² mole^-2 sec^-1 at 25˚C. The activation energy, obtained from measurements between 25˚C and 50˚C, was -1.197 (±0.02) kilocalories.

The presence of water vapor during the oxidation caused the formation of nitrous acid vapor when nitric oxide, nitrogen dioxide and water vapor combined. By measuring the difference between the concentrations of nitrogen dioxide during the wet and dry oxidations, the rate of formation of nitrous acid vapor was found. The third-order rate constant for the formation of nitrous acid vapor was equal to 1.5 (± 0.5) x 10⁵ liter² mole^-2 sec^-1 at 40˚C. The reaction rate did not change measurably when the temperature was increased to 50˚C. The formation of nitric acid vapor was prevented by keeping the concentration of nitrogen dioxide low.

Surface effects were appreciable for the wet tests. Below 35˚C, the rate of appearance of nitrogen dioxide increased with increasing surface. Above 40˚C, the effect of surface was small.

Analysis of dispersion compensation for position-dependence in externally modulated CATV lightwave systems by using chirped fiber grating

The dispersion compensation characteristics of the chirped fiber grating (CFG) for different dispersion compensation positions are analyzed in externally modulated cable television (CATV) lightwave system and the analytic expression of the composite second order (CSO) distortion is derived. The analyses give a reasonable explanation for the position-dependent effect of CFG dispersion compensator, which was found in practical systems. Moreover, the theoretical result is also verified by an experiment. It is believed that the theory will be helpful in designing optical CATV fiber links with nodes at proper positions both for intensity amplification and dispersion compensation.

Range of validity of the method of averaging

Sufficient conditions are derived for the validity of approximate periodic solutions of a class of second order ordinary nonlinear differential equations. An approximate solution is defined to be valid if an exact solution exists in a neighborhood of the approximation.

Two classes of validity criteria are developed. Existence is obtained using the contraction mapping principle in one case, and the Schauder-Leray fixed point theorem in the other. Both classes of validity criteria make use of symmetry properties of periodic functions, and both classes yield an upper bound on a norm of the difference between the approximate and exact solution. This bound is used in a procedure which establishes sufficient stability conditions for the approximated solution.

Application to a system with piecewise linear restoring force (bilinear system) reveals that the approximate solution obtained by the method of averaging is valid away from regions where the response exhibits vertical tangents. A narrow instability region is obtained near one-half the natural frequency of the equivalent linear system. Sufficient conditions for the validity of resonant solutions are also derived, and two term harmonic balance approximate solutions which exhibit ultraharmonic and subharmonic resonances are studied.

Some experiments on solving multistage stochastic mixed 0-1 programs with time stochastic dominance constraints

In this work we extend to the multistage case two recent risk averse measures for two-stage stochastic programs based on first- and second-order stochastic dominance constraints induced by mixed-integer linear recourse. Additionally, we consider Time Stochastic Dominance (TSD) along a given horizon. Given the dimensions of medium-sized problems augmented by the new variables and constraints required by those risk measures, it is unrealistic to solve the problem up to optimality by plain use of MIP solvers in a reasonable computing time, at least. Instead of it, decomposition algorithms of some type should be used. We present an extension of our Branch-and-Fix Coordination algorithm, so named BFC-TSD, where a special treatment is given to cross scenario group constraints that link variables from different scenario groups. A broad computational experience is presented by comparing the risk neutral approach and the tested risk averse strategies. The performance of the new version of the BFC algorithm versus the plain use of a state-of-the-artMIP solver is also reported.

Optimum linear coding for additive noise systems using information feedback

The feedback coding problem for Gaussian systems in which the noise is neither white nor statistically independent between channels is formulated in terms of arbitrary linear codes at the transmitter and at the receiver. This new formulation is used to determine a number of feedback communication systems. In particular, the optimum linear code that satisfies an average power constraint on the transmitted signals is derived for a system with noiseless feedback and forward noise of arbitrary covariance. The noisy feedback problem is considered and signal sets for the forward and feedback channels are obtained with an average power constraint on each. The general formulation and results are valid for non-Gaussian systems in which the second order statistics are known, the results being applicable to the determination of error bounds via the Chebychev inequality.

Energy loss of light ions channeling in silicon

Experimental measurements of rate of energy loss were made for protons of energy .5 to 1.6 MeV channeling through 1 μm thick silicon targets along the <110>, <111>, and <211> axial directions, and the {100}, {110}, {111}, and {211} planar directions. A .05% resolution automatically controlled magnetic spectrometer was used. The data are presented graphically along with an extensive summary of data in the literature. The data taken cover a wider range of channels than has previously been examined, and are in agreement with the data of F. Eisen, et al., Radd. Eff. 13, 93 (1972).

The theory in the literature for channeling energy loss due to interaction with local electrons, core electrons, and distant valence electrons of the crystal atoms is summarized. Straggling is analyzed, and a computer program which calculates energy loss and straggling using this theory and the Moliere approximation to the Thomas Fermi potential, V_TF, and the detailed silicon crystal structure is described. Values for the local electron density Z_loc in each of the channels listed above are extracted from the data by graphical matching of the experimental and computer results.

Zeroth and second order contributions to Z_loc as a function of distance from the center of the channel were computed from ∇²V_TF = 4πρ for various channels in silicon. For data taken in this work and data of F. Eisen, et al., Rad. Eff. 13, 93 (1972), the calculated zeroth order contribution to Z_loc lies between the experimentally extracted Z_loc values obtained by using the peak and the leading edge of the transmission spectra, suggesting that the observed straggling is due both to statistical fluctuations and to path variation.

Theoretical investigation of turbulent boundary layer over a wavy surface

The important features of the two-dimensional incompressible turbulent flow over a wavy surface of wavelength comparable with the boundary layer thickness are analyzed.

A turbulent field method using model equation for turbulent shear stress similar to the scheme of Bradshaw, Ferriss and Atwell (1967) is employed with suitable modification to cover the viscous sublayer. The governing differential equations are linearized based on the small but finite amplitude to wavelength ratio. An orthogonal wavy coordinate system, accurate to the second order in the amplitude ratio, is adopted to avoid the severe restriction to the validity of linearization due to the large mean velocity gradient near the wall. Analytic solution up to the second order is obtained by using the method of matched-asymptotic-expansion based on the large Reynolds number and hence the small skin friction coefficient.

In the outer part of the layer, the perturbed flow is practically "inviscid." Solutions for the velocity, Reynolds stress and also the wall pressure distributions agree well with the experimental measurement. In the wall region where the perturbed Reynolds stress plays an important role in the process of momentum transport, only a qualitative agreement is obtained. The results also show that the nonlinear second-order effect is negligible for amplitude ratio of 0.03. The discrepancies in the detailed structure of the velocity, shear stress, and skin friction distributions near the wall suggest modifications to the model are required to describe the present problem.

Energy inequalities and error estimates for axisymmetric torsion of thin elastic shells of revolution

The problem motivating this investigation is that of pure axisymmetric torsion of an elastic shell of revolution. The analysis is carried out within the framework of the three-dimensional linear theory of elastic equilibrium for homogeneous, isotropic solids. The objective is the rigorous estimation of errors involved in the use of approximations based on thin shell theory.

The underlying boundary value problem is one of Neumann type for a second order elliptic operator. A systematic procedure for constructing pointwise estimates for the solution and its first derivatives is given for a general class of second-order elliptic boundary-value problems which includes the torsion problem as a special case.

The method used here rests on the construction of “energy inequalities” and on the subsequent deduction of pointwise estimates from the energy inequalities. This method removes certain drawbacks characteristic of pointwise estimates derived in some investigations of related areas.

Special interest is directed towards thin shells of constant thickness. The method enables us to estimate the error involved in a stress analysis in which the exact solution is replaced by an approximate one, and thus provides us with a means of assessing the quality of approximate solutions for axisymmetric torsion of thin shells.

Finally, the results of the present study are applied to the stress analysis of a circular cylindrical shell, and the quality of stress estimates derived here and those from a previous related publication are discussed.

Biossorção de tório com emprego de Sargassum filipendula

Com o aumento do rigor para descarte de efluentes contaminado com metais pesados, as pesquisas têm se intensificado na busca de métodos de remoção, que tragam bons resultados de captação dos metais, aliado a um baixo custo. O uso de biomassas como bactérias, fungos e algas marinhas como material adsorvente, tem se apresentado como uma alternativa, principalmente para soluções com baixo teor de metais. Neste trabalho, a alga marinha Sargassum filipendula foi avaliada na sua capacidade de remoção do metal tório de uma solução sintética e do efluente dos laboratórios de análises ambientais do IRD. A cinética desta biossorção foi estudada em regime de batelada e o equilíbrio foi alcançado com 180 min de reação. Dois modelos cinéticos foram utilizados nesta avaliação, um de primeira ordem e um modelo de pseudo segunda ordem, tendo o modelo de segunda ordem apresentado um melhor ajuste dos dados. Na avaliação da capacidade máxima de captação do tório pela biomassa marinha em regime de batelada, foi construída a isoterma que apresentou um perfil crescente na captação alcançando um máximo de 2,59 mol/g. Os modelos de Langmuir e Freundlich foram utilizados para ajustar os dados da isoterma, tendo apresentado maior correlação com os dados o modelo de Langmuir, resultando num valor de captação máxima calculado pelo modelo de 2,92 mol/g. A capacidade de remoção do metal tório da alga Sargassum filipendula também foi avaliada em regime contínuo. Um estudo de altura crítica de leito foi realizado preenchendo-se uma coluna com diferentes massas de alga correspondendo a diferentes alturas de leito. A concentração de tório residual foi quantificada na solução de saída e a menor concentração na saída foi alcançada com 40 cm de leito ou 96 g de biomassa. Após este estudo um sistema contínuo com duas colunas, com 96 g de biomassa cada, e uma bomba peristáltica foi utilizado para o tratamento do efluente real do IRD, contendo não apenas o metal tório, mas outros metais como urânio, cálcio, cromo, ferro, chumbo, etc. Cento e cinco litros de efluente foram tratados numa concentração de 6 mg/L, a concentração do efluente de saída foi de 3,75 mg/L. A caracterização do efluente quanto aos metais presentes foi realizada em ICP-MS, os resultados demonstraram que não houve competição pelos sítios ligantes da biomassa entre o tório e os demais metais. Além disso, alguns metais como cálcio, ferro e magnésio, tiveram um aumento na concentração de saída indicando a presença do mecanismo de troca iônica na biossorção do tório por Sargassum filipendula

Analysis of dispersion compensation by chirped fiber grating in directly and externally modulated CATV systems

The dispersion compensation effect of the chirped fiber grating (CFG) is analyzed theoretically, and analytic expressions are derived for composite second-order (CSO) distortion in analog modulated sub-carrier multiplexed (AM-SCM) cable television (CATV) systems with externally and directly modulated transmitters. Simulations are given for the two kinds of modulations and for standard single mode fiber and non-zero dispersion shift fiber (NZDSF) systems. The results show that CFG could be used as a dispersion compensator in directly modulated systems, but its dispersion coefficient should be adjusted much more precisely than the externally modulated system. The requirements for the NZDSF system could be loosened much. It is proposed that directly modulated source may be used as a transmitter in CATV systems combined with tunable CFG dispersion compensator being adjusted precisely, which may be more cost-effective than externally modulation technology. (c) 2006 Elsevier GmbH. All rights reserved.

Synthesis and characterization of photorefractive polymer based on chemically hybridized CdS-PVK nanocomposite with a new azo chromophore

We report an organic/inorganic polymer composite based on the chemically hybridized photoconductor CdS-PVK nanocomposite doped with a new second-order optically nonlinear chromophore 1-n-butoxy-2-methyl-(4-p-nitrophenylazo)benzene (BMNPAB) and plasticizer 9-ethyl-carbazole (ECZ) to manifest a photorefractive (PR) effect. A detailed description of the synthesis and characterization of BMNPAB is presented. The poled film including PVK-10-CdS nanocomposite and BMNPAB exhibits a high second harmonic generation (SHG) coefficient of 31 pm/V The photoconductivity of PVK-CdS nanocomposite also was studied here. Two-beam coupling experiment clearly indicated an asymmetric optical energy exchange between two beams on the polymer composite at zero electrical field, and the two-beam coupling gain of 50.0 cm(-1) and diffraction efficiencv of 4.2% were obtained at 647.1 nm wavelength. (C) 2007 Elsevier Ltd. All rights reserved.