Rapid evaluation of two-dimensional retarded time integrals

Bizzozero, D. A.; Ellison, J. A.; Heinemann, K.; Lau, S. R.

doi:10.1016/j.cam.2017.04.007

Title: Rapid evaluation of two-dimensional retarded time integrals

Abstract

We present two methods for rapid evaluation of two-dimensional retarded time integrals. For example, such integrals arise as the z = 0 trace U(t, x, y, 0) of a solution U(t, x, y, z) to 3 + 1 wave equation U = –2f(t, x, y)δ(z) forced by a ‘‘sheet source’’ at z = 0. The spatial Fourier transform of a two-dimensional retarded time integral involves a temporal convolution with the zeroth order Bessel function J₀(t). Appealing to work by Alpert, Greengard, and Hagstrom and by Xu and Jiang on rational approximation in the Laplace-transform domain, our first method relies on approximation of J₀(t) as a sum of exponential functions. We achieve approximations with double precision accuracy near t ≃ 0, and maintain single precision accuracy out to T ≃ 10⁸. Our second method involves evolution of the 3 + 1 wave equation in a ‘‘thin block’’ above the sheet, adopting the radiation boundary conditions of Hagstrom, Warburton, and Givoli based on complete plane wave expansions. We review their technique, present its implementation for our problem, and present new results on the nonlocal spacetime form of radiation boundary conditions. Our methods are relevant for the sheet-bunch formulation of the Vlasov–Maxwell system,more »« less

Authors:

Bizzozero, D. A. ^[1]; Ellison, J. A. ^[1]; Heinemann, K. ^[1]; Lau, S. R. ^[1]

Univ. of New Mexico, Albuquerque, NM (United States)

Publication Date:: 2017-04-12

Research Org.:: Univ. of New Mexico, Albuquerque, NM (United States); Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)

Sponsoring Org.:: USDOE Office of Science (SC), High Energy Physics (HEP) (SC-25)

OSTI Identifier:: 1511845

Alternate Identifier(s):: OSTI ID: 1416621; OSTI ID: 1530304

Grant/Contract Number:: FG02-99ER41104; FG03-99ER41104; FG-99ER41104; AC02-05CH11231

Resource Type:: Accepted Manuscript

Journal Name:: Journal of Computational and Applied Mathematics

Additional Journal Information:: Journal Volume: 324; Journal Issue: C; Journal ID: ISSN 0377-0427

Publisher:: Elsevier

Country of Publication:: United States

Language:: English

Subject:: 98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL PROTECTION; Retarded time integral; Rational approximation; Radiation boundary conditions; Initial boundary value problem; Vlasov–Maxwell system; Accelerator beam physics

Citation Formats


                    Bizzozero, D. A., Ellison, J. A., Heinemann, K., and Lau, S. R. Rapid evaluation of two-dimensional retarded time integrals.  United States: N. p., 2017. 
        Web.  doi:10.1016/j.cam.2017.04.007.

Copy to clipboard


                    Bizzozero, D. A., Ellison, J. A., Heinemann, K., & Lau, S. R. Rapid evaluation of two-dimensional retarded time integrals.  United States.  doi:https://doi.org/10.1016/j.cam.2017.04.007

Copy to clipboard


                    Bizzozero, D. A., Ellison, J. A., Heinemann, K., and Lau, S. R. Wed .  
        "Rapid evaluation of two-dimensional retarded time integrals".  United States.  doi:https://doi.org/10.1016/j.cam.2017.04.007.  https://www.osti.gov/servlets/purl/1511845.

