Computation of confined swirling flow using RNG and k-epsilon turbulence models

In this work we numerically model isothermal turbulent swirling flow in a cylindrical burner. Three versions of the RNG k-epsilon model are assessed against performance of the standard k-epsilon model. Sensitivity of numerical predictions to grid refinement, differing convective differencing schemes and choice of (unknown) inlet dissipation rate, were closely scrutinised to ensure accuracy. Particular attention is paid to modelling the inlet conditions to within the range of uncertainty of the experimental data, as model predictions proved to be significantly sensitive to relatively small changes in upstream flow conditions. We also examine the characteristics of the swirl--induced recirculation zone predicted by the models over an extended range of inlet conditions. Our main findings are: - (i) the standard k-epsilon model performed best compared with experiment; - (ii) no one inlet specification can simultaneously optimize the performance of the models considered; - (iii) the RNG models predict both single-cell and double-cell IRZ characteristics, the latter both with and without additional internal stagnation points. The first finding indicates that the examined RNG modifications to the standard k-e model do not result in an improved eddy viscosity based model for the prediction of swirl flows. The second finding suggests that tuning established models for optimal performance in swirl flows a priori is not straightforward. The third finding indicates that the RNG based models exhibit a greater variety of structural behaviour, despite being of the same level of complexity as the standard k-e model. The plausibility of the predicted IRZ features are discussed in terms of known vortex breakdown phenomena.

Computation of turbulent flow in an IFRF quarl burner

A computational model for isothermal axisymmetric turbulent flow in a quarl burner is set up using the CFD package FLUENT, and numerical solutions obtained from the model are compared with available experimental data. A standard k-e model and and two versions of the RNG k-e model are used to model the turbulence. One of the aims of the computational study is to investigate whether the RNG based k-e turbulence models are capable of yielding improved flow predictions compared with the standard k-e turbulence model. A difficulty is that the flow considered here features a confined vortex breakdown which can be highly sensitive to flow behaviour both upstream and downstream of the breakdown zone. Nevertheless, the relatively simple confining geometry allows us to undertake a systematic study so that both grid-independent and domain-independent results can be reported. The systematic study includes a detailed investigation of the effects of upstream and downstream conditions on the predictions, in addition to grid refinement and other tests to ensure that numerical error is not significant. Another important aim is to determine to what extent the turbulence model predictions can provide us with new insights into the physics of confined vortex breakdown flows. To this end, the computations are discussed in detail with reference to known vortex breakdown phenomena and existing theories. A major conclusion is that one of the RNG k-e models investigated here is able to correctly capture the complex forward flow region inside the recirculating breakdown zone. This apparently pathological result is in stark contrast to the findings of previous studies, most of which have concluded that either algebraic or differential Reynolds stress modelling is needed to correctly predict the observed flow features. Arguments are given as to why an isotropic eddy-viscosity turbulence model may well be able to capture the complex flow structure within the recirculating zone for this flow setup. With regard to the flow physics, a major finding is that the results obtained here are more consistent with the view that confined vortex breakdown is a type of axisymmetric boundary layer separation, rather than a manifestation of a subcritical flow state.

A Case Study of the Conformation of Poly(alpha-aminoisobutyric acid): alpha- or 310-Helix

The relative stabilities of a- and Blo-helical structures for polymers of a-aminoisobutyric acid (Aib) have been worked out, using the classical potential energy functions. To make a comparative study, we have used Buckingham "6-exp" and Kitaigorodsky's potential functions. Conformational analysis of the dipeptide segment with Aib residue indicates the necessity for nonplanar distortion of the peptide unit, which is a common feature in the observed crystal structures with Aib residues. In the range of Aw -10 to +loo studied, a-helical conformations are preferred in the region -3" < Aw < +loo, and Blo-helical conformations are preferred in the region -3" > Aw > -10'. Minimum energy conformations for right-handed structures are found in the +ue region of Aw and correspondingly for left-handed structures in the -ue region of Aw. For Aw - 6", a-helical structures have four- or near fourfold symmetry with h - 1.5 A. Such a helix with n = 4 and h = 1.5 A is termed an a'-helix. This structure is found to be consistent with the electron diffraction data of Malcolm3 and energetically more favorable than the standard 310-helix.

