Some nonstandard error analysis of discontinuous Galerkin methods for elliptic problems

An a priori error analysis of discontinuous Galerkin methods for a general elliptic problem is derived under a mild elliptic regularity assumption on the solution. This is accomplished by using some techniques from a posteriori error analysis. The model problem is assumed to satisfy a GAyenrding type inequality. Optimal order L (2) norm a priori error estimates are derived for an adjoint consistent interior penalty method.

Comparative study of different formulations and solution methods for two phase flow in saturated porous media

Six models (Simulators) are formulated and developed with all possible combinations of pressure and saturation of the phases as primary variables. A comparative study between six simulators with two numerical methods, conventional simultaneous and modified sequential methods are carried out. The results of the numerical models are compared with the laboratory experimental results to study the accuracy of the model especially in heterogeneous porous media. From the study it is observed that the simulator using pressure and saturation of the wetting fluid (PW, SW formulation) is the best among the models tested. Many simulators with nonwetting phase as one of the primary variables did not converge when used along with simultaneous method. Based on simulator 1 (PW, SW formulation), a comparison of different solution methods such as simultaneous method, modified sequential and adaptive solution modified sequential method are carried out on 4 test problems including heterogeneous and randomly heterogeneous problems. It is found that the modified sequential and adaptive solution modified sequential methods could save the memory by half and as also the CPU time required by these methods is very less when compared with that using simultaneous method. It is also found that the simulator with PNW and PW as the primary variable which had problem of convergence using the simultaneous method, converged using both the modified sequential method and also using adaptive solution modified sequential method. The present study indicates that pressure and saturation formulation along with adaptive solution modified sequential method is the best among the different simulators and methods tested.

Latitudinal gradient in aerosol properties over the Indian and Southern Oceans during the austral summer

Extensive measurements of columnar aerosol optical depth (AOD), composite (M-T) and black carbon aerosol mass (M-B) concentrations were made over the tropical Indian and Southern Oceans as a part of the Pilot Expedition to the Southern Ocean during the boreal winter. The AOD, M-T and M-B show large latitudinal gradient towards south up to ITCZ. Beyond ITCZ, up to 56 degrees S, AOD and M-B show very low and steady values. However M-T shows large variations in the Southern Ocean due to the enhanced production of sea salt aerosols associated with high sea surface winds. The short wave aerosol radiative forcing at the surface over north of equator is in the range - 10 to -23 W m(-2), whereas that over the Southern Ocean was in the range -4 to -5 W m(-2). The corresponding atmospheric forcing was in the range of 6-13 W m(-2) and 0.8-1.4 W m(-2). This large north south change in the aerosol radiative forcing has important implications to the meridional circulation and hence to climate.

Integral method for the calculation of incompressible two-dimensional transitional boundary layers

The method proposed here considers the mean flow in the transition zone as a linear combination of the laminar and turbulent boundary layer in proportions determined by the transitional intermittency, the component flows being calculated by approximate integral methods. The intermittency distribution adopted takes into account the possibility of subtransitions within the zone in the presence of strong pressure gradients. A new nondimensional spot formation rate, whose value depends on the pressure gradient, is utilized to estimate the extent of the transition zone. Onset location is determined by a correlation that takes into account freestream turbulence and facility-specific residual disturbances in test data. Extensive comparisons with available experimental results in strong pressure gradients show that the proposed method performs at least as well as differential models, in many cases better, and is always faster.

Gradient solid-electrolyte for use with dissimilar gas electrodes

The EMF of a solid-state cell, incorporating a composite solid-electrolyte with gradual variation in composition, and dissimilar gas electrodes, has been studied as a function of temperature and partial pressures at the electrodes. The cell with the configuration: Pt, CO2' + O2' parallel-to Na2CO3\Na(SO4)x(CO3)1-x\Na2SO4 parallel-to SO3'' + SO2'' + O2'', Pt x=0 x=1 was investigated in the temperature range 973 to 1079 K. The solid-electrolyte surface exposed to SO3 + SO2 + O2 gas mixture was doped-Na2SO4, whereas the CO2 + O2 gas mixture was in contact with pure Na2CO3. The composition of the solid solution between the carbonate and sulfate, with hexagonal structure, was varied gradually between the boundary values. It has been found that the EMF of the cell is close to that calculated from thermodynamic data, assuming unit transport number for Na+ ions. The gradient in the concentration of sulfate and carbonate ions in the electrolyte does not give rise to a significant diffusion potential.

