Composite scheme using localized relaxation with non-standard finite difference method for hyperbolic conservation laws

Non-standard finite difference methods (NSFDM) introduced by Mickens [Non-standard Finite Difference Models of Differential Equations, World Scientific, Singapore, 1994] are interesting alternatives to the traditional finite difference and finite volume methods. When applied to linear hyperbolic conservation laws, these methods reproduce exact solutions. In this paper, the NSFDM is first extended to hyperbolic systems of conservation laws, by a novel utilization of the decoupled equations using characteristic variables. In the second part of this paper, the NSFDM is studied for its efficacy in application to nonlinear scalar hyperbolic conservation laws. The original NSFDMs introduced by Mickens (1994) were not in conservation form, which is an important feature in capturing discontinuities at the right locations. Mickens [Construction and analysis of a non-standard finite difference scheme for the Burgers–Fisher equations, Journal of Sound and Vibration 257 (4) (2002) 791–797] recently introduced a NSFDM in conservative form. This method captures the shock waves exactly, without any numerical dissipation. In this paper, this algorithm is tested for the case of expansion waves with sonic points and is found to generate unphysical expansion shocks. As a remedy to this defect, we use the strategy of composite schemes [R. Liska, B. Wendroff, Composite schemes for conservation laws, SIAM Journal of Numerical Analysis 35 (6) (1998) 2250–2271] in which the accurate NSFDM is used as the basic scheme and localized relaxation NSFDM is used as the supporting scheme which acts like a filter. Relaxation schemes introduced by Jin and Xin [The relaxation schemes for systems of conservation laws in arbitrary space dimensions, Communications in Pure and Applied Mathematics 48 (1995) 235–276] are based on relaxation systems which replace the nonlinear hyperbolic conservation laws by a semi-linear system with a stiff relaxation term. The relaxation parameter (λ) is chosen locally on the three point stencil of grid which makes the proposed method more efficient. This composite scheme overcomes the problem of unphysical expansion shocks and captures the shock waves with an accuracy better than the upwind relaxation scheme, as demonstrated by the test cases, together with comparisons with popular numerical methods like Roe scheme and ENO schemes.

Structure of polyamidoamide dendrimers up to limiting generations: A mesoscale description

The polyamidoamide (PAMAM) class of dendrimers was one of the first dendrimers synthesized by Tomalia and co-workers at Dow. Since its discovery the PAMAMs have stimulated many discussions on the structure and dynamics of such hyperbranched polymers. Many questions remain open because the huge conformation disorder combined with very similar local symmetries have made it difficult to characterize experimentally at the atomistic level the structure and dynamics of PAMAM dendrimers. The higher generation dendrimers have also been difficult to characterize computationally because of the large size (294852 atoms for generation 11) and the huge number of conformations. To help provide a practical means of atomistic computational studies, we have developed an atomistically informed coarse-grained description for the PAMAM dendrimer. We find that a two-bead per monomer representation retains the accuracy of atomistic simulations for predicting size and conformational complexity, while reducing the degrees of freedom by tenfold. This mesoscale description has allowed us to study the structural properties of PAMAM dendrimer up to generation 11 for time scale of up to several nanoseconds. The gross properties such as the radius of gyration compare very well with those from full atomistic simulation and with available small angle x-ray experiment and small angle neutron scattering data. The radial monomer density shows very similar behavior with those obtained from the fully atomistic simulation. Our approach to deriving the coarse-grain model is general and straightforward to apply to other classes of dendrimers.

Insight into the early stages of thermal unfolding of peanut agglutinin by molecular dynamics simulations

Peanut agglutinin is a homotetrameric nonglycosylated protein. The protein has a unique open quaternary structure. Molecular dynamics simulations have been employed follow the atomistic details of its unfolding at different temperatures. The early events of the deoligomerization of the protein have been elucidated in the present study. Simulation trajectories of the monomer as well as those of the tetramer have been compared and the tetramer is found to be substantially more stable than its monomeric counterpart. The tetramer shows retention of most of its.. secondary structure but considerable loss of the tertiary structure at high temperature. e generation of a This observation impies the molten globule-like intermediate in the later stages of deoligomerization. The quaternary structure of the protein has weakened to a large extent, but none of the subunits are separated. In addition, the importance of the metal-binding to the stability of the protein structure has also been investigated. Binding of the metal ions not only enhances the local stability of the metal-ion binding loop, but also imparts a global stability to the overall structure. The dynamics of different interfaces vary significantly as probed through interface clusters. The differences are substantially enhanced at higher temperatures. The dynamics and the stability of the interfaces have been captured mainly by cluster analysis, which has provided detailed information on the thermal deoligomerization of the protein.