Conformational analysis of cyclo (l-cystine

Conformational analysis of cyclo(L-cystine) shows that the diketopiperazine ring has to exist only in the boat form. With this geometry, the molecule can adopt two distinct forms differing mainly in the chirality of the disulphide bridge. In both the P- and M-models, corresponding to dihedral angles of nearly + 90° and —90° respectively about the S-S bond, the molecule displays an approximate two-fold symmetry. According to our semi-empirical energy calculations, the minimum energy of the M-model is —9.2 kcal/mol, only 0.3 kcal/mol lower than that of the P-model. Because the difference between the two minima is so small, neither form is clearly superior to the other. However, the number of low energy conformations of the M-model in the allowed conformational space is significantly larger than that of the P-model by a ratio of 3 to 1, and therefore the former is likely to be thermodynamically favoured.

Zeeman effect of NQR in 4,6-dichloropyrimidine and 6-chloro 2,4-dimethoxypyrimidine

Zeeman (35Cl) NQR studies in polycrystalline samples of 4,6-dichloropyrimidine and 6 chloro 2,4 dimethoxypyrimidine show that the asymmetry at the four chemically inequivalent chlorine sites in the former is about 10%, while in the latter (one line) the asymmetry is almost zero. Using a valence-bond picture, C-Cl bonds in 4,6-dichloropyrimidine have been characterised, and the results are also compared with those in a corresponding benzene compound using a simple molecular orbital calculation. The axial symmetry of C-Cl bond in the second compound has been attributed to mesomeric effects.

Significance of Branching Angles in Streamer Discharges

The phenomenon of branching at specific angles in streamer breakdown studies is found to be more universal than was ever thought. The angles measured in the breakdown of gases show that the coefficient of field distortion, K, lies in the range from 1 to less than 0.1. The values of K, so obtained agree well with those envisaged in the criterion of the streamer mechanism. It is hoped that branching angles observed in various types of breakdown may, possibly, be explained similarly.

ESR study of copper diethyldithiophosphate

ESR investigations are reported in single crystals of copper diethyldithiophosphate, magnetically diluted with the corresponding diamagnetic nickel complex. The spectrum at normal gain shows hyperfine components from 63Cu, 65Cu, and 31P nuclei. At much higher gain, hyperfine interaction from 33S nuclei in the ligand is detected. The spin Hamiltonian parameters relating to copper show tetragonal symmetry. The measured parameters are g|| = 2.085, g[perpendicular]=2.025, A63Cu = 149.6 × 10−4 cm−1, A65Cu = 160.8 × 10−4 cm−1, BCu = 32.5 × 10−4 cm−1 and QCu [infinity] 5.5 × 10−4cm−1. The 31P interaction is isotropic with a coupling constant AP = 9.6 × 10−4 cm−1. Angular variation of the 33S lines shows two different hyperfine tensors indicating the presence of two chemically inequivalent Cu[Single Bond]S bonds. The experimentally determined hyperfine constants are A ₁^s=34.9×10−4 cm−1, B ₁^s=26.1×10−4 cm−1, A ₂^s=60.4×10−4 cm−1, B₂^s=55.5×10−4 cm−1. The hyperfine parameters show that the hybridization of the ligand orbitals is very sensitive to the symmetry around the ligand. The g values and Cu hyperfine parameters are not much affected by the distortions occurring in the ligand. The energies of the d-d transitions are determined by optical absorption measurements on Cu diethyldithiophosphate in solution. Using the spin Hamiltonian parameters together with optical absorption results, the MO parameters for the complex are calculated. It is found that in addition to the sigma bond, the pi bonds are also strongly covalent. ©1973 The American Institute of Physics.

Social evolution - Has nature ever rewound the tape?

Social insects such as ants, bees, wasps and termites exhibit extreme forms of altruism where some individuals remain sterile and assist other individuals in reproduction. Hamilton's inclusive fitness theory provides a powerful framework for investigating the evolution of such altruism. Using the paper wasp Ropalidia marginata, we have quantified and delineated the role of ecological, physiological, genetic and demographic factors in social evolution. An interesting feature of the models we have developed is their symmetry so that either altruism or selfishness can evolve, depending on the numerical values of various parameters. This suggests that selfish/solitary behaviour must occasionally re-emerge even from the eusocial state, It is useful to contemplate expected intermediate states during such potential reversals. We can perhaps envisage three successive steps in such a hypothetical process: i) workers revolt against the hegemony of the queen and challenge her status as the sole reproductive, ii) workers stop producing queens and one or more of them function as egg layers (functional queen/s) capable of producing both haploid as well as diploid offspring and iii) social evolution reverses completely so that a eusocial species becomes solitary, at least facultatively. It appears that the third step, namely transition from eusociality to the solitary state, is rare and has been restricted to transitions from the primitively eusocial state only. The absence of transitions from the highly eusocial state to the solitary state may be attributed to a number of 'preventing mechanisms' such as (a) queen control of workers (b) loss of spermathecae and ability to mate (c) morphological specialization (d) caste polyethism and (e) homeostasis, which must each make the transition difficult and, taken together, perhaps very difficult. However, the discovery of a transition from the highly eusocial to the solitary state can hardly he ruled out, given that little or no effort has gone into its detection. In this paper I discuss social evolution and its possible reversal and cite potential examples of stages in the transition from the social to the solitary.

Thin nanoporous metal-insulator-metal membranes

Insulating nanoporous materials are promising platforms for soft-ionizing membranes; however, improvement in fabrication processes and the quality and high breakdown resistance of the thin insulator layers are needed for high integration and performance. Here, scalable fabrication of highly porous, thin, silicon dioxide membranes with controlled thickness is demonstrated using plasma-enhanced chemical-vapor-deposition. The fabricated membranes exhibit good insulating properties with a breakdown voltage of 1 × 107 V/cm. Our calculations suggest that the average electric field inside a nanopore of the membranes can be as high as 1 × 106 V/cm; sufficient for ionization of wide range of molecules. These metal–insulator–metal nanoporous arrays are promising for applications such soft ionizing membranes for mass spectroscopy.

The electronic factor in Bergman cyclization

Through-bond interactions in 1,4-dehydrobenzene preferentially stabilize the out-of-phase combination of the radical hydrids, The resultant splitting between the frontier orbitals is crucial in making Bergman cyclization a symmetry-allowed process. Orbital symmetry also inhibits the radical centers from forming a C-C bond, enabling the biradical to survive as a local minimum capable of intermolecular hydrogen abstraction, Both these factors, which are important in the design of DNA cleaving molecules, are confirmed through calculations on biradicals formed from diynes in which through-bond interactions stabilize the in-phase combination of hybrids at the radical centers.

Structure of bis(thiocyanato-N)bis(thiosemicarbazide-N1,S)nickel(II): a redetermination

[Ni(NCS)2(CHsN3S)2], Mr = 356.7, monoclinic, P21/c , a = 5-297 (1), b = 7.869 (1), c - 16-078 (2) A,/3 = 91.53 (1) °, V-= 669.9 A 3, Z= 2, Om = 1"76, Dx = 1"771 g cm -3, A(Mo Ka) = 0-71069 ]k, /.~ = 19"9 cm-l, F(000) = 364, T = 295 K, final R = 0.026 for 1576 significant [F > 10g(F)] reflections. The complex lies on a crystallographic centre of symmetry. The Ni atom is octahedrally coordinated by two thiocyanates (through N atoms) and by two thiosemicarbazide molecules (through hydrazinic N and S atoms). The crystal structure is stabilized by N--H...S hydrogen bonds. Early work on this structure [Garaj & Dunaj-Jurco (1968). Chem. Commun. p. 518] used photographic data and was refined to R = 0-13 for 512 reflections.

