The aim of the work is to study the electromagnetic field theory, finite element method and the interaction of electromagnetic field with tissues. Gained knowledge is then used to calculate spreading of the electromagnetic field in the microwave field and to create a temperature profile of spreading the electromagnetic fields in human tissue. The finite element method was implemented in the Matlab programming environment, where the 1D model was created in the frequency and time domain and a simple 2D model created in time domain. The program was developed to analyze spreading electromagnetic wave. Another part of work was done in the programming environment of COMSOL Multiphysics. In this case was the human leg exposed to electromagnetic fields. The analysis determined the changes of temperature in these biological tissues for six minutes.

This thesis concerns the applications and analysis of the Heterogeneous Multiscale methods (HMM) for Multiscale Elliptic and Hyperbolic Partial Differential Equations. We have gathered the main contributions in two papers. The first paper deals with the cell-boundary error which is present in multi-scale algorithms for elliptic homogenization problems. Typical multi-scale methods have two essential components: a macro and a micro model. The micro model is used to upscale parameter values which are missing in the macro model. Solving the micro model requires, on the other hand, imposing boundary conditions on the boundary of the microscopic domain. Imposing a naive boundary condition leads to $O(\varepsilon/\eta)$ error in the computation, where $\varepsilon$ is the size of the microscopic variations in the media and $\eta$ is the size of the micro-domain. Until now, strategies were proposed to improve the convergence rate up to fourth-order in $\varepsilon/\eta$ at best. However, the removal of this error in multi-scale algorithms still remains an important open problem. In this paper, we present an approach with a time-dependent model which is general in terms of dimension. With this approach we are able to obtain $O((\varepsilon/\eta)^q)$ and $O((\varepsilon/\eta)^q  + \eta^p)$ convergence rates in periodic and locally-periodic media respectively, where $p,q$ can be chosen arbitrarily large.      In the second paper, we analyze a multi-scale method developed under the Heterogeneous Multi-Scale Methods (HMM) framework for numerical approximation of wave propagation problems in periodic media. In particular, we are interested in the long time $O(\varepsilon^{-2})$ wave propagation. In the method, the microscopic model uses the macro solutions as initial data. In short-time wave propagation problems a linear interpolant of the macro variables can be used as the initial data for the micro-model. However, in long-time multi-scale wave problems the linear data does not suffice and one has to use a third-degree interpolant of the coarse data to capture the $O(1)$ dispersive effects apperaing in the long time. In this paper, we prove that through using an initial data consistent with the current macro state, HMM captures this dispersive effects up to any desired order of accuracy in terms of $\varepsilon/\eta$. We use two new ideas, namely quasi-polynomial solutions of periodic problems and local time averages of solutions of periodic hyperbolic PDEs. As a byproduct, these ideas naturally reveal the role of consistency for high accuracy approximation of homogenized quantities. / <p>QC 20130926</p>

Dans cette thèse, on se propose de participer à la popularisation des méthodes de résolution de problèmes de propagation d'onde basées sur des formulations intégrales en fournissant des indicateurs d'erreur a posteriori utilisable dans le cadre d'algorithmes de raffinement autoadaptatif. Le développement de tels indicateurs est complexe du fait de la non-localité des normes associées aux espaces de Sobolev et des opérateurs entrant en jeu. Des indicateurs de la littérature sont étendus au cas de la propagation d'une onde acoustique. On étend les preuves de convergence quasi-optimale (de la littérature) des algorithmes autoadaptatifs associés dans ce cas. On propose alors une nouvelle approche par rapport à la littérature qui consiste à utiliser une technique de localisation des normes, non pas basée sur des inégalités inverses, mais sur l'utilisation d'un opérateur Λ de localisation bien choisi.On peut alors construire des indicateurs d'erreur a posteriori fiables, efficaces, locaux et asymptotiquement exacts par rapport à la norme de Galerkin de l'erreur. On donne ensuite une méthode pour la construction de tels indicateurs. Les applications numériques sur des géométries 2D et 3D confirment l'exactitude asymptotique ainsi que l'optimalité du guidage de l'algorithme autoadaptatif.On étend ensuite ces indicateurs au cas de la propagation d'une onde électromagnétique. Plus précisément, on s'intéresse au cas de l'EFIE. On propose des généralisations des indicateurs de la littérature. On effectue la preuve de convergence quasi-optimale dans le cas d'un indicateur basé sur une localisation de la norme du résidu. On utilise le principe du Λ pour obtenir le premier indicateur d'erreur fiable, efficace et local pour cette équation. On en propose une seconde forme qui est également, théoriquement asymptotiquement exacte. / The aim of this work is to participate to the popularization of methods for the resolution of wave propagation problems based on integral equations formulations by developping a posteriori error estimates in the context of autoadaptive mesh refinement strategies. The development of such estimates is difficult because of the non-locality of the norms associated to the Sobolev spaces and of the involved integral operators. Estimates from the literature are extended in the case of the propagation of an acoustic wave. The proofs of quasi-optimal convergence of the autoadaptive algorithms are established. We then introduce a new approach with respect to the literature which is based on a new norm-localization technique based on the use of a well-chosen Λ operator and not on inverse inequalities as it was the case previously.We then establish new a posteriori error estimates which are reliable, efficient, local and asymptotically exact with respect to the Galerkin norm of the error. We give a method for the construction of such estimates. Numerical applications on 2D and 3D geometries confirm the asymptotic exactness and the optimality of the autoadaptive algorithm.These estimates are extended in the case of the propagation of an electromagnetic wave. More precisely, we are interested in the EFIE. We suggest generalization of the estimates of the literature. A proof for quasi-optimal convergence is given for an estimate based on a localization of the norm of the residual. The principle of Λ is used to construct the first reliable, efficient, local error estimate for this equation. We give a second forme which is eventually theoretically asymptotically exact.

