Experimentelle und numerische Untersuchungen zur Dispersion von Dichteströmungen in einem stochastischen Modellaquifer

Es ist bekannt, dass die Dichte eines gelösten Stoffes die Richtung und die Stärke seiner Bewegung im Untergrund entscheidend bestimmen kann. Eine Vielzahl von Untersuchungen hat gezeigt, dass die Verteilung der Durchlässigkeiten eines porösen Mediums diese Dichteffekte verstärken oder abmindern kann. Wie sich dieser gekoppelte Effekt auf die Vermischung zweier Fluide auswirkt, wurde in dieser Arbeit untersucht und dabei das experimentelle sowohl mit dem numerischen als auch mit dem analytischen Modell gekoppelt. Die auf der Störungstheorie basierende stochastische Theorie der macrodispersion wurde in dieser Arbeit für den Fall der transversalen Makodispersion. Für den Fall einer stabilen Schichtung wurde in einem Modelltank (10m x 1.2m x 0.1m) der Universität Kassel eine Serie sorgfältig kontrollierter zweidimensionaler Experimente an einem stochastisch heterogenen Modellaquifer durchgeführt. Es wurden Versuchsreihen mit variierenden Konzentrationsdifferenzen (250 ppm bis 100 000 ppm) und Strömungsgeschwindigkeiten (u = 1 m/ d bis 8 m/d) an drei verschieden anisotrop gepackten porösen Medien mit variierender Varianzen und Korrelationen der lognormal verteilten Permeabilitäten durchgeführt. Die stationäre räumliche Konzentrationsausbreitung der sich ausbreitenden Salzwasserfahne wurde anhand der Leitfähigkeit gemessen und aus der Höhendifferenz des 84- und 16-prozentigen relativen Konzentrationsdurchgang die Dispersion berechnet. Parallel dazu wurde ein numerisches Modell mit dem dichteabhängigen Finite-Elemente-Strömungs- und Transport-Programm SUTRA aufgestellt. Mit dem kalibrierten numerischen Modell wurden Prognosen für mögliche Transportszenarien, Sensitivitätsanalysen und stochastische Simulationen nach der Monte-Carlo-Methode durchgeführt. Die Einstellung der Strömungsgeschwindigkeit erfolgte - sowohl im experimentellen als auch im numerischen Modell - über konstante Druckränder an den Ein- und Auslauftanks. Dabei zeigte sich eine starke Sensitivität der räumlichen Konzentrationsausbreitung hinsichtlich lokaler Druckvariationen. Die Untersuchungen ergaben, dass sich die Konzentrationsfahne mit steigendem Abstand von der Einströmkante wellenförmig einem effektiven Wert annähert, aus dem die Makrodispersivität ermittelt werden kann. Dabei zeigten sich sichtbare nichtergodische Effekte, d.h. starke Abweichungen in den zweiten räumlichen Momenten der Konzentrationsverteilung der deterministischen Experimente von den Erwartungswerten aus der stochastischen Theorie. Die transversale Makrodispersivität stieg proportional zur Varianz und Korrelation der lognormalen Permeabilitätsverteilung und umgekehrt proportional zur Strömungsgeschwindigkeit und Dichtedifferenz zweier Fluide. Aus dem von Welty et al. [2003] mittels Störungstheorie entwickelten dichteabhängigen Makrodispersionstensor konnte in dieser Arbeit die stochastische Formel für die transversale Makrodispersion weiter entwickelt und - sowohl experimentell als auch numerisch - verifiziert werden.

Is the Helmholtz equation really sign-indefinite?

The usual variational (or weak) formulations of the Helmholtz equation are sign-indefinite in the sense that the bilinear forms cannot be bounded below by a positive multiple of the appropriate norm squared. This is often for a good reason, since in bounded domains under certain boundary conditions the solution of the Helmholtz equation is not unique at wavenumbers that correspond to eigenvalues of the Laplacian, and thus the variational problem cannot be sign-definite. However, even in cases where the solution is unique for all wavenumbers, the standard variational formulations of the Helmholtz equation are still indefinite when the wavenumber is large. This indefiniteness has implications for both the analysis and the practical implementation of finite element methods. In this paper we introduce new sign-definite (also called coercive or elliptic) formulations of the Helmholtz equation posed in either the interior of a star-shaped domain with impedance boundary conditions, or the exterior of a star-shaped domain with Dirichlet boundary conditions. Like the standard variational formulations, these new formulations arise just by multiplying the Helmholtz equation by particular test functions and integrating by parts.

A GRADIENT-LIKE NONAUTONOMOUS EVOLUTION PROCESS

