Helmholtz Equation Research Articles

Direct Electromagnetic Wave Scattering Calculation Using Methods of Moments through Layered Rough Surface

This thesis focuses on the direct calculation of electromagnetic wave scattering through layered rough surfaces using the Method of Moments. The study aims to contribute to existing knowledge by addressing the problems associated with accurately predicting scattered electromagnetic fields from rough surfaces with layered structures. By using the Method of Moments, this research seeks to provide insights into the behavior of electromagnetic waves when interacting with multi-layered rough surfaces, offering valuable implications for practical engineering applications and the development of advanced computational tools for electromagnetic wave analysis. In practical terms, the surface is typically divided into small patches, and the electromagnetic field at each patch is determined using MoM. This involves solving integral equations describing the interaction between the induced currents on the surface and the incident field. The layered structure may require additional equations to model the reflection of waves and the transmission at the boundaries between the layers. The direct problem is also examined in a two-dimensional picture, where the rough surface profile and the medium are known parameters, and the electromagnetic waves scattered from the surface to the entire space are the unknowns. The rough surface is considered to have a finite length and can be deterministic or arise from a stationary stochastic random process, with a focus on Gaussian-type rough surfaces. The roughness is considered as two basic types: as a perfect electric conductor (PEC) in free space and as a boundary between two dielectric media. The equations of integration obtained are based on the scalar solution of Green's theory, Helmholtz equations, and using Hankel type green functions as kernels due to the one-dimensional nature. The MoM is extensively used to solve these integral equations and is compared to asymptotic approaches and analytical. Various simulations are conducted to test the feasibility of the algorithms, and the results are thoroughly discussed. The study of electromagnetic wave scattering in one-dimensional rough surfaces involves the direct problem, where the field distribution in the whole space is determined using the equation of integration. These equations take into account the surface roughness and medium parameters they are solved using the Method of Moments (MoM) by expressing the field distribution on the surface using basis functions. The direct calculation of electromagnetic wave scattering through a layered rough surface using the Method of Moments (MoM) is a computationally intensive process and complex and. It contains simulating the interaction of EM waves with a complex surface, which presents challenges due to the multiple layers with different properties, and a roughly surface causing and diffraction of incident waves and scattering. When working with a layered rough surface increases the electromagnetic wave interaction complexity. Additional multiple layers introduce considerations, such as layer interfaces and accounting for different material abilities. This requires careful modeling and analysis to exactly capture the behavior of electromagnetic waves as they traverse through and interact with layers. MoM, a commonly used numerical method for solving electromagnetic scattering problems, discreteness the surface into small segments and solves integral equations to determine the fields that is scattered. This approach allows for a maximum analysis of electromagnetic interactions at a grained level.

Enhanced Bolometric Detection of THz Signals by a Resonant Structure for Inclined Radiation Incidence

In this work, we consider the possibility of enhancing terahertz bolometric detection efficiency using resonant structures in the case of an inclined incidence of radiation. The structures are made of a sequence of doped and undoped semiconductors, including epsilon-near-zero areas. Undoped regions act as electromagnetic resonators, thus ensuring resonant signal penetration through the opaque (doped) regions of the structure. A set of epsilon-near-zero areas can ensure substantial enhancements to the electric field in the material. In the doped regions, absorption occurs. The structure described above can provide efficient resonant energy absorption for a wide range of angles of incidence. The numerical calculations based on the solution of the Helmholtz equation have shown that the studied resonant structures ensure the absorption of up to 50% of the incident radiation energy for a 60-degree incidence.

Time-independent, paraxial and time-dependent Madelung trajectories near zeros

The Madelung trajectories associated with a wavefunction are the integral curves (streamlines) of its phase gradient, interpretable in terms of the local velocity (momentum) vector field. The pattern of trajectories provides an immediately visualisable representation of the wave. The patterns can be completely different when the same wave equation describes different physical contexts. For the time-independent Schrödinger or Helmholtz equation, trajectories circulate around the phase singularities (zeros) of the wavefunction; and in the paraxially approximate wave, streamlines spiral slowly in or out of the zeros as well as circulating. But if the paraxial wave equation is reinterpreted as the time-dependent Schrödinger equation, its Madelung trajectories do not circulate around the zeros in spacetime: they undulate while avoiding them, except for isolated trajectories that encounter each zero in a cusp singularity. The different local trajectory geometries are illustrated with two examples; a local model explains the spacetime cusps.