Quality assessment of CdZnTe (Zn similar to 4 %) crystals

We have studied the as grown and annealed CdZnTe (Zn similar to 4 %) crystals for the assessment of their crystalline quality. As grown crystals suffer from tellurium precipitates and cadmium vacancies, which are inherent, due to retrograde solid solubility curve in the phase diagram. This is reflected in the Fourier transform infrared (FTIR) spectra over the 400 - 4500 cm(-1) range by a strong absorption around 2661 cm(-1) which corresponds to the band gap of tellurium confirming their presence, where-as a monotonic decrease in the transmission with the decrease in wave number indicates the presence of cadmium vacancies. Obviously the presence of Cd vacancies lead to the formation of tellurium precipitates confirming their presence. Annealed samples under cadmium + zinc ambient at 650 degrees C for 6 hours show an improvement in the transmission over the same range. This can be attributed to thermo-migration of tellurium precipitates and hence bonding with Cd or Zn to form CdZnTe. This is further supported by the reduced full width at half maximum in the X-ray diffraction rocking curve of these CdZnTe crystals. Cadmium annealing although can passivate Cd vacancy related defects and reduce the Te precipitates, as is observed in our low temperature Photoluminescence (PL) spectra, alone may not be sufficient possibly due to the loss of Zn. Vacuum annealing at 650 degrees C for 6 hours further deteriorated the material quality as is reflected in the low temperature PL spectra by the introduction of a new defect band around 0.85 eV and reduced IR transmission.

The effects of manganese content and mould size on abrasion and slurry erosion behaviour of chromium-manganese iron systems investigated by positron lifetime spectroscopy

Abrasion and slurry erosion behaviour of chromium-manganese iron samples with chromium (Cr) in the range similar to 16-19% and manganese (Mn) at 5 and 10% levels have been characterized for hardness followed by microstructural examination using optical and scanning electron microscopy. Positron lifetime studies have been conducted to understand the defects/microporosity influence on the microstructure. The samples were heat treated and characterized to understand the structural transformations in the matrix. The data reveals that hardness decreased with increase in Mn content from 5 to 10% in the first instance and then increase in the section size in the other case, irrespective of the sample conditions. The abrasion and slurry erosion losses show increase with increase in the section size as well as with increase in Mn content. The positron results show that as hardness increases from as-cast to heat treated sample, the positron trapping rate and hence defect concentration showed opposite trend as expected. So a good correlation between defects concentration and the hardness has been observed. These findings also corroborate well with the microstructural features obtained from optical and scanning electron microscopy. (C) 2009 Elsevier B. V. All rights reserved.

Enhanced Tetrahedral Ordering of Water Molecules in Minor Grooves of DNA: Relative Role of DNA Rigidity, Nanoconfinement, and Surface Specific Interactions

Confinement and Surface specific interactions call induce Structures otherwise unstable at that temperature and pressure. Here we Study the groove specific water dynamics ill the nucleic acid sequences, poly-AT and poly-GC, in long B-DNA duplex chains by large scale atomistic molecular dynamics simulations, accompanied by thermodynamic analysis. While water dynamics in the major groove remains insensitive to the sequence differences, exactly the opposite is true for the minor groove water. Much slower water dynamics observed in the minor grooves (especially in the AT minor) call be attributed to all enhanced tetrahedral ordering (< t(h)>) of water. The largest value of < t(h)> in the AT minor groove is related to the spine of hydration found in X-ray Structure. The calculated configurational entropy (S-C) of the water molecules is found to be correlated with the self-diffusion coefficient of water in different region via Adam-Gibbs relation D = A exp(-B/TSC), and also with < t(h)>.

