An Optimal Order Interior Penalty Discontinuous Galerkin Discretization of the Compressible Navier-Stokes Equations

In this article we propose a new symmetric version of the interior penalty discontinuous Galerkin finite element method for the numerical approximation of the compressible Navier-Stokes equations. Here, particular emphasis is devoted to the construction of an optimal numerical method for the evaluation of certain target functionals of practical interest, such as the lift and drag coefficients of a body immersed in a viscous fluid. With this in mind, the key ingredients in the construction of the method include: (i) An adjoint consistent imposition of the boundary conditions; (ii) An adjoint consistent reformulation of the underlying target functional of practical interest; (iii) Design of appropriate interior-penalty stabilization terms. Numerical experiments presented within this article clearly indicate the optimality of the proposed method when the error is measured in terms of both the L_2-norm, as well as for certain target functionals. Computational comparisons with other discontinuous Galerkin schemes proposed in the literature, including the second scheme of Bassi & Rebay, cf. [11], the standard SIPG method outlined in [25], and an NIPG variant of the new scheme will be undertaken.

Emergent Antiferromagnetism Out Of The Hidden-order" State In Uru2si2: High Magnetic Field Nuclear Magnetic Resonance To 40 T."

Very high field (29)Si-NMR measurements using a fully (29)Si-enriched URu(2)Si(2) single crystal were carried out in order to microscopically investigate the hidden order (HO) state and adjacent magnetic phases in the high field limit. At the lowest measured temperature of 0.4 K, a clear anomaly reflecting a Fermi surface instability near 22 T inside the HO state is detected by the (29)Si shift, (29)K(c). Moreover, a strong enhancement of (29)K(c) develops near a critical field H(c) ≃ 35.6 T, and the ^{29}Si-NMR signal disappears suddenly at H(c), indicating the total suppression of the HO state. Nevertheless, a weak and shifted (29)Si-NMR signal reappears for fields higher than H(c) at 4.2 K, providing evidence for a magnetic structure within the magnetic phase caused by the Ising-type anisotropy of the uranium ordered moments.

SOLITON SOLUTIONS TO SYSTEMS OF COUPLED SCHRODINGER EQUATIONS OF HAMILTONIAN TYPE

We study the existence of positive solutions of Hamiltonian-type systems of second-order elliptic PDE in the whole space. The systems depend on a small parameter and involve a potential having a global well structure. We use dual variational methods, a mountain-pass type approach and Fourier analysis to prove positive solutions exist for sufficiently small values of the parameter.

Graviton emission in the bulk by a simply rotating black hole

In this work, we study the emission of tensor-type gravitational degrees of freedom from a higher-dimensional, simply rotating black hole in the bulk. The decoupled radial part of the corresponding field equation is first solved analytically in the limit of low-energy emitted particles and low-angular momentum of the black hole in order to derive the absorption probability. Both the angular and radial equations are then solved numerically, and the comparison of the analytical and numerical results shows a very good agreement in the low and intermediate energy regimes. By using our exact, numerical results we compute the energy and angular-momentum emission rates and their dependence on the spacetime parameters such as the number of additional spacelike dimensions and the angular momentum of the black hole. Particular care is given to the convergence of our results in terms of the number of modes taken into account in the calculation and the multiplicity of graviton tensor modes that correspond to the same angular-momentum numbers.

Stability of higher dimensional Reissner-Nordstrom-anti-de Sitter black holes

We investigate stability of the D-dimensional Reissner-Nordstrom-anti-de Sitter metrics as solutions of the Einstein-Maxwell equations. We have shown that asymptotically anti-de Sitter (AdS) black holes are dynamically stable for all values of charge and anti-de Sitter radius in D=5,6...11 dimensional space-times. This does not contradict dynamical instability of RNAdS black holes found by Gubser in N=8 gauged supergravity, because the latter instability comes from the tachyon mode of the scalar field, coupled to the system. Asymptotically AdS black holes are known to be thermodynamically unstable for some region of parameters, yet, as we have shown here, they are stable against gravitational perturbations.

Third-order nonlinear spectra and optical limiting of lead oxifluoroborate glasses

We have determined two-photon absorption and nonlinear refraction spectra of the 50BO(1.5) - (50-x)PbF(2) - xPbO glasses (with x = 25, 35, 50 cationic %) at the range of the 470 and 1550 nm. The replacement of fluor atoms by oxygen leads to an increase in the third-order susceptibility, due to the formation of non-bridging oxygens (NBO). The nonlinear index of refraction is one order of magnitude higher than the one for fused silica, and it increases almost twice for the sample with x = 50. This sample has also shown promising features for all-optical switching as well as for optical limiting. (C) 2011 Optical Society of America

Verification of a mixed high-order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high-order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier-Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity-velocity formulation for a two-dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high-order compact and non-compact finite-differences from fourth-order to sixth-order of accuracy. The other scheme used spectral methods instead of finite-difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed-order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright (C) 2009 John Wiley & Sons, Ltd.