NDDO-based CI methods for the prediction of electronic spectra and sum-over-states molecular hyperpolarization

NDDO-based (AM1) configuration interaction (CI) calculations have been used to calculate the wavelength and oscillator strengths of electronic absorptions in organic molecules and the results used in a sum-over-states treatment to calculate second-order-hyperpolarizabilities. The results for both spectra and hyperpolarizabilities are of acceptable quality as long as a suitable CI-expansion is used. We have found that using an active space of eight electrons in eight orbitals and including all single and pair-double excitations in the CI leads to results that agree well with experiment and that do not change significantly with increasing active space for most organic molecules. Calculated second-order hyperpolarizabilities using this type of CI within a sum-over-states calculation appear to be of useful accuracy.

A study of the synthetic methods and properties of graphenes

Graphenes with varying number of layers can be synthesized by using different strategies. Thus, single-layer graphene is prepared by micromechanical cleavage, reduction of single-layer graphene oxide, chemical vapor deposition and other methods. Few-layer graphenes are synthesized by conversion of nanodiamond, arc discharge of graphite and other methods. In this article, we briefly overview the various synthetic methods and the surface, magnetic and electrical properties of the produced graphenes. Few-layer graphenes exhibit ferromagnetic features along with antiferromagnetic properties, independent of the method of preparation. Aside from the data on electrical conductivity of graphenes and graphene-polymer composites, we also present the field-effect transistor characteristics of graphenes. Only single-layer reduced graphene oxide exhibits ambipolar properties. The interaction of electron donor and acceptor molecules with few-layer graphene samples is examined in detail.

Weighted subspace methods and spatial smoothing: analysis and comparison

The effect of using a spatially smoothed forward-backward covariance matrix on the performance of weighted eigen-based state space methods/ESPRIT, and weighted MUSIC for direction-of-arrival (DOA) estimation is analyzed. Expressions for the mean-squared error in the estimates of the signal zeros and the DOA estimates, along with some general properties of the estimates and optimal weighting matrices, are derived. A key result is that optimally weighted MUSIC and weighted state-space methods/ESPRIT have identical asymptotic performance. Moreover, by properly choosing the number of subarrays, the performance of unweighted state space methods can be significantly improved. It is also shown that the mean-squared error in the DOA estimates is independent of the exact distribution of the source amplitudes. This results in a unified framework for dealing with DOA estimation using a uniformly spaced linear sensor array and the time series frequency estimation problems.

Finite data performance of MUSIC and minimum norm methods

In the direction of arrival (DOA) estimation problem, we encounter both finite data and insufficient knowledge of array characterization. It is therefore important to study how subspace-based methods perform in such conditions. We analyze the finite data performance of the multiple signal classification (MUSIC) and minimum norm (min. norm) methods in the presence of sensor gain and phase errors, and derive expressions for the mean square error (MSE) in the DOA estimates. These expressions are first derived assuming an arbitrary array and then simplified for the special case of an uniform linear array with isotropic sensors. When they are further simplified for the case of finite data only and sensor errors only, they reduce to the recent results given in [9-12]. Computer simulations are used to verify the closeness between the predicted and simulated values of the MSE.

Measurement of complex optical flow with use of an augmented generalized gradient scheme