Melting of an Anchored Bilayer: Molecular Dynamics Simulations of the Structural Transition in (CnH2n+1NH3)(2)PbI4 (n=12, 14, 16, 18)

The thermally driven Structural phase transition in the organic-inorganic hybrid perovskite (CnH2n+1NH3)(2)PbI4 has been investigated using molecular dynamics (MD) simulations. This system consists of positively charged alkyl-amine chains anchored to a rigid negatively charged PbI4 sheet with the chains organized as bilayers with a herringbone arrangement. Atomistic simulations were performed using ail isothermal-isobaric ensemble over a wide temperature range from 65 to 665 K for different alkyl chain lengths, n = 12, 14, 16, and 18. The simulations are able to reproduce the essential Features of the experimental observations of this system, including the existence of a transition, the linear variation of the transition temperature with alkyl chain length, and the expansion of the bilayer thickness at the transition. By use of the distance fluctuation Criteria, it is Shown that the transition is associated With a Melting of the alkyl chains of the anchored bilayer. Ail analysis of the conformation of the alkyl chains shows increased disorder in the form of gauche defects above due melting transition. Simulations also show that the melting transition is characterized by the complete disappearance of all-trans alkyl chains in the anchored bilayer, in agreement with experimental observations. A conformationally disordered chain has a larger effective cross-sectional area, and above due transition a uniformly tilted arrangement of the anchored chains call no longer be Sustained. At the melt the angular distribution of the orientation of the chains are 110 longer uniform; the chains are splayed allowing for increased space for individual chains of the anchored bilayer. This is reflected in a sharp rise in the ratio of the mean head-to-head to tail-to-tail distance of the chains of the bilayer at the transition resulting in in expansion of the bilayer thickness. The present MD simulations provide a simple explanation as to how changes in conformation of individual alkyl-chains gives rise to the observed increase in the interlayer lattice spacing of (CnH2n+1NH3)(2)PbI4 at the melting transition.

Overview No.144 - Mechanical behavior of amorphous alloys

The mechanical properties of amorphous alloys have proven both scientifically unique and of potential practical interest, although the underlying deformation physics of these materials remain less firmly established as compared with crystalline alloys. In this article, we review recent advances in understanding the mechanical behavior of metallic glasses, with particular emphasis on the deformation and fracture mechanisms. Atomistic as well as continuum modeling and experimental work on elasticity, plastic flow and localization, fracture and fatigue are all discussed, and theoretical developments are connected, where possible, with macroscopic experimental responses. The role of glass structure on mechanical properties, and conversely, the effect of deformation upon glass structure, are also described. The mechanical properties of metallic glass-derivative materials – including in situ and ex situ composites, foams and nanocrystal-reinforced glasses – are reviewed as well. Finally, we identify a number of important unresolved issues for the field.

Macro fiber composite (MFC) as a delamination sensor in antisymmetric laminates

The change in extension-twist Coupling due to delamination in antisymmetric laminates is experimentally measured. Experimental results are compared with the results from analytical expression existing in literature and finite element analysis. The application of the Macro-Fiber Composite (MFC) developed at the NASA Langley Research Center for sensing the delamination in the laminates is investigated. While many applications have been reported in the literature using the MFC as an actuator, here its use as a twist sensor has been studied. The real-life application envisaged is structural health monitoring of laminated composite flexbeams taking advantage of the symmetry in the structure. Apart from the defect detection under symmetric conditions, other methods of health monitoring for the same structure are reported for further validation. Results show that MFC works well as a sensor.

Optimization of friction stir welding parameters for dissimilar aluminum alloys

Assembly consisting of cast and wrought aluminum alloys has wide spread application in defense and aero space industries. For the efficacious use of the transition joints, the weld should have adequate strength and formability. In the present investigation, A356 and 6061 aluminum alloys were friction stir welded under tool rotational speed of 1000-1400 rpm and traversing speed of 80-240 mm/min, keeping other parameters same. The variable process window is responsible for the change in total heat input and cooling rate during welding. Structural characterization of the bonded assemblies exhibits recovery-recrystallization in the stirring zone and breaking of coarse eutectic network of Al-Si. Dispersion of fine Si rich particles, refinement of 6061 grain size, low residual stress level and high defect density within weld nugget contribute towards the improvement in bond strength. Lower will be the tool rotational and traversing speed, more dominant will be the above phenomena. Therefore, the joint fabricated using lowest tool traversing and rotational speed, exhibits substantial improvement in bond strength (similar to 98% of that of 6061 alloy), which is also maximum with respect to others. (C) 2010 Elsevier Ltd. All rights reserved.

