3-D Relativistic MHD Simulations
NASA Astrophysics Data System (ADS)
Nishikawa, K.-I.; Frank, J.; Koide, S.; Sakai, J.-I.; Christodoulou, D. M.; Sol, H.; Mutel, R. L.
1998-12-01
We present 3-D numerical simulations of moderately hot, supersonic jets propagating initially along or obliquely to the field lines of a denser magnetized background medium with Lorentz factors of W = 4.56 and evolving in a four-dimensional spacetime. The new results are understood as follows: Relativistic simulations have consistently shown that these jets are effectively heavy and so they do not suffer substantial momentum losses and are not decelerated as efficiently as their nonrelativistic counterparts. In addition, the ambient magnetic field, however strong, can be pushed aside with relative ease by the beam, provided that the degrees of freedom associated with all three spatial dimensions are followed self-consistently in the simulations. This effect is analogous to pushing Japanese ``noren'' or vertical Venetian blinds out of the way while the slats are allowed to bend in 3-D space rather than as a 2-D slab structure.
3-D Relativistic MHD Simulations
NASA Astrophysics Data System (ADS)
Nishikaw, K.-I.; Frank, J.; Christodoulou, D. M.; Koide, S.; Sakai, J.-I.; Sol, H.; Mutel, R. L.
1998-12-01
We present 3-D numerical simulations of moderately hot, supersonic jets propagating initially along or obliquely to the field lines of a denser magnetized background medium with Lorentz factors of W=4.56 and evolving in a four-dimensional spacetime. The new results are understood as follows: Relativistic simulations have consistently shown that these jets are effectively heavy and so they do not suffer substantial momentum losses and are not decelerated as efficiently as their nonrelativistic counterparts. In addition, the ambient magnetic field, however strong, can be pushed aside with relative ease by the beam, provided that the degrees of freedom associated with all three spatial dimensions are followed self-consistently in the simulations. This effect is analogous to pushing Japanese ``noren'' or vertical Venetian blinds out of the way while the slats are allowed to bend in 3-D space rather than as a 2-D slab structure. We also simulate jets with the more realistic initial conditions for injecting jets for helical mangetic field, perturbed density, velocity, and internal energy, which are supposed to be caused in the process of jet generation. Three possible explanations for the observed variability are (i) tidal disruption of a star falling into the black hole, (ii) instabilities in the relativistic accretion disk, and (iii) jet-related PRocesses. New results will be reported at the meeting.
3D Hydrodynamic Simulations of Relativistic Jets
NASA Astrophysics Data System (ADS)
Hughes, P. A.; Miller, M. A.; Duncan, G. C.; Swift, C. M.
1998-12-01
We present the results of validation runs and the first extragalactic jet simulations performed with a 3D relativistic numerical hydrodynamic code employing a solver of the RHLLE type and using adaptive mesh refinement (AMR; Duncan & Hughes, 1994, Ap. J., 436, L119). Test problems include the shock tube, blast wave and spherical shock reflection (implosion). Trials with the code show that as a consequence of AMR it is viable to perform exploratory runs on workstation class machines (with no more than 128Mb of memory) prior to production runs. In the former case we achieve a resolution not much less than that normally regarded as the minimum needed to capture the essential physics of a problem, which means that such runs can provide valuable guidance allowing the optimum use of supercomputer resources. We present initial results from a program to explore the 3D stability properties of flows previously studied using a 2D axisymmetric code, and our first attempt to explore the structure and morphology of a relativistic jet encountering an ambient density gradient that mimics an ambient inhomogeneity or cloud.
3D Electromagnetic inversion using conjugate gradients
Newman, G.A.; Alumbaugh, D.L.
1997-06-01
In large scale 3D EM inverse problems it may not be possible to directly invert a full least-squares system matrix involving model sensitivity elements. Thus iterative methods must be employed. For the inverse problem, we favor either a linear or non-linear (NL) CG scheme, depending on the application. In a NL CG scheme, the gradient of the objective function is required at each relaxation step along with a univariate line search needed to determine the optimum model update. Solution examples based on both approaches will be presented.
3-D Relativistic MHD Simulations of Extragalactic Jets
NASA Astrophysics Data System (ADS)
Nishikawa, K.-I.; Koide, S.; Sakai, J.-I.; Frank, J.; Christodoulou, D. M.; Sol, H.; Mutel, R. L.
1997-12-01
We present the numerical simulations of relativistic jets propagating initially oblique to the field lines of a magnetized ambient medium. Our simulations incorporate relativistic MHD in a four-dimensional spacetime and clearly show that (a) relatively weak, oblique fields (at 1/16 of the equipartition value) have only a negligible influence on the propagating jet and they are passively pushed away by the relativistically moving head; (b) oblique fields in equipartition with the ambient plasma provide more resistance and cause bending at the jet head, but the magnitude of this deflection and the associated backflow are small compared to those identified by previous studies with a 2-D slab model. The new results are understood as follows: Relativistic simulations have consistently shown that these jets are effectively heavy and so they do not suffer substantial momentum losses and are not decelerated as efficiently as their nonrelativistic counterparts. In addition, the ambient magnetic field, however strong, can be pushed aside with relative ease by the beam, provided that the degrees of freedom associated with all three spatial dimensions are followed self-consistently during the simulations. The effect is analogous to pushing Japanese ``noren'' or vertical Venetian blinds out of the way while the slats are allowed to bend in 3-D space rather than as a 2-D slab structure. Applied to relativistic extragalactic jets from blazars, the new results are encouraging since superluminal outflows exhibit bending near their sources and their environments are profoundly magnetized---but observations do not provide support for irregular kinematics such as large-scale vortical motions and pronounced reverse flows near the points of origin.
Photospheric Emission from Collapsar Jets in 3D Relativistic Hydrodynamics
NASA Astrophysics Data System (ADS)
Ito, Hirotaka; Matsumoto, Jin; Nagataki, Shigehiro; Warren, Donald C.; Barkov, Maxim V.
2015-12-01
We explore the photospheric emission from a relativistic jet breaking out from a massive stellar envelope based on relativistic hydrodynamical simulations and post-process radiation transfer calculations in three dimensions. To investigate the impact of three-dimensional (3D) dynamics on the emission, two models of injection conditions are considered for the jet at the center of the progenitor star: one with periodic precession and another without precession. We show that structures developed within the jet due to the interaction with the stellar envelope, as well as due to the precession, have a significant imprint on the resulting emission. Particularly, we find that the signature of precession activity by the central engine is not smeared out and can be directly observed in the light curve as a periodic signal. We also show that non-thermal features, which can account for observations of gamma-ray bursts, are produced in the resulting spectra even though only thermal photons are injected initially and the effect of non-thermal particles is not considered.
3-D electromagnetic modeling of wakefields in accelerator components
Poole, B.R.; Caporaso, G.J.; Ng, Wang C.; Shang, C.C.; Steich, D.
1996-09-18
We discuss the use of 3-D finite-difference time-domain (FDTD) electromagnetic codes for modeling accelerator components. Computational modeling of cylindrically symmetric structures such as induction accelerator cells has been very successful in predicting the wake potential and wake impedances of these structures, but full 3-D modeling of complex structures has been limited due to substantial computer resources required for a full 3-D model. New massively parallel 3-D time domain electromagnetic codes now under development using conforming unstructured meshes allow a substantial increase in the geometric fidelity of the structures being modeled. Development of these new codes are discussed in context of applicability to accelerator problems. Various 3-D structures are tested with an existing cubical cell FDTD code and wake impedances compared with simple analytic models for the structures; results will be used as benchmarks for testing the new time time domain codes. Structures under consideration include a stripline beam position monitor as well as circular and elliptical apertures in circular waveguides. Excellent agreement for monopole and dipole impedances with models were found for these structures below the cutoff frequency of the beam line.
3-D electromagnetic modeling of wakefields in accelerator components
Poole, Brian R.; Caporaso, George J.; Ng, Wang C.; Shang, Clifford C.; Steich, David
1997-02-01
We discuss the use of 3-D finite-difference time-domain (FDTD) electromagnetic codes for the modeling of accelerator components. Computational modeling of cylindrically symmetric structures such as induction accelerator cells has been extremely successful in predicting the wake potential and wake impedances of these structures, but fully 3-D modeling of complex structures has been limited due to the substantial computer resources required for a fully 3-D model. New massively parallel 3-D time domain electromagnetic codes now under development using conforming unstructured meshes allow a substantial increase in the geometric fidelity of the structures being modeled. Development of these new codes will be discussed in the context of their applicability to accelerator problems. A variety of 3-D structures are tested with an existing cubical cell FDTD code and the wake impedances are compared with simple analytic models for the structures. These results will provide a set of benchmarks for testing the new time domain codes. Structures under consideration include a stripline beam position monitor as well as circular and elliptical apertures in circular waveguides. Excellent agreement for the monopole and dipole impedances with the models are found for these structures below the cutoff frequency of the beam line.
Advances in 3D electromagnetic finite element modeling
Nelson, E.M.
1997-08-01
Numerous advances in electromagnetic finite element analysis (FEA) have been made in recent years. The maturity of frequency domain and eigenmode calculations, and the growth of time domain applications is briefly reviewed. A high accuracy 3D electromagnetic finite element field solver employing quadratic hexahedral elements and quadratic mixed-order one-form basis functions will also be described. The solver is based on an object-oriented C++ class library. Test cases demonstrate that frequency errors less than 10 ppm can be achieved using modest workstations, and that the solutions have no contamination from spurious modes. The role of differential geometry and geometrical physics in finite element analysis is also discussed.
Electromagnetic processes in relativistic heavy ion collisions
NASA Astrophysics Data System (ADS)
Bertulani, C. A.; Baur, G.
1986-10-01
Electromagnetic effects in relativistic heavy ion collisions with impact parameter larger than the sum of the nuclear radii are studied using the virtual photon method. With increasing value of the relativistic parameter γ the hardness of the virtual photon spectrum increases. This leads to interesting new effects which will also have to be considered in the design of future relativistic heavy ion machines and experiments. The excitation of high-lying giant E1 and E2 multipole resonances is calculated as well as electromagnetic pion production. Coulomb bremsstrahlung is calculated and compared to the bremsstrahlung emitted in the more violent central nuclear collisions. K-shell ionization and electron-positron pair production is studied. The latter process has a very large cross section for heavy ions and contributes significantly to the stopping power of relativistic heavy ions in a dense medium.
Relativistic diffusive motion in random electromagnetic fields
NASA Astrophysics Data System (ADS)
Haba, Z.
2011-08-01
We show that the relativistic dynamics in a Gaussian random electromagnetic field can be approximated by the relativistic diffusion of Schay and Dudley. Lorentz invariant dynamics in the proper time leads to the diffusion in the proper time. The dynamics in the laboratory time gives the diffusive transport equation corresponding to the Jüttner equilibrium at the inverse temperature β-1 = mc2. The diffusion constant is expressed by the field strength correlation function (Kubo's formula).
Relativistic diffusive motion in thermal electromagnetic fields
NASA Astrophysics Data System (ADS)
Haba, Z.
2013-04-01
We discuss relativistic dynamics in a random electromagnetic field which can be considered as a high temperature limit of the quantum electromagnetic field in a heat bath (cavity) moving with a uniform velocity w. We derive a diffusion approximation for the particle’s dynamics generalizing the diffusion of Schay and Dudley. It is shown that the Jüttner distribution is the equilibrium state of the diffusion.
Image Appraisal for 2D and 3D Electromagnetic Inversion
Alumbaugh, D.L.; Newman, G.A.
1999-01-28
Linearized methods are presented for appraising image resolution and parameter accuracy in images generated with two and three dimensional non-linear electromagnetic inversion schemes. When direct matrix inversion is employed, the model resolution and posterior model covariance matrices can be directly calculated. A method to examine how the horizontal and vertical resolution varies spatially within the electromagnetic property image is developed by examining the columns of the model resolution matrix. Plotting the square root of the diagonal of the model covariance matrix yields an estimate of how errors in the inversion process such as data noise and incorrect a priori assumptions about the imaged model map into parameter error. This type of image is shown to be useful in analyzing spatial variations in the image sensitivity to the data. A method is analyzed for statistically estimating the model covariance matrix when the conjugate gradient method is employed rather than a direct inversion technique (for example in 3D inversion). A method for calculating individual columns of the model resolution matrix using the conjugate gradient method is also developed. Examples of the image analysis techniques are provided on 2D and 3D synthetic cross well EM data sets, as well as a field data set collected at the Lost Hills Oil Field in Central California.
Relativistically strong electromagnetic radiation in a plasma
NASA Astrophysics Data System (ADS)
Bulanov, S. V.; Esirkepov, T. Zh.; Kando, M.; Kiriyama, H.; Kondo, K.
2016-03-01
Physical processes in a plasma under the action of relativistically strong electromagnetic waves generated by high-power lasers have been briefly reviewed. These processes are of interest in view of the development of new methods for acceleration of charged particles, creation of sources of bright hard electromagnetic radiation, and investigation of macroscopic quantum-electrodynamical processes. Attention is focused on nonlinear waves in a laser plasma for the creation of compact electron accelerators. The acceleration of plasma bunches by the radiation pressure of light is the most efficient regime of ion acceleration. Coherent hard electromagnetic radiation in the relativistic plasma is generated in the form of higher harmonics and/or electromagnetic pulses, which are compressed and intensified after reflection from relativistic mirrors created by nonlinear waves. In the limit of extremely strong electromagnetic waves, radiation friction, which accompanies the conversion of radiation from the optical range to the gamma range, fundamentally changes the behavior of the plasma. This process is accompanied by the production of electron-positron pairs, which is described within quantum electrodynamics theory.
Relativistic MHD simulations of core-collapse GRB jets: 3D instabilities and magnetic dissipation
NASA Astrophysics Data System (ADS)
Bromberg, Omer; Tchekhovskoy, Alexander
2016-02-01
Relativistic jets are associated with extreme astrophysical phenomena, like the core collapse of massive stars in gamma-ray bursts (GRBs) and the accretion on to supermassive black holes in active galactic nuclei. It is generally accepted that these jets are powered electromagnetically, by the magnetized rotation of a central compact object (black hole or neutron star). However, how the jets produce the observed emission and survive the propagation for many orders of magnitude in distance without being disrupted by current-driven instabilities is the subject of active debate. We carry out time-dependent 3D relativistic magnetohydrodynamic (MHD) simulations of relativistic, Poynting-flux-dominated jets. The jets are launched self-consistently by the rotation of a strongly magnetized central object. This determines the natural degree of azimuthal magnetic field winding, a crucial factor that controls jet stability. We find that the jets are susceptible to two types of instability: (i) a global, external kink mode that grows on long time-scales. It bodily twists the jet, reducing its propagation velocity. We show analytically that in flat density profiles, like the ones associated with galactic cores, the external mode grows and may stall the jet. In the steep profiles of stellar envelopes the external kink weakens as the jet propagates outward. (ii) a local, internal kink mode that grows over short time-scales and causes small-angle magnetic reconnection and conversion of about half of the jet electromagnetic energy flux into heat. We suggest that internal kink instability is the main dissipation mechanism responsible for powering GRB prompt emission.
3D synthetic aperture for controlled-source electromagnetics
NASA Astrophysics Data System (ADS)
Knaak, Allison
Locating hydrocarbon reservoirs has become more challenging with smaller, deeper or shallower targets in complicated environments. Controlled-source electromagnetics (CSEM), is a geophysical electromagnetic method used to detect and derisk hydrocarbon reservoirs in marine settings, but it is limited by the size of the target, low-spatial resolution, and depth of the reservoir. To reduce the impact of complicated settings and improve the detecting capabilities of CSEM, I apply synthetic aperture to CSEM responses, which virtually increases the length and width of the CSEM source by combining the responses from multiple individual sources. Applying a weight to each source steers or focuses the synthetic aperture source array in the inline and crossline directions. To evaluate the benefits of a 2D source distribution, I test steered synthetic aperture on 3D diffusive fields and view the changes with a new visualization technique. Then I apply 2D steered synthetic aperture to 3D noisy synthetic CSEM fields, which increases the detectability of the reservoir significantly. With more general weighting, I develop an optimization method to find the optimal weights for synthetic aperture arrays that adapts to the information in the CSEM data. The application of optimally weighted synthetic aperture to noisy, simulated electromagnetic fields reduces the presence of noise, increases detectability, and better defines the lateral extent of the target. I then modify the optimization method to include a term that minimizes the variance of random, independent noise. With the application of the modified optimization method, the weighted synthetic aperture responses amplifies the anomaly from the reservoir, lowers the noise floor, and reduces noise streaks in noisy CSEM responses from sources offset kilometers from the receivers. Even with changes to the location of the reservoir and perturbations to the physical properties, synthetic aperture is still able to highlight targets
Relativistic spin operators in various electromagnetic environments
NASA Astrophysics Data System (ADS)
Bauke, Heiko; Ahrens, Sven; Keitel, Christoph H.; Grobe, Rainer
2014-05-01
Different operators have been suggested in the literature to describe the electron's spin degree of freedom within the relativistic Dirac theory. We compare concrete predictions of the various proposed relativistic spin operators in different physical situations. In particular, we investigate the so-called Pauli, Foldy-Wouthuysen, Czachor, Frenkel, Chakrabarti, Pryce, and Fradkin-Good spin operators. We demonstrate that when a quantum system interacts with electromagnetic potentials the various spin operators predict different expectation values. This is explicitly illustrated for the scattering dynamics at a potential step and in a standing laser field and also for energy eigenstates of hydrogenic ions. Therefore, one may distinguish between the proposed relativistic spin operators experimentally.
Solution accelerators for large scale 3D electromagnetic inverse problems
Newman, Gregory A.; Boggs, Paul T.
2004-04-05
We provide a framework for preconditioning nonlinear 3D electromagnetic inverse scattering problems using nonlinear conjugate gradient (NLCG) and limited memory (LM) quasi-Newton methods. Key to our approach is the use of an approximate adjoint method that allows for an economical approximation of the Hessian that is updated at each inversion iteration. Using this approximate Hessian as a preconditoner, we show that the preconditioned NLCG iteration converges significantly faster than the non-preconditioned iteration, as well as converging to a data misfit level below that observed for the non-preconditioned method. Similar conclusions are also observed for the LM iteration; preconditioned with the approximate Hessian, the LM iteration converges faster than the non-preconditioned version. At this time, however, we see little difference between the convergence performance of the preconditioned LM scheme and the preconditioned NLCG scheme. A possible reason for this outcome is the behavior of the line search within the LM iteration. It was anticipated that, near convergence, a step size of one would be approached, but what was observed, instead, were step lengths that were nowhere near one. We provide some insights into the reasons for this behavior and suggest further research that may improve the performance of the LM methods.
Image appraisal for 2D and 3D electromagnetic inversion
Alumbaugh, D.L.; Newman, G.A.
1998-04-01
Linearized methods are presented for appraising image resolution and parameter accuracy in images generated with two and three dimensional non-linear electromagnetic inversion schemes. When direct matrix inversion is employed, the model resolution and model covariance matrices can be directly calculated. The columns of the model resolution matrix are shown to yield empirical estimates of the horizontal and vertical resolution throughout the imaging region. Plotting the square root of the diagonal of the model covariance matrix yields an estimate of how the estimated data noise maps into parameter error. When the conjugate gradient method is employed rather than a direct inversion technique (for example in 3D inversion), an iterative method can be applied to statistically estimate the model covariance matrix, as well as a regularization covariance matrix. The latter estimates the error in the inverted results caused by small variations in the regularization parameter. A method for calculating individual columns of the model resolution matrix using the conjugate gradient method is also developed. Examples of the image analysis techniques are provided on a synthetic cross well EM data set.
3D electromagnetic simulation of spatial autoresonance acceleration of electron beams
NASA Astrophysics Data System (ADS)
Dugar-Zhabon, V. D.; González, J. D.; Orozco, E. A.
2016-02-01
The results of full electromagnetic simulations of the electron beam acceleration by a TE 112 linear polarized electromagnetic field through Space Autoresonance Acceleration mechanism are presented. In the simulations, both the self-sustaned electric field and selfsustained magnetic field produced by the beam electrons are included into the elaborated 3D Particle in Cell code. In this system, the space profile of the magnetostatic field maintains the electron beams in the acceleration regime along their trajectories. The beam current density evolution is calculated applying the charge conservation method. The full magnetic field in the superparticle positions is found by employing the trilinear interpolation of the mesh node data. The relativistic Newton-Lorentz equation presented in the centered finite difference form is solved using the Boris algorithm that provides visualization of the beam electrons pathway and energy evolution. A comparison between the data obtained from the full electromagnetic simulations and the results derived from the motion equation depicted in an electrostatic approximation is carried out. It is found that the self-sustained magnetic field is a factor which improves the resonance phase conditions and reduces the beam energy spread.
3D Relativistic Magnetohydrodynamic Simulations of Magnetized Spine-Sheath Relativistic Jets
NASA Technical Reports Server (NTRS)
Mizuno, Yosuke; Hardee, Philip; Nishikawa, Ken-Ichi
2006-01-01
Numerical simulations of weakly magnetized and strongly magnetized relativistic jets embedded in a weakly magnetized and strongly magnetized stationary or weakly relativistic (v = c/2) sheath have been performed. A magnetic field parallel to the flow is used in these simulations performed by the new GRMHD numerical code RAISHIN used in its RMHD configuration. In the numerical simulations the Lorentz factor gamma = 2.5 jet is precessed to break the initial equilibrium configuration. In the simulations sound speeds are less than or equal to c/the square root of 3 in the weakly magnetized simulations and less than or equal to 0.56 c in the strongly magnetized simulations. The Alfven wave speed is less than or equal to 0.07 c in the weakly magnetized simulations and less than or equal to 0.56 c in the strongly magnetized simulations. The results of the numerical simulations are compared to theoretical predictions from a normal mode analysis of the linearized relativistic magnetohydrodynamic (RMHD) equations capable of describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized relativistically moving sheath. The theoretical dispersion relation allows investigation of effects associated with maximum possible sound speeds, Alfven wave speeds near light speed and relativistic sheath speeds. The prediction of increased stability of the weakly magnetized system resulting from c/2 sheath speeds and the stabilization of the strongly magnetized system resulting from c/2 sheath speeds is verified by the numerical simulation results.
Current loop coalescence studied by 3-D electromagnetic particle code
NASA Technical Reports Server (NTRS)
Nishikawa, Ken-Ichi; Sakai, Jun-Ichi; Koide, Shinji; Buneman, O.; Neubert, T.
1993-01-01
Solar flare plasma data from the Yohkoh satellite is analyzed. The interactions of current loops were observed in the active regions on the Sun. This observation pointed out the importance of the idea that the solar flare is generated by the coalescence of current loops. The three dimensional electromagnetic particle simulations are to help in understanding the global interaction between two current loops including the evolution of the twist of loops due to instabilities. Associated rapid dynamics of current loop coalescence such as reconnection, shock waves and associated kinetic processes such as energy transfer, acceleration of particles, and electromagnetic emissions are to be studied by the code to complement analytical theories and magnetohydrodynamic simulations of the current loop coalescence. The simulation results show the strong interactions between two current loops, beam and whistler instabilities, and associated parallel and perpendicular particle heating.
Exact Relativistic Ideal Hydrodynamical Solutions in (1+3)D with Longitudinal and Transverse Flows
Liao, Jinfeng; Koch, Volker
2009-05-20
A new method for solving relativistic ideal hydrodynamics in (1+3)D is developed. Longitudinal and transverse radial flows are explicitly embedded into the ansatz for velocity field and the hydrodynamic equations are reduced to a single equation for the transverse velocity field only, which is analytically more tractable as compared with the full hydrodynamic equations. As an application we use the method to find analytically all possible solutions whose transverse velocity fields have power dependence on proper time and transverse radius. Possible application to the Relativistic Heavy Ion Collisions and possible generalizations of the method are discussed.
Development of 3D electromagnetic modeling tools for airborne vehicles
NASA Technical Reports Server (NTRS)
Volakis, John L.
1992-01-01
The main goal of this report is to advance the development of methodologies for scattering by airborne composite vehicles. Although the primary focus continues to be the development of a general purpose computer code for analyzing the entire structure as a single unit, a number of other tasks are also being pursued in parallel with this effort. One of these tasks discussed within is on new finite element formulations and mesh termination schemes. The goal here is to decrease computation time while retaining accuracy and geometric adaptability.The second task focuses on the application of wavelets to electromagnetics. Wavelet transformations are shown to be able to reduce a full matrix to a band matrix, thereby reducing the solutions memory requirements. Included within this document are two separate papers on finite element formulations and wavelets.
Fast Simulation of 3D Electromagnetic Problems Using Potentials
NASA Astrophysics Data System (ADS)
Haber, E.; Ascher, U. M.; Aruliah, D. A.; Oldenburg, D. W.
2000-09-01
We consider solving three-dimensional electromagnetic problems in parameter regimes where the quasi-static approximation applies, the permeability is constant, the conductivity may vary significantly, and the range of frequencies is moderate. The difficulties encountered include handling solution discontinuities across interfaces and accelerating convergence of traditional iterative methods for the solution of the linear systems of algebraic equations that arise when discretizing Maxwell's equations in the frequency domain. We use a potential-current formulation (A, φ ,Ĵ) with a Coulomb gauge. The potentials A and φ decompose the electric field E into components in the active and null spaces of the ∇×operator. We develop a finite volume discretization on a staggered grid that naturally employs harmonic averages for the conductivity at cell faces. After discretization, we eliminate the current and the resulting large, sparse, linear system of equations has a block structure that is diagonally dominant, allowing an efficient solution with preconditioned Krylov space methods. A particularly efficient algorithm results from the combination of BICGSTAB and an incomplete LU-decomposition. We demonstrate the efficacy of our method in several numerical experiments.
Time efficient 3-D electromagnetic modeling on massively parallel computers
Alumbaugh, D.L.; Newman, G.A.
1995-08-01
A numerical modeling algorithm has been developed to simulate the electromagnetic response of a three dimensional earth to a dipole source for frequencies ranging from 100 to 100MHz. The numerical problem is formulated in terms of a frequency domain--modified vector Helmholtz equation for the scattered electric fields. The resulting differential equation is approximated using a staggered finite difference grid which results in a linear system of equations for which the matrix is sparse and complex symmetric. The system of equations is solved using a preconditioned quasi-minimum-residual method. Dirichlet boundary conditions are employed at the edges of the mesh by setting the tangential electric fields equal to zero. At frequencies less than 1MHz, normal grid stretching is employed to mitigate unwanted reflections off the grid boundaries. For frequencies greater than this, absorbing boundary conditions must be employed by making the stretching parameters of the modified vector Helmholtz equation complex which introduces loss at the boundaries. To allow for faster calculation of realistic models, the original serial version of the code has been modified to run on a massively parallel architecture. This modification involves three distinct tasks; (1) mapping the finite difference stencil to a processor stencil which allows for the necessary information to be exchanged between processors that contain adjacent nodes in the model, (2) determining the most efficient method to input the model which is accomplished by dividing the input into ``global`` and ``local`` data and then reading the two sets in differently, and (3) deciding how to output the data which is an inherently nonparallel process.
Electromagnetic wave equations for relativistically degenerate quantum magnetoplasmas.
Masood, Waqas; Eliasson, Bengt; Shukla, Padma K
2010-06-01
A generalized set of nonlinear electromagnetic quantum hydrodynamic (QHD) equations is derived for a magnetized quantum plasma, including collisional, electron spin- 1/2, and relativistically degenerate electron pressure effects that are relevant for dense astrophysical systems, such as white dwarfs. For illustrative purposes, linear dispersion relations are derived for one-dimensional magnetoacoustic waves for a collisionless nonrelativistic degenerate gas in the presence of the electron spin- 1/2 contribution and for magnetoacoustic waves in a plasma containing relativistically degenerate electrons. It is found that both the spin and relativistic degeneracy at high densities tend to slow down the magnetoacoustic wave due to the Pauli paramagnetic effect and relativistic electron mass increase. The present study outlines the theoretical framework for the investigation of linear and nonlinear behaviors of electromagnetic waves in dense astrophysical systems. The results are applied to calculate the magnetoacoustic speeds for both the nonrelativistic and relativistic electron degeneracy cases typical for white dwarf stars. PMID:20866534
Kinetic theory of electromagnetic ion waves in relativistic plasmas
Marklund, Mattias; Shukla, Padma K.
2006-09-15
A kinetic theory for electromagnetic ion waves in a cold relativistic plasma is derived. The kinetic equation for the broadband electromagnetic ion waves is coupled to the slow density response via an acoustic equation driven by a ponderomotive force-like term linear in the electromagnetic field amplitude. The modulational instability growth rate is derived for an arbitrary spectrum of waves. The monochromatic and random phase cases are studied.
HERO - A 3D general relativistic radiative post-processor for accretion discs around black holes
NASA Astrophysics Data System (ADS)
Zhu, Yucong; Narayan, Ramesh; Sadowski, Aleksander; Psaltis, Dimitrios
2015-08-01
HERO (Hybrid Evaluator for Radiative Objects) is a 3D general relativistic radiative transfer code which has been tailored to the problem of analysing radiation from simulations of relativistic accretion discs around black holes. HERO is designed to be used as a post-processor. Given some fixed fluid structure for the disc (i.e. density and velocity as a function of position from a hydrodynamic or magnetohydrodynamic simulation), the code obtains a self-consistent solution for the radiation field and for the gas temperatures using the condition of radiative equilibrium. The novel aspect of HERO is that it combines two techniques: (1) a short-characteristics (SC) solver that quickly converges to a self-consistent disc temperature and radiation field, with (2) a long-characteristics (LC) solver that provides a more accurate solution for the radiation near the photosphere and in the optically thin regions. By combining these two techniques, we gain both the computational speed of SC and the high accuracy of LC. We present tests of HERO on a range of 1D, 2D, and 3D problems in flat space and show that the results agree well with both analytical and benchmark solutions. We also test the ability of the code to handle relativistic problems in curved space. Finally, we discuss the important topic of ray defects, a major limitation of the SC method, and describe our strategy for minimizing the induced error.
3D inversion of airborne electromagnetic data using a moving footprint
NASA Astrophysics Data System (ADS)
Cox, Leif H.; Wilson, Glenn A.; Zhdanov, Michael S.
2010-12-01
It is often argued that 3D inversion of entire airborne electromagnetic (AEM) surveys is impractical, and that 1D methods provide the only viable option for quantitative interpretation. However, real geological formations are 3D by nature and 3D inversion is required to produce accurate images of the subsurface. To that end, we show that it is practical to invert entire AEM surveys to 3D conductivity models with hundreds of thousands if not millions of elements. The key to solving a 3D AEM inversion problem is the application of a moving footprint approach. We have exploited the fact that the area of the footprint of an AEM system is significantly smaller than the area of an AEM survey, and developed a robust 3D inversion method that uses a moving footprint. Our implementation is based on the 3D integral equation method for computing data and sensitivities, and uses the re-weighted regularised conjugate gradient method for minimising the objective functional. We demonstrate our methodology with the 3D inversion of AEM data acquired for salinity mapping over the Bookpurnong Irrigation District in South Australia. We have inverted 146 line km of RESOLVE data for a 3D conductivity model with ~310000 elements in 45min using just five processors of a multi-processor workstation.
Relativistic transport equations for electromagnetic scalar, and pseudoscalar potentials
Shin, G.R.; Rafelski, J.
1995-10-01
The authors propose a particular form of relativistic transport equations arising from the classical limit of single-time Wigner function for Dirac particles evolving in the presence of scalar, pseudoscalar, and electromagnetic fields. These relativistic Vlasov-type equations for the particle and the antiparticle sector of the Fock space can be also obtained assuming the validity of the Liouville`s equation given a suitable classical Hamiltonian and the associated force. 11 refs.
NASA Astrophysics Data System (ADS)
Narayan, Ramesh; Zhu, Yucong; Psaltis, Dimitrios; Saḑowski, Aleksander
2016-03-01
We describe Hybrid Evaluator for Radiative Objects Including Comptonization (HEROIC), an upgraded version of the relativistic radiative post-processor code HERO described in a previous paper, but which now Includes Comptonization. HEROIC models Comptonization via the Kompaneets equation, using a quadratic approximation for the source function in a short characteristics radiation solver. It employs a simple form of accelerated lambda iteration to handle regions of high scattering opacity. In addition to solving for the radiation field, HEROIC also solves for the gas temperature by applying the condition of radiative equilibrium. We present benchmarks and tests of the Comptonization module in HEROIC with simple 1D and 3D scattering problems. We also test the ability of the code to handle various relativistic effects using model atmospheres and accretion flows in a black hole space-time. We present two applications of HEROIC to general relativistic magnetohydrodynamics simulations of accretion discs. One application is to a thin accretion disc around a black hole. We find that the gas below the photosphere in the multidimensional HEROIC solution is nearly isothermal, quite different from previous solutions based on 1D plane parallel atmospheres. The second application is to a geometrically thick radiation-dominated accretion disc accreting at 11 times the Eddington rate. Here, the multidimensional HEROIC solution shows that, for observers who are on axis and look down the polar funnel, the isotropic equivalent luminosity could be more than 10 times the Eddington limit, even though the spectrum might still look thermal and show no signs of relativistic beaming.
Relativistic 3D precessing jet simulations for the X-ray binary SS433
NASA Astrophysics Data System (ADS)
Monceau-Baroux, Rémi; Porth, Oliver; Meliani, Zakaria; Keppens, Rony
2014-01-01
Context. Modern high-resolution radio observations allow us a closer look into the objects that power relativistic jets. This is especially the case for SS433, an X-ray binary that emits a precessing jet that is observed down to the subparsec scale. Aims: We aim to study full 3D dynamics of relativistic jets associated with active galactic nuclei or X-ray binaries (XRB). In particular, we incorporate the precessing motion of a jet into a model for the jet associated with the XRB SS433. Our study of the jet dynamics in this system focuses on the subparsec scales. We investigate the impact of jet precession and the variation of the Lorentz factor of the injected matter on the general 3D jet dynamics and its energy transfer to the surrounding medium. After visualizing and quantifying jet dynamics, we aim to realize synthetic radio mapping of the data, to compare our results with observations. Methods: For our study we used a block-tree adaptive mesh refinement scheme and an inner time-dependent boundary prescription to inject precessing bipolar supersonic jets. Parameters extracted from observations were used. Different 3D jet realizations that match the kinetic flux of the SS433 jet were intercompared, which vary in density contrast and jet beam velocity. We tracked the energy content deposited in different regions of the domain affected by the jet. Our code allows us to follow the adiabatic cooling of a population of relativistic particles injected by the jet. This evolving energy spectrum of accelerated electrons, using a pressure-based proxy for the magnetic field, allowed us to obtain the radio emission from our simulation. Results: We find a higher energy transfer for a precessing jet than for standing jets with otherwise identical parameters as a result of the effectively increased interaction area. We obtain synthetic radio maps for all jets, from which one can see that dynamical flow features are clearly linked with enhanced emission sites. Conclusions: The
Preparation for a 3D Electromagnetic inversion-Application to GREATEM data
NASA Astrophysics Data System (ADS)
Abd allah, S.; Mogi, T.; Kim, H.; Fomenko, E.
2013-12-01
Previous studies conducted by the Grounded Electrical-Source Airborne Transient Electromagnetic (GREATEM) have shown that, this system is a promising method for modelling 3D resistivity structures in coastal areas. To expand the application of the GREATEM system in the future for studying hazardous wastes, sea water incursion and hydrocarbon exploration, a 3D-resistivity modelling that considers large lateral resistivity variations is required in case of large resistivity contrasts between land and sea in surveys of coastal areas where 1D resistivity model that assumes a horizontally layered structure might be inaccurate. In this abstract we present the preparation for developing a consistent three dimensional electromagnetic inversion algorithm to calculate the EM response over arbitrary 3D conductivity structure using GREATEM system. In forward modelling the second order partial differential equations for scalar and vector potential are discretized on a staggered-grid using the finite difference method (Fomenko and Mogi, 2002, Mogi et al., 2011). In the inversion method the 3D model discretized into a large number of rectangular cells of constant conductivity and the final solution is obtained by minimizing a global objective function composed of the model objective function and data misfit. To deal with a huge number of grids and wide range of frequencies in air borne data sets, a method for approximating sensitivities is introduced for the efficient 3-D inversion. Approximate sensitivities are derived by replacing adjoint secondary electric fields with those computed in the previous iteration. These sensitivities can reduce the computation time, without significant loss of accuracy when constructing a full sensitivity matrix for 3-D inversion, based on the Gauss-Newton method (N. Han et al., 2008). Now, we tested the algorithm in the frequency domain electromagnetic response of synthetic model considering a 3D conductor. Frequency-domain computation is executed
Performance improvement study of a relativistic magnetron using MAGIC-3D
Maurya, S.; Singh, V.V.P.; Jain, P.K.
2011-07-01
A three dimensional particle-in-cell (PlC) code, MAGIC3D, is used to examine the performance improvement in a relativistic magnetron by perturbing technique. Asymmetrical metal rods of different length have been used to perturb the magnetic field in the annular sector of the resonant system. Enhancement up to 45% in the radiated output power has been obtained in the perturbed magnetic field case over the unperturbed one. It has also been found in the simulation that oscillation start up time is reduced by 16 %, and the amplitude of the nearest competing mode goes down 9dB compared to unperturbed case. Perturbed magnetic field also reduces the end caps current improving the efficiency. (author)
Dynamic Analysis of 2D Electromagnetic Resonant Optical Scanner Using 3D Finite Element Method
NASA Astrophysics Data System (ADS)
Hirata, Katsuhiro; Hong, Sara; Maeda, Kengo
The optical scanner is a scanning device in which a laser beam is reflected by a mirror that can be rotated or oscillated. In this paper, we propose a new 2D electromagnetic resonant optical scanner that employs electromagnets and leaf springs. Torque characteristics and resonance characteristics of the scanner are analyzed using the 3D finite element method. The validity of the analysis is shown by comparing the characteristics inferred from the analysis with the characteristics of the prototype. Further, 2D resonance is investigated by introducing a superimposed-frequency current in a single coil.
Using the CAVE virtual-reality environment as an aid to 3-D electromagnetic field computation
Turner, L.R.; Levine, D.; Huang, M.; Papka, M; Kettunen, L.
1995-08-01
One of the major problems in three-dimensional (3-D) field computation is visualizing the resulting 3-D field distributions. A virtual-reality environment, such as the CAVE, (CAVE Automatic Virtual Environment) is helping to overcome this problem, thus making the results of computation more usable for designers and users of magnets and other electromagnetic devices. As a demonstration of the capabilities of the CAVE, the elliptical multipole wiggler (EMW), an insertion device being designed for the Advanced Photon Source (APS) now being commissioned at Argonne National Laboratory (ANL), wa made visible, along with its fields and beam orbits. Other uses of the CAVE in preprocessing and postprocessing computation for electromagnetic applications are also discussed.
Electromagnetic Response Inversion for a 3D Distribution of Conductivity/Dielect
Energy Science and Technology Software Center (ESTSC)
2001-10-24
NLCGCS inverts electromagnetic responses for a 3D distribution of electrical conductivity and dielectric permittivity within the earth for geophysical applications using single processor computers. The software comes bundled with a graphical user interface to aid in model construction and analysis and viewing of earth images. The solution employs both dipole and finite size source configurations for harmonic oscillatory sources. A new nonlinear preconditioner is included in the solution to speed up solution convergence.
Quantitative 3D electromagnetic field determination of 1D nanostructures from single projection.
Phatak, C; de Knoop, L; Houdellier, F; Gatel, C; Hÿtch, M J; Masseboeuf, A
2016-05-01
One-dimensional (1D) nanostructures have been regarded as the most promising building blocks for nanoelectronics and nanocomposite material systems as well as for alternative energy applications. Although they result in confinement of a material, their properties and interactions with other nanostructures are still very much three-dimensional (3D) in nature. In this work, we present a novel method for quantitative determination of the 3D electromagnetic fields in and around 1D nanostructures using a single electron wave phase image, thereby eliminating the cumbersome acquisition of tomographic data. Using symmetry arguments, we have reconstructed the 3D magnetic field of a nickel nanowire as well as the 3D electric field around a carbon nanotube field emitter, from one single projection. The accuracy of quantitative values determined here is shown to be a better fit to the physics at play than the value obtained by conventional analysis. Moreover the 3D reconstructions can then directly be visualized and used in the design of functional 3D architectures built using 1D nanostructures. PMID:26998702
Relativistic particle motion in nonuniform electromagnetic waves
NASA Technical Reports Server (NTRS)
Schmidt, G.; Wilcox, T.
1973-01-01
It is shown that a charged particle moving in a strong nonuniform electromagnetic wave suffers a net acceleration in the direction of the negative intensity gradient of the wave. Electrons will be expelled perpendicularly from narrow laser beams and various instabilities can result.
Relativistic particle motion in nonuniform electromagnetic waves
NASA Technical Reports Server (NTRS)
Schmidt, G.; Wilcox, T.
1973-01-01
A charged particle moving in a strong nonuniform electromagnetic wave which suffers a net acceleration in the direction of the negative intensity gradient of the wave was investigated. Electrons will be expelled perpendicularly from narrow laser beams and various instabilities result.
Models the Electromagnetic Response of a 3D Distribution using MP COMPUTERS
Energy Science and Technology Software Center (ESTSC)
1999-05-01
EM3D models the electromagnetic response of a 3D distribution of conductivity, dielectric permittivity and magnetic permeability within the earth for geophysical applications using massively parallel computers. The simulations are carried out in the frequency domain for either electric or magnetic sources for either scattered or total filed formulations of Maxwell''s equations. The solution is based on the method of finite differences and includes absorbing boundary conditions so that responses can be modeled up into themore » radar range where wave propagation is dominant. Recent upgrades in the software include the incorporation of finite size sources, that in addition to dipolar source fields, and a low induction number preconditioner that can significantly reduce computational run times. A graphical user interface (GUI) is bundled with the software so that complicated 3D models can be easily constructed and simulated with the software. The GUI also allows for plotting of the output.« less
Weakly relativistic electromagnetic solitons in warm plasmas
NASA Astrophysics Data System (ADS)
Sundar, Sita
2016-06-01
For slowly propagating electromagnetic solitons, validity of the cold plasma model is addressed using a more realistic model involving effects arising due to temperature as well as ion dynamics. Small amplitude single peak structures which are quasineutral are studied, and different regions of existence of bright and dark classes of solitons are delineated. Influence of temperature on spectral characteristics of the solitary structures is presented.
Pegoraro, F.
2009-11-10
A tutorial presentation is given describing the nature and the effects of relativistic nonlinearities in a plasma and indicating how they can be exploited in order to manipulate and enhance locally the energy density of the electromagnetic fields. The mathematical formulation and the examples presented are chosen from results available in the scientific literature.
Relativistic electromagnetic waves in an electron-ion plasma
NASA Technical Reports Server (NTRS)
Chian, Abraham C.-L.; Kennel, Charles F.
1987-01-01
High power laser beams can drive plasma particles to relativistic energies. An accurate description of strong waves requires the inclusion of ion dynamics in the analysis. The equations governing the propagation of relativistic electromagnetic waves in a cold electron-ion plasma can be reduced to two equations expressing conservation of energy-momentum of the system. The two conservation constants are functions of the plasma stream velocity, the wave velocity, the wave amplitude, and the electron-ion mass ratio. The dynamic parameter, expressing electron-ion momentum conversation in the laboratory frame, can be regarded as an adjustable quantity, a suitable choice of which will yield self-consistent solutions when other plasma parameters were specified. Circularly polarized electromagnetic waves and electrostatic plasma waves are used as illustrations.
Parametric instability of a relativistically strong electromagnetic wave.
NASA Technical Reports Server (NTRS)
Max, C. E.
1973-01-01
The stability of a circularly polarized electromagnetic wave that is strong enough to make plasma electrons, but not ions, relativistic is studied. Small perturbations are considered which propagate parallel to the large-amplitude driver. A relativistically strong wave can be unstable on time scales as short as twice its own oscillation period, and decays into a forward-going plasma oscillation and either one or two electromagnetic waves. Ion motion introduces an additional instability which can be important at short perturbation wavelengths, where the driver would otherwise be stable. The unstable ion and electron modes both have potential for producing anomalously large acceleration of relativistic particles, as well as significant amounts of backscattered light. These effects may be important in two applications: (1) the use of intense lasers to heat or compress plasma, and (2) the plasma surrounding a pulsar, if the pulsar is losing energy by radiation of electromagnetic waves at its rotation frequency. Instability persists in the nonrelativistic regime, reducing to stimulated Raman scattering as a special case.
Truesdell invariance in relativistic electromagnetic fields
NASA Astrophysics Data System (ADS)
Walwadkar, B. B.; Virkar, K. V.
1984-01-01
The Truesdell derivative of a contravariant tensor fieldX ab is defined with respect to a null congruencel a analogous to the Truesdell stress rate in classical continuum mechanics. The dynamical consequences of the Truesdell invariance with respect to a timelike vectoru a of the stress-energy tensor characterizing a charged perfect fluid with null conductivity are the conservation of pressure (p), charged density (e) an expansion-free flow, constancy of the Maxwell scalars, and vanishing spin coefficientsα+¯β = ¯σ - λ = τ = 0 (assuming freedom conditionsk = λ = ɛ ψ + ¯γ = 0). The electromagnetic energy momentum tensor for the special subcases of Ruse-Synge classification for typesA andB are described in terms of the spin coefficients introduced by Newman-Penrose.
3D electromagnetic modelling of a TTI medium and TTI effects in inversion
NASA Astrophysics Data System (ADS)
Jaysaval, Piyoosh; Shantsev, Daniil; de la Kethulle de Ryhove, Sébastien
2016-04-01
We present a numerical algorithm for 3D electromagnetic (EM) forward modelling in conducting media with general electric anisotropy. The algorithm is based on the finite-difference discretization of frequency-domain Maxwell's equations on a Lebedev grid, in which all components of the electric field are collocated but half a spatial step staggered with respect to the magnetic field components, which also are collocated. This leads to a system of linear equations that is solved using a stabilized biconjugate gradient method with a multigrid preconditioner. We validate the accuracy of the numerical results for layered and 3D tilted transverse isotropic (TTI) earth models representing typical scenarios used in the marine controlled-source EM method. It is then demonstrated that not taking into account the full anisotropy of the conductivity tensor can lead to misleading inversion results. For simulation data corresponding to a 3D model with a TTI anticlinal structure, a standard vertical transverse isotropic inversion is not able to image a resistor, while for a 3D model with a TTI synclinal structure the inversion produces a false resistive anomaly. If inversion uses the proposed forward solver that can handle TTI anisotropy, it produces resistivity images consistent with the true models.
Nguyen, B.T.; Hutchinson, S.A.
1995-07-01
The upwind leapfrog scheme for electromagnetic scattering is briefly described. Its application to the 3D Maxwell`s time domain equations is shown in detail. The scheme`s use of upwind characteristic variables and a narrow stencil result in a smaller demand in communication overhead, making it ideal for implementation on distributed memory parallel computers. The algorithm`s implementation on two message passing computers, a 1024-processor nCUBE 2 and a 1840-processor Intel Paragon, is described. Performance evaluation demonstrates that the scheme performs well with both good scaling qualities and high efficiencies on these machines.
3D modelling of the electromagnetic response of geophysical targets using the FDTD method
Debroux, P.S.
1996-05-01
A publicly available and maintained electromagnetic finite-difference time domain (FDTD) code has been applied to the forward modelling of the response of 1D, 2D and 3D geophysical targets to a vertical magnetic dipole excitation. The FDTD method is used to analyze target responses in the 1 MHz to 100 MHz range, where either conduction or displacement currents may have the controlling role. The response of the geophysical target to the excitation is presented as changes in the magnetic field ellipticity. The results of the FDTD code compare favorably with previously published integral equation solutions of the response of 1D targets, and FDTD models calculated with different finite-difference cell sizes are compared to find the effect of model discretization on the solution. The discretization errors, calculated as absolute error in ellipticity, are presented for the different ground geometry models considered, and are, for the most part, below 10% of the integral equation solutions. Finally, the FDTD code is used to calculate the magnetic ellipticity response of a 2D survey and a 3D sounding of complicated geophysical targets. The response of these 2D and 3D targets are too complicated to be verified with integral equation solutions, but show the proper low- and high-frequency responses.
Yi, Hee-Gyeong; Kang, Kyung Shin; Hong, Jung Min; Jang, Jinah; Park, Moon Nyeo; Jeong, Young Hun; Cho, Dong-Woo
2016-07-01
In cartilage tissue engineering, electromagnetic field (EMF) therapy has been reported to have a modest effect on promoting cartilage regeneration. However, these studies were conducted using different frequencies of EMF to stimulate chondrocytes. Thus, it is necessary to investigate the effect of EMF frequency on cartilage formation. In addition to the stimulation, a scaffold is required to satisfy the characteristics of cartilage such as its hydrated and dense extracellular matrix, and a mechanical resilience to applied loads. Therefore, we 3D-printed a composite construct composed of a polymeric framework and a chondrocyte-laden hydrogel. Here, we observed frequency-dependent positive and negative effects on chondrogenesis using a 3D cell-printed cartilage tissue. We found that a frequency of 45 Hz promoted gene expression and secretion of extracellular matrix molecules of chondrocytes. In contrast, a frequency of 7.5 Hz suppressed chondrogenic differentiation in vitro. Additionally, the EMF-treated composite constructs prior to implantation showed consistent results with those of in vitro, suggesting that in vitro pre-treatment with different EMF frequencies provides different capabilities for the enhancement of cartilage formation in vivo. This correlation between EMF frequency and 3D-printed chondrocytes suggests the necessity for optimization of EMF parameters when this physical stimulus is applied to engineered cartilage. © 2016 Wiley Periodicals, Inc. J Biomed Mater Res Part A: 104A: 1797-1804, 2016. PMID:26991030
Resonance of relativistic electrons with electromagnetic ion cyclotron waves
Denton, R. E.; Jordanova, V. K.; Bortnik, J.
2015-06-29
Relativistic electrons have been thought to more easily resonate with electromagnetic ion cyclotron EMIC waves if the total density is large. We show that, for a particular EMIC mode, this dependence is weak due to the dependence of the wave frequency and wave vector on the density. A significant increase in relativistic electron minimum resonant energy might occur for the H band EMIC mode only for small density, but no changes in parameters significantly decrease the minimum resonant energy from a nominal value. The minimum resonant energy depends most strongly on the thermal velocity associated with the field line motionmore » of the hot ring current protons that drive the instability. High density due to a plasmasphere or plasmaspheric plume could possibly lead to lower minimum resonance energy by causing the He band EMIC mode to be dominant. We demonstrate these points using parameters from a ring current simulation.« less
Resonance of relativistic electrons with electromagnetic ion cyclotron waves
Denton, R. E.; Jordanova, V. K.; Bortnik, J.
2015-06-29
Relativistic electrons have been thought to more easily resonate with electromagnetic ion cyclotron EMIC waves if the total density is large. We show that, for a particular EMIC mode, this dependence is weak due to the dependence of the wave frequency and wave vector on the density. A significant increase in relativistic electron minimum resonant energy might occur for the H band EMIC mode only for small density, but no changes in parameters significantly decrease the minimum resonant energy from a nominal value. The minimum resonant energy depends most strongly on the thermal velocity associated with the field line motion of the hot ring current protons that drive the instability. High density due to a plasmasphere or plasmaspheric plume could possibly lead to lower minimum resonance energy by causing the He band EMIC mode to be dominant. We demonstrate these points using parameters from a ring current simulation.
NASA Astrophysics Data System (ADS)
Ullmann, A.; Scheunert, M.; Afanasjew, M.; Börner, R.-U.; Siemon, B.; Spitzer, K.
2016-07-01
As a standard procedure, multi-frequency helicopter-borne electromagnetic (HEM) data are inverted to conductivity-depth models using 1-D inversion methods, which may, however, fail in areas of strong lateral conductivity contrasts (so-called induction anomalies). Such areas require more realistic multi-dimensional modelling. Since the full 3-D inversion of an entire HEM data set is still extremely time consuming, our idea is to combine fast 1-D and accurate but numerically expensive 3-D inversion of HEM data in such a way that the full 3-D inversion is only carried out for those parts of a HEM survey which are affected by induction anomalies. For all other parts, a 1-D inversion method is sufficient. We present a newly developed algorithm for identification, selection, and extraction of induction anomalies in HEM data sets and show how the 3-D inversion model of the anomalous area is re-integrated into the quasi-1-D background. Our proposed method is demonstrated to work properly on a synthetic and a field HEM data set from the Cuxhaven tunnel valley in Germany. We show that our 1-D/3-D approach yields better results compared to 1-D inversions in areas where 3-D effects occur.
Electromagnetic fluctuations and normal modes of a drifting relativistic plasma
Ruyer, C.; Gremillet, L.; Bénisti, D.; Bonnaud, G.
2013-11-15
We present an exact calculation of the power spectrum of the electromagnetic fluctuations in a relativistic equilibrium plasma described by Maxwell-Jüttner distribution functions. We consider the cases of wave vectors parallel or normal to the plasma mean velocity. The relative contributions of the subluminal and supraluminal fluctuations are evaluated. Analytical expressions of the spatial fluctuation spectra are derived in each case. These theoretical results are compared to particle-in-cell simulations, showing a good reproduction of the subluminal fluctuation spectra.
Charged and Electromagnetic Fields from Relativistic Quantum Geometry
NASA Astrophysics Data System (ADS)
Arcodía, Marcos; Bellini, Mauricio
2016-06-01
In the Relativistic Quantum Geometry (RQG) formalism recently introduced, was explored the possibility that the variation of the tensor metric can be done in a Weylian integrable manifold using a geometric displacement, from a Riemannian to a Weylian integrable manifold, described by the dynamics of an auxiliary geometrical scalar field $\\theta$, in order that the Einstein tensor (and the Einstein equations) can be represented on a Weyl-like manifold. In this framework we study jointly the dynamics of electromagnetic fields produced by quantum complex vector fields, which describes charges without charges. We demonstrate that complex fields act as a source of tetra-vector fields which describe an extended Maxwell dynamics.
NASA Technical Reports Server (NTRS)
Nishikawa, K.-I.; Mizuno, Y.; Hardee, P.; Hededal, C. B.; Fishman, G. J.
2006-01-01
Recent PIC simulations using injected relativistic electron-ion (electro-positron) jets into ambient plasmas show that acceleration occurs in relativistic shocks. The Weibel instability created in shocks is responsible for particle acceleration, and generation and amplification of highly inhomogeneous, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection in relativistic jets. The "jitter" radiation from deflected electrons has different properties than the synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understand the complex time evolution and spectral structure in relativistic jets and gamma-ray bursts. We will present recent PIC simulations which show particle acceleration and magnetic field generation. We will also calculate associated self-consistent emission from relativistic shocks.
Large-scale 3D inversion of frequency domain controlled-source electromagnetic data
NASA Astrophysics Data System (ADS)
Miller, C. R.; Routh, P. S.; Donaldson, P.; Oldenburg, D. W.
2005-05-01
these targets. In this paper we will present a strategy to invert a large-scale EM data set. Haber, E., Ascher, U. and Oldenburg, D., 2004, Inversion of 3D electromagnetic data in frequency and time domain using an inexact all-at-once approach: Geophysics, Soc. of Expl. Geophys., 69, 1216-1228.
NASA Astrophysics Data System (ADS)
Boivin, A.; Lai, P.; Samson, C.; Cloutis, E.; Holladay, S.; Monteiro Santos, F. A.
2013-07-01
The Mars Methane Analogue Mission simulates a micro-rover mission whose purpose is to detect, analyze, and determine the source of methane emissions on the planet's surface. As part of this project, both an electromagnetic induction sounder (EMIS) and a high-resolution triangulation-based 3D laser scanner were tested at the Jeffrey open-pit asbestos mine to identify and characterize geological environments favourable to the occurrence of methane. The presence of serpentinite in the form of chrysotile (asbestos), magnesium carbonate, and iron oxyhydroxides make the mine a likely location for methane production. The EMIS clearly delineated the contacts between the two geological units found at the mine, peridotite and slate, which are separated by a shear zone. Both the peridotite and slate units have low and uniform apparent electrical conductivity and magnetic susceptibility, while the shear zone has much higher conductivity and susceptibility, with greater variability. The EMIS data were inverted and the resulting model captured lateral conductivity variations through the different bedrock geological units buried beneath a gravel road. The 3D point cloud data acquired by the laser scanner were fitted with triangular meshes where steeply dipping triangles were plotted in dark grey to accentuate discontinuities. The resulting images were further processed using Sobel edge detection to highlight networks of fractures which are potential pathways for methane seepage.
NASA Astrophysics Data System (ADS)
Rezzolla, Luciano; Ahmedov, Bobomurat J.
2016-07-01
An important issue in the asteroseismology of compact and magnetized stars is the determination of the dissipation mechanism which is most efficient in damping the oscillations when these are produced. In a linear regime and for low-multipolarity modes, these mechanisms are confined to either gravitational-wave or electromagnetic losses. We here consider the latter and compute the energy losses in the form of Poynting fluxes, Joule heating and Ohmic dissipation in a relativistic oscillating spherical star with a dipolar magnetic field in vacuum. While this approach is not particularly realistic for rapidly rotating stars, it has the advantage that it is fully analytic and that it provides expressions for the electric and magnetic fields produced by the most common modes of oscillation both in the vicinity of the star and far away from it. In this way, we revisit and extend to a relativistic context the classical estimates of McDermott et al. Overall, we find that general-relativistic corrections lead to electromagnetic damping time-scales that are at least one order of magnitude smaller than in Newtonian gravity. Furthermore, with the only exception of g (gravity) modes, we find that f (fundamental), p (pressure), i (interface) and s (shear) modes are suppressed more efficiently by gravitational losses than by electromagnetic ones.
External Electromagnetic Fields of Slowly Rotating Relativistic Magnetized NUT Stars
NASA Astrophysics Data System (ADS)
Ahmedov, B. J.; Khugaev, A. V.
2006-08-01
Analytic general relativistic expressions for the electromagnetic fields external to a slowly-rotating magnetized NUT star with non-vanishing gravitomagnetic charge have been presented. Solutions for the electric and magnetic fields have been found after separating the Maxwell equations in the external background spacetime of a slowly rotating NUT star into angular and radial parts in the lowest order approximation in specific angular momentum and NUT parameter . The relativistic star is considered isolated and in vacuum, with different models for stellar magnetic field: i) monopolar magnetic field and ii) dipolar magnetic field aligned with the axis of rotation. It has been shown that the general relativistic corrections due to the dragging of reference frames and gravitomagnetic charge are not present in the form of the magnetic fields but emerge only in the form of the electric fields. In particular, it has been obtained that the frame-dragging and gravitomagnetic charge provide an additional induced electric field which is analogous to the one introduced by the rotation of the star in the flat spacetime limit.
Parallel goal-oriented adaptive finite element modeling for 3D electromagnetic exploration
NASA Astrophysics Data System (ADS)
Zhang, Y.; Key, K.; Ovall, J.; Holst, M.
2014-12-01
We present a parallel goal-oriented adaptive finite element method for accurate and efficient electromagnetic (EM) modeling of complex 3D structures. An unstructured tetrahedral mesh allows this approach to accommodate arbitrarily complex 3D conductivity variations and a priori known boundaries. The total electric field is approximated by the lowest order linear curl-conforming shape functions and the discretized finite element equations are solved by a sparse LU factorization. Accuracy of the finite element solution is achieved through adaptive mesh refinement that is performed iteratively until the solution converges to the desired accuracy tolerance. Refinement is guided by a goal-oriented error estimator that uses a dual-weighted residual method to optimize the mesh for accurate EM responses at the locations of the EM receivers. As a result, the mesh refinement is highly efficient since it only targets the elements where the inaccuracy of the solution corrupts the response at the possibly distant locations of the EM receivers. We compare the accuracy and efficiency of two approaches for estimating the primary residual error required at the core of this method: one uses local element and inter-element residuals and the other relies on solving a global residual system using a hierarchical basis. For computational efficiency our method follows the Bank-Holst algorithm for parallelization, where solutions are computed in subdomains of the original model. To resolve the load-balancing problem, this approach applies a spectral bisection method to divide the entire model into subdomains that have approximately equal error and the same number of receivers. The finite element solutions are then computed in parallel with each subdomain carrying out goal-oriented adaptive mesh refinement independently. We validate the newly developed algorithm by comparison with controlled-source EM solutions for 1D layered models and with 2D results from our earlier 2D goal oriented
NASA Astrophysics Data System (ADS)
Yang, Dikun; Oldenburg, Douglas W.; Haber, Eldad
2014-03-01
Airborne electromagnetic (AEM) methods are highly efficient tools for assessing the Earth's conductivity structures in a large area at low cost. However, the configuration of AEM measurements, which typically have widely distributed transmitter-receiver pairs, makes the rigorous modelling and interpretation extremely time-consuming in 3-D. Excessive overcomputing can occur when working on a large mesh covering the entire survey area and inverting all soundings in the data set. We propose two improvements. The first is to use a locally optimized mesh for each AEM sounding for the forward modelling and calculation of sensitivity. This dedicated local mesh is small with fine cells near the sounding location and coarse cells far away in accordance with EM diffusion and the geometric decay of the signals. Once the forward problem is solved on the local meshes, the sensitivity for the inversion on the global mesh is available through quick interpolation. Using local meshes for AEM forward modelling avoids unnecessary computing on fine cells on a global mesh that are far away from the sounding location. Since local meshes are highly independent, the forward modelling can be efficiently parallelized over an array of processors. The second improvement is random and dynamic down-sampling of the soundings. Each inversion iteration only uses a random subset of the soundings, and the subset is reselected for every iteration. The number of soundings in the random subset, determined by an adaptive algorithm, is tied to the degree of model regularization. This minimizes the overcomputing caused by working with redundant soundings. Our methods are compared against conventional methods and tested with a synthetic example. We also invert a field data set that was previously considered to be too large to be practically inverted in 3-D. These examples show that our methodology can dramatically reduce the processing time of 3-D inversion to a practical level without losing resolution
Scattering and coupling effects of electromagnetic waves in 3D networks of spheres
NASA Astrophysics Data System (ADS)
Defos Du Rau, M.; Pessan, F.; Ruffie, G.; Vignéras-Lefebvre, V.; Parneix, J. P.
1998-01-01
In this paper, the problem of electromagnetic scattering from a 3D system of spheres is considered and an iterative solution that accounts for multiple scattering is proposed. The Mie formalism used for a single sphere is extended to account for multiple scattered fields between several particles. The translational addition theorems for spherical wave functions are used to express the electromagnetic field scattered by a sphere S_i in terms of an incident field for a sphere S_k in a spherical coordinates system attached to the sphere S_k. In this work, the numerical convergence of the method is discussed and associated computational times are given. Numerical computations including Radar Cross Section (RCS) and radiation patterns for various 3D configurations are presented. Some of them are compared with free-space measurements made in the 8 to 100 GHz frequency band using vectorial network analyzers. 11.55.-m S-matrix theory; analytic structure of amplitudes Cet article étudie la diffusion des ondes électromagnétiques par des réseaux tridimensionnels de sphères et propose une méthode itérative pour prendre en compte les effets de multidiffusion. Le formalisme de Mie utilisé dans le cas d'une sphère est étendu pour calculer les champs "multidiffusés" entre plusieurs particules. Les théorèmes d'addition et de translation des fonctions d'onde sphériques sont utilisés pour exprimer le champ diffusé par une sphère S_i comme étant incident sur une sphère S_k, dans un système de coordonnées sphériques lié au centre de S_k. La convergence numérique de la méthode est discutée et des temps de calcul sont donnés. Des résultats numériques tels que des Surfaces Équivalentes Radar (SER) et des diagrammes de rayonnement pour différentes configurations tridimensionnelles sont montrés. Certains d'entre eux sont comparés à des mesures en espace libre faites à l'aide d'analyseurs de réseaux vectoriels dans la bande de fréquence 8{-}100 GHz.
NASA Astrophysics Data System (ADS)
Tandon, K.; Egbert, G.; Siripunvaraporn, W.
2003-12-01
We are developing a modular system for three-dimensional inversion of electromagnetic (EM) induction data, using an object oriented programming approach. This approach allows us to modify the individual components of the inversion scheme proposed, and also reuse the components for variety of problems in earth science computing howsoever diverse they might be. In particular, the modularity allows us to (a) change modeling codes independently of inversion algorithm details; (b) experiment with new inversion algorithms; and (c) modify the way prior information is imposed in the inversion to test competing hypothesis and techniques required to solve an earth science problem. Our initial code development is for EM induction equations on a staggered grid, using iterative solution techniques in 3D. An example illustrated here is an experiment with the sensitivity of 3D magnetotelluric inversion to uncertainties in the boundary conditions required for regional induction problems. These boundary conditions should reflect the large-scale geoelectric structure of the study area, which is usually poorly constrained. In general for inversion of MT data, one fixes boundary conditions at the edge of the model domain, and adjusts the earth?s conductivity structure within the modeling domain. Allowing for errors in specification of the open boundary values is simple in principle, but no existing inversion codes that we are aware of have this feature. Adding a feature such as this is straightforward within the context of the modular approach. More generally, a modular approach provides an efficient methodology for setting up earth science computing problems to test various ideas. As a concrete illustration relevant to EM induction problems, we investigate the sensitivity of MT data near San Andreas Fault at Parkfield (California) to uncertainties in the regional geoelectric structure.
Parametric decay of an extraordinary electromagnetic wave in relativistic plasma
Dorofeenko, V. G.; Krasovitskiy, V. B.; Turikov, V. A.
2015-03-15
Parametric instability of an extraordinary electromagnetic wave in plasma preheated to a relativistic temperature is considered. A set of self-similar nonlinear differential equations taking into account the electron “thermal” mass is derived and investigated. Small perturbations of the parameters of the heated plasma are analyzed in the linear approximation by using the dispersion relation determining the phase velocities of the fast and slow extraordinary waves. In contrast to cold plasma, the evanescence zone in the frequency range above the electron upper hybrid frequency vanishes and the asymptotes of both branches converge. Theoretical analysis of the set of nonlinear equations shows that the growth rate of decay instability increases with increasing initial temperature of plasma electrons. This result is qualitatively confirmed by numerical simulations of plasma heating by a laser pulse injected from vacuum.
RELATIVISTIC MERGERS OF SUPERMASSIVE BLACK HOLES AND THEIR ELECTROMAGNETIC SIGNATURES
Bode, Tanja; Haas, Roland; Laguna, Pablo; Shoemaker, Deirdre; Bogdanovic, Tamara
2010-06-01
Coincident detections of electromagnetic (EM) and gravitational wave (GW) signatures from coalescence events of supermassive black holes (SMBHs) are the next observational grand challenge. Such detections will provide the means to study cosmological evolution and accretion processes associated with these gargantuan compact objects. More generally, the observations will enable testing general relativity in the strong, nonlinear regime and will provide independent cosmological measurements to high precision. Understanding the conditions under which coincidences of EM and GW signatures arise during SMBH mergers is therefore of paramount importance. As an essential step toward this goal, we present results from the first fully general relativistic, hydrodynamical study of the late inspiral and merger of equal-mass, spinning SMBH binaries in a gas cloud. We find that variable EM signatures correlated with GWs can arise in merging systems as a consequence of shocks and accretion combined with the effect of relativistic beaming. The most striking EM variability is observed for systems where spins are aligned with the orbital axis and where orbiting black holes form a stable set of density wakes, but all systems exhibit some characteristic signatures that can be utilized in searches for EM counterparts. In the case of the most massive binaries observable by the Laser Interferometer Space Antenna, calculated luminosities imply that they may be identified by EM searches to z {approx} 1, while lower mass systems and binaries immersed in low density ambient gas can only be detected in the local universe.
A non-linear 3D printed electromagnetic vibration energy harvester
NASA Astrophysics Data System (ADS)
Constantinou, P.; Roy, S.
2015-12-01
This paper describes a novel electromagnetic energy harvester that exploits the low flexural modulus of ABS and comprises of a nonlinear mechanism to enhance the generated power and bandwidth. The device is printed using desktop additive manufacturing techniques (3D printing) that use thermoplastics. It has a ‘V’ spring topology and exhibits a softening spring non-linearity introduced through the magnetic arrangement, which introduces a monostable potential well. A model is presented and measurements correspond favourably. The produced prototype generates a peak power of approximately 2.5mW at a frame acceleration of 1g and has a power bandwidth of approximately 1.2→1.5Hz and 3.5→3.9Hz during up and down sweeps respectively. The device has a power density of 0.4mW/cm3 at a frame acceleration of 1g and a density of 0.04mW/cm3 from a generated power of 25μW at 0.1g.
Propagation of Electromagnetic Waves in 3D Opal-based Magnetophotonic Crystals
NASA Astrophysics Data System (ADS)
Pardavi-Horvath, Martha; Makeeva, Galina S.; Golovanov, Oleg A.; Rinkevich, Anatolii B.
2013-03-01
Opals, a class of self-organized 3D nanostructures, are typical representatives of photonic bandgap structures. The voids inside of the opal structure of close packed SiO2 spheres can be infiltrated by a magnetic material, creating magnetically tunable magnetophotonic crystals with interesting and potentially useful properties at GHz and THz frequencies. The propagation of electromagnetic waves at microwave frequencies was investigated numerically in SiO2 opal based magnetic nanostructures, using rigorous mathematical models to solve Maxwell's equations complemented by the Landau-Lifshitz equation with electrodynamic boundary conditions. The numerical approach is based on Galerkin's projection method using the decomposition algorithm on autonomous blocks with Floquet channels. The opal structure consists of SiO2 nanospheres, with inter-sphere voids infiltrated with nanoparticles of Ni-Zn ferrites. Both the opal matrix and the ferrite are assumed to be lossy. A model, taking into account the real structure of the ferrite particles in the opal's voids was developed to simulate the measured FMR lineshape of the ferrite infiltrated opal. The numerical technique shows an excellent agreement when applied to model recent experimental data on similar ferrite opals.
FDFD: A 3D Finite-Difference Frequency-Domain Code for Electromagnetic Induction Tomography
NASA Astrophysics Data System (ADS)
Champagne, Nathan J.; Berryman, James G.; Buettner, H. Michael
2001-07-01
A new 3D code for electromagnetic induction tomography with intended applications to environmental imaging problems has been developed. The approach consists of calculating the fields within a volume using an implicit finite-difference frequency-domain formulation. The volume is terminated by an anisotropic perfectly matched layer region that simulates an infinite domain by absorbing outgoing waves. Extensive validation of this code has been done using analytical and semianalytical results from other codes, and some of those results are presented in this paper. The new code is written in Fortran 90 and is designed to be easily parallelized. Finally, an adjoint field method of data inversion, developed in parallel for solving the fully nonlinear inverse problem for electrical conductivity imaging (e.g., for mapping underground conducting plumes), uses this code to provide solvers for both forward and adjoint fields. Results obtained from this inversion method for high-contrast media are encouraging and provide a significant improvement over those obtained from linearized inversion methods.
Towards 3D and Multilayer Electromagnetic Metamaterials Structures in the THz Range
NASA Astrophysics Data System (ADS)
Casse, B. D. F.; Moser, H. O.; Lee, J. W.; Inglis, S.; Bahou, M.; Jian, L. K.
2007-03-01
V. G. Veselago predicted that left-handed materials would exhibit a plethora of unusual effects such as a negative index of refraction as used in Snell's law, a reverse Doppler and Cerenkov effect. This novel class of materials, following Pendry's recipes, can potentially restore evanescent waves to focus subwavelength details in an image. Micron-size electromagnetic metamaterials (EM^3) structures which exist so far were produced and characterized as single layer structures. Furthermore the structures were produced with a low yield. In the first part of the talk, we will present techniques for producing copious amount of EM^3 chips via the LIGA process using Synchrotron radiation and demonstrate assembly of the first multilayer THz EM^3 structures. The planar micro- or nanoEM^3 produced so far are also highly anisotropic. Recently, we proposed schemes to produce more isotropic structures, within the same matrix, via tilted X-ray exposures that were introduced in the LIGA process years ago. In the second part of the talk, we will show the results of microfabrication of nearly 3D EM^3 structures for the THz range.
NASA Technical Reports Server (NTRS)
Nishikawa, K.-I.
2006-01-01
Nonthermal radiation observed from astrophysical systems containing (relativistic) jets and shocks, e.g., supernova remnants, active galactic nuclei (AGNs), gamma-ray bursts (GRBs), and Galactic microquasar systems usually have power-law emission spectra. Fermi acceleration is the mechanism usually assumed for the acceleration of particles in astrophysical environments. Recent PIC simulations using injected relativistic electron-ion (electro-positron) jets show that acceleration occurs within the downstream jet, rather than by the scattering of particles back and forth across the shock as in Fermi acceleration. Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, other two-streaming instability, and the Weibel instability) created in the .shocks are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The "jitter" radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants. We will review recent PIC simulations which show particle acceleration in jets.
Finite-element 3D simulation tools for high-current relativistic electron beams
NASA Astrophysics Data System (ADS)
Humphries, Stanley; Ekdahl, Carl
2002-08-01
The DARHT second-axis injector is a challenge for computer simulations. Electrons are subject to strong beam-generated forces. The fields are fully three-dimensional and accurate calculations at surfaces are critical. We describe methods applied in OmniTrak, a 3D finite-element code suite that can address DARHT and the full range of charged-particle devices. The system handles mesh generation, electrostatics, magnetostatics and self-consistent particle orbits. The MetaMesh program generates meshes of conformal hexahedrons to fit any user geometry. The code has the unique ability to create structured conformal meshes with cubic logic. Organized meshes offer advantages in speed and memory utilization in the orbit and field solutions. OmniTrak is a versatile charged-particle code that handles 3D electric and magnetic field solutions on independent meshes. The program can update both 3D field solutions from the calculated beam space-charge and current-density. We shall describe numerical methods for orbit tracking on a hexahedron mesh. Topics include: 1) identification of elements along the particle trajectory, 2) fast searches and adaptive field calculations, 3) interpolation methods to terminate orbits on material surfaces, 4) automatic particle generation on multiple emission surfaces to model space-charge-limited emission and field emission, 5) flexible Child law algorithms, 6) implementation of the dual potential model for 3D magnetostatics, and 7) assignment of charge and current from model particle orbits for self-consistent fields.
Finite-Difference Algorithm for Simulating 3D Electromagnetic Wavefields in Conductive Media
NASA Astrophysics Data System (ADS)
Aldridge, D. F.; Bartel, L. C.; Knox, H. A.
2013-12-01
Electromagnetic (EM) wavefields are routinely used in geophysical exploration for detection and characterization of subsurface geological formations of economic interest. Recorded EM signals depend strongly on the current conductivity of geologic media. Hence, they are particularly useful for inferring fluid content of saturated porous bodies. In order to enhance understanding of field-recorded data, we are developing a numerical algorithm for simulating three-dimensional (3D) EM wave propagation and diffusion in heterogeneous conductive materials. Maxwell's equations are combined with isotropic constitutive relations to obtain a set of six, coupled, first-order partial differential equations governing the electric and magnetic vectors. An advantage of this system is that it does not contain spatial derivatives of the three medium parameters electric permittivity, magnetic permeability, and current conductivity. Numerical solution methodology consists of explicit, time-domain finite-differencing on a 3D staggered rectangular grid. Temporal and spatial FD operators have order 2 and N, where N is user-selectable. We use an artificially-large electric permittivity to maximize the FD timestep, and thus reduce execution time. For the low frequencies typically used in geophysical exploration, accuracy is not unduly compromised. Grid boundary reflections are mitigated via convolutional perfectly matched layers (C-PMLs) imposed at the six grid flanks. A shared-memory-parallel code implementation via OpenMP directives enables rapid algorithm execution on a multi-thread computational platform. Good agreement is obtained in comparisons of numerically-generated data with reference solutions. EM wavefields are sourced via point current density and magnetic dipole vectors. Spatially-extended inductive sources (current carrying wire loops) are under development. We are particularly interested in accurate representation of high-conductivity sub-grid-scale features that are common
Robust and scalable 3-D geo-electromagnetic modelling approach using the finite element method
NASA Astrophysics Data System (ADS)
Grayver, Alexander V.; Bürg, Markus
2014-07-01
We present a robust and scalable solver for time-harmonic Maxwell's equations for problems with large conductivity contrasts, wide range of frequencies, stretched grids and locally refined meshes. The solver is part of the fully distributed adaptive 3-D electromagnetic modelling scheme which employs the finite element method and unstructured non-conforming hexahedral meshes for spatial discretization using the open-source software deal.II. We use the complex-valued electric field formulation and split it into two real-valued equations for which we utilize an optimal block-diagonal pre-conditioner. Application of this pre-conditioner requires the solution of two smaller real-valued symmetric problems. We solve them by using either a direct solver or the conjugate gradient method pre-conditioned with the recently introduced auxiliary space technique. The auxiliary space pre-conditioner reformulates the original problem in form of several simpler ones, which are then solved using highly efficient algebraic multigrid methods. In this paper, we consider the magnetotelluric case and verify our numerical scheme by using COMMEMI 3-D models. Afterwards, we run a series of numerical experiments and demonstrate that the solver converges in a small number of iterations for a wide frequency range and variable problem sizes. The number of iterations is independent of the problem size, but exhibits a mild dependency on frequency. To test the stability of the method on locally refined meshes, we have implemented a residual-based a posteriori error estimator and compared it with uniform mesh refinement for problems up to 200 million unknowns. We test the scalability of the most time consuming parts of our code and show that they fulfill the strong scaling assumption as long as each MPI process possesses enough degrees of freedom to alleviate communication overburden. Finally, we refer back to a direct solver-based pre-conditioner and analyse its complexity in time. The results show
Poulin, Eric; Racine, Emmanuel; Beaulieu, Luc; Binnekamp, Dirk
2015-03-15
Purpose: In high dose rate brachytherapy (HDR-B), current catheter reconstruction protocols are relatively slow and error prone. The purpose of this technical note is to evaluate the accuracy and the robustness of an electromagnetic (EM) tracking system for automated and real-time catheter reconstruction. Methods: For this preclinical study, a total of ten catheters were inserted in gelatin phantoms with different trajectories. Catheters were reconstructed using a 18G biopsy needle, used as an EM stylet and equipped with a miniaturized sensor, and the second generation Aurora{sup ®} Planar Field Generator from Northern Digital Inc. The Aurora EM system provides position and orientation value with precisions of 0.7 mm and 0.2°, respectively. Phantoms were also scanned using a μCT (GE Healthcare) and Philips Big Bore clinical computed tomography (CT) system with a spatial resolution of 89 μm and 2 mm, respectively. Reconstructions using the EM stylet were compared to μCT and CT. To assess the robustness of the EM reconstruction, five catheters were reconstructed twice and compared. Results: Reconstruction time for one catheter was 10 s, leading to a total reconstruction time inferior to 3 min for a typical 17-catheter implant. When compared to the μCT, the mean EM tip identification error was 0.69 ± 0.29 mm while the CT error was 1.08 ± 0.67 mm. The mean 3D distance error was found to be 0.66 ± 0.33 mm and 1.08 ± 0.72 mm for the EM and CT, respectively. EM 3D catheter trajectories were found to be more accurate. A maximum difference of less than 0.6 mm was found between successive EM reconstructions. Conclusions: The EM reconstruction was found to be more accurate and precise than the conventional methods used for catheter reconstruction in HDR-B. This approach can be applied to any type of catheters and applicators.
Relativistic Laser Pulse Intensification with 3D Printed Micro-Tube Plasma Target
NASA Astrophysics Data System (ADS)
Ji, Liangliang; Snyder, Joseph; Pukhov, Alexander; Akli, Kramer
2015-11-01
The potential and applications of laser-plasma interactions (LPI) are restricted by the parameter space of existing lasers and targets. Advancing the laser intensity to the extreme regime is motivated by the production of energetic particle beams and by the quest to explore the exotic regimes of light-matter interaction. Target density and dimensions can always be varied to optimize the outcome. Here, we propose to create another degree of freedom in the parameter space of LPI using recent advances in 3D printing of materials. Fine structures at nm scale with high repetition and accuracy can nowadays be manufactured, allowing for a full precise control of the target. We demonstrate, via particle-in-cell (PIC) simulations, that 3D-printed micro-tube plasma (MTP) targets yield an intensity enhancement factor of 2-5. The novel MTP targets not only act as a plasma optical device to reach the 1023W/cm2 threshold based on today's intensities, but can also boost the generation of secondary particle and radiation sources. This work demonstrates that the combination of high contrast high power lasers and nano-3D printing techniques opens new paths in the intensity frontier and LPI micro-engineering.
Electromagnetic lepton-pair production in relativistic collisions
NASA Astrophysics Data System (ADS)
Albert, C. J.; Ernst, D. J.; Strayer, M. R.; Bottcher, C.
Electromagnetic lepton-pair production in relativistic collisions is studied in an ab initio approach with no free parameters. After a semi-classical approximation to the relative motion of the two incident particles is made, the resulting second-order diagram is calculated using a Monte Carlo technique to evaluate the resulting seven-dimensional integral. We examine the case of electron-positron pair production in pi(-)p collisions at p(sub pi) = 17 GeV. We find that a significant fraction of the measured pairs in this reaction are produced via the magnetic spin-flip current of the proton. Approaches, such as the equivalent photon approximation, which neglect this part of the current predict much too small a cross section. This feature is traced to the cuts imposed in taking the experimental data. Lepton-pair production in the scattering of He-3, He-4 and He-4, He-4 is proposed as a clean way of experimentally separating the spin-flip and non-flip processes; predictions are made for these systems.
Chenot, Jean-Loup; Bay, Francois
2007-04-07
The different stages of metal forming processes often involve - beyond the mechanical deformations processes - other physical coupled problems, such as heat transfer, electromagnetism or metallurgy. The purpose of this paper is to focus on problems involving electromagnetic couplings. After a brief recall on electromagnetic modeling, we shall then focus on induction heating processes and present some results regarding heat transfer, as well as mechanical couplings. A case showing coupling for metallurgic microstructure evolution will conclude this paper.
NASA Technical Reports Server (NTRS)
Mizuno, Yosuke; Lyubarsky, Yuri; ishikawa, Ken-Ichi; Hardee, Philip E.
2010-01-01
We have investigated the development of current-driven (CD) kink instability through three-dimensional relativistic MHD simulations. A static force-free equilibrium helical magnetic configuration is considered in order to study the influence of the initial configuration on the linear and nonlinear evolution of the instability. We found that the initial configuration is strongly distorted but not disrupted by the kink instability. The instability develops as predicted by linear theory. In the non-linear regime the kink amplitude continues to increase up to the terminal simulation time, albeit at different rates, for all but one simulation. The growth rate and nonlinear evolution of the CD kink instability depends moderately on the density profile and strongly on the magnetic pitch profile. The growth rate of the kink mode is reduced in the linear regime by an increase in the magnetic pitch with radius and the non-linear regime is reached at a later time than for constant helical pitch. On the other hand, the growth rate of the kink mode is increased in the linear regime by a decrease in the magnetic pitch with radius and reaches the non-linear regime sooner than the case with constant magnetic pitch. Kink amplitude growth in the non-linear regime for decreasing magnetic pitch leads to a slender helically twisted column wrapped by magnetic field. On the other hand, kink amplitude growth in the non-linear regime nearly ceases for increasing magnetic pitch.
NASA Astrophysics Data System (ADS)
Li, Ji-Jiao; Li, Bo; Peng, Qin-Mei; Zhou, Ji; Li, Long-Tu
2014-09-01
Three groups of three-dimensional (3D) TiO2 woodpile electromagnetic gap materials with tailed rheological properties were developed for direct-written fabrication. Appropriate amount of polyethyleneimine (PEI) dispersants allow the preparation of TiO2 inks with a high solid content of 42 vol.%, which enables them to flow through the nozzles easily. The inks exhibit pseudoplastic behavior. The measured microwave characteristics of the results agree well with simulations based on plane wave expansion (PWE).
Xie, G.; Li, J.; Majer, E.; Zuo, D.
1998-07-01
This paper describes a new 3D parallel GILD electromagnetic (EM) modeling and nonlinear inversion algorithm. The algorithm consists of: (a) a new magnetic integral equation instead of the electric integral equation to solve the electromagnetic forward modeling and inverse problem; (b) a collocation finite element method for solving the magnetic integral and a Galerkin finite element method for the magnetic differential equations; (c) a nonlinear regularizing optimization method to make the inversion stable and of high resolution; and (d) a new parallel 3D modeling and inversion using a global integral and local differential domain decomposition technique (GILD). The new 3D nonlinear electromagnetic inversion has been tested with synthetic data and field data. The authors obtained very good imaging for the synthetic data and reasonable subsurface EM imaging for the field data. The parallel algorithm has high parallel efficiency over 90% and can be a parallel solver for elliptic, parabolic, and hyperbolic modeling and inversion. The parallel GILD algorithm can be extended to develop a high resolution and large scale seismic and hydrology modeling and inversion in the massively parallel computer.
NASA Astrophysics Data System (ADS)
Zhdanov, M. S.; Cuma, M.; Black, N.; Wilson, G. A.
2009-12-01
The marine controlled source electromagnetic (MCSEM) method has become widely used in offshore oil and gas exploration. Interpretation of MCSEM data is still a very challenging problem, especially if one would like to take into account the realistic 3D structure of the subsurface. The inversion of MCSEM data is complicated by the fact that the EM response of a hydrocarbon-bearing reservoir is very weak in comparison with the background EM fields generated by an electric dipole transmitter in complex geoelectrical structures formed by a conductive sea-water layer and the terranes beneath it. In this paper, we present a review of the recent developments in the area of large-scale 3D EM forward modeling and inversion. Our approach is based on using a new integral form of Maxwell’s equations allowing for an inhomogeneous background conductivity, which results in a numerically effective integral representation for 3D EM field. This representation provides an efficient tool for the solution of 3D EM inverse problems. To obtain a robust inverse model of the conductivity distribution, we apply regularization based on a focusing stabilizing functional which allows for the recovery of models with both smooth and sharp geoelectrical boundaries. The method is implemented in a fully parallel computer code, which makes it possible to run large-scale 3D inversions on grids with millions of inversion cells. This new technique can be effectively used for active EM detection and monitoring of the subsurface targets.
3D Airborne Electromagnetic Inversion: A case study from the Musgrave Region, South Australia
NASA Astrophysics Data System (ADS)
Cox, L. H.; Wilson, G. A.; Zhdanov, M. S.; Sunwall, D. A.
2012-12-01
Geophysicists know and accept that geology is inherently 3D, and is resultant from complex, overlapping processes related to genesis, metamorphism, deformation, alteration, weathering, and/or hydrogeology. Yet, the geophysics community has long relied on qualitative analysis, conductivity depth imaging (CDIs), 1D inversion, and/or plate modeling. There are many reasons for this deficiency, not the least of which has been the lack of capacity for historic 3D AEM inversion algorithms to invert entire surveys so as to practically affect exploration decisions. Our recent introduction of a moving sensitivity domain (footprint) methodology has been a paradigm shift in AEM interpretation. The basis of this method is that one needs only to calculate the responses and sensitivities for that part of the 3D earth model that is within the AEM system's sensitivity domain (footprint), and then superimpose all sensitivity domains into a single, sparse sensitivity matrix for the entire 3D earth model which is then updated in a regularized inversion scheme. This has made it practical to rigorously invert entire surveys with thousands of line kilometers of AEM data to mega-cell 3D models in hours using multi-processor workstations. Since 2010, over eighty individual projects have been completed for Aerodat, AEROTEM, DIGHEM, GEOTEM, HELITEM, HoisTEM, MEGATEM, RepTEM, RESOLVE, SkyTEM, SPECTREM, TEMPEST, and VTEM data from Australia, Brazil, Canada, Finland, Ghana, Peru, Tanzania, the US, and Zambia. Examples of 3D AEM inversion have been published for a variety of applications, including mineral exploration, oil sands exploration, salinity, permafrost, and bathymetry mapping. In this paper, we present a comparison of 3D inversions for SkyTEM, SPECTREM, TEMPET and VTEM data acquired over the same area in the Musgrave region of South Australia for exploration under cover.
3D relaxation MHD modeling with FOI-PERFECT code for electromagnetically driven HED systems
NASA Astrophysics Data System (ADS)
Wang, Ganghua; Duan, Shuchao; Xie, Weiping; Kan, Mingxian; Institute of Fluid Physics Collaboration
2015-11-01
One of the challenges in numerical simulations of electromagnetically driven high energy density (HED) systems is the existence of vacuum region. The electromagnetic part of the conventional model adopts the magnetic diffusion approximation (magnetic induction model). The vacuum region is approximated by artificially increasing the resistivity. On one hand the phase/group velocity is superluminal and hence non-physical in the vacuum region, on the other hand a diffusion equation with large diffusion coefficient can only be solved by implicit scheme. Implicit method is usually difficult to parallelize and converge. A better alternative is to solve the full electromagnetic equations for the electromagnetic part. Maxwell's equations coupled with the constitutive equation, generalized Ohm's law, constitute a relaxation model. The dispersion relation is given to show its transition from electromagnetic propagation in vacuum to resistive MHD in plasma in a natural way. The phase and group velocities are finite for this system. A better time stepping is adopted to give a 3rd full order convergence in time domain without the stiff relaxation term restriction. Therefore it is convenient for explicit & parallel computations. Some numerical results of FOI-PERFECT code are also given. Project supported by the National Natural Science Foundation of China (Grant No. 11172277,11205145).
NASA Astrophysics Data System (ADS)
Vidal, A.; San-Blas, A. A.; Quesada-Pereira, F. D.; Pérez-Soler, J.; Gil, J.; Vicente, C.; Gimeno, B.; Boria, V. E.
2015-07-01
A novel technique for the full-wave analysis of 3-D complex waveguide devices is presented. This new formulation, based on the Boundary Integral-Resonant Mode Expansion (BI-RME) method, allows the rigorous full-wave electromagnetic characterization of 3-D arbitrarily shaped metallic structures making use of extremely low CPU resources (both time and memory). The unknown electric current density on the surface of the metallic elements is represented by means of Rao-Wilton-Glisson basis functions, and an algebraic procedure based on a singular value decomposition is applied to transform such functions into the classical solenoidal and nonsolenoidal basis functions needed by the original BI-RME technique. The developed tool also provides an accurate computation of the electromagnetic fields at an arbitrary observation point of the considered device, so it can be used for predicting high-power breakdown phenomena. In order to validate the accuracy and efficiency of this novel approach, several new designs of band-pass waveguides filters are presented. The obtained results (S-parameters and electromagnetic fields) are successfully compared both to experimental data and to numerical simulations provided by a commercial software based on the finite element technique. The results obtained show that the new technique is specially suitable for the efficient full-wave analysis of complex waveguide devices considering an integrated coaxial excitation, where the coaxial probes may be in contact with the metallic insets of the component.
NASA Astrophysics Data System (ADS)
Gansen, A.; El Hachemi, M.; Belouettar, S.; Hassan, O.; Morgan, K.
2015-12-01
In computational electromagnetics, the advantages of the standard Yee algorithm are its simplicity and its low computational costs. However, because of the accuracy losses resulting from the staircased representation of curved interfaces, it is normally not the method of choice for modelling electromagnetic interactions with objects of arbitrary shape. For these problems, an unstructured mesh finite volume time domain method is often employed, although the scheme does not satisfy the divergence free condition at the discrete level. In this paper, we generalize the standard Yee algorithm for use on unstructured meshes and solve the problem concerning the loss of accuracy linked to staircasing, while preserving the divergence free nature of the algorithm. The scheme is implemented on high quality primal Delaunay and dual Voronoi meshes. The performance of the approach was validated in previous work by simulating the scattering of electromagnetic waves by spherical 3D PEC objects in free space. In this paper we demonstrate the performance of this scheme for penetration problems in lossy dielectrics using a new averaging technique for Delaunay and Voronoi edges at the interface. A detailed explanation of the implementation of the method, and a demonstration of the quality of the results obtained for transmittance and scattering simulations by 3D objects of arbitrary shapes, are presented.
Modeling 3-D objects with planar surfaces for prediction of electromagnetic scattering
NASA Technical Reports Server (NTRS)
Koch, M. B.; Beck, F. B.; Cockrell, C. R.
1992-01-01
Electromagnetic scattering analysis of objects at resonance is difficult because low frequency techniques are slow and computer intensive, and high frequency techniques may not be reliable. A new technique for predicting the electromagnetic backscatter from electrically conducting objects at resonance is studied. This technique is based on modeling three dimensional objects as a combination of flat plates where some of the plates are blocking the scattering from others. A cube is analyzed as a simple example. The preliminary results compare well with the Geometrical Theory of Diffraction and with measured data.
NASA Astrophysics Data System (ADS)
Wong, Kenneth H.; Tang, Jonathan; Zhang, Hui J.; Varghese, Emmanuel; Cleary, Kevin R.
2005-04-01
An effective treatment method for organs that move with respiration (such as the lungs, pancreas, and liver) is a major goal of radiation medicine. In order to treat such tumors, we need (1) real-time knowledge of the current location of the tumor, and (2) the ability to adapt the radiation delivery system to follow this constantly changing location. In this study, we used electromagnetic tracking in a swine model to address the first challenge, and to determine if movement of a marker attached to the skin could accurately predict movement of an internal marker embedded in an organ. Under approved animal research protocols, an electromagnetically tracked needle was inserted into a swine liver and an electromagnetically tracked guidewire was taped to the abdominal skin of the animal. The Aurora (Northern Digital Inc., Waterloo, Canada) electromagnetic tracking system was then used to monitor the position of both of these sensors every 40 msec. Position readouts from the sensors were then tested to see if any of the movements showed correlation. The strongest correlations were observed between external anterior-posterior motion and internal inferior-superior motion, with many other axes exhibiting only weak correlation. We also used these data to build a predictive model of internal motion by taking segments from the data and using them to derive a general functional relationship between the internal needle and the external guidewire. For the axis with the strongest correlation, this model enabled us to predict internal organ motion to within 1 mm.
NASA Astrophysics Data System (ADS)
Kim, B.; Byun, J.; Seol, S. J.; Jeong, S.; Chung, Y.; Kwon, T.
2015-12-01
For many decades, gas hydrates have been received great attention as a potential source of natural gas. Therefore, the detailed information of structures of buried gas hydrates and their concentrations are prerequisite for the production for the gas hydrate as a reliable source of alternate energy. Recently, for this reason, a lot of gas hydrate assessment methods have been proposed by many researchers. However, it is still necessary to establish as new method for the further improvement of the accuracy of the 3D gas hydrate distribution. In this study, we present a 3D joint inversion method that provides superior quantitative information of gas hydrate distributions using 3D seismic data obtained by ocean-bottom cable (OBC) and marine controlled-source electromagnetic (CSEM) data. To verify our inversion method, we first built the general 3D gas hydrate model containing vertical methane-flow pathways. With the described model, we generated synthetic 3D OBC data and marine CSEM data using finite element modeling algorithms, respectively. In the joint inversion process, to obtain the high-resolution volumetric P-wave velocity structure, we applied the 3D full waveform inversion algorithm to the acquired OBC data. After that, the obtained P-wave velocity model is used as the structure constraint to compute cross-gradients with the updated resistivity model in the EM inversion process. Finally, petrophysical relations were applied to estimate volumetric gas hydrate concentrations. The proposed joint inversion process makes possible to obtain more precise quantitative gas hydrate assessment than inversion processes using only seismic or EM data. This technique can be helpful for accurate decision-making in gas hydrate development as well as in their production monitoring.
NASA Astrophysics Data System (ADS)
Pare, Pascal; Gribenko, Alexander V.; Cox, Leif H.; Čuma, Martin; Wilson, Glenn A.; Zhdanov, Michael S.; Legault, Jean; Smit, Jaco; Polome, Louis
2012-04-01
Geological, geochemical, and geophysical surveys have been conducted in the area of the Pebble Cu-Au-Mo porphyry deposit in south-west Alaska since 1985. This case study compares three-dimensional (3D) inversion results from Anglo American's proprietary SPECTREM 2000 fixed-wing time-domain airborne electromagnetic (AEM) and Geotech's ZTEM airborne audio-frequency magnetics (AFMAG) systems flown over the Pebble deposit. Within the commonality of their physics, 3D inversions of both SPECTREM and ZTEM recover conductivity models consistent with each other and the known geology. Both 3D inversions recover conductors coincident with alteration associated with both Pebble East and Pebble West. The high grade CuEqn 0.6% ore shell is not consistently following the high conductive trend, suggesting that the SPECTREM and ZTEM responses correspond in part to the sulphide distribution, but not directly with the ore mineralization. As in any exploration project, interpretation of both surveys has yielded an improved understanding of the geology, alteration and mineralization of the Pebble system and this will serve well for on-going exploration activities. There are distinct practical advantages to the use of both SPECTREM and ZTEM, so we draw no recommendation for either system. We can conclude however, that 3D inversion of both AEM and ZTEM surveys is now a practical consideration and that it has added value to exploration at Pebble.
Xiong, Z.; Tripp, A.C.
1994-12-31
This paper presents an integral equation algorithm for 3D EM modeling at high frequencies for applications in engineering an environmental studies. The integral equation method remains the same for low and high frequencies, but the dominant roles of the displacements currents complicate both numerical treatments and interpretations. With singularity extraction technique they successively extended the application of the Hankel filtering technique to the computation of Hankel integrals occurring in high frequency EM modeling. Time domain results are calculated from frequency domain results via Fourier transforms. While frequency domain data are not obvious for interpretations, time domain data show wave-like pictures that resemble seismograms. Both 1D and 3D numerical results show clearly the layer interfaces.
NASA Astrophysics Data System (ADS)
Abubakirov, É. B.; Denisenko, A. N.; Konyushkov, A. P.; Soluyanov, E. I.; Yastrebov, V. V.
2014-10-01
We study operation of a relativistic backward-wave oscillator driven by an external electromagnetic signal. Such operation regimes as hard excitation of self-oscillations and amplification of the external signal are implemented experimentally. The conditions for possible synchronization of the relativistic backward-wave oscillator by an external signal is discussed. The possibility of accelerating the onset of oscillations by the action of an external signal is confirmed experimentally. The conditions of realization of the amplification regime and the main effects, which limit the amplification coefficient, are determined. The obtained results can be used to optimize the parameters of generators and amplifiers based on relativistic backward-wave oscillators.
Electromagnetic Scattering of Finite and Infinite 3D Lattices in Polarizable Backgrounds
Gallinet, Benjamin; Martin, Olivier J. F.
2009-10-07
A novel method is elaborated for the electromagnetic scattering from periodical arrays of scatterers embedded in a polarizable background. A dyadic periodic Green's function is introduced to calculate the scattered electric field in a lattice of dielectric or metallic objects. The method exhibits strong advantages: discretization and computation of the field are restricted to the volume of the scatterers in the unit cell, open and periodic boundary conditions for the electric field are included in the Green's tensor, and finally both near and far-fields physics are directly revealed, without any additional computational effort. Promising applications include the design of periodic structures such as frequency-selective surfaces, photonic crystals and metamaterials.
Electromagnetic and absorbing property of CIPs/resin composite using the 3D forming process
NASA Astrophysics Data System (ADS)
Xu, Yonggang; Liang, Zichang; Wang, Xiaobing; Yuan, Liming; Li, Xinghao
2016-08-01
The absorbing composite filled with the flaky carbonyl iron particles (CIPs) were prepared using a three-dimensional (3D) forming process, in which the forming powder was fabricated using a milling process. The surface morphology was characterized by the scanning electron microscopy, the static magnetic property was evaluated on a vibrating sample magnetometer, and X-ray diffraction (XRD) patterns were done to analyze the particle crystal grain structure. The complex permittivity and permeability were measured using a vector network analyzer in the frequency range of 4-18 GHz. With the variable thickness was set, the reflection loss (RL) was simulated to analyze the absorbing property of the composite. The results showed that the forming powder was uniformly dispersed in the absorber, and the saturation magnetization and the grain structure of the CIPs in the forming powder nearly did not change in the milling process. With the same volume content CIPs added, the average permittivity and the imaginary permeability of the samples added the powder was smaller than the directly mixing sample due to the aggregation effect. The RL results showed that the absorbing composites using the 3D forming process with thickness 6 or 8 mm had an better absorbing property (minimum RL -13.58 and -21.85 dB) in 4-18 GHz.
Liu, Jia; Cao, Mao-Sheng; Luo, Qiang; Shi, Hong-Long; Wang, Wen-Zhong; Yuan, Jie
2016-08-31
We fabricated the nickel chains by a facile wet chemical method. The morphology of nickel chains were tailored by adjusting the amount of PVP during the synthesis process. Both the complex permittivity and permeability of the three-dimensional (3D) nets constructed by nickel chains present strong dependences on temperature in the frequency range of 8.2-12.4 GHz and temperature range of 323-573 K. The peaks in imaginary component of permittivity and permeability mainly derive from interfacial polarizations and resonances, devoting to dielectric and magnetic loss, respectively. The effect from both dielectric and magnetism contribute to enhancing the microwave absorption. The maximum absorption value of the 3D nickel chain nets is approximately -50 dB at 8.8 GHz and 373 K with a thickness of 1.8 mm, and the bandwidth less than -10 dB almost covers the whole investigated frequency band. These are encouraging findings, which provide the potential advantages of magnetic transition metal-based materials for microwave absorption application at elevated temperature. PMID:27509241
Relativistic quark model for the Omega- electromagnetic form factors
G. Ramalho, K. Tsushima, Franz Gross
2009-08-01
We compute the Omega- electromagnetic form factors and the decuplet baryon magnetic moments using a quark model application of the Covariant Spectator Theory. Our predictions for the Omega- electromagnetic form factors can be tested in the future by lattice QCD simulations at the physical strange quark mass.
NASA Astrophysics Data System (ADS)
Xu, Fujun; Yao, Lan; Zhao, Da; Jiang, Muwen; Qiu, Yipping
2013-10-01
A three-dimensionally integrated microstrip antenna (3DIMA) is a microstrip antenna woven into the three-dimensional woven composite for load bearing while functioning as an antenna. In this study, the effect of weaving direction of conductive yarns on electromagnetic performance of 3DIMAs are investigated by designing, simulating and experimental testing of two microstrip antennas with different weaving directions of conductive yarns: one has the conductive yarns along the antenna feeding direction (3DIMA-Exp1) and the other has the conductive yarns perpendicular the antenna feeding direction (3DIMA-Exp2). The measured voltage standing wave ratio (VSWR) of 3DIMA-Exp1 was 1.4 at the resonant frequencies of 1.39 GHz; while that of 3DIMA-Exp2 was 1.2 at the resonant frequencies of 1.35 GHz. In addition, the measured radiation pattern of the 3DIMA-Exp1 has smaller back lobe and higher gain value than those of the 3DIMA-Exp2. This result indicates that the waving direction of conductive yarns may have a significant impact on electromagnetic performance of textile structural antennas.
Advanced 3D electromagnetic and particle-in-cell modeling on structured/unstructured hybrid grids
Seidel, D.B.; Pasik, M.F.; Kiefer, M.L.; Riley, D.J.; Turner, C.D.
1998-01-01
New techniques have been recently developed that allow unstructured, free meshes to be embedded into standard 3-dimensional, rectilinear, finite-difference time-domain grids. The resulting hybrid-grid modeling capability allows the higher resolution and fidelity of modeling afforded by free meshes to be combined with the simplicity and efficiency of rectilinear techniques. Integration of these new methods into the full-featured, general-purpose QUICKSILVER electromagnetic, Particle-In-Cell (PIC) code provides new modeling capability for a wide variety of electromagnetic and plasma physics problems. To completely exploit the integration of this technology into QUICKSILVER for applications requiring the self-consistent treatment of charged particles, this project has extended existing PIC methods for operation on these hybrid unstructured/rectilinear meshes. Several technical issues had to be addressed in order to accomplish this goal, including the location of particles on the unstructured mesh, adequate conservation of charge, and the proper handling of particles in the transition region between structured and unstructured portions of the hybrid grid.
NASA Astrophysics Data System (ADS)
Woodbury, D.; Kubota, S.; Johnson, I.
2014-10-01
Computer simulations of electromagnetic wave propagation in magnetized plasmas are an important tool for both plasma heating and diagnostics. For active millimeter-wave and microwave diagnostics, accurately modeling the evolution of the beam parameters for launched, reflected or scattered waves in a toroidal plasma requires that calculations be done using the full 3-D geometry. Previously, we reported on the application of GPGPU (General-Purpose computing on Graphics Processing Units) to a 3-D vacuum Maxwell code using the FDTD (Finite-Difference Time-Domain) method. Tests were done for Gaussian beam propagation with a hard source antenna, utilizing the parallel processing capabilities of the NVIDIA K20M. In the current study, we have modified the 3-D code to include a soft source antenna and an induced current density based on the cold plasma approximation. Results from Gaussian beam propagation in an inhomogeneous anisotropic plasma, along with comparisons to ray- and beam-tracing calculations will be presented. Additional enhancements, such as advanced coding techniques for improved speedup, will also be investigated. Supported by U.S. DoE Grant DE-FG02-99-ER54527 and in part by the U.S. DoE, Office of Science, WDTS under the Science Undergraduate Laboratory Internship program.
NASA Astrophysics Data System (ADS)
Bai, Xianchen; Yang, Jianhua; Zhang, Jiande
2012-08-01
By using an electromagnetic particle-in-cell (PIC) code, an S-band two-cavity wide-gap klystron amplifier (WKA) loaded with washers/rods structure is designed and investigated for high power injection application. Influences of the washers/rods structure on the high frequency characteristics and the basic operation of the amplifier are presented. Generally, the rod structure has great impacts on the space-charge potential depression and the resonant frequency of the cavities. Nevertheless, if only the resonant frequency is tuned to the desired operation frequency, effects of the rod size on the basic operation of the amplifier are expected to be very weak. The 3-dimension (3-D) PIC simulation results show an output power of 0.98 GW corresponding to an efficiency of 33% for the WKA, with a 594 keV, 5 kA electron beam guided by an external magnetic field of 1.5 Tesla. Moreover, if a conductive plane is placed near the output gap, such as the electron collector, the beam potential energy can be further released, and the RF power can be increased to about 1.07 GW with the conversion efficiency of about 36%.
Bai Xianchen; Yang Jianhua; Zhang Jiande
2012-08-15
By using an electromagnetic particle-in-cell (PIC) code, an S-band two-cavity wide-gap klystron amplifier (WKA) loaded with washers/rods structure is designed and investigated for high power injection application. Influences of the washers/rods structure on the high frequency characteristics and the basic operation of the amplifier are presented. Generally, the rod structure has great impacts on the space-charge potential depression and the resonant frequency of the cavities. Nevertheless, if only the resonant frequency is tuned to the desired operation frequency, effects of the rod size on the basic operation of the amplifier are expected to be very weak. The 3-dimension (3-D) PIC simulation results show an output power of 0.98 GW corresponding to an efficiency of 33% for the WKA, with a 594 keV, 5 kA electron beam guided by an external magnetic field of 1.5 Tesla. Moreover, if a conductive plane is placed near the output gap, such as the electron collector, the beam potential energy can be further released, and the RF power can be increased to about 1.07 GW with the conversion efficiency of about 36%.
Hewett, D.W.; DiPeso, G.; Gibbons, M; Lambert, M.; Tung, L.S.
1997-02-03
Low-frequency electromagnetic simulation models have a wide range of industrial applications. We have built several models, differentiated by slightly different physics approximations or computational solution methods, that have proven quite useful in a variety of applications. Our models been used to investigate beam plasma interactions in ICF targets, antenna plasma coupling in plasma processing, and magnetic implosion drive in Z-pinch pulsed power generators. The common feature of these models is that they retain inductive effects but implicitly ignore computationally intensive, fully electromagnetic effects. However, the preponderance of our work has been limited to only two dimensions. We have made significant progress modeling low-frequency electromagnetic physics with a new model in 2-D that is now capable of modeling antenna structures in 3-D. Although LLNL`s interest in plasma processing has diminished, we have certainly added to LLNL`s capabilities. Interestingly, we have already found another application, the magnetic behavior of read/write heads in the magnetic storage industry, that can make use of many of the computational methods described here, rewarding us again for maintaining a strong core competency in low-frequency EM plasmas.
Fall, Mandiaye; Boutami, Salim; Glière, Alain; Stout, Brian; Hazart, Jerome
2013-06-01
A combination of the multilevel fast multipole method (MLFMM) and boundary element method (BEM) can solve large scale photonics problems of arbitrary geometry. Here, MLFMM-BEM algorithm based on a scalar and vector potential formulation, instead of the more conventional electric and magnetic field formulations, is described. The method can deal with multiple lossy or lossless dielectric objects of arbitrary geometry, be they nested, in contact, or dispersed. Several examples are used to demonstrate that this method is able to efficiently handle 3D photonic scatterers involving large numbers of unknowns. Absorption, scattering, and extinction efficiencies of gold nanoparticle spheres, calculated by the MLFMM, are compared with Mie's theory. MLFMM calculations of the bistatic radar cross section (RCS) of a gold sphere near the plasmon resonance and of a silica coated gold sphere are also compared with Mie theory predictions. Finally, the bistatic RCS of a nanoparticle gold-silver heterodimer calculated with MLFMM is compared with unmodified BEM calculations. PMID:24323115
3D electromagnetic optimization of the front face of the ITER ICRF antenna
NASA Astrophysics Data System (ADS)
Louche, F.; Dumortier, P.; Messiaen, A.; Durodié, F.
2011-10-01
In the framework of the ion cyclotron resonance heating (ICRH) antenna development for ITER, a design based on an external matching concept has been proposed [1]. We present in this work a series of electromagnetic simulations of this design performed with the commercial code CST Microwave Studio [2]. On the one hand, we explore how various geometrical modifications of some parts of the antenna (the straps and the four-port junction) can practically double the RF power coupled to the plasma. This optimization is supported by transmission line analysis. On the other hand, we treat the important question of the opportunity to tilt the straps in the toroidal direction to follow the plasma curvature as close as possible. We show that a configuration with two toroidal segments is sufficient and that further segmentation is not necessary. This work also underlines significant progress in the realism of ICRH antenna modelling and the importance of considering realistic load shaping in the models.
The 3-D reconstruction of medieval wetland reclamation through electromagnetic induction survey.
De Smedt, Philippe; Van Meirvenne, Marc; Herremans, Davy; De Reu, Jeroen; Saey, Timothy; Meerschman, Eef; Crombé, Philippe; De Clercq, Wim
2013-01-01
Studies of past human-landscape interactions rely upon the integration of archaeological, biological and geological information within their geographical context. However, detecting the often ephemeral traces of human activities at a landscape scale remains difficult with conventional archaeological field survey. Geophysical methods offer a solution by bridging the gap between point finds and the surrounding landscape, but these surveys often solely target archaeological features. Here we show how simultaneous mapping of multiple physical soil properties with a high resolution multi-receiver electromagnetic induction (EMI) survey permits a reconstruction of the three-dimensional layout and pedological setting of a medieval reclaimed landscape in Flanders (Belgium). Combined with limited and directed excavations, the results offer a unique insight into the way such marginal landscapes were reclaimed and occupied during the Middle Ages. This approach provides a robust foundation for unravelling complex historical landscapes and will enhance our understanding of past human-landscape interactions. PMID:23519060
The 3-D reconstruction of medieval wetland reclamation through electromagnetic induction survey
De Smedt, Philippe; Van Meirvenne, Marc; Herremans, Davy; De Reu, Jeroen; Saey, Timothy; Meerschman, Eef; Crombé, Philippe; De Clercq, Wim
2013-01-01
Studies of past human-landscape interactions rely upon the integration of archaeological, biological and geological information within their geographical context. However, detecting the often ephemeral traces of human activities at a landscape scale remains difficult with conventional archaeological field survey. Geophysical methods offer a solution by bridging the gap between point finds and the surrounding landscape, but these surveys often solely target archaeological features. Here we show how simultaneous mapping of multiple physical soil properties with a high resolution multi-receiver electromagnetic induction (EMI) survey permits a reconstruction of the three-dimensional layout and pedological setting of a medieval reclaimed landscape in Flanders (Belgium). Combined with limited and directed excavations, the results offer a unique insight into the way such marginal landscapes were reclaimed and occupied during the Middle Ages. This approach provides a robust foundation for unravelling complex historical landscapes and will enhance our understanding of past human-landscape interactions. PMID:23519060
Morgan, F. Dale; Sogade, John
2004-12-14
This project was designed as a broad foundational study of spectral induced polarization (SIP) for characterization of contaminated sites. It encompassed laboratory studies of the effects of chemistry on induced polarization, development of 3D forward modeling and inversion codes, and investigations of inductive and capacitive coupling problems. In the laboratory part of the project a physico-chemical model developed in this project was used to invert laboratory IP spectra for the grain size and the effective grain size distribution of the sedimentary rocks as well as the formation factor, porosity, specific surface area, and the apparent fractal dimension. Furthermore, it was established that the IP response changed with the solution chemistry, the concentration of a given solution chemistry, valence of the constituent ions, and ionic radius. In the field part of the project, a 3D complex forward and inverse model was developed. It was used to process data acquired at two frequencies (1/16 Hz and 1/ 4Hz) in a cross-borehole configuration at the A-14 outfall area of the Savannah River Site (SRS) during March 2003 and June 2004. The chosen SRS site was contaminated with Tetrachloroethylene (TCE) and Trichloroethylene (PCE) that were disposed in this area for several decades till the 1980s. The imaginary conductivity produced from the inverted 2003 data correlated very well with the log10 (PCE) concentration derived from point sampling at 1 ft spacing in five ground-truth boreholes drilled after the data acquisition. The equivalent result for the 2004 data revealed that there were significant contaminant movements during the period March 2003 and June 2004, probably related to ground-truth activities and nearby remediation activities. Therefore SIP was successfully used to develop conceptual models of volume distributions of PCE/TCE contamination. In addition, the project developed non-polarizing electrodes that can be deployed in boreholes for years. A total of 28
Octet Baryon Electromagnetic Form Factors in a Relativistic Quark Model
Gilberto Ramalho, Kazuo Tsushima
2011-09-01
We study the octet baryon electromagnetic properties by applying the covariant spectator quark model, and provide covariant parametrization that can be used to study baryon electromagnetic reactions. While we use the lattice QCD data in the large pion mass regime (small pion cloud effects) to determine the parameters of the model in the valence quark sector, we use the nucleon physical and octet baryon magnetic moment data to parameterize the pion cloud contributions. The valence quark contributions for the octet baryon electromagnetic form factors are estimated by extrapolating the lattice parametrization in the large pion mass regime to the physical regime. As for the pion cloud contributions, we parameterize them in a covariant, phenomenological manner, combined with SU(3) symmetry. We also discuss the impact of the pion cloud effects on the octet baryon electromagnetic form factors and their radii.
NASA Technical Reports Server (NTRS)
Wu, Jianshi; Derrickson, J. H.; Parnell, T. A.; Strayer, M. R.
1999-01-01
We study the screening effects of the atomic electrons in the electromagnetic production of electron-positron pairs in relativistic nucleus-atom collisions for fixed target experiments. Our results are contrasted with those obtained in bare collisions, with particular attention given to its dependence on the beam energy and the target atom.
Efficient computational methods for electromagnetic imaging with applications to 3D magnetotellurics
NASA Astrophysics Data System (ADS)
Kordy, Michal Adam
The motivation for this work is the forward and inverse problem for magnetotellurics, a frequency domain electromagnetic remote-sensing geophysical method used in mineral, geothermal, and groundwater exploration. The dissertation consists of four papers. In the first paper, we prove the existence and uniqueness of a representation of any vector field in H(curl) by a vector lying in H(curl) and H(div). It allows us to represent electric or magnetic fields by another vector field, for which nodal finite element approximation may be used in the case of non-constant electromagnetic properties. With this approach, the system matrix does not become ill-posed for low-frequency. In the second paper, we consider hexahedral finite element approximation of an electric field for the magnetotelluric forward problem. The near-null space of the system matrix for low frequencies makes the numerical solution unstable in the air. We show that the proper solution may obtained by applying a correction on the null space of the curl. It is done by solving a Poisson equation using discrete Helmholtz decomposition. We parallelize the forward code on multicore workstation with large RAM. In the next paper, we use the forward code in the inversion. Regularization of the inversion is done by using the second norm of the logarithm of conductivity. The data space Gauss-Newton approach allows for significant savings in memory and computational time. We show the efficiency of the method by considering a number of synthetic inversions and we apply it to real data collected in Cascade Mountains. The last paper considers a cross-frequency interpolation of the forward response as well as the Jacobian. We consider Pade approximation through model order reduction and rational Krylov subspace. The interpolating frequencies are chosen adaptively in order to minimize the maximum error of interpolation. Two error indicator functions are compared. We prove a theorem of almost always lucky failure in the
NASA Astrophysics Data System (ADS)
Smirnov, M. Yu.; Korja, T.; Pedersen, L. B.
2009-04-01
Two electromagnetic arrays are used in the EMMA project to study conductivity structure of the Archaean lithosphere in the Fennoscandian Shield. The first array was operated during almost one year, while the second one was running only during the summer time. Twelve 5-components magnetotelluric instruments with fluxgate magnetometers recorded simultaneously time variations of Earth's natural electromagnetic field at the sites separated by c. 30 km. To better control the source field and to obtain galvanic distortion free responses we have applied horizontal spatial gradient (HSG) technique to the data. The study area is highly inhomogeneous, thus classical HSG might give erroneous results. The method was extended to include anomalous field effects by implementing multivariate analysis. The HSG transfer functions were then used to control static shift distortions of apparent resistivities. During the BEAR experiment 1997-2002, the conductance map of entire Fennoscandia was assembled and finally converted into 3D volume resistivity model. We have used the model, refined it to get denser grid around measurement area and calculated MT transfer functions after 3D modeling. We have used trial-and-error method in order to further improve the model. The data set was also inverted using 3D code of Siripunvaraporn (2005). In the first stage we have used homogeneous halfspace as starting model for the inversion. In the next step we have used final 3D forward model as apriori model. The usage of apriori information significantly stabilizes the inverse solution, especially in case of a limited amount of data available. The results show that in the Archaean Domain a conductive layer is found in the upper/middle crust on contrary to previous results from other regions of the Archaean crust in the Fennoscandian Shield. Data also suggest enhanced conductivity at the depth of c. 100 km. Conductivity below the depth of 200-250 km is lower than that of the laboratory based estimates
Modelling of 3-D electromagnetic responses using the time-wavenumber method
Lee, S.
1991-12-01
The diffusion of electromagnetic fields in time and the three spatial dimensions can be modelled using a new numerical algorithm that is tailored for geophysical applications. The novel feature of the algorithm is that a large part of the computation is done in the wavenumber domain. Here, the spatial Fourier transforms of the vertical magnetic field and the vertical current density are used to define two scalar potentials. For either a vertical electric or a vertical magnetic dipole source at the subsurface these wavenumber potentials can be represented by a simple Gaussian distribution function. In the air, the fields satisfy the Laplace equation. The flow of this algorithm is as follows: the potentials are defined in the wavenumber domain as an initial condition depending on the source configuration, the vector current density J in space is obtained from the potentials using the inverse Fourier transform, the vector electric field E is obtained by multiplying J by resistivity, the updated potentials are then obtained from the forward Fourier transform of E. Using the updated potential as a subsequent initial condition these steps are repeated until the solution reaches the final time. Since spatial derivatives can be exactly evaluated in the wavenumber domain by simple multiplications, this algorithm requires far less memory than the conventional finite difference (FD) method. The conventional FD method needs finer discretization in space in order to minimize the numerical dispersion caused by numerical differentiation in space. The conductivity distribution for this algorithm is piece-wise continuous and bounded in the wavenumber domain.
Aydmer, A.A.; Chew, W.C.; Cui, T.J.; Wright, D.L.; Smith, D.V.; Abraham, J.D.
2001-01-01
A simple and efficient method for large scale three-dimensional (3-D) subsurface imaging of inhomogeneous background is presented. One-dimensional (1-D) multifrequency distorted Born iterative method (DBIM) is employed in the inversion. Simulation results utilizing synthetic scattering data are given. Calibration of the very early time electromagnetic (VETEM) experimental waveforms is detailed along with major problems encountered in practice and their solutions. This discussion is followed by the results of a large scale application of the method to the experimental data provided by the VETEM system of the U.S. Geological Survey. The method is shown to have a computational complexity that is promising for on-site inversion.
Brosten, T.R.; Day-Lewis, F. D.; Schultz, G.M.; Curtis, G.P.; Lane, J.W.
2011-01-01
Electromagnetic induction (EMI) instruments provide rapid, noninvasive, and spatially dense data for characterization of soil and groundwater properties. Data from multi-frequency EMI tools can be inverted to provide quantitative electrical conductivity estimates as a function of depth. In this study, multi-frequency EMI data collected across an abandoned uranium mill site near Naturita, Colorado, USA, are inverted to produce vertical distribution of electrical conductivity (EC) across the site. The relation between measured apparent electrical conductivity (ECa) and hydraulic conductivity (K) is weak (correlation coefficient of 0.20), whereas the correlation between the depth dependent EC obtained from the inversions, and K is sufficiently strong to be used for hydrologic estimation (correlation coefficient of -0.62). Depth-specific EC values were correlated with co-located K measurements to develop a site-specific ln(EC)-ln(K) relation. This petrophysical relation was applied to produce a spatially detailed map of K across the study area. A synthetic example based on ECa values at the site was used to assess model resolution and correlation loss given variations in depth and/or measurement error. Results from synthetic modeling indicate that optimum correlation with K occurs at ~0.5m followed by a gradual correlation loss of 90% at 2.3m. These results are consistent with an analysis of depth of investigation (DOI) given the range of frequencies, transmitter-receiver separation, and measurement errors for the field data. DOIs were estimated at 2.0??0.5m depending on the soil conductivities. A 4-layer model, with varying thicknesses, was used to invert the ECa to maximize available information within the aquifer region for improved correlations with K. Results show improved correlation between K and the corresponding inverted EC at similar depths, underscoring the importance of inversion in using multi-frequency EMI data for hydrologic estimation. ?? 2011.
NASA Astrophysics Data System (ADS)
Brosten, Troy R.; Day-Lewis, Frederick D.; Schultz, Gregory M.; Curtis, Gary P.; Lane, John W., Jr.
2011-04-01
Electromagnetic induction (EMI) instruments provide rapid, noninvasive, and spatially dense data for characterization of soil and groundwater properties. Data from multi-frequency EMI tools can be inverted to provide quantitative electrical conductivity estimates as a function of depth. In this study, multi-frequency EMI data collected across an abandoned uranium mill site near Naturita, Colorado, USA, are inverted to produce vertical distribution of electrical conductivity ( EC) across the site. The relation between measured apparent electrical conductivity ( ECa) and hydraulic conductivity ( K) is weak (correlation coefficient of 0.20), whereas the correlation between the depth dependent EC obtained from the inversions, and K is sufficiently strong to be used for hydrologic estimation (correlation coefficient of - 0.62). Depth-specific EC values were correlated with co-located K measurements to develop a site-specific ln( EC)-ln( K) relation. This petrophysical relation was applied to produce a spatially detailed map of K across the study area. A synthetic example based on ECa values at the site was used to assess model resolution and correlation loss given variations in depth and/or measurement error. Results from synthetic modeling indicate that optimum correlation with K occurs at ~ 0.5 m followed by a gradual correlation loss of 90% at 2.3 m. These results are consistent with an analysis of depth of investigation (DOI) given the range of frequencies, transmitter-receiver separation, and measurement errors for the field data. DOIs were estimated at 2.0 ± 0.5 m depending on the soil conductivities. A 4-layer model, with varying thicknesses, was used to invert the ECa to maximize available information within the aquifer region for improved correlations with K. Results show improved correlation between K and the corresponding inverted EC at similar depths, underscoring the importance of inversion in using multi-frequency EMI data for hydrologic estimation.
NASA Astrophysics Data System (ADS)
Hu, Y.; Ji, Y.; Egbert, G. D.
2015-12-01
The fictitious time domain method (FTD), based on the correspondence principle for wave and diffusion fields, has been developed and used over the past few years primarily for marine electromagnetic (EM) modeling. Here we present results of our efforts to apply the FTD approach to land and airborne TEM problems which can reduce the computer time several orders of magnitude and preserve high accuracy. In contrast to the marine case, where sources are in the conductive sea water, we must model the EM fields in the air; to allow for topography air layers must be explicitly included in the computational domain. Furthermore, because sources for most TEM applications generally must be modeled as finite loops, it is useful to solve directly for the impulse response appropriate to the problem geometry, instead of the point-source Green functions typically used for marine problems. Our approach can be summarized as follows: (1) The EM diffusion equation is transformed to a fictitious wave equation. (2) The FTD wave equation is solved with an explicit finite difference time-stepping scheme, with CPML (Convolutional PML) boundary conditions for the whole computational domain including the air and earth , with FTD domain source corresponding to the actual transmitter geometry. Resistivity of the air layers is kept as low as possible, to compromise between efficiency (longer fictitious time step) and accuracy. We have generally found a host/air resistivity contrast of 10-3 is sufficient. (3)A "Modified" Fourier Transform (MFT) allow us recover system's impulse response from the fictitious time domain to the diffusion (frequency) domain. (4) The result is multiplied by the Fourier transformation （FT） of the real source current avoiding time consuming convolutions in the time domain. (5) The inverse FT is employed to get the final full waveform and full time response of the system in the time domain. In general, this method can be used to efficiently solve most time-domain EM
Steady-state solutions for relativistically strong electromagnetic waves in plasmas.
NASA Technical Reports Server (NTRS)
Max, C. E.
1973-01-01
New steady-state solutions are derived which describe electromagnetic waves strong enough to make plasma ions and electrons relativistic. A two-fluid model is used throughout. The following solutions are studied: (1) linearly polarized waves with phase velocity much greater than c; (2) arbitrarily polarized waves with phase velocity near c, in a cold uniform plasma; (3) circularly polarized waves in a uniform plasma characterized by a scalar pressure tensor. All of these waves are capable of propagating in normally overdense plasmas, due to nonlinearities introduced by relativistic effects. The propagation of relativistically strong waves in a density gradient is examined, for the example of a circularly polarized wave strong enough to make electrons but not ions relativistic. It is shown that such a wave propagates at constant energy flux despite the nonlinearity of the system.
NASA Astrophysics Data System (ADS)
Yan, X.; Cai, D.; Nishikawa, K.; Lembege, B.
2004-12-01
We made our efforts to parallelize the global 3D HPF Electromagnetic particle model (EMPM) for several years and have also reported our meaningful simulation results that revealed the essential physics involved in interaction of the solar wind with the Earth's magnetosphere using this EMPM (Nishikawa et al., 1995; Nishikawa, 1997, 1998a, b, 2001, 2002) in our PC cluster and supercomputer(D.S. Cai et al., 2001, 2003). Sash patterns and related phenomena have been observed and reported in some satellite observations (Fujumoto et al. 1997; Maynard, 2001), and have motivated 3D MHD simulations (White and al., 1998). We also investigated it with our global 3D parallelized HPF EMPM with dawnward IMF By (K.-I. Nishikawa, 1998) and recently new simulation with dusk-ward IMF By was accomplished in the new VPP5000 supercomputer. In the new simulations performed on the new VPP5000 supercomputer of Tsukuba University, we used larger domain size, 305×205×205, smaller grid size (Δ ), 0.5R E(the radium of the Earth), more total particle number, 220,000,000 (about 8 pairs per cell). At first, we run this code until we get the so-called quasi-stationary status; After the quasi-stationary status was established, we applied a northward IMF (B z=0.2), and then wait until the IMF arrives around the magnetopuase. After the arrival of IMF, we begin to change the IMF from northward to duskward (IMF B y=-0.2). The results revealed that the groove structure at the day-side magnetopause, that causes particle entry into inner magnetosphere and the cross structure or S-structure at near magneto-tail are formed. Moreover, in contrast with MHD simulations, kinetic characteristic of this event is also analyzed self-consistently with this simulation. The new simulation provides new and more detailed insights for the observed sash event.
Electromagnetic dissociation of relativistic {sup 28}Si by nucleon emission
Sonnadara, U.J.
1992-12-01
A detailed study of the electromagnetic dissociation of {sup 28}Si by nucleon emission at E{sub lab}/A = 14.6 (GeV/nucleon was carried out with {sup 28}Si beams interacting on {sup 208}Pb). {sup 120}Sn. {sup 64}C targets. The measurements apparatus consists of detectors in the target area which measure the energy and charged multiplicity, and a forward spectrometer which measures the position, momentum and energy of the reaction fragments. The exclusive electromagnetic dissociation cross sections for decay channels having multiple nucleons in the final state have been measured which enables the selection of events produced in pure electromagnetic interactions. The measured cross sections agree well with previous measurements obtained for the removal of a few nucleons as well as with measurements on total charge removal cross sections from other experiments. The dependence of the integrated cross sections on the target charge Z{sub T} and the target mass AT confirms that for higher Z targets the excitation is largely electromagnetic. Direct measurements of the excitation energy for the electromagnetic dissociation of {sup 28}Si {yields} p+{sup 27}Al and {sup 28}Si {yields} n+{sup 27}Si have been obtained through a calculation of the invariant mass in kinematically, reconstructed events. The excitation energy spectrum for all targets peak near the isovector giant dipole resonance in {sup 28}Si. These distributions are well reproduced by combining the photon spectrum calculated using the Weizsaecker-Williams approximation with the experimental data on the photonuclear {sup 28}Si({sub {gamma},p}){sup 27}Al and {sup 28}Si({sub {gamma},n}){sup 27}Si. The possibilities of observing double giant dipole resonance excitations in {sup 28}Si have been investigated with cross section measurements as well as with excitation energy reconstruction.
Electromagnetic dissociation of relativistic [sup 28]Si by nucleon emission
Sonnadara, U.J.
1992-12-01
A detailed study of the electromagnetic dissociation of [sup 28]Si by nucleon emission at E[sub lab]/A = 14.6 (GeV/nucleon was carried out with [sup 28]Si beams interacting on [sup 208]Pb). [sup 120]Sn. [sup 64]C targets. The measurements apparatus consists of detectors in the target area which measure the energy and charged multiplicity, and a forward spectrometer which measures the position, momentum and energy of the reaction fragments. The exclusive electromagnetic dissociation cross sections for decay channels having multiple nucleons in the final state have been measured which enables the selection of events produced in pure electromagnetic interactions. The measured cross sections agree well with previous measurements obtained for the removal of a few nucleons as well as with measurements on total charge removal cross sections from other experiments. The dependence of the integrated cross sections on the target charge Z[sub T] and the target mass AT confirms that for higher Z targets the excitation is largely electromagnetic. Direct measurements of the excitation energy for the electromagnetic dissociation of [sup 28]Si [yields] p+[sup 27]Al and [sup 28]Si [yields] n+[sup 27]Si have been obtained through a calculation of the invariant mass in kinematically, reconstructed events. The excitation energy spectrum for all targets peak near the isovector giant dipole resonance in [sup 28]Si. These distributions are well reproduced by combining the photon spectrum calculated using the Weizsaecker-Williams approximation with the experimental data on the photonuclear [sup 28]Si([sub [gamma],p])[sup 27]Al and [sup 28]Si([sub [gamma],n])[sup 27]Si. The possibilities of observing double giant dipole resonance excitations in [sup 28]Si have been investigated with cross section measurements as well as with excitation energy reconstruction.
NASA Astrophysics Data System (ADS)
Scheunert, M.; Ullmann, A.; Afanasjew, M.; Börner, R.-U.; Siemon, B.; Spitzer, K.
2016-06-01
We present an inversion concept for helicopter-borne frequency-domain electromagnetic (HEM) data capable of reconstructing 3-D conductivity structures in the subsurface. Standard interpretation procedures often involve laterally constrained stitched 1-D inversion techniques to create pseudo-3-D models that are largely representative for smoothly varying conductivity distributions in the subsurface. Pronounced lateral conductivity changes may, however, produce significant artifacts that can lead to serious misinterpretation. Still, 3-D inversions of entire survey data sets are numerically very expensive. Our approach is therefore based on a cut-&-paste strategy whereupon the full 3-D inversion needs to be applied only to those parts of the survey where the 1-D inversion actually fails. The introduced 3-D Gauss-Newton inversion scheme exploits information given by a state-of-the-art (laterally constrained) 1-D inversion. For a typical HEM measurement, an explicit representation of the Jacobian matrix is inevitable which is caused by the unique transmitter-receiver relation. We introduce tensor quantities which facilitate the matrix assembly of the forward operator as well as the efficient calculation of the Jacobian. The finite difference forward operator incorporates the displacement currents because they may seriously affect the electromagnetic response at frequencies above 100. Finally, we deliver the proof of concept for the inversion using a synthetic data set with a noise level of up to 5%.
NASA Astrophysics Data System (ADS)
Deca, J.; Lapenta, G.; Divin, A. V.; Lembege, B.; Markidis, S.
2013-12-01
Unlike the Earth and Mercury, our Moon has no global magnetic field and is therefore not shielded from the impinging solar wind by a magnetosphere. However, lunar magnetic field measurements made by the Apollo missions provided direct evidence that the Moon has regions of small-scale crustal magnetic fields, ranging up to a few 100km in scale size with surface magnetic field strengths up to hundreds of nanoTeslas. More recently, the Lunar Prospector spacecraft has provided high-resolution observations allowing to construct magnetic field maps of the entire Moon, confirming the earlier results from Apollo, but also showing that the lunar plasma environment is much richer than earlier believed. Typically the small-scale magnetic fields are non-dipolar and rather tiny compared to the lunar radius and mainly clustered on the far side of the moon. Using iPic3D we present the first 3D fully kinetic and electromagnetic Particle-in-Cell simulations of the solar wind interaction with lunar magnetic anomalies. We study the behaviour of a dipole model with variable surface magnetic field strength under changing solar wind conditions and confirm that lunar crustal magnetic fields may indeed be strong enough to stand off the solar wind and form a mini-magnetosphere, as suggested by MHD and hybrid simulations and spacecraft observations. 3D-PIC simulations reveal to be very helpful to analyze the diversion/braking of the particle flux and the characteristics of the resulting particles accumulation. The particle flux to the surface is significantly reduced at the magnetic anomaly, surrounded by a region of enhanced density due to the magnetic mirror effect. Second, the ability of iPic3D to resolve all plasma components (heavy ions, protons and electrons) allows to discuss in detail the electron physics leading to the highly non-adiabatic interactions expected as well as the implications for solar wind shielding of the lunar surface, depending on the scale size (solar wind protons
THz electromagnetic radiation driven by intense relativistic electron beam based on ion focus regime
NASA Astrophysics Data System (ADS)
Zhou, Qing; Yang, Shengpeng; Xu, Jin; Zhang, Wenchao; Tang, Changjian; Duan, Zhaoyun; Gong, Yubin
2016-06-01
The simulation study finds that the relativistic electron beam propagating through the plasma background can produce electromagnetic (EM) radiation. With the propagation of the electron beam, the oscillations of the beam electrons in transverse and longitudinal directions have been observed simultaneously, which provides the basis for the electromagnetic radiation. The simulation results clearly show that the electromagnetic radiation frequency can reach up to terahertz (THz) wave band which may result from the filter-like property of plasma background, and the electromagnetic radiation frequency closely depends on the plasma density. To understand the above simulation results physically, the dispersion relation of the beam-plasma system has been derived using the field-matching method, and the dispersion curves show that the slow wave modes can couple with the electron beam effectively in THz wave band, which is an important theoretical evidence of the EM radiation.
NASA Astrophysics Data System (ADS)
Koldan, Jelena; Puzyrev, Vladimir; de la Puente, Josep; Houzeaux, Guillaume; Cela, José María
2014-06-01
We present an elaborate preconditioning scheme for Krylov subspace methods which has been developed to improve the performance and reduce the execution time of parallel node-based finite-element (FE) solvers for 3-D electromagnetic (EM) numerical modelling in exploration geophysics. This new preconditioner is based on algebraic multigrid (AMG) that uses different basic relaxation methods, such as Jacobi, symmetric successive over-relaxation (SSOR) and Gauss-Seidel, as smoothers and the wave front algorithm to create groups, which are used for a coarse-level generation. We have implemented and tested this new preconditioner within our parallel nodal FE solver for 3-D forward problems in EM induction geophysics. We have performed series of experiments for several models with different conductivity structures and characteristics to test the performance of our AMG preconditioning technique when combined with biconjugate gradient stabilized method. The results have shown that, the more challenging the problem is in terms of conductivity contrasts, ratio between the sizes of grid elements and/or frequency, the more benefit is obtained by using this preconditioner. Compared to other preconditioning schemes, such as diagonal, SSOR and truncated approximate inverse, the AMG preconditioner greatly improves the convergence of the iterative solver for all tested models. Also, when it comes to cases in which other preconditioners succeed to converge to a desired precision, AMG is able to considerably reduce the total execution time of the forward-problem code-up to an order of magnitude. Furthermore, the tests have confirmed that our AMG scheme ensures grid-independent rate of convergence, as well as improvement in convergence regardless of how big local mesh refinements are. In addition, AMG is designed to be a black-box preconditioner, which makes it easy to use and combine with different iterative methods. Finally, it has proved to be very practical and efficient in the
Interaction of relativistically strong electromagnetic waves with a layer of overdense plasma
Korzhimanov, A. V.; Eremin, V. I. Kim, A. V.; Tushentsov, M. R.
2007-10-15
Plasma-field structures that arise under the interaction between a relativistically strong electromagnetic wave and a layer of overdense plasma are considered within a quasistationary approximation. It is shown that, together with known solutions, which are nonlinear generalizations of skin-layer solutions, multilayer structures containing cavitation regions with completely removed electrons (ion layers) can be excited when the amplitude of the incident field exceeds a certain threshold value. Under symmetric irradiation, these cavitation regions, which play the role of self-consistent resonators, may amplify the field and accumulate electromagnetic energy.
DiPerna-Lions Flow for Relativistic Particles in an Electromagnetic Field
NASA Astrophysics Data System (ADS)
Jabin, P.-E.; Masmoudi, N.
2015-09-01
We show the existence and uniqueness of a DiPerna-Lions flow for relativistic particles subject to a Lorentz force in an electromagnetic field. The electric and magnetic fields solve the linear Maxwell system in the vacuum but for singular initial conditions which are only in the physical energy space. As the corresponding force field is only in L 2, we have to perform a careful analysis of the cancellations over a trajectory.
NASA Astrophysics Data System (ADS)
Petrov, P.; Newman, G. A.
2010-12-01
-Fourier domain we had developed 3D code for full-wave field simulation in the elastic media which take into account nonlinearity introduced by free-surface effects. Our approach is based on the velocity-stress formulation. In the contrast to conventional formulation we defined the material properties such as density and Lame constants not at nodal points but within cells. This second order finite differences method formulated in the cell-based grid, generate numerical solutions compatible with analytical ones within the range errors determinate by dispersion analysis. Our simulator will be embedded in an inversion scheme for joint seismic- electromagnetic imaging. It also offers possibilities for preconditioning the seismic wave propagation problems in the frequency domain. References. Shin, C. & Cha, Y. (2009), Waveform inversion in the Laplace-Fourier domain, Geophys. J. Int. 177(3), 1067- 1079. Shin, C. & Cha, Y. H. (2008), Waveform inversion in the Laplace domain, Geophys. J. Int. 173(3), 922-931. Commer, M. & Newman, G. (2008), New advances in three-dimensional controlled-source electromagnetic inversion, Geophys. J. Int. 172(2), 513-535. Newman, G. A., Commer, M. & Carazzone, J. J. (2010), Imaging CSEM data in the presence of electrical anisotropy, Geophysics, in press.
Samluk, Jesse P.; Geiger, Cathleen A.; Weiss, Chester J.; Kolodzey, James
2015-10-01
In this article we explore simulated responses of electromagnetic (EM) signals relative to in situ field surveys and quantify the effects that different values of conductivity in sea ice have on the EM fields. We compute EM responses of ice types with a three-dimensional (3-D) finite-volume discretization of Maxwell's equations and present 2-D sliced visualizations of their associated EM fields at discrete frequencies. Several interesting observations result: First, since the simulator computes the fields everywhere, each gridcell acts as a receiver within the model volume, and captures the complete, coupled interactions between air, snow, sea ice and sea water as a function of their conductivity; second, visualizations demonstrate how 1-D approximations near deformed ice features are violated. But the most important new finding is that changes in conductivity affect EM field response by modifying the magnitude and spatial patterns (i.e. footprint size and shape) of current density and magnetic fields. These effects are demonstrated through a visual feature we define as 'null lines'. Null line shape is affected by changes in conductivity near material boundaries as well as transmitter location. Our results encourage the use of null lines as a planning tool for better ground-truth field measurements near deformed ice types.
Samluk, Jesse P.; Geiger, Cathleen A.; Weiss, Chester J.; Kolodzey, James
2015-10-01
In this article we explore simulated responses of electromagnetic (EM) signals relative to in situ field surveys and quantify the effects that different values of conductivity in sea ice have on the EM fields. We compute EM responses of ice types with a three-dimensional (3-D) finite-volume discretization of Maxwell's equations and present 2-D sliced visualizations of their associated EM fields at discrete frequencies. Several interesting observations result: First, since the simulator computes the fields everywhere, each gridcell acts as a receiver within the model volume, and captures the complete, coupled interactions between air, snow, sea ice and sea water asmore » a function of their conductivity; second, visualizations demonstrate how 1-D approximations near deformed ice features are violated. But the most important new finding is that changes in conductivity affect EM field response by modifying the magnitude and spatial patterns (i.e. footprint size and shape) of current density and magnetic fields. These effects are demonstrated through a visual feature we define as 'null lines'. Null line shape is affected by changes in conductivity near material boundaries as well as transmitter location. Our results encourage the use of null lines as a planning tool for better ground-truth field measurements near deformed ice types.« less
Electromagnetic field evolution in relativistic heavy-ion collisions
Voronyuk, V.; Toneev, V. D.; Cassing, W.; Bratkovskaya, E. L.; Konchakovski, V. P.; Voloshin, S. A.
2011-05-15
The hadron string dynamics (HSD) model is generalized to include the creation and evolution of retarded electromagnetic fields as well as the influence of the magnetic and electric fields on the quasiparticle propagation. The time-space structure of the fields is analyzed in detail for noncentral Au + Au collisions at {radical}(s{sub NN})=200 GeV. It is shown that the created magnetic field is highly inhomogeneous, but in the central region of the overlapping nuclei it changes relatively weakly in the transverse direction. For the impact parameter b=10 fm, the maximal magnetic field - perpendicularly to the reaction plane - is obtained of order eB{sub y}/m{sub {pi}}{sup 2}{approx}5 for a very short time {approx}0.2 fm/c, which roughly corresponds to the time of a maximal overlap of the colliding nuclei. We find that at any time, the location of the maximum in the eB{sub y} distribution correlates with that of the energy density of the created particles. In contrast, the electric field distribution, being also highly inhomogeneous, has a minimum in the center of the overlap region. Furthermore, the field characteristics are presented as a function of the collision energy and the centrality of the collisions. To explore the effect of the back reaction of the fields on hadronic observables, a comparison of HSD results with and without fields is exemplified. Our actual calculations show no noticeable influence of the electromagnetic fields--created in heavy-ion collisions--on the effect of the electric charge separation with respect to the reaction plane.
NASA Astrophysics Data System (ADS)
Vujanovic, Gojko; Paquet, Jean-François; Denicol, Gabriel S.; Luzum, Matthew; Jeon, Sangyong; Gale, Charles
2016-07-01
The penetrating nature of electromagnetic signals makes them suitable probes to explore the properties of the strongly interacting medium created in relativistic nuclear collisions. We examine the effects of the initial conditions and shear relaxation time on the spectra and flow coefficients of electromagnetic probes, using an event-by-event 3+1-dimensional viscous hydrodynamic simulation (music).
A statistical model for relativistic quantum fluids interacting with an intense electromagnetic wave
NASA Astrophysics Data System (ADS)
Mahajan, Swadesh M.; Asenjo, Felipe A.
2016-05-01
A statistical model for relativistic quantum fluids interacting with an arbitrary amplitude circularly polarized electromagnetic wave is developed in two steps. First, the energy spectrum and the wave function for a quantum particle (Klein Gordon and Dirac) embedded in the electromagnetic wave are calculated by solving the appropriate eigenvalue problem. The energy spectrum is anisotropic in the momentum K and reflects the electromagnetic field through the renormalization of the rest mass m to M =√{m2+q2A2 } . Based on this energy spectrum of this quantum particle plus field combination (QPF), a statistical mechanics model of the quantum fluid made up of these weakly interacting QPF is developed. Preliminary investigations of the formalism yield highly interesting results—a new scale for temperature, and fundamental modification of the dispersion relation of the electromagnetic wave. It is expected that this formulation could, inter alia, uniquely advance our understanding of laboratory as well as astrophysical systems where one encounters arbitrarily large electromagnetic fields.
Standing electromagnetic solitons in hot ultra-relativistic electron-positron plasmas
Heidari, E.; Aslaninejad, M.; Eshraghi, H.; Rajaee, L.
2014-03-15
Using a one-dimensional self-consistent fluid model, we investigate standing relativistic bright solitons in hot electron-positron plasmas. The positron dynamics is taken into account. A set of nonlinear coupled differential equations describing the evolution of electromagnetic waves in fully relativistic two-fluid plasma is derived analytically and solved numerically. As a necessary condition for the existence of standing solitons the system should be relativistic. For the case of ultra-relativistic plasma, we investigate non-drifting bright solitary waves. Detailed discussions of the acceptable solutions are presented. New single hump non-trivial symmetric solutions for the scalar potential were found, and single and multi-nodal symmetric and anti-symmetric solutions for the vector potential are presented. It is shown that for a fixed value of the fluid velocity excited modes with more zeros in the profile of the vector potential show a higher magnitude for the scalar potential. An increase in the plasma fluid velocity also increases the magnitude of the scalar potential. Furthermore, the Hamiltonian and the first integral of the system are given.
Relativistic Two-Boson System in Presence of Electromagnetic Plane Wave
NASA Astrophysics Data System (ADS)
Droz-Vincent, Ph.
2016-09-01
The relativistic two-body problem is considered for spinless particles subject to an external electromagnetic field. When this field is made of the monochromatic superposition of two counter-propagating plane waves (and provided the mutual interaction between particles is known), it is possible to write down explicitly a pair of coupled wave equations (corresponding to a pair of mass-shell constraints) which takes into account also the field contribution. These equations are manifestly covariant; constants of the motion are exhibited, so one ends up with a reduced problem involving five degrees of freedom.
NASA Astrophysics Data System (ADS)
Paschalidis, Vasileios; Etienne, Zachariah B.; Shapiro, Stuart L.
2013-07-01
We perform the first general relativistic force-free simulations of neutron star magnetospheres in orbit about spinning and nonspinning black holes. We find promising precursor electromagnetic emission: typical Poynting luminosities at, e.g., an orbital separation of r=6.6RNS are LEM˜6×1042(BNS,p/1013G)2(MNS/1.4M⊙)2erg/s. The Poynting flux peaks within a broad beam of ˜40° in the azimuthal direction and within ˜60° from the orbital plane, establishing a possible lighthouse effect. Our calculations, though preliminary, preview more detailed simulations of these systems that we plan to perform in the future.
NASA Astrophysics Data System (ADS)
Artemyev, A. V.; Mourenas, D.; Agapitov, O. V.; Vainchtein, D. L.; Mozer, F. S.; Krasnoselskikh, V.
2015-08-01
In the present paper, we investigate the trapping of relativistic electrons by intense whistler-mode waves or electromagnetic ion cyclotron waves in the Earth's radiation belts. We consider the non-resonant impact of additional, lower amplitude magnetic field fluctuations on the stability of electron trapping. We show that such additional non-resonant fluctuations can break the adiabatic invariant corresponding to trapped electron oscillations in the effective wave potential. This destruction results in a diffusive escape of electrons from the trapped regime of motion and thus can lead to a significant reduction of the efficiency of electron acceleration. We demonstrate that when energetic electrons are trapped by intense parallel or very oblique whistler-mode waves, non-resonant magnetic field fluctuations in the whistler-mode frequency range with moderate amplitudes around 3 -15 pT (much less intense than the primary waves) can totally disrupt the trapped motion. However, the trapping of relativistic electrons by electromagnetic ion cyclotron waves is noticeably more stable. We also discuss how the proposed approach can be used to estimate the effects of wave amplitude modulations on the motion of trapped particles.
NASA Astrophysics Data System (ADS)
Deffenbaugh, Paul Issac
3D printing has garnered immense attention from many fields including in-office rapid prototyping of mechanical parts, outer-space satellite replication, garage functional firearm manufacture, and NASA rocket engine component fabrication. 3D printing allows increased design flexibility in the fabrication of electronics, microwave circuits and wireless antennas and has reached a level of maturity which allows functional parts to be printed. Much more work is necessary in order to perfect the processes of 3D printed electronics especially in the area of automation. Chapter 1 shows several finished prototypes of 3D printed electronics as well as newly developed techniques in fabrication. Little is known about the RF and microwave properties and applications of the standard materials which have been developed for 3D printing. Measurement of a wide variety of materials over a broad spectrum of frequencies up to 10 GHz using a variety of well-established measurement methods is performed throughout chapter 2. Several types of high frequency RF transmission lines are fabricated and valuable model-matched data is gathered and provided in chapter 3 for future designers' use. Of particular note is a fully 3D printed stripline which was automatically fabricated in one process on one machine. Some core advantages of 3D printing RF/microwave components include rapid manufacturing of complex, dimensionally sensitive circuits (such as antennas and filters which are often iteratively tuned) and the ability to create new devices that cannot be made using standard fabrication techniques. Chapter 4 describes an exemplary fully 3D printed curved inverted-F antenna.
NASA Technical Reports Server (NTRS)
Sakai, J. I.; Zhao, J.; Nishikawa, K.-I.
1994-01-01
We have shown that a current-carrying plasma loop can be heated by magnetic pinch driven by the pressure imbalance between inside and outside the loop, using a 3-dimensional electromagnetic (EM) particle code. Both electrons and ions in the loop can be heated in the direction perpendicular to the ambient magnetic field, therefore the perpendicular temperature can be increased about 10 times compared with the parallel temperature. This temperature anisotropy produced by the magnetic pinch heating can induce a plasma instability, by which high-frequency electromagnetic waves can be excited. The plasma current which is enhanced by the magnetic pinch can also excite a kinetic kink instability, which can heat ions perpendicular to the magnetic field. The heating mechanism of ions as well as the electromagnetic emission could be important for an understanding of the coronal loop heating and the electromagnetic wave emissions from active coronal regions.
NASA Astrophysics Data System (ADS)
Kholmetskii, Alexander; Missevitch, Oleg; Yarman, Tolga
2016-02-01
We address to the Poynting theorem for the bound (velocity-dependent) electromagnetic field, and demonstrate that the standard expressions for the electromagnetic energy flux and related field momentum, in general, come into the contradiction with the relativistic transformation of four-vector of total energy-momentum. We show that this inconsistency stems from the incorrect application of Poynting theorem to a system of discrete point-like charges, when the terms of self-interaction in the product {\\varvec{j}} \\cdot {\\varvec{E}} (where the current density {\\varvec{j}} and bound electric field {\\varvec{E}} are generated by the same source charge) are exogenously omitted. Implementing a transformation of the Poynting theorem to the form, where the terms of self-interaction are eliminated via Maxwell equations and vector calculus in a mathematically rigorous way (Kholmetskii et al., Phys Scr 83:055406, 2011), we obtained a novel expression for field momentum, which is fully compatible with the Lorentz transformation for total energy-momentum. The results obtained are discussed along with the novel expression for the electromagnetic energy-momentum tensor.
Electromagnetic heavy-lepton pair production in relativistic heavy-ion collisions
NASA Astrophysics Data System (ADS)
Şengül, M. Y.; Güçlü, M. C.; Mercan, Ö.; Karakuş, N. G.
2016-08-01
We calculate the cross sections of electromagnetic productions of muon- and tauon-pair productions from the ultra-relativistic heavy ion collisions. Since the Compton wavelengths of muon and tauon are comparable to the radius of the colliding ions, nuclear form factors play important roles for calculating the cross sections. Recent measurement (Abrahamyan et al., Phys Rev Lett 108:112502, 2012) indicates that the neutrons are differently distributed from the protons; therefore this affects the cross section of the heavy-lepton pair production. In order to see the effects of the neutron distributions in the nucleus, we used analytical expression of the Fourier transforms of the Wood-Saxon distribution. Cross section calculations show that the Wood-Saxon distribution function is more sensitive to the parameter R compared to the parameter a.
NASA Astrophysics Data System (ADS)
Ryutov, D. D.
2008-11-01
The physics of relativistic electron particle beams propagating through the plasma is of a significant interest for laboratory astrophysics, fast ignition, and Z-pinch research. Most attention has been directed towards the analysis of electromagnetic filamentation instabilities. On the other hand, there exists a broad class of very powerful electrostatic instabilities, e.g., the Buneman instability. The author considers in a unified fashion linear theory for both types of instabilities under conditions where there is no magnetic field in an unperturbed state (i.e., the beam current is fully neutralized by the plasma current). The following factors are taken into account: the beam energy and angular spread; plasma non-uniformity; particle collisions in the background plasma. The areas of the parameter domain where particular instability is prevalent are identified; the results are presented in the form of several easy-to-use diagrams. Prepared by LLNL under contract DE-AC52-07NA27344.
NASA Astrophysics Data System (ADS)
Gill, Tarsem Singh; Kaur, Ravinder; Mahajan, Ranju
2010-09-01
This paper presents an analysis of self-consistent, steady-state, theoretical model, which explains the ring formation in a Gaussian electromagnetic beam propagating in a magnetoplasma, characterized by relativistic nonlinearity. Higher order terms (up to r4) in the expansion of the dielectric function and the eikonal have been taken into account. The condition for the formation of a dark and bright ring derived earlier by Misra and Mishra [J. Plasma Phys. 75, 769 (2009)] has been used to study focusing/defocusing of the beam. It is seen that inclusion of higher order terms does significantly affect the dependence of the beam width on the distance of propagation. Further, the effect of the magnetic field and the nature of nonlinearity on the ring formation and self-focusing of the beam have been explored.
NASA Astrophysics Data System (ADS)
Cuperman, S.; Heristchi, D.
1992-08-01
The transcendental dispersion equation for electromagnetic waves propagating in the slow mode in sheared non-neutral relativistic cylindrical electron beams in strong applied magnetic fields is solved exactly. Thus, rather than truncated power series for the modified Bessel functions involved, use is made of modern algorithms able to compute such functions up to 18-digit accuracy. Consequently, new and significantly more important branches of the velocity shear instability are found. When the shear-factor and/or the geometrical parameter a/b (pipe-to-beam radius ratio) are increased, the unstable branches join, and the higher-frequency, larger-wavenumber modes are significantly enhanced. Since analytical solutions of the exact dispersion relation cannot be obtained, it is suggested that in all similar cases the methods proposed and demonstrated here should be used to carry out a rigorous stability analysis.
Hodgdon, M.L.; Oona, H.; Martinez, A.R.; Salon, S.; Wendling, P.; Krahenbuhl, L.; Nicolas, A.; Nicolas, L.
1989-01-01
We present herein the results of three electromagnetic field problems for compressed magnetic field generators and their associated power flow channels. The first problem is the computation of the transient magnetic field in a two-dimensional model of helical generator during loading. The second problem is the three-dimensional eddy current patterns in a section of an armature beneath a bifurcation point of a helical winding. Our third problem is the calculation of the three-dimensional electrostatic fields in a region known as the post-hole convolute in which a rod connects the inner and outer walls of a system of three concentric cylinders through a hole in the middle cylinder. While analytic solutions exist for many electromagnetic field problems in cases of special and ideal geometries, the solutions of these and similar problems for the proper analysis and design of compressed magnetic field generators and their related hardware require computer simulations. In earlier studies, computer models have been proposed, several based on research oriented hydrocodes to which uncoupled or partially coupled Maxwell's equations solvers are added. Although the hydrocode models address the problem of moving, deformable conductors, they are not useful for electromagnetic analysis, nor can they be considered design tools. For our studies, we take advantage of the commercial, electromagnetic computer-aided design software packages FLUX2D nd PHI3D that were developed for motor manufacturers and utilities industries. 4 refs., 6 figs.
Cremaschini, Claudio Stuchlík, Zdeněk; Tessarotto, Massimo
2014-03-15
The kinetic description of relativistic plasmas in the presence of time-varying and spatially non-uniform electromagnetic (EM) fields is a fundamental theoretical issue both in astrophysics and plasma physics. This refers, in particular, to the treatment of collisionless and strongly-magnetized plasmas in the presence of intense radiation sources. In this paper, the problem is investigated in the framework of a covariant gyrokinetic treatment for Vlasov–Maxwell equilibria. The existence of a new class of kinetic equilibria is pointed out, which occur for spatially-symmetric systems. These equilibria are shown to exist in the presence of non-uniform background EM fields and curved space-time. In the non-relativistic limit, this feature permits the determination of kinetic equilibria even for plasmas in which particle energy is not conserved due to the occurrence of explicitly time-dependent EM fields. Finally, absolute stability criteria are established which apply in the case of infinitesimal symmetric perturbations that can be either externally or internally produced.
NASA Astrophysics Data System (ADS)
Manuwong, Theerapatt; Zhang, Wei; Kazinczi, Peter Lobo; Bodey, Andrew J.; Rau, Christoph; Mi, Jiawei
2015-07-01
A novel programmable electromagnetic pulse device was developed and used to study the solidification of Al-15 pct Cu and Al-35 pct Cu alloys. The pulsed magnetic fluxes and Lorentz forces generated inside the solidifying melts were simulated using finite element methods, and their effects on the solidification microstructures were characterized using electron microscopy and synchrotron X-ray tomography. Using a discharging voltage of 120 V, a pulsed magnetic field with the peak Lorentz force of ~1.6 N was generated inside the solidifying Al-Cu melts which were showed sufficiently enough to disrupt the growth of the primary Al dendrites and the Al2Cu intermetallic phases. The microstructures exhibit a strong correlation to the characteristics of the applied pulse, forming a periodical pattern that resonates the frequency of the applied electromagnetic field.
NASA Astrophysics Data System (ADS)
Gubchenko, V. M.
2015-12-01
In part I of the work, the physical effects responsible for the formation of low-speed flows in plasma coronas, coupled with formation of coronas magnetosphere-like structures, are described qualitatively. Coronal domain structures form if we neglect scales of spatial plasma dispersion: high-speed flows are accumulated in magnetic tubes of the open domains, while magnetic structures and low-speed flows are concentrated within boundaries of domains. The inductive electromagnetic process occurring in flows of the hot collisionless plasma is shown to underlie the formation of magnetosphere-like structures. Depending on the form of the velocity distribution function of particles (PDF), a hot flow differently reveals its electromagnetic properties, which are expressed by the induction of resistive and diamagnetic scales of spatial dispersion. These determine the magnetic structure scales and structure reconstruction. The inductive electromagnetic process located in lines of the plasma nontransparency and absorption, in which the structures of excited fields are spatially aperiodic and skinned to the magnetic field sources. The toroidal and dipole magnetic sources of different configurations are considered for describing the corona structures during the solar maximum and solar minimum.
NASA Astrophysics Data System (ADS)
Schultz, A.
2010-12-01
3D forward solvers lie at the core of inverse formulations used to image the variation of electrical conductivity within the Earth's interior. This property is associated with variations in temperature, composition, phase, presence of volatiles, and in specific settings, the presence of groundwater, geothermal resources, oil/gas or minerals. The high cost of 3D solutions has been a stumbling block to wider adoption of 3D methods. Parallel algorithms for modeling frequency domain 3D EM problems have not achieved wide scale adoption, with emphasis on fairly coarse grained parallelism using MPI and similar approaches. The communications bandwidth as well as the latency required to send and receive network communication packets is a limiting factor in implementing fine grained parallel strategies, inhibiting wide adoption of these algorithms. Leading Graphics Processor Unit (GPU) companies now produce GPUs with hundreds of GPU processor cores per die. The footprint, in silicon, of the GPU's restricted instruction set is much smaller than the general purpose instruction set required of a CPU. Consequently, the density of processor cores on a GPU can be much greater than on a CPU. GPUs also have local memory, registers and high speed communication with host CPUs, usually through PCIe type interconnects. The extremely low cost and high computational power of GPUs provides the EM geophysics community with an opportunity to achieve fine grained (i.e. massive) parallelization of codes on low cost hardware. The current generation of GPUs (e.g. NVidia Fermi) provides 3 billion transistors per chip die, with nearly 500 processor cores and up to 6 GB of fast (DDR5) GPU memory. This latest generation of GPU supports fast hardware double precision (64 bit) floating point operations of the type required for frequency domain EM forward solutions. Each Fermi GPU board can sustain nearly 1 TFLOP in double precision, and multiple boards can be installed in the host computer system. We
Campione, Salvatore; Capolino, Filippo
2016-01-25
In this study, we investigate the effect on wave propagation of array packing and electromagnetic coupling between spheres in a three-dimensional (3D) lattice of microspheres with large permittivity that exhibit strong magnetic polarizability. We report on the complex wavenumber of Bloch waves in the lattice when each sphere is assumed to possess both electric and magnetic dipoles and full electromagnetic coupling is accounted for. While for small material-filling fractions we always determine one dominant mode with low attenuation constant, the same does not happen for large filling fractions, when electromagnetic coupling is included. In the latter case we peculiarly observemore » two dominant modes with low attenuation constant, dominant in different frequency ranges. The filling fraction threshold for which two dominant modes appear varies for different metamaterial constituents, as proven by considering spheres made by either titanium dioxide or lead telluride. As further confirmation of our findings, we retrieve the complex propagation constant of the dominant mode(s) via a field fitting procedure employing two sets of waves (direct and reflected) pertaining to two distinct modes, strengthening the presence of the two distinct dominant modes for increasing filling fractions. However, given that one mode only, with transverse polarization, at any given frequency, is dominant and able to propagate inside the lattice, we are able to accurately treat the metamaterial that is known to exhibit artificial magnetism as a homogeneous material with effective parameters, such as the refractive index. Results clearly show that the account of both electric and magnetic scattering processes in evaluating all electromagnetic intersphere couplings is essential for a proper description of the electromagnetic propagation in lattices.« less
NASA Astrophysics Data System (ADS)
Zheleznyakov, V. V.; Bespalov, P. A.
2016-04-01
In part I of this work [1], we study the dispersion characteristics of low-frequency waves in a relativistic electron-positron plasma. In part II, we examine the electromagnetic wave instability in this plasma caused by an admixture of nonrelativistic protons with energy comparable with the energy of relativistic low-mass particles. The instability occurs in the frequency band between the fundamental harmonic of proton gyrofrequency and the fundamental harmonic of relativistic electron gyrofrequency. The results can be used for the interpretation of known observations of the pulsar emissions obtained with a high time and frequency resolution. The considered instability can probably be the initial stage of the microwave radio emission nanoshots typical of the pulsar in the Crab Nebula.
NASA Astrophysics Data System (ADS)
Zheleznyakov, V. V.; Bespalov, P. A.
2016-05-01
In part I of this work [1], we study the dispersion characteristics of low-frequency waves in a relativistic electron-positron plasma. In part II, we examine the electromagnetic wave instability in this plasma caused by an admixture of nonrelativistic protons with energy comparable with the energy of relativistic low-mass particles. The instability occurs in the frequency band between the fundamental harmonic of proton gyrofrequency and the fundamental harmonic of relativistic electron gyrofrequency. The results can be used for the interpretation of known observations of the pulsar emissions obtained with a high time and frequency resolution. The considered instability can probably be the initial stage of the microwave radio emission nanoshots typical of the pulsar in the Crab Nebula.
NASA Astrophysics Data System (ADS)
Yi, Myeong-Jong; Sasaki, Yutaka
2015-11-01
Frequency-domain loop-loop electromagnetic (EM) methods are sensitive to the magnetic susceptibility of the Earth as well as its resistivity. Thus, inversion techniques have been used to simultaneously reconstruct both resistivity and susceptibility models from EM data. However, to take full advantage of inversion methods, calibration errors must be assessed and removed because ignoring them can result in misleading models. We present a multidimensional inversion method that jointly inverts EM and direct current (DC) resistivity data to derive offset errors as well as resistivity and susceptibility models, assuming that calibration errors can be represented by in-phase and quadrature offsets at each frequency. Addition of independent data such as DC data is effective for more accurately estimating the offsets, resulting in more reliable subsurface models. Synthetic examples involving small-loop EM data show that simultaneous inversion for resistivity and susceptibility is not stable, because of strong correlations between in-phase offset parameters and background susceptibility, but that the offsets are well determined when the data misfit is reduced rapidly in the early iteration step. Improvements achieved by joint inversion are mainly on the resistivity model. For airborne electromagnetic (AEM) data, the inversion process is stable, because AEM data are acquired using more loop-loop geometries and a wider range of frequencies. As a result, both the resistivity and susceptibility models are significantly improved by joint inversion.
Xie, G.; Li, J.
1997-05-01
A new 3D electromagnetic modeling and nonlinear inversion algorithm is presented based on global integral and local differential equations decomposition (GILD). The GILD parallel nonlinear inversion algorithm consists of five parts: (1) the domain is decomposed into subdomain SI and subdomain SII; (2) a new global magnetic integral equation in SI and the local magnetic differential equations IN SII will be used together to obtain the magnetic field in the modeling step; (3) the new global magnetic integral Jacobian equation in SI and the local magnetic differential Jacobian equations in SII will be used together to update the electric conductivity and permittivity from the magnetic field data in the inversion step; (4) the subdomain SII can naturally and uniformly be decomposed into 2{sup n} smaller sub-cubic-domains; the sparse matrix in each sub-cubic-domain can be eliminated separately, in parallel; (5) a new parallel multiple hierarchy substructure algorithm will be used to solve the smaller full matrices in SI, in parallel. The applications of the new 3D parallel GILD EM modeling and nonlinear inversion algorithm and software are: (1) to create high resolution controlled-source electric conductivity and permittivity imaging for interpreting electromagnetic field data acquired from cross hole, surface to borehole, surface to surface, single hole, and multiple holes; (2) to create the magnetotelluric high resolution imaging from the surface impedance and field data. The new GILD parallel nonlinear inversion will be a 3D/2.5D powerful imaging tool for the oil geophysical exploration and environmental remediation and monitoring.
NASA Astrophysics Data System (ADS)
Shemelya, Corey M.; Rivera, Armando; Perez, Angel Torrado; Rocha, Carmen; Liang, Min; Yu, Xiaoju; Kief, Craig; Alexander, David; Stegeman, James; Xin, Hao; Wicker, Ryan B.; MacDonald, Eric; Roberson, David A.
2015-08-01
Material-extrusion three-dimensional (3D) printing has recently attracted much interest because of its process flexibility, rapid response to design alterations, and ability to create structures "on-the-go". For this reason, 3D printing has possible applications in rapid creation of space-based devices, for example cube satellites (CubeSat). This work focused on fabrication and characterization of tungsten-doped polycarbonate polymer matrix composites specifically designed for x-ray radiation-shielding applications. The polycarbonate-tungsten polymer composite obtained intentionally utilizes low loading levels to provide x-ray shielding while limiting effects on other properties of the material, for example weight, electromagnetic functionality, and mechanical strength. The fabrication process, from tungsten functionalization to filament extrusion and material characterization, is described, including printability, determination of x-ray attenuation, tensile strength, impact resistance, and gigahertz permittivity, and failure analysis. The proposed materials are uniquely advantageous when implemented in 3D printed structures, because even a small volume fraction of tungsten has been shown to substantially alter the properties of the resulting composite.
Commer, Michael; Helwig, Stefan, L.; Hordt, Andreas; Scholl,Carsten; Tezkan, Bulent
2006-06-14
Three long-offset transient electromagnetic (LOTEM) surveyswerecarried out at the active volcano Merapi in Central Java (Indonesia)during the years 1998, 2000, and 2001. The measurements focused on thegeneral resistivity structure of the volcanic edifice at depths of 0.5-2km and the further investigation of a southside anomaly. The measurementswere insufficient for a full 3D inversion scheme, which could enable theimaging of finely discretized resistivity distributions. Therefore, astable, damped least-squares joint-inversion approach is used to optimize3D models with a limited number of parameters. The mode ls feature therealistic simulation of topography, a layered background structure, andadditional coarse 3D blocks representing conductivity anomalies.Twenty-eight LOTEM transients, comprising both horizontal and verticalcomponents of the magnetic induction time derivative, were analyzed. Inview of the few unknowns, we were able to achieve reasonable data fits.The inversion results indicate an upwelling conductor below the summit,suggesting hydrothermal activity in the central volcanic complex. Ashallow conductor due to a magma-filled chamber, at depths down to 1 kmbelow the summit, suggested by earlier seismic studies, is not indicatedby the inversion results. In conjunction with an anomalous-density model,derived from arecent gravity study, our inversion results provideinformation about the southern geological structure resulting from amajor sector collapse during the Middle Merapi period. The density modelallows to assess a porosity range andthus an estimated vertical salinityprofile to explain the high conductivities on a larger scale, extendingbeyond the foothills of Merapi.
Milant'ev, V. P. Castillo, A. J.
2013-04-15
Averaged relativistic equations of motion of a charged particle in the field of intense electromagnetic radiation have been obtained in the geometrical optics approximation using the Bogoliubov method. Constraints are determined under which these equations are valid. Oscillating additions to the smoothed dynamical variables of the particle have been found; they are reduced to known expressions in the case of the circularly and linearly polarized plane waves. It has been shown that the expressions for the averaged relativistic force in both cases contain new additional small terms weakening its action. The known difference between the expressions for the ponderomotive force in the cases of circularly and linearly polarized waves has been confirmed.
Bauer, Carl A.; Werner, Gregory R.; Cary, John R.
2011-03-01
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.
Colavita, E.; Hacyan, S.
2014-03-15
We analyze the solutions of the Klein–Gordon and Dirac equations describing a charged particle in an electromagnetic plane wave combined with a magnetic field parallel to the direction of propagation of the wave. It is shown that the Klein–Gordon equation admits coherent states as solutions, while the corresponding solutions of the Dirac equation are superpositions of coherent and displaced-number states. Particular attention is paid to the resonant case in which the motion of the particle is unbounded. -- Highlights: •We study a relativistic electron in a particular electromagnetic field configuration. •New exact solutions of the Klein–Gordon and Dirac equations are obtained. •Coherent and displaced number states can describe a relativistic particle.
Boutaleb, Samir; Fillion, Olivier; Bonillas, Antonio; Hautvast, Gilion; Binnekamp, Dirk; Beaulieu, Luc
2015-01-01
Purpose Accurate insertion and overall needle positioning are key requirements for effective brachytherapy treatments. This work aims at demonstrating the accuracy performance and the suitability of the Aurora® V1 Planar Field Generator (PFG) electromagnetic tracking system (EMTS) for real-time treatment assistance in interstitial brachytherapy procedures. Material and methods The system's performance was characterized in two distinct studies. First, in an environment free of EM disturbance, the boundaries of the detection volume of the EMTS were characterized and a tracking error analysis was performed. Secondly, a distortion analysis was conducted as a means of assessing the tracking accuracy performance of the system in the presence of potential EM disturbance generated by the proximity of standard brachytherapy components. Results The tracking accuracy experiments showed that positional errors were typically 2 ± 1 mm in a zone restricted to the first 30 cm of the detection volume. However, at the edges of the detection volume, sensor position errors of up to 16 mm were recorded. On the other hand, orientation errors remained low at ± 2° for most of the measurements. The EM distortion analysis showed that the presence of typical brachytherapy components in vicinity of the EMTS had little influence on tracking accuracy. Position errors of less than 1 mm were recorded with all components except with a metallic arm support, which induced a mean absolute error of approximately 1.4 mm when located 10 cm away from the needle sensor. Conclusions The Aurora® V1 PFG EMTS possesses a great potential for real-time treatment assistance in general interstitial brachytherapy. In view of our experimental results, we however recommend that the needle axis remains as parallel as possible to the generator surface during treatment and that the tracking zone be restricted to the first 30 cm from the generator surface. PMID:26622231
Poulin, E; Racine, E; Beaulieu, L; Binnekamp, D
2014-06-15
Purpose: In high dose rate brachytherapy (HDR-B), actual catheter reconstruction protocols are slow and errors prompt. The purpose of this study was to evaluate the accuracy and robustness of an electromagnetic (EM) tracking system for improved catheter reconstruction in HDR-B protocols. Methods: For this proof-of-principle, a total of 10 catheters were inserted in gelatin phantoms with different trajectories. Catheters were reconstructed using a Philips-design 18G biopsy needle (used as an EM stylet) and the second generation Aurora Planar Field Generator from Northern Digital Inc. The Aurora EM system exploits alternating current technology and generates 3D points at 40 Hz. Phantoms were also scanned using a μCT (GE Healthcare) and Philips Big Bore clinical CT system with a resolution of 0.089 mm and 2 mm, respectively. Reconstructions using the EM stylet were compared to μCT and CT. To assess the robustness of the EM reconstruction, 5 catheters were reconstructed twice and compared. Results: Reconstruction time for one catheter was 10 seconds or less. This would imply that for a typical clinical implant of 17 catheters, the total reconstruction time would be less than 3 minutes. When compared to the μCT, the mean EM tip identification error was 0.69 ± 0.29 mm while the CT error was 1.08 ± 0.67 mm. The mean 3D distance error was found to be 0.92 ± 0.37 mm and 1.74 ± 1.39 mm for the EM and CT, respectively. EM 3D catheter trajectories were found to be significantly more accurate (unpaired t-test, p < 0.05). A mean difference of less than 0.5 mm was found between successive EM reconstructions. Conclusion: The EM reconstruction was found to be faster, more accurate and more robust than the conventional methods used for catheter reconstruction in HDR-B. This approach can be applied to any type of catheters and applicators. We would like to disclose that the equipments, used in this study, is coming from a collaboration with Philips Medical.
NASA Astrophysics Data System (ADS)
von Hebel, Christian; Rudolph, Sebastian; Huisman, Johan A.; van der Kruk, Jan; Vereecken, Harry
2013-04-01
three different coil offsets in HCP and VCP measurement modes. This resulted in six high spatial resolution data sets of approximately 60000 measurements with different sensing depths. A 5 m block-kriging was applied to all six data sets to re-grid the sampling points on the same regular grid. For each grid node, the six measured apparent conductivities were used in a three-layer inversion. The three-layer inversion results of electrical conductivity thus obtained were used to derive a three-dimensional (3D) model of subsurface heterogeneity, which clearly indicated lateral and vertical conductivity changes of the subsurface that are related to changes in soil texture and soil water content.
Hatt, Charles R.; Jain, Ameet K.; Parthasarathy, Vijay; Lang, Andrew; Raval, Amish N.
2014-01-01
Myocardial infarction (MI) is one of the leading causes of death in the world. Small animal studies have shown that stem-cell therapy offers dramatic functional improvement post-MI. An endomyocardial catheter injection approach to therapeutic agent delivery has been proposed to improve efficacy through increased cell retention. Accurate targeting is critical for reaching areas of greatest therapeutic potential while avoiding a life-threatening myocardial perforation. Multimodal image fusion has been proposed as a way to improve these procedures by augmenting traditional intra-operative imaging modalities with high resolution pre-procedural images. Previous approaches have suffered from a lack of real-time tissue imaging and dependence on X-ray imaging to track devices, leading to increased ionizing radiation dose. In this paper, we present a new image fusion system for catheter-based targeted delivery of therapeutic agents. The system registers real-time 3D echocardiography, magnetic resonance, X-ray, and electromagnetic sensor tracking within a single flexible framework. All system calibrations and registrations were validated and found to have target registration errors less than 5 mm in the worst case. Injection accuracy was validated in a motion enabled cardiac injection phantom, where targeting accuracy ranged from 0.57 to 3.81 mm. Clinical feasibility was demonstrated with in-vivo swine experiments, where injections were successfully made into targeted regions of the heart. PMID:23561056
NASA Astrophysics Data System (ADS)
Balsara, Dinshaw S.; Amano, Takanobu; Garain, Sudip; Kim, Jinho
2016-08-01
In various astrophysics settings it is common to have a two-fluid relativistic plasma that interacts with the electromagnetic field. While it is common to ignore the displacement current in the ideal, classical magnetohydrodynamic limit, when the flows become relativistic this approximation is less than absolutely well-justified. In such a situation, it is more natural to consider a positively charged fluid made up of positrons or protons interacting with a negatively charged fluid made up of electrons. The two fluids interact collectively with the full set of Maxwell's equations. As a result, a solution strategy for that coupled system of equations is sought and found here. Our strategy extends to higher orders, providing increasing accuracy. The primary variables in the Maxwell solver are taken to be the facially-collocated components of the electric and magnetic fields. Consistent with such a collocation, three important innovations are reported here. The first two pertain to the Maxwell solver. In our first innovation, the magnetic field within each zone is reconstructed in a divergence-free fashion while the electric field within each zone is reconstructed in a form that is consistent with Gauss' law. In our second innovation, a multidimensionally upwinded strategy is presented which ensures that the magnetic field can be updated via a discrete interpretation of Faraday's law and the electric field can be updated via a discrete interpretation of the generalized Ampere's law. This multidimensional upwinding is achieved via a multidimensional Riemann solver. The multidimensional Riemann solver automatically provides edge-centered electric field components for the Stokes law-based update of the magnetic field. It also provides edge-centered magnetic field components for the Stokes law-based update of the electric field. The update strategy ensures that the electric field is always consistent with Gauss' law and the magnetic field is always divergence-free. This
Andreev, A A; Galkin, A L; Kalashnikov, M P; Korobkin, V V; Romanovsky, Mikhail Yu; Shiryaev, O B
2011-08-31
We study the motion of an electron and emission of electromagnetic waves by an electron in the field of a relativistically intense laser pulse. The dynamics of the electron is described by the Newton equation with the Lorentz force in the right-hand side. It is shown that the electrons may be ejected from the interaction region with high energy. The energy spectrum of these electrons and the technique of using the spectrum to assess the maximal intensity in the focus are analysed. It is found that electromagnetic radiation of an electron moving in an intense laser field occurs within a small angle around the direction of the electron trajectory tangent. The tangent quickly changes its direction in space; therefore, electromagnetic radiation of the electron in the far-field zone in a certain direction in the vicinity of the tangent is a short pulse with a duration as short as zeptoseconds. The calculation of the temporary and spectral distribution of the radiation field is carried out. (superintense laser fields)
Relativistic electron precipitation events driven by electromagnetic ion-cyclotron waves
Khazanov, G. Sibeck, D.; Tel'nikhin, A.; Kronberg, T.
2014-08-15
We adopt a canonical approach to describe the stochastic motion of relativistic belt electrons and their scattering into the loss cone by nonlinear EMIC waves. The estimated rate of scattering is sufficient to account for the rate and intensity of bursty electron precipitation. This interaction is shown to result in particle scattering into the loss cone, forming ∼10 s microbursts of precipitating electrons. These dynamics can account for the statistical correlations between processes of energization, pitch angle scattering, and relativistic electron precipitation events, that are manifested on large temporal scales of the order of the diffusion time ∼tens of minutes.
High frequency electromagnetic modes in a weakly magnetized relativistic electron plasma
Abbas, Gohar; Murtaza, G.; Kingham, R. J.
2010-07-15
Using the linearized Vlasov-Maxwell model, the polarization tensor for a weakly magnetized electron plasma is derived. For isotropic relativistic Maxwellian velocity distribution function, dispersion relations are obtained for both parallel and perpendicular propagations. The integrals (called Meijer G functions) that arise due to relativistic effects are examined in various limits and dispersion relations are derived for the nonrelativistic, weakly, strongly, and ultrarelativistic Maxwellian velocity distributions. It is generally observed that the propagation domains of the modes are enlarged as one proceeds from the nonrelativistic to the highly relativistic regime. Resultantly, due to the relativistic effects, the Whistler mode is suppressed in the R-wave, the nonpropagation band of X-mode is reduced, and the X-mode itself approaches the O-mode. Further, the results derived in the ultra- and nonrelativistic limits found to be in agreement with the earlier calculations [G. Abbas et al. Phys. Scr. 76, 649 (2007); F. F. Chen, Introduction to Plasma Physics and Controlled Fusion (Plenum, New York, 1984), Vol. 1].
NASA Astrophysics Data System (ADS)
Shebeko, A. V.; Shirokov, M. I.
2000-04-01
A growing interest in the method of unitary transformations (UT's) in the quantum theory of particles and nuclei has been seen during the last years. We express the total Hamiltonian H of interacting fields through new operators for particle creation and destruction and show that this can be understood as a UT of H. The respective particles may be called “cothed”. They are identified with the physical particles. The Hamiltonian in the new form turns out to be dependent on the renormalized particle masses and not the “bare” ones. Forms of the same kind are derived for all the Poincaré group generators. By using this new form of the Hamiltonian we suggest an approach to the bound state problem in relativistic quantum field theories (RQFT's). We also discuss applications of the developed formalism for constructing effective electromagnetic (EM) currents in the theory of photonuclear reactions.
Gillani, S. S. A.; Shah, H. A.; Tsintsadze, N. L.; Razzaq, M.
2010-08-15
It is shown that the interaction of the superstrong laser radiation with an isotropic plasma leads to the generation of low frequency electromagnetic (EM) waves and in particular a quasistationary magnetic field. When the relativistic circularly polarized transverse EM wave propagates along z-axis, it creates a ponderomotive force, which affects the motion of particles along the direction of its propagation. On the other hand, motion of the particles across the direction of propagation is defined by the ponderomotive potential. The dispersion relation for the transverse EM wave using a special distribution function, which has an anisotropic form, is derived. The dispersion relation is subsequently investigated for a number of special cases. In general, it is shown that the growth rate of the EM wave strongly depends upon its intensity.
Particle Production in Strong Electromagnetic Fields in Relativistic Heavy-Ion Collisions
Tuchin, Kirill
2013-01-01
I reviewmore » the origin and properties of electromagnetic fields produced in heavy-ion collisions. The field strength immediately after a collision is proportional to the collision energy and reaches ~ m π 2 at RHIC and ~ 10 m π 2 at LHC. I demonstrate by explicit analytical calculation that after dropping by about one-two orders of magnitude during the first fm/c of plasma expansion, it freezes out and lasts for as long as quark-gluon plasma lives as a consequence of finite electrical conductivity of the plasma. Magnetic field breaks spherical symmetry in the direction perpendicular to the reaction plane, and therefore all kinetic coefficients are anisotropic. I examine viscosity of QGP and show that magnetic field induces azimuthal anisotropy on plasma flow even in spherically symmetric geometry. Very strong electromagnetic field has an important impact on particle production. I discuss the problem of energy loss and polarization of fast fermions due to synchrotron radiation, consider photon decay induced by magnetic field, elucidate J / ψ dissociation via Lorentz ionization mechanism, and examine electromagnetic radiation by plasma. I conclude that all processes in QGP are affected by strong electromagnetic field and call for experimental investigation.« less
An electromagnetic calorimeter for the solenoidal tracker at the Relativistic Heavy Ion Collider
Westfall, G.D.; Llope, W.J.; Underwood, D.G.
1993-10-01
In this document, we outline a proposal to the National Science Foundation (NSF) for the construction of an electromagnetic calorimeter for STAR that fulfills these requirements. This proposal creates the opportunity for the NSF to make a major impact on the experimental program at RHIC by providing a crucial, but defensibly omitted, component of the STAR experiment as approved.
The Relativistic Transformation for an Electromagnetic Plane Wave with General Time Dependence
ERIC Educational Resources Information Center
Smith, Glenn S.
2012-01-01
In special relativity, the transformation between inertial frames for an electromagnetic plane wave is usually derived for the time-harmonic case (the field is a sinusoid of infinite duration), even though all practical waves are of finite duration and may not even contain a dominant sinusoid. This paper presents an alternative derivation in which…
NASA Astrophysics Data System (ADS)
Cooray, G. V.; Cooray, G. K.
2011-12-01
Gurevich et al. [1] postulated that the source of narrow bipolar pulses, a class of high energy pulses that occur during thunderstorms, could be a runaway electron avalanche driven by the intense electric fields of a thunderstorm. Recently, Watson and Marshall [2] used the modified transmission line model to test the mechanism of the source of narrow bipolar pulses. In a recent paper, Cooray and Cooray [3] demonstrated that the electromagnetic fields of accelerating charges could be used to evaluate the electromagnetic fields from electrical discharges if the temporal and spatial variation of the charges in the discharge is known. In the present study, those equations were utilized to evaluate the electromagnetic fields generated by a relativistic electron avalanche. In the analysis it is assumed that all the electrons in the avalanche are moving with the same speed. In other words, the growth or the decay of the number of electrons takes place only at the head of the avalanche. It is shown that the radiation is emanating only from the head of the avalanche where electrons are being accelerated. It is also shown that an analytical expression for the radiation field of the avalanche at any distance can be written directly in terms of the e-folding length of the avalanche. This makes it possible to extract directly the spatial variation of the e-folding length of the avalanche from the measured radiation fields. In the study this model avalanche was used to investigate whether it can be used to describe the measured electromagnetic fields of narrow bipolar pulses. The results obtained are in reasonable agreement with the two station data of Eack [4] for speeds of propagation around (2 - 2.5) x 10^8 m/s and when the propagation effects on the electric fields measured at the distant station is taken into account. [1] Gurevich et al. (2004), Phys. Lett. A., 329, pp. 348 -361. [2] Watson, S. S. and T. C. Marshall (2007), Geophys. Res. Lett., Vol. 34, L04816, doi: 10
Mankiewicz, L. ); Sawicki, M. )
1989-11-15
Within a relativistically correct yet analytically solvable model of light-front quantum mechanics we construct the electromagnetic form factor of the two-body bound state and we study the validity of the static approximation to the full form factor. Upon comparison of full form factors calculated for different values of binding energy we observe an unexpected effect that for very strongly bound states further increase in binding leads to an increase in the size of the bound system. A similar effect is found for another quantum-mechanical model of relativistic dynamics.
Electromagnetic dissociation of relativistic sup 28 Si into p + sup 27 Al
Barrette, J.; Bellwied, R.; Braun-Munzinger, P.; Cleland, W.E.; David, G.; Dee, J.; Dietzsch, O.; Duek, E.; Fatyga, M.; Fox, D.; Greene, S.V.; Hall, J.R.; Hemmick, T.K.; Herrmann, N.; Hogue, R.W.; Hong, B.; Jayananda, K.; Kraus, D.; Shiva Kumar, B.; Lacasse, R.; Lissauer, D.; Llope, W.J.; Ludlam, T.; Majka, R.; Makowiecki, D.; Mark, S.K.; McCorkle, S.; Mitchell, J.T.; Muthuswamy, M.; O'Brien, E.; Polychronakos, V.; Pruneau, C.; Rotondo, F.S.; Sandweiss, J.; Simon-Gillo, J.; Sonnadara, U.; Stachel, J.; Takai, H.; Takagui, E.M.; Throwe, T.G.; Waters, L.; Willis, W.J.; Winter, C.; Wolf, K.; Wolfe, D.; Woody, C.L.; Xu, N.; Zhang, Y.; Zhang, Z.; Zou, Z. European Organization for Nuclear Research , Geneva Los Alamos National Laboratory, Los Alamos, New Mexico 87545 McGill University, Montreal, Quebec, H3A 2K6 University of New Mexico, Albuquerque, New Mexico 87131
1992-05-01
We report a direct measurement of the final-state energy spectrum in the electromagnetic dissociation of {sup 28}Si into {ital p}+ {sup 27}Al at an energy of 14.6 GeV/nucleon. The final-state energy is obtained through a calculation of the {ital p}-{sup 27}Al invariant mass in kinematically reconstructed events. The final-state energy spectrum for all targets is peaked near the isovector giant-dipole resonance in {sup 28}Si and the dependence of the magnitude of the cross section on target charge confirms that the excitation is largely electromagnetic. By exploiting the expected scaling behavior on target {ital Z} and {ital A}, the background from nuclear interactions is evaluated and subtracted, leaving a pure electromagnetic dissociation final-state energy distribution. This distribution is well reproduced by simulated events, in which the photon spectrum calculated in the Weiszaecker-Williams approximation is combined with experimental data on the photonuclear reaction {sup 28}Si({gamma},{ital p}) {sup 27}Al, and slight differences are observed only at low final-state energy.
Beam Optics Analysis - An Advanced 3D Trajectory Code
Ives, R. Lawrence; Bui, Thuc; Vogler, William; Neilson, Jeff; Read, Mike; Shephard, Mark; Bauer, Andrew; Datta, Dibyendu; Beal, Mark
2006-01-03
Calabazas Creek Research, Inc. has completed initial development of an advanced, 3D program for modeling electron trajectories in electromagnetic fields. The code is being used to design complex guns and collectors. Beam Optics Analysis (BOA) is a fully relativistic, charged particle code using adaptive, finite element meshing. Geometrical input is imported from CAD programs generating ACIS-formatted files. Parametric data is inputted using an intuitive, graphical user interface (GUI), which also provides control of convergence, accuracy, and post processing. The program includes a magnetic field solver, and magnetic information can be imported from Maxwell 2D/3D and other programs. The program supports thermionic emission and injected beams. Secondary electron emission is also supported, including multiple generations. Work on field emission is in progress as well as implementation of computer optimization of both the geometry and operating parameters. The principle features of the program and its capabilities are presented.
Beam Optics Analysis — An Advanced 3D Trajectory Code
NASA Astrophysics Data System (ADS)
Ives, R. Lawrence; Bui, Thuc; Vogler, William; Neilson, Jeff; Read, Mike; Shephard, Mark; Bauer, Andrew; Datta, Dibyendu; Beal, Mark
2006-01-01
Calabazas Creek Research, Inc. has completed initial development of an advanced, 3D program for modeling electron trajectories in electromagnetic fields. The code is being used to design complex guns and collectors. Beam Optics Analysis (BOA) is a fully relativistic, charged particle code using adaptive, finite element meshing. Geometrical input is imported from CAD programs generating ACIS-formatted files. Parametric data is inputted using an intuitive, graphical user interface (GUI), which also provides control of convergence, accuracy, and post processing. The program includes a magnetic field solver, and magnetic information can be imported from Maxwell 2D/3D and other programs. The program supports thermionic emission and injected beams. Secondary electron emission is also supported, including multiple generations. Work on field emission is in progress as well as implementation of computer optimization of both the geometry and operating parameters. The principle features of the program and its capabilities are presented.
NASA Technical Reports Server (NTRS)
Gabriel, G. J.
1980-01-01
A new theorem on a class of four-dimensional skew-symmetric tensors is demonstrated. Coupled with the relativistic covariant form of Maxwell's equations, this theorem consolidates the classifications of guided waves by combining the three types - TE, TM, TEM - under a uniform condition applied to the generating four-potential which is Lorentz invariant. Each type corresponds to a potential of which a pair of the four components vanishes in a particular frame. Through appropriate normalization conditions, the resulting time-domain equations for the field amplitudes are readily reduced to modified telegraphist equations, which in turn lead to distributed network representations for each of the three types. The ambiguity of distributed network formalisms in general is elucidated and the concept of network parameter densities such as traditionally employed in TEM transmission line theory is questioned.
Relativistic Plasma Polarizer: Impact of Temperature Anisotropy on Relativistic Transparency.
Stark, David J; Bhattacharjee, Chinmoy; Arefiev, Alexey V; Toncian, Toma; Hazeltine, R D; Mahajan, S M
2015-07-10
3D particle-in-cell simulations demonstrate that the enhanced transparency of a relativistically hot plasma is sensitive to how the energy is partitioned between different degrees of freedom. For an anisotropic electron distribution, propagation characteristics, like the critical density, will depend on the polarization of the electromagnetic wave. Despite the onset of the Weibel instability in such plasmas, the anisotropy can persist long enough to affect laser propagation. This plasma can then function as a polarizer or a wave plate to dramatically alter the pulse polarization. PMID:26207473
NASA Technical Reports Server (NTRS)
Gamayunov, K. V.; Khazanov, G. V.
2007-01-01
We consider the effect of oblique EMIC waves on relativistic electron scattering in the outer radiation belt using simultaneous observations of plasma and wave parameters from CRRES. The main findings can be s ummarized as follows: 1. In 1comparison with field-aligned waves, int ermediate and highly oblique distributions decrease the range of pitc h-angles subject to diffusion, and reduce the local scattering rate b y an order of magnitude at pitch-angles where the principle absolute value of n = 1 resonances operate. Oblique waves allow the absolute va lue of n > 1 resonances to operate, extending the range of local pitc h-angle diffusion down to the loss cone, and increasing the diffusion at lower pitch angles by orders of magnitude; 2. The local diffusion coefficients derived from CRRES data are qualitatively similar to the local results obtained for prescribed plasma/wave parameters. Conseq uently, it is likely that the bounce-averaged diffusion coefficients, if estimated from concurrent data, will exhibit the dependencies similar to those we found for model calculations; 3. In comparison with f ield-aligned waves, intermediate and highly oblique waves decrease th e bounce-averaged scattering rate near the edge of the equatorial lo ss cone by orders of magnitude if the electron energy does not excee d a threshold (approximately equal to 2 - 5 MeV) depending on specified plasma and/or wave parameters; 4. For greater electron energies_ ob lique waves operating the absolute value of n > 1 resonances are more effective and provide the same bounce_averaged diffusion rate near the loss cone as fiel_aligned waves do.
Attosecond Electro-Magnetic Forces Acting on Metal Nanospheres Induced By Relativistic Electrons
NASA Astrophysics Data System (ADS)
Lagos, M. J.; Batson, P. E.; Reyes-Coronado, A.; Echenique, P. M.; Aizpurua, J.
2014-03-01
Swift electron scattering near nanoscale materials provides information about light-matter behavior, including induced forces. We calculate time-dependent electromagnetic forces acting on 1-1.5 nm metal nanospheres induced by passing swift electrons, finding both impulse-like and oscillatory response forces. Initially, impulse-like forces are generated by a competition between attractive electric forces and repulsive magnetic forces, lasting a few attoseconds (5-10 as). Oscillatory, plasmonic response forces take place later in time, last a few femtoseconds (1- 5 fs), and apparently rely on photon emission by decay of the electron-induced surface plasmons. A comparison of the strength of these two forces suggests that the impulse-like behavior dominates the process, and can transfer significant linear momentum to the sphere. Our results advance understanding of the physics behind the observation of both attractive and repulsive behavior of gold nano-particles induced by electron beams in aberration-corrected electron microscopy. Work supported under DOE, Award # DE-SC0005132, Basque Gov. project ETORTEK inano, Spanish Ministerio de Ciencia e Innovacion, No. FIS2010-19609-C02-01.
Electromagnetic fluctuations in magnetized plasmas. I. The rigorous relativistic kinetic theory
Schlickeiser, R. E-mail: yoonp@umd.edu; Yoon, P. H. E-mail: yoonp@umd.edu
2015-07-15
Using the system of the Klimontovich and Maxwell equations, the general linear fluctuation theory for magnetized plasmas is developed. General expressions for the electromagnetic fluctuation spectra (electric and magnetic fields) from uncorrelated plasma particles in plasmas with a uniform magnetic field are derived, which are covariantly correct within the theory of special relativity. The general fluctuation spectra hold for plasmas of arbitrary composition, arbitrary momentum dependences of the plasma particle distribution functions, and arbitrary orientations of the wave vector with respect to the uniform magnetic field. Moreover, no restrictions on the values of the real and the imaginary parts of the frequency are made. The derived fluctuation spectra apply to both non-collective fluctuations and collective plasma eigenmodes in magnetized plasmas. In the latter case, kinetic equations for the components of fluctuating electric and magnetic fields in magnetized plasmas are derived that include the effect of spontaneous emission and absorption. In the limiting case of an unmagnetized plasmas, the general fluctuation spectra correctly reduce to the unmagnetized fluctuation spectra derived before.
NASA Astrophysics Data System (ADS)
Tiberi, Gianluigi; Fontana, Nunzia; Monorchio, Agostino; Stara, Riccardo; Retico, Alessandra; Tosetti, Michela
2015-12-01
A procedure for evaluating radio-frequency electromagnetic fields in anatomical human models for any matching and coupling conditions is introduced. The procedure resorts to the extraction of basis functions: such basis functions, which represent the fields produced by each individual port without any residual coupling, are derived through an algebraic procedure which uses the S parameter matrix and the fields calculated in one (only) full-wave simulation. The basis functions are then used as building-blocks for calculating the fields for any other S parameter matrix. The proposed approach can be used both for volume coil driven in quadrature and for parallel transmission configuration.
Three-dimensional parallel UNIPIC-3D code for simulations of high-power microwave devices
NASA Astrophysics Data System (ADS)
Wang, Jianguo; Chen, Zaigao; Wang, Yue; Zhang, Dianhui; Liu, Chunliang; Li, Yongdong; Wang, Hongguang; Qiao, Hailiang; Fu, Meiyan; Yuan, Yuan
2010-07-01
This paper introduces a self-developed, three-dimensional parallel fully electromagnetic particle simulation code UNIPIC-3D. In this code, the electromagnetic fields are updated using the second-order, finite-difference time-domain method, and the particles are moved using the relativistic Newton-Lorentz force equation. The electromagnetic field and particles are coupled through the current term in Maxwell's equations. Two numerical examples are used to verify the algorithms adopted in this code, numerical results agree well with theoretical ones. This code can be used to simulate the high-power microwave (HPM) devices, such as the relativistic backward wave oscillator, coaxial vircator, and magnetically insulated line oscillator, etc. UNIPIC-3D is written in the object-oriented C++ language and can be run on a variety of platforms including WINDOWS, LINUX, and UNIX. Users can use the graphical user's interface to create the complex geometric structures of the simulated HPM devices, which can be automatically meshed by UNIPIC-3D code. This code has a powerful postprocessor which can display the electric field, magnetic field, current, voltage, power, spectrum, momentum of particles, etc. For the sake of comparison, the results computed by using the two-and-a-half-dimensional UNIPIC code are also provided for the same parameters of HPM devices, the numerical results computed from these two codes agree well with each other.
Electric dipole transitions for 3d64s-3d64p in Mn I
NASA Astrophysics Data System (ADS)
Kabakçı, Selda; Özdemir, Leyla; Usta, Betül Karaçoban
2015-10-01
We have calculated the logarithmic weighted oscillator strengths and transition probabilities (or rates) for 3d64s-3d64p electric dipole transitions in neutral manganese (Mn I, Z=25) by using two configuration interaction methods (the multiconfiguration Hartree-Fock (MCHF) method within the framework of Breit-Pauli relativistic corrections developed by Fischer and Cowan's relativistic Hartree-Fock (HFR) method). Results obtained have been compared with other calculations and experiments.
Zhou Jun; Sebastian, Evelyn; Mangona, Victor; Yan Di
2013-02-15
Purpose: In order to increase the accuracy and speed of catheter reconstruction in a high-dose-rate (HDR) prostate implant procedure, an automatic tracking system has been developed using an electromagnetic (EM) device (trakSTAR, Ascension Technology, VT). The performance of the system, including the accuracy and noise level with various tracking parameters and conditions, were investigated. Methods: A direct current (dc) EM transmitter (midrange model) and a sensor with diameter of 1.3 mm (Model 130) were used in the trakSTAR system for tracking catheter position during HDR prostate brachytherapy. Localization accuracy was assessed under both static and dynamic analyses conditions. For the static analysis, a calibration phantom was used to investigate error dependency on operating room (OR) table height (bottom vs midposition vs top), sensor position (distal tip of catheter vs connector end of catheter), direction [left-right (LR) vs anterior-posterior (AP) vs superior-inferior (SI)], sampling frequency (40 vs 80 vs 120 Hz), and interference from OR equipment (present vs absent). The mean and standard deviation of the localization offset in each direction and the corresponding error vectors were calculated. For dynamic analysis, the paths of five straight catheters were tracked to study the effects of directions, sampling frequency, and interference of EM field. Statistical analysis was conducted to compare the results in different configurations. Results: When interference was present in the static analysis, the error vectors were significantly higher at the top table position (3.3 {+-} 1.3 vs 1.8 {+-} 0.9 mm at bottom and 1.7 {+-} 1.0 mm at middle, p < 0.001), at catheter end position (3.1 {+-} 1.1 vs 1.4 {+-} 0.7 mm at the tip position, p < 0.001), and at 40 Hz sampling frequency (2.6 {+-} 1.1 vs 2.4 {+-} 1.5 mm at 80 Hz and 1.8 {+-} 1.1 at 160 Hz, p < 0.001). So did the mean offset errors in the LR direction (-1.7 {+-} 1.4 vs 0.4 {+-} 0.5 mm in AP and 0
NASA Astrophysics Data System (ADS)
Cooray, Vernon; Cooray, Gerald; Marshall, Thomas; Arabshahi, Shahab; Dwyer, Joseph; Rassoul, Hamid
2014-11-01
In the present study, electromagnetic fields of accelerating charges were utilized to evaluate the electromagnetic fields generated by a relativistic electron avalanche. In the analysis it is assumed that all the electrons in the avalanche are moving with the same speed. In other words, the growth or the decay of the number of electrons takes place only at the head of the avalanche. It is shown that the radiation is emanating only from the head of the avalanche where electrons are being accelerated. It is also shown that an analytical expression for the radiation field of the avalanche at any distance can be written directly in terms of the e-folding length of the avalanche. This model of the avalanche was utilized to test the idea whether the source of the lightning signatures known as narrow bipolar pulses could be relativistic avalanches. The idea was tested by using the simultaneously measured electric fields of narrow bipolar pulses at two distances, one measured far away from the source and the other in the near vicinity. The avalanche parameters were extracted from the distant field and they are used to evaluate the close field. The results show that the source of the NBP can be modeled either as a single or a multiple burst of relativistic avalanches with speed of avalanches in the range of 2-3 × 108 m/s. The multiple avalanche model agrees better with the experimental data in that it can also generate the correct signature of the time derivatives and the HF and VHF radiation bursts of NBP.
Parallel 3-D particle-in-cell modelling of charged ultrarelativistic beam dynamics
NASA Astrophysics Data System (ADS)
Boronina, Marina A.; Vshivkov, Vitaly A.
2015-12-01
> ) in supercolliders. We use the 3-D set of Maxwell's equations for the electromagnetic fields, and the Vlasov equation for the distribution function of the beam particles. The model incorporates automatically the longitudinal effects, which can play a significant role in the cases of super-high densities. We present numerical results for the dynamics of two focused ultrarelativistic beams with a size ratio 10:1:100. The results demonstrate high efficiency of the proposed computational methods and algorithms, which are applicable to a variety of problems in relativistic plasma physics.
NASA Astrophysics Data System (ADS)
Pletinckx, D.
2011-09-01
The current 3D hype creates a lot of interest in 3D. People go to 3D movies, but are we ready to use 3D in our homes, in our offices, in our communication? Are we ready to deliver real 3D to a general public and use interactive 3D in a meaningful way to enjoy, learn, communicate? The CARARE project is realising this for the moment in the domain of monuments and archaeology, so that real 3D of archaeological sites and European monuments will be available to the general public by 2012. There are several aspects to this endeavour. First of all is the technical aspect of flawlessly delivering 3D content over all platforms and operating systems, without installing software. We have currently a working solution in PDF, but HTML5 will probably be the future. Secondly, there is still little knowledge on how to create 3D learning objects, 3D tourist information or 3D scholarly communication. We are still in a prototype phase when it comes to integrate 3D objects in physical or virtual museums. Nevertheless, Europeana has a tremendous potential as a multi-facetted virtual museum. Finally, 3D has a large potential to act as a hub of information, linking to related 2D imagery, texts, video, sound. We describe how to create such rich, explorable 3D objects that can be used intuitively by the generic Europeana user and what metadata is needed to support the semantic linking.
NASA Astrophysics Data System (ADS)
Bulanov, Sergei V.; Esirkepov, Timur Z.; Hayashi, Yukio; Kando, Masaki; Kiriyama, Hiromitsu; Koga, James K.; Kondo, Kiminori; Kotaki, Hideyuki; Pirozhkov, Alexander S.; Bulanov, Stepan S.; Zhidkov, Alexei G.; Chen, Pisin; Neely, David; Kato, Yoshiaki; Narozhny, Nikolay B.; Korn, Georg
2011-06-01
The critical electric field of quantum electrodynamics, called also the Schwinger field, is so strong that it produces electron-positron pairs from vacuum, converting the energy of light into matter. Since the dawn of quantum electrodynamics, there has been a dream on how to reach it on Earth. With the rise of laser technology this field has become feasible through the construction of extremely high power lasers or/and with the sophisticated use of nonlinear processes in relativistic plasmas. This is one of the most attractive motivations for extremely high power laser development, i.e. producing matter from vacuum by pure light in fundamental process of quantum electrodynamics in the nonperturbative regime. Recently it has been realized that a laser with intensity well below the Schwinger limit can create an avalanche of electron-positron pairs similar to a discharge before attaining the Schwinger field. It has also been realized that the Schwinger limit can be reached using an appropriate configuration of laser beams. In experiments on the collision of laser light and high intensity electromagnetic pulses generated by relativistic flying mirrors, with electron bunches produced by a conventional accelerator and with laser wake field accelerated electrons the studying of extreme field limits in the nonlinear interaction of electromagnetic waves is proposed. The regimes of dominant radiation reaction, which completely changes the electromagnetic wave-matter interaction, will be revealed. This will result in a new powerful source of high brightness gamma-rays. A possibility of the demonstration of the electronpositron pair creation in vacuum via multi-photon processes can be realized. This will allow modeling under terrestrial laboratory conditions neutron star magnetospheres, cosmological gamma ray bursts and the Leptonic Era of the Universe.
3d-3d correspondence revisited
NASA Astrophysics Data System (ADS)
Chung, Hee-Joong; Dimofte, Tudor; Gukov, Sergei; Sułkowski, Piotr
2016-04-01
In fivebrane compactifications on 3-manifolds, we point out the importance of all flat connections in the proper definition of the effective 3d {N}=2 theory. The Lagrangians of some theories with the desired properties can be constructed with the help of homological knot invariants that categorify colored Jones polynomials. Higgsing the full 3d theories constructed this way recovers theories found previously by Dimofte-Gaiotto-Gukov. We also consider the cutting and gluing of 3-manifolds along smooth boundaries and the role played by all flat connections in this operation.
An, Zhenguo; Zhang, Jingjie; Pan, Shunlong
2010-04-14
Novel 3D urchin-like glass/Ni-Ni(3)P/Co(2)P(2)O(7) core/shell/shell composite hollow structures are fabricated for the first time by controlled stepwise assembly of granular Ni-Ni(3)P alloy and ribbon-like Co(2)P(2)O(7) nanocrystals on hollow glass spheres in aqueous solutions at mild conditions. It is found that the shell structure and the overall morphology of the products can be tailored by properly tuning the annealing temperature. The as-obtained composite core/shell/shell products possess low density (ca. 1.18 g cm(-3)) and shape-dependent magnetic and microwave absorbing properties, and thus may have some promising applications in the fields of low-density magnetic materials, microwave absorbers, etc. Based on a series of contrast experiments, the probable formation mechanism of the core/shell/shell hierarchical structures is proposed. This work provides an additional strategy to prepare core/shell composite spheres with tailored shell morphology and electromagnetic properties. PMID:20379530
NASA Astrophysics Data System (ADS)
Meulien Ohlmann, Odile
2013-02-01
Today the industry offers a chain of 3D products. Learning to "read" and to "create in 3D" becomes an issue of education of primary importance. 25 years professional experience in France, the United States and Germany, Odile Meulien set up a personal method of initiation to 3D creation that entails the spatial/temporal experience of the holographic visual. She will present some different tools and techniques used for this learning, their advantages and disadvantages, programs and issues of educational policies, constraints and expectations related to the development of new techniques for 3D imaging. Although the creation of display holograms is very much reduced compared to the creation of the 90ies, the holographic concept is spreading in all scientific, social, and artistic activities of our present time. She will also raise many questions: What means 3D? Is it communication? Is it perception? How the seeing and none seeing is interferes? What else has to be taken in consideration to communicate in 3D? How to handle the non visible relations of moving objects with subjects? Does this transform our model of exchange with others? What kind of interaction this has with our everyday life? Then come more practical questions: How to learn creating 3D visualization, to learn 3D grammar, 3D language, 3D thinking? What for? At what level? In which matter? for whom?
NASA Technical Reports Server (NTRS)
Norbury, John W.
1989-01-01
Previous analyses of the comparison of Weizsacker-Williams (WW) theory to experiment for nucleon emission via electromagnetic (EM) excitations in nucleus-nucleus collisions were not definitive because of different assumptions concerning the value of the minimum impact parameter. This situation is corrected by providing criteria that allows definitive statements to be made concerning agreement or disagreement between WW theory and experiment.
ERIC Educational Resources Information Center
Hastings, S. K.
2002-01-01
Discusses 3 D imaging as it relates to digital representations in virtual library collections. Highlights include X-ray computed tomography (X-ray CT); the National Science Foundation (NSF) Digital Library Initiatives; output peripherals; image retrieval systems, including metadata; and applications of 3 D imaging for libraries and museums. (LRW)
Radiation from Relativistic Jets in Turbulent Magnetic Fields
NASA Technical Reports Server (NTRS)
Nishikawa, K.-I.; Mizuno, Y.; Medvedev, M.; Zhang, B.; Hardee, P.; Niemiec, J.; Nordlund, A.; Frederiksen, J.; Mizuno, Y.; Sol, H.; Fishman, G. J.
2008-01-01
Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we have investigated long-term particle acceleration associated with an relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. The simulations have been performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. The acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value. Behind the bow shock in the jet shock strong electromagnetic fields are generated. These fields may lead to the afterglow emission. We have calculated the time evolution of the spectrum from two electrons propagating in a uniform parallel magnetic field to verify the technique.
Crandall, K.R.
1987-08-01
TRACE 3-D is an interactive beam-dynamics program that calculates the envelopes of a bunched beam, including linear space-charge forces, through a user-defined transport system. TRACE 3-D provides an immediate graphics display of the envelopes and the phase-space ellipses and allows nine types of beam-matching options. This report describes the beam-dynamics calculations and gives detailed instruction for using the code. Several examples are described in detail.
NASA Astrophysics Data System (ADS)
Haba, Z.
2009-02-01
We discuss relativistic diffusion in proper time in the approach of Schay (Ph.D. thesis, Princeton University, Princeton, NJ, 1961) and Dudley [Ark. Mat. 6, 241 (1965)]. We derive (Langevin) stochastic differential equations in various coordinates. We show that in some coordinates the stochastic differential equations become linear. We obtain momentum probability distribution in an explicit form. We discuss a relativistic particle diffusing in an external electromagnetic field. We solve the Langevin equations in the case of parallel electric and magnetic fields. We derive a kinetic equation for the evolution of the probability distribution. We discuss drag terms leading to an equilibrium distribution. The relativistic analog of the Ornstein-Uhlenbeck process is not unique. We show that if the drag comes from a diffusion approximation to the master equation then its form is strongly restricted. The drag leading to the Tsallis equilibrium distribution satisfies this restriction whereas the one of the Jüttner distribution does not. We show that any function of the relativistic energy can be the equilibrium distribution for a particle in a static electric field. A preliminary study of the time evolution with friction is presented. It is shown that the problem is equivalent to quantum mechanics of a particle moving on a hyperboloid with a potential determined by the drag. A relation to diffusions appearing in heavy ion collisions is briefly discussed.
Haba, Z
2009-02-01
We discuss relativistic diffusion in proper time in the approach of Schay (Ph.D. thesis, Princeton University, Princeton, NJ, 1961) and Dudley [Ark. Mat. 6, 241 (1965)]. We derive (Langevin) stochastic differential equations in various coordinates. We show that in some coordinates the stochastic differential equations become linear. We obtain momentum probability distribution in an explicit form. We discuss a relativistic particle diffusing in an external electromagnetic field. We solve the Langevin equations in the case of parallel electric and magnetic fields. We derive a kinetic equation for the evolution of the probability distribution. We discuss drag terms leading to an equilibrium distribution. The relativistic analog of the Ornstein-Uhlenbeck process is not unique. We show that if the drag comes from a diffusion approximation to the master equation then its form is strongly restricted. The drag leading to the Tsallis equilibrium distribution satisfies this restriction whereas the one of the Jüttner distribution does not. We show that any function of the relativistic energy can be the equilibrium distribution for a particle in a static electric field. A preliminary study of the time evolution with friction is presented. It is shown that the problem is equivalent to quantum mechanics of a particle moving on a hyperboloid with a potential determined by the drag. A relation to diffusions appearing in heavy ion collisions is briefly discussed. PMID:19391727
Sur, Chiranjib; Chaudhuri, Rajat K.
2007-07-15
We report the forbidden 6s{sup 2} {sup 1}S{sub 0}{yields}6s5d {sup 3}D{sub 1} magnetic-dipole transition amplitude computed using multireference Fock-space coupled-cluster theory. Our computed transition matrix element (1.34x10{sup -4}){mu}{sub B} is in excellent agreement with the experimental value (1.33x10{sup -4}){mu}{sub B}. This value in combination with other known quantities will be helpful in determining the parity-nonconserving amplitude for the 6s{sup 2} {sup 1}S{sub 0}{yields}6s5d {sup 3}D{sub 1} transition in atomic Yb. To our knowledge, this calculation is the most accurate to date, and can be very important in the search for physics beyond the standard model. We further report the 6s6p {sup 3}P{sub 0}{yields}6s6p {sup 1}P{sub 1} and 6s5d {sup 3}D{sub 1}{yields}6s6p {sup 3}P{sub 0} transition matrix elements, which are also in good agreement with the earlier theoretical estimates.
NASA Astrophysics Data System (ADS)
Oldham, Mark
2015-01-01
Radiochromic materials exhibit a colour change when exposed to ionising radiation. Radiochromic film has been used for clinical dosimetry for many years and increasingly so recently, as films of higher sensitivities have become available. The two principle advantages of radiochromic dosimetry include greater tissue equivalence (radiologically) and the lack of requirement for development of the colour change. In a radiochromic material, the colour change arises direct from ionising interactions affecting dye molecules, without requiring any latent chemical, optical or thermal development, with important implications for increased accuracy and convenience. It is only relatively recently however, that 3D radiochromic dosimetry has become possible. In this article we review recent developments and the current state-of-the-art of 3D radiochromic dosimetry, and the potential for a more comprehensive solution for the verification of complex radiation therapy treatments, and 3D dose measurement in general.
NASA Astrophysics Data System (ADS)
Yu, Peicheng; Xu, Xinlu; Davidson, Asher; Tableman, Adam; Dalichaouch, Thamine; Li, Fei; Meyers, Michael D.; An, Weiming; Tsung, Frank S.; Decyk, Viktor K.; Fiuza, Frederico; Vieira, Jorge; Fonseca, Ricardo A.; Lu, Wei; Silva, Luis O.; Mori, Warren B.
2016-07-01
When modeling laser wakefield acceleration (LWFA) using the particle-in-cell (PIC) algorithm in a Lorentz boosted frame, the plasma is drifting relativistically at βb c towards the laser, which can lead to a computational speedup of ∼ γb2 = (1 - βb2)-1. Meanwhile, when LWFA is modeled in the quasi-3D geometry in which the electromagnetic fields and current are decomposed into a limited number of azimuthal harmonics, speedups are achieved by modeling three dimensional (3D) problems with the computational loads on the order of two dimensional r - z simulations. Here, we describe a method to combine the speedups from the Lorentz boosted frame and quasi-3D algorithms. The key to the combination is the use of a hybrid Yee-FFT solver in the quasi-3D geometry that significantly mitigates the Numerical Cerenkov Instability (NCI) which inevitably arises in a Lorentz boosted frame due to the unphysical coupling of Langmuir modes and EM modes of the relativistically drifting plasma in these simulations. In addition, based on the space-time distribution of the LWFA data in the lab and boosted frame, we propose to use a moving window to follow the drifting plasma, instead of following the laser driver as is done in the LWFA lab frame simulations, in order to further reduce the computational loads. We describe the details of how the NCI is mitigated for the quasi-3D geometry, the setups for simulations which combine the Lorentz boosted frame, quasi-3D geometry, and the use of a moving window, and compare the results from these simulations against their corresponding lab frame cases. Good agreement is obtained among these sample simulations, particularly when there is no self-trapping, which demonstrates it is possible to combine the Lorentz boosted frame and the quasi-3D algorithms when modeling LWFA. We also discuss the preliminary speedups achieved in these sample simulations.
NASA Astrophysics Data System (ADS)
Iliesiu, Luca; Kos, Filip; Poland, David; Pufu, Silviu S.; Simmons-Duffin, David; Yacoby, Ran
2016-03-01
We study the conformal bootstrap for a 4-point function of fermions < ψψψψ> in 3D. We first introduce an embedding formalism for 3D spinors and compute the conformal blocks appearing in fermion 4-point functions. Using these results, we find general bounds on the dimensions of operators appearing in the ψ × ψ OPE, and also on the central charge C T . We observe features in our bounds that coincide with scaling dimensions in the GrossNeveu models at large N . We also speculate that other features could coincide with a fermionic CFT containing no relevant scalar operators.
NASA Technical Reports Server (NTRS)
Kocharov, L. G.; Kovaltsov, G. A.; Kocharov, G. E.; Chuikin, E. I.; Usoskin, I. G.; Shea, M. A.; Smart, D. F.; Melnikov, V. F.; Podstrigach, T. S.; Armstrong, T. P.
1994-01-01
Data on X-, gamma-ray, optical and radio emission from the 1991 June 15 solar flare are considered. We have calculated the spectrum of protons that produces gamma-rays during the gradual phase of the flare. The primary proton spectrum can be described as a Bessel-function-type up to 0.8 GeV and a power law with the spectral index approximately equal to 3 from 0.8 up to 10 GeV or above. We have also analyzed data on energetic particles near the Earth. Their spectrum differed from that of primary protons producing gamma-ray line emission. In the gradual phase of the flare additional pulses of energy release occurred and the time profiles of cm-radio emission and gamma-rays in the 0.8-10 MeV energy band and above 50 MeV coincided. A continuous and simultaneous stochastic acceleration of the protons and relativistic electrons at the gradual phase of the flare is considered as a natural explanation of the data.
NASA Astrophysics Data System (ADS)
Lindroth, Eva; Selstø, Sølve; Bengtsson, Jakob
2007-06-01
We solve the time dependent Dirac equation in a basis consisting of the eigenstates of the field free Hamiltonian H0=c α .p + V(r)+ m c^2 β. The dynamics induced by the time dependent perturbation, H'v=c α .A(r,t) in the velocity gauge, or equivalently H'l=r .E(r,t) (1-α .k) in the length gauge, is resolved by solving the first order differential equation arising from the expansion. The number of continuum states needed in order to get converged results is reduced by complex scaling of the coordinates, r ->r e^i θ. We investigate the importance of relativistic effects for various field strengths, E0, and nuclear charges by comparing the solutions of the Dirac equation with those of the Schr"odinger equation. The dynamics is described both within and beyond the dipole approximation (A(r,t) A(t)). For high E0 and high photon energy, φ, it is clear that this approximation breaks down. However, for increasing Z, the electron is increasingly tightly bound, which to some extent reduces the importance of non-dipole effects.
Relativistic Plasma Polarizer: Impact of Temperature Anisotropy on Relativistic Transparency
NASA Astrophysics Data System (ADS)
Hazeltine, R. D.; Stark, David J.; Bhattacharjee, Chinmoy; Arefiev, Alexey V.; Toncian, Toma; Mahajan, S. M.
2015-11-01
3D particle-in-cell simulations demonstrate that the enhanced transparency of a relativistically hot plasma is sensitive to how the energy is partitioned between different degrees of freedom. We consider here the simplest problem: the propagation of a low amplitude pulse through a preformed relativistically hot anisotropic electron plasma to explore its intrinsic dielectric properties. We find that: 1) the critical density for propagation depends strongly on the pulse polarization, 2) two plasmas with the same density and average energy per electron can exhibit profoundly different responses to electromagnetic pulses, 3) the anisotropy-driven Weibel instability develops as expected; the timescales of the growth and back reaction (on anisotropy), however, are long enough that sufficient anisotropy persists for the entire duration of the simulation. This plasma can then function as a polarizer or a wave plate to dramatically alter the pulse polarization. This work was supported by the U.S. DOE Contract Nos. DE-FG02-04ER54742 and DE-AC05-06OR23100 (D. J. S.) and NNSA Contract No. DE-FC52-08NA28512.
NASA Astrophysics Data System (ADS)
Kurkin, S. A.; Hramov, A. E.; Koronovskii, A. A.
2013-07-01
The study of the output power of the electromagnetic radiation of the relativistic electron beam (REB) with virtual cathode in the presence of external magnetic field has been found out. The typical dependencies of the output microwave power of the vircator versus external magnetic field have been analyzed by means of 3D electromagnetic simulation. It has been shown that the power of vircator demonstrates several maxima with external magnetic field growth. The characteristic features of the power behavior are determined by the conditions of the virtual cathode formation in the presence of the external transversal magnetic field and the REB self-magnetic fields.
Kurkin, S. A.; Hramov, A. E.; Koronovskii, A. A.; Saratov State Technical University, Politechnicheskaja 77, Saratov 410028
2013-07-22
The study of the output power of the electromagnetic radiation of the relativistic electron beam (REB) with virtual cathode in the presence of external magnetic field has been found out. The typical dependencies of the output microwave power of the vircator versus external magnetic field have been analyzed by means of 3D electromagnetic simulation. It has been shown that the power of vircator demonstrates several maxima with external magnetic field growth. The characteristic features of the power behavior are determined by the conditions of the virtual cathode formation in the presence of the external transversal magnetic field and the REB self-magnetic fields.
Ghizzo, A.; Reveille, T.; Bertrand, P.; Albrecht-Marc, M.; Johnston, T. W.
2006-10-15
Particle simulations on a flat-topped somewhat underdense (typically n{sub 0}/n{sub c}=0.6) plasma slab by Nikolic et al. [Phys. Rev. E 66, 036404 (2002)] were seen to give transient stimulated scattering behavior with frequency shift [{omega}{sub 0}-{omega}{sub s}({approx_equal}{omega}{sub p})] considerably less than the plasma frequency {omega}{sub p}. This has been linked to the electron acoustic wave (EAW) and the scattering was thus seen as another example of stimulated electron acoustic scattering inferred by Montgomery et al. [Phys. Rev. Lett. 87, 155001 (2001)] from experiments on low-density plasmas. Montgomery et al. had noted the difficulty of how one could have a very narrow observed scattering from a wave whose damping was at least initially very high. Our Vlasov-Maxwell simulations for such somewhat underdense (n{sub 0}/n{sub c}{>=}0.25) plasmas show that the simulation resonance was in fact determined by the beating of the pump with a new 'radiating pseudocavity' electromagnetic mode for the slab at a frequency close to {omega}{sub p} with relatively low loss. This allows the initial narrow-band excitation of the kinetic electrostatic electron nonlinear (KEEN) waves (the nonlinear 'cousins' of EAWs) at a well-defined frequency ({omega}{sub K}{approx_equal}{omega}{sub 0}-{omega}{sub p}<{omega}{sub p}) which is not necessarily the value given by the EAW dispersion relation. (The KEEN wave characteristics have been discussed by Afeyan et al. [33rd AAAC (2003), no. 238, IFSA 2003].) The consideration of such a mechanism is relevant to moderately underdense hot plasmas.
NASA Astrophysics Data System (ADS)
Iizuka, Keigo
2008-02-01
In order to circumvent the fact that only one observer can view the image from a stereoscopic microscope, an attachment was devised for displaying the 3D microscopic image on a large LCD monitor for viewing by multiple observers in real time. The principle of operation, design, fabrication, and performance are presented, along with tolerance measurements relating to the properties of the cellophane half-wave plate used in the design.
Angular description for 3D scattering centers
NASA Astrophysics Data System (ADS)
Bhalla, Rajan; Raynal, Ann Marie; Ling, Hao; Moore, John; Velten, Vincent J.
2006-05-01
The electromagnetic scattered field from an electrically large target can often be well modeled as if it is emanating from a discrete set of scattering centers (see Fig. 1). In the scattering center extraction tool we developed previously based on the shooting and bouncing ray technique, no correspondence is maintained amongst the 3D scattering center extracted at adjacent angles. In this paper we present a multi-dimensional clustering algorithm to track the angular and spatial behaviors of 3D scattering centers and group them into features. The extracted features for the Slicy and backhoe targets are presented. We also describe two metrics for measuring the angular persistence and spatial mobility of the 3D scattering centers that make up these features in order to gather insights into target physics and feature stability. We find that features that are most persistent are also the most mobile and discuss implications for optimal SAR imaging.
Generation of radiation by intense plasma and electromagnetic undulators
Joshi, C.
1991-10-01
We examine the characteristics of the classical radiation emission resulting from the interaction of a relativistic electron beam that propagates perpendicularly through a large amplitude relativistic plasma wave. Such a study is useful for evaluating the feasibility of using relativistic plasma waves as extremely short wavelength undulators for generating short wavelength radiation. The electron trajectories in a plasma wave undulator and in an ac FEL undulator are obtained using perturbation techniques. The spontaneous radiation frequency spectrum and angular distribution emitted by a single electron oscillating in these two undulators are then calculated. The radiation gain of a copropagating electromagnetic wave is calculated. The approximate analytic results for the trajectories, spontaneous radiation and gain are compared with 3-D simulation results. The characteristics of the plasma wave undulator are compared with the ac FEL undulator and linearly polarized magnetic undulator. 50 refs., 26 figs., 3 tabs.
NASA Astrophysics Data System (ADS)
Kostrzewski, Andrew A.; Aye, Tin M.; Kim, Dai Hyun; Esterkin, Vladimir; Savant, Gajendra D.
1998-09-01
Physical Optics Corporation has developed an advanced 3-D virtual reality system for use with simulation tools for training technical and military personnel. This system avoids such drawbacks of other virtual reality (VR) systems as eye fatigue, headaches, and alignment for each viewer, all of which are due to the need to wear special VR goggles. The new system is based on direct viewing of an interactive environment. This innovative holographic multiplexed screen technology makes it unnecessary for the viewer to wear special goggles.
NASA Technical Reports Server (NTRS)
1992-01-01
Ames Research Center research into virtual reality led to the development of the Convolvotron, a high speed digital audio processing system that delivers three-dimensional sound over headphones. It consists of a two-card set designed for use with a personal computer. The Convolvotron's primary application is presentation of 3D audio signals over headphones. Four independent sound sources are filtered with large time-varying filters that compensate for motion. The perceived location of the sound remains constant. Possible applications are in air traffic control towers or airplane cockpits, hearing and perception research and virtual reality development.
Cevidanes, Lucia; Tucker, Scott; Styner, Martin; Kim, Hyungmin; Chapuis, Jonas; Reyes, Mauricio; Proffit, William; Turvey, Timothy; Jaskolka, Michael
2009-01-01
This paper discusses the development of methods for computer-aided jaw surgery. Computer-aided jaw surgery allows us to incorporate the high level of precision necessary for transferring virtual plans into the operating room. We also present a complete computer-aided surgery (CAS) system developed in close collaboration with surgeons. Surgery planning and simulation include construction of 3D surface models from Cone-beam CT (CBCT), dynamic cephalometry, semi-automatic mirroring, interactive cutting of bone and bony segment repositioning. A virtual setup can be used to manufacture positioning splints for intra-operative guidance. The system provides further intra-operative assistance with the help of a computer display showing jaw positions and 3D positioning guides updated in real-time during the surgical procedure. The CAS system aids in dealing with complex cases with benefits for the patient, with surgical practice, and for orthodontic finishing. Advanced software tools for diagnosis and treatment planning allow preparation of detailed operative plans, osteotomy repositioning, bone reconstructions, surgical resident training and assessing the difficulties of the surgical procedures prior to the surgery. CAS has the potential to make the elaboration of the surgical plan a more flexible process, increase the level of detail and accuracy of the plan, yield higher operative precision and control, and enhance documentation of cases. Supported by NIDCR DE017727, and DE018962 PMID:20816308
TE/TM alternating direction scheme for wake field calculation in 3D
NASA Astrophysics Data System (ADS)
Zagorodnov, Igor; Weiland, Thomas
2006-03-01
In the future, accelerators with very short bunches will be used. It demands developing new numerical approaches for long-time calculation of electromagnetic fields in the vicinity of relativistic bunches. The conventional FDTD scheme, used in MAFIA, ABCI and other wake and PIC codes, suffers from numerical grid dispersion and staircase approximation problem. As an effective cure of the dispersion problem, a numerical scheme without dispersion in longitudinal direction can be used as it was shown by Novokhatski et al. [Transition dynamics of the wake fields of ultrashort bunches, TESLA Report 2000-03, DESY, 2000] and Zagorodnov et al. [J. Comput. Phys. 191 (2003) 525]. In this paper, a new economical conservative scheme for short-range wake field calculation in 3D is presented. As numerical examples show, the new scheme is much more accurate on long-time scale than the conventional FDTD approach.
NASA Astrophysics Data System (ADS)
Gil, José J.; San José, Ignacio
2010-11-01
From our previous definition of the indices of polarimetric purity for 3D light beams [J.J. Gil, J.M. Correas, P.A. Melero and C. Ferreira, Monogr. Semin. Mat. G. de Galdeano 31, 161 (2004)], an analysis of their geometric and physical interpretation is presented. It is found that, in agreement with previous results, the first parameter is a measure of the degree of polarization, whereas the second parameter (called the degree of directionality) is a measure of the mean angular aperture of the direction of propagation of the corresponding light beam. This pair of invariant, non-dimensional, indices of polarimetric purity contains complete information about the polarimetric purity of a light beam. The overall degree of polarimetric purity is obtained as a weighted quadratic average of the degree of polarization and the degree of directionality.
Caspi, S.; Helm, M.; Laslett, L.J.
1991-03-30
We have developed an harmonic representation for the three dimensional field components within the windings of accelerator magnets. The form by which the field is presented is suitable for interfacing with other codes that make use of the 3D field components (particle tracking and stability). The field components can be calculated with high precision and reduced cup time at any location (r,{theta},z) inside the magnet bore. The same conductor geometry which is used to simulate line currents is also used in CAD with modifications more readily available. It is our hope that the format used here for magnetic fields can be used not only as a means of delivering fields but also as a way by which beam dynamics can suggest correction to the conductor geometry. 5 refs., 70 figs.
NASA Technical Reports Server (NTRS)
2004-01-01
The Mars Exploration Rover Spirit took this 3-D navigation camera mosaic of the crater called 'Bonneville' after driving approximately 13 meters (42.7 feet) to get a better vantage point. Spirit's current position is close enough to the edge to see the interior of the crater, but high enough and far enough back to get a view of all of the walls. Because scientists and rover controllers are so pleased with this location, they will stay here for at least two more martian days, or sols, to take high resolution panoramic camera images of 'Bonneville' in its entirety. Just above the far crater rim, on the left side, is the rover's heatshield, which is visible as a tiny reflective speck.
Shim3d Helmholtz Solution Package
Energy Science and Technology Software Center (ESTSC)
2009-01-29
This suite of codes solves the Helmholtz Equation for the steady-state propagation of single-frequency electromagnetic radiation in an arbitrary 2D or 3D dielectric medium. Materials can be either transparent or absorptive (including metals) and are described entirely by their shape and complex dielectric constant. Dielectric boundaries are assumed to always fall on grid boundaries and the material within a single grid cell is considered to be uniform. Input to the problem is in the formmore » of a Dirichlet boundary condition on a single boundary, and may be either analytic (Gaussian) in shape, or a mode shape computed using a separate code (such as the included eigenmode solver vwave20), and written to a file. Solution is via the finite difference method using Jacobi iteration for 3D problems or direct matrix inversion for 2D problems. Note that 3D problems that include metals will require different iteration parameters than described in the above reference. For structures with curved boundaries not easily modeled on a rectangular grid, the auxillary codes helmholtz11(2D), helm3d (semivectoral), and helmv3d (full vectoral) are provided. For these codes the finite difference equations are specified on a topological regular triangular grid and solved using Jacobi iteration or direct matrix inversion as before. An automatic grid generator is supplied.« less
NASA Astrophysics Data System (ADS)
Froger, Etienne
1993-05-01
A description of the electromagnetic behavior of a satellite subjected to an electric discharge is given using a specially developed numerical code. One of the particularities of vacuum discharges, obtained by irradiation of polymers, is the intense emission of electrons into the spacecraft environment. Electromagnetic radiation, associated with the trajectories of the particles around the spacecraft, is considered as the main source of the interference observed. In the absence of accurate orbital data and realistic ground tests, the assessment of these effects requires numerical simulation of the interaction between this electron source and the spacecraft. This is done by the GEODE particle code which is applied to characteristic configurations in order to estimate the spacecraft response to a discharge, which is simulated from a vacuum discharge model designed in laboratory. The spacecraft response to a current injection is simulated by the ALICE numerical three dimensional code. The comparison between discharge and injection effects, from the results given by the two codes, illustrates the representativity of electromagnetic susceptibility tests and the main parameters for their definition.
NASA Astrophysics Data System (ADS)
Speck, Jared R.
The two hyperbolic systems of PDEs we consider in this work are the source-free Maxwell-Born-Infeld (MBI) field equations and the Euler-Nordstrom system for gravitationally self-interacting fluids. The former system plays a central role in Kiessling's recently proposed self-consistent model of classical electrodynamics with point charges, a model that does not suffer from the infinities found in the classical Maxwell-Maxwell model with point charges. The latter system is a scalar gravity caricature of the incredibly more complex Euler-Einstein system. The primary original contributions of the thesis can be summarized as follows: (1) We give a sharp non-local criterion for the formation of singularities in plane-symmetric solutions to the source-free MBI field equations. We also use a domain of dependence argument to show that 3-d initial data agreeing with certain plane-symmetric data on a large enough ball lead to solutions that form singularities in finite time. This work is an extension of a theorem of Brenier, who studied singularity formation in periodic plane-symmetric solutions. (2) We prove well-posedness for the Euler-Nordstrom system with a cosmological constant (ENkappa) for initial data that are an HN perturbation (not necessarily small) of a uniform, quiet fluid, for N ≥ 3. The method of proof relies on the framework of energy currents that has been recently developed by Christodoulou. We turn to this method out of necessity: two common frameworks for showing well-posedness in HN, namely symmetric hyperbolicity and strict hyperbolicity, do not apply to the ENkappa system, while Christodoulou's techniques apply to all hyperbolic systems derivable from a Lagrangian, of which the ENkappa system is an example. (3) We insert the speed of light c as a parameter into the ENkappa system (and designate the family of systems ENck ) in order to study the non-relativistic limit c to infinity. Taking the formal limit in the equations gives the Euler
NASA Technical Reports Server (NTRS)
1997-01-01
Many prominent rocks near the Sagan Memorial Station are featured in this image, taken in stereo by the Imager for Mars Pathfinder (IMP) on Sol 3. 3D glasses are necessary to identify surface detail. Wedge is at lower left; Shark, Half-Dome, and Pumpkin are at center. Flat Top, about four inches high, is at lower right. The horizon in the distance is one to two kilometers away.
Mars Pathfinder is the second in NASA's Discovery program of low-cost spacecraft with highly focused science goals. The Jet Propulsion Laboratory, Pasadena, CA, developed and manages the Mars Pathfinder mission for NASA's Office of Space Science, Washington, D.C. JPL is an operating division of the California Institute of Technology (Caltech). The Imager for Mars Pathfinder (IMP) was developed by the University of Arizona Lunar and Planetary Laboratory under contract to JPL. Peter Smith is the Principal Investigator.
Click below to see the left and right views individually. [figure removed for brevity, see original site] Left [figure removed for brevity, see original site] Right
NASA Technical Reports Server (NTRS)
2004-01-01
This 3-D, microscopic imager mosaic of a target area on a rock called 'Diamond Jenness' was taken after NASA's Mars Exploration Rover Opportunity ground into the surface with its rock abrasion tool for a second time.
Opportunity has bored nearly a dozen holes into the inner walls of 'Endurance Crater.' On sols 177 and 178 (July 23 and July 24, 2004), the rover worked double-duty on Diamond Jenness. Surface debris and the bumpy shape of the rock resulted in a shallow and irregular hole, only about 2 millimeters (0.08 inch) deep. The final depth was not enough to remove all the bumps and leave a neat hole with a smooth floor. This extremely shallow depression was then examined by the rover's alpha particle X-ray spectrometer.
On Sol 178, Opportunity's 'robotic rodent' dined on Diamond Jenness once again, grinding almost an additional 5 millimeters (about 0.2 inch). The rover then applied its Moessbauer spectrometer to the deepened hole. This double dose of Diamond Jenness enabled the science team to examine the rock at varying layers. Results from those grindings are currently being analyzed.
The image mosaic is about 6 centimeters (2.4 inches) across.
NASA Technical Reports Server (NTRS)
Mitchell, T. P.
1973-01-01
The motion of a charged particle in electromagnetic fields of various geometric configurations and arising from a variety of sources is of intrinsic interest in electromagnetic theory. The particular configuration consisting of a plane wave propagating in the presence of a static uniform magnetic field whose direction is parallel to the wave normal is examined. The analysis presented here is treated within the context of classical electromagnetic theory. A numerical solution - at least to the approximate Lorentz-Dirac equation - is obtained.
NASA Astrophysics Data System (ADS)
Mediavilla, Evencio; Arribas, Santiago; Roth, Martin; Cepa-Nogué, Jordi; Sánchez, Francisco
2011-09-01
Preface; Acknowledgements; 1. Introductory review and technical approaches Martin M. Roth; 2. Observational procedures and data reduction James E. H. Turner; 3. 3D Spectroscopy instrumentation M. A. Bershady; 4. Analysis of 3D data Pierre Ferruit; 5. Science motivation for IFS and galactic studies F. Eisenhauer; 6. Extragalactic studies and future IFS science Luis Colina; 7. Tutorials: how to handle 3D spectroscopy data Sebastian F. Sánchez, Begona García-Lorenzo and Arlette Pécontal-Rousset.
3D Elevation Program—Virtual USA in 3D
Lukas, Vicki; Stoker, J.M.
2016-01-01
The U.S. Geological Survey (USGS) 3D Elevation Program (3DEP) uses a laser system called ‘lidar’ (light detection and ranging) to create a virtual reality map of the Nation that is very accurate. 3D maps have many uses with new uses being discovered all the time.
NASA Astrophysics Data System (ADS)
Baraka, S. M.; Ben-Jaffel, L. B.
2014-12-01
We use particle-in-cell PIC 3D Electromagnetic, relativistic global code to address large-scale problems in magnetosphere electrodynamics. Terrestrial bow shock is simulated as an example. 3D Magnetohydrodynamics model ,MHD GUMICS in CCMC project, have been used in parallel with PIC under same scaled Solar wind (SW) and IMF conditions. We report new results from the coupling between the two models. Further investigations are required for confirmations of these results. In both codes the Earth's bow shock position is found at ~14.8 RE along the Sun-Earth line, and ~29 RE on the dusk side which is consistent with past in situ observation. Both simulations reproduce the theoretical jump conditions at the shock. However, PIC code density and temperature distributions are inflated and slightly shifted sunward when compared to MHD results. Reflected ions upstream of the bow shock may cause this sunward shift for density and temperature. Distribution of reflected ions and electrons are shown in the foreshock region, within the transition of the shock and in the downstream. The current version of PIC code can be run under modest computing facilities and resources. Additionally, existing MHD simulations should be useful to calibrate scaled properties of plasma resulting from PIC simulations for comparison with observations. Similarities and drawbacks of the results obtained by the two models are listed. The ultimate goal of using these different models in a complimentary manner rather than competitive is to better understand the macrostructure of the magnetosphere
MT3D was first developed by Chunmiao Zheng in 1990 at S.S. Papadopulos & Associates, Inc. with partial support from the U.S. Environmental Protection Agency (USEPA). Starting in 1990, MT3D was released as a pubic domain code from the USEPA. Commercial versions with enhanced capab...
NASA Technical Reports Server (NTRS)
1977-01-01
A market study of a proposed version of a 3-D eyetracker for initial use at NASA's Ames Research Center was made. The commercialization potential of a simplified, less expensive 3-D eyetracker was ascertained. Primary focus on present and potential users of eyetrackers, as well as present and potential manufacturers has provided an effective means of analyzing the prospects for commercialization.
Energy Science and Technology Software Center (ESTSC)
2013-10-01
Earth3D is a computer code designed to allow fast calculation of seismic rays and travel times through a 3D model of the Earth. LLNL is using this for earthquake location and global tomography efforts and such codes are of great interest to the Earth Science community.
[3-D ultrasound in gastroenterology].
Zoller, W G; Liess, H
1994-06-01
Three-dimensional (3D) sonography represents a development of noninvasive diagnostic imaging by real-time two-dimensional (2D) sonography. The use of transparent rotating scans, comparable to a block of glass, generates a 3D effect. The objective of the present study was to optimate 3D presentation of abdominal findings. Additional investigations were made with a new volumetric program to determine the volume of selected findings of the liver. The results were compared with the estimated volumes of 2D sonography and 2D computer tomography (CT). For the processing of 3D images, typical parameter constellations were found for the different findings, which facilitated processing of 3D images. In more than 75% of the cases examined we found an optimal 3D presentation of sonographic findings with respect to the evaluation criteria developed by us for the 3D imaging of processed data. Great differences were found for the estimated volumes of the findings of the liver concerning the three different techniques applied. 3D ultrasound represents a valuable method to judge morphological appearance in abdominal findings. The possibility of volumetric measurements enlarges its potential diagnostic significance. Further clinical investigations are necessary to find out if definite differentiation between benign and malign findings is possible. PMID:7919882
2013-10-30
This video provides an overview of the Sandia National Laboratories developed 3-D World Model Building capability that provides users with an immersive, texture rich 3-D model of their environment in minutes using a laptop and color and depth camera.
None
2014-02-26
This video provides an overview of the Sandia National Laboratories developed 3-D World Model Building capability that provides users with an immersive, texture rich 3-D model of their environment in minutes using a laptop and color and depth camera.
NASA Astrophysics Data System (ADS)
Grigoryan, G. V.; Grigoryan, R. P.
1995-03-01
The canonical quantization of a ( D=2n)-dimensional Dirac particle with spin in an arbitrary external electromagnetic field is performed in a gauge that makes it possible to describe simultaneously particles and antiparticles (both massive and massiess) already at the classical level. A pseudoclassical Foldy-Wouthuysen transformation is used to find the canonical (Newton-Wigner) coordinates. The connection between this quantization scheme and Blount's picture describing the behavior of a Dirac particle in an external electromagnetic field is discussed.
NASA Astrophysics Data System (ADS)
Walsh, J. R.
2004-02-01
The Euro3D RTN is an EU funded Research Training Network to foster the exploitation of 3D spectroscopy in Europe. 3D spectroscopy is a general term for spectroscopy of an area of the sky and derives its name from its two spatial + one spectral dimensions. There are an increasing number of instruments which use integral field devices to achieve spectroscopy of an area of the sky, either using lens arrays, optical fibres or image slicers, to pack spectra of multiple pixels on the sky (``spaxels'') onto a 2D detector. On account of the large volume of data and the special methods required to reduce and analyse 3D data, there are only a few centres of expertise and these are mostly involved with instrument developments. There is a perceived lack of expertise in 3D spectroscopy spread though the astronomical community and its use in the armoury of the observational astronomer is viewed as being highly specialised. For precisely this reason the Euro3D RTN was proposed to train young researchers in this area and develop user tools to widen the experience with this particular type of data in Europe. The Euro3D RTN is coordinated by Martin M. Roth (Astrophysikalisches Institut Potsdam) and has been running since July 2002. The first Euro3D science conference was held in Cambridge, UK from 22 to 23 May 2003. The main emphasis of the conference was, in keeping with the RTN, to expose the work of the young post-docs who are funded by the RTN. In addition the team members from the eleven European institutes involved in Euro3D also presented instrumental and observational developments. The conference was organized by Andy Bunker and held at the Institute of Astronomy. There were over thirty participants and 26 talks covered the whole range of application of 3D techniques. The science ranged from Galactic planetary nebulae and globular clusters to kinematics of nearby galaxies out to objects at high redshift. Several talks were devoted to reporting recent observations with newly
NASA Technical Reports Server (NTRS)
Walatka, Pamela P.; Buning, Pieter G.; Pierce, Larry; Elson, Patricia A.
1990-01-01
PLOT3D is a computer graphics program designed to visualize the grids and solutions of computational fluid dynamics. Seventy-four functions are available. Versions are available for many systems. PLOT3D can handle multiple grids with a million or more grid points, and can produce varieties of model renderings, such as wireframe or flat shaded. Output from PLOT3D can be used in animation programs. The first part of this manual is a tutorial that takes the reader, keystroke by keystroke, through a PLOT3D session. The second part of the manual contains reference chapters, including the helpfile, data file formats, advice on changing PLOT3D, and sample command files.
Dawood, A; Marti Marti, B; Sauret-Jackson, V; Darwood, A
2015-12-01
3D printing has been hailed as a disruptive technology which will change manufacturing. Used in aerospace, defence, art and design, 3D printing is becoming a subject of great interest in surgery. The technology has a particular resonance with dentistry, and with advances in 3D imaging and modelling technologies such as cone beam computed tomography and intraoral scanning, and with the relatively long history of the use of CAD CAM technologies in dentistry, it will become of increasing importance. Uses of 3D printing include the production of drill guides for dental implants, the production of physical models for prosthodontics, orthodontics and surgery, the manufacture of dental, craniomaxillofacial and orthopaedic implants, and the fabrication of copings and frameworks for implant and dental restorations. This paper reviews the types of 3D printing technologies available and their various applications in dentistry and in maxillofacial surgery. PMID:26657435
PLOT3D/AMES, APOLLO UNIX VERSION USING GMR3D (WITH TURB3D)
NASA Technical Reports Server (NTRS)
Buning, P.
1994-01-01
PLOT3D is an interactive graphics program designed to help scientists visualize computational fluid dynamics (CFD) grids and solutions. Today, supercomputers and CFD algorithms can provide scientists with simulations of such highly complex phenomena that obtaining an understanding of the simulations has become a major problem. Tools which help the scientist visualize the simulations can be of tremendous aid. PLOT3D/AMES offers more functions and features, and has been adapted for more types of computers than any other CFD graphics program. Version 3.6b+ is supported for five computers and graphic libraries. Using PLOT3D, CFD physicists can view their computational models from any angle, observing the physics of problems and the quality of solutions. As an aid in designing aircraft, for example, PLOT3D's interactive computer graphics can show vortices, temperature, reverse flow, pressure, and dozens of other characteristics of air flow during flight. As critical areas become obvious, they can easily be studied more closely using a finer grid. PLOT3D is part of a computational fluid dynamics software cycle. First, a program such as 3DGRAPE (ARC-12620) helps the scientist generate computational grids to model an object and its surrounding space. Once the grids have been designed and parameters such as the angle of attack, Mach number, and Reynolds number have been specified, a "flow-solver" program such as INS3D (ARC-11794 or COS-10019) solves the system of equations governing fluid flow, usually on a supercomputer. Grids sometimes have as many as two million points, and the "flow-solver" produces a solution file which contains density, x- y- and z-momentum, and stagnation energy for each grid point. With such a solution file and a grid file containing up to 50 grids as input, PLOT3D can calculate and graphically display any one of 74 functions, including shock waves, surface pressure, velocity vectors, and particle traces. PLOT3D's 74 functions are organized into
PLOT3D/AMES, APOLLO UNIX VERSION USING GMR3D (WITHOUT TURB3D)
NASA Technical Reports Server (NTRS)
Buning, P.
1994-01-01
PLOT3D is an interactive graphics program designed to help scientists visualize computational fluid dynamics (CFD) grids and solutions. Today, supercomputers and CFD algorithms can provide scientists with simulations of such highly complex phenomena that obtaining an understanding of the simulations has become a major problem. Tools which help the scientist visualize the simulations can be of tremendous aid. PLOT3D/AMES offers more functions and features, and has been adapted for more types of computers than any other CFD graphics program. Version 3.6b+ is supported for five computers and graphic libraries. Using PLOT3D, CFD physicists can view their computational models from any angle, observing the physics of problems and the quality of solutions. As an aid in designing aircraft, for example, PLOT3D's interactive computer graphics can show vortices, temperature, reverse flow, pressure, and dozens of other characteristics of air flow during flight. As critical areas become obvious, they can easily be studied more closely using a finer grid. PLOT3D is part of a computational fluid dynamics software cycle. First, a program such as 3DGRAPE (ARC-12620) helps the scientist generate computational grids to model an object and its surrounding space. Once the grids have been designed and parameters such as the angle of attack, Mach number, and Reynolds number have been specified, a "flow-solver" program such as INS3D (ARC-11794 or COS-10019) solves the system of equations governing fluid flow, usually on a supercomputer. Grids sometimes have as many as two million points, and the "flow-solver" produces a solution file which contains density, x- y- and z-momentum, and stagnation energy for each grid point. With such a solution file and a grid file containing up to 50 grids as input, PLOT3D can calculate and graphically display any one of 74 functions, including shock waves, surface pressure, velocity vectors, and particle traces. PLOT3D's 74 functions are organized into
Stanton, M M; Samitier, J; Sánchez, S
2015-08-01
Three-dimensional (3D) bioprinting has recently emerged as an extension of 3D material printing, by using biocompatible or cellular components to build structures in an additive, layer-by-layer methodology for encapsulation and culture of cells. These 3D systems allow for cell culture in a suspension for formation of highly organized tissue or controlled spatial orientation of cell environments. The in vitro 3D cellular environments simulate the complexity of an in vivo environment and natural extracellular matrices (ECM). This paper will focus on bioprinting utilizing hydrogels as 3D scaffolds. Hydrogels are advantageous for cell culture as they are highly permeable to cell culture media, nutrients, and waste products generated during metabolic cell processes. They have the ability to be fabricated in customized shapes with various material properties with dimensions at the micron scale. 3D hydrogels are a reliable method for biocompatible 3D printing and have applications in tissue engineering, drug screening, and organ on a chip models. PMID:26066320
Unassisted 3D camera calibration
NASA Astrophysics Data System (ADS)
Atanassov, Kalin; Ramachandra, Vikas; Nash, James; Goma, Sergio R.
2012-03-01
With the rapid growth of 3D technology, 3D image capture has become a critical part of the 3D feature set on mobile phones. 3D image quality is affected by the scene geometry as well as on-the-device processing. An automatic 3D system usually assumes known camera poses accomplished by factory calibration using a special chart. In real life settings, pose parameters estimated by factory calibration can be negatively impacted by movements of the lens barrel due to shaking, focusing, or camera drop. If any of these factors displaces the optical axes of either or both cameras, vertical disparity might exceed the maximum tolerable margin and the 3D user may experience eye strain or headaches. To make 3D capture more practical, one needs to consider unassisted (on arbitrary scenes) calibration. In this paper, we propose an algorithm that relies on detection and matching of keypoints between left and right images. Frames containing erroneous matches, along with frames with insufficiently rich keypoint constellations, are detected and discarded. Roll, pitch yaw , and scale differences between left and right frames are then estimated. The algorithm performance is evaluated in terms of the remaining vertical disparity as compared to the maximum tolerable vertical disparity.
Arena3D: visualization of biological networks in 3D
Pavlopoulos, Georgios A; O'Donoghue, Seán I; Satagopam, Venkata P; Soldatos, Theodoros G; Pafilis, Evangelos; Schneider, Reinhard
2008-01-01
Background Complexity is a key problem when visualizing biological networks; as the number of entities increases, most graphical views become incomprehensible. Our goal is to enable many thousands of entities to be visualized meaningfully and with high performance. Results We present a new visualization tool, Arena3D, which introduces a new concept of staggered layers in 3D space. Related data – such as proteins, chemicals, or pathways – can be grouped onto separate layers and arranged via layout algorithms, such as Fruchterman-Reingold, distance geometry, and a novel hierarchical layout. Data on a layer can be clustered via k-means, affinity propagation, Markov clustering, neighbor joining, tree clustering, or UPGMA ('unweighted pair-group method with arithmetic mean'). A simple input format defines the name and URL for each node, and defines connections or similarity scores between pairs of nodes. The use of Arena3D is illustrated with datasets related to Huntington's disease. Conclusion Arena3D is a user friendly visualization tool that is able to visualize biological or any other network in 3D space. It is free for academic use and runs on any platform. It can be downloaded or lunched directly from . Java3D library and Java 1.5 need to be pre-installed for the software to run. PMID:19040715
NASA Astrophysics Data System (ADS)
Otis, Collin; Ferrero, Pietro; Candler, Graham; Givi, Peyman
2013-11-01
The scalar filtered mass density function (SFMDF) methodology is implemented into the computer code US3D. This is an unstructured Eulerian finite volume hydrodynamic solver and has proven very effective for simulation of compressible turbulent flows. The resulting SFMDF-US3D code is employed for large eddy simulation (LES) on unstructured meshes. Simulations are conducted of subsonic and supersonic flows under non-reacting and reacting conditions. The consistency and the accuracy of the simulated results are assessed along with appraisal of the overall performance of the methodology. The SFMDF-US3D is now capable of simulating high speed flows in complex configurations.
Chilcoat, S.R. Hildebrand, S.T.
1995-12-31
Travel time computation in inhomogeneous media is essential for pre-stack Kirchhoff imaging in areas such as the sub-salt province in the Gulf of Mexico. The 2D algorithm published by Vinje, et al, has been extended to 3D to compute wavefronts in complicated inhomogeneous media. The 3D wavefront construction algorithm provides many advantages over conventional ray tracing and other methods of computing travel times in 3D. The algorithm dynamically maintains a reasonably consistent ray density without making a priori guesses at the number of rays to shoot. The determination of caustics in 3D is a straight forward geometric procedure. The wavefront algorithm also enables the computation of multi-valued travel time surfaces.
NASA Astrophysics Data System (ADS)
Yang, Xu; Zhang, Yong; Yang, Chenghua; Xu, Lu; Wang, Qiang; Zhao, Yuan
2016-06-01
Conventional three dimensional (3D) ghost imaging measures range of target based on pulse fight time measurement method. Due to the limit of data acquisition system sampling rate, range resolution of the conventional 3D ghost imaging is usually low. In order to take off the effect of sampling rate to range resolution of 3D ghost imaging, a heterodyne 3D ghost imaging (HGI) system is presented in this study. The source of HGI is a continuous wave laser instead of pulse laser. Temporal correlation and spatial correlation of light are both utilized to obtain the range image of target. Through theory analysis and numerical simulations, it is demonstrated that HGI can obtain high range resolution image with low sampling rate.
Combinatorial 3D Mechanical Metamaterials
NASA Astrophysics Data System (ADS)
Coulais, Corentin; Teomy, Eial; de Reus, Koen; Shokef, Yair; van Hecke, Martin
2015-03-01
We present a class of elastic structures which exhibit 3D-folding motion. Our structures consist of cubic lattices of anisotropic unit cells that can be tiled in a complex combinatorial fashion. We design and 3d-print this complex ordered mechanism, in which we combine elastic hinges and defects to tailor the mechanics of the material. Finally, we use this large design space to encode smart functionalities such as surface patterning and multistability.
Grigoryan, G.V.; Grigoryan, R.P.
1995-09-01
The canonical quantization of a (D=2n)-dimensional Dirac particle with spin in an arbitrary external electromagnetic field is performed in a gauge that makes it possible to describe simultaneously particles and antiparticles (both massive and massless) already at the classical level. A pseudoclassical Foldy-Wouthuysen transformation is used to find the canonical (Newton-Wigner) coordinates. The connection between this quantization scheme and Blount`s picture describing the behavior of a Dirac particle in an external electromagnetic field is discussed.
3-D magnetic field calculations for wiggglers using MAGNUS-3D
Pissanetzky, S.; Tompkins, P.
1988-01-01
The recent but steady trend toward increased magnetic and geometric complexity in the design of wigglers and undulators, of which tapered wigglers, hybrid structures, laced electromagnetic wigglers, magnetic cladding, twisters and magic structures are examples, has caused a need for reliable 3-D computer models and a better understanding of the behavior of magnetic systems in three dimensions. The capabilities of the MAGNUS-3D Group of Programs are ideally suited to solve this class of problems and provide insight into 3-D effects. MAGNUS-3D can solve any problem of Magnetostatics involving permanent magnets, linear or nonlinear ferromagnetic materials and electric conductors of any shape in space. The magnetic properties of permanent magnets are described by the complete nonlinear demagnetization curve as provided by the manufacturer, or, at the user's choice, by a simpler approximation involving the coercive force, the residual induction and the direction of magnetization. The ferromagnetic materials are described by a magnetization table and an accurate interpolation relation. An internal library with properties of common industrial steels is available. The conductors are independent of the mesh and are described in terms of conductor elements from an internal library.
Enhancements to the opera-3d suite
NASA Astrophysics Data System (ADS)
Riley, Christopher P.
1997-02-01
The OPERA-3D suite of programs has been enhanced to include 2 additional 3 dimensional finite element based solvers, with complimentary features in the pre- and postprocessing. SOPRANO computes electromagnetic fields at high frequency including displacement current effects. It has 2 modules—a deterministic solution at a user defined frequency and an eigenvalue solution for modal analysis. It is suitable for designing microwave structures and cavities found in particle accelerators. SCALA computes electrostatic fields in the presence of space charge from charged particle beams. The user may define the emission characteristics of electrodes or plasma surfaces and compute the resultant space charge limited beams, including the presence of magnetic fields. Typical applications in particle accelerators are electron guns and ion sources. Other enhancements to the suite include additional capabilities in TOSCA and ELEKTRA, the static and dynamic solvers.