In this paper we consider a dissipative damped wave equation with nonautonomous damping of the form u(tt) + beta(t)u(t) - Delta u + f(u) (1) in a bounded smooth domain Omega subset of R(n) with Dirichlet boundary conditions, where f is a dissipative smooth nonlinearity and the damping beta : R -> (0, infinity) is a suitable function. We prove, if (1) has finitely many equilibria, that all global bounded solutions of (1) are backwards and forwards asymptotic to equilibria. Thus, we give a class of examples of nonautonomous evolution processes for which the structure of the pullback attractors is well understood. That complements the results of [Carvalho & Langa, 2009] on characterization of attractors, where it was shown that a small nonautonomous perturbation of an autonomous gradient-like evolution process is also gradient-like. Note that the evolution process associated to (1) is not a small nonautonomous perturbation of any autonomous gradient-like evolution processes. Moreover, we are also able to prove that the pullback attractor for (1) is also a forwards attractor and that the rate of attraction is exponential.

Continuity of attractors for parabolic problems with localized large diffusion

In this paper we study the continuity of asymptotics of semilinear parabolic problems of the form u(t) - div(p(x)del u) + lambda u =f(u) in a bounded smooth domain ohm subset of R `` with Dirichlet boundary conditions when the diffusion coefficient p becomes large in a subregion ohm(0) which is interior to the physical domain ohm. We prove, under suitable assumptions, that the family of attractors behave upper and lower semicontinuously as the diffusion blows up in ohm(0). (c) 2006 Elsevier Ltd. All rights reserved.

The two-loop massless (lambda/4!)phi(4) model in nontranslational invariant domain

We study the (lambda/4!)phi(4) massless scalar field theory in a four-dimensional Euclidean space, where all but one of the coordinates are unbounded. We are considering Dirichlet boundary conditions in two hyperplanes, breaking the translation invariance of the system. We show how to implement the perturbative renormalization up to two-loop level of the theory. First, analyzing the full two and four-point functions at the one-loop level, we show that the bulk counterterms are sufficient to render the theory finite. Meanwhile, at the two-loop level, we must also introduce surface counterterms in the bare Lagrangian in order to make finite the full two and also four-point Schwinger functions. (c) 2006 American Institute of Physics.

Effects of varying curvature and width on the electronic states of GaAs quantum rings

Stationary states of an electron in thin GaAs elliptical quantum rings are calculated within the effective-mass approximation. The width of the ring varies smoothly along the centerline, which is an ellipse. The solutions of the Schrödinger equation with Dirichlet boundary conditions are approximated by a product of longitudinal and transversal wave functions. The ground-state probability density shows peaks: (i) where the curvature is larger in a constant-with ring, and (ii) in thicker parts of a circular ring. For rings of typical dimensions, it is shown that the effects of a varying width may be stronger than those of the varying curvature. Also, a width profile which compensates the main localization effects of the varying curvature is obtained.

Interpolating EFGM for computing continuous and discontinuous electromagnetic fields

Purpose - This paper proposes an interpolating approach of the element-free Galerkin method (EFGM) coupled with a modified truncation scheme for solving Poisson's boundary value problems in domains involving material non-homogeneities. The suitability and efficiency of the proposed implementation are evaluated for a given set of test cases of electrostatic field in domains involving different material interfaces.Design/methodology/approach - the authors combined an interpolating approximation with a modified domain truncation scheme, which avoids additional techniques for enforcing the Dirichlet boundary conditions and for dealing with material interfaces usually employed in meshfree formulations.Findings - the local electric potential and field distributions were correctly described as well as the global quantities like the total potency and resistance. Since, the treatment of the material interfaces becomes practically the same for both the finite element method (FEM) and the proposed EFGM, FEM-oriented programs can, thus, be easily extended to provide EFGM approximations.Research limitations/implications - the robustness of the proposed formulation became evident from the error analyses of the local and global variables, including in the case of high-material discontinuity.Practical implications - the proposed approach has shown to be as robust as linear FEM. Thus, it becomes an attractive alternative, also because it avoids the use of additional techniques to deal with boundary/interface conditions commonly employed in meshfree formulations.Originality/value - This paper reintroduces the domain truncation in the EFGM context, but by using a set of interpolating shape functions the authors avoided the use of Lagrange multipliers as well Mathematics in Engineering high-material discontinuity.

Quantum rings of arbitrary shape and non-uniform width in a threading magnetic field

The electronic states of quantum rings with centerlines of arbitrary shape and non-uniform width in a threading magnetic field are calculated. The solutions of the Schrodinger equation with Dirichlet boundary conditions are obtained by a variational separation of variables in curvilinear coordinates. We obtain a width profile that compensates for the main effects of the curvature variations in the centerline. Numerical results are shown for circular, elliptical, and limacon-shaped quantum rings. We also show that smooth and tiny variations in the width may strongly affect the Aharonov-Bohm oscillations.

Aplicação da transformada integral e da transformação conforme na solução de uma classe de problemas difusivo-convectivos em domínios de geometrias não-convencionais

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

Sistemas ecológicos modelados por equações de reação-difusão

Pós-graduação em Física - IFT

Densidade de energia e força de radiação sobre fronteiras em movimento para estados arbitrários do campo: aplicação ao estado coerente

No presente trabalho foi considerado um campo escalar real não massivo em um espaço-tempo bidimensional, satisfazendo à condição de fronteira Dirichlet ou Neumann na posição instantânea de uma fronteira em movimento. Para uma lei de movimento relativística, foi mostrado que as condiçõoes de fronteira Dirichlet e Neumann produzem a mesma força de radiação sobre um espelho em movimento quando o estado inicial do campo é invariante sobre translações temporais. Obtemos as fórmulas exatas para a densidade de energia do campo e da força de radiação na fronteira para os estados de vácuo, coerente e comprimido. No limite não-relativistico, os resultados obtidos coincidem com os encontrados na literatura. Também foi investigado o campo dentro de uma cavidade oscilante. Considerando as condiçõoes de fronteira Neumann e Dirichlet, escreveu-se a fórmulas exata para a densidade de energia dentro de uma cavidade não-estática, para um estado inicial arbitrário do campo. Tomando como base a equação de Moore, nós calculamos recursivamente a densidade de energia e investigamos a evolução temporal da densidade de energia para o estado coerente do campo.

Stability analysis of Crank-Nicolson and Euler schemes for time-dependent diffusion equations

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

Isospectrality and Heat Content

We present examples of isospectral operators that do not have the same heat content. Several of these examples are planar polygons that are isospectral for the Laplace operator with Dirichlet boundary conditions. These include examples with infinitely many components. Other planar examples have mixed Dirichlet and Neumann boundary conditions. We also consider Schrodinger operators acting in L-2[0,1] with Dirichlet boundary conditions, and show that an abundance of isospectral deformations do not preserve the heat content.

hp-dGFEM for Second-Order Elliptic Problems in Plyhedra II: Exponential Convergence

The goal of this paper is to establish exponential convergence of $hp$-version interior penalty (IP) discontinuous Galerkin (dG) finite element methods for the numerical approximation of linear second-order elliptic boundary-value problems with homogeneous Dirichlet boundary conditions and piecewise analytic data in three-dimensional polyhedral domains. More precisely, we shall analyze the convergence of the $hp$-IP dG methods considered in [D. Schötzau, C. Schwab, T. P. Wihler, SIAM J. Numer. Anal., 51 (2013), pp. 1610--1633] based on axiparallel $\sigma$-geometric anisotropic meshes and $\bm{s}$-linear anisotropic polynomial degree distributions.

Simulación numérica computacional de fenómenos dinámicos

En este proyecto se trata la simulación numérica de un fenómeno dinámico, basado en el comportamiento de una onda transmitida a lo largo de una cuerda elástica de un instrumento musical, cuyos extremos se encuentran anclados. El fenómeno físico, se desarrolla utilizando una ecuación en derivadas parciales hiperbólicas con variables espacial y temporal, acompañada por unas condiciones de contorno tipo Dirichlet en los extremos y por más condiciones iniciales que dan comienzo al proceso. Posteriormente se han generado algoritmos para el método numérico empleado (Diferencias finitas centrales y progresivas) y la programación del problema aproximado con su consistencia, estabilidad y convergencia, obteniéndose unos resultados acordes con la solución analítica del problema matemático. La programación y salida de resultados se ha realizado con Visual Studio 8.0. y la programación de objetos con Visual Basic .Net In this project the topic is the numerical simulation of a dynamic phenomenon, based on the behavior of a transmitted wave along an elastic string of a musical instrument, whose ends are anchored. The physical phenomenon is developed using a hyperbolic partial differential equation with spatial and temporal variables, accompanied by a Dirichlet boundary conditions at the ends and more initial conditions that start the process. Subsequently generated algorithms for the numerical method used (central and forward finite differences) and the programming of the approximate problem with consistency, stability and convergence, yielding results in line with the analytical solution of the mathematical problem. Programming and output results has been made with Visual Studio 8.0. and object programming with Visual Basic. Net