Influence Of Rayleigh Effect Combined With Marangoni Effect On The Onset Of Convection In A Liquid Layer Overlying A Porous Layer

The coupling mechanism of Rayleigh effect and Marangoni effect in a liquid-porous system is investigated using a linear stability analysis. The eigenvalue problem is solved by means of a Chebyshev tau method. Results indicate that there are three coupling modes between the Rayleigh effect and the Marangoni effect for different depth ratios. (C) 2008 Elsevier Ltd. All rights reserved.

Teleseismic array analysis of upper mantle compressional velocity structure

Large quantities of teleseismic short-period seismograms recorded at SCARLET provide travel time, apparent velocity and waveform data for study of upper mantle compressional velocity structure. Relative array analysis of arrival times from distant (30° < Δ < 95°) earthquakes at all azimuths constrains lateral velocity variations beneath southern California. We compare dT/dΔ back azimuth and averaged arrival time estimates from the entire network for 154 events to the same parameters derived from small subsets of SCARLET. Patterns of mislocation vectors for over 100 overlapping subarrays delimit the spatial extent of an east-west striking, high-velocity anomaly beneath the Transverse Ranges. Thin lens analysis of the averaged arrival time differences, called 'net delay' data, requires the mean depth of the corresponding lens to be more than 100 km. Our results are consistent with the PKP-delay times of Hadley and Kanamori (1977), who first proposed the high-velocity feature, but we place the anomalous material at substantially greater depths than their 40-100 km estimate.

Detailed analysis of travel time, ray parameter and waveform data from 29 events occurring in the distance range 9° to 40° reveals the upper mantle structure beneath an oceanic ridge to depths of over 900 km. More than 1400 digital seismograms from earthquakes in Mexico and Central America yield 1753 travel times and 58 dT/dΔ measurements as well as high-quality, stable waveforms for investigation of the deep structure of the Gulf of California. The result of a travel time inversion with the tau method (Bessonova et al., 1976) is adjusted to fit the p(Δ) data, then further refined by incorporation of relative amplitude information through synthetic seismogram modeling. The application of a modified wave field continuation method (Clayton and McMechan, 1981) to the data with the final model confirms that GCA is consistent with the entire data set and also provides an estimate of the data resolution in velocity-depth space. We discover that the upper mantle under this spreading center has anomalously slow velocities to depths of 350 km, and place new constraints on the shape of the 660 km discontinuity.

Seismograms from 22 earthquakes along the northeast Pacific rim recorded in southern California form the data set for a comparative investigation of the upper mantle beneath the Cascade Ranges-Juan de Fuca region, an ocean-continent transit ion. These data consist of 853 seismograms (6° < Δ < 42°) which produce 1068 travel times and 40 ray parameter estimates. We use the spreading center model initially in synthetic seismogram modeling, and perturb GCA until the Cascade Ranges data are matched. Wave field continuation of both data sets with a common reference model confirms that real differences exist between the two suites of seismograms, implying lateral variation in the upper mantle. The ocean-continent transition model, CJF, features velocities from 200 and 350 km that are intermediate between GCA and T7 (Burdick and Helmberger, 1978), a model for the inland western United States. Models of continental shield regions (e.g., King and Calcagnile, 1976) have higher velocities in this depth range, but all four model types are similar below 400 km. This variation in rate of velocity increase with tectonic regime suggests an inverse relationship between velocity gradient and lithospheric age above 400 km depth.

Métodos espectrais aplicados à relatividade numérica: determinação dos dados iniciais

Neste trabalho aplicamos métodos espectrais para a determinação da configuração inicial de três espaços-tempos contendo buracos negros. Para isto apresentamos primeiro a foliação do espaço-tempo em hipersuperfícies tridimensionais espaciais parametrizadas pela função temporal t. Este processo é chamado de decomposição 3+1 [2] [5]. O resultado deste processo são dois conjuntos de equações classificadas em equações de vínculo e evolução [4]. As equações de vínculo podem ser divididas em vínculos Hamiltoniano e dos momentos. Para a obtenção dos dados iniciais dos problemas estudados aqui, apenas a equação de vínculo Hamiltoniano será resolvida numericamente, pois as equações de vínculo dos momentos possuem solução analítica nestes casos. Uma pequena descrição dos métodos espectrais é apresentada, destacando-se os método de Galerkin, método pseudoespectral ou de colocação e método de Tau, que são empregados na resolução das equações de vínculo Hamiltoniano dos problemas estudados. Verificamos que os resultados obtidos neste trabalho superam aqueles produzidos por Kidder e Finn [15], devido a uma escolha diferente das funções de base, que aqui satisfazem uma das condições de contorno.

On a generalized Laguerre operational matrix of fractional integration

A new operationalmatrix of fractional integration of arbitrary order for generalized Laguerre polynomials is derived.The fractional integration is described in the Riemann-Liouville sense.This operational matrix is applied together with generalized Laguerre tau method for solving general linearmultitermfractional differential equations (FDEs).Themethod has the advantage of obtaining the solution in terms of the generalized Laguerre parameter. In addition, only a small dimension of generalized Laguerre operational matrix is needed to obtain a satisfactory result. Illustrative examples reveal that the proposedmethod is very effective and convenient for linear multiterm FDEs on a semi-infinite interval.

