Study of three elastic-plastic constitutive models by non-proportional finite deformations of OFHC copper

Experimental stress-strain data of OFHC copper first under torsion to 13% and then under torsion-tension to about 10% are used to study the characteristics of three elastic-plastic constitutive models: Chaboche's super-positional nonlinear model, Dafalias and Popov's two surface model and Watanabe and Atluri's version of the endochronic model. The three models, originally oriented for infinitesimal deformation, have been extended for finite deformation. The results show (a) the Mises-type yield surface used in the three models brings about significant departure of the predictions from the experimental data; (b) Chaboche's and Dafalias' models are easier than Watanabe and Atluri's model in determining the material parameters in them, and (c) Chaboche's and Watanabe & Atluri's models produce almost the same prediction to the data, while Dafalias' model cannot accurately predict the plastic deformations when a loading path changes in its direction. Copyright (C) 1996 Elsevier Science Ltd

Influence of liquid bridge volume on the onset of oscillation in floating-zone convection .2. Numerical-simulation

The onset of oscillation in the floating zone convection driven by the gradient of surface tension was studied numerically for an unsteady and two-dimensional model, and studies were concentrated on the influence of liquid bridge volume on the onset of oscillation in comparison with the experimental results in the Paper I. The numerical results agree with the experimental ones presented in the previous paper, in which the distributions of critical applied temperature difference depending on the volume of liquid bridge and a gap range of liquid volume in marginal stability curve were obtained.

Influence of convections on dopant segregation in floating-zone crystal-growth system

The steady and axisymmetric crystal growth process of floating zone model was studied numerically to concern with the influence of convection and phase change on effective segregation. An iteration method of numerical simulation considering both thermocapillary and buoyancy effects for GaAs crystal growth gave the effective segregation coefficient, which was compared with the space experiment of GaAs on board the Chinese recoverable satellite. The calculated segregation coefficient of a two-dimensional model was found to be smaller than the one suggested by space experiment with the simplified assumption of an one-dimensional model.

Mechanical modeling of grain boundary sliding of polycrystals

Using a variational method, a general three-dimensional solution to the problem of a sliding spherical inclusion embedded in an infinite anisotropic medium is presented in this paper. The inclusion itself is also a general anisotropic elastic medium. The interface is treated as a thin interface layer with interphase anisotropic properties. The displacements in the matrix and the inclusion are expressed as polynomial series of the cartesian coordinate components. Using the virtual work principle, a set of linear algebraic equations about unknown coefficients are obtained. Then the general sliding spherical inclusion problem is accurately solved. Based on this solution, a self-consistent method for sliding polycrystals is proposed. Combining this with a two-dimensional model of an aggregate polycrystal, a systematic analysis of the mechanical behaviour of sliding polycrystals is given in detail. Numerical results are given to show the significant effect of grain boundary sliding on the overall mechanical properties of aggregate polycrystals.

Asymmetric structure of an isolated flux arch

This paper extends two-dimensional model of symmetric magnetostatic flux arches confined in stratified atmospheres (Zhang and Hu, 1992, 1993) to asymmetric models. Numerical results show that the flux structure is influenced greatly by the boundary condition of magnetic field, the force-free factor, the atmospheric pressure distribution and the position of footpoints (especially the width ratio of outlet to entrance, which differs from symmetric case).

Monte Carlo Valuation of natural gas investments

This paper deals with the valuation of energy assets related to natural gas. In particular, we evaluate a baseload Natural Gas Combined Cycle (NGCC) power plant and an ancillary instalation, namely a Liquefied Natural Gas (LNG) facility, in a realistic setting; specifically, these investments enjoy a long useful life but require some non-negligible time to build. Then we focus on the valuation of several investment options again in a realistic setting. These include the option to invest in the power plant when there is uncertainty concerning the initial outlay, or the option's time to maturity, or the cost of CO2 emission permits, or when there is a chance to double the plant size in the future. Our model comprises three sources of risk. We consider uncertain gas prices with regard to both the current level and the long-run equilibrium level; the current electricity price is also uncertain. They all are assumed to show mean reversion. The two-factor model for natural gas price is calibrated using data from NYMEX NG futures contracts. Also, we calibrate the one-factor model for electricity price using data from the Spanish wholesale electricity market, respectively. Then we use the estimated parameter values alongside actual physical parameters from a case study to value natural gas plants. Finally, the calibrated parameters are also used in a Monte Carlo simulation framework to evaluate several American-type options to invest in these energy assets. We accomplish this by following the least squares MC approach.

