A Finite Variable Difference Relaxation Scheme for hyperbolic-parabolic equations

Using the framework of a new relaxation system, which converts a nonlinear viscous conservation law into a system of linear convection-diffusion equations with nonlinear source terms, a finite variable difference method is developed for nonlinear hyperbolic-parabolic equations. The basic idea is to formulate a finite volume method with an optimum spatial difference, using the Locally Exact Numerical Scheme (LENS), leading to a Finite Variable Difference Method as introduced by Sakai [Katsuhiro Sakai, A new finite variable difference method with application to locally exact numerical scheme, journal of Computational Physics, 124 (1996) pp. 301-308.], for the linear convection-diffusion equations obtained by using a relaxation system. Source terms are treated with the well-balanced scheme of Jin [Shi Jin, A steady-state capturing method for hyperbolic systems with geometrical source terms, Mathematical Modeling Numerical Analysis, 35 (4) (2001) pp. 631-645]. Bench-mark test problems for scalar and vector conservation laws in one and two dimensions are solved using this new algorithm and the results demonstrate the efficiency of the scheme in capturing the flow features accurately.

Graphical simulation of Petri Nets

An interactive graphics package for modeling with Petri Nets has been implemented. It uses the VT-11 graphics terminal supported on the PDP-11/35 computer to draw, execute, analyze, edit and redraw a Petri Net. Each of the above mentioned tasks can be performed by selecting appropriate items from a menu displayed on the screen. Petri Nets with a reasonably large number of nodes can be created and analyzed using this package. The number of nodes supported may be increased by making simple changes in the program. Being interactive, the program seeks information from the user after displaying appropriate messages on the terminal. After completing the Petri Net, it may be executed step by step and the changes in the number of tokens may be observed on the screen, at each place. Some properties of Petri Nets like safety, boundedness, conservation and redundancy can be checked using this package. This package can be used very effectively for modeling asynchronous (concurrent) systems with Petri Nets and simulating the model by “graphical execution.”

Response of a turbulent boundary layer to a step change in surface heat flux

Measurements of both the velocity and the temperature field have been made in the thermal layer that grows inside a turbulent boundary layer which is subjected to a small step change in surface heat flux. Upstream of the step, the wall heat flux is zero and the velocity boundary layer is nearly self-preserving. The thermal-layer measurements are discussed in the context of a self-preserving analysis for the temperature disturbance which grows underneath a thick external turbulent boundary layer. A logarithmic mean temperature profile is established downstream of the step but the budget for the mean-square temperature fluctuations shows that, in the inner region of the thermal layer, the production and dissipation of temperature fluctuations are not quite equal at the furthest downstream measurement station. The measurements for both the mean and the fluctuating temperature field indicate that the relaxation distance for the thermal layer is quite large, of the order of 1000θ0, where θ0 is the momentum thickness of the boundary layer at the step. Statistics of the thermal-layer interface and conditionally sampled measurements with respect to this interface are presented. Measurements of the temperature intermittency factor indicate that the interface is normally distributed with respect to its mean position. Near the step, the passive heat contaminant acts as an effective marker of the organized turbulence structure that has been observed in the wall region of a boundary layer. Accordingly, conditional averages of Reynolds stresses and heat fluxes measured in the heated part of the flow are considerably larger than the conventional averages when the temperature intermittency factor is small.

CFC Alternatives: A Fresh Look

This article examines, through a molecular perspective, the 'ozone-friendly' refrigerants R-134a and R-123 vis-à-vis R-12 and R-11, which are targeted to be phased out under the Montreal Protocol on Substances that Deplete the Ozone Layer, Final Act (1987). It appears that the molecular weight, size parameter, and dipole moment, of R-134a and R-123, may induce a pronounced effect on the chemical equilibrium of ice particles in the polar stratospheric clouds and subsequently influence the photochemical reactions therein. Non-polar, high-molecular-weight perfluoropropane (R-218), could be a better substitute for R-12, while R-134, which is a non-polar HFC of the ethane family, could also be a candidate although its molecular weight is lower than that of R-12. The search for a good substitute for R-11, however, must continue.

The fully implicit stochastic-alpha method for stiff stochastic differential equations