Computational Accuracy and Efficiency of Room Acoustics Simulation Using a Frequency Domain FEM with Air Absorption: 2D Study

Recently, to simulate sound propagation inside architectural spaces at high frequencies, the application of computationally expensive wave-based numerical methods to room acoustics simulation is increasing gradually. Generally, standard room acoustics simulations in the frequency domain are performed based on the lossless Helmholtz equation. However, for acoustics simulation at high frequencies, consideration of the sound attenuation effect caused by air absorption is an aspect to increase the reliability of predictions. Although a simple approach based on the lossy Helmholtz equation is available to include the air absorption effect in the frequency domain, its accuracy and efficiency are still not discussed well. This paper presents an accuracy and efficiency estimation of FEM based on the lossy Helmholtz equation via two numerical problems in two dimensions: a plane wave propagation problem up to 20 kHz in a long duct with 1 km length and a sound propagation problem in a real-scale office up to 6 kHz. Results revealed that the lossy Helmholtz equation-based FEM can include the air absorption effect accurately up to 20 kHz. Moreover, a possibility of providing a higher computational efficiency at higher frequencies is suggested when the magnitude of the pure-tone sound attenuation coefficient is large.

Three-Dimensional Forward Modeling of Transient Electromagnetic Method Considering Induced Polarization Effect Based on Spectral Element Method

The transient electromagnetic method (TEM) is widely used in the exploration of mineral, petroleum, and geothermal resources due to its sensitivity to low-resistivity bodies, limited site constraints, and strong resistance to interference. In practical applications, the TEM often uses a long wire source instead of an idealized horizontal electric dipole (HED) source as the excitation source. This is due to the complex external conditions and the relatively large distance between the receiving zone and the transmitter source. Compared to the HED, the long wire source can provide a larger excitation current, generating stronger signals to meet the requirements of a higher signal-to-noise ratio or deeper exploration. It also produces longer-duration signals, thereby providing better resolution. Additionally, for the interpretation of TEM data, three-dimensional forward modeling plays a crucial role. However, the mature traditional TEM forward method is based on a simple, sometimes inappropriate model, as it is well established that the induced polarization (IP) effect is widely present in the deep earth, especially in oil and gas reservoirs. The presence of the IP effect results in negative responses in field data that do not conform to the traditional theoretical decay law of TEM, which can significantly impact data processing and inversion results. To address this issue, a TEM forward modeling method considering the IP effect based on the spectral element method (SEM) has been developed in this study. Firstly, starting from the Helmholtz equation satisfied by the time domain electric field, we introduce the Debye model with polarization information into the forward modeling by utilizing the differential form of Ohm’s law. As a result, we derive the boundary value problem for the time domain electric field that considers the induced polarization effect. Using Gauss–Lobatto–Legendre (GLL) polynomials as the basis functions, the SEM is employed to discretize the governing equations at each time step and obtain spectral element discretization equations. Then, temporal discretization equations are derived using the second-order backward Euler formula, and the linear system of equations is solved using the Pardiso direct solver. Finally, the electromagnetic responses at any time channel are obtained via SEM interpolation and numerical integration, thereby achieving three-dimensional TEM forward modeling considering the IP effect. The results indicate that this method can effectively reflect the spatial distribution of polarizable subsurface media. It provides valuable references for studying the polarization parameters of subsurface media and performing a three-dimensional inversion of TEM data considering the induced polarization effect.

A highly accurate indirect boundary integral equation solution for three dimensional elastic scattering problem

In this paper, we propose a numerical formula for three dimensional elastic wave scattering problem. Different from the classical boundary integral equation method, the three dimensional elastic scattering problem is transformed into a coupled boundary value problem including a scalar Helmholtz equation and a vector Helmholtz equation by the decomposition of the displacement. The coupled boundary value problem can be solved by the boundary integral equations based on the Helmholtz equation. A denseness result is given to support this algorithm. A comparative study with the classical method of fundamental solutions getting from the boundary integral equation method is given to check the accuracy with different types of boundary conditions. From the numerical results, the present method can give more accurate results than the classical method of fundamental solutions method.

Methodological derivation of the eikonal equation

