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.

Pressure shocks in relativistic viscous heat-conducting fluids

In this paper we have studied the propagation of pressure shocks in viscous, heat-conducting, relativistic fluids. Velocities of wave fronts and growth equations for the strength of the waves are obtained in the case of low and high temperatures with variable transport coefficients. On the basis of numerical integrations the growth equation results have been discussed. In the case of constant transport coefficients and for all admissible values of ratio of specific heats of the fluid, an analytical solution for the velocity of the wave as a function of distance along the normal trajectory to the wave front, has been obtained.

Some Exact-Solutions Describing Unsteady Plane Gas-Flows With Shocks

A new class of exact solutions of plane gasdynamic equations is found which describes piston-driven shocks into non-uniform media. The governing equations of these flows are taken in the coordinate system used earlier by Ustinov, and their similarity form is determined by the method of infinitesimal transformations. The solutions give shocks with velocities which either decay or grown in a finite or infinite time depending on the density distribution in the ambient medium, although their strength remains constant. The results of the present study are related to earlier investigations describing the propagation of shocks of constant strength into non-uniform media.

Influence of the roughness of aggregate surface on the interface bond strength

An experimental investigation on the bond strength of the interface between mortar and aggregate is reported. Composite compact specimens were used for applying Mode I and Mode 11 loading effects. The influence of the type of mortar and type of aggregate and its roughness on the bond strength of the interface has been studied. It has been observed that the bond strength of the interface in tension is significantly low, though the mortars exhibited higher strength. The highest tensile bond strength values have been observed with rough concrete surface with M-13 mortar. The bond strength of the interface in Mode I load depends on the type of aggregate surface and its roughness, and the type of mortar, The bond strength of the interface between mortar M-13 cast against rough concrete in direct tension seems to be about one third of the strength of the mortar. However, it is about 1/20th to 1/10th with the mortar M-12 in sandwiched composite specimens. The bond strength of the interface in shear (Mode IT) significantly increases as the roughness and the phase angle of the aggregate surface increase. The strength of mortar on the interface bond strength has been very significant. The sandwiched composite specimens show relatively low bond strength in Mode I loading. The behavior of the interface in both Mode I and Mode 11 loading effects has been brittle, indicating catastrophic failure. (C) 2002 Elsevier Science Ltd. All rights reserved.

Aggregate download throughput for TCP-controlled long file transfers in a WLAN with multiple STA-AP association rates

We consider several WLAN stations associated at rates r(1), r(2), ... r(k) with an Access Point. Each station (STA) is downloading a long file from a local server, located on the LAN to which the Access Point (AP) is attached, using TCP. We assume that a TCP ACK will be produced after the reception of d packets at an STA. We model these simultaneous TCP-controlled transfers using a semi-Markov process. Our analytical approach leads to a procedure to compute aggregate download, as well as per-STA throughputs, numerically, and the results match simulations very well. (C) 2012 Elsevier B.V. All rights reserved.

Efficient external memory structures for range-aggregate queries

We present external memory data structures for efficiently answering range-aggregate queries. The range-aggregate problem is defined as follows: Given a set of weighted points in R-d, compute the aggregate of the weights of the points that lie inside a d-dimensional orthogonal query rectangle. The aggregates we consider in this paper include COUNT, sum, and MAX. First, we develop a structure for answering two-dimensional range-COUNT queries that uses O(N/B) disk blocks and answers a query in O(log(B) N) I/Os, where N is the number of input points and B is the disk block size. The structure can be extended to obtain a near-linear-size structure for answering range-sum queries using O(log(B) N) I/Os, and a linear-size structure for answering range-MAX queries in O(log(B)(2) N) I/Os. Our structures can be made dynamic and extended to higher dimensions. (C) 2012 Elsevier B.V. All rights reserved.

Transonic flutter, shocks, and parameter effects

There is a drop in the flutter boundary of an aeroelastic system placed in a transonic flow due to compressibility effects and is known as the transonic dip. Viscous effects can shift the lo-cation of the shock and depending on the shock strength the boundary layer may separate leading to changes in the flutter speed. An unsteady Euler flow solver coupled with the structural dynamic equations is used to understand the effect of shock on the transonic dip. The effect of various system parameters such as mass ratio, location of the center of mass, position of the elastic axis, ratio of uncoupled natural frequencies in heave and pitch are also studied. Steady turbulent flow results are presented to demonstrate the effect of viscosity on the location and strength of the shock.

