The effect of the equivalent-weights particle filter on dynamical balance in a primitive equation model

The disadvantage of the majority of data assimilation schemes is the assumption that the conditional probability density function of the state of the system given the observations [posterior probability density function (PDF)] is distributed either locally or globally as a Gaussian. The advantage, however, is that through various different mechanisms they ensure initial conditions that are predominantly in linear balance and therefore spurious gravity wave generation is suppressed. The equivalent-weights particle filter is a data assimilation scheme that allows for a representation of a potentially multimodal posterior PDF. It does this via proposal densities that lead to extra terms being added to the model equations and means the advantage of the traditional data assimilation schemes, in generating predominantly balanced initial conditions, is no longer guaranteed. This paper looks in detail at the impact the equivalent-weights particle filter has on dynamical balance and gravity wave generation in a primitive equation model. The primary conclusions are that (i) provided the model error covariance matrix imposes geostrophic balance, then each additional term required by the equivalent-weights particle filter is also geostrophically balanced; (ii) the relaxation term required to ensure the particles are in the locality of the observations has little effect on gravity waves and actually induces a reduction in gravity wave energy if sufficiently large; and (iii) the equivalent-weights term, which leads to the particles having equivalent significance in the posterior PDF, produces a change in gravity wave energy comparable to the stochastic model error. Thus, the scheme does not produce significant spurious gravity wave energy and so has potential for application in real high-dimensional geophysical applications.

Patch-size and isolation effects in the Fisher-Kolmogorov equation

We examine the classical problem of the existence of a threshold size for a patch to allow for survival of a given population in the case where the patch is not completely isolated. The surrounding habitat matrix is characterized by a non-zero carrying capacity. We show that a critical patch size cannot be strictly defined in this case. We also obtain the saturation density in such a patch as a function of the size of the patch and the relative carrying capacity of the outer region. We argue that this relative carrying capacity is a measure of the isolation of the patch. Our results are then compared with conclusions drawn from observations of the population dynamics of understorey birds in fragments of the Amazonian forest and shown to qualitatively agree with them, offering an explanation for the importance of dispersal and isolation in these observations. Finally, we show that a generalized critical patch size can be introduced resorting to threshold densities for the observation of a given species.

Generalised single-equation load flow method for unbalanced distribution systems

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Renormalization of the ladder light-front Bethe-Salpeter equation in the Yukawa model

The reduction of the two-fermion Bethe-Salpeter equation in the framework of light-front dynamics is studied for the Yukawa model. It yields auxiliary three-dimensional quantities for the transition matrix and the bound state. The arising effective interaction can be perturbatively expanded in powers of the coupling constant gs allowing a defined number of boson exchanges; it is divergent and needs renormalization; it also includes the instantaneous term of the Dirac propagator. One possible solution of the renormalization problem of the boson exchanges is shown to be provided by expanding the effective interaction beyond single boson exchange. The effective interaction in ladder approximation up to order g4 s is discussed in detail. It is shown that the effective interaction naturally yields the box counterterm required to be introduced ad hoc previously. The covariant results of the Bethe-Salpeter equation can be recovered from the corresponding auxiliary three-dimensional quantities.

Fokker-Planck Equation for a Metastable Time Dependent Potential

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

The Fokker-Planck equation for a bistable potential

The Fokker-Planck equation is studied through its relation to a Schrodinger-type equation. The advantage of this combination is that we can construct the probability distribution of the Fokker-Planck equation by using well-known solutions of the Schrodinger equation. By making use of such a combination, we present the solution of the Fokker-Planck equation for a bistable potential related to a double oscillator. Thus, we can observe the temporal evolution of the system describing its dynamic properties such as the time tau to overcome the barrier. By calculating the rates k = 1/tau as a function of the inverse scaled temperature 1/D, where D is the diffusion coefficient, we compare the aspect of the curve k x 1/D, with the ones obtained from other studies related to four different kinds of activated process. We notice that there are similarities in some ranges of the scaled temperatures, where the different processes follow the Arrhenius behavior. We propose that the type of bistable potential used in this study may be used, qualitatively, as a simple model, whose rates share common features with the rates of some single rate-limited thermally activated processes. (C) 2014 Elsevier B.V. All rights reserved.

Three dimensional Lifshitz black hole and the Korteweg-de Vries equation

We consider a solution of three dimensional New Massive Gravity with a negative cosmological constant and use the AdS/CTF correspondence to inquire about the equivalent two dimensional model at the boundary. We conclude that there should be a close relation of the theory with the Korteweg-de Vries equation. (C) 2012 Elsevier B.V..All rights reserved.