Electron cryo-microscopy studies of bacteriophage phi8 and archaeal virus SH1

Symmetry is a key principle in viral structures, especially the protein capsid shells. However, symmetry mismatches are very common, and often correlate with dynamic functionality of biological significance. The three-dimensional structures of two isometric viruses, bacteriophage phi8 and the archaeal virus SH1 were reconstructed using electron cryo-microscopy. Two image reconstruction methods were used: the classical icosahedral method yielded high resolution models for the symmetrical parts of the structures, and a novel asymmetric in-situ reconstruction method allowed us to resolve the symmetry mismatches at the vertices of the viruses. Evidence was found that the hexameric packaging enzyme at the vertices of phi8 does not rotate relative to the capsid. The large two-fold symmetric spikes of SH1 were found not to be responsible for infectivity. Both virus structures provided insight into the evolution of viruses. Comparison of the phi8 polymerase complex capsid with those of phi6 and other dsRNA viruses suggests that the quaternary structure in dsRNA bacteriophages differs from other dsRNA viruses. SH1 is unusual because there are two major types of capsomers building up the capsid, both of which seem to be composed mainly of single beta-barrels perpendicular to the capsid surface. This indicates that the beta-barrel may be ancestral to the double beta-barrel fold.

The hydrogen bond: a molecular beam microwave spectroscopist's view with a universal appeal