Usually, when working with the eikonal equation, reference is made to its derivation in the monograph by Born and Wolf. The derivation of this equation was done rather carelessly. Understanding this derivation requires a certain number of implicit assumptions. For a better understanding of the eikonal approximation and for methodological purposes, the authors decided to repeat the derivation of the eikonal equation, explicating all possible assumptions. Methodically, the following algorithm for deriving the eikonal equation is proposed. The wave equation is derived from Maxwell’s equation. In this case, all conditions are explicitly introduced under which it is possible to do this. Further, from the wave equation, the transition to the Helmholtz equation is carried out. From the Helmholtz equation, with the application of certain assumptions, a transition is made to the eikonal equation. After analyzing all the assumptions and steps, the transition from the Maxwell’s equations to the eikonal equation is actually implemented. When deriving the eikonal equation, several formalisms are used. The standard formalism of vector analysis is used as the first formalism. Maxwell’s equations and the eikonal equation are written as three-dimensional vectors. After that, both the Maxwell’s equations and the eikonal equation use the covariant 4-dimensional formalism. The result of the work is a methodically consistent description of the eikonal equation.

On the Fredholm Property and Solvability of a System of Integral Equations in the Transmission Problem for the Helmholtz Equation

A scalar three-dimensional boundary value problem of wave diffraction for the Helmholtz equation with transmission conditions that assume the presence of an infinitely thin material at the media interface is considered. Uniqueness and existence theorems for solutions are proved. The original problem is reduced to a system of integral equations over the media interface. Calculation formulas for the system of linear algebraic equations obtained after applying the collocation method and numerical results for solving the problem when the domain is a ball with certain transmission conditions are given.

Method of virtual sources using on-surface radiation conditions for the Helmholtz equation

We develop a novel method of virtual sources to formulate boundary integral equations for exterior wave propagation problems. However, by contrast to classical boundary integral formulations, we displace the singularity of the Green’s function by a small distance h>0. As a result, the discretization can be performed on the actual physical boundary with continuous kernels so that any naive quadrature scheme can be used to approximate integral operators. Using on-surface radiation conditions, we combine single- and double-layer potential representations of the solution to arrive at a well-conditioned system upon discretization. The virtual displacement parameter h controls the conditioning of the discrete system. We provide mathematical guidance to choose h, in terms of the wavelength and mesh refinements, in order to strike a balance between accuracy and stability. Proof-of-concept implementations are presented, including piecewise linear and isogeometric element formulations in two- and three-dimensional settings. We observe exceptionally well-behaved spectra, and solve the corresponding systems using matrix-free GMRES iterations. The method is implemented for problems with known analytical solutions for comparison. We conclude that the proposed method leads to accurate solutions and good stability for a wide range of wavelengths and mesh refinements.

Estimation of the stress intensity factor of an interface crack by the scattered SH-far field

The diffraction of the plane elastic SH-wave from a semi-infinite interface crack on junction of two elastic semi-infinite media is studied. The crack is modelled by the mathematical cut with no stress on its faces. The displacement and the stress fields are continuous outside of the crack. The wave diffraction problem is reduced to the solution of the mixed boundary value problem for Helmholtz equation. We search the solution which satisfies the Neumann boundary condition on the crack faces and the continuity condition for the stress and displacement fields outside of the crack. The radiation condition at the infinity and the Meixners condition at the crack tip must be also satisfied. Using the Fourier integral transform, the mixed boundary value problem is reduced to the Wiener-Hopf functional equation which is valid in the given strip of regularity in the complex plane. The method of factorization and decomposition, as well as the Liouvilles theorem, are used to solve this equation. Its kernel function is factorized and represented as a product of two split functions that are regular in the overlapping half-planes. These ones allow for the simple poles outside of the regularity regions. The solution of the Wiener-Hopf equation is presented in analytical form. The scattered displacement field is found for an arbitrary frequency and sounding angle by applying an inverse Fourier integral transform to the solution. The asympŹtotic formula for the stress intensity factor (SIF) at the crack tip is obtained. The correlation between complex amplitude of the SH-wave far scattered displacement field and the SIF caused by this field is obtained for an arbitrary radiation angle, frequency and medium parameters. This allows us to express the SIF through the amplitude of the far field if the radiation frequency, sounding angle as well as the physical characteristics of materials are known. It is shown that in the plane that is normal to the tip of the crack the ratio of SIF for the given junction and two fixed values of the sounding angle or two sounding frequencies is proportional to the ratio of the scattered fields. It is found that under above mentioned condition the proportionality rate does not depend on material properties. The obtained relations can be applied for estimation of the SIF with changing frequency and sounding angle.