Numerical solution of a PDE system with non-linear steady state conditions that translates the air stripping pollutants removal

This work deals with the numerical simulation of air stripping process for the pre-treatment of groundwater used in human consumption. The model established in steady state presents an exponential solution that is used, together with the Tau Method, to get a spectral approach of the solution of the system of partial differential equations associated to the model in transient state.

Spectral solution for the air stripping pollutants removal dynamic model with non linear steady state conditions

This work deals with the numerical simulation of air stripping process for the pre-treatment of groundwater used in human consumption. The model established in steady state presents an exponential solution that is used, together with the Tau Method, to get a spectral approach of the solution of the system of partial differential equations associated to the model in transient state.

A Review of Operational Matrices and Spectral Techniques for Fractional Calculus

Recently, operational matrices were adapted for solving several kinds of fractional differential equations (FDEs). The use of numerical techniques in conjunction with operational matrices of some orthogonal polynomials, for the solution of FDEs on finite and infinite intervals, produced highly accurate solutions for such equations. This article discusses spectral techniques based on operational matrices of fractional derivatives and integrals for solving several kinds of linear and nonlinear FDEs. More precisely, we present the operational matrices of fractional derivatives and integrals, for several polynomials on bounded domains, such as the Legendre, Chebyshev, Jacobi and Bernstein polynomials, and we use them with different spectral techniques for solving the aforementioned equations on bounded domains. The operational matrices of fractional derivatives and integrals are also presented for orthogonal Laguerre and modified generalized Laguerre polynomials, and their use with numerical techniques for solving FDEs on a semi-infinite interval is discussed. Several examples are presented to illustrate the numerical and theoretical properties of various spectral techniques for solving FDEs on finite and semi-infinite intervals.

Generalized binomial Tau-leap method for biochemical kinetics incorporating both delay and intrinsic noise

The delay stochastic simulation algorithm (DSSA) by Barrio et al. [Plos Comput. Biol.2, 117–E (2006)] was developed to simulate delayed processes in cell biology in the presence of intrinsic noise, that is, when there are small-to-moderate numbers of certain key molecules present in a chemical reaction system. These delayed processes can faithfully represent complex interactions and mechanisms that imply a number of spatiotemporal processes often not explicitly modeled such as transcription and translation, basic in the modeling of cell signaling pathways. However, for systems with widely varying reaction rate constants or large numbers of molecules, the simulation time steps of both the stochastic simulation algorithm (SSA) and the DSSA can become very small causing considerable computational overheads. In order to overcome the limit of small step sizes, various τ-leap strategies have been suggested for improving computational performance of the SSA. In this paper, we present a binomial τ- DSSA method that extends the τ-leap idea to the delay setting and avoids drawing insufficient numbers of reactions, a common shortcoming of existing binomial τ-leap methods that becomes evident when dealing with complex chemical interactions. The resulting inaccuracies are most evident in the delayed case, even when considering reaction products as potential reactants within the same time step in which they are produced. Moreover, we extend the framework to account for multicellular systems with different degrees of intercellular communication. We apply these ideas to two important genetic regulatory models, namely, the hes1 gene, implicated as a molecular clock, and a Her1/Her 7 model for coupled oscillating cells.

Distinguishing SUSY scenarios using tau polarisation and neutralino Dark Matter

