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.

Parallel parsing of programming languages

A new method of specifying the syntax of programming languages, known as hierarchical language specifications (HLS), is proposed. Efficient parallel algorithms for parsing languages generated by HLS are presented. These algorithms run on an exclusive-read exclusive-write parallel random-access machine. They require O(n) processors and O(log2n) time, where n is the length of the string to be parsed. The most important feature of these algorithms is that they do not use a stack.

A comment on the consistency of truncated nonlinear integral equation based theories of freezing

We report the results of two studies of aspects of the consistency of truncated nonlinear integral equation based theories of freezing: (i) We show that the self-consistent solutions to these nonlinear equations are unfortunately sensitive to the level of truncation. For the hard sphere system, if the Wertheim–Thiele representation of the pair direct correlation function is used, the inclusion of part but not all of the triplet direct correlation function contribution, as has been common, worsens the predictions considerably. We also show that the convergence of the solutions found, with respect to number of reciprocal lattice vectors kept in the Fourier expansion of the crystal singlet density, is slow. These conclusions imply great sensitivity to the quality of the pair direct correlation function employed in the theory. (ii) We show the direct correlation function based and the pair correlation function based theories of freezing can be cast into a form which requires solution of isomorphous nonlinear integral equations. However, in the pair correlation function theory the usual neglect of the influence of inhomogeneity of the density distribution on the pair correlation function is shown to be inconsistent to the lowest order in the change of density on freezing, and to lead to erroneous predictions. The Journal of Chemical Physics is copyrighted by The American Institute of Physics.

Studies in Nonlinear Optical Materials: Structure of 3-Menthyl 2-[Bis(dimethylarnino)methylene]-3-oxobutyrate Hemihydrate

C 19Ha4N203.~xH 2 O, Mr= 347.5, monoclinic, C2, a = 15.473 (3), b = 6.963 (2), c = 20.708 (4) ]1, //=108.2(2) ° , V=2119(2)A 3, Z=4, Ox= 1.089 Mg m -3, ,~(Cu Ktx) = 1.5418 ]1, p = 0.523 mm -~, F(000) = 760.0, T= 293 K, R = 0.068 for 1967 unique reflections. The C=C bond length is 1-447 (6)]1, significantly longer than in ethylene, 1.336 (2)]1. The crystal structure is stabilized by O-H...O hydrogen bonding. Explanation for the observed low second-harmonic-generation efficiency (0.5 times that of urea) is provided.

Studies in Nonlinear Optical Materials: Structure of Methyl 2-(4-Ethyl-5-methyl- 2-thioxo-2,3-dihydro-l,3-thiazol- 3-yl) propionate

CIoH15NO282, Mr=245"0, orthorhombic, P21212 ~, a = 6.639 (2), b = 8.205 (2), c = 22.528(6)A, V= I227.2(6)A 3, z=4, Dm= 1.315, Dx= 1.326gem -3, MoKa, 2=0.7107A, 12= 3.63 cm -1, F(000) = 520, T= 293 K, R = 0.037 for 1115 significant reflections. The second-harmonicgeneration (SHG) efficiency of this compound is only 1/10th of the urea standard. The observed low second-order nonlinear response may be attributed to the unfavourable packing of the molecules in the crystal lattice.

A Nonlinear System Under Combined Periodic and Random Excitation

The anharmonic oscillator under combined sinusoidal and white noise excitation is studied using the Gaussian closure approximation. The mean response and the steady-state variance of the system is obtained by the WKBJ approximation and also by the Fokker Planck equation. The multiple steadystate solutions are obtained and their stability analysis is presented. Numerical results are obtained for a particular set of system parameters. The theoretical results are compared with a digital simulation study to bring out the usefulness of the present approximate theory.

Nonlinear I-V characteristics of senmiconducting paints

The oxide-based high resistivity semiconducting Paints can be used on high voltage insulators for improved performance. They were found to exhibit nonlinear I-V characteristics which are reported here. The paints exhibited power law relationship and a sharp transition in merhanism of conduction at a critical voltage.

3-D kinematical conservation laws (KCL): Evolution of a surface in R-3 - in particular propagation of a nonlinear wavefront

