Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ELECTROMAGNETS”

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 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Electromagnetic Detection and Real-Time DMLC Adaptation to Target Rotation During Radiotherapy

Purpose: Intrafraction rotation of more than 45 Degree-Sign and 25 Degree-Sign has been observed for lung and prostate tumors, respectively. Such rotation is not routinely adapted to during current radiotherapy, which may compromise tumor dose coverage. The aim of the study was to investigate the geometric and dosimetric performance of an electromagnetically guided real-time dynamic multileaf collimator (DMLC) tracking system to adapt to intrafractional tumor rotation. Materials/Methods: Target rotation was provided by changing the treatment couch angle. The target rotation was measured by a research Calypso system integrated with a real-time DMLC tracking system employed on a Varian linac. The geometric beam-target rotational alignment difference was measured using electronic portal images. The dosimetric accuracy was quantified using a two-dimensional ion chamber array. For each beam, the following five delivery modes were tested: 1) nonrotated target (reference); 2) fixed rotated target with tracking; 3) fixed rotated target without tracking; 4) actively rotating target with tracking; and 5) actively rotating target without tracking. Dosimetric performance of the latter four modes was measured and compared to the reference dose distribution using a 3 mm/3% {gamma}-test. Results: Geometrically, the beam-target rotational alignment difference was 0.3 Degree-Sign {+-} 0.6 Degree-Sign for fixed rotation and 0.3 Degree-Sign {+-} 1.3 Degree-Sign for active rotation. Dosimetrically, the average failure rate for the {gamma}-test for a fixed rotated target was 11% with tracking and 36% without tracking. The average failure rate for an actively rotating target was 9% with tracking and 35% without tracking. Conclusions: For the first time, real-time target rotation has been accurately detected and adapted to during radiation delivery via DMLC tracking. The beam-target rotational alignment difference was mostly within 1 Degree-Sign . Dose distributions to fixed and actively rotating targets with DMLC tracking were significantly superior to those without tracking.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Adaptive Radiation Therapy for Postprostatectomy Patients Using Real-Time Electromagnetic Target Motion Tracking During External Beam Radiation Therapy

Purpose: Using real-time electromagnetic (EM) transponder tracking data recorded by the Calypso 4D Localization System, we report inter- and intrafractional target motion of the prostate bed, describe a strategy to evaluate treatment adequacy in postprostatectomy patients receiving intensity modulated radiation therapy (IMRT), and propose an adaptive workflow. Methods and Materials: Tracking data recorded by Calypso EM transponders was analyzed for postprostatectomy patients that underwent step-and-shoot IMRT. Rigid target motion parameters during beam delivery were calculated from recorded transponder positions in 16 patients with rigid transponder geometry. The delivered doses to the clinical target volume (CTV) were estimated from the planned dose matrix and the target motion for the first 3, 5, 10, and all fractions. Treatment adequacy was determined by comparing the delivered minimum dose (D{sub min}) with the planned D{sub min} to the CTV. Treatments were considered adequate if the delivered CTV D{sub min} is at least 95% of the planned CTV D{sub min}. Results: Translational target motion was minimal for all 16 patients (mean: 0.02 cm; range: −0.12 cm to 0.07 cm). Rotational motion was patient-specific, and maximum pitch, yaw, and roll were 12.2, 4.1, and 10.5°, respectively. We observed inadequate treatments in 5 patients. In these treatments, we observed greater target rotations along with large distances between the CTV centroid and transponder centroid. The treatment adequacy from the initial 10 fractions successfully predicted the overall adequacy in 4 of 5 inadequate treatments and 10 of 11 adequate treatments. Conclusion: Target rotational motion could cause underdosage to partial volume of the postprostatectomy targets. Our adaptive treatment strategy is applicable to post-prostatectomy patients receiving IMRT to evaluate and improve radiation therapy delivery.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Respiratory-Induced Prostate Motion Using Wavelet Decomposition of the Real-Time Electromagnetic Tracking Signal

Purpose: The objective of this work is to characterize and quantify the impact of respiratory-induced prostate motion. Methods and Materials: Real-time intrafraction motion is observed with the Calypso 4-dimensional nonradioactive electromagnetic tracking system (Calypso Medical Technologies, Inc. Seattle, Washington). We report the results from a total of 1024 fractions from 31 prostate cancer patients. Wavelet transform was used to decompose the signal to extract and isolate the respiratory-induced prostate motion from the total prostate displacement. Results: Our results show that the average respiratory motion larger than 0.5 mm can be observed in 68% of the fractions. Fewer than 1% of the patients showed average respiratory motion of less than 0.2 mm, whereas 99% of the patients showed average respiratory-induced motion ranging between 0.2 and 2 mm. The maximum respiratory range of motion of 3 mm or greater was seen in only 25% of the fractions. In addition, about 2% patients showed anxiety, indicated by a breathing frequency above 24 times per minute. Conclusions: Prostate motion is influenced by respiration in most fractions. Real-time intrafraction data are sensitive enough to measure the impact of respiration by use of wavelet decomposition methods. Although the average respiratory amplitude observed in this study is small, this technique provides a tool that can be useful if one moves to smaller treatment margins (≤5 mm). This also opens ups the possibility of being able to develop patient specific margins, knowing that prostate motion is not unpredictable.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2018