Human thermal comfort: an irreversibility-based approach emulating empirical clothed-body correlations and the conceptual energy balance equation

Exergetic analysis can provide useful information as it enables the identification of irreversible phenomena bringing about entropy generation and, therefore, exergy losses (also referred to as irreversibilities). As far as human thermal comfort is concerned, irreversibilities can be evaluated based on parameters related to both the occupant and his surroundings. As an attempt to suggest more insights for the exergetic analysis of thermal comfort, this paper calculates irreversibility rates for a sitting person wearing fairly light clothes and subjected to combinations of ambient air and mean radiant temperatures. The thermodynamic model framework relies on the so-called conceptual energy balance equation together with empirical correlations for invoked thermoregulatory heat transfer rates adapted for a clothed body. Results suggested that a minimum irreversibility rate may exist for particular combinations of the aforesaid surrounding temperatures. By separately considering the contribution of each thermoregulatory mechanism, the total irreversibility rate rendered itself more responsive to either convective or radiative clothing-influenced heat transfers, with exergy losses becoming lower if the body is able to transfer more heat (to the ambient) via convection.

Spectrum reconstruction from a scattering measurement using the adjoint Boltzmann transport equation for photons

Quality control of medical radiological systems is of fundamental importance, and requires efficient methods for accurately determine the X-ray source spectrum. Straightforward measurements of X-ray spectra in standard operating require the limitation of the high photon flux, and therefore the measure has to be performed in a laboratory. However, the optimal quality control requires frequent in situ measurements which can be only performed using a portable system. To reduce the photon flux by 3 magnitude orders an indirect technique based on the scattering of the X-ray source beam by a solid target is used. The measured spectrum presents a lack of information because of transport and detection effects. The solution is then unfolded by solving the matrix equation that represents formally the scattering problem. However, the algebraic system is ill-conditioned and, therefore, it is not possible to obtain a satisfactory solution. Special strategies are necessary to circumvent the ill-conditioning. Numerous attempts have been done to solve this problem by using purely mathematical methods. In this thesis, a more physical point of view is adopted. The proposed method uses both the forward and the adjoint solutions of the Boltzmann transport equation to generate a better conditioned linear algebraic system. The procedure has been tested first on numerical experiments, giving excellent results. Then, the method has been verified with experimental measurements performed at the Operational Unit of Health Physics of the University of Bologna. The reconstructed spectra have been compared with the ones obtained with straightforward measurements, showing very good agreement.

Analytische Ableitungen der Energie für den "Equation-of-Motion-Coupled-Cluster"-Ansatz zur Berechnung von Ionisationspotentialen

In der vorliegenden Arbeit wird die Theorie der analytischen zweiten Ableitungen für die EOMIP-CCSD-Methode formuliert sowie die durchgeführte Implementierung im Quantenchemieprogramm CFOUR beschrieben. Diese Ableitungen sind von Bedeutung bei der Bestimmung statischer Polarisierbarkeiten und harmonischer Schwingungsfrequenzen und in dieser Arbeit wird die Genauigkeit des EOMIP-CCSD-Ansatzes bei der Berechnung dieser Eigenschaften für verschiedene radikalische Systeme untersucht. Des Weiteren können mit Hilfe der ersten und zweiten Ableitungen vibronische Kopplungsparameter berechnet werden, welche zur Simulation von Molekülspektren in Kombination mit dem Köppel-Domcke-Cederbaum (KDC)-Modell - in der Arbeit am Beispiel des Formyloxyl (HCO2)-Radikals demonstriert - benötigt werden.rnrnDer konzeptionell einfache EOMIP-CC-Ansatz wurde gewählt, da hier die Wellenfunktion eines Radikalsystems ausgehend von einem stabilen geschlossenschaligen Zustand durch die Entfernung eines Elektrons gebildet wird und somit die Problematik der Symmetriebrechung umgangen werden kann. Im Rahmen der Implementierung wurden neue Programmteile zur Lösung der erforderlichen Gleichungen für die gestörten EOMIP-CC-Amplituden und die gestörten Lagrange-Multiplikatoren zeta zum Quantenchemieprogramm CFOUR hinzugefügt. Die unter Verwendung des Programms bestimmten Eigenschaften werden hinsichtlich ihrer Leistungsfähigkeit im Vergleich zu etablierten Methoden wie z.B. CCSD(T) untersucht. Bei der Berechnung von Polarisierbarkeiten und harmonischen Schwingungsfrequenzen liefert die EOMIP-CCSD-Theorie meist gute Resultate, welche nur wenig von den CCSD(T)-Ergebnissen abweichen. Einzig bei der Betrachtung von Radikalen, für die die entsprechenden Anionen nicht stabil sind (z.B. NH2⁻ und CH3⁻), liefert der EOMIP-CCSD-Ansatz aufgrund methodischer Nachteile keine aussagekräftige Beschreibung. rnrnDie Ableitungen der EOMIP-CCSD-Energie lassen sich auch zur Simulation vibronischer Kopplungen innerhalb des KDC-Modells einsetzen.rnZur Kopplung verschiedener radikalischer Zustände in einem solchen Modellpotential spielen vor allem die Ableitungen von Übergangsmatrixelementen eine wichtige Rolle. Diese sogenannten Kopplungskonstanten können in der EOMIP-CC-Theorie besonders leicht definiert und berechnet werden. Bei der Betrachtung des Photoelektronenspektrums von HCO2⁻ werden zwei Alternativen untersucht: Die vertikale Bestimmung an der Gleichgewichtsgeometrie des HCO2⁻-Anions und die Ermittlung adiabatischer Kraftkonstanten an den Gleichgewichtsgeometrien des Radikals. Lediglich das adiabatische Modell liefert bei Beschränkung auf harmonische Kraftkonstanten eine qualitativ sinnvolle Beschreibung des Spektrums. Erweitert man beide Modelle um kubische und quartische Kraftkonstanten, so nähern sich diese einander an und ermöglichen eine vollständige Zuordnung des gemessenen Spektrums innerhalb der ersten 1500 cm⁻¹. Die adiabatische Darstellung erreicht dabei nahezu quantitative Genauigkeit.