In this manuscript, we propose a criterion for a weakly bound complex formed in a supersonic beam to be characterized as a `hydrogen bonded complex'. For a `hydrogen bonded complex', the zero point energy along any large amplitude vibrational coordinate that destroys the orientational preference for the hydrogen bond should be significantly below the barrier along that coordinate so that there is at least one bound level. These are vibrational modes that do not lead to the breakdown of the complex as a whole. If the zero point level is higher than the barrier, the `hydrogen bond' would not be able to stabilize the orientation which favors it and it is no longer sensible to characterize a complex as hydrogen bonded. Four complexes, Ar-2-H2O, Ar-2-H2S, C2H4-H2O and C2H4-H2S, were chosen for investigations. Zero point energies and barriers for large amplitude motions were calculated at a reasonable level of calculation, MP2(full)/aug-cc-pVTZ, for all these complexes. Atoms in molecules (AIM) theoretical analyses of these complexes were carried out as well. All these complexes would be considered hydrogen bonded according to the AIM theoretical criteria suggested by Koch and Popelier for C-H center dot center dot center dot O hydrogen bonds (U. Koch and P. L. A. Popelier, J. Phys. Chem., 1995, 99, 9747), which has been widely and, at times, incorrectly used for all types of contacts involving H. It is shown that, according to the criterion proposed here, the Ar-2-H2O/H2S complexes are not hydrogen bonded even at zero kelvin and C2H4-H2O/H2S complexes are. This analysis can naturally be extended to all temperatures. It can explain the recent experimental observations on crystal structures of H2S at various conditions and the crossed beam scattering studies on rare gases with H2O and H2S.

Metal-Nucleotide Interactions: Crystal Structures of Alkali (Li+, Na+, K+) and Alkaline Earth (Ca2+, Mg2+) Metal Complexes of Adenosine 2'-Monophosphate

