A fast two-grid and finite section method for a class of integral equations on the real line with application to an acoustic scattering problem in the half-plane

We consider the numerical treatment of second kind integral equations on the real line of the form ∅(s) = ∫_(-∞)^(+∞)▒〖κ(s-t)z(t)ϕ(t)dt,s=R〗 (abbreviated ϕ= ψ+K_z ϕ) in which K ϵ L_1 (R), z ϵ L_∞ (R) and ψ ϵ BC(R), the space of bounded continuous functions on R, are assumed known and ϕ ϵ BC(R) is to be determined. We first derive sharp error estimates for the finite section approximation (reducing the range of integration to [-A, A]) via bounds on (1-K_z )^(-1)as an operator on spaces of weighted continuous functions. Numerical solution by a simple discrete collocation method on a uniform grid on R is then analysed: in the case when z is compactly supported this leads to a coefficient matrix which allows a rapid matrix-vector multiply via the FFT. To utilise this possibility we propose a modified two-grid iteration, a feature of which is that the coarse grid matrix is approximated by a banded matrix, and analyse convergence and computational cost. In cases where z is not compactly supported a combined finite section and two-grid algorithm can be applied and we extend the analysis to this case. As an application we consider acoustic scattering in the half-plane with a Robin or impedance boundary condition which we formulate as a boundary integral equation of the class studied. Our final result is that if z (related to the boundary impedance in the application) takes values in an appropriate compact subset Q of the complex plane, then the difference between ϕ(s)and its finite section approximation computed numerically using the iterative scheme proposed is ≤C_1 [kh log⁡〖(1⁄kh)+(1-Θ)^((-1)⁄2) (kA)^((-1)⁄2) 〗 ] in the interval [-ΘA,ΘA](Θ<1) for kh sufficiently small, where k is the wavenumber and h the grid spacing. Moreover this numerical approximation can be computed in ≤C_2 N log⁡N operations, where N = 2A/h is the number of degrees of freedom. The values of the constants C1 and C2 depend only on the set Q and not on the wavenumber k or the support of z.

Newtonian perturbations on models with matter creation

Creation of cold dark matter (CCDM) can macroscopically be described by a negative pressure, and, therefore, the mechanism is capable to accelerate the Universe, without the need of an additional dark energy component. In this framework, we discuss the evolution of perturbations by considering a Neo-Newtonian approach where, unlike in the standard Newtonian cosmology, the fluid pressure is taken into account even in the homogeneous and isotropic background equations (Lima, Zanchin, and Brandenberger, MNRAS 291, L1, 1997). The evolution of the density contrast is calculated in the linear approximation and compared to the one predicted by the Lambda CDM model. The difference between the CCDM and Lambda CDM predictions at the perturbative level is quantified by using three different statistical methods, namely: a simple chi(2)-analysis in the relevant space parameter, a Bayesian statistical inference, and, finally, a Kolmogorov-Smirnov test. We find that under certain circumstances, the CCDM scenario analyzed here predicts an overall dynamics (including Hubble flow and matter fluctuation field) which fully recovers that of the traditional cosmic concordance model. Our basic conclusion is that such a reduction of the dark sector provides a viable alternative description to the accelerating Lambda CDM cosmology.

Generalized Coupled Algebraic Riccati Equations for Discrete-time Markov Jump with Multiplicative Noise Systems

In this paper we consider the existence of the maximal and mean square stabilizing solutions for a set of generalized coupled algebraic Riccati equations (GCARE for short) associated to the infinite-horizon stochastic optimal control problem of discrete-time Markov jump with multiplicative noise linear systems. The weighting matrices of the state and control for the quadratic part are allowed to be indefinite. We present a sufficient condition, based only on some positive semi-definite and kernel restrictions on some matrices, under which there exists the maximal solution and a necessary and sufficient condition under which there exists the mean square stabilizing solution fir the GCARE. We also present a solution for the discounted and long run average cost problems when the performance criterion is assumed be composed by a linear combination of an indefinite quadratic part and a linear part in the state and control variables. The paper is concluded with a numerical example for pension fund with regime switching.

Least energy solutions of a critical Neumann problem with a weight

We investigate the effect of the coefficient of the critical nonlinearity for the Neumann problem on the existence of least energy solutions. As a by-product we establish a Sobolev inequality with interior norm.

On a semilinear Schrodinger equation with critical Sobolev exponent

