Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Maxwell equations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Equations of motion for a free-electron laser with an electromagnetic pump field and an axial electrostatic field

The equations of motion for a free-electron laser (FEL) with an electromagnetic pump field and a static axial electric field are derived using a Hamiltonian formalism. Equations governing the energy transfer between the electron beam and each of the electromagnetic fields are given, and the phase shift for each of the electromagnetic fields is derived from a linearized Maxwell wave equation. The relation between the static axial electric field and the resonant phase is given. Laser gain and the fraction of the electron energy converted to photon energy are determined using a simplified resonant particle model. These results are compared to those of a more exact particle simulation code.

Hiddleston, H. R.↗

Linear analysis of an axially grooved rectangular gyrotron for harmonic operation

In an axially grooved rectangular waveguide the linearized Vlasov equation is solved to find the perturbed distribution function resulting from the electromagnetic forces on the electrons. The resulting beam current and the propagating electromagnetic waves of the cold tube are used in the inhomogeneous Maxwell's equation to derive the general dispersion equation. This equation is then transformed into the electron beam frame and the resulting linear growth rate of amplification is calculated. By maximizing the linear growth rate, the operational parameters of the gyrotron are then optimized.

Ferendeci, A. M.↗

Analytical evaluation of a surface integral expressing the coupling between interior and exterior volumes in a FE-IE approach

Recently an approach which combines the finite element technique and an integral equation to determine the fields scattered by inhomogeneous bodies of complicated shape has been proposed. Basically, a mathematical surface which encloses the scatterers is introduced, thus dividing the space into an interior and an exterior volume, in which the finite element technique and an integral equation for EM scattering, respectively, are applied. The integral equation is set up for the tangential components of the fields at the surface, while the interior volume the unknowns are the total fields. Continuity of the tangential fields at the boundary, as required by Maxwell's equations, is imposed, thus coupling the two methods to obtain a consistent solution. The coupling term is expressed by a surface integral formed by the dot product of a FE basis function and an IE testing function, or viceversa. By choosing the boundary to be a surface of revolution and by making a convenient selection of IE basis (testing) functions, it is possible to evaluate the integrals analytically on surfaces such as curved triangles, curved quadrilaterals and curved pentagons. We will illustrate the salient steps involved in setting up and carrying out these integrals and discuss what class of basis (testing) functions and analytic surfaces of revolution they are applicable to. Analytic calculations offer the advantage of better accuracy than purely numerical ones, and, when combined with them, often shed light on issues of numerical convergence and limiting values. Furthermore, they may reduce computation time and storage requirements.

Zuffada, C.↗

A Formulation of Asymptotic and Exact Boundary Conditions Using Local Operators

In this paper we describe a systematic approach for constructing asymptotic boundary conditions for isotropic wave-like equations using local operators. The conditions take a recursive form with increasing order of accuracy. In three dimensions the recursion terminates and the resulting conditions are exact for solutions which are described by finite combinations of angular spherical harmonics. First, we develop the expansion for the two-dimensional wave equation and construct a sequence of easily implementable boundary conditions. We show that in three dimensions and analogous conditions are again easily implementable in addition to being exact. Also, we provide extensions of these ideas to hyperbolic systems. Namely, Maxwell's equations for TM waves are used to demonstrate the construction. Finally, we provide numerical examples to demonstrate the effectiveness of these conditions for a model problem governed by the wave equation.

Hagstrom, T.↗

Three-Dimensional Electron Optics Model Developed for Traveling-Wave Tubes