Crystal structures of lithium, sodium, potassium, calcium and magnesium salts of adenosine 2'-monophosphate (2'-AMP) have been obtained at atomic resolution by X-ray crystallographic methods. 2'-AMP.Li belongs to the monoclinic space group P21 with a = 7.472(3)Å, b = 26.853(6) Å, c = 9.184(1)Å, b = 113.36(1)Å and Z= 4. 2'-AMP.Na and 2'-AMP.K crystallize in the trigonal space groups P31 and P3121 with a = 8.762(1)Å, c = 34.630(5)Å, Z= 6 and a = 8.931(4), Åc = 34.852(9)Å and Z= 6 respectively while 2'-AMP.Ca and 2'-AMP.Mg belong to space groups P6522 and P21 with cell parameters a = 9.487(2), c = 74.622(13), Z = 12 and a = 4.973(1), b = 10.023(2), c = 16.506(2), beta = 91.1(0) and Z = 2 respectively. All the structures were solved by direct methods and refined by full matrix least-squares to final R factors of 0.033, 0.028, 0.075, 0.069 and 0.030 for 2'-AMP.Li, 2'-AMP.Na, 2'- AMP.K, 2'-AMP.Ca and 2'-AMP.Mg, respectively. The neutral adenine bases in all the structures are in syn conformation stabilized by the O5'-N3 intramolecular hydrogen bond as in free acid and ammonium complex reported earlier. In striking contrast, the adenine base is in the anti geometry (cCN = -156.4(2)°) in 2'-AMP.Mg. Ribose moieties adopt C2'-endo puckering in 2'-AMP.Li and 2'-AMP.Ca, C2'-endo-C3'-exo twist puckering in 2'-AMP.Na and 2'-AMP.K and a C3'-endo-C2'-exo twist puckering in 2'-AMP.Mg structure. The conformation about the exocyclic C4'-C5' bond is the commonly observed gauche-gauche (g+) in all the structures except the gauche- trans (g-) conformation observed in 2'-AMP.Mg structure. Lithium ions coordinate with water, ribose and phosphate oxygens at distances 1.88 to 1.99Å. Na+ ions and K+ ions interact with phosphate and ribose oxygens directly and with N7 indirectly through a water oxygen. A distinct feature of 2'-AMP.Na and 2'-AMP.K structures is the involvement of ribose O4' in metal coordination. The calcium ion situated on a two-fold axis coordinates directly with three oxygens OW1, OW2 and O2 and their symmetry mates at distances 2.18 to 2.42Å forming an octahedron. A classic example of an exception to the existence of the O5'-N3 intramolecular hydorgen bond is the 2'-AMP.Mg strucure. Magnesium ion forms an octahedral coordination with three water and three phosphate oxygens at distances ranging from 2.02 to 2.11Å. A noteworthy feature of its coordination is the indirect link with N3 through OW3 oxygen resulting in macrochelation between the base and the phosphate group. Greater affnity of metal clays towards 5' compared to 2' and 3' nucleotides (J. Lawless, E. Edelson, and L. Manring, Am. Chem. Soc. Northwest Region Meeting, Seattle. 1978) due to macrochelation infered from solution studies (S. S. Massoud, H. Sigel, Eur. J. Biochem. 179, 451-458 (1989)) and interligand hydrogen bonding induced by metals postulated from metal-nucleotide structures in solid state (V. Swaminathan and M. Sundaralingam, CRC. Crit. Rev. Biochem. 6, 245-336 (1979)) are borne out by our structures also. The stacking patterns of adenine bases of both 2'-AMP.Na and 2'-AMP.K structures resemble the 2'-AMP.NH4 structure reported in the previous article. 2'-AMP.Li, 2'-AMP.Ca and 2'-AMP.Mg structures display base-ribose O4' stacking. An overview of interaction of monovalent and divalent cations with 2' and 5'-nucleotides has been presented.

Theoretical calculations of base-base interactions in nucleic acids: I Stacking interactions in free bases

Stacking interactions in free bases were computed on the basis of molecular association. The results of the calculations were compared with the stacking patterns observed in a few single crystals of nucleic acid components as examples. The following are the conclusions: (i) there can be two types of stacking pattern classified as normal and inverted types for any two interacting bases and both can be energetically favourable (ii) in both the types the stacking interaction is a combined effect of the overlap of the interacting bases and relative positions and orientations of the atomic centres of the two bases (iii) crystal symmetry and H-bonding interaction may influence stacking patterns.