Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electromagnetic field calculations”

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 91 records · Page 5

Computational modeling of femtosecond optical solitons from Maxwell's equations

An algorithm is developed that permits the direct time integration of full-vector nonlinear Maxwell's equations. This capability permits the modeling of both linear and nonlinear instantaneous and dispersive effects in the electric polarization in material media. The modeling of the optical carrier is retained. The fundamental innovation is to notice that it is possible to treat the linear and nonlinear convolution integrals, which describe the dispersion, as new dependent variables. A coupled system of nonlinear second-order ordinary differential equations can then be derived for the linear and nonlinear convolution integrals, by differentiating them in the time domain. These equations, together with Maxwell's equations, are solved to determine the electromagnetic fields in nonlinear dispersive media. Results are presented of calculations in one dimension of the propagation and collision of femtosecond electromagnetic solitons that retain the optical carrier, taking into account as the Kerr and Raman interactions.

Goorjian, Peter M.↗

High-Power Fiber Lasers Using Photonic Band Gap Materials

High-power fiber lasers (HPFLs) would be made from photonic band gap (PBG) materials, according to the proposal. Such lasers would be scalable in the sense that a large number of fiber lasers could be arranged in an array or bundle and then operated in phase-locked condition to generate a superposition and highly directed high-power laser beam. It has been estimated that an average power level as high as 1,000 W per fiber could be achieved in such an array. Examples of potential applications for the proposed single-fiber lasers include welding and laser surgery. Additionally, the bundled fibers have applications in beaming power through free space for autonomous vehicles, laser weapons, free-space communications, and inducing photochemical reactions in large-scale industrial processes. The proposal has been inspired in part by recent improvements in the capabilities of single-mode fiber amplifiers and lasers to produce continuous high-power radiation. In particular, it has been found that the average output power of a single strand of a fiber laser can be increased by suitably changing the doping profile of active ions in its gain medium to optimize the spatial overlap of the electromagnetic field with the distribution of active ions. Such optimization minimizes pump power losses and increases the gain in the fiber laser system. The proposal would expand the basic concept of this type of optimization to incorporate exploitation of the properties (including, in some cases, nonlinearities) of PBG materials to obtain power levels and efficiencies higher than are now possible. Another element of the proposal is to enable pumping by concentrated sunlight. Somewhat more specifically, the proposal calls for exploitation of the properties of PBG materials to overcome a number of stubborn adverse phenomena that have impeded prior efforts to perfect HPFLs. The most relevant of those phenomena is amplified spontaneous emission (ASE), which causes saturation of gain and power at undesirably low levels, and scattering of light from dopants. In designing a given fiber laser for reduced ASE, care must be taken to maintain a correct fiber structure for eventual scaling to an array of many such lasers such that the interactions among all the members of the array would cause them to operate in phase lock. Hence, the problems associated with improving a single-fiber laser are not entirely separate from the bundling problem, and some designs for individual fiber lasers may be better than others if the fibers are to be incorporated into bundles. Extensive calculations, expected to take about a year, must be performed in order to determine design parameters before construction of prototype individual and fiber lasers can begin. The design effort can be expected to include calculations to optimize overlaps between the electromagnetic modes and the gain media and calculations of responses of PBG materials to electromagnetic fields. Design alternatives and physical responses that may be considered include simple PBG fibers with no intensity-dependent responses, PBG fibers with intensity- dependent band-gap shifting (see figure), and broad-band pumping made possible by use of candidate broad-band pumping media in place of the air or vacuum gaps used in prior PBG fibers.

DiDomenico, Leo↗

The effects of self-generated and applied magnetic fields on the computation of flow over a Mars return aerobrake

A CFD technique is developed to calculate the electromagnetic phenomena simultaneously with the fluid flow in the shock layer over an axisymmetric blunt body in a thermal-equilibrium chemical-nonequilibrium environment. The flowfield is solved using an explicit time-marching, first-order spatially accurate scheme. The electromagnetic phenomena are coupled to the real-gas flow solver through an iterative procedure. The electromagnetic terms introduce a strong stiffness, which was overcome by using significantly smaller time steps for the electromagnetic conservation equation. The technique is applied in calculating the flow over a Mars return aerobrake vehicle entering the Earth's atmosphere. For the case where no external field is applied, the electromagnetic effects have little impact on the flowfield.

Palmer, Grant↗

Auxiliary field diffusion Monte Carlo calculations of magnetic moments of light nuclei with chiral effective field theory interactions