Recursive analytical quadrature for the close evaluation of Laplace and Helmholtz layer potentials over flat boundary elements in [formula omitted

A method for the analytical evaluation of singular and nearly singular layer potentials arising in the collocation boundary element method for the Laplace and Helmholtz equation is developed for flat boundary elements with polynomial shape functions. The method is based on dimension-reduction via the divergence theorem and a Recursive scheme for evaluating the resulting line Integrals for Polynomial Elements (RIPE). It is used to evaluate single layer, double layer, adjoint double layer, and hypersingular potentials, for both the Laplace and the Helmholtz kernels. It naturally supports nearly singular, singular, and hypersingular integrals under a single framework. The developed recursive algorithm allows accurate evaluation of layer potentials associated with O(p2) density functions used in a pth order boundary element in O(p3) time for the Laplace case.

Application of Adomian Decomposition Method for Solving Helmholtz Equation

The Adomian decomposition method (ADM) is effective for solving linear and nonlinear ordinary and partial differential equations, integral equations, algebraic equations and integro-differential equations and it leads to an analytical solution in the form of an infinite power series. We will solve the Helmholtz equation in two dimensional by the Adomian decomposition method.
 Key words: Application of Adomian Decomposition Method, Helmholtz Equation, Adomian Decomposition Method. 

Far field broadband approximate cloaking for the Helmholtz equation with a Drude-Lorentz refractive index

This paper concerns the analysis of a passive, broadband approximate cloaking scheme for the Helmholtz equation in Rd for d=2 or d=3. Using ideas from transformation optics, we construct an approximate cloak by “blowing up” a small ball of radius ϵ>0 to one of radius 1. In the anisotropic cloaking layer resulting from the “blow-up” change of variables, we incorporate a Drude-Lorentz-type model for the index of refraction, and we assume that the cloaked object is a soft (perfectly conducting) obstacle. We first show that (for any fixed ϵ) there are no real transmission eigenvalues associated with the inhomogeneity representing the cloak, which implies that the cloaking devices we have created will not yield perfect cloaking at any frequency, even for a single incident time harmonic wave. Secondly, we establish estimates on the scattered field due to an arbitrary time harmonic incident wave. These estimates show that, as ϵ approaches 0, the L2-norm of the scattered field outside the cloak, and its far field pattern, approach 0 uniformly over any bounded band of frequencies. In other words: our scheme leads to broadband approximate cloaking for arbitrary incident time harmonic waves.

A Modified Residual Power Series Method for the Approximate Solution of Two-Dimensional Fractional Helmholtz Equations

In this paper, we suggest a modification for the residual power series method that is used to solve fractional-order Helmholtz equations, which is called the Shehu-transform residual power series method (ST-RPSM). This scheme uses a combination of the Shehu transform (ST) and the residual power series method (RPSM). The fractional derivatives are taken with respect to Caputo order. The novelty of this approach is that it does not restrict the fractional order and reduces the need for heavy computational work. The results were obtained using an iterative series that led to an exact solution. The 3D graphical plots for different values of fractional orders are shown to compare ST-RPSM results with exact solutions.

Revising the Boundary Element Method for Thermoviscous Acoustics: An Iterative Approach via Schur Complement

The Helmholtz equation is a reliable model for acoustics in inviscid fluids. Real fluids, however, experience viscous and thermal dissipation that impact the sound propagation dynamics. The viscothermal losses primarily arise in the boundary region between the fluid and solid, the acoustic boundary layers. To preserve model accuracy for structures housing acoustic cavities of comparable size to the boundary layer thickness, meticulous consideration of these losses is essential. Recent research efforts aim to integrate viscothermal effects into acoustic boundary element methods (BEM). While the reduced discretization of BEM is advantageous over finite element methods, it results in fully populated system matrices whose conditioning deteriorates when extended with additional degrees of freedom to account for viscothermal dissipation. Solving such a linear system of equations becomes prohibitively expensive for large-scale applications, as only direct solvers can be used. This work proposes a revised formulation for the viscothermal BEM employing the Schur complement and a change of basis for the boundary coupling. We demonstrate that static condensation significantly improves the conditioning of the coupled problem. When paired with an iterative solution scheme, the approach lowers the algorithmic complexity and thus reduces the computational costs in terms of runtime and storage requirements. The results demonstrate the favorable performance of the new method, indicating its usability for applications of practical relevance in thermoviscous acoustics.

Simulation of filamentation dynamics of microwave discharge in nitrogen

The present study deals with numerical simulations of microwave (MW) discharges in nitrogen based on extended fluid-dynamic model. The set of governing equations for non-equilibrium gas-discharge plasma includes conservation equations for species number densities, electron energy density, Poisson equation for the electric field coupled to the multi-temperature Navier–Stokes fluid-dynamic equations taking into account thermal nonequilibrium; the power transmitted from MW radiation to electrons is determined from the Helmholtz equation. The kinetic scheme includes 61 reactions involving neutral molecules and atoms in the ground and electronically excited states, ions and electrons. The set of equations is solved for a two-dimensional problem under conditions of experiments at a pressure of 40 and 50 Torr and different electromagnetic wave frequencies and pulse duration. The dynamics of discharge formation and transition from the diffuse to the filament form is studied. The results are compared with experimental data, and a good agreement is shown for the time larger than 10 µs. The possible reasons for discrepancies at a shorter time are discussed and the effect of small oxygen impurities on the quantitative characteristics of the discharge are evaluated. The presence of a small oxygen impurity and seed electrons in the region of discharge formation yields a better agreement between numerical and experimental data.

A review of the application of spectral methods in computational ocean acoustics

Spectral methods are a new and niche numerical discretization method with the main characteristic of high accuracy. Spectral methods have been applied in many fields of engineering numerical simulation. Since their introduction into ocean acoustics in 1993, spectral methods have made significant progress in computational ocean acoustics. This article systematically introduces the basic principles of spectral methods, their applicable conditions, and their applications and developments in the normal mode model, the wavenumber integration model, the parabolic equation model, and acoustic Helmholtz equation simulations. At the same time, this article points out the shortcomings of the current application of spectral methods in computational ocean acoustics and potential research directions in the future. The aim is to provide a comprehensive research foundation for subsequent researchers and to promote the application of spectral methods in computational ocean acoustics to go further and deeper.

Hyperparameter Tuning of Neural Network for High-Dimensional Problems in the Case of Helmholtz Equation

Hyperparameter Tuning of Neural Network for High-Dimensional Problems in the Case of Helmholtz Equation

Level set-based topology optimization for programmable acoustic structures

Controlling acoustic waves is essential in noise-reduction applications. Acoustic structures, which can efficiently manipulate acoustic waves based on their structural design, have been used to achieve the desired wave propagation behaviors. To meet various demands on sound propagation, it is crucial to design acoustic structures with multiple functionalities. In this research, we propose a level set-based topology optimization method for acoustic structures with desired acoustic performances that depend on their deformation. Firstly, we introduce a mathematical model that represents acoustic wave propagation based on the Helmholtz equation and the coordinate transformation. Next, we formulate the optimization problem based on the framework of a level set-based topology optimization method. A two-dimensional numerical example is provided to demonstrate the validity of the proposed method. As an example, we optimize a sound reflective structure with the desired acoustic responses that are dependent on the deformation of the optimized structure.

On the use of a tailored fluid-fluid Green's function to predict scattering from two-phase fluid interfaces

In naval applications, the precise knowledge of the acoustical behaviour of two-phase fluids is of interest. Indeed, they appear in many configurations as bubbles curtains, cavitation (along propellers and pipes) or two-phase turbulent boundary layers and wakes (gas exhaust). To model this problem, we focus on a Boundary Element/ Boundary Element coupling for low to moderate frequencies. We propose two approaches: (a) solving the Helmholtz equations in each phase by introducing the free field Green's function and (b) determining a tailored Green's function, taking into account the presence of the two-phase fluid. The determination of a tailored Green's function has two main advantages: (i) it allows a reduction of the numerical model, since it contains all the information on the fluid-fluid coupling, in particular the transmission conditions and resonances and (ii) it is source-independent and thus it gives directly the answer to any source. Numerical tests on the evaluation of radiated noise are performed to determine the efficiency of the approach based on a tailored Green's function compared to the free-field based one.