A new estimator to minimize the error due to breathing in the measurement of respiratory impedance

This paper presents a new respiratory impedance estimator to minimize the error due to breathing. Its practical reliability was evaluated in a simulation using realistic signals. These signals were generated by superposing pressure and flow records obtained in two conditions: 1) when applying forced oscillation to a resistance- inertance- elastance (RIE) mechanical model; 2) when healthy subjects breathed through the unexcited forced oscillation generator. Impedances computed (4-32 Hz) from the simulated signals with the new estimator resulted in a mean value which was scarcely biased by the added breathing (errors less than 1 percent in the mean R, I , and E ) and had a small variability (coefficients of variation of R, I, and E of 1.3, 3.5, and 9.6 percent, respectively). Our results suggest that the proposed estimator reduces the error in measurement of respiratory impedance without appreciable extracomputational cost.

Local and global error models to improve uncertainty quantification

In groundwater applications, Monte Carlo methods are employed to model the uncertainty on geological parameters. However, their brute-force application becomes computationally prohibitive for highly detailed geological descriptions, complex physical processes, and a large number of realizations. The Distance Kernel Method (DKM) overcomes this issue by clustering the realizations in a multidimensional space based on the flow responses obtained by means of an approximate (computationally cheaper) model; then, the uncertainty is estimated from the exact responses that are computed only for one representative realization per cluster (the medoid). Usually, DKM is employed to decrease the size of the sample of realizations that are considered to estimate the uncertainty. We propose to use the information from the approximate responses for uncertainty quantification. The subset of exact solutions provided by DKM is then employed to construct an error model and correct the potential bias of the approximate model. Two error models are devised that both employ the difference between approximate and exact medoid solutions, but differ in the way medoid errors are interpolated to correct the whole set of realizations. The Local Error Model rests upon the clustering defined by DKM and can be seen as a natural way to account for intra-cluster variability; the Global Error Model employs a linear interpolation of all medoid errors regardless of the cluster to which the single realization belongs. These error models are evaluated for an idealized pollution problem in which the uncertainty of the breakthrough curve needs to be estimated. For this numerical test case, we demonstrate that the error models improve the uncertainty quantification provided by the DKM algorithm and are effective in correcting the bias of the estimate computed solely from the MsFV results. The framework presented here is not specific to the methods considered and can be applied to other combinations of approximate models and techniques to select a subset of realizations

Front and domain growth in the presence of gravity

Front and domain growth of a binary mixture in the presence of a gravitational field is studied. The interplay of bulk- and surface-diffusion mechanisms is analyzed. An equation for the evolution of interfaces is derived from a time-dependent Ginzburg-Landau equation with a concentration-dependent diffusion coefficient. Scaling arguments on this equation give the exponents of a power-law growth. Numerical integrations of the Ginzburg-Landau equation corroborate the theoretical analysis.

Gauge group and reality conditions in Ashtekar's complex formulation of canonical gravity

We discuss reality conditions and the relation between spacetime diffeomorphisms and gauge transformations in Ashtekars complex formulation of general relativity. We produce a general theoretical framework for the stabilization algorithm for the reality conditions, which is different from Diracs method of stabilization of constraints. We solve the problem of the projectability of the diffeomorphism transformations from configuration-velocity space to phase space, linking them to the reality conditions. We construct the complete set of canonical generators of the gauge group in the phase space which includes all the gauge variables. This result proves that the canonical formalism has all the gauge structure of the Lagrangian theory, including the time diffeomorphisms.

Accurate and efficient simulation of multiphase flow in a heterogeneous reservoir by using error estimate and control in the multiscale finite-volume framework

The multiscale finite-volume (MSFV) method is designed to reduce the computational cost of elliptic and parabolic problems with highly heterogeneous anisotropic coefficients. The reduction is achieved by splitting the original global problem into a set of local problems (with approximate local boundary conditions) coupled by a coarse global problem. It has been shown recently that the numerical errors in MSFV results can be reduced systematically with an iterative procedure that provides a conservative velocity field after any iteration step. The iterative MSFV (i-MSFV) method can be obtained with an improved (smoothed) multiscale solution to enhance the localization conditions, with a Krylov subspace method [e.g., the generalized-minimal-residual (GMRES) algorithm] preconditioned by the MSFV system, or with a combination of both. In a multiphase-flow system, a balance between accuracy and computational efficiency should be achieved by finding a minimum number of i-MSFV iterations (on pressure), which is necessary to achieve the desired accuracy in the saturation solution. In this work, we extend the i-MSFV method to sequential implicit simulation of time-dependent problems. To control the error of the coupled saturation/pressure system, we analyze the transport error caused by an approximate velocity field. We then propose an error-control strategy on the basis of the residual of the pressure equation. At the beginning of simulation, the pressure solution is iterated until a specified accuracy is achieved. To minimize the number of iterations in a multiphase-flow problem, the solution at the previous timestep is used to improve the localization assumption at the current timestep. Additional iterations are used only when the residual becomes larger than a specified threshold value. Numerical results show that only a few iterations on average are necessary to improve the MSFV results significantly, even for very challenging problems. Therefore, the proposed adaptive strategy yields efficient and accurate simulation of multiphase flow in heterogeneous porous media.

Stochastic semiclassical gravity

In the first part of this paper, we show that the semiclassical Einstein-Langevin equation, introduced in the framework of a stochastic generalization of semiclassical gravity to describe the back reaction of matter stress-energy fluctuations, can be formally derived from a functional method based on the influence functional of Feynman and Vernon. In the second part, we derive a number of results for background solutions of semiclassical gravity consisting of stationary and conformally stationary spacetimes and scalar fields in thermal equilibrium states. For these cases, fluctuation-dissipation relations are derived. We also show that particle creation is related to the vacuum stress-energy fluctuations and that it is enhanced by the presence of stochastic metric fluctuations.