We consider the semilinear Schrodinger equation -Deltau+V(x)u= K(x) \u \ (2*-2 u) + g(x; u), u is an element of W-1,W-2 (R-N), where N greater than or equal to4, V, K, g are periodic in x(j) for 1 less than or equal toj less than or equal toN, K>0, g is of subcritical growth and 0 is in a gap of the spectrum of -Delta +V. We show that under suitable hypotheses this equation has a solution u not equal 0. In particular, such a solution exists if K equivalent to 1 and g equivalent to 0.

On the nonlinear Neumann problem with critical and supercritical nonlinearities

We investigate the solvability of the Neumann problem (1.1) involving a critical Sobolev exponent. In the first part of this work it is assumed that the coeffcients Q and h are at least continuous. Moreover Q is positive on overline Omega and lambda > 0 is a parameter. We examine the common effect of the mean curvature and the shape of the graphs of the coeffcients Q and h on the existence of low energy solutions. In the second part of this work we consider the same problem with Q replaced by - Q. In this case the problem can be supercritical and the existence results depend on integrability conditions on Q and h.

A demonstration of artificial dissipation effects in some finite volume solution methods for the Navier-Stokes equations

The artificial dissipation effects in some solutions obtained with a Navier-Stokes flow solver are demonstrated. The solvers were used to calculate the flow of an artificially dissipative fluid, which is a fluid having dissipative properties which arise entirely from the solution method itself. This was done by setting the viscosity and heat conduction coefficients in the Navier-Stokes solvers to zero everywhere inside the flow, while at the same time applying the usual no-slip and thermal conducting boundary conditions at solid boundaries. An artificially dissipative flow solution is found where the dissipation depends entirely on the solver itself. If the difference between the solutions obtained with the viscosity and thermal conductivity set to zero and their correct values is small, it is clear that the artificial dissipation is dominating and the solutions are unreliable.

Truncation errors in finite difference models for solute transport equation with first-order reaction

The truncation errors associated with finite difference solutions of the advection-dispersion equation with first-order reaction are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation and a temporally and spatially weighted parametric approach is used for differentiating among the various finite difference schemes. The numerical truncation errors are defined using Peclet and Courant numbers and a new Sink/Source dimensionless number. It is shown that all of the finite difference schemes suffer from truncation errors. Tn particular it is shown that the Crank-Nicolson approximation scheme does not have second order accuracy for this case. The effects of these truncation errors on the solution of an advection-dispersion equation with a first order reaction term are demonstrated by comparison with an analytical solution. The results show that these errors are not negligible and that correcting the finite difference scheme for them results in a more accurate solution. (C) 1999 Elsevier Science B.V. All rights reserved.

Quantum analysis of the nondegenerate optical parametric oscillator with injected signal

In this paper we study the nondegenerate optical parametric oscillator with injected signal, both analytically and numerically. We develop a perturbation approach which allows us to find approximate analytical solutions, starting from the full equations of motion in the positive-P representation. We demonstrate the regimes of validity of our approximations via comparison with the full stochastic results. We find that, with reasonably low levels of injected signal, the system allows for demonstrations of quantum entanglement and the Einstein-Podolsky-Rosen paradox. In contrast to the normal optical parametric oscillator operating below threshold, these features are demonstrated with relatively intense fields.

Robust algorithms for solving stochastic partial differential equations

A robust semi-implicit central partial difference algorithm for the numerical solution of coupled stochastic parabolic partial differential equations (PDEs) is described. This can be used for calculating correlation functions of systems of interacting stochastic fields. Such field equations can arise in the description of Hamiltonian and open systems in the physics of nonlinear processes, and may include multiplicative noise sources. The algorithm can be used for studying the properties of nonlinear quantum or classical field theories. The general approach is outlined and applied to a specific example, namely the quantum statistical fluctuations of ultra-short optical pulses in chi((2)) parametric waveguides. This example uses a non-diagonal coherent state representation, and correctly predicts the sub-shot noise level spectral fluctuations observed in homodyne detection measurements. It is expected that the methods used wilt be applicable for higher-order correlation functions and other physical problems as well. A stochastic differencing technique for reducing sampling errors is also introduced. This involves solving nonlinear stochastic parabolic PDEs in combination with a reference process, which uses the Wigner representation in the example presented here. A computer implementation on MIMD parallel architectures is discussed. (C) 1997 Academic Press.

Validation of plasma clearance of (51)Cr-EDTA in adult renal transplant recipients: comparison with inulin renal clearance