Suppression of spinel formation to induce reversible thermal behavior in the layered double hydroxides (LDHs) of Co with Al, Fe, Ga, and In

The layered double hydroxides (LDHs) of Co with trivalent cations decompose irreversibly to yield oxides with the spinel structure. Spinel formation is aided by the oxidation of Co(II) to Co(III) in the ambient atmosphere. When the decomposition is carried out under N-2, the oxidation of Co(II) is suppressed, and the resulting oxide has the rock salt structure. Thus, the Co-Al-CO32-/Cl- LDHs yield oxides of the type Co1- Al-x(2x/3)rectangle O-x/3, which are highly metastable, given the large defect concentration. This defect oxide rapidly reverts back to the original hydroxide on soaking in a Na2CO3 solution. Interlayer NO3- anions, on the other hand, decompose generating a highly oxidizing atmosphere, whereby the Co-Al-NO3- LDH decomposes to form the spinel phase even in a N-2 atmosphere. The oxide with the defect rock salt structure formed by the thermal decomposition of the Co-Fe-CO32- LDH under N2, on soaking in a Na2CO3 solution, follows a different kinetic pathway and undergoes a solution transformation into the inverse spinel Co(Co, Fe)(2)O-4. Fe3+ has a low octahedral crystal field stabilization energy and therefore prefers the tetrahedral coordination offered by the structure of the inverse spinel rather than the octahedral coordination of the parent LDH. Similar considerations do not hold in the case of Ga- and In-containing LDHs, given the considerable barriers to the diffusion of M3+ (M=Ga, In) from octahedral to tetrahedral sites owing to their large size. Consequently, the In-containing oxide residue reverts back to the parent hydroxide, whereas this reconstruction is partial in the case of the Ga-containing oxide. These studies show that the reversible thermal behavior offers a competing kinetic pathway to spinel formation. Suppression of the latter induces the reversible behavior in an LDH that otherwise decomposes irreversibly to the spinel.

Localised modes due to impurities in III–V semiconductors

Local mode frequencies due to substitutional impurities in some III–V semiconductors are calculated using Green functions on the mass defect approximation and compared with experimental results.

Phase transformation, oxidation and hydration of raseodymium sesquioxide

The C→A transformation of Pr2O3 has been examined. The transformation is sluggish and takes place at and above 750°C. C-Pr2O3 is oxidized at a lower temperature than the A-form and oxidation proceeds in two stages in both cases. A hydrate Pr 2O3 . H2O, formed by the interaction of C-Pr2O3 and water at relatively high temperature, is described and its behaviour compared with that of Pr(OH)3. The C-form undergoes hydration at 40°C while the A-form does not. All the observations have been explained in terms of the defect structure of C-Pr2O 3.

A nonlinear theory of wake development

A simple new series, using an expansion of the velocity profile in parabolic cylinder functions, has been developed to describe the nonlinear evolution of a steady, laminar, incompressible wake from a given arbitrary initial profile. The first term in this series is itself found to provide a very satisfactory prediction of the decay of the maximum velocity defect in the wake behind a flat plate or aft of the recirculation zone behind a symmetric blunt body. A detailed analysis, including higher order terms, has been made of the flat plate wake with a Blasius profile at the trailing edge. The same method yields, as a special case, complete results for the development of linearized wakes with arbitrary initial profile under the influence of arbitrary pressure gradients. Finally, for purposes of comparison, a simple approximate solution is obtained using momentum integral methods, and found to predict satisfactorily the decay of the maximum velocity defect. © 1970 Wolters-Noordhoff Publishing.

Salmonella typhimurium mutants affecting establishment of lysogeny

Salmonella typhimurium mutants have been isolated in which phage P22 failsto establish lysogeny. These appear to be defective in cAMP metabolism. A phage mutation overcoming the bacterial defect has been mapped between gene c 1 and gene