Modeling Time-Varying Aggregate Interference in Cognitive Radio Systems, and Application to Primary Exclusive Zone Design

Accurately characterizing the time-varying interference caused to the primary users is essential in ensuring a successful deployment of cognitive radios (CR). We show that the aggregate interference at the primary receiver (PU-Rx) from multiple, randomly located cognitive users (CUs) is well modeled as a shifted lognormal random process, which is more accurate than the lognormal and the Gaussian process models considered in the literature, even for a relatively dense deployment of CUs. It also compares favorably with the asymptotically exact stable and symmetric truncated stable distribution models, except at high CU densities. Our model accounts for the effect of imperfect spectrum sensing, which depends on path-loss, shadowing, and small-scale fading of the link from the primary transmitter to the CU; the interweave and underlay modes or CR operation, which determine the transmit powers of the CUs; and time-correlated shadowing and fading of the links from the CUs to the PU-Rx. It leads to expressions for the probability distribution function, level crossing rate, and average exceedance duration. The impact of cooperative spectrum sensing is also characterized. We validate the model by applying it to redesign the primary exclusive zone to account for the time-varying nature of interference.

Riemann shock tube: 1D normal shocks in air, simulations and experiments

This paper presents numerical simulation of the evolution of one-dimensional normal shocks, their propagation, reflection and interaction in air using a single diaphragm Riemann shock tube and validate them using experimental results. Mathematical model is derived for one-dimensional compressible flow of viscous and conducting medium. Dimensionless form of the mathematical model is used to construct space-time finite element processes based on minimization of the space-time residual functional. The space-time local approximation functions for space-time p-version hierarchical finite elements are considered in higher order GRAPHICS] spaces that permit desired order of global differentiability of local approximations in space and time. The resulting algebraic systems from this approach yield unconditionally positive-definite coefficient matrices, hence ensure unique numerical solution. The evolution is computed for a space-time strip corresponding to a time increment Delta t and then time march to obtain the evolution up to any desired value of time. Numerical studies are designed using recently invented hand-driven shock tube (Reddy tube) parameters, high/low side density and pressure values, high- and low-pressure side shock tube lengths, so that numerically computed results can be compared with actual experimental measurements.

A stabilization approach for mesh-free simulations of systems developing shocks or extreme strain localizations

A new stabilization scheme, based on a stochastic representation of the discretized field variables, is proposed with a view to reduce or even eliminate unphysical oscillations in the mesh-free numerical simulations of systems developing shocks or exhibiting localized bands of extreme deformation in the response. The origin of the stabilization scheme may be traced to nonlinear stochastic filtering and, consistent with a class of such filters, gain-based additive correction terms are applied to the simulated solution of the system, herein achieved through the element-free Galerkin method, in order to impose a set of constraints that help arresting the spurious oscillations. The method is numerically illustrated through its Applications to inviscid Burgers' equations, wherein shocks may develop as a result of intersections of the characteristics, and to a gradient plasticity model whose response is often characterized by a developing shear band as the external load is gradually increased. The potential of the method in stabilized yet accurate numerical simulations of such systems involving extreme gradient variations in the response is thus brought forth. (C) 2014 Elsevier Ltd. All rights reserved.

Molecular Origin of the Intrinsic Bending Force for Helical Morphology Observed in Chiral Amphiphilic Assemblies: Concentration and Size Dependence

The formation of the helical morphology in monolayers and bilayers of chiral amphiphilic assemblies is believed to be driven at least partly by the interactions at the chiral centers of the amphiphiles. However, a detailed microscopic understanding of these interactions and their relation with the helix formation is still not clear. In this article a study of the molecular origin of the chirality-driven helix formation is presented by calculating, for the first time, the effective pair potential between a pair of chiral molecules. This effective potential depends on the relative sizes of the groups attached to the two chiral centers, on the orientation of the amphiphile molecules, and also on the distance between them. We find that for the mirror-image isomers (in the racemic modification) the minimum energy conformation is a nearly parallel alignment of the molecules. On the other hand, the same for a pair of molecules of one kind of enantiomer favors a tilt angle between them, thus leading to the formation of a helical morphology of the aggregate. The tilt angle is determined by the size of the groups attached to the chiral centers of the pair of molecules considered and in many cases predicted it to be close to 45 degrees. The present study, therefore, provides a molecular origin of the intrinsic bending force, suggested by Helfrich (J. Chem. Phys. 1986, 85, 1085-1087), to be responsible for the formation of helical structure. This effective potential may explain many of the existing experimental results, such as the size and the concentration dependence of the formation of helical morphology. It is further found that the elastic forces can significantly modify the pitch predicted by the chiral interactions alone and that the modified real pitch is close to the experimentally observed value. The present study is expected to provide a starting point for future microscopic studies.

