The transient response of an anisotropic rectangular plate subjected to impact is described through a Chebyshev collocation multidomain discretization of the Reissner-Mindlin plate equations. The trapezoidal rule is used for time-integration. The spatial collocation derivative operators are represented by matrices, and the subdomains are patched by natural and essential conditions. At each time level the resulting governing matrix equation is reduced by two consecutive block Gaussian eliminations, so that an equation for the variables at the subdomain corners has to be solved. Back-substitution gives the variables at all other collocation points. The time history as represented by computed contour plots has been compared with analytical results and with photos produced by holographic interferometry. The agreements are satisfactory

Here are considered time-harmonic electromagnetic waves in a quadratic waveguide consisting of a periodic dielectric core enclosed by conducting walls. The permittivity function may be smooth or have jumps. The electromagnetic field is given by a magnetic vector potential in Lorenz gauge, and defined on a Floquet cell. The Helmholtz operator is approximated by a Chebyshev collocation, Fourier–Galerkin method. Laurent's rule and the inverse rule are employed for the representation of Fourier coefficients of products of functions. The computations yield, for known wavenumbers, values of the first few eigenfrequencies of the field. In general, the dispersion curves exhibit band gaps. Field patterns are identified as transverse electric, TE, transverse magnetic, TM, or hybrid modes. Maxwell's equations are fulfilled. A few trivial solutions appear when the permittivity varies in the guiding direction and across it. The results of the present method are consistent with exact results and with those obtained by a low-order finite element software. The present method is more efficient than the low-order finite element method

The transient response of a tube subjected to impact is described through Fourier-Galerkin and Chebyshev collocation multidomain discretizations of the equations of linear elasticity. The trapezoidal rule is used for time integration. For each Fourier mode the spatial collocation derivative operators are represented by matrices, and the subdomains are patched by natural and essential conditions. At each time level the resulting governing matrix equation is reduced by two consecutive block Gaussian eliminations, so that an equation for the complex Fourier coefficients at die subdomain corners has to be solved. Back-substitution gives the coefficients at all other collocation points. An inverse discrete Fourier transform generates, at optional time levels, the three components of the displacement field. Through this method the long-term evolution of the field may be calculated, provided the impact time is long enough. The time history as represented by computed contour plots has been compared with photos produced by holographic interferometry. The agreements are satisfactory.

The behaviour is compared of two solvers for the discrete equations arising from the discretization using Chebyshev collocation of a second-order linear partial differential equation on a square. The alternative solvers considered are a direct solver and an iterative solver based on preconditioning with the matrix arising from finite-difference discretization of the governing equation. The total error of the collocation derivatives and the separate contributions from round-off and discretization error are examined. The efficiency of the two solvers is compared. The iterative solver is more efficient than the direct solver on fine grids for equations similar to the Poisson equation, provided that there are Dirichlet boundary conditions on at least three of the sides of the square.

Here are described four solvers for time-harmonic electromagnetic fields in checkerboard patterns. A pattern is built by four squares with constant permittivity, inline image or inline image. It is enclosed by conducting walls or is a unit cell of a periodic structure. The field is represented in two ways: by inline image, the transverse component of the magnetic induction, and by inline image, the magnetic vector potential in Lorenz gauge. inline image and inline image satisfy Helmholtz equations in each square as well as transmission and boundary conditions (BCs). These governing equations yield eigensolutions inline image and inline image, which are found to be inline image at worst. Variational versions of the governing equations are introduced. The weak formulations for inline image are standard, while those for inline image are new. They imply that the derivative transmission and BCs are satisfied weakly on interfaces between regions with different permittivity. Eigenpairs are computed approximately by spectral element methods. They yield mutually consistent eigenpairs. However, only about half of the eigenpairs (inline image) correspond to eigenpairs (inline image). For each set of BCs, the first few eigenfrequencies inline image are given by tables, and some of the eigenfunctions are presented by contour plots.

This paper presents a viscous compressible flow problem to which an equilibrium solution, in terms of density and velocity, can be given implicitly by elementary functions. The corresponding initial boundary value problem is solved by time discretization by the Crank-Nicolson method, Newton linearization and space discretization using multidomain Chebyshev collocation techniques. The physical interval is covered by subintervals of equal length. Each subinterval utilizes the same number of collocation points and each interface consists of one or two points. Six ways of patching are tested. All of them yield solutions with spectral accuracy for a few time steps, but only three are stable in the long run. Details of the density evolution are illustrated

The electromagnetic initial-boundary value problem for a cavity enclosed by perfectly conducting walls is considered. The cavity medium is defined by its permittivity and permeability which vary continuously in space. The electromagnetic field comes from a source in the cavity. The field is described by a magnetic vector potential A satisfying a wave equation with initial-boundary conditions. This description through A is rigorously shown to give a unique solution of the problem and is the starting point for numerical computations. A Chebyshev collocation solver has been implemented for a cubic cavity, and it has been compared to a standard finite element solver. The results obtained are consistent while the collocation solver performs substantially faster. Some time histories and spectra are computed.