We calculate the magnetic moments of light nuclei $( A < 20 )$ using the auxiliary field diffusion Monte Carlo method and local two- and three-nucleon forces with electromagnetic currents from chiral effective field theory. For all nuclei under consideration, we also calculate the ground-state energies and charge radii. We generally find a good agreement with experimental values for all of these observables. For the electromagnetic currents, we explore the impact of employing two different power counting schemes, and study theoretical uncertainties stemming from the truncation of the chiral expansion order by order for select nuclei within these two approaches. In conclusion, we find that it is crucial to employ consistent power counting schemes for interactions and currents to achieve a systematic order-by-order convergence.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Predictive Capability Maturity Model Demonstration for Cylindrical Cavity Coupling Using Gemma in the Next Generation Workflow

The predictive capability maturity model (PCMM) uses the expert elicitation process to generate credibility evidence for a particular analysis. To ensure Gemma has the capability to efficiently produce this credibility evidence, next generation workflows (NGW) are created for the solution verification, calibration/validation, and input uncertainty quantification portions of the PCMM assessment. These workflows are then used on the Higgins cylinder problem, which is representative of applications involving external-to-internal electromagnetic field coupling through a slot. The uncertainties calculated using these workflows are then used to calculate the validation comparison error and the validation uncertainty for the model following the American Society of Mechanical Engineers (ASME) verification and validation (V&V) 20 standard. These workflows will enable analysts to iterate each element of PCMM more efficiently than if completed without using a NGW workflow. An example of this iterative process is shown in Section 7.2.

42 ENGINEERING↗

Horizontal fields generated by return strokes

Horizontal fields generated by return strokes play an important role in the interaction of lightning generated electric fields with power lines. In many of the recent investigations on the interaction of lightning electromagnetic fields with power lines, the horizontal field was calculated by employing the expression for the tilt of the electric field of a plane wave propagating over finitely conducting earth. The method is suitable for calculating horizontal fields generated by return strokes at distances as close as 200m. At these close ranges, the use of the wavetilt expression can cause large errors.

Cooray, Vernon↗

Sensitivity analysis in multipole-accelerated panel methods for potential flow

In the design of an airframe, the effect of changing the geometry on resulting computations is necessary for design optimization. The geometry is defined in terms of a series of design variables, including design variables to define the wing planform, tail, canard, pylon, and nacelle. Design optimization in this research is based on how these design variable affect the potential flow. The potential flow is computed as a function of the geometry and location of a series of panels describing the airframe, which are in turn a function of the design variables. Multipole accelerated panel methods improve the computational complexity of the problem and thus are an attractive approach. To utilize the methods in design optimization, it was necessary to define the appropriate sensitivity derivatives. The overhead incurred from finding the sensitivity derivatives in conjunction with the original computation should be small. This research developed the background for multipole-accelerated panel methods and the framework for finding sensitivity derivatives in the methods. Potential flow panel codes are commonly used for powered-lift aerodynamic predictions for three dimensional geometries. Given an airframe which has been discretized into a series of panels to define the airframe geometry, potential is computed as a function of the influence of all panels on all other panels. This is a computationally intensive problem for which efficient solutions are desired to improve the computational time and to allow greater resolution by use of more panels. One such solution is the use of hierarchical multipole methods which entail approximations of the effects of far-field terms. Hierarchical multipole methods have become prevalent in molecular dynamics and gravitational physics, and have been introduced into the fields of capacitance calculations, computational fluid dynamics, and electromagnetics. The methods utilize multipole expansions to describe the effect of bodies (i.e. particles, astrophysical bodies, panels, etc.) within a sphere on points distant from the sphere, where the influence diminishes as a function of distance. The expansions are exact with infinite series, however, for practical computations, the series are truncated and accuracy is selected based on the number of terms retained in the expansions. A hierarchical tree structure groups bodies together based on proximity to allow definition of multipole expansions for each group. The multipole expansions are then used to compute the effect of the bodies in a group on distant bodies.

Leathrum, James F., Jr.↗

Transformation of Spinors in Accelerators and Beam Transfer Lines