We discuss rst a method of measuring polarisation at the ILC using the 1{prong hadronic decays of the . We then show in this contribution how a study of the ~sector and particularly use of decay polarisation can oer a very good handle for distinguishing between mSUGRA and a SUSY-GUTs scenario, both of which can give rise to appropriate Dark Matter.

A novel method of determining semiconductor parameters in EBIC and SEBIV modes of SEM

We present a novel method for determining semiconductor parameters such as diffusion length L, lifetime tau and surface recombination velocity S of minority carriers by employing scanning electron microscopy (SEM). This new method is applicable to both electron beam induced current (EBIC and surface electron beam induced voltage (SEBIV) modes in SEM. The quantitative descriptions for EBIC and SEBIV signals are derived. The parameters L, S and tau can be directly extracted from the expressions for EBIC or SEBIV signals and their relaxation characteristics in experiment. As an example, the values of L, S and tau for n-p junction and p-Si crystal are determined by using the novel method in EBIC or SEBIV mode. The carrier diffusion length of a p-Si crystal is determined to be 8.74 mum in SEBIV mode. It is very close to the normal diffusion length of 7.41 mum of this sample. The novel method is proved to be very helpful for the quantitative characterization of semiconductor materials and devices. Especially, the SEBIV mode in SEM shows great potential for investigating semiconductor structures nondestructively.

Study of the conductivity and side chain mobility of a comb-like polymer electrolyte using the dynamic mechanical method

Using poly(styrene-co-maleic anhydride) as a backbone and poly(ethylene glycol) methyl ether (PEGME) with different molecular weights as side chains, three comb-like polymers and their Li salt complexes were synthesized. The dynamic mechanical properties and conductivities were investigated. Results showed that the polymer electrolytes possess two glass transitions: alpha -transition and beta -transition, and the temperature dependence of the ionic conductivity shows WLF (Williams-Landel-Ferry) behavior. Based on the time-temperature equivalence principle, a master curve was constructed by selecting T-beta as reference temperature. The values of the WLF parameters (C-1 and C-2) were obtained and were found to be almost independent of the length of the PEGME side chain and the content of Li salt. By reference to T-0 = 50 degreesC. the relation between log tau (c) and c was found to be linear. The master curves are displaced progressively to higher frequencies as the molecular weight of the side chain is increased. The relation between log tau (n) and the molecular weight of the side chain is also linear. (C) 2001 Elsevier Science B.V. All rights reserved.

Measurement of intrinsic properties of amyloid fibrils by the peak force QNM method

We report the investigation of the mechanical properties of different types of amyloid fibrils by the peak force quantitative nanomechanical (PF-QNM) technique. We demonstrate that this technique correctly measures the Young’s modulus independent of the polymorphic state and the cross-sectional structural details of the fibrils, and we show that values for amyloid fibrils assembled from heptapeptides, a-synuclein, Ab(1–42), insulin, b-lactoglobulin,lysozyme, ovalbumin, Tau protein and bovine serum albumin all fall in the range of 2–4 GPa.

Ab initio calculation of the BCC Fe-Al-Mo (Iron-Aluminum-Molybdenum) phase diagram: Implications for the nature of the tau(2) phase

The metastable phase diagram of the BCC-based ordering equilibria in the Fe-Al-Mo system has been calculated via a truncated cluster expansion, through the combination of Full-Potential-Linear augmented Plane Wave (FP-LAPW) electronic structure calculations and of Cluster Variation Method (CVM) thermodynamic calculations in the irregular tetrahedron approximation. Four isothermal sections at 1750 K, 2000 K, 2250 K and 2500 K are calculated and correlated with recently published experimental data on the system. The results confirm that the critical temperature for the order-disorder equilibrium between Fe(3)Al-D0(3) and FeAl-B2 is increased by Mo additions, while the critical temperature for the FeAl-B2/A2 equilibrium is kept approximately invariant with increasing Mo contents. The stabilization of the Al-rich A2 phase in equilibrium with overstoichiometric B2-(Fe,Mo)Al is also consistent with the attribution of the A2 structure to the tau(2) phase, stable at high temperatures in overstoichiometric B2-FeAl. (C) 2009 Elsevier Ltd. All rights reserved.

THE REDUCTIVE PERTURBATION METHOD AND THE KORTEWEG-DE VRIES HIERARCHY

By using the reductive perturbation method of Taniuti with the introduction of an infinite sequence of slow time variables tau(1), tau(3), tau(5), ..., we study the propagation of long surface-waves in a shallow inviscid fluid. The Korteweg-de Vries (KdV) equation appears as the lowest order amplitude equation in slow variables. In this context, we show that, if the lowest order wave amplitude zeta(0) satisfies the KdV equation in the time tau(3), it must satisfy the (2n+1)th order equation of the KdV hierarchy in the time tau(2n+1), With n = 2, 3, 4,.... AS a consequence of this fact, we show with an explicit example that the secularities of the evolution equations for the higher-order terms (zeta(1), zeta(2),...) of the amplitude can be eliminated when zeta(0) is a solitonic solution to the KdV equation. By reversing this argument, we can say that the requirement of a secular-free perturbation theory implies that the amplitude zeta(0) satisfies the (2n+1)th order equation of the KdV hierarchy in the time tau(2n+1) This essentially means that the equations of the KdV hierarchy do play a role in perturbation theory. Thereafter, by considering a solitary-wave solution, we show, again with an explicit, example that the elimination of secularities through the use of the higher order KdV hierarchy equations corresponds, in the laboratory coordinates, to a renormalization of the solitary-wave velocity. Then, we conclude that this procedure of eliminating secularities is closely related to the renormalization technique developed by Kodama and Taniuti.

Tau-functions and dressing transformations for zero-curvature affine integrable equations

The solutions of a large class of hierarchies of zero-curvature equations that includes Toda- and KdV-type hierarchies are investigated. All these hierarchies are constructed from affine (twisted or untwisted) Kac-Moody algebras g. Their common feature is that they have some special vacuum solutions corresponding to Lax operators lying in some Abelian (up to the central term) subalgebra of g; in some interesting cases such subalgebras are of the Heisenberg type. Using the dressing transformation method, the solutions in the orbit of those vacuum solutions are constructed in a uniform way. Then, the generalized tau-functions for those hierarchies are defined as an alternative set of variables corresponding to certain matrix elements evaluated in the integrable highest-weight representations of g. Such definition of tau-functions applies for any level of the representation, and it is independent of its realization (vertex operator or not). The particular important cases of generalized mKdV and KdV hierarchies as well as the Abelian and non-Abelian affine Toda theories are discussed in detail. © 1997 American Institute of Physics.