Chemical equilibrium analysis and CFD simulation of coal combustion and NOx formation

A full two-fluid model of reacting gas-particle flows with an algebraic unified second-order moment (AUSM) turbulence-chemistry model is used to simulate Beijing coal combustion and NOx formation. The sub-models are the k-epsilon-kp two-phase turbulence model, the EBU-Arrhenius volatile and CO combustion model, the six-flux radiation model, coal devolatilization model and char combustion model. The blocking effect on NOx formation is discussed. In addition, the chemical equilibrium analysis is used to predict NOx concentration at different temperature. Results of CID simulation and chemical equilibrium analysis show that, optimizing air dynamic parameters can delay the NOx formation and decrease NOx emission, but it is effective only in a restricted range. In order to decrease NOx emission near to zero, the re-burning or other chemical methods must be used.

Influence of Evaporation on Marangoni-Benard Instability

Evaporative convection and instability give rise to both scientific and technological interests. Practically, a number of the industrial applications such as thin-film evaporators, boiling technologies and heat pipes concern with the evaporation process of which through the vapor-liquid interface the heat and mass transfer occur. From a physical viewpoint, one of interesting questions is the mechanisms of convection instability in thin-liquid layers induced by the coupling of evaporation phenomenon and Marangoni effect at the mass exchanged interface. Classical theories, including Rayleigh’s and Pearson’s, have only successfully explained convection in a liquid layer heated from below without evaporation. However these theories are unable to explain the convection in an evaporating thin layer, especially liquid layer is cooled from below. In present paper, a new two-sided model is put forward rather than the one-sided model in previous works. In previous works, the vapor is treated as passive gas and dynamics of vapor has been ignored. In this case, the vapor liquid system can be described by one-sided model. In our two-sided model, the dynamics of vapor should be considered. Linear instability analysis of the Marangoni-Bénard convection in the two-layer system with an evaporation interface is performed. We define a new evaporating Biot number which is different from the Biot number in one-sided model and obtain the curves of critical Marangoni number versus wave number. In our theoretical results, the Biot number and the evaporating velocity play a major role in the stability of the vapor-liquid system.

New Numerical Algorithms in SUPER CE/SE and Their Applications in Explosion Mechanics

The new numerical algorithms in SUPER/CESE and their applications in explosion mechanics are studied. The researched algorithms and models include an improved CE/SE (space-time Conservation Element and Solution Element) method, a local hybrid particle level set method, three chemical reaction models and a two-fluid model. Problems of shock wave reflection over wedges, explosive welding, cellular structure of gaseous detonations and two-phase detonations in the gas-droplet system are simulated by using the above-mentioned algorithms and models. The numerical results reveal that the adopted algorithms have many advantages such as high numerical accuracy, wide application field and good compatibility. The numerical algorithms presented in this paper may be applied to the numerical research of explosion mechanics.

Growth, movement, and attrition of northern bluefin tuna, Thunnus thynnus, in the Pacific Ocean, as determined by tagging