A fully implicit integration method for stochastic differential equations with significant multiplicative noise and stiffness in both the drift and diffusion coefficients has been constructed, analyzed and illustrated with numerical examples in this work. The method has strong order 1.0 consistency and has user-selectable parameters that allow the user to expand the stability region of the method to cover almost the entire drift-diffusion stability plane. The large stability region enables the method to take computationally efficient time steps. A system of chemical Langevin equations simulated with the method illustrates its computational efficiency.

Pi-turns in proteins and peptides: Classification, conformation, occurrence, hydration and sequence.

The i + 5-->i hydrogen bonded turn conformation (pi-turn) with the fifth residue adopting alpha L conformation is frequently found at the C-terminus of helices in proteins and hence is speculated to be a "helix termination signal." An analysis of the occurrence of i + 5-->i hydrogen bonded turn conformation at any general position in proteins (not specifically at the helix C-terminus), using coordinates of 228 protein crystal structures determined by X-ray crystallography to better than 2.5 A resolution is reported in this paper. Of 486 detected pi-turn conformations, 367 have the (i + 4)th residue in alpha L conformation, generally occurring at the C-terminus of alpha-helices, consistent with previous observations. However, a significant number (111) of pi-turn conformations occur with (i + 4)th residue in alpha R conformation also, generally occurring in alpha-helices as distortions either at the terminii or at the middle, a novel finding. These two sets of pi-turn conformations are referred to by the names pi alpha L and pi alpha R-turns, respectively, depending upon whether the (i + 4)th residue adopts alpha L or alpha R conformations. Four pi-turns, named pi alpha L'-turns, were noticed to be mirror images of pi alpha L-turns, and four more pi-turns, which have the (i + 4)th residue in beta conformation and denoted as pi beta-turns, occur as a part of hairpin bend connecting twisted beta-strands. Consecutive pi-turns occur, but only with pi alpha R-turns. The preference for amino acid residues is different in pi alpha L and pi alpha R-turns. However, both show a preference for Pro after the C-termini. Hydrophilic residues are preferred at positions i + 1, i + 2, and i + 3 of pi alpha L-turns, whereas positions i and i + 5 prefer hydrophobic residues. Residue i + 4 in pi alpha L-turns is mainly Gly and less often Asn. Although pi alpha R-turns generally occur as distortions in helices, their amino acid preference is different from that of helices. Poor helix formers, such as His, Tyr, and Asn, also were found to be preferred for pi alpha R-turns, whereas good helix former Ala is not preferred. pi-Turns in peptides provide a picture of the pi-turn at atomic resolution. Only nine peptide-based pi-turns are reported so far, and all of them belong to pi alpha L-turn type with an achiral residue in position i + 4. The results are of importance for structure prediction, modeling, and de novo design of proteins.

Three-dimensional hybrid networks based on aspartic acid

Three-dimensional achiral coordination polymers of the general formula M2(D, l-NHCH (COO)CH2COO)2·C4H4N2 where M = Ni and Co and pyrazine acts as the linker molecule have been prepared under hydrothermal conditions starting with [M(L-NHCH(COO)CH2COO)·3H2O] possessing a helical chain structure. A three-dimensional hybrid compound of the formula Pb2.5[N{CH(COO) CH2COO}22H2O] has also been prepared hydrothermally starting with aspartic acid and Pb(NO3)2. In this lead compound, where a secondary amine formed by the dimerisation of aspartic acid acts as the ligand, there is two-dimensional inorganic connectivity and one-dimensional organic connectivity.

Negative nuclear Overhauser effects as probes of macromolecular structure

Nuclear Overhauser effects (NOE) and circular dichroism (CD) techniques have been used to probe @-turn conformations in acyclic and cyclic peptides containingPro-Xsequences. The model peptides studied are of the type Piv-Pro-X-NHMe (X = Aib, D-Ala, Gly, Val, and Leu) and Boc-Cys-Pro-X-C s NHMe (X = Aib, L-Ala, D-Ala, Gly, and Leu). In the acyclic series, observation of NOES between Pro C"H and X-NH, together with solvent and temperature dependence of NH chemical shifts, establishes a 4 - 1 hydrogen bond stabilized type I1 @-turn in the Gly, D-Ala, and Aib peptides, in CDC13 and (CD3)2S0. A positive n-r* CD band at -225-230 nm appears to be characteristic of this structure. For the acyclic Pro-Leu peptide the observation of NOE's for both Pro and Leu C"H resonances on saturation of Leu NH is compatible with a type V bend or consecutive y-turn conformation. In the cyclic disulfide series the Pro-Aib and Pro-D-Ala peptides favor type I1 @-turns, whereas all other peptides adopt type I (111) conformations. All the cyclic disulfides exhibit an intense negative CD band at -228-230 nm. The results suggest thatgeneralcorrelations between CD spectral type and specific 0-turn conformations may not be obtained. Evidence for solvent-dependent structural changes in the Pro-Aib sequence in both cyclic and acyclic peptides is presented.