A Prospective Cohort Study of Gated Stereotactic Liver Radiation Therapy Using Continuous Internal Electromagnetic Motion Monitoring

Intrafraction motion can compromise the treatment accuracy in liver stereotactic body radiation therapy (SBRT). Respiratory gating can improve treatment delivery; however, gating based on external motion surrogates is inaccurate. The present study reports the use of Calypso-based internal electromagnetic motion monitoring for gated liver SBRT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

Simulating the effects of warm O + ions on the growth of electromagnetic ion cyclotron (EMIC) waves

Electromagnetic ion cyclotron (EMIC) waves are believed to play a crucial role in the dynamics of the Earth’s magnetosphere. It has been widely accepted that plasma compositions can influence the growth rate of EMIC waves in the inner magnetosphere, but how warm O + ions change the wave growth rate is not well known. In this study, we investigate the impact of ring current O + concentration on the EMIC growth rate during a specific storm when the O + ion flux significantly increases. We calculate the growth rate of EMIC waves by using physics-based ion distributions output from the Ring current–Atmosphere interactions Model with Self-Consistent magnetic (B) and Electric (E) fields (RAM-SCBE) model. The percentage of warm O + is parametrically varied to examine how the maximum EMIC growth rate changes over time and how the global distribution of the maximum EMIC growth rate is affected. We found that the maximum growth rate of H-band appears in the dusk-to-midnight sector near the plasmapause, while O-band is excited at a slightly outer region. The maximum growth rate of He-band is closely related to the cold plasma density. With the increase of warm O + density, the maximum growth rate of H-band and He-band EMIC wave is reduced, while that of O-band EMIC wave is increased, and the region with this wave excitation is widened. Furthermore, such variation in the maximum EMIC growth rate implies a potential impact on the associated wave–particle interactions and change of the decay rate in the ring current.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

A second-order 3D electromagnetics algorithm for curved interfaces between anisotropic dielectrics on a Yee mesh

A new frequency-domain electromagnetics algorithm is developed for simulating curved interfaces between anisotropic dielectrics embedded in a Yee mesh with second-order error in resonant frequencies. The algorithm is systematically derived using the finite integration formulation of Maxwell's equations on the Yee mesh. Second-order convergence of the error in resonant frequencies is achieved by guaranteeing first-order error on dielectric boundaries and second-order error in bulk (possibly anisotropic) regions. Convergence studies, conducted for an analytically solvable problem and for a photonic crystal of ellipsoids with anisotropic dielectric constant, both show second-order convergence of frequency error; the convergence is sufficiently smooth that Richardson extrapolation yields roughly third-order convergence. The convergence of electric fields near the dielectric interface for the analytic problem is also presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Perfectly matched layers in a divergence preserving ADI scheme for electromagnetics

For numerical simulations of highly relativistic and transversely accelerated charged particles including radiation fast algorithms are needed. While the radiation in particle accelerators has wavelengths in the order of 100 {mu}m the computational domain has dimensions roughly five orders of magnitude larger resulting in very large mesh sizes. The particles are confined to a small area of this domain only. To resolve the smallest scales close to the particles subgrids are envisioned. For reasons of stability the alternating direction implicit (ADI) scheme by Smithe et al. [D.N. Smithe, J.R. Cary, J.A. Carlsson, Divergence preservation in the ADI algorithms for electromagnetics, J. Comput. Phys. 228 (2009) 7289-7299] for Maxwell equations has been adopted. At the boundary of the domain absorbing boundary conditions have to be employed to prevent reflection of the radiation. In this paper we show how the divergence preserving ADI scheme has to be formulated in perfectly matched layers (PML) and compare the performance in several scenarios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

A fast multigrid-based electromagnetic eigensolver for curved metal boundaries on the Yee mesh