ENGLISH: The growth of northern bluefin tuna is described by a two-stanza model. For fish between 191 and 564 mm in length the Gompertz curve, with values of 581 mm and 4.32 for Loo and K (annual), respectively, is used. The fish between 564 and 1530 mm grow linearly, at the rate of 0.709 mm per day. Age-O fish tagged and released in the western Pacific Ocean have been recaptured in the western, central, and eastern Pacific. The minimum time between release in the western Pacific and recapture in the eastern Pacific is 215 days. Older fish, mostly Land 2-year olds, tagged and released in the eastern Pacific have been recaptured in the eastern and western Pacific. The minimum time between release in eastern Pacific and recapture in the western Pacific is 674 days. The coefficient of natural mortality is estimated from data on growth and ambient temperature to be 0.276 on an annual basis, with 90-percent confidence limits of 0.161 and 0.47L Spawning of northern bluefin takes place only in the western Pacific. Some of the juveniles migrate to the eastern Pacific, where they reside for several months to several years before returning to the western Pacific. The portion of fish which migrate to the eastern Pacific varies among years, and this appears to be an important cause of the annual variation in the catches in the eastern Pacific Ocean. SPANISH: El crecimiento del atún aleta azul del norte es descrito por un modelo de dos estadios. Para los peces de entre 191 y 564 mm de talla se usa la curva de Gompertz, con valores de 581 mm y 4.32 para Loo y K (anual), respectivamente. Los peces de entre 564 y 1530 mm crecen de forma lineal, a 0.709 mm por día. Peces de edad Omarcados y liberados en el Pacífico occidental han sido recapturados en el Pacífico occidental, central, y oriental. La demora mínima entre la liberación en el Pacífico occidental y la recaptura en el Pacífico oriental es de 215 días. Peces mayores, principalmente de 1 ó 2 años de edad, marcados y liberados en el Pacífico oriental han sido re capturados en el Pacífico occidental y oriental. La demora mínima entre la liberación en el Pacífico oriental y la recaptura en el Pacífico occidental es de 674 días. Se estima el coeficiente de mortalidad natural a partir de los datos de crecimiento y temperatura ambiental en un 0.276 anual, con límites de confianza al 90% de 0.161 y 0.471. El aleta azul del norte desova únicamente en el Pacífico occidental. Algunos de los juveniles migran al Pacífico oriental, donde permanecen entre varios meses y varios años antes de regresar al Pacífico occidental. La porción de los peces que migran al Pacífico oriental varía entre años, y ésto parece ser una causa importante de la variación anual en las capturas en el Océano Pacífico oriental. (PDF contains 94 pages.)

Instabilities of a liquid film flowing down an inclined porous plane

The problem of a film flowing down an inclined porous layer is considered. The fully developed basic flow is driven by gravitation. A careful linear instability analysis is carried out. We use Darcy's law to describe the porous layer and solve the coupling equations of the fluid and the porous medium rather than the decoupled equations of the one-sided model used in previous works. The eigenvalue problem is solved by means of a Chebyshev collocation method. We compare the instability of the two-sided model with the results of the one-sided model. The result reveals a porous mode instability which is completely neglected in previous works. For a falling film on an inclined porous plane there are three instability modes, i.e., the surface mode, the shear mode, and the porous mode. We also study the influences of the depth ratio d, the Darcy number delta, and the Beavers-Joseph coefficient alpha(BJ) on the instability of the system.

Properties of α-monoclinic selenium

The growth of bulky and platelet shaped α-monoclinic crystals is discussed. A simple method is devised for identifying and orienting them.

The density, previously in disagreement with the value calculated from x-ray studies, is carefully redetermined, and found to be in good agreement with the latter.

The relative dielectric constant is determined, an effort being made to eliminate errors inherent in previous measurements, which have not been in agreement. A two parameter model is derived which explains the anisotropy in the relative dielectric constant of orthorhombic sulfur, which is also composed of 8-atom puckered ring molecules. The model works less well for α-monoclinic selenium. The relative dielectric constant anisotropy is quite noticeable, being 6.06 along the crystal b axis, and 8.52-8.93 normal to the axis.

Thin crystal platelets of α-monoclinic selenium (less than 1µ thick) are used to extend optical transmission measurements up to 4.5eV. Previously the measurements extended up to 2.1 eV, limited by the thickness of the available crystals. The absorption edge is at 2.20 eV, with changes in slope of the absorption coefficient occurring at 2.85 eV and 3.8 eV. Measurement of transmission through solutions of selenium in CS_2 and trichlorethylene yield an absorption edge of 2.75 eV. There is evidence the selenium exists in solution partly as Se_8 rings, the building block of monoclinic selenium. Transmission is measured at low temperatures (80°K and 10°K) using the platelets. The absorption edge is at 2.38 eV and 2.39 eV, respectively, for the two temperatures. Measurements at low temperatures with polarized and unpolarized light reveal interesting absorption anisotropy near 2.65 eV.

基于U型塑料光纤和邻啡咯啉钌的溶解氧传感器

本文提出了一种基于U型塑料光纤和邻啡咯啉钌的溶解氧传感器。光纤传感头放在样品池中，测量了不同氮氧体积比情况下荧光强度和寿命。研究了邻啡咯啉钌的浓度、传感头的退火时间和U型光纤的弯曲半径对系统灵敏度的影响。我们提出了双层模型来解释实验结果不符合S-V方程的实验现象，并用双层模型分析了弯曲半径对系统灵敏度的影响。通过对溶解氧浓度和荧光寿命亚线性曲线的拟合，我们得出敏感层的厚度和S-V系数分别为0.59和0.61。