We present a method for measuring the local velocities and first-order variations in velocities in a timevarying image. The scheme is an extension of the generalized gradient model that encompasses the local variation of velocity within a local patch of the image. Motion within a patch is analyzed in parallel by 42 different spatiotemporal filters derived from 6 linearly independent spatiotemporal kernels. No constraints are imposed on the image structure, and there is no need for smoothness constraints on the velocity field. The aperture problem does not arise so long as there is some two-dimensional structure in the patch being analyzed. Among the advantages of the scheme is that there is no requirement to calculate second or higher derivatives of the image function. This makes the scheme robust in the presence of noise. The spatiotemporal kernels are of simple form, involving Gaussian functions, and are biologically plausible receptive fields. The validity of the scheme is demonstrated by application to both synthetic and real video images sequences and by direct comparison with another recently published scheme Biol. Cybern. 63, 185 (1990)] for the measurement of complex optical flow.

Fuel drop size measurement in a carburetted induction system with 2 throttling methods

A new throttling system far SI engines is examined. The SMD of the fuel droplets in the induction system is measured to evaluate the performance of the new device with respect to the conventional throttle plate arrangement. The measurements are conducted at steady now conditions. A forward angular scattering technique with a He-Ne laser beam is used for droplet size measurement. The experiments are carried out with different mixture strength, stream velocity and throttle positions. It is observed that A/F ratio has no effect on SMD. However, stream velocity and throttle position have a significant influence on SMD. The new throttling method is found to be more effective in reducing the SMD, particularly at low throttle opening and high stream velocity compared to the conventional throttle plate.

Measurement of complex optical flow with use of an augmented generalized gradient scheme

We present a method for measuring the local velocities and first-order variations in velocities in a time-varying image. The scheme is an extension of the generalized gradient model that encompasses the local variation of velocity within a local patch of the image. Motion within a patch is analyzed in parallel by 42 different spatiotemporal filters derived from 6 linearly independent spatiotemporal kernels. No constraints are imposed on the image structure, and there is no need for smoothness constraints on the velocity field. The aperture problem does not arise so long as there is some two-dimensional structure in the patch being analyzed. Among the advantages of the scheme is that there is no requirement to calculate second or higher derivatives of the image function. This makes the scheme robust in the presence of noise. The spatiotemporal kernels are of simple form, involving Gaussian functions, and are biologically plausible receptive fields. The validity of the scheme is demonstrated by application to both synthetic and real video images sequences and by direct comparison with another recently published scheme [Biol. Cybern. 63, 185 (1990)] for the measurement of complex optical flow.

Influence of poling methods on the orientational dynamics of 2-methyl-4-nitro-aniline in poly(methyl methacrylate) studied by second harmonic generation

The magnitude and stability of the induced dipolar orientation of 2-methyl-4-nitroaniline (MNA)/poly(methyl methacrylate) (PMMA) guest/host system is investigated. The chromophores are aligned using both the corona discharge and contact electrode poling techniques. The magnitude of order parameter (also an indicator for the second order nonlinear susceptibility) is measured by recording absorbances of the poled (by the two different techniques) and unpoled PMMA films at different concentrations of MNA. Under the same conditions the corona poling technique creates a higher alignment of molecules along the field direction. The time dependence of the second harmonic intensity of the MNA/PMMA film prepared by the two techniques can be described by a Kohlrausch-Williams-Watts stretched exponential. The temperature dependence of the decay time constant is found to generally follow a modified Williams-Landel-Ferry (WLF) or Vogel-Tamann-Fulcher (VTF) equation. The glass transition temperature seems to be the single most important parameter for determining the relaxation time tau(T).

On generating optimal signal probabilities for random tests: A genetic approach

Genetic Algorithms are robust search and optimization techniques. A Genetic Algorithm based approach for determining the optimal input distributions for generating random test vectors is proposed in the paper. A cost function based on the COP testability measure for determining the efficacy of the input distributions is discussed, A brief overview of Genetic Algorithms (GAs) and the specific details of our implementation are described. Experimental results based on ISCAS-85 benchmark circuits are presented. The performance pf our GA-based approach is compared with previous results. While the GA generates more efficient input distributions than the previous methods which are based on gradient descent search, the overheads of the GA in computing the input distributions are larger. To account for the relatively quick convergence of the gradient descent methods, we analyze the landscape of the COP-based cost function. We prove that the cost function is unimodal in the search space. This feature makes the cost function amenable to optimization by gradient-descent techniques as compared to random search methods such as Genetic Algorithms.