Plasma clearance of (51)Cr-EDTA ((51)Cr-EDTA-Cl) is an alternative method to evaluate glomerular filtration rate (GFR). This study aimed to investigate the concordance between (51)Cr-EDTA-Cl and renal inulin clearance (In-Cl) in renal transplant recipients as well to determine the repeatability of (51)Cr-EDTA-Cl in kidney donors. Forty four kidney recipients and 22 kidney donors were enrolled. Simultaneous measurements of (51)Cr-EDTA-Cl and In-Cl were performed. A single dose of 3.7MBq of (51)Cr-EDTA was injected and the plasma disappearance curve was created by taking blood samples at 2, 4, 6 and 8 h after injection. Bland and Altman statistical approach was used to quantify the agreement between In-Cl and (51)Cr-EDTA-Cl and to determine the better concordance between all possibilities of measure for the (51)Cr-EDTA-Cl. The mean of In-Cl was 44.5 +/- 17.9 ml/min/1.73 m(2). There was a positive correlation between In-Cl and all possible measurements of (51)Cr-EDTA-Cl. (51)Cr-EDTA-Cl with two samples taken at 4 and 8 h or at 4 and 6 h presenting the narrow limits of agreement and a difference (bias) of 2.8 and 2.7 ml/min, respectively. Two plasma sampling for (51)Cr-EDTA-Cl was a reliable method to measure GFR compared with In-Cl and comprises a suitable method to be used in kidney transplanted patients.

Comparison of renal transplantation outcomes in children with and without bladder dysfunction. A customized approach equals the difference

Purpose: We examined the development of urological abnormalities in a group of pediatric renal transplant recipients. Materials and Methods: We reviewed 211 patients younger than 19 years who underwent 226 renal transplants. Three groups of patients were studied-136 children with end stage renal disease due to a nonurological cause (group 1), 56 children with a urological disorder but with an adequate bladder (group 2a) and 19 children with lower urinary tract dysfunction and/or inadequate bladder drainage (group 2b). A total of 15 children in group 2b underwent bladder augmentation (ureterocystoplasty in 6, enterocystoplasty in 9), 2 underwent continent urinary diversion, 1 underwent autoaugmentation and 1 underwent a Mitrofanoff procedure at the bladder for easier drainage. Kidney transplantation was performed in the classic manner by extraperitoneal access, and whenever possible the ureter was reimplanted using an antireflux procedure. Results: At a mean followup of 75 months 13 children had died, 59 grafts were lost and 15 children had received a second transplant. Two patients in group 2a required a complementary urological procedure to preserve renal function (1 enterocystoplasty, 1 vesicostomy). A total of 12 major surgical complications occurred in 226 kidney transplants (5.3%), with a similar incidence in all groups. The overall graft survival at 5 years was 75%, 74% and 84%, respectively, in groups 1, 2a and 2b. Conclusions: With individualized treatment children with severely inferior lower urinary tract function may undergo renal transplantation with a safe and adequate outcome.

Comparison of measured sleeping metabolic rate and predicted basal metabolic rate during the first year of life: evidence of a bias changing with increasing metabolic rate

Objective: To compare measurements of sleeping metabolic rate (SMR) in infancy with predicted basal metabolic rate (BMR) estimated by the equations of Schofield. Methods: Some 104 serial measurements of SMR by indirect calorimetry were performed in 43 healthy infants at 1.5, 3, 6, 9 and 12 months of age. Predicted BMR was calculated using the weight only (BMR-wo) and weight and height (BMR-wh) equations of Schofield for 0-3-y-olds. Measured SMR values were compared with both predictive values by means of the Bland-Altman statistical test. Results: The mean measured SMR was 1.48 MJ/day. The mean predicted BMR values were 1.66 and 1.47 MJ/day for the weight only and weight and height equations, respectively. The Bland-Altman analysis showed that BMR-wo equation on average overestimated SMR by 0.18 MJ/day (11%) and the BMR-wh equation underestimated SMR by 0.01 MJ/day (1%). However the 95% limits of agreement were wide: - 0.64 to - 0.28MJ/day (28%) for the former equation and - 0.39 to +0.41 MJ/day (27%) for the latter equation. Moreover there was a significant correlation between the mean of the measured and predicted metabolic rate and the difference between them. Conclusions: The wide variation seen in the difference between measured and predicted metabolic rate and the bias probably with age indicates there is a need to measure actual metabolic rate for individual clinical care in this age group.