The Regulation of Sleep and Circadian Rhythms: The Role of Melatonin and Adenosine in Zebrafish

Sleep is a highly conserved behavioral state whose regulation is still unclear. In this thesis I initially briefly introduce the known sleep circuitry and regulation in vertebrates, and why zebrafish is seen as a good model to study sleep-regulation. I describe the existing two-process model of sleep regulation, which posits that the two processes C (circadian) and S (homeostatic) control timing of sleep-wake behavior. I then study the role melatonin plays in the circadian regulation of sleep using zebrafish. Firstly, we find that the absence of melatonin results in a reduction of sleep at night, establishing that endogenous melatonin is required for sleep at night. Secondly, melatonin mutants show a reduction in sleep in animals with no functional behavioral rhythms suggesting that melatonin does not require intact circadian rhythms for its effect on sleep. Thirdly, melatonin mutants do not exhibit any changes in circadian rhythms, suggesting that the circadian clock does not require melatonin for its function. Fourthly, we find that in the absence of melatonin, there is no rhythmic expression of sleep, suggesting that melatonin is the output molecule of process C. Lastly, we describe a connection between adenosine signaling (output molecules of process S), and melatonin. Following this we proceed to study the role adenosine signaling plays in sleep-wake behavior. We find that firstly, adenosine receptor A1 and A2 are involved in sleep- wake behavior in zebrafish, based on agonist/antagonist behavioral results. Secondly, we find that several brain regions such as PACAP cells in the rostral midbrain, GABAergic cells in the forebrain and hindbrain, Dopamine and serotonin cells in the caudal hypothalamus and sox2 cells lining the hindbrain ventricle are activated in response to the A1 antagonist and VMAT positive cells are activated in response to the A2A agonist, suggesting these areas are involved in adenosine signaling in zebrafish. Thirdly, we find that knocking out the zebrafish adenosine receptors has no effect on sleep architecture. Lastly, we find that while the A1 agonist phenotype requires the zfAdora1a receptor, the antagonist and the A2A agonist behavioral phenotypes are not mediated by the zfAdora1a, zfAdora1b and zfAdoraA2Aa, zfAdora2Ab receptors respectively.

A molecular dynamics study of the microscopic properties of simple dense fluids

The microscopic properties of a two-dimensional model dense fluid of Lennard-Jones disks have been studied using the so-called "molecular dynamics" method. Analyses of the computer-generated simulation data in terms of "conventional" thermodynamic and distribution functions verify the physical validity of the model and the simulation technique.

The radial distribution functions g(r) computed from the simulation data exhibit several subsidiary features rather similar to those appearing in some of the g(r) functions obtained by X-ray and thermal neutron diffraction measurements on real simple liquids. In the case of the model fluid, these "anomalous" features are thought to reflect the existence of two or more alternative configurations for local ordering.

Graphical display techniques have been used extensively to provide some intuitive insight into the various microscopic phenomena occurring in the model. For example, "snapshots" of the instantaneous system configurations for different times show that the "excess" area allotted to the fluid is collected into relatively large, irregular, and surprisingly persistent "holes". Plots of the particle trajectories over intervals of 2.0 to 6.0 x 10^-12 sec indicate that the mechanism for diffusion in the dense model fluid is "cooperative" in nature, and that extensive diffusive migration is generally restricted to groups of particles in the vicinity of a hole.

A quantitative analysis of diffusion in the model fluid shows that the cooperative mechanism is not inconsistent with the statistical predictions of existing theories of singlet, or self-diffusion in liquids. The relative diffusion of proximate particles is, however, found to be retarded by short-range dynamic correlations associated with the cooperative mechanism--a result of some importance from the standpoint of bimolecular reaction kinetics in solution.

A new, semi-empirical treatment for relative diffusion in liquids is developed, and is shown to reproduce the relative diffusion phenomena observed in the model fluid quite accurately. When incorporated into the standard Smoluchowski theory of diffusion-controlled reaction kinetics, the more exact treatment of relative diffusion is found to lower the predicted rate of reaction appreciably.

Finally, an entirely new approach to an understanding of the liquid state is suggested. Our experience in dealing with the simulation data--and especially, graphical displays of the simulation data--has led us to conclude that many of the more frustrating scientific problems involving the liquid state would be simplified considerably, were it possible to describe the microscopic structures characteristic of liquids in a concise and precise manner. To this end, we propose that the development of a formal language of partially-ordered structures be investigated.