Mesh-free approximations via the error reproducing kernel method and applications to nonlinear systems developing shocks

Error estimates for the error reproducing kernel method (ERKM) are provided. The ERKM is a mesh-free functional approximation scheme [A. Shaw, D. Roy, A NURBS-based error reproducing kernel method with applications in solid mechanics, Computational Mechanics (2006), to appear (available online)], wherein a targeted function and its derivatives are first approximated via non-uniform rational B-splines (NURBS) basis function. Errors in the NURBS approximation are then reproduced via a family of non-NURBS basis functions, constructed using a polynomial reproduction condition, and added to the NURBS approximation of the function obtained in the first step. In addition to the derivation of error estimates, convergence studies are undertaken for a couple of test boundary value problems with known exact solutions. The ERKM is next applied to a one-dimensional Burgers equation where, time evolution leads to a breakdown of the continuous solution and the appearance of a shock. Many available mesh-free schemes appear to be unable to capture this shock without numerical instability. However, given that any desired order of continuity is achievable through NURBS approximations, the ERKM can even accurately approximate functions with discontinuous derivatives. Moreover, due to the variation diminishing property of NURBS, it has advantages in representing sharp changes in gradients. This paper is focused on demonstrating this ability of ERKM via some numerical examples. Comparisons of some of the results with those via the standard form of the reproducing kernel particle method (RKPM) demonstrate the relative numerical advantages and accuracy of the ERKM.

Set theoretic operations on polygons using the scan-grid approach

This paper describes an algorithm to compute the union, intersection and difference of two polygons using a scan-grid approach. Basically, in this method, the screen is divided into cells and the algorithm is applied to each cell in turn. The output from all the cells is integrated to yield a representation of the output polygon. In most cells, no computation is required and thus the algorithm is a fast one. The algorithm has been implemented for polygons but can be extended to polyhedra as well. The algorithm is shown to take O(N) time in the average case where N is the total number of edges of the two input polygons.

A simulation study of the mixing, atomic flow and velocity profiles of crystalline materials during sliding

Self-contained Non-Equilibrium Molecular Dynamics (NEMD) simulations using Lennard-Jones potentials were performed to identify the origin and mechanisms of atomic scale interfacial behavior between sliding metals. The mixing sequence and velocity profiles were compared via MD simulations for three cases, viz.: sell-mated, similar and hard-softvcrystal pairs. The results showed shear instability, atomic scale mixing, and generation of eddies at the sliding interface. Vorticity at the interface suggests that atomic flow during sliding is similar to fluid flow under Kelvin-Helmholtz instability and this is supported by velocity profiles from the simulations. The initial step-function velocity profile spreads during sliding. However the velocity profile does not change much at later stages of the simulation and it eventually stops spreading. The steady state friction coefficient during simulation was monitored as a function of sliding velocity. Frictional behavior can be explained on the basis of plastic deformation and adiabatic effects. The mixing layer growth kinetics was also investigated.

Stress-induced phase transformation and pseudo-elastic/pseudo-plastic recovery in intermetallic Ni-Al nanowires

Extensive molecular dynamics (MD) simulations have been performed in a B2-NiAl nanowire using an embedded atom method (EAM) potential. We show a stress induced B2 -> body-centered-tetragonal (BCT) phase transformation and a novel temperature and cross-section dependent pseudo-elastic/pseudo-plastic recovery from such an unstable BCT phase with a recoverable strain of similar to 30% as compared to 5-8% in polycrystalline materials. Such a temperature and cross-section dependent pseudo-elastic/pseudo-plastic strain recovery can be useful in various interesting applications of shape memory and strain sensing in nanoscale devices. Effects of size, temperature, and strain rate on the structural and mechanical properties have also been analyzed in detail. For a given size of the nanowire the yield stress of both the B2 and the BCT phases is found to decrease with increasing temperature, whereas for a given temperature and strain rate the yield stress of both the B2 and the BCT phase is found to increase with increase in the cross-sectional dimensions of the nanowire. A constant elastic modulus of similar to 80 GPa of the B2 phase is observed in the temperature range of 200-500 K for nanowires of cross-sectional dimensions in the range of 17.22-28.712 angstrom, whereas the elastic modulus of the BCT phase shows a decreasing trend with an increase in the temperature.