The stable spin direction and the spin tune are two important concepts used when polarized beams are transferred along beam lines or circulating in an accelerator. One of the methods to calculate the stable spin direction and the spin tune of the reference particle of a polarized beam bunch is by integrating simultaneously the Lorentz Force $\vec{F}$ = q($\vec{E}$ +$\vec{v}$ x $\vec{B}$) and the BMT equation of the reference particle as it is transported by the magnets of a beam line or accelerator. A modified version of the Raytrace code was used to integrate the Lorentz equation and the BMT equation. Using the spin motion which is governed by the BMT equation one can compute both, the stable spin direction of the reference particle at any point along the reference orbit of an accelerator or a beam line, and the spin tune of the closed orbit in an accelerator. The simultaneous integration of the Lorentz force and of the BMT equation is the most accurate method to calculate the trajectory and the spin motion of a charged particle moving in an electromagnetic field, however this method requires the accurate description of the electromagnetic fields of the elements in which the charged particles are moving. An alternative method of calculating the stable spin direction and the spin tune of spin ½ charged particles along the reference trajectory of a beam line or an accelerator is the transformation of the spinor which is the wavefunction of spin ½ particles. This method requires the knowledge of the stable spin direction of each element of the beam line or the accelerator. This technical note describes the “spinor-method” to calculate the stable-spin- direction and the spin-tune of a spin ½ particle moving along the reference trajectory of a beam line or an accelerator.

43 PARTICLE ACCELERATORS↗

Diffraction Analysis of Antennas With Mesh Surfaces

Strip-aperture model replaces wire-grid model. Far-field radiation pattern of antenna with mesh reflector calculated more accurately with new strip-aperture model than with wire-grid model of reflector surface. More adaptable than wire-grid model to variety of practical configurations and decidedly superior for reflectors in which mesh-cell width exceeds mesh thickness. Satisfies reciprocity theorem. Applied where mesh cells are no larger than tenth of wavelength. Small cell size permits use of simplifying approximation that reflector-surface current induced by electromagnetic field is present even in apertures. Approximation useful in calculating far field.

Rahmat-Samii, Yahya↗

Field–potential finite-difference time-domain (FiPo FDTD) technique for computational electromagnetics

Modeling light–matter interactions at the nanoscale requires accurate handling of coupled quantum and electromagnetic systems. This coupling requires information about the electric scalar potential Φ and the magnetic vector potential A, which are not typically calculated in standard computational electromagnetics implementations. To that end, we have developed a field–potential finite-difference time-domain (FiPo FDTD) algorithm, which solves a set of first-order equations for Φ and A alongside equations for the electric and magnetic fields E and H. The FiPo Basic code is essentially conventional FDTD, but with an added module that calculates the potentials. The FiPo Hybrid code self-consistently calculates both fields and potentials and is particularly suitable for coupling with quantum electronic transport solvers because it can be sourced by the potentials themselves. To terminate the domain and mimic infinite space, we have derived and implemented a convolutional perfectly matched layer (CPML) absorbing boundary condition for FiPo FDTD whose performance is on par with state-of-the-art CPMLs for standard FDTD. We present FiPo simulation results on several example systems.

