Korteweg-De Vries Equation for Non-Linear Ion-Acoustic Waves in a Plasma with Negative Ions

The theoretical analysis, based on the perturbation technique, of ion-acoustic waves in the vicinity of a Korteweg-de Vries (K-dV) equation derived in a plasma with some negative ions has been made. The investigation shows that the negative ions in plasma with isothermal electrons introduced a critical concentration at which the ion-acoustic wave plays an important role of wave-breaking and forming a precursor while the plasma with non-isothermal electrons has no such singular behaviour of the wave. These two distinct features of ion waves lead to an overall different approach of present study of ion-waves. A distinct feature of non-uniform transition from the nonisothermal case to isothermal case has been shown. Few particular plasma models have been chosen to show the characteristics behaviour of the ion-waves existing in different cases

Limit analysis of rectangular R.C. slabs on ground

A method is presented for obtaining lower bound on the carrying capacity of reinforced concrete foundation slab-structures subject to non-uniform contact pressure distributions. Functional approach suggested by Vallance for simply supported square slabs subject to uniform pressure distribution has been extended to simply supported rectangular slabs subject to symmetrical non-uniform pressure distributions. Radial solutions, ideally suited for rotationally symmetric problems, are shown to be adoptable for regular polygonal slabs subject to contact pressure paraboloids with constant edge pressures. The functional approach has been shown to be well suited even when the pressure is varying along the edges.

A Smooth Finite Element Method Based on Reproducing Kernel DMS-Splines