Transduction in mycobacterium smegmatis

SEVERAL people have studied lysogeny and lysogenic conversion in the advanced prokaryote, Mycobacterium1-3, and presented evidence for reciprocal genetic exchange from the lysogenic to the non-lysogenic state. Biochemical and physiological changes induced by lysogenization of M. phlei with mycobacteriophage B2 have been demonstrated by Juhasz4. We wish to report the isolation of a mycobacteriophage capable of mediating transduction in Mycobacterium smegmatis strain SN2, thus providing an opportunity to study the genetics of mycobacteria.

A class of one-dimensional steady-state flows in radiation-gas-dynamics

The study of steady-state flows in radiation-gas-dynamics, when radiation pressure is negligible in comparison with gas pressure, can be reduced to the study of a single first-order ordinary differential equation in particle velocity and radiation pressure. The class of steady flows, determined by the fact that the velocities in two uniform states are real, i.e. the Rankine-Hugoniot points are real, has been discussed in detail in a previous paper by one of us, when the Mach number M of the flow in one of the uniform states (at x=+∞) is greater than one and the flow direction is in the negative direction of the x-axis. In this paper we have discussed the case when M is less than or equal to one and the flow direction is still in the negative direction of the x-axis. We have drawn the various phase planes and the integral curves in each phase plane give various steady flows. We have also discussed the appearance of discontinuities in these flows.

Evidence for two distinct orthorhombic structures associated with different Tc regimes in LnBa2Cu3O7 − δ (Ln = Nd, Eu, Gd and Dy): A study of the dependence of superconductivity on oxygen stoichiometry

Superconductivity in LnBa2Cu3O7 − δ with Ln = Nd, Eu, Gdand Dy has been investigated as a function of δ, closely following the accompanying changes in crystal structure. Orthorhombic GdBa2Cu3O7 − δ and DyBa2Cu3O7 − δ show a Tc of ≈ 90 K up to δ = 0.2 and a lower Tc plateau (40–50 K) in the δ range 02 to 0.4, similar to that found in YBa2Cu3O7 − δ. The orthorhombic structure II in the lower Tc regions is different from the structure I in the 90 K Tc (low δ) region. The unit cell parameters of the orthorhombic I structure in the high Tc region bear the relationship of a a ≠ b not, vert, similar c/3. This relationship is not seen in the low Tc plateau. The low Tc plateau region does not distinctly manifest itself in NdBa2Cu3O7 − δ just as in LaBa2Cu3O7 − δ.

A series of metallic oxides of the formula La3LnBaCu5O13+δ (Ln = rare earth or Y)

Oxides of the formula La3LnBaCu5O13+δ (Ln = Nd, Sm, Gd, Dy, or Y) exhibiting metallic resistivity have been prepared and characterized. In the case of yttrium, a composition close to La2Y2BaCu5O13+δ, which is also metallic, could be prepared.

A comparative study of superconducting LaBa2Cu3O7−δ and La3−xBa3+xCu6O14+δ

Based on X-ray diffraction and electron microscopy it is shown that oxides of the general formula La Ba2Cu3O7−δ become tetragonal when δ deviates slightly from 0. This tetragonal structure is similar to that of La3−xBa3+xCu6O14+δ, with a cubic perovskite subcell and triple periodicity. Electron micrographs of these tetragonal oxides show 90° microdomains. Orthorhombic LaBa2Cu3O7−δ with high Tc (not, vert, similar77 K) is found only when δ reverse similar, equals 0; this sample is subject to formation of twins. Fluorine substitution seems to favor superconductivity.