Shock propagation in gas dynamics: explicit form of higher order compatibility conditions

With the use of tensor analysis and the method of singular surfaces, an infinite system of equations can be derived to study the propagation of curved shocks of arbitrary strength in gas dynamics. The first three of these have been explicitly given here. This system is further reduced to one involving scalars only. The choice of dependent variables in the infinite system is quite important, it leads to coefficients free from singularities for all values of the shock strength.

Synthesis of some copper(II)-chelating (dialkylamino)pyridine amphiphiles and evaluation of their esterolytic capacities in cationic micellar media

Three new (dialkylamino)pyridine (DAAP)-based ligand amphiphiles 3-5 have been synthesized. All of the compounds possess a metal ion binding subunit in the form of a 2,6-disubstituted DAAP moiety. In addition, at least one ortho-CH2OH substituent is present in all the ligands. Complex formation by these ligands with various metal ions were examined under micellar conditions, but only complexes with Cu(II) ions showed kinetically potent esterolytic capacities under micellar conditions. Complexes with Cu(II) were prepared in host comicellar cetyltrimethylammonium bromide (CTABr) media at pH 7.6. Individual complexes were characterized by UV-visible absorption spectroscopy and electron paramagnetic resonance spectroscopy. These metallomicelles speed the cleavage of the substrates p-nitrophenyl hexanoate or p-nitrophenyl diphenyl phosphate. To ascertain the nature of the active esterolytic species, the stoichiometries of the respective Cu(II) complexes were determined from the kinetic version of Job's plot. In all the instances, 2:1 complex ligand/Cu(II) ion are the most kinetically competent species. The apparent pK(a) values of the Cu(II)-coordinated hydroxyl groups of the ligands 3, 4, and 5, in the comicellar aggregate, are 7.8, 8.0, and 8.0, respectively, as estimated from the rate constant vs pH: profiles of the ester cleavage reactions. The nucleophilic metallomicellar reagents and the second-order "catalytic" rate constants toward esterolysis of the substrate p-nitrophenyl hexanoate (at 25 degrees C, pH 7.6) are 37.5 for 3, 11.4 for 4, and 13.8 for 5. All catalytic systems comprising the coaggregates of 3, 4, or 5 and CTABr demonstrate turnover behavior in the presence of excess substrate.

X-ray studies on crystalline complexes involving amino acids and peptides. XXXIII. Crystal structures of L- and DL-arginine complexed with oxalic acid and a comparative study of amino acid oxalic acid complexes

The DL- and L-arginine complexes of oxalic acid are made up of zwitterionic positively charged amino acid molecules and semi-oxalate ions. The dissimilar molecules aggregate into separate alternating layers in the former. The basic unit in the arginine layer is a centrosymmetric dimer, while the semi-oxalate ions form hydrogen-bonded strings in their layer. In the L-arginine complex each semi-oxalate ion is surrounded by arginine molecules and the complex can be described as an inclusion compound. The oxalic acid complexes of basic amino acids exhibit a variety of ionization states and stoichiometry. They illustrate the effect of aggregation and chirality on ionization state and stoichiometry, and that of molecular properties on aggregation. The semi-oxalate/oxalate ions tend to be planar, but large departures from planarity are possible. The amino acid aggregation in the different oxalic acid complexes do not resemble one another significantly, but the aggregation of a particular amino acid in its oxalic acid complex tends to have similarities with its aggregation in other structures. Also, semi-oxalate ions aggregate into similar strings in four of the six oxalic acid complexes. Thus, the intrinsic aggregation propensities of individual molecules tend to be retained in the complexes.