3-D KCL are equations of evolution of a propagating surface (or a wavefront) Omega(t), in 3-space dimensions and were first derived by Giles, Prasad and Ravindran in 1995 assuming the motion of the surface to be isotropic. Here we discuss various properties of these 3-D KCL.These are the most general equations in conservation form, governing the evolution of Omega(t) with singularities which we call kinks and which are curves across which the normal n to Omega(t) and amplitude won Omega(t) are discontinuous. From KCL we derive a system of six differential equations and show that the KCL system is equivalent to the ray equations of 2, The six independent equations and an energy transport equation (for small amplitude waves in a polytropic gas) involving an amplitude w (which is related to the normal velocity m of Omega(t)) form a completely determined system of seven equations. We have determined eigenvalues of the system by a very novel method and find that the system has two distinct nonzero eigenvalues and five zero eigenvalues and the dimension of the eigenspace associated with the multiple eigenvalue 0 is only 4. For an appropriately defined m, the two nonzero eigenvalues are real when m > 1 and pure imaginary when m < 1. Finally we give some examples of evolution of weakly nonlinear wavefronts.

Nonlinear transport in semiconducting polymers at high carrier densities

Conducting and semiconducting polymers are important materials in the development of printed, flexible, large-area electronics such as flat-panel displays and photovoltaic cells. There has been rapid progress in developing conjugated polymers with high transport mobility required for high-performance field-effect transistors (FETs), beginning(1) with mobilities around 10(-4) cm(2) V-1 s(-1) to a recent report(2) of 1 cm(2) V-1 s(-1) for poly(2,5-bis(3-tetradecylthiophen-2-yl) thieno[3,2-b] thiophene) (PBTTT). Here, the electrical properties of PBTTT are studied at high charge densities both as the semiconductor layer in FETs and in electrochemically doped films to determine the transport mechanism. We show that data obtained using a wide range of parameters (temperature, gate-induced carrier density, source-drain voltage and doping level) scale onto the universal curve predicted for transport in the Luttinger liquid description of the one-dimensional `metal'.

Nonlinear conduction and electrical switching in one-dimensional conductors

Many one-dimensional conductors show pronounced nonlinear electrical conduction. Some of them show very interesting electrical switching from a low conducting state to a high conducting state. Such electrical switching is often associated with memory. These are discussed with particular emphasis on charge transfer complexestmbine-tcnq, tmpd-tcnq, Cs2(tcnq)3,tea-(tcnq) 2 ando-tolidine-iodine.

Optical and nonlinear absorption properties of Na doped ZnO nanoparticle dispersions

We report linear and nonlinear optical properties of the biologically important Na doped ZnO nanoparticle dispersions. Interesting morphological changes involving a spherical to flowerlike transition have been observed with Na doping. Optical absorption measurements show an exciton absorption around 368 nm. Photoluminescence measurements reveal exciton recombination emission, along with shallow and deep trap emissions. The increased intensity of shallow trap emission with Na doping is attributed to oxygen deficiency and shape changes associated with doping. Nonlinear optical measurements show a predominantly two-photon induced, excited state absorption, when excited with 532 nm, 5 ns laser pulses, indicating potential optical limiting applications.

Functional programming systems revisited

Functional Programming (FP) systems are modified and extended to form Nondeterministic Functional Programming (NFP) systems in which nondeterministic programs can be specified and both deterministic and nondeterministic programs can be verified essentially within the system. It is shown that the algebra of NFP programs has simpler laws in comparison with the algebra of FP programs. "Regular" forms are introduced to put forward a disciplined way of reasoning about programs. Finally, an alternative definition of "linear" forms is proposed for reasoning about recursively defined programs. This definition, when used to test the linearity of forms, results in simpler verification conditions than those generated by the original definition of linear forms.

Observability for a Class of Nonlinear Systems

t - N m and sufficient computable conditions are obtained for the obsemabii of systems with linear state equations and polgwmIal outputs. Based on these, initial state reconstmctors are also described.

Some direct numerical approaches to the solution of the singular nonlinear transonic integral equation

The nonlinear singular integral equation of transonic flow is examined, noting that standard numerical techniques are not applicable in solving it. The difficulties in approximating the integral term in this expression were solved by special methods mitigating the inaccuracies caused by standard approximations. It was shown how the infinite domain of integration can be reduced to a finite one; numerical results were plotted demonstrating that the methods proposed here improve accuracy and computational economy.

Coupled amplitude theory of nonlinear surface acoustic waves

The nonlinear propagation characteristics of surface acoustic waves on an isotropic elastic solid have been studied in this paper. The solution of the harmonic boundary value problem for Rayleigh waves is obtained as a generalized Fourier series whose coefficients are proportional to the slowly varying amplitudes of the various harmonics. The infinite set of coupled equations for the amplitudes when solved exhibit an oscillatory slow variation signifying a continuous transfer of energy back and forth among the various harmonics. A conservation relation is derived among all the harmonic amplitudes.