## New Taylor-expansion method for solving a general class of wave equations

JOSA A, Vol. 15, Issue 8, pp. 2182-2198 (1998)

http://dx.doi.org/10.1364/JOSAA.15.002182

Enhanced HTML Acrobat PDF (1097 KB)

### Abstract

We derive a general method for solving second-order wave equations on the basis of a fourth-order Taylor expansion, with respect to time, of the field and its time derivative. The method requires explicit evaluation of space derivatives up to fourth order of the field and its time derivative, and the time step is subject to a stability criterion, similar to the one routinely applied in conjunction with the time-domain–finite-difference method. The Taylor-series method is generalizable to vector wave equations, to wave equations involving dissipation, and to the Helmholtz wave equation. The method is equally applicable to situations involving large or small variations in the refractive index. Sample numerical solutions in one and two space dimensions are described.

© 1998 Optical Society of America

**OCIS Codes**

(000.4430) General : Numerical approximation and analysis

(130.2790) Integrated optics : Guided waves

(230.3120) Optical devices : Integrated optics devices

(260.2110) Physical optics : Electromagnetic optics

(350.5500) Other areas of optics : Propagation

(350.7420) Other areas of optics : Waves

**History**

Original Manuscript: September 16, 1997

Revised Manuscript: January 7, 1998

Manuscript Accepted: February 6, 1998

Published: August 1, 1998

**Citation**

Joseph A. Fleck, "New Taylor-expansion method for solving a general class of wave equations," J. Opt. Soc. Am. A **15**, 2182-2198 (1998)

http://www.opticsinfobase.org/josaa/abstract.cfm?URI=josaa-15-8-2182

Sort: Year | Journal | Reset

### References

- K. S. Yee, “Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media,” IEEE Trans. Antennas Propag. AP-14, 302–307 (1966).
- A. Taflove, “Review of the formulation and applications of the finite-difference time-domain method for numerical modeling of electromagnetic wave interactions with arbitrary structures,” Wave Motion 10, 547–582 (1988). [CrossRef]
- W. C. Chew, Waves and Fields in Inhomogeneous Media (Van Nostrand Reinhold, New York, 1990).
- A. Taflove, K. R. Umashanker, “The finite-difference time-domain method for numerical modeling of electromagnetics wave interactions with arbitrary structures,” in Finite Element and Finite Difference Methods in Electromagnetic Scattering. (Elsevier, New York, 1990).
- S. T. Chu, W. P. Huang, S. K. Chaudari, “Simulation and analysis of waveguide based optical integrated circuits,” Comput. Phys. Pub. 68, 451–484 (1991). [CrossRef]
- W. P. Huang, S. T. Chu, A. Goss, S. K. Chaudhuri, “A scalar finite-difference time-domain approach for guided-wave optics,” IEEE Photonics Technol. Lett. 3, 524–526 (1991). [CrossRef]
- D. M. Young, R. T. Gregory, A Survey of Numerical Methods (Addison-Wesley, Reading, Mass., 1971), Vol. 1.
- M. D. Feit, J. A. Fleck, “Simple spectral method for solving propagation problems in cylindrical geometry with fast Fourier transforms,” Opt. Lett. 14, 662–664 (1989). [CrossRef] [PubMed]
- A. Splett, M. Majd, K. Petermann, “A novel beam propagation method for large refractive index steps and large propagation distances,” IEEE Photonics Technol. Lett. 3, 466–468 (1991). [CrossRef]
- R. P. Ratowsky, J. A. Fleck, “Treatment of angular derivatives in the Schrödinger equation by the finite Fourier series method,” J. Comput. Phys. 93, 376–387 (1991). [CrossRef]
- The split operator beam-propagation method takes advantage of the ease and accuracy of computing spatial derivatives by the FFT method. See, for example, J. A. Fleck, J. R. Morris, M. D. Feit, “Time-dependent propagation of high energy laser beams through the atmosphere,” Appl. Phys. 10, 129–160 (1976);J. A. Fleck, J. R. Morris, E. S. Bliss, “Small-scale self-focusing effects in a high power glass laser amplifier,” IEEE J. Quantum Electron. QE-14, 353–363 (1978); andM. D. Feit, J. A. Fleck, “Light propagation in graded index fibers,” Appl. Opt. 17, 3990–3998 (1978). [CrossRef] [PubMed]
- For a discussion of the sampling theorem and its implications on the accuracy of band-limited functions and their derivatives see, for example, E. O. Brigham, The Fast Fourier Transform (Prentice-Hall, Englewood Cliffs, N.J., 1974).
- M. D. Feit, J. A. Fleck, A. Steiger, “Solution of the Schrödinger equation by a spectral method,” J. Comput. Phys. 47, 412–433 (1982). [CrossRef]
- M. R. Hermann, J. A. Fleck, “Split operator spectral method for solving the time-dependent Schrödinger equation in spherical coordinates,” Phys. Rev. A 38, 6000–6012 (1988). [CrossRef] [PubMed]
- J. F. Claerbout, Fundamentals of Geophysical Data Processing with Applications to Petroleum Prospecting (McGraw-Hill, New York, 1976).
- D. Y. Lee, Y. Saad, M. H. Schultz, “An efficient method for solving the three-dimensional wide angle wave equation,” in Computational Acoustics, Vol 1—Wave Propagation, D. Lee, R. I. Sternberg, M. H. Schultz, eds. (North-Holland, Amsterdam, 1988).
- R. P. Ratowsky, J. A. Fleck, M. D. Feit, “Accurate solution of the Helmholtz equation by Lanczos orthoganalization for media with loss or gain,” Opt. Lett. 17, 10–12 (1992). [CrossRef] [PubMed]
- R. P. Ratowsky, J. A. Fleck, M. D. Feit, “Accurate description of ultra wide-angle beam propation in homogeneous media by Lanczos orthoganalization,” Opt. Lett. 19, 1284–1286 (1994). [CrossRef] [PubMed]
- J. A. Fleck “Solution of the scalar Helmholtz wave equation by Lanczos reduction,” in Methods for Modeling and Simulation of Guided-Wave Optoelectronic Devices: Part II: Waves and Interactions, W. P. Huang, ed. (EMW Publishing, Cambridge, Mass., 1995).
- A. Taflove, M. E. Baldwin, “Numerical solution of steady-state electromagnetic scattering problems using the time-dependent Maxwell’s equations,” IEEE Trans. Microwave Theory Tech. 23, 623–630 (1975). [CrossRef]
- See, for example, M. Born, E. Wolf, Principles of Optics (Pergamon, New York, 1959), pp. 22–23.
- See, for example, M. D. Feit, J. A. Fleck, “Filament formation and beam breakup in the self-focusing of optical beams,” J. Opt. Soc. Am. B 5, 633–644 (1988) and references therein. [CrossRef]

## Cited By |
Alert me when this paper is cited |

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

### Figures

Fig. 1 |
Fig. 2 |
Fig. 3 |

Fig. 4 |
Fig. 5 |
Fig. 6 |

Fig. 7 |
Fig. 8 |
Fig. 9 |

Fig. 10 |
Fig. 11 |
Fig. 12 |

Fig. 13 |
Fig. 14 |
Fig. 15 |

Fig. 16 |
Fig. 17 |
Fig. 18 |

Fig. 19 |
Fig. 20 |
Fig. 21 |

Fig. 22 |
||

« Previous Article | Next Article »

OSA is a member of CrossRef.