Avazpour, L. [University of Wisconsin-Madison, WI ↗

Magnetic Field Solver

The Magnetic Field Solver computer program calculates the magnetic field generated by a group of collinear, cylindrical axisymmetric electromagnet coils. Given the current flowing in, and the number of turns, axial position, and axial and radial dimensions of each coil, the program calculates matrix coefficients for a finite-difference system of equations that approximates a two-dimensional partial differential equation for the magnetic potential contributed by the coil. The program iteratively solves these finite-difference equations by use of the modified incomplete Cholesky preconditioned-conjugate-gradient method. The total magnetic potential as a function of axial (z) and radial (r) position is then calculated as a sum of the magnetic potentials of the individual coils, using a high-accuracy interpolation scheme. Then the r and z components of the magnetic field as functions of r and z are calculated from the total magnetic potential by use of a high-accuracy finite-difference scheme. Notably, for the finite-difference calculations, the program generates nonuniform two-dimensional computational meshes from nonuniform one-dimensional meshes. Each mesh is generated in such a way as to minimize the numerical error for a benchmark one-dimensional magnetostatic problem.

Ilin, Andrew V.↗

Radio-Frequency Rapid Thermal Processing Enabling Spatial Phase Transformation and Nanocrystallization of Soft Magnetic Amorphous Alloys

Thermal processing of soft magnetic amorphous and nanocrystalline alloys is explored under the influence of radio-frequency induction-heating techniques. Direct induction-heating concepts based on longitudinal and transverse flux heating are examined and the details of electromagnetic fields interaction with metallic strips are discussed by analytical calculations as well as finite element analysis. Initial experimental results confirming spatial control of phase transformations and nanocrystallization within a single strip of Finemet Fe-based amorphous ribbons are reported. Here, the degree to which primary and secondary crystallization temperature are achieved depends on the spacing between the ribbon relative to the induction coil as well as the coil design and configuration. For transverse coil configurations, the local temperature and therefore microstructural evolution is different across the lateral dimension of processed ribbons, with reduced gap sizes producing enhanced peak temperatures and larger temperature distributions with greater spatial variation in microstructure. In addition, indirect susceptor-based induction heating under tension is performed and the impact of microstructure is demonstrated. Herein, potential for exploiting spatially optimized phase transformations is illustrated through electromagnetic field–assisted processing in a scalable manufacturing process with amorphous and nanocrystalline soft magnetic alloys.

36 MATERIALS SCIENCE↗

The physical meaning of the magnetic scalar potential and its use in the design of hermetic electromagnetic coils

The magnetic scalar potential U(r⃗) has an intuitive physical interpretation, representing the electric current I = ΔU, which must be directed between any pair of isocontours along the boundary of a region, differing in potential by ΔU, in order to generate the corresponding magnetic intensity H⃗=−∇⃗U inside the region with no tangential component outside the boundary. This physical significance is exploited to invert the design process of hermetic (fringeless) electromagnetic coils from standard iterative techniques (calculating the magnetic field of refined current distributions) to a procedure for calculating the exact surface currents required to generate a given field configuration. A practical construction algorithm is given to produce the prescribed field in a “Target” region (constrained only by Maxwell’s equations) while confining the fringe field to a specified hermetic “Return” region. Example coils are analyzed along with a discussion of the limits of precision of the constructed field.

Crawford, C. B. (ORCID:0000000219324334)↗

A wave guide model of lightning currents and their electromagnetic field

Lightning channels are considered as resonant wave guides in which only standing resonant wave modes can be excited. Two types of discharging currents develop. Type 1 is an aperiodic wave; type 2 is a damped oscillation. The electromagnetic radiation field of both types of currents is calculated and compared with the observation.

Volland, H.↗

Theory of high-harmonic rectangular gyrotron for TE(mn) modes

In this paper, the linearized relativistic Vlasov equation is solved and a dispersion relation is calculated for the interaction between a relativistic electron beam and the electromagnetic fields of a rectangular waveguide supporting TE(mn) modes. The dispersion relation is simplified to the special case of a frame of reference moving with the electrons and the resulting coupling coefficient epsilon(mn)super l is calculated for various TE(mn) modes. The dependence of epsilon(mn)super l on the parity of the harmonic number l and various mn modes of the waveguide are discussed.

Ferendeci, A. M.↗

A VLSI Processor Engine to Support Calculation of the Vector Wave Equation

The development of a VLSI processor engine tailored to support computation of the near-field vector wave equation (VWE) is presented. Calculation of the near field is essential for evaluation of the electromagnetic fields generated by a complex radiating structure. The processor engine presented is optimized such that it takes advantage of the extensive parallelism inherent in the computational algorithm. The utilization of very large scale integration (VLSI) circuit technology and a customized reduced instruction set computer (RISC) architecture support development of a system of these processor engines, each of which is capable of performing the VWE algorithm calculation. By taking maximum advantage of the parallelism present in the calculation and dividing the problem among an array of these processor engines, the calculation can be computed in a realistic time frame.

J W Degroat↗

Ladder symmetries and Love numbers of Reissner-Nordström black holes

It is well known that asymptotically flat black holes in general relativity have vanishing tidal Love numbers. In the case of Schwarzschild and Kerr black holes, this property has been shown to be a consequence of a hidden structure of ladder symmetries for the perturbations. In this work, we extend the ladder symmetries to non-rotating charged black holes in general relativity. As opposed to previous works in this context, we adopt a more general definition of Love numbers, including quadratic operators that mix gravitational and electromagnetic perturbations in the point-particle effective field theory. We show that the calculation of a subset of those couplings in full general relativity is affected by an ambiguity in the split between source and response, which we resolve through an analytic continuation. As a result, we derive a novel master equation that unifies scalar, electromagnetic and gravitational perturbations around Reissner-Nordström black holes. The equation is hypergeometric and can be obtained from previous formulations via nontrivial field redefinitions, which allow to systematically remove some of the singularities and make the presence of the ladder symmetries more manifest.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