Use of high-resolution geophysical data to characterize heterogeneous aquifers: influence of data integration method on hydrological predictions

The integration of geophysical data into the subsurface characterization problem has been shown in many cases to significantly improve hydrological knowledge by providing information at spatial scales and locations that is unattainable using conventional hydrological measurement techniques. The investigation of exactly how much benefit can be brought by geophysical data in terms of its effect on hydrological predictions, however, has received considerably less attention in the literature. Here, we examine the potential hydrological benefits brought by a recently introduced simulated annealing (SA) conditional stochastic simulation method designed for the assimilation of diverse hydrogeophysical data sets. We consider the specific case of integrating crosshole ground-penetrating radar (GPR) and borehole porosity log data to characterize the porosity distribution in saturated heterogeneous aquifers. In many cases, porosity is linked to hydraulic conductivity and thus to flow and transport behavior. To perform our evaluation, we first generate a number of synthetic porosity fields exhibiting varying degrees of spatial continuity and structural complexity. Next, we simulate the collection of crosshole GPR data between several boreholes in these fields, and the collection of porosity log data at the borehole locations. The inverted GPR data, together with the porosity logs, are then used to reconstruct the porosity field using the SA-based method, along with a number of other more elementary approaches. Assuming that the grid-cell-scale relationship between porosity and hydraulic conductivity is unique and known, the porosity realizations are then used in groundwater flow and contaminant transport simulations to assess the benefits and limitations of the different approaches.

A possible way to relate the covariantization and the negative dimensional integration method in the light-cone gauge

Here we present a possible way to relate the method of covariantizing the gauge-dependent pole and the negative dimensional integration method for computing Feynman integrals pertinent to the light-cone gauge fields. Both techniques are applicable to the algebraic light-cone gauge and dispense with prescriptions to treat the characteristic poles.

Evaluating residues and integrals through negative dimensional integration method (NDIM)

