A spectral multiscale method for wave propagation analysis: Atomistic-continuum coupled simulation

In this paper, we present a new multiscale method which is capable of coupling atomistic and continuum domains for high frequency wave propagation analysis. The problem of non-physical wave reflection, which occurs due to the change in system description across the interface between two scales, can be satisfactorily overcome by the proposed method. We propose an efficient spectral domain decomposition of the total fine scale displacement along with a potent macroscale equation in the Laplace domain to eliminate the spurious interfacial reflection. We use Laplace transform based spectral finite element method to model the macroscale, which provides the optimum approximations for required dynamic responses of the outer atoms of the simulated microscale region very accurately. This new method shows excellent agreement between the proposed multiscale model and the full molecular dynamics (MD) results. Numerical experiments of wave propagation in a 1D harmonic lattice, a 1D lattice with Lennard-Jones potential, a 2D square Bravais lattice, and a 2D triangular lattice with microcrack demonstrate the accuracy and the robustness of the method. In addition, under certain conditions, this method can simulate complex dynamics of crystalline solids involving different spatial and/or temporal scales with sufficient accuracy and efficiency. (C) 2014 Elsevier B.V. All rights reserved.

Space-Time Coded Spatial Modulated Physical Layer Network Coding for Two-Way Relaying

Using the spatial modulation approach, where only one transmit antenna is active at a time, we propose two transmission schemes for two-way relay channel using physical layer network coding with space time coding using coordinate interleaved orthogonal designs (CIODs). It is shown that using two uncorrelated transmit antennas at the nodes, but using only one RF transmit chain and space-time coding across these antennas can give a better performance without using any extra resources and without increasing the hardware implementation cost and complexity. In the first transmission scheme, two antennas are used only at the relay, adaptive network coding (ANC) is employed at the relay and the relay transmits a CIOD space time block code (STBC). This gives a better performance compared to an existing ANC scheme for two-way relay channel which uses one antenna each at all the three nodes. It is shown that for this scheme at high SNR the average end-to-end symbol error probability (SEP) is upper bounded by twice the SEP of a point-to-point fading channel. In the second transmission scheme, two transmit antennas are used at all the three nodes, CIOD STBCs are transmitted in multiple access and broadcast phases. This scheme provides a diversity order of two for the average end-to-end SEP with an increased decoding complexity of O(M-3) for an arbitrary signal set and O(M-2 root M) for square QAM signal set. Simulation results show that the proposed schemes performs better than the existing ANC schemes under perfect and imperfect channel state information.

Experimentally validated Raman Monte Carlo simulation for a cuboid object to obtain Raman spectroscopic signatures for hidden material

In conventional Raman spectroscopic measurements of liquids or surfaces the preferred geometry for detection of the Raman signal is the backscattering (or reflection) mode. For non-transparent layered materials, sub-surface Raman signals have been retrieved using spatially offset Raman spectroscopy (SORS), usually with light collection in the same plane as the point of excitation. However, as a result of multiple scattering in a turbid medium, Raman photons will be emitted in all directions. In this study, Monte Carlo simulations for a three-dimensional layered sample with finite geometry have been performed to confirm the detectability of Raman signals at all angles and at all sides of the object. We considered a non-transparent cuboid container (high density polyethylene) with explosive material (ammonium nitrate) inside. The simulation results were validated with experimental Raman intensities. Monte Carlo simulation results reveal that the ratio of sub-surface to surface signals improves at geometries other than backscattering. In addition, we demonstrate through simulations the effects of the absorption and scattering coefficients of the layers, and that of the diameter of the excitation beam. The advantage of collecting light from all possible 4 angles, over other collection modes, is that this technique is not geometry specific and molecular identification of layers underneath non-transparent surfaces can be obtained with minimal interference from the surface layer. To what extent all sides of the object will contribute to the total signal will depend on the absorption and scattering coefficients and the physical dimensions. Copyright (c) 2015 John Wiley & Sons, Ltd.

Influence of fuel hydrogen fraction on syngas fueled SI engine: Fuel thermo-physical property analysis and in-cylinder experimental investigations