The element-based piecewise smooth functional approximation in the conventional finite element method (FEM) results in discontinuous first and higher order derivatives across element boundaries Despite the significant advantages of the FEM in modelling complicated geometries, a motivation in developing mesh-free methods has been the ease with which higher order globally smooth shape functions can be derived via the reproduction of polynomials There is thus a case for combining these advantages in a so-called hybrid scheme or a `smooth FEM' that, whilst retaining the popular mesh-based discretization, obtains shape functions with uniform C-p (p >= 1) continuity One such recent attempt, a NURBS based parametric bridging method (Shaw et al 2008b), uses polynomial reproducing, tensor-product non-uniform rational B-splines (NURBS) over a typical FE mesh and relies upon a (possibly piecewise) bijective geometric map between the physical domain and a rectangular (cuboidal) parametric domain The present work aims at a significant extension and improvement of this concept by replacing NURBS with DMS-splines (say, of degree n > 0) that are defined over triangles and provide Cn-1 continuity across the triangle edges This relieves the need for a geometric map that could precipitate ill-conditioning of the discretized equations Delaunay triangulation is used to discretize the physical domain and shape functions are constructed via the polynomial reproduction condition, which quite remarkably relieves the solution of its sensitive dependence on the selected knotsets Derivatives of shape functions are also constructed based on the principle of reproduction of derivatives of polynomials (Shaw and Roy 2008a) Within the present scheme, the triangles also serve as background integration cells in weak formulations thereby overcoming non-conformability issues Numerical examples involving the evaluation of derivatives of targeted functions up to the fourth order and applications of the method to a few boundary value problems of general interest in solid mechanics over (non-simply connected) bounded domains in 2D are presented towards the end of the paper

The effect of boundaries on the plane Couette flow of granular materials: a bifurcation analysis

The tendency of granular materials in rapid shear flow to form non-uniform structures is well documented in the literature. Through a linear stability analysis of the solution of continuum equations for rapid shear flow of a uniform granular material, performed by Savage (1992) and others subsequently, it has been shown that an infinite plane shearing motion may be unstable in the Lyapunov sense, provided the mean volume fraction of particles is above a critical value. This instability leads to the formation of alternating layers of high and low particle concentrations oriented parallel to the plane of shear. Computer simulations, on the other hand, reveal that non-uniform structures are possible even when the mean volume fraction of particles is small. In the present study, we have examined the structure of fully developed layered solutions, by making use of numerical continuation techniques and bifurcation theory. It is shown that the continuum equations do predict the existence of layered solutions of high amplitude even when the uniform state is linearly stable. An analysis of the effect of bounding walls on the bifurcation structure reveals that the nature of the wall boundary conditions plays a pivotal role in selecting that branch of non-uniform solutions which emerges as the primary branch. This demonstrates unequivocally that the results on the stability of bounded shear how of granular materials presented previously by Wang et al. (1996) are, in general, based on erroneous base states.

The effect of boundaries on the plane Couette flow of granular materials: a bifurcation analysis

The tendency of granular materials in rapid shear ow to form non-uniform structures is well documented in the literature. Through a linear stability analysis of the solution of continuum equations for rapid shear flow of a uniform granular material, performed by Savage (1992) and others subsequently, it has been shown that an infinite plane shearing motion may be unstable in the Lyapunov sense, provided the mean volume fraction of particles is above a critical value. This instability leads to the formation of alternating layers of high and low particle concentrations oriented parallel to the plane of shear. Computer simulations, on the other hand, reveal that non-uniform structures are possible even when the mean volume fraction of particles is small. In the present study, we have examined the structure of fully developed layered solutions, by making use of numerical continuation techniques and bifurcation theory. It is shown that the continuum equations do predict the existence of layered solutions of high amplitude even when the uniform state is linearly stable. An analysis of the effect of bounding walls on the bifurcation structure reveals that the nature of the wall boundary conditions plays a pivotal role in selecting that branch of non-uniform solutions which emerges as the primary branch. This demonstrates unequivocally that the results on the stability of bounded shear flow of granular materials presented previously by Wang et al. (1996) are, in general, based on erroneous base states.

Two-level Mapping Based Cache Index Selection for Packet Forwarding Engines

Packet forwarding is a memory-intensive application requiring multiple accesses through a trie structure. The efficiency of a cache for this application critically depends on the placement function to reduce conflict misses. Traditional placement functions use a one-level mapping that naively partitions trie-nodes into cache sets. However, as a significant percentage of trie nodes are not useful, these schemes suffer from a non-uniform distribution of useful nodes to sets. This in turn results in increased conflict misses. Newer organizations such as variable associativity caches achieve flexibility in placement at the expense of increased hit-latency. This makes them unsuitable for L1 caches.We propose a novel two-level mapping framework that retains the hit-latency of one-level mapping yet incurs fewer conflict misses. This is achieved by introducing a secondlevel mapping which reorganizes the nodes in the naive initial partitions into refined partitions with near-uniform distribution of nodes. Further as this remapping is accomplished by simply adapting the index bits to a given routing table the hit-latency is not affected. We propose three new schemes which result in up to 16% reduction in the number of misses and 13% speedup in memory access time. In comparison, an XOR-based placement scheme known to perform extremely well for general purpose architectures, can obtain up to 2% speedup in memory access time.

Discrimination of Nuclear Grade Steel Samples Subjected to Different Heat Treatment Procedures Based on Acoustic Emission (AE) Data Profiling at Crack Initiation

The present work is an attempt to study crack initiation in nuclear grade, 9Cr-1Mo ferritic steel using AE as an online NDE tool. Laboratory experiments were conducted on 5 heat treated Compact Tension (CT) specimens made out of nuclear grade 9Cr-1Mo ferritic steel by subjecting them to cyclic tensile load. The CT Specimens were of 12.5 mm thickness. The Acoustic emission test system was setup to acquire the data continuously during the test by mounting AE sensor on one of the surfaces of the specimen. This was done to characterize AE data pertaining to crack initiation and then discriminate the samples in terms of their heat treatment processes based on AE data. The AE signatures at crack initiation could conclusively bring to fore the heat treatment distinction on a sample to sample basis in a qualitative sense.Thus, the results obtained through these investigations establish a step forward in utilizing AE technique as an on-line measurement tool for accurate detection and understanding of crack initiation and its profile in 9Cr-1Mo nuclear grade steel subjected to different processes of heat treatment.

On-Chip Memory Architecture Exploration Framework for DSP Processor-Based Embedded System on Chip

Today's SoCs are complex designs with multiple embedded processors, memory subsystems, and application specific peripherals. The memory architecture of embedded SoCs strongly influences the power and performance of the entire system. Further, the memory subsystem constitutes a major part (typically up to 70%) of the silicon area for the current day SoC. In this article, we address the on-chip memory architecture exploration for DSP processors which are organized as multiple memory banks, where banks can be single/dual ported with non-uniform bank sizes. In this paper we propose two different methods for physical memory architecture exploration and identify the strengths and applicability of these methods in a systematic way. Both methods address the memory architecture exploration for a given target application by considering the application's data access characteristics and generates a set of Pareto-optimal design points that are interesting from a power, performance and VLSI area perspective. To the best of our knowledge, this is the first comprehensive work on memory space exploration at physical memory level that integrates data layout and memory exploration to address the system objectives from both hardware design and application software development perspective. Further we propose an automatic framework that explores the design space identifying 100's of Pareto-optimal design points within a few hours of running on a standard desktop configuration.

On the capacity of quantized gaussian MAC channels with finite input alphabet

In this paper, we investigate the achievable rate region of Gaussian multiple access channels (MAC) with finite input alphabet and quantized output. With finite input alphabet and an unquantized receiver, the two-user Gaussian MAC rate region was studied. In most high throughput communication systems based on digital signal processing, the analog received signal is quantized using a low precision quantizer. In this paper, we first derive the expressions for the achievable rate region of a two-user Gaussian MAC with finite input alphabet and quantized output. We show that, with finite input alphabet, the achievable rate region with the commonly used uniform receiver quantizer has a significant loss in the rate region compared. It is observed that this degradation is due to the fact that the received analog signal is densely distributed around the origin, and is therefore not efficiently quantized with a uniform quantizer which has equally spaced quantization intervals. It is also observed that the density of the received analog signal around the origin increases with increasing number of users. Hence, the loss in the achievable rate region due to uniform receiver quantization is expected to increase with increasing number of users. We, therefore, propose a novel non-uniform quantizer with finely spaced quantization intervals near the origin. For a two-user Gaussian MAC with a given finite input alphabet and low precision receiver quantization, we show that the proposed non-uniform quantizer has a significantly larger rate region compared to what is achieved with a uniform quantizer.

Binaural Signal Processing Motivated Generalized Analytic Signal Construction and AM-FM Demodulation

Binaural hearing studies show that the auditory system uses the phase-difference information in the auditory stimuli for localization of a sound source. Motivated by this finding, we present a method for demodulation of amplitude-modulated-frequency-modulated (AM-FM) signals using a ignal and its arbitrary phase-shifted version. The demodulation is achieved using two allpass filters, whose impulse responses are related through the fractional Hilbert transform (FrHT). The allpass filters are obtained by cosine-modulation of a zero-phase flat-top prototype halfband lowpass filter. The outputs of the filters are combined to construct an analytic signal (AS) from which the AM and FM are estimated. We show that, under certain assumptions on the signal and the filter structures, the AM and FM can be obtained exactly. The AM-FM calculations are based on the quasi-eigenfunction approximation. We then extend the concept to the demodulation of multicomponent signals using uniform and non-uniform cosine-modulated filterbank (FB) structures consisting of flat bandpass filters, including the uniform cosine-modulated, equivalent rectangular bandwidth (ERB), and constant-Q filterbanks. We validate the theoretical calculations by considering application on synthesized AM-FM signals and compare the performance in presence of noise with three other multiband demodulation techniques, namely, the Teager-energy-based approach, the Gabor's AS approach, and the linear transduction filter approach. We also show demodulation results for real signals.

Isotropic finite volume discretization

Finite volume methods traditionally employ dimension by dimension extension of the one-dimensional reconstruction and averaging procedures to achieve spatial discretization of the governing partial differential equations on a structured Cartesian mesh in multiple dimensions. This simple approach based on tensor product stencils introduces an undesirable grid orientation dependence in the computed solution. The resulting anisotropic errors lead to a disparity in the calculations that is most prominent between directions parallel and diagonal to the grid lines. In this work we develop isotropic finite volume discretization schemes which minimize such grid orientation effects in multidimensional calculations by eliminating the directional bias in the lowest order term in the truncation error. Explicit isotropic expressions that relate the cell face averaged line and surface integrals of a function and its derivatives to the given cell area and volume averages are derived in two and three dimensions, respectively. It is found that a family of isotropic approximations with a free parameter can be derived by combining isotropic schemes based on next-nearest and next-next-nearest neighbors in three dimensions. Use of these isotropic expressions alone in a standard finite volume framework, however, is found to be insufficient in enforcing rotational invariance when the flux vector is nonlinear and/or spatially non-uniform. The rotationally invariant terms which lead to a loss of isotropy in such cases are explicitly identified and recast in a differential form. Various forms of flux correction terms which allow for a full recovery of rotational invariance in the lowest order truncation error terms, while preserving the formal order of accuracy and discrete conservation of the original finite volume method, are developed. Numerical tests in two and three dimensions attest the superior directional attributes of the proposed isotropic finite volume method. Prominent anisotropic errors, such as spurious asymmetric distortions on a circular reaction-diffusion wave that feature in the conventional finite volume implementation are effectively suppressed through isotropic finite volume discretization. Furthermore, for a given spatial resolution, a striking improvement in the prediction of kinetic energy decay rate corresponding to a general two-dimensional incompressible flow field is observed with the use of an isotropic finite volume method instead of the conventional discretization. (C) 2014 Elsevier Inc. All rights reserved.

Adaptive executions of hyperbolic block-structured AMR applications on GPU systems

A block-structured adaptive mesh refinement (AMR) technique has been used to obtain numerical solutions for many scientific applications. Some block-structured AMR approaches have focused on forming patches of non-uniform sizes where the size of a patch can be tuned to the geometry of a region of interest. In this paper, we develop strategies for adaptive execution of block-structured AMR applications on GPUs, for hyperbolic directionally split solvers. While effective hybrid execution strategies exist for applications with uniform patches, our work considers efficient execution of non-uniform patches with different workloads. Our techniques include bin-packing work units to load balance GPU computations, adaptive asynchronism between CPU and GPU executions using a knapsack formulation, and scheduling communications for multi-GPU executions. Our experiments with synthetic and real data, for single-GPU and multi-GPU executions, on Tesla S1070 and Fermi C2070 clusters, show that our strategies result in up to a 3.23 speedup in performance over existing strategies.

Mechanism on Progressive Failure of a Faulted Rock Slope Due to Slip-Weakening

Slip-weakening is one of the characteristics of geological materials under certain loadings. Non-uniform rock structure may exist in the vicinity of the slip surface for a rock slope. Some portion of the slip surface may be penetrated but the other not. For the latter case, the crack or the fault surface will undergo shear deformation before it becomes a successive surface under a certain loading. As the slipped portion advances,slip-weakening occurs over a distance behind the crack tip. In the weakening zone, the shear strength will decrease from its peak value to residual friction level. The stress will redistribute along the surface of crack and in the weakening zone. Thus the changed local stress concentration leads the crack to extend and the ratio of penetration of the slip surface to increase. From the view of large-scale for the whole slip surface, the shear strength will decrease due to the damage of interior rock structure, and the faulted rock behaves as a softening material. Such a kind of mechanism performs in a large number of practical landslides in the zones experienced strong earthquakes. It should be noted that the mechanism mentioned above is different from that of the breakage of structural clay,in which the geological material is regarded as a medium containing structural lumps and structural bands. In this paper, the softening behavior of a faulted rock should be regarded as a comprehensive result of the whole complicated process including slip-weakening, redistribution of stress, extension of crack tip, and the penetration of the slip surface. This process is accompanied by progressive failure and abrupt structural damage. The size of slip-weakening zone is related to the undergoing strain. Once the relative slide is initiated (local or integrated), the effect of slip-weakening will behave in a certain length behind the crack tip until the formation of the whole slip surface.

On applications of high-frequency asymptotics in aeroacoustics.

The aim of this paper is to survey a range of applications of high-frequency asymptotic methods in aeroacoustics. Specifically, we are concerned with problems associated with noise generation, propagation and scattering as found in large modern aeroengines. With regard to noise generation, we consider the interaction between high-frequency vortical waves and thin aerofoils, with particular emphasis being placed on the way in which the vortical waves act on the non-uniform mean flow around the aerofoil. A ray-theoretic description of the resulting sound as it propagates along the engine intake is then presented, followed by consideration of the diffraction of these rays by the (possibly asymmetric) intake lip to produce sound in the far field. A range of more detailed possible extensions is also presented.

Effect of imperfections on the yielding of two-dimensional foams

The influence of each of the six different types of morphological imperfection - waviness, non-uniform cell wall thickness, cell-size variations, fractured cell walls, cell-wall misalignments, and missing cells - on the yielding of 2D cellular solids has been studied systematically for biaxial loading. Emphasis is placed on quantifying the knock-down effect of these defects on the hydrostatic yield strength and upon understanding the associated deformation mechanisms. The simulations in the present study indicate that the high hydrostatic strength, characteristic of ideal honeycombs, is reduced to a level comparable with the deviatoric strength by several types of defect. The common source of this large knock-down is a switch in deformation mode from cell wall stretching to cell wall bending under hydrostatic loading. Fractured cell edges produce the largest knock-down effect on the yield strength of 2D foams, followed in order by missing cells, wavy cell edges, cell edge misalignments, Γ Voronoi cells, δ Voronoi cells, and non-uniform wall thickness. A simple elliptical yield function with two adjustable material parameters successfully fits the numerically predicted yield surfaces for the imperfect 2D foams, and shows potential as a phenomenological constitutive law to guide the design of structural components made from metallic foams.