The standard way of evaluating residues and some real integrals through the residue theorem (Cauchy's theorem) is well-known and widely applied in many branches of Physics. Herein we present an alternative technique based on the negative dimensional integration method (NDIM) originally developed to handle Feynman integrals. The advantage of this new technique is that we need only to apply Gaussian integration and solve systems of linear algebraic equations, with no need to determine the poles themselves or their residues, as well as obtaining a whole class of results for differing orders of poles simultaneously.

Power conversion and signal transmission integration method based on dual modulation of DC-DC converters

For the development of communication systems such as Internet of Things, integrating communication with power supplies is an attractive solution to reduce supply cost. This paper presents a novel method of power/signal dual modulation (PSDM), by which signal transmission is integrated with power conversion. This method takes advantage of the intrinsic ripple initiated in switch mode power supplies as signal carriers, by which cost-effective communications can be realized. The principles of PSDM are discussed, and two basic dual modulation methods (specifically PWM/FSK and PWM/PSK) are concluded. The key points of designing a PWM/FSK system, including topology selection, carrier shape, and carrier frequency, are discussed to provide theoretical guidelines. A practical signal modulation-demodulation method is given, and a prototype system provides experimental results to verify the effectiveness of the proposed solution.

Finite integration methods for isospectral flows

In this paper we consider the approximate computation of isospectral flows based on finite integration methods( FIM) with radial basis functions( RBF) interpolation,a new algorithm is developed. Our method ensures the symmetry of the solutions. Numerical experiments demonstrate that the solutions have higher accuracy by our algorithm than by the second order Runge- Kutta( RK2) method.

Parallel Monte Carlo approach for integration of the rendering equation

This paper is addressed to the numerical solving of the rendering equation in realistic image creation. The rendering equation is integral equation describing the light propagation in a scene accordingly to a given illumination model. The used illumination model determines the kernel of the equation under consideration. Nowadays, widely used are the Monte Carlo methods for solving the rendering equation in order to create photorealistic images. In this work we consider the Monte Carlo solving of the rendering equation in the context of the parallel sampling scheme for hemisphere. Our aim is to apply this sampling scheme to stratified Monte Carlo integration method for parallel solving of the rendering equation. The domain for integration of the rendering equation is a hemisphere. We divide the hemispherical domain into a number of equal sub-domains of orthogonal spherical triangles. This domain partitioning allows to solve the rendering equation in parallel. It is known that the Neumann series represent the solution of the integral equation as a infinity sum of integrals. We approximate this sum with a desired truncation error (systematic error) receiving the fixed number of iteration. Then the rendering equation is solved iteratively using Monte Carlo approach. At each iteration we solve multi-dimensional integrals using uniform hemisphere partitioning scheme. An estimate of the rate of convergence is obtained using the stratified Monte Carlo method. This domain partitioning allows easy parallel realization and leads to convergence improvement of the Monte Carlo method. The high performance and Grid computing of the corresponding Monte Carlo scheme are discussed.

A Nyström method for a class of integral equations on the real line with applications to scattering by diffraction gratings and rough surfaces

We propose a Nystr¨om/product integration method for a class of second kind integral equations on the real line which arise in problems of two-dimensional scalar and elastic wave scattering by unbounded surfaces. Stability and convergence of the method is established with convergence rates dependent on the smoothness of components of the kernel. The method is applied to the problem of acoustic scattering by a sound soft one-dimensional surface which is the graph of a function f, and superalgebraic convergence is established in the case when f is infinitely smooth. Numerical results are presented illustrating this behavior for the case when f is periodic (the diffraction grating case). The Nystr¨om method for this problem is stable and convergent uniformly with respect to the period of the grating, in contrast to standard integral equation methods for diffraction gratings which fail at a countable set of grating periods.

Negative-dimensional integration revisited

Feynman diagrams are the best tool we have to study perturbative quantum field theory. For this very reason the development of any new technique that allows us to compute Feynman integrals is welcome. By the middle of the 1980s, Halliday and Ricotta suggested the possibility of using negative-dimensional integrals to tackle the problem. The aim of this work is to revisit the technique as such and check on its possibilities. For this purpose, we take a box diagram integral contributing to the photon-photon scattering amplitude in quantum electrodynamics using the negative-dimensional integration method. Our approach enables us to quickly reproduce the known results as well as six other solutions as yet unknown in the literature. These six new solutions arise quite naturally in the context of negative-dimensional integration method, revealing a promising technique to handle Feynman integrals.

Integration Methods Used in Numerical Simulations of Electromagnetic Transients

This paper made an analysis of some numerical integration methods that can be used in electromagnetic transient simulations. Among the existing methods, we analyzed the trapezoidal integration method (or Heun formula), Simpson's Rule and Runge-Kutta. These methods were used in simulations of electromagnetic transients in power systems, resulting from switching operations and maneuvers that occur in transmission lines. Analyzed the characteristics such as accuracy, computation time and robustness of the methods of integration.

Field-theory Two-point Functions via Negative Dimension Integration

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Negative dimensional integration: Lab testing at two loops

Negative dimensional integration method (NDIM) is a technique to deal with D-dimensional Feynman loop integrals. Since most of the physical quantities in perturbative Quantum Field Theory (pQFT) require the ability of solving them, the quicker and easier the method to evaluate them the better. The NDIM is a novel and promising technique, ipso facto requiring that we put it to test in different contexts and situations and compare the results it yields with those that we already know by other well-established methods. It is in this perspective that we consider here the calculation of an on-shell two-loop three point function in a massless theory. Surprisingly this approach provides twelve non-trivial results in terms of double power series. More astonishing than this is the fact that we can show these twelve solutions to be different representations for the same well-known single result obtained via other methods. It really comes to us as a surprise that the solution for the particular integral we are dealing with is twelvefold degenerate.

Probing negative-dimensional integration: Two-loop covariant vertex and one-loop light-cone integrals

The negative-dimensional integration method (NDIM) seems to be a very promising technique for evaluating massless and/or massive Feynman diagrams. It is unique in the sense that the method gives solutions in different regions of external momenta simultaneously. Moreover, it is a technique whereby the difficulties associated with performing parametric integrals in the standard approach are transferred to a simpler solving of a system of linear algebraic equations, thanks to the polynomial character of the relevant integrands. We employ this method to evaluate a scalar integral for a massless two-loop three-point vertex with all the external legs off-shell, and consider several special cases for it, yielding results, even for distinct simpler diagrams. We also consider the possibility of NDIM in non-covariant gauges such as the light-cone gauge and do some illustrative calculations, showing that for one-degree violation of covariance (i.e. one external, gauge-breaking, light-like vector n μ) the ensuing results are concordant with the ones obtained via either the usual dimensional regularization technique, or the use of the principal value prescription for the gauge-dependent pole, while for two-degree violation of covariance - i.e. two external, light-like vectors n μ, the gauge-breaking one, and (its dual) n * μ - the ensuing results are concordant with the ones obtained via causal constraints or the use of the so-called generalized Mandelstam-Leibbrandt prescription. © 1999 Elsevier Science B.V.

Feynman integrals with tensorial structure in the negative dimensional integration scheme

The negative-dimensional integration method (NDIM) is revealing itself as a very useful technique for computing massless and/or massive Feynman integrals, covariant and noncovanant alike. Up until now however, the illustrative calculations done using such method have been mostly covariant scalar integrals/without numerator factors. We show here how those integrals with tensorial structures also can be handled straightforwardly and easily. However, contrary to the absence of significant features in the usual approach, here the NDIM also allows us to come across surprising unsuspected bonuses. Toward this end, we present two alternative ways of working out the integrals and illustrate them by taking the easiest Feynman integrals in this category that emerge in the computation of a standard one-loop self-energy diagram. One of the novel and heretofore unsuspected bonuses is that there are degeneracies in the way one can express the final result for the referred Feynman integral.

Simulation of dynamic pinion course using Runge-Kutta's method and impact modeling

Once defined the relationship between the Starter Motor components and their functions, it is possible to develop a mathematical model capable to predict the Starter behavior during operation. One important aspect is the engagement system behavior. The development of a mathematical tool capable of predicting it is a valuable step in order to reduce the design time, cost and engineering efforts. A mathematical model, represented by differential equations, can be developed using physics laws, evaluating force balance and energy flow through the systems degrees of freedom. Another important physical aspect to be considered in this modeling is the impact conditions (particularly on the pinion and ring-gear contact). This work is a report of those equations application on available mathematical software and the resolution of those equations by Runge-Kutta's numerical integration method, in order to build an accessible engineering tool. Copyright © 2011 SAE International.

A new method for the estimation of the reverberation time from measured room impulse responses with application to Italian opera houses

The Schroeder's backward integration method is the most used method to extract the decay curve of an acoustic impulse response and to calculate the reverberation time from this curve. In the literature the limits and the possible improvements of this method are widely discussed. In this work a new method is proposed for the evaluation of the energy decay curve. The new method has been implemented in a Matlab toolbox. Its performance has been tested versus the most accredited literature method. The values of EDT and reverberation time extracted from the energy decay curves calculated with both methods have been compared in terms of the values themselves and in terms of their statistical representativeness. The main case study consists of nine Italian historical theatres in which acoustical measurements were performed. The comparison of the two extraction methods has also been applied to a critical case, i.e. the structural impulse responses of some building elements. The comparison underlines that both methods return a comparable value of the T30. Decreasing the range of evaluation, they reveal increasing differences; in particular, the main differences are in the first part of the decay, where the EDT is evaluated. This is a consequence of the fact that the new method returns a “locally" defined energy decay curve, whereas the Schroeder's method accumulates energy from the tail to the beginning of the impulse response. Another characteristic of the new method for the energy decay extraction curve is its independence on the background noise estimation. Finally, a statistical analysis is performed on the T30 and EDT values calculated from the impulse responses measurements in the Italian historical theatres. The aim of this evaluation is to know whether a subset of measurements could be considered representative for a complete characterization of these opera houses.