Hydrogen, either in pure form or as a gaseous fuel mixture specie enhances the fuel conversion efficiency and reduce emissions in an internal combustion engine. This is due to the reduction in combustion duration attributed to higher laminar flame speeds. Hydrogen is also expected to increase the engine convective heat flux, attributed (directly or indirectly) to parameters like higher adiabatic flame temperature, laminar flame speed, thermal conductivity and diffusivity and lower flame quenching distance. These factors (adversely) affect the thermo-kinematic response and offset some of the benefits. The current work addresses the influence of mixture hydrogen fraction in syngas on the engine energy balance and the thermo-kinematic response for close to stoichiometric operating conditions. Four different bio-derived syngas compositions with fuel calorific value varying from 3.14 MJ/kg to 7.55 MJ/kg and air fuel mixture hydrogen fraction varying from 7.1% to 14.2% by volume are used. The analysis comprises of (a) use of chemical kinetics simulation package CHEMKIN for quantifying the thermo-physical properties (b) 0-D model for engine in-cylinder analysis and (c) in-cylinder investigations on a two-cylinder engine in open loop cooling mode for quantifying the thermo-kinematic response and engine energy balance. With lower adiabatic flame temperature for Syngas, the in-cylinder heat transfer analysis suggests that temperature has little effect in terms of increasing the heat flux. For typical engine like conditions (700 K and 25 bar at CR of 10), the laminar flame speed for syngas exceeds that of methane (55.5 cm/s) beyond mixture hydrogen fraction of 11% and is attributed to the increase in H based radicals. This leads to a reduction in the effective Lewis number and laminar flame thickness, potentially inducing flame instability and cellularity. Use of a thermodynamic model to assess the isolated influence of thermal conductivity and diffusivity on heat flux suggests an increase in the peak heat flux between 2% and 15% for the lowest (0.420 MW/m(2)) and highest (0.480 MW/m(2)) hydrogen containing syngas over methane (0.415 MW/m(2)) fueled operation. Experimental investigations indicate the engine cooling load for syngas fueled engine is higher by about 7% and 12% as compared to methane fueled operation; the losses are seen to increase with increasing mixture hydrogen fraction. Increase in the gas to electricity efficiency is observed from 18% to 24% as the mixture hydrogen fraction increases from 7.1% to 9.5%. Further increase in mixture hydrogen fraction to 14.2% results in the reduction of efficiency to 23%; argued due to the changes in the initial and terminal stages of combustion. On doubling of mixture hydrogen fraction, the flame kernel development and fast burn phase duration decrease by about 7% and 10% respectively and the terminal combustion duration, corresponding to 90%-98% mass burn, increases by about 23%. This increase in combustion duration arises from the cooling of the near wall mixture in the boundary layer attributed to the presence of hydrogen. The enhancement in engine cooling load and subsequent reduction in the brake thermal efficiency with increasing hydrogen fraction is evident from the engine energy balance along with the cumulative heat release profiles. Copyright (C) 2015, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.

Physical insights into salicylic acid release from poly(anhydrides)

Salicylic acid (SA) based biodegradable polyanhydrides (PAHs) are of great interest for drug delivery in a variety of diseases and disorders owing to the multi-utility of SA. There is a need for the design of SA-based PAHs for tunable drug release, optimized for the treatment of different diseases. In this study, we devised a simple strategy for tuning the release properties and erosion kinetics of a family of PAHs. PAHs incorporating SA were derived from related aliphatic diacids, varying only in the chain length, and prepared by simple melt condensation polymerization. Upon hydrolysis induced erosion, the polymer degrades into cytocompatible products, including the incorporated bioactive SA and diacid. The degradation follows first order kinetics with the rate constant varying by nearly 25 times between the PAH obtained with adipic acid and that with dodecanedioic acid. The release profiles have been tailored from 100% to 50% SA release in 7 days across the different PAHs. The release rate constants of these semi-crystalline, surface eroding PAHs decreased almost linearly with an increase in the diacid chain length, and varied by nearly 40 times between adipic acid and dodecanedioic acid PAH. The degradation products with SA concentration in the range of 30-350 ppm were used to assess cytocompatibility and showed no cytotoxicity to HeLa cells. This particular strategy is expected to (a) enable synthesis of application specific PAHs with tunable erosion and release profiles; (b) encompass a large number of drugs that may be incorporated into the PAH matrix. Such a strategy can potentially be extended to the controlled release of other drugs that may be incorporated into the PAH backbone and has important implications for the rational design of drug eluting bioactive polymers.

Investigations on the physical origin of lateral photovoltage in PbS-colloidal quantum dot/Si heterojunctions

Restricted area heterojunctions, an array of lead sulfide colloidal quantum dots (PbS-CQDs) and crystalline silicon, are studied with a non-destructive remote contact light beam induced current (RC-LBIC) technique. As well as getting good quality active area images we observed an anomalous unipolar signal response for the PbS-CQD/n-Si devices and a conventionally expected bipolar signal profile for the PbS-CQD/p-Si devices. Interestingly, our simulation results consistently yielded a unipolar and bipolar nature in the signals related to the PbSCQD/n-Si and PbS-CQD/p-Si heterostructures, respectively. In order to explain the physical mechanism involved in the unipolar signal response of the PbS-CQD/n-Si devices, we propose a model based on the band alignment in the heterojunctions, in addition to the distribution of photo-induced excess majority carriers across the junction. Given that the RC-LBIC technique is well suited to this context, the presence of these two distinct mechanisms (the bipolar and unipolar nature of the signals) needs to be considered in order to have a better interpretation of the data in the characterization of an array of homo/heterojunctions.