A new nitrosyl ruthenium complex nitric oxide donor presents higher efficacy than sodium nitroprusside on relaxation of airway smooth muscle

Nitric oxide (NO) has been demonstrated to be the primary agent in relaxing airways in humans and animals. We investigated the mechanisms involved in the relaxation induced by NO-donors, ruthenium complex [Ru(terpy)(bdq)NO(+)](3+) (TERPY) and sodium nitroprusside (SNP) in isolated trachea of rats contracted with carbachol in an isolated organs chamber. For instance, we verified the contribution of K(+) channels, the importance of sGC/cGMP pathway, the influence of the extra and intracellular Ca(2+) sources and the contribution of the epithelium on the relaxing response. Additionally, we have used confocal microscopy in order to analyze the action of the NO-donors on cytosolic Ca(2+) concentration. The results demonstrated that both compounds led to the relaxation of trachea in a dependent-concentration way. However, the maximum effect (E(max)) of TERPY is higher than the SNP. The relaxation induced by SNP (but not TERPY) was significantly reduced by pretreatment with ODQ (sGC inhibitor). Only TERPY-induced relaxation was reduced by tetraethylammonium (K(+) channels blocker) and by pre-contraction with 75 mM KCl (membrane depolarization). The response to both NO-donors was not altered by the presence of thapsigargin (sarcoplasmic reticulum Ca(2+)-ATPase inhibitor). The epithelium removal has reduced the relaxation only to SNP, and it has no effect on TERPY. The both NO-donors reduced the contraction evoked by Ca(2+) influx, while TERPY have shown a higher inhibitory effect on contraction. Moreover, the TERPY was more effective than SNP in reducing the cytosolic Ca(2+) concentration measured by confocal microscopy. In conclusion, these results show that TERPY induces airway smooth muscle relaxation by cGMP-independent mechanisms, it involves the fluxes of Ca(2+) and K(+) across the membrane, it is more effective in reducing cytosolic Ca(2+) concentration and inducing relaxation in the rat trachea than the standard drug, SNP. (C) 2011 Elsevier B.V. All rights reserved.

Existence of Multiple Solutions for Second Order Boundary Value Problems

We give conditions on f involving pairs of lower and upper solutions which lead to the existence of at least three solutions of the two point boundary value problem y" + f(x, y, y') = 0, x epsilon [0, 1], y(0) = 0 = y(1). In the special case f(x, y, y') = f(y) greater than or equal to 0 we give growth conditions on f and apply our general result to show the existence of three positive solutions. We give an example showing this latter result is sharp. Our results extend those of Avery and of Lakshmikantham et al.

Periodic or Bounded Solutions of Carathéodory Systems of Ordinary Differential Equations

In this paper we extend the guiding function approach to show that there are periodic or bounded solutions for first order systems of ordinary differential equations of the form x1 =f(t,x), a.e. epsilon[a,b], where f satisfies the Caratheodory conditions. Our results generalize recent ones of Mawhin and Ward.

Fixed bed reactors with catalyst poisoning - first order kinetics.

The long performance of an isothermal fixed bed reactor undergoing catalyst poisoning is theoretically analyzed using the dispersion model. First order reaction with dth order deactivation is assumed and the model equations are solved by matched asymptotic expansions for large Peclet number. Simple closed-form solutions, uniformly valid in time, are obtained.

Subcycling first- and second-order generalizations of the trapezoidal rule

Subcycling algorithms which employ multiple timesteps have been previously proposed for explicit direct integration of first- and second-order systems of equations arising in finite element analysis, as well as for integration using explicit/implicit partitions of a model. The author has recently extended this work to implicit/implicit multi-timestep partitions of both first- and second-order systems. In this paper, improved algorithms for multi-timestep implicit integration are introduced, that overcome some weaknesses of those proposed previously. In particular, in the second-order case, improved stability is obtained. Some of the energy conservation properties of the Newmark family of algorithms are shown to be preserved in the new multi-timestep extensions of the Newmark method. In the first-order case, the generalized trapezoidal rule is extended to multiple timesteps, in a simple way that permits an implicit/implicit partition. Explicit special cases of the present algorithms exist. These are compared to algorithms proposed previously. (C) 1998 John Wiley & Sons, Ltd.

Critical values for unit root and cointegration test statistics - the use of response surface equations

This note considers the value of surface response equations which can be used to calculate critical values for a range of unit root and cointegration tests popular in applied economic research.

Existence results for a damped second order abstract functional differential equation with impulses

In this work we study the existence and regularity of mild solutions for a damped second order abstract functional differential equation with impulses. The results are obtained using the cosine function theory and fixed point criterions. (C) 2009 Elsevier Ltd. All rights reserved.

On the loss of regularity for a class of weakly hyperbolic operators

In this paper we determine bounds for the optimal loss of regularity in the Sobolev scale for a class of weakly hyperbolic operators. (C) 2009 Elsevier Inc. All rights reserved.