No description available.

Low frequency acoustic room behaviour can be approximated using numerical methods. Traditionally, music studio control rooms are built with complex geometries, making their eigenmodes difficult to predict mathematically. Hence, a summation-by-parts method with simultaneous-approximation-terms is derived to approximate the time dependent acoustic wave equation in three dimensions. The derived model is limited to rectangular prismatic rooms but planned to be expanded to handle complex geometries in the future. Semi-reflecting boundary conditions are used, corresponding to tabulated reflection and absorption properties of real. walls. Two speakers are modeled as omnidirectional point sources placed on a boundary, to mimic common studio setups. Through tests and examination of eigenvalues of the matrices in the method, conditions for stability and reflection coefficients are derived. Simulations of sound pressure distribution produced by the model correlate well to room mode theory, suggesting the model to be accurate in the application of predicting low frequency acoustic room behaviour. However, the convergence rate of the model turns out to be lower than expected when point sources are introduced. Future steps towards applying the model to real music studio control rooms include modeling the walls as changes in density and wave speed rather than boundaries of the domain. This would potentially allow more complex geometries to be modeled within a larger, rectangular domain.

京都大学 / 新制・課程博士 / 博士(農学) / 甲第23527号 / 農博第2474号 / 新制||農||1087(附属図書館) / 学位論文||R3||N5358(農学部図書室) / 京都大学大学院農学研究科地域環境科学専攻 / (主査)教授 藤原 正幸, 教授 中村 公人, 准教授 宇波 耕一 / 学位規則第4条第1項該当 / Doctor of Agricultural Science / Kyoto University / DGAM

A stochastic method is developed, implemented and investigated here for solving Laplace, Poisson's, and standard parabolic wave equations. This method is based on the properties of random walk, diffusion process, Ito formula, Dynkin formula and Monte Carlo simulations. The developed method is a local method i:e: it gives the value of the solution directly at an arbitrary point rather than extracting its value from complete field solution and thus is inherently parallel. Field computation by this method is demonstrated for electrostatic and electrodynamic propagation problems by considering simple examples and numerical results are presented to validate this method. Numerical investigations are carried out to understand efficacy and limitations of this method and to provide qualitative understanding of various parameters involved in this method.

In this work the scattering of an incident plane wave propagating along a plane perpendicular to the xy-plane is studied. The wave is scattered from an infinitely long cylindrical object of arbitrary cross-section. Due to the arbitrary geometry of the obstacle, a finite differences numerical method is employed to approximate the solution of the scattering problems. The wave equation is expressed in terms of generalized curvilinear coordinates. Boundary conforming grids are generated using elliptic grid generators. Then, a explicit marching in time scheme is implemented over these grids. It is found that as time grows the numerical solution converges to a wave with harmonic time dependence. The amplitude of these waves is analyzed and graphed over generalized grids for different geometries. An important physical measure of the energy scattered, the differential scattering cross section, is also obtained. In particular, the method is applied to a circular cylindrical obstacle. For this case, the analytical solution can also be obtained by traditional spectral techniques. The method is validated by comparing this exact solution with the numerical approximation obtained from the application of it.

Wind power is a rapidly increasing resource of electrical power world-wide. With the increasing number of wind turbines installed one major concern is the noise they generate. Sometimes already built wind turbines have to be put down or down-regulated, when certain noise levels are exceeded, resulting in economical and environmental losses. Therefore, accurate sound propagation calculations would be beneficial already in a planning stage of a wind farm. A model that can account for varying wind speeds and complex terrains could therefore be of great importance when future wind farms are planned. In this report an extended version of the classical wave equation that allows for variations in wind speed and terrain is derived which can be used to solve complex terrain and wind settings. The equation are solved with the use of Fourier transforms and Chebyshev polynomials and a numerical code is developed. The numerical code is evaluated against test cases where analytical and simple solutions exist. Tests with no wind for both totally free propagation and with a ground surface is evaluated in both 2D and 3D settings. For these simple cases the developed code shows good agreement to analytical solutions if the computational domain is sufficiently large. More advanced test cases with wind and terrain is not evaluated in this report and needs further validation. If the sound pressure needs to be calculated for a large area, and if the frequency is high, the developed model has problems regarding computational time and memory. These problems could be solved by further development of the numerical code or by using other solution methods.

No description available.