Virtual-bound, filamentary and layered states in a box-shaped quantum dot of square potential form the exact numerical solution of the effective mass Schrodinger equation

The effective mass Schrodinger equation of a QD of parallelepipedic shape with a square potential well is solved by diagonalizing the exact Hamiltonian matrix developed in a basis of separation-of-variables wavefunctions. The expected below bandgap bound states are found not to differ very much from the former approximate calculations. In addition, the presence of bound states within the conduction band is confirmed. Furthermore, filamentary states bounded in two dimensions and extended in one dimension and layered states with only one dimension bounded, all within the conduction band which are similar to those originated in quantum wires and quantum wells coexist with the ordinary continuum spectrum of plane waves. All these subtleties are absent in spherically shaped quantum dots, often used for modeling.

A new equation for macroscopic description of capillary rise in porous media

Capillary rise in porous media is frequently modeled using the Washburn equation. Recent accurate measurements of advancing fronts clearly illustrate its failure to describe the phenomenon in the long term. The observed underprediction of the position of the front is due to the neglect of dynamic saturation gradients implicit in the formulation of the Washburn equation. We consider an approximate solution of the governing macroscopic equation, which retains these gradients, and derive new analytical formulae for the position of the advancing front, its speed of propagation, and the cumulative uptake. The new solution properly describes the capillary rise in the long term, while the Washburn equation may be recovered as a special case. (C) 2004 Elsevier Inc. All rights reserved.

On an indirect integral equation approach for stationary heat transfer in semi-infinite layered domains in R3 with cavities

Regions containing internal boundaries such as composite materials arise in many applications.We consider a situation of a layered domain in IR3 containing a nite number of bounded cavities. The model is stationary heat transfer given by the Laplace equation with piecewise constant conductivity. The heat ux (a Neumann condition) is imposed on the bottom of the layered region and various boundary conditions are imposed on the cavities. The usual transmission (interface) conditions are satised at the interface layer, that is continuity of the solution and its normal derivative. To eciently calculate the stationary temperature eld in the semi-innite region, we employ a Green's matrix technique and reduce the problem to boundary integral equations (weakly singular) over the bounded surfaces of the cavities. For the numerical solution of these integral equations, we use Wienert's approach [20]. Assuming that each cavity is homeomorphic with the unit sphere, a fully discrete projection method with super-algebraic convergence order is proposed. A proof of an error estimate for the approximation is given as well. Numerical examples are presented that further highlights the eciency and accuracy of the proposed method.

Null Condition for Semilinear Wave Equation with Variable Coefficients

∗The author was partially supported by M.U.R.S.T. Progr. Nazionale “Problemi Non Lineari...”

Fermat's Equation in Matrices

The Fermat equation is solved in integral two by two matrices of determinant one as well as in finite order integral three by three matrices.