Controlled time integration for the numerical simulation of meteor radar reflections

We model meteoroids entering the Earth[U+05F3]s atmosphere as objects surrounded by non-magnetized plasma, and consider efficient numerical simulation of radar reflections from meteors in the time domain. Instead of the widely used finite difference time domain method (FDTD), we use more generalized finite differences by applying the discrete exterior calculus (DEC) and non-uniform leapfrog-style time discretization. The computational domain is presented by convex polyhedral elements. The convergence of the time integration is accelerated by the exact controllability method. The numerical experiments show that our code is efficiently parallelized. The DEC approach is compared to the volume …

On a numerical solution of the Maxwell equations by discrete exterior calculus

Generalized finite difference schemes with higher order Whitney forms

Finite difference kind of schemes are popular in approximating wave propagation problems in finite dimensional spaces. While Yee’s original paper on the finite difference method is already from the sixties, mathematically there still remains questions which are not yet satisfactorily covered. In this paper, we address two issues of this kind. Firstly, in the literature Yee’s scheme is constructed separately for each particular type of wave problem. Here, we explicitly generalize the Yee scheme to a class of wave problems that covers at large physics field theories. For this we introduce Yee’s scheme for all problems of a class characterised on a Minkowski manifold by (i) a pair of first ord…

Using Wave Propagation Simulations and Convolutional Neural Networks to Retrieve Thin Film Thickness from Hyperspectral Images

Ill-posed inversion problems are one of the major challenges when there is a need to combine measurements with the theory and numerical model. In this study, we demonstrate the use of wave propagation simulations to train a convolutional neural network (CNN) for retrieving sub-wavelength thickness profiles of thin film coatings from hyperspectral images. The simulations are produced by solving numerically one-dimensional wave equation with a method based on Discrete Exterior Calculus (DEC). This approach provides a powerful tool to produce large sets of training data for the neural network. CNN was verified by simulated verification sets and measured reflectance spectra, both of which showe…

Ray optics for absorbing particles with application to ice crystals at near-infrared wavelengths

Abstract Light scattering by particles large compared to the wavelength of incident light is traditionally solved using ray optics which considers absorption inside the particle approximately, along the ray paths. To study the effects rising from this simplification, we have updated the ray-optics code SIRIS to take into account the propagation of light as inhomogeneous plane waves inside an absorbing particle. We investigate the impact of this correction on traditional ray-optics computations in the example case of light scattering by ice crystals through the extended near-infrared (NIR) wavelength regime. In this spectral range, ice changes from nearly transparent to opaque, and therefore…

Systematisation of Systems Solving Physics Boundary Value Problems

A general conservation law that defines a class of physical field theories is constructed. First, the notion of a general field is introduced as a formal sum of differential forms on a Minkowski manifold. By the action principle the conservation law is defined for such a general field. By construction, particular field notions of physics, e.g., magnetic flux, electric field strength, stress, strain etc. become instances of the general field. Hence, the differential equations that constitute physical field theories become also instances of the general conservation law. The general field and the general conservation law together correspond to a large class of relativistic hyperbolic physical …

Three-dimensional splitting dynamics of giant vortices in Bose-Einstein condensates

We study the splitting dynamics of giant vortices in dilute Bose-Einstein condensates by numerically integrating the three-dimensional Gross-Pitaevskii equation in time. By taking advantage of tetrahedral tiling in the spatial discretization, we decrease the error and increase the reliability of the numerical method. An extensive survey of vortex splitting symmetries is presented for different aspect ratios of the harmonic trapping potential. The symmetries of the splitting patterns observed in the simulated dynamics are found to be in good agreement with predictions obtained by solving the dominant dynamical instabilities from the corresponding Bogoliubov equations. Furthermore, we observe…

Generalized wave propagation problems and discrete exterior calculus

We introduce a general class of second-order boundary value problems unifying application areas such as acoustics, electromagnetism, elastodynamics, quantum mechanics, and so on, into a single framework. This also enables us to solve wave propagation problems very efficiently with a single software system. The solution method precisely follows the conservation laws in finite-dimensional systems, whereas the constitutive relations are imposed approximately. We employ discrete exterior calculus for the spatial discretization, use natural crystal structures for three-dimensional meshing, and derive a “discrete Hodge” adapted to harmonic wave. The numerical experiments indicate that the cumulat…

Efficient Time Integration of Maxwell's Equations with Generalized Finite Differences

We consider the computationally efficient time integration of Maxwell’s equations using discrete exterior calculus (DEC) as the computational framework. With the theory of DEC, we associate the degrees of freedom of the electric and magnetic fields with primal and dual mesh structures, respectively. We concentrate on mesh constructions that imitate the geometry of the close packing in crystal lattices that is typical of elemental metals and intermetallic compounds. This class of computational grids has not been used previously in electromagnetics. For the simulation of wave propagation driven by time-harmonic source terms, we provide an optimized Hodge operator and a novel time discretizati…

How much is enough? : The convergence of finite sample scattering properties to those of infinite media

We study the scattering properties of a cloud of particles. The particles are spherical, close to the incident wavelength in size, have a high albedo, and are randomly packed to 20% volume density. We show, using both numerically exact methods for solving the Maxwell equations and radiative-transfer-approximation methods, that the scattering properties of the cloud converge after about ten million particles in the system. After that, the backward-scattered properties of the system should estimate the properties of a macroscopic, practically infinite system. (C) 2021 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license (http://creativecommons.o…

High-quality discretizations for microwave simulations

We apply high-quality discretizations to simulate electromagnetic microwaves. Instead of the vector field presentations, we focus on differential forms and discretize the model in the spatial domain using the discrete exterior calculus. At the discrete level, both the Hodge operators and the time discretization are optimized for time-harmonic simulations. Non-uniform spatial and temporal discretization are applied in problems in which the wavelength is highly-variable and geometry contains sub-wavelength structures. peerReviewed