Heat pipes-concepts, materials and applications

The heat pipe is an innovative engineering structure characterized by its capacity to transfer large quantities of heat through relatively small cross-sectional areas with very small temperature differences; it also possesses high thermal conductance and low thermal impedance. In recent times, heat pipes in various forms and designs have found a wide variety of applications. This paper briefly presents the basic concepts of heat pipes, recent innovations in design and their applications.

Effects of Isothermal and Adiabatic Thermal Loadings on Size and Strain Rate Dependence of Copper Nanowire

In the present paper, the size and strain rate effects on ultra-thin < 100 >/{100} Cu nanowires at an initial temperature of 10 K have been discussed. Extensive molecular dynamics (MD) simulations have been performed using Embedded atom method (EAM) to investigate the structural behaviours and properties under high strain rate. Velocity-Verlet algorithm has been used to solve the equation of motions. Two different thermal loading cases have been considered: (i) Isothermal loading, in which Nose-Hoover thermostat is used to maintain the constant system temperature, and (ii) Adiabatic loading, i.e., without any thermostat. Five different wire cross-sections were considered ranging from 0.723 x 0.723 nm(2) to 2.169 x 2.169 nm(2) The strain rates used in the present study were 1 x 10(9) s(-1), 1 x 10(8) s(-1), and 1 x 10(7) s(-1). The effect of strain rate on the mechanical properties of copper nanowires was analysed, which shows that elastic properties are independent of thermal loading for a given strain rate and cross-sectional dimension of nanowire. It showed a decreasing yield stress and yield strain with decreasing strain rate for a given cross- section. Also, a decreasing yield stress and increasing yield strain were observed for a given strain rate with increasing cross-sectional area. Elastic modulus was found to be similar to 100 GPa, which was independent of processing temperature, strain rate, and size for a given initial temperature. Reorientation of < 100 >/{100} square cross-sectional copper nanowire into a series of stable ultra-thin Pentagon copper nanobridge structures with dia of similar to 1 nm at 10 K was observed under high strain rate tensile loading. The effect of isothermal and adiabatic loading on the formation of such pentagonal nanobridge structure has been discussed.

Investigation of Deep Defects Due to -Particle Irradiation in n-Silicon

Electrical properties of deep defects induced in n-silicon by -particles of about 10 MeV energy at a dose of 1014 and 1015 cm-2 are studied by DLTS. The levels at Ec -0.18 eV, Ec -0.26 eV, and Ec -0.48 eV are identified as A center, V2 (=/-) and V2 (-/0) on the basis of activation energy, electron capture cross section, and annealing behavior. Two other irradiation related levels at Ec -0.28 eV and Ec -0.51 eV could not be related to any known center.

Complexity analysis of EEG in patients with schizophrenia using fractal dimension

We computed Higuchi's fractal dimension (FD) of resting, eyes closed EEG recorded from 30 scalp locations in 18 male neuroleptic-naive, recent-onset schizophrenia (NRS) subjects and 15 male healthy control (HC) subjects, who were group-matched for age. Schizophrenia patients showed a diffuse reduction of FD except in the bilateral temporal and occipital regions, with the reduction being most prominent bifrontally. The positive symptom (PS) schizophrenia subjects showed FD values similar to or even higher than HC in the bilateral temporo-occipital regions, along with a co-existent bifrontal FD reduction as noted in the overall sample of NRS. In contrast, this increase in FD values in the bilateral temporo-occipital region was absent in the negative symptom (NS) subgroup. The regional differences in complexity suggested by these findings may reflect the aberrant brain dynamics underlying the pathophysiology of schizophrenia and its symptom dimensions. Higuchi's method of measuring FD directly in the time domain provides an alternative for the more computationally intensive nonlinear methods of estimating EEG complexity.