Mode decomposition and renormalization in semiclassical gravity

We compute the influence action for a system perturbatively coupled to a linear scalar field acting as the environment. Subtleties related to divergences that appear when summing over all the modes are made explicit and clarified. Being closely connected with models used in the literature, we show how to completely reconcile the results obtained in the context of stochastic semiclassical gravity when using mode decomposition with those obtained by other standard functional techniques.

Cosmic rays as probes of large extra dimensions and TeV gravity

If there are large extra dimensions and the fundamental Planck scale is at the TeV scale, then the question arises of whether ultrahigh energy cosmic rays might probe them. We study the neutrino-nucleon cross section in these models. The elastic forward scattering is analyzed in some detail, hoping to clarify earlier discussions. We also estimate the black hole production rate. We study energy loss from graviton mediated interactions and conclude that they cannot explain the cosmic ray events above the GZK energy limit. However, these interactions could start horizontal air showers with characteristic profile and at a rate higher than in the standard model.

Induced quantum metric fluctuations and the validity of semiclassical gravity

We propose a criterion for the validity of semiclassical gravity (SCG) which is based on the stability of the solutions of SCG with respect to quantum metric fluctuations. We pay special attention to the two-point quantum correlation functions for the metric perturbations, which contain both intrinsic and induced fluctuations. These fluctuations can be described by the Einstein-Langevin equation obtained in the framework of stochastic gravity. Specifically, the Einstein-Langevin equation yields stochastic correlation functions for the metric perturbations which agree, to leading order in the large N limit, with the quantum correlation functions of the theory of gravity interacting with N matter fields. The homogeneous solutions of the Einstein-Langevin equation are equivalent to the solutions of the perturbed semiclassical equation, which describe the evolution of the expectation value of the quantum metric perturbations. The information on the intrinsic fluctuations, which are connected to the initial fluctuations of the metric perturbations, can also be retrieved entirely from the homogeneous solutions. However, the induced metric fluctuations proportional to the noise kernel can only be obtained from the Einstein-Langevin equation (the inhomogeneous term). These equations exhibit runaway solutions with exponential instabilities. A detailed discussion about different methods to deal with these instabilities is given. We illustrate our criterion by showing explicitly that flat space is stable and a description based on SCG is a valid approximation in that case.

Cosmological perturbations from stochastic gravity

In inflationary cosmological models driven by an inflaton field the origin of the primordial inhomogeneities which are responsible for large-scale structure formation are the quantum fluctuations of the inflaton field. These are usually calculated using the standard theory of cosmological perturbations, where both the gravitational and the inflaton fields are linearly perturbed and quantized. The correlation functions for the primordial metric fluctuations and their power spectrum are then computed. Here we introduce an alternative procedure for calculating the metric correlations based on the Einstein-Langevin equation which emerges in the framework of stochastic semiclassical gravity. We show that the correlation functions for the metric perturbations that follow from the Einstein-Langevin formalism coincide with those obtained with the usual quantization procedures when the scalar field perturbations are linearized. This method is explicitly applied to a simple model of chaotic inflation consisting of a Robertson-Walker background, which undergoes a quasi-de Sitter expansion, minimally coupled to a free massive quantum scalar field. The technique based on the Einstein-Langevin equation can, however, deal naturally with the perturbations of the scalar field even beyond the linear approximation, as is actually required in inflationary models which are not driven by an inflaton field, such as Starobinsky¿s trace-anomaly driven inflation or when calculating corrections due to nonlinear quantum effects in the usual inflaton driven models.

Stability of de Sitter spacetime under isotropic perturbations in semiclassical gravity

A spatially flat Robertson-Walker spacetime driven by a cosmological constant is nonconformally coupled to a massless scalar field. The equations of semiclassical gravity are explicitly solved for this case, and a self-consistent de Sitter solution associated with the Bunch-Davies vacuum state is found (the effect of the quantum field is to shift slightly the effective cosmological constant). Furthermore, it is shown that the corrected de Sitter spacetime is stable under spatially isotropic perturbations of the metric and the quantum state. These results are independent of the free renormalization parameters.

Gravity in the Randall-Sundrum brane world

We discuss the weak gravitational field created by isolated matter sources in the Randall-Sundrum brane world. For the case of a single wall of positive tension, the field stays localized near the wall if the source is stationary. We calculate the leading Kaluza-Klein corrections to the linearized gravitational field of a nonrelativistic spherical object, which is different from the Schwarzschild solution at large distances. In the case of two branes of opposite tension, linearized Brans-Dicke (BD) gravity is recovered on either wall, with different BD parameters. On the wall with positive tension the BD parameter is larger than 3000 provided that the separation between walls is larger than 4 times the AdS radius. The gravitational field due to shadow matter is also considered.

Quantization of AdS3 black holes in external fields: semiclassical results from pure gravity

(2+1)-dimensional anti-de Sitter (AdS) gravity is quantized in the presence of an external scalar field. We find that the coupling between the scalar field and gravity is equivalently described by a perturbed conformal field theory at the boundary of AdS3. This allows us to perform a microscopic computation of the transition rates between black hole states due to absorption and induced emission of the scalar field. Detailed thermodynamic balance then yields Hawking radiation as spontaneous emission, and we find agreement with the semiclassical result, including greybody factors. This result also has application to four and five-dimensional black holes in supergravity.