A three-dimensional traveling-wave tube (TWT) electron beam optics model including periodic permanent magnet (PPM) focusing has been developed at the NASA Glenn Research Center at Lewis Field. This accurate model allows a TWT designer to develop a focusing structure while reducing the expensive and time-consuming task of building the TWT and hot-testing it (with the electron beam). In addition, the model allows, for the first time, an investigation of the effect on TWT operation of the important azimuthally asymmetric features of the focusing stack. The TWT is a vacuum device that amplifies signals by transferring energy from an electron beam to a radiofrequency (RF) signal. A critically important component is the focusing structure, which keeps the electron beam from diverging and intercepting the RF slow wave circuit. Such an interception can result in excessive circuit heating and decreased efficiency, whereas excessive growth in the beam diameter can lead to backward wave oscillations and premature saturation, indicating a serious reduction in tube performance. The most commonly used focusing structure is the PPM stack, which consists of a sequence of cylindrical iron pole pieces and opposite-polarity magnets. Typically, two-dimensional electron optics codes are used in the design of magnetic focusing devices. In general, these codes track the beam from the gun downstream by solving equations of motion for the electron beam in static-electric and magnetic fields in an azimuthally symmetric structure. Because these two-dimensional codes cannot adequately simulate a number of important effects, the simulation code MAFIA (solution of Maxwell's equations by the Finite-Integration-Algorithm) was used at Glenn to develop a three-dimensional electron optics model. First, a PPM stack was modeled in three dimensions. Then, the fields obtained using the magnetostatic solver were loaded into a particle-in-cell solver where the fully three-dimensional behavior of the beam was simulated in the magnetic focusing field. For the first time, the effects of azimuthally asymmetric designs and critical azimuthally asymmetric characteristics of the focusing stack (such as shunts, C-magnets, or magnet misalignment) on electron beam behavior have been investigated. A cutaway portion of a simulated electron beam focused by a PPM stack is illustrated.

Kory, Carol L.↗

Reaction formulation for radiation and scattering from plates, corner reflectors and dielectric-coated cylinders

The reaction concept is employed to formulate an integral equation for radiation and scattering from plates, corner reflectors, and dielectric-coated conducting cylinders. The surface-current density on the conducting surface is expanded with subsectional bases. The dielectric layer is modeled with polarization currents radiating in free space. Maxwell's equation and the boundary conditions are employed to express the polarization-current distribution in terms of the surface-current density on the conducting surface. By enforcing reaction tests with an array of electric test sources, the moment method is employed to reduce the integral equation to a matrix equation. Inversion of the matrix equation yields the current distribution, and the scattered field is then obtained by integrating the current distribution. The theory, computer program and numerical results are presented for radiation and scattering from plates, corner reflectors, and dielectric-coated conducting cylinders.

Wang, N. N.↗

Transport equations for multicomponent anisotropic space plasmas - A review

An attempt is made to present a unified approach to the study of transport phenomena in multicomponent anisotropic space plasmas. In particular, a system of generalized transport equations is presented that can be applied to widely different plasma flow conditions. The generalized transport equations can describe subsonic and supersonic flows, collision-dominated and collisionless flows, plasma flows in rapidly changing magnetic field configurations, multicomponent plasma flows with large temperature differences between the interacting species, and plasma flows that contain anisotropic temperature distributions. In addition, if Maxwell's equations of electricity and magnetism are added to the system of transport equations, they can be used to model electrostatic shocks, double layers, and magnetic merging processes. These transport equations also contain terms which act to regulate both the heat flow and temperature anisotropy, processes which appear to be operating in the solar wind.

Barakat, A. R.↗

Bounded Error Schemes for the Wave Equation on Complex Domains

This paper considers the application of the method of boundary penalty terms ("SAT") to the numerical solution of the wave equation on complex shapes with Dirichlet boundary conditions. A theory is developed, in a semi-discrete setting, that allows the use of a Cartesian grid on complex geometries, yet maintains the order of accuracy with only a linear temporal error-bound. A numerical example, involving the solution of Maxwell's equations inside a 2-D circular wave-guide demonstrates the efficacy of this method in comparison to others (e.g. the staggered Yee scheme) - we achieve a decrease of two orders of magnitude in the level of the L2-error.

Abarbanel, Saul↗

Maxwell's second- and third-order equations of transfer for non-Maxwellian gases

Condensed algebraic forms for Maxwell's second- and third-order equations of transfer are developed for the case of molecules described by either elastic hard spheres, inverse-power potentials, or by Bird's variable hard-sphere model. These hardly reduced, yet exact, equations provide a new point of origin, when using the moment method, in seeking approximate solutions in the kinetic theory of gases for molecular models that are physically more realistic than that provided by the Maxwell model. An important by-product of the analysis when using these second- and third-order relations is that a clear mathematical connection develops between Bird's variable hard-sphere model and that for the inverse-power potential.

Baganoff, D.↗

Characterization of thunderstorm induced Maxwell current densities in the middle atmosphere

Middle atmospheric transient Maxwell current densities generated by lightning induced charge perturbations are investigated via a simulation of Maxwell's equations. A time domain finite element analysis is employed for the simulations. The atmosphere is modeled as a region contained within a right circular cylinder with a height of 110 km and radius of 80 km. A composite conductivity profile based on measured data is used when charge perturbations are centered about the vertical axis at altitudes of 6 and 10 km. The simulations indicate that the temporal structure of the Maxwell current density is relatively insensitive to altitude variation within the region considered. It is also shown that the electric field and Maxwell current density are not generally aligned.

Baginski, Michael Edward↗

Progressive wave expansions and open boundary problems

In this paper we construct progressive wave expansions and asymptotic boundary conditions for wave-like equations in exterior domains, including applications to electromagnetics, compressible flows and aero-acoustics. The development of the conditions will be discussed in two parts. The first part will include derivations of asymptotic conditions based on the well-known progressive wave expansions for the two-dimensional wave equations. A key feature in the derivations is that the resulting family of boundary conditions involves a single derivative in the direction normal to the open boundary. These conditions are easy to implement and an application in electromagnetics will be presented. The second part of the paper will discuss the theory for hyperbolic systems in two dimensions. Here, the focus will be to obtain the expansions in a general way and to use them to derive a class of boundary conditions that involve only time derivatives or time and tangential derivatives. Maxwell's equations and the compressible Euler equations are used as examples. Simulations with the linearized Euler equations are presented to validate the theory.

Hagstrom, T.↗

SYNTHESIS OF NOVEL ALL-DIELECTRIC GRATING FILTERS USING GENETIC ALGORITHMS

We are concerned with the design of inhomogeneous, all dielectric (lossless) periodic structures which act as filters. Dielectric filters made as stacks of inhomogeneous gratings and layers of materials are being used in optical technology, but are not common at microwave frequencies. The problem is then finding the periodic cell's geometric configuration and permittivity values which correspond to a specified reflectivity/transmittivity response as a function of frequency/illumination angle. This type of design can be thought of as an inverse-source problem, since it entails finding a distribution of sources which produce fields (or quantities derived from them) of given characteristics. Electromagnetic sources (electric and magnetic current densities) in a volume are related to the outside fields by a well known linear integral equation. Additionally, the sources are related to the fields inside the volume by a constitutive equation, involving the material properties. Then, the relationship linking the fields outside the source region to those inside is non-linear, in terms of material properties such as permittivity, permeability and conductivity. The solution of the non-linear inverse problem is cast here as a combination of two linear steps, by explicitly introducing the electromagnetic sources in the computational volume as a set of unknowns in addition to the material unknowns. This allows to solve for material parameters and related electric fields in the source volume which are consistent with Maxwell's equations. Solutions are obtained iteratively by decoupling the two steps. First, we invert for the permittivity only in the minimization of a cost function and second, given the materials, we find the corresponding electric fields through direct solution of the integral equation in the source volume. The sources thus computed are used to generate the far fields and the synthesized triter response. The cost function is obtained by calculating the deviation between the synthesized value of reflectivity/transmittivity and the desired one. Solution geometries for the periodic cell are sought as gratings (ensembles of columns of different heights and widths), or combinations of homogeneous layers of different dielectric materials and gratings. Hence the explicit unknowns of the inversion step are the material permittivities and the relative boundaries separating homogeneous parcels of the periodic cell.

Zuffada, Cinzia↗

An Introduction to Atmospheric Radiation: Review for the Bulletin of AMS

Whether you like a certain geophysical book or not, largely depends on your background. The field of radiative transfer and atmospheric radiation, in particular, combines people with a wide range of mathematical skills: from theoretical astrophysicists and nuclear physicists to meteorologists and ecologists. There is always a delicate balance between physical explanations and their mathematical interpretations. This balance is very personal and is based on your background. I came to the field of atmospheric radiative transfer as a mathematician with little knowledge of atmospheric physics. After being in the field for more than a decade, I still have gaps in my atmospheric science education. Thus I assess a radiative transfer book fi-om two main criteria: how well does it describe the material that is familiar to me (the radiative transfer equation and its numerical solutions) and how well does it help me to fill the gaps in my personal knowledge. So I present this review fi-om the perspective of a former mathematician working in the field of atmospheric radiation. . After being asked to review the book, my first intention was to compare the new edition with the previous one (Liou, 1980). In doing so, you can clearly follow the progress made in the field of atmospheric radiation over the past two decades. If there are few changes (as in Fundamental Radiative Transfer) or no changes at all (as in the Maxwell s equations), then the field has not seen much development. To the contrary, many differences between the two editions illustrate areas of major progress in the field, such as evidenced in Thermal Ineared Radiative Transfer and even in the creations of completely new fields like Three-Dimensional Radiative Transfer or Light Scattering by Nonspherical Particles. Obviously, the major changes happened not in the theory, which is at least half a century old, but in data quality and completely new measurements (mostly due to new satellite data) with higher accuracy and more reliability. The new edition illustrates this progress well.

Marshak, Alexander↗

Direct Solve of Electrically Large Integral Equations for Problem Sizes to 1M Unknowns

Matrix methods for solving integral equations via direct solve LU factorization are presently limited to weeks to months of very expensive supercomputer time for problems sizes of several hundred thousand unknowns. This report presents matrix LU factor solutions for electromagnetic scattering problems for problem sizes to one million unknowns with thousands of right hand sides that run in mere days on PC level hardware. This EM solution is accomplished by utilizing the numerical low rank nature of spatially blocked unknowns using the Adaptive Cross Approximation for compressing the rank deficient blocks of the system Z matrix, the L and U factors, the right hand side forcing function and the final current solution. This compressed matrix solution is applied to a frequency domain EM solution of Maxwell's equations using standard Method of Moments approach. Compressed matrix storage and operations count leads to orders of magnitude reduction in memory and run time.

Shaeffer, John↗

Scattering and Diffraction of Electromagnetic Radiation: An Effective Probe to Material Structure

Scattered electromagnetic waves from material bodies of different forms contain, in an intricate way, precise information on the intrinsic, geometrical and physical properties of the objects. Scattering theories, ever deepening, aim to provide dependable interpretation and prediction to the complicated interaction of electromagnetic radiation with matter. There are well-established multiple-scattering formulations based on classical electromagnetic theories. An example is the Generalized Multi-particle Mie-solution (GMM), which has recently been extended to a special version ̶ the GMM-PA approach, applicable to finite periodic arrays consisting of a huge number (e.g., >>106) of identical scattering centers [1]. The framework of the GMM-PA is nearly complete. When the size of the constituent unit scatterers becomes considerably small in comparison with incident wavelength, an appropriate array of such small element volumes may well be a satisfactory representation of a material entity having an arbitrary structure. X-ray diffraction is a powerful characterization tool used in a variety of scientific and technical fields, including material science. A diffraction pattern is nothing more than the spatial distribution of scattered intensity, determined by the distribution of scattering matter by way of its Fourier transform [1]. Since all linear dimensions entered into Maxwell's equations are normalized by wavelength, an analogy exists between optical and X-ray diffraction patterns. A large set of optical diffraction patterns experimentally obtained can be found in the literature [e.g., 2,3]. Theoretical results from the GMM-PA have been scrutinized using a large collection of publically accessible, experimentally obtained Fraunhofer diffraction patterns. As far as characteristic structures of the patterns are concerned, theoretical and experimental results are in uniform agreement; no exception has been found so far. Closely connected with the spatial distribution of scattered intensities are cross sections, such as for extinction, scattering, absorption, and radiation pressure, as a critical type of key quantity addressed in most theoretical and experimental studies of radiative scattering. Cross sections predicted from different scattering theories are supposed to be in general agreement. For objects of irregular shape, the GMM-PA solutions can be compared with the highly flexible Discrete Dipole Approximation (DDA) [4,5] when dividing a target to no more than ~106 unit cells. Also, there are different ways to calculate the cross sections in the GMM-PA, providing an additional means to examine the accuracy of the numerical solutions and to unveil potential issues concerning the theoretical formulations and numerical aspects. To solve multiple scattering by an assembly of material volumes through classical theories such as the GMM-PA, the radiative properties of the component scatterers, the complex refractive index in particular, must be provided as input parameters. When using a PA to characterize a material body, this involves the use of an adequate theoretical tool, an effective medium theory, to connect Maxwell's phenomenogical theory with the atomistic theory of matter. In the atomic theory, one regards matter as composed of interacting particles (atoms and molecules) embedded in the vacuum [6]. However, the radiative properties of atomic-scaled particles are known to be substantially different from bulk materials. Intensive research efforts in the fields of cluster science and nanoscience attempt to bridge the gap between bulk and atom and to understand the transition from classical to quantum physics. The GMM-PA calculations, which place virtually no restriction on the component-particle size, might help to gain certain insight into the transition.

Xu, Yu-Lin↗

An electromagnetic Thirring problem.

A neutral rotating mass shell surrounds a concentric stationary electrically charged insulating shell. The dipole-like magnetic field induced by (and proportional to) the rotation of the neutral shell is calculated on the basis of the coupled linearized Einstein-Maxwell field equations of general relativity. This field is apparently at variance with a conjecture made on Machian grounds, for which a possible explanation is suggested. The corresponding induced quadrupolar electric field is calculated for the region within the charged shell, and the potential is given for this field everywhere. Though understandable on mutually inconsistent elementary grounds, we regard this field as a useful example of a solution of linearized general relativity.-

Ehlers, J.↗