Atomistic Simulation of Stacked Nucleosome Core Particles: Tail Bridging, the H4 Tail, and Effect of Hydrophobic Forces

We report the first atomistic simulation of two stacked nucleosome core particles (NCPs), with an aim to understand, in molecular detail, how they interact, the effect of salt concentration, and how different histone tails contribute to their interaction, with a special emphasis on the H4 tail, known to have the largest stabilizing effect on the NCP-NCP interaction. We do not observe specific K16-mediated interaction between the H4 tail and the H2A-H2B acidic patch, in contrast with the findings from crystallographic studies, but find that the stacking was stable even in the absence of this interaction. We perform simulations with the H4 tail (partially/completely) removed and find that the region between LYS-16 and LYS-20 of the H4 tail holds special importance in mediating the inter-NCP interaction. Performing similar tail-clipped simulations with the H3 tail removed, we compare the roles of the H3 and H4 tails in maintaining the stacking. We discuss the relevance of our simulation results to the bilayer and other liquid-crystalline phases exhibited by NCPs in vitro and, through an analysis of the histone-histone interface, identify the interactions that could possibly stabilize the inter-NCP interaction in these columnar mesophases. Through the mechanical disruption of the stacked nucleosome system using steered molecular dynamics, we quantify the strength of inter-NCP stacking in the presence and absence of salt. We disrupt the stacking at some specific sites of internucleosomal tail-DNA contact and perform a comparative quantification of the binding strengths of various tails in stabilizing the stacking. We also examine how hydrophobic interactions may contribute to the overall stability of the stacking and find a marked difference in the role of hydrophobic forces as compared with electrostatic forces in determining the stability of the stacked nucleosome system.

Intensity Distribution Design for Laser-Induced Thermal Loading Based on Numerical Simulation

To accomplish laser-induced thermal loading simulation tests for pistons,the Gaussian beam was modulated into multi-circular beam with specific intensity distribution.A reverse method was proposed to design the intensity distribution for the laser-induced thermal loading based on finite element(FE) analysis.Firstly,the FE model is improved by alternating parameters of boundary conditions and thermal-physical properties of piston material in a reasonable range,therefore it can simulate the experimental resul...

Numerical Simulation of the Flow Field and Concentration Distribution in the Bulk Growth of Silicon Carbide Crystals

The physical vapor transport (PVT) method is being widely used to grow large-size single SiC crystals. The growth process is associated with heat and mass transport in the growth chamber, chemical reactions among multiple species as well as phase change at the crystal/gas interface. The current paper aims at studying and verifying the transport mechanism and growth kinetics model by demonstrating the flow field and species concentration distribution in the growth system. We have developed a coupled model, which takes into account the mass transport and growth kinetics. Numerical simulation is carried out by employing an in-house developed software based on finite volume method. The results calculated are in good agreement with the experimental observation.

Simulation of the Load-Unload Response Ratio and Critical Sensitivity in the Lattice Solid Model

The Load-Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It involves calculating the ratio of a specified energy release measure during loading and unloading where loading and unloading periods are determined from the earth tide induced perturbations in the Coulomb Failure Stress on optimally oriented faults. In the lead-up to large earthquakes, high LURR values are frequently observed a few months or years prior to the event. These signals may have a similar origin to the observed accelerating seismic moment release (AMR) prior to many large earthquakes or may be due to critical sensitivity of the crust when a large earthquake is imminent. As a first step towards studying the underlying physical mechanism for the LURR observations, numerical studies are conducted using the particle based lattice solid model (LSM) to determine whether LURR observations can be reproduced. The model is initialized as a heterogeneous 2-D block made up of random-sized particles bonded by elastic-brittle links. The system is subjected to uniaxial compression from rigid driving plates on the upper and lower edges of the model. Experiments are conducted using both strain and stress control to load the plates. A sinusoidal stress perturbation is added to the gradual compressional loading to simulate loading and unloading cycles and LURR is calculated. The results reproduce signals similar to those observed in earthquake prediction practice with a high LURR value followed by a sudden drop prior to macroscopic failure of the sample. The results suggest that LURR provides a good predictor for catastrophic failure in elastic-brittle systems and motivate further research to study the underlying physical mechanisms and statistical properties of high LURR values. The results provide encouragement for earthquake prediction research and the use of advanced simulation models to probe the physics of earthquakes.

Semi-insulating layer for novel high voltage polysilicon thin film transistors