For embedded boundary electromagnetics using the Dey–Mittra (Dey and Mittra, 1997) [1] algorithm, a special grad–div matrix constructed in this work allows use of multigrid methods for efficient inversion of Maxwell’s curl–curl matrix. Efficient curl–curl inversions are demonstrated within a shift-and-invert Krylov-subspace eigensolver (open-sourced at ([ofortt]https://github.com/bauerca/maxwell[cfortt])) on the spherical cavity and the 9-cell TESLA superconducting accelerator cavity. The accuracy of the Dey–Mittra algorithm is also examined: frequencies converge with second-order error, and surface fields are found to converge with nearly second-order error. In agreement with previous work (Nieter et al., 2009) [2], neglecting some boundary-cut cell faces (as is required in the time domain for numerical stability) reduces frequency convergence to first-order and surface-field convergence to zeroth-order (i.e. surface fields do not converge). Additionally and importantly, neglecting faces can reduce accuracy by an order of magnitude at low resolutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Identification of subsurface structures using electromagnetic data and shape priors

We consider the inverse problem of identifying large-scale subsurface structures using the controlled source electromagnetic method. To identify structures in the subsurface where the contrast in electric conductivity can be small, regularization is needed to bias the solution towards preserving structural information. We propose to combine two approaches for regularization of the inverse problem. In the first approach we utilize a model-based, reduced, composite representation of the electric conductivity that is highly flexible, even for a moderate number of degrees of freedom. With a low number of parameters, the inverse problem is efficiently solved using a standard, second-order gradient-based optimization algorithm. Further regularization is obtained using structural prior information, available, e.g., from interpreted seismic data. The reduced conductivity representation is suitable for incorporation of structural prior information. Such prior information cannot, however, be accurately modeled with a gaussian distribution. To alleviate this, we incorporate the structural information using shape priors. The shape prior technique requires the choice of kernel function, which is application dependent. We argue for using the conditionally positive definite kernel which is shown to have computational advantages over the commonly applied gaussian kernel for our problem. Numerical experiments on various test cases show that the methodology is able to identify fairly complex subsurface electric conductivity distributions while preserving structural prior information during the inversion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A curvilinear, fully implicit, conservative electromagnetic PIC algorithm in multiple dimensions

Here, we extend a recently proposed fully implicit PIC algorithm for the Vlasov–Darwin model in multiple dimensions (Chen and Chacón (2015) [1]) to curvilinear geometry. As in the Cartesian case, the approach is based on a potential formulation (Φ, A), and overcomes many difficulties of traditional semi-implicit Darwin PIC algorithms. Conservation theorems for local charge and global energy are derived in curvilinear representation, and then enforced discretely by a careful choice of the discretization of field and particle equations. Additionally, the algorithm conserves canonical-momentum in any ignorable direction, and preserves the Coulomb gauge ∇ • A = 0 exactly. An asymptotically well-posed fluid preconditioner allows efficient use of large cell sizes, which are determined by accuracy considerations, not stability, and can be orders of magnitude larger than required in a standard explicit electromagnetic PIC simulation. We demonstrate the accuracy and efficiency properties of the algorithm with numerical experiments in mapped meshes in 1D-3V and 2D-3V.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

GPU-advanced 3D electromagnetic simulations of superconductors in the Ginzburg–Landau formalism

Ginzburg–Landau theory is one of the most powerful phenomenological theories in physics, with particular predictive value in superconductivity. The formalism solves coupled nonlinear differential equations for both the electronic and magnetic responsiveness of a given superconductor to external electromagnetic excitations. With order parameter varying on the short scale of the coherence length, and the magnetic field being long-range, the numerical handling of 3D simulations becomes extremely challenging and time-consuming for realistic samples. Here we show precisely how one can employ graphics-processing units (GPUs) for this type of calculations, and obtain physics answers of interest in a reasonable time-frame – with speedup of over 100× compared to best available CPU implementations of the theory on a 256{sup 3} grid.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

A fully implicit, conservative, non-linear, electromagnetic hybrid particle-ion/fluid-electron algorithm

The quasi-neutral hybrid model with kinetic ions and fluid electrons is a promising approach for bridging the inherent multi-scale nature of many problems in space and laboratory plasmas. A novel, implicit, particle-in-cell based scheme for the hybrid model is derived for fully 3D electromagnetic problems with multiple ion species, which features global mass, momentum and energy conservation. The scheme includes sub-cycling and orbit-averaging for the ions, with cell-centered finite differences and implicit midpoint time advance. To reduce discrete particle noise, the scheme allows arbitrary-order shape functions for the particle-mesh interpolations and the application of conservative binomial smoothing. The algorithm is verified for a number of test problems to demonstrate the correctness of the implementation, the unique conservation properties, and the favorable stability properties of the new scheme. In particular, there is no indication of unstable growth of the finite-grid instability for a population of cold ions drifting through a uniform spatial mesh, in a set-up where several commonly used non-conservative schemes are highly unstable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

A hierarchical butterfly LU preconditioner for two-dimensional electromagnetic scattering problems involving open surfaces

This paper introduces a hierarchical interpolative decomposition butterfly-LU factorization (H-IDBF-LU) preconditioner for solving two-dimensional electric-field integral equations (EFIEs) in electromagnetic scattering problems of perfect electrically conducting objects with open surfaces. H-IDBF-LU leverages the interpolative decomposition butterfly factorization (IDBF) to compress dense blocks of the discretized EFIE operator to expedite its application; this compressed operator also serves as an approximate LU factorization of the EFIE operator leading to an efficient preconditioner in iterative solvers. Both the memory requirement and computational cost of the H-IDBF-LU solver scale as $O$($N$ log 2 $N$) in one iteration; the total number of iterations required for a reasonably good accuracy scales as $O$(1) to O(log 2 $N$) in all of our numerical tests. The efficacy and accuracy of the proposed preconditioned iterative solver are demonstrated via its application to a broad range of scatterers involving up to 100 million unknowns.

Liu, Yang↗
DOE OSTI2020

An energy-conserving and asymptotic-preserving charged-particle orbit implicit time integrator for arbitrary electromagnetic fields

Here we present a new implicit asymptotic preserving time integration scheme for charged-particle orbit computation in arbitrary electromagnetic fields. The scheme is built on the Crank-Nicolson integrator and continues to recover full-orbit motion in the small time-step limit, but also recovers all the first-order guiding center drifts as well as the correct gyroradius when stepping over the gyration time-scale. In contrast to previous efforts in this direction, the new scheme also features exact energy conservation. In the derivation of the scheme, we find that a new numerical time-scale is introduced. This scale is analyzed and the resulting restrictions on time-step are derived. Based on this analysis, we develop an adaptive time-stepping strategy the respects these constraints while stepping over the gyration scale when physically justified. It is shown through numerical tests on single-particle motion that the scheme's energy conservation property results in tremendous improvements in accuracy, and that the scheme is able to transition smoothly between magnetized and unmagnetized regimes as a result of the adaptive time-stepping.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Asymptotic-preserving gyrokinetic implicit particle-orbit integrator for arbitrary electromagnetic fields

We extend the asymptotic preserving and energy conserving time integrator for charged-particle motion developed in Ricketson and Chacón (2020) to include finite Larmor-radius (FLR) effects in the presence of electric-field length-scales comparable to the particle gyro-radius (the gyro-kinetic limit). We introduce two modifications to the earlier scheme. The first is the explicit gyro-averaging of the electric field at the half time-step, along with an analogous modification to the current deposition, which we show preserves total energy conservation in implicit PIC schemes. The number of gyrophase samples is chosen adaptively, ensuring proper averaging for large timesteps and the recovery of full-orbit dynamics in the small time-step limit. The second modification is an alternating large and small time-step strategy that ensures the particle trajectory samples gyrophases evenly. We show that this strategy relaxes the time-step restrictions on the scheme, allowing even larger speed-ups than previously achievable. We demonstrate the new method with several single-particle motion tests in a variety of electromagnetic field configurations featuring gyro-scale variation in the electric field. Finally, the results demonstrate the advertised ability to capture FLR effects accurately even when significantly stepping over the gyration time-scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Precision current regulator for NMR electromagnets

When using electromagnets or non-persistent superconducting magnets, the field stability can be directly proportional to the current stability. We present in this report a simple design for a current regulator based upon a high-precision, dc-dc fluxgate transformer. The output current and NMR frequency are stabilized by about a factor of 50, over times of order one hour. The method can be applied to essentially any current supply.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

High-pressure crystal growth and electromagnetic properties of 5d double-perovskite Ca₃OsO₆

Single crystals of the osmium-containing compound Ca₃OsO₆ have been successfully grown under high-pressure conditions, for the first time. The crystal structure of Ca₃OsO₆ atoms being fully ordered at the perovskite B-site. The electromagnetic analysis shows that the crystal exhibits a semiconductor-like behavior below 300 K and undergoes an antiferromagnetic transition at 50 K. - Graphical Abstract: Schematic image of crystal structure of Ca₃OsO₆ as determined by X-ray diffraction, where the gray and black octahedrons are occupied by Ca and Os, respectively. Top inset reveals an optic image of a typical Ca₃OsO₆ single crystal. Highlights: • Single crystals of Ca₃OsO₆ have been successfully grown under high-pressure. • Ca₃OsO₆ crystalizes into an ordered double-perovskite structure. • The Ca₃OsO₆ undergoes an antiferromagnetic transition at 50 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