Copy to clipboard


                    
@article{osti_1511845,

  title        = {Rapid evaluation of two-dimensional retarded time integrals},

  author       = {Bizzozero, D. A. and Ellison, J. A. and Heinemann, K. and Lau, S. R.},

  abstractNote = {We present two methods for rapid evaluation of two-dimensional retarded time integrals. For example, such integrals arise as the z = 0 trace U(t, x, y, 0) of a solution U(t, x, y, z) to 3 + 1 wave equation U = –2f(t, x, y)δ(z) forced by a ‘‘sheet source’’ at z = 0. The spatial Fourier transform of a two-dimensional retarded time integral involves a temporal convolution with the zeroth order Bessel function J0(t). Appealing to work by Alpert, Greengard, and Hagstrom and by Xu and Jiang on rational approximation in the Laplace-transform domain, our first method relies on approximation of J0(t) as a sum of exponential functions. We achieve approximations with double precision accuracy near t ≃ 0, and maintain single precision accuracy out to T ≃ 108. Our second method involves evolution of the 3 + 1 wave equation in a ‘‘thin block’’ above the sheet, adopting the radiation boundary conditions of Hagstrom, Warburton, and Givoli based on complete plane wave expansions. We review their technique, present its implementation for our problem, and present new results on the nonlocal spacetime form of radiation boundary conditions. Our methods are relevant for the sheet-bunch formulation of the Vlasov–Maxwell system, although here we only test methods on a model problem, a Gaussian source following an elliptical orbit. Here, our concluding section discusses the complexity of both methods in comparison with naive evaluation of a retarded-time integral.},

  doi          = {10.1016/j.cam.2017.04.007},

  journal      = {Journal of Computational and Applied Mathematics},

  number       = C,

  volume       = 324,

  place        = {United States},

  year         = {2017},

  month        = {4}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (Publisher)

Accepted Manuscript (DOE)

Publisher's Version of Record

DOI: https://doi.org/10.1016/j.cam.2017.04.007

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 1 work

Citation information provided by
Web of Science

Figures / Tables:

Figure 1: Retarded-time integration in terms of bunch history. In the right-hand figure the $z$-coordinate normal to the sheet is suppressed.

All figures and tables (16 total)

Save / Share:

Export Metadata

Save to My Library

Figures / Tables found in this record:

Figures/Tables have been extracted from DOE-funded journal article accepted manuscripts.

Similar Records in DOE PAGES and OSTI.GOV collections:

Stable, high-order discretization for evolution of the wave equation in 1+1 dimensions

Journal Article Visher, John ; Wandzura, S ; White, Amanda M - Journal of Computational Physics, 194(2):395-408

We carry forward the approach of Alpert, Greengard, and Hagstrom to construct stable high-order explicit discretizations for the wave equation in one space and one time dimension. They presented their scheme as an integral form of the Lax-Wendroff method. Our perspective is somewhat different from theirs; our focus is on the discretization of the evolution formula rather than on its form (integral, differential, etc.). A key feature of our approach is the independent computation of three discretizations, one for bulk (away from boundaries) propagation, one for propagation near boundaries, and a projection operator to enforce boundary conditions. This is donemore »« less
DOI: https://doi.org/10.1016/j.jcp.2003.09.028
Analytic solution to leading order coupled DGLAP evolution equations: A new perturbative QCD tool

Journal Article Block, Martin M ; Durand, Loyal ; Ha, Phuoc ; ... - Physical Review. D, Particles Fields

We have analytically solved the LO perturbative QCD singlet DGLAP equations [V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15, 438 (1972)][G. Altarelli and G. Parisi, Nucl. Phys. B126, 298 (1977)][Y. L. Dokshitzer, Sov. Phys. JETP 46, 641 (1977)] using Laplace transform techniques. Newly developed, highly accurate, numerical inverse Laplace transform algorithms [M. M. Block, Eur. Phys. J. C 65, 1 (2010)][M. M. Block, Eur. Phys. J. C 68, 683 (2010)] allow us to write fully decoupled solutions for the singlet structure function F{sub s}(x,Q{sup 2}) and G(x,Q{sup 2}) as F{sub s}(x,Q{sup 2})=F{sub s}(F{sub s0}(x{sub 0}),G{sub 0}(x{submore »« less
DOI: https://doi.org/10.1103/PHYSREVD.83.054009
Measuring the Higgs-vector boson couplings at linear e{sup +}e{sup -} collider

Journal Article Dutta, Sukanta ; KEK Theory Division, 1-1 Oho Tsukuba, 305-0801 ; Hagiwara, Kaoru ; ... - Physical Review. D, Particles Fields

We estimate the accuracy with which the coefficient of the CP even dimension-six operators involving Higgs and two vector bosons (HVV) can be measured at linear e{sup +}e{sup -} colliders. Using the optimal observables method for the kinematic distributions, our analysis is based on the five different processes. First is the WW fusion process in the t-channel (e{sup +}e{sup -}{yields}{nu}{sub e}{nu}{sub e}H), where we use the rapidity y and the transverse momentum p{sub T} of the Higgs boson as observables. Second is the ZH pair production process in the s channel, where we use the scattering angle of the Zmore »« less
DOI: https://doi.org/10.1103/PHYSREVD.78.115016
Superconducting instability in the infinite-[ital U] Anderson lattice in the presence of crystal electric fields

Journal Article Trees, B R - Physical Review, B: Condensed Matter; (United States)

We have studied the mean-field quasiparticle states and quasiparticle interactions in the infinite-[ital U] Anderson lattice in the presence of crystal electric field of cubic symmetry. We assume a lattice of 4[ital f] sites, each with a total angular momentum of [ital J]=5/2 that is split into a low-lying doublet of [Gamma][sub 7] symmetry and an excited quartet of [Gamma][sub 8] symmetry. Slave bosons on the 4[ital f] sites create and destroy 4[ital f][sup 0] configurations and Lagrange multipliers at each 4[ital f] site enforce the occupancy constraint due to the infinite Coulomb repulsion. The quasiparticle states are strongly peakedmore »« less
DOI: https://doi.org/10.1103/PhysRevB.51.470
The minimal set of electroweak precision parameters

Journal Article Cacciapaglia, G ; Csaki, C ; Marandella, G ; ... - Physical Review. D, Particles Fields

We present a simple method for analyzing the impact of precision electroweak data above and below the Z peak on flavor-conserving, heavy new physics. We find that experiments have probed about ten combinations of new-physics effects, which to a good approximation can be condensed into the effective oblique parameters S-circumflex, T-circumflex, U-circumflex, V, X, W, Y (we prove positivity constraints W,Y{>=}0) and three combinations of quark couplings (including a distinct parameter for the bottom). We apply our method to generic extra Z{sup '} vectors.
DOI: https://doi.org/10.1103/PHYSREVD.74.033011

Similar Records