This work describes the deposition and characterisation of semi-insulating oxygen-doped silicon films for the development of high voltage polycrystalline silicon (poly-Si) circuitry on glass. The performance of a novel poly-Si High Voltage Thin Film Transistor (HVTFT) structure, incorporating a layer of semi-insulating material, has been investigated using a two dimensional device simulator. The semi-insulating layer increases the operating voltage of the HVTFT structure by linearising the potential distribution in the device offset region. A glass compatible semi-insulating layer, suitable for HVTFT applications, has been deposited by the Plasma Enhanced Chemical Vapour Deposition (PECVD) technique from silane (SiH4), nitrous oxide (N2O) and helium (He) gas mixtures. The as-deposited films are furnace annealed at 600°C which is the maximum process temperature. By varying the N2O/SiH4 ratio the conductivity of the annealed films can be accurately controlled up to a maximum of around 10-7 Ω-1cm-1. Helium dilution of the reactant gases improves both film uniformity and reproducibility. Raman analysis shows the as-deposited and annealed films to be completely amorphous. A model for the microstructure of these Semi-Insulating Amorphous Oxygen-Doped Silicon (SIAOS) films is proposed to explain the observed physical and electrical properties.

Low-temperature (≤600 °C) semi-insulating oxygen-doped silicon films by the PECVD technique for large-area power applications

This work describes the deposition, annealing and characterisation of semi-insulating oxygen-doped silicon films at temperatures compatible with polysilicon circuitry on glass. The semi-insulating layers are deposited by the plasma enhanced chemical vapour deposition technique from silane (SiH4), nitrous oxide (N2O) and helium (He) gas mixtures at a temperature of 350 °C. The as-deposited films are then furnace annealed at 600 °C which is the maximum process temperature. Raman analysis shows the as-deposited and annealed films to be completely amorphous. The most important deposition variable is the N2O SiH4 gas ratio. By varying the N2O SiH4 ratio the conductivity of the annealed films can be accurately controlled, for the first time, down to a minimum of ≈10-7Ω-1cm-1 where they exhibit a T -1 4 temperature dependence indicative of a hopping conduction mechanism. Helium dilution of the reactant gases is shown to improve both film uniformity and reproducibility. A model for the microstructure of these semi-insulating amorphous oxygen-doped silicon films is proposed to explain the observed physical and electrical properties. © 1995.

Low temperature (≤ 600°C) semi-insulating oxygen-doped silicon films by the PECVD technique for large area power applications

This work describes the annealing and characterisation of semi-insulating oxygen-doped silicon films deposited by the Plasma Enhanced Chemical Vapour Deposition (PECVD) technique from silane (SiH4), nitrous oxide (N2O) and helium (He) gas mixtures. The maximum process temperature is chosen to be compatible with large area polycrystalline silicon (poly-Si) circuitry on glass. The most important deposition variable is shown to be the N2O SiH4 gas ratio. Helium dilution results in improved film uniformity and reproducibility. Raman analysis shows the 'as-deposited' and annealed films to be completely amorphous. A model for the microstructure of these Semi-Insulating Amorphous Oxygen-doped Silicon (SIAOS) films is proposed to explain the observed physical and electrical properties. © 1995.

Thickness And Component Distributions Of Yttrium-Titanium Alloy Films In Electron-Beam Physical Vapor Deposition

Thickness and component distributions of large-area thin films are an issue of international concern in the field of material processing. The present work employs experiments and direct simulation Monte Carlo (DSMC) method to investigate three-dimensional low-density, non-equilibrium jets of yttrium and titanium vapor atoms in an electron-beams physical vapor deposition (EBPVD) system furnished with two or three electron-beams, and obtains their deposition thickness and component distributions onto 4-inch and 6-inch mono-crystal silicon wafers. The DSMC results are found in excellent agreement with our measurements, such as evaporation rates of yttrium and titanium measured in-situ by quartz crystal resonators, deposited film thickness distribution measured by Rutherford backscattering spectrometer (RBS) and surface profilometer and deposited film molar ratio distribution measured by RBS and inductively coupled plasma atomic emission spectrometer (ICP-AES). This can be taken as an indication that a combination of DSMC method with elaborate measurements may be satisfactory for predicting and designing accurately the transport process of EBPVD at the atomic level.

Numerical simulation of Richtmyer-Meshkov instability

The compressible Navier-Stokes equations discretized with a fourth order accurate compact finite difference scheme with group velocity control are used to simulate the Richtmyer-Meshkov (R-M) instability problem produced by cylindrical shock-cylindrical material interface with shock Mach number Ms = 1.2 and density ratio 1:20 (interior density/outer density). Effect of shock refraction, reflection, interaction of the reflected shock with the material interface, and effect of initial perturbation modes on R-M instability are investigated numerically. It is noted that the shock refraction is a main physical mechanism of the initial phase changing of the material surface. The multiple interactions of the reflected shock from the origin with the interface and the R-M instability near the material interface are the reason for formation of the spike-bubble structures. Different viscosities lead to different spike-bubble structure characteristics. The vortex pairing phenomenon is found in the initial double mode simulation. The mode interaction is the main factor of small structures production near the interface.