Spelling suggestions: "subject:"helmholtzequation"" "subject:"helmholtzekvation""
81 |
Stabilized finite element methods for convection-diffusion-reaction, helmholtz and stokes problemsNadukandi, Prashanth 13 May 2011 (has links)
We present three new stabilized finite element (FE) based Petrov-Galerkin methods for the convection-diffusionreaction
(CDR), the Helmholtz and the Stokes problems, respectively. The work embarks upon a priori analysis of a
consistency recovery procedure for some stabilization methods belonging to the Petrov- Galerkin framework. It was
ound that the use of some standard practices (e.g. M-Matrices theory) for the design of essentially non-oscillatory
numerical methods is not appropriate when consistency recovery methods are employed. Hence, with respect to
convective stabilization, such recovery methods are not preferred. Next, we present the design of a high-resolution
Petrov-Galerkin (HRPG) method for the CDR problem. The structure of the method in 1 D is identical to the consistent
approximate upwind (CAU) Petrov-Galerkin method [doi: 10.1016/0045-7825(88)90108-9] except for the definitions of
he stabilization parameters. Such a structure may also be attained via the Finite Calculus (FIC) procedure [doi:
10.1 016/S0045-7825(97)00119-9] by an appropriate definition of the characteristic length. The prefix high-resolution is
used here in the sense popularized by Harten, i.e. second order accuracy for smooth/regular regimes and good
shock-capturing in non-regular re9jmes. The design procedure in 1 D embarks on the problem of circumventing the
Gibbs phenomenon observed in L projections. Next, we study the conditions on the stabilization parameters to
ircumvent the global oscillations due to the convective term. A conjuncture of the two results is made to deal with the
problem at hand that is usually plagued by Gibbs, global and dispersive oscillations in the numerical solution. A multi
dimensional extension of the HRPG method using multi-linear block finite elements is also presented.
Next, we propose a higher-order compact scheme (involving two parameters) on structured meshes for the Helmholtz
equation. Making the parameters equal, we recover the alpha-interpolation of the Galerkin finite element method
(FEM) and the classical central finite difference method. In 1 D this scheme is identical to the alpha-interpolation
method [doi: 10.1 016/0771 -050X(82)90002-X] and in 2D choosing the value 0.5 for both the parameters, we recover
he generalized fourth-order compact Pade approximation [doi: 10.1 006/jcph.1995.1134, doi: 10.1016/S0045-
7825(98)00023-1] (therein using the parameter V = 2). We follow [doi: 10.1 016/0045-7825(95)00890-X] for the
analysis of this scheme and its performance on square meshes is compared with that of the quasi-stabilized FEM [doi:
10.1016/0045-7825(95)00890-X]. Generic expressions for the parameters are given that guarantees a dispersion
accuracy of sixth-order should the parameters be distinct and fourth-order should they be equal. In the later case, an
expression for the parameter is given that minimizes the maximum relative phase error in 2D. A Petrov-Galerkin
ormulation that yields the aforesaid scheme on structured meshes is also presented. Convergence studies of the
error in the L2 norm, the H1 semi-norm and the I ~ Euclidean norm is done and the pollution effect is found to be small. / Presentamos tres nuevos metodos estabilizados de tipo Petrov- Galerkin basado en elementos finitos (FE) para los
problemas de convecci6n-difusi6n- reacci6n (CDR), de Helmholtz y de Stokes, respectivamente. El trabajo comienza
con un analisis a priori de un metodo de recuperaci6n de la consistencia de algunos metodos de estabilizaci6n que
pertenecen al marco de Petrov-Galerkin. Hallamos que el uso de algunas de las practicas estandar (por ejemplo, la
eoria de Matriz-M) para el diserio de metodos numericos esencialmente no oscilatorios no es apropiado cuando
utilizamos los metodos de recu eraci6n de la consistencia. Por 10 tanto, con res ecto a la estabilizaci6n de
conveccion, no preferimos tales metodos de recuperacion . A continuacion, presentamos el diser'io de un metodo de
Petrov-Galerkin de alta-resolucion (HRPG) para el problema CDR. La estructura del metodo en 10 es identico al
metodo CAU [doi: 10.1016/0045-7825(88)90108-9] excepto en la definicion de los parametros de estabilizacion. Esta
estructura tambien se puede obtener a traves de la formulacion del calculo finito (FIC) [doi: 10.1 016/S0045-
7825(97)00119-9] usando una definicion adecuada de la longitud caracteristica. El prefijo de "alta-resolucion" se
utiliza aqui en el sentido popularizado por Harten, es decir, tener una solucion con una precision de segundo orden
en los regimenes suaves y ser esencialmente no oscilatoria en los regimenes no regulares. El diser'io en 10 se
embarca en el problema de eludir el fenomeno de Gibbs observado en las proyecciones de tipo L2. A continuacion,
estudiamos las condiciones de los parametros de estabilizacion para evitar las oscilaciones globales debido al
ermino convectivo. Combinamos los dos resultados (una conjetura) para tratar el problema COR, cuya solucion
numerica sufre de oscilaciones numericas del tipo global, Gibbs y dispersiva. Tambien presentamos una extension
multidimensional del metodo HRPG utilizando los elementos finitos multi-lineales.
fa. continuacion, proponemos un esquema compacto de orden superior (que incluye dos parametros) en mallas
estructuradas para la ecuacion de Helmholtz. Haciendo igual ambos parametros, se recupera la interpolacion lineal
del metodo de elementos finitos (FEM) de tipo Galerkin y el clasico metodo de diferencias finitas centradas. En 10
este esquema es identico al metodo AIM [doi: 10.1 016/0771 -050X(82)90002-X] y en 20 eligiendo el valor de 0,5 para
ambos parametros, se recupera el esquema compacto de cuarto orden de Pade generalizada en [doi:
10.1 006/jcph.1 995.1134, doi: 10.1 016/S0045-7825(98)00023-1] (con el parametro V = 2). Seguimos [doi:
10.1 016/0045-7825(95)00890-X] para el analisis de este esquema y comparamos su rendimiento en las mallas
uniformes con el de "FEM cuasi-estabilizado" (QSFEM) [doi: 10.1016/0045-7825 (95) 00890-X]. Presentamos
expresiones genericas de los para metros que garantiza una precision dispersiva de sexto orden si ambos parametros
son distintos y de cuarto orden en caso de ser iguales. En este ultimo caso, presentamos la expresion del parametro
que minimiza el error maxima de fase relativa en 20. Tambien proponemos una formulacion de tipo Petrov-Galerkin
~ue recupera los esquemas antes mencionados en mallas estructuradas. Presentamos estudios de convergencia del
error en la norma de tipo L2, la semi-norma de tipo H1 y la norma Euclidiana tipo I~ y mostramos que la perdida de
estabilidad del operador de Helmholtz ("pollution effect") es incluso pequer'ia para grandes numeros de onda.
Por ultimo, presentamos una coleccion de metodos FE estabilizado para el problema de Stokes desarrollados a
raves del metodo FIC de primer orden y de segundo orden. Mostramos que varios metodos FE de estabilizacion
existentes y conocidos como el metodo de penalizacion, el metodo de Galerkin de minimos cuadrados (GLS) [doi:
10.1016/0045-7825(86)90025-3], el metodo PGP (estabilizado a traves de la proyeccion del gradiente de presion)
[doi: 10.1 016/S0045-7825(96)01154-1] Y el metodo OSS (estabilizado a traves de las sub-escalas ortogonales) [doi:
10.1016/S0045-7825(00)00254-1] se recuperan del marco general de FIC. Oesarrollamos una nueva familia de
metodos FE, en adelante denominado como PLS (estabilizado a traves del Laplaciano de presion) con las formas no
lineales y consistentes de los parametros de estabilizacion. Una caracteristica distintiva de la familia de los metodos
PLS es que son no lineales y basados en el residuo, es decir, los terminos de estabilizacion dependera de los
residuos discretos del momento y/o las ecuaciones de incompresibilidad. Oiscutimos las ventajas y desventajas de
estas tecnicas de estabilizaci6n y presentamos varios ejemplos de aplicacion
|
82 |
Fast algorithms for frequency domain wave propagationTsuji, Paul Hikaru 22 February 2013 (has links)
High-frequency wave phenomena is observed in many physical settings, most notably in acoustics, electromagnetics, and elasticity. In all of these fields, numerical simulation and modeling of the forward propagation problem is important to the design and analysis of many systems; a few examples which rely on these computations are the development of
metamaterial technologies and geophysical prospecting for natural resources. There are two modes of modeling the forward problem: the frequency domain and the time domain. As the title states, this work is concerned with the former regime.
The difficulties of solving the high-frequency wave propagation problem accurately lies in the large number of degrees of freedom required. Conventional wisdom in the computational electromagnetics commmunity suggests that about 10 degrees of freedom per wavelength be used in each coordinate direction to resolve each oscillation. If K is the width of the domain in wavelengths, the number of unknowns N grows at least by O(K^2) for surface discretizations and O(K^3) for volume discretizations in 3D. The memory requirements and asymptotic complexity estimates of direct algorithms such as the multifrontal method are too costly for such problems. Thus, iterative solvers must be used. In this dissertation, I will present fast algorithms which, in conjunction with GMRES, allow the solution of the forward problem in O(N) or O(N log N) time. / text
|
83 |
Stability of finite element solutions to Maxwell's equations in frequency domainSchwarzbach, Christoph 12 October 2009 (has links) (PDF)
Eine Standardformulierung der Randwertaufgabe für die Beschreibung zeitharmonischer elektromagnetischer Phänomene hat die Vektor-Helmholtzgleichung für das elektrische Feld zur Grundlage. Bei niedrigen Frequenzen führt der große Nullraum des Rotationsoperators zu einem instabilen Lösungsverhalten. Wird die Randwertaufgabe zum Beispiel mit Hilfe der Methode der Finiten Elemente in ein lineares Gleichungssystem überführt, äußert sich die Instabilität in einer schlechten Konditionszahl ihrer Koeffizientenmatrix. Eine stabilere Formulierung wird durch die explizite Berücksichtigung der Kontinuitätsgleichung erreicht. Zur numerischen Lösung der Randwertaufgaben wurde eine Finite-Elemente-Software erstellt. Sie berücksichtigt unter anderem unstrukturierte Gitter, räumlich variable, anisotrope Materialparameter sowie die Erweiterung der Maxwell-Gleichungen durch Perfectly Matched Layers. Die Software wurde anhand von Anwendungen in der marinen Geophysik erfolgreich getestet. Insbesondere demonstriert die Einbeziehung von Seebodentopographie in Form einer stetigen Oberflächentriangulierung die geometrische Flexibilität der Software. / The physics of time-harmonic electromagnetic phenomena can be mathematically described by boundary value problems. A standard approach is based on the vector Helmholtz equation in terms of the electric field. The curl operator involved has a large, non-trivial kernel which leads to an instable solution behaviour at low frequencies. If the boundary value problem is solved approximately using, e. g., the
finite element method, the instability expresses itself by a badly conditioned coefficient matrix of the ensuing system of linear equations. A stable formulation is obtained by taking the continuity equation explicitly into account. In order to solve the boundary value problem numerically a finite element software package has been implemented. Its features comprise, amongst others, the treatment of
unstructured meshes and piecewise polynomial, anisotropic constitutive parameters as well as the extension of Maxwell’s equations to the Perfectly Matched Layer. Successful application of the software is demonstrated with examples from marine geophysics. In particular, the incorporation of seafloor topography by a continuous
surface triangulation illustrates the geometric flexibility of the software.
|
84 |
Développement d'outils statistiques pour l'amélioration de dispositifs d'imagerie acoustique et micro-ondeDiong, Mouhamadou 09 December 2015 (has links)
L'un des enjeux majeurs pour les systèmes d'imagerie par diffraction acoustique et micro-onde, est l'amélioration des performances obtenues au moment de la reconstruction des objets étudiés. Cette amélioration peut s'effectuer par la recherche d'algorithmes d'imagerie plus performants d'une part et par la recherche d'une meilleure configuration de mesures d'autre part. La première approche (recherche d'algorithmes) permet d'améliorer le processus d'extraction de l'information présente dans un échantillon de mesures donné. Néanmoins, la qualité des résultats d'imagerie reste limitée par la quantité d'information initialement disponible. La seconde approche consiste à choisir la configuration de mesures de manière à augmenter la quantité d'information disponible dans les données. Pour cette approche, il est nécessaire de quantifier la quantité d'information dans les données. En théorie de l'estimation, ceci équivaut à quantifier la performance du système. Dans cette thèse, nous utilisons la Borne de Cramer Rao comme mesure de performance, car elle permet d'analyser la performance des systèmes de mesures sans être influencé par le choix de la méthode d'inversion utilisée. Deux analyses sont proposées dans ce manuscrit. La première consiste en l'évaluation de l'influence des facteurs expérimentaux sur la performance d'inversion. Cette analyse a été effectuée pour différents objets le tout sous une hypothèse de configuration bidimensionnelle. La seconde analyse consiste à comparer les performances de l'estimateur obtenu avec l'approximation de Born aux valeurs de la borne de Cramer Rao (BCR); l'objectif étant d'illustrer d'autres applications possibles de la BCR. / Improving the performance of diffraction based imaging systems constitutes a major issue in both acoustic and electromagnetic scattering. To solve this problem, two main approaches can be explored. The first one consists in improving the inversion algorithms used in diffraction based imaging. However, while this approach generally leads to a significant improvement of the performance of the imaging system, it remains limited by the initial amount of information available within the measurements. The second one consists in improving the measurement system in order to maximize the amount of information within the experimental data. This approach does require a quantitative mean of measuring the amount of information available. In estimation problems, the {appraisal of the} performance of the system is often used for that purpose. In this Ph.D. thesis, we use the Cramer Rao bound to assess the performance of the imaging system. In fact, this quantity has the advantage of providing an assessment which is independent from the inversion algorithm used. Two main analysis are discussed in this thesis. The first analysis explores the influence on the system's performance, of several experimental conditions such as the antennas positions, the radiation pattern of the source, the properties of the background medium, etc. Two classes of objects are considered: 2D homogeneous circular cylindrical objects and 2D cylindrical objects with defect. The second analysis studies the performance of an estimator based on Born approximation with the Cramer Rao Bound as reference. The aim of this second analysis is to showcase other possible applications for the Cramer Rao Bound.
|
85 |
Stability of finite element solutions to Maxwell's equations in frequency domainSchwarzbach, Christoph 10 August 2009 (has links)
Eine Standardformulierung der Randwertaufgabe für die Beschreibung zeitharmonischer elektromagnetischer Phänomene hat die Vektor-Helmholtzgleichung für das elektrische Feld zur Grundlage. Bei niedrigen Frequenzen führt der große Nullraum des Rotationsoperators zu einem instabilen Lösungsverhalten. Wird die Randwertaufgabe zum Beispiel mit Hilfe der Methode der Finiten Elemente in ein lineares Gleichungssystem überführt, äußert sich die Instabilität in einer schlechten Konditionszahl ihrer Koeffizientenmatrix. Eine stabilere Formulierung wird durch die explizite Berücksichtigung der Kontinuitätsgleichung erreicht. Zur numerischen Lösung der Randwertaufgaben wurde eine Finite-Elemente-Software erstellt. Sie berücksichtigt unter anderem unstrukturierte Gitter, räumlich variable, anisotrope Materialparameter sowie die Erweiterung der Maxwell-Gleichungen durch Perfectly Matched Layers. Die Software wurde anhand von Anwendungen in der marinen Geophysik erfolgreich getestet. Insbesondere demonstriert die Einbeziehung von Seebodentopographie in Form einer stetigen Oberflächentriangulierung die geometrische Flexibilität der Software. / The physics of time-harmonic electromagnetic phenomena can be mathematically described by boundary value problems. A standard approach is based on the vector Helmholtz equation in terms of the electric field. The curl operator involved has a large, non-trivial kernel which leads to an instable solution behaviour at low frequencies. If the boundary value problem is solved approximately using, e. g., the
finite element method, the instability expresses itself by a badly conditioned coefficient matrix of the ensuing system of linear equations. A stable formulation is obtained by taking the continuity equation explicitly into account. In order to solve the boundary value problem numerically a finite element software package has been implemented. Its features comprise, amongst others, the treatment of
unstructured meshes and piecewise polynomial, anisotropic constitutive parameters as well as the extension of Maxwell’s equations to the Perfectly Matched Layer. Successful application of the software is demonstrated with examples from marine geophysics. In particular, the incorporation of seafloor topography by a continuous
surface triangulation illustrates the geometric flexibility of the software.
|
86 |
Méthodes d'ordre élevé et méthodes de décomposition de domaine efficaces pour les équations de Maxwell en régime harmonique / Efficient high order and domain decomposition methods for the time-harmonic Maxwell's equationsBonazzoli, Marcella 11 September 2017 (has links)
Les équations de Maxwell en régime harmonique comportent plusieurs difficultés lorsque la fréquence est élevée. On peut notamment citer le fait que leur formulation variationnelle n’est pas définie positive et l’effet de pollution qui oblige à utiliser des maillages très fins, ce qui rend problématique la construction de solveurs itératifs. Nous proposons une stratégie de solution précise et rapide, qui associe une discrétisation par des éléments finis d’ordre élevé à des préconditionneurs de type décomposition de domaine. La conception, l’implémentation et l’analyse des deux méthodes sont assez difficiles pour les équations de Maxwell. Les éléments finis adaptés à l’approximation du champ électrique sont les éléments finis H(rot)-conformes ou d’arête. Ici nous revisitons les degrés de liberté classiques définis par Nédélec, afin d’obtenir une expression plus pratique par rapport aux fonctions de base d’ordre élevé choisies. De plus, nous proposons une technique pour restaurer la dualité entre les fonctions de base et les degrés de liberté. Nous décrivons explicitement une stratégie d’implémentation qui a été appliquée dans le langage open source FreeFem++. Ensuite, nous nous concentrons sur les techniques de préconditionnement du système linéaire résultant de la discrétisation par éléments finis. Nous commençons par la validation numérique d’un préconditionneur à un niveau, de type Schwarz avec recouvrement, avec des conditions de transmission d’impédance entre les sous-domaines. Enfin, nous étudions comment des préconditionneurs à deux niveaux, analysés récemment pour l’équation de Helmholtz, se comportent pour les équations de Maxwell, des points de vue théorique et numérique. Nous appliquons ces méthodes à un problème à grande échelle qui découle de la modélisation d’un système d’imagerie micro-onde, pour la détection et le suivi des accidents vasculaires cérébraux. La précision et la vitesse de calcul sont essentielles dans cette application. / The time-harmonic Maxwell’s equations present several difficulties when the frequency is large, such as the sign-indefiniteness of the variational formulation, the pollution effect and the problematic construction of iterative solvers. We propose a precise and efficient solution strategy that couples high order finite element (FE) discretizations with domain decomposition (DD) preconditioners. High order FE methods make it possible for a given precision to reduce significantly the number of unknowns of the linear system to be solved. DD methods are then used as preconditioners for the iterative solver: the problem defined on the global domain is decomposed into smaller problems on subdomains, which can be solved concurrently and using robust direct solvers. The design, implementation and analysis of both these methods are particularly challenging for Maxwell’s equations. FEs suited for the approximation of the electric field are the curl-conforming or edge finite elements. Here, we revisit the classical degrees of freedom (dofs) defined by Nédélec to obtain a new more friendly expression in terms of the chosen high order basis functions. Moreover, we propose a general technique to restore duality between dofs and basis functions. We explicitly describe an implementation strategy, which we embedded in the open source language FreeFem++. Then we focus on the preconditioning of the linear system, starting with a numerical validation of a one-level overlapping Schwarz preconditioner, with impedance transmission conditions between subdomains. Finally, we investigate how two-level preconditioners recently analyzed for the Helmholtz equation work in the Maxwell case, both from the theoretical and numerical points of view. We apply these methods to the large scale problem arising from the modeling of a microwave imaging system, for the detection and monitoring of brain strokes. In this application accuracy and computing speed are indeed of paramount importance.
|
87 |
Better imaging for landmine detection : an exploration of 3D full-wave inversion for ground-penetrating radarWatson, Francis Maurice January 2016 (has links)
Humanitarian clearance of minefields is most often carried out by hand, conventionally using a a metal detector and a probe. Detection is a very slow process, as every piece of detected metal must treated as if it were a landmine and carefully probed and excavated, while many of them are not. The process can be safely sped up by use of Ground-Penetrating Radar (GPR) to image the subsurface, to verify metal detection results and safely ignore any objects which could not possibly be a landmine. In this thesis, we explore the possibility of using Full Wave Inversion (FWI) to improve GPR imaging for landmine detection. Posing the imaging task as FWI means solving the large-scale, non-linear and ill-posed optimisation problem of determining the physical parameters of the subsurface (such as electrical permittivity) which would best reproduce the data. This thesis begins by giving an overview of all the mathematical and implementational aspects of FWI, so as to provide an informative text for both mathematicians (perhaps already familiar with other inverse problems) wanting to contribute to the mine detection problem, as well as a wider engineering audience (perhaps already working on GPR or mine detection) interested in the mathematical study of inverse problems and FWI.We present the first numerical 3D FWI results for GPR, and consider only surface measurements from small-scale arrays as these are suitable for our application. The FWI problem requires an accurate forward model to simulate GPR data, for which we use a hybrid finite-element boundary-integral solver utilising first order curl-conforming N\'d\'{e}lec (edge) elements. We present a novel `line search' type algorithm which prioritises inversion of some target parameters in a region of interest (ROI), with the update outside of the area defined implicitly as a function of the target parameters. This is particularly applicable to the mine detection problem, in which we wish to know more about some detected metallic objects, but are not interested in the surrounding medium. We may need to resolve the surrounding area though, in order to account for the target being obscured and multiple scattering in a highly cluttered subsurface. We focus particularly on spatial sensitivity of the inverse problem, using both a singular value decomposition to analyse the Jacobian matrix, as well as an asymptotic expansion involving polarization tensors describing the perturbation of electric field due to small objects. The latter allows us to extend the current theory of sensitivity in for acoustic FWI, based on the Born approximation, to better understand how polarization plays a role in the 3D electromagnetic inverse problem. Based on this asymptotic approximation, we derive a novel approximation to the diagonals of the Hessian matrix which can be used to pre-condition the GPR FWI problem.
|
88 |
Field reconstructions and range tests for acoustics and electromagnetics in homogeneous and layered media / Feld-Rekonstruktionen und Range Tests für Akustik und Elektromagnetik in homogenen und geschichteten MedienSchulz, Jochen 04 December 2007 (has links)
No description available.
|
89 |
Partikelmodellierung der Strukturbildung akustischer Kavitationsblasen in Wechselwirkung mit dem Schalldruckfeld / Particle modeling of acoustic cavitation bubble structure formation and interaction with the acoustic pressure fieldKoch, Philipp 29 August 2006 (has links)
No description available.
|
Page generated in 0.0848 seconds