skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A New Method for Solving the Difference Equations with Variable Coefficients

Abstract

Very-effective method of solving of difference equations is the operator's one, where the operators are the translational ones. It is demonstrated how this method gives the solution of the second order equation with variable function. The solution is applicable to investigation of Frenkel's excitons behavior in thin molecular film.

Authors:
 [1];  [2];  [3]
  1. Technical Faculty 'M. Pupin', Zrenjanin (Serbia)
  2. Vojvodina Academy of Sciences and Arts, Novi Sad (Serbia)
  3. Department of Physics, Faculty of Sciences, Novi Sad (Serbia)
Publication Date:
OSTI Identifier:
21057259
Resource Type:
Journal Article
Resource Relation:
Journal Name: AIP Conference Proceedings; Journal Volume: 899; Journal Issue: 1; Conference: 6. international conference of the Balkan Physical Union, Istanbul (Turkey), 22-26 Aug 2006; Other Information: DOI: 10.1063/1.2733450; (c) 2007 American Institute of Physics; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; EQUATIONS; EXCITONS; FRENKEL DEFECTS; MATHEMATICAL OPERATORS; MATHEMATICAL SOLUTIONS; MOLECULES; THIN FILMS

Citation Formats

Sajfert, V. D., Tosic, B. S., and Setrajcic, J. P.. A New Method for Solving the Difference Equations with Variable Coefficients. United States: N. p., 2007. Web. doi:10.1063/1.2733450.
Sajfert, V. D., Tosic, B. S., & Setrajcic, J. P.. A New Method for Solving the Difference Equations with Variable Coefficients. United States. doi:10.1063/1.2733450.
Sajfert, V. D., Tosic, B. S., and Setrajcic, J. P.. Mon . "A New Method for Solving the Difference Equations with Variable Coefficients". United States. doi:10.1063/1.2733450.
@article{osti_21057259,
title = {A New Method for Solving the Difference Equations with Variable Coefficients},
author = {Sajfert, V. D. and Tosic, B. S. and Setrajcic, J. P.},
abstractNote = {Very-effective method of solving of difference equations is the operator's one, where the operators are the translational ones. It is demonstrated how this method gives the solution of the second order equation with variable function. The solution is applicable to investigation of Frenkel's excitons behavior in thin molecular film.},
doi = {10.1063/1.2733450},
journal = {AIP Conference Proceedings},
number = 1,
volume = 899,
place = {United States},
year = {Mon Apr 23 00:00:00 EDT 2007},
month = {Mon Apr 23 00:00:00 EDT 2007}
}
  • We present a method for solving Poisson and heat equations with discon- tinuous coefficients in two- and three-dimensions. It uses a Cartesian cut-cell/embedded boundary method to represent the interface between materi- als, as described in Johansen& Colella (1998). Matching conditions across the interface are enforced using an approximation to fluxes at the boundary. Overall second order accuracy is achieved, as indicated by an array of tests using non-trivial interface geometries. Both the elliptic and heat solvers are shown to remain stable and efficient for material coefficient contrasts up to 106, thanks in part to the use of geometric multigrid. Amore » test of accuracy when adaptive mesh refinement capabilities are utilized is also performed. An example problem relevant to nuclear reactor core simulation is presented, demonstrating the ability of the method to solve problems with realistic physical parameters.« less
  • We present a method for solving Poisson and heat equations with discontinuous coefficients in two- and three-dimensions. It uses a Cartesian cut-cell/embedded boundary method to represent the interface between materials, as described in Johansen and Colella (1998). Matching conditions across the interface are enforced using an approximation to fluxes at the boundary. Overall second order accuracy is achieved, as indicated by an array of tests using non-trivial interface geometries. Both the elliptic and heat solvers are shown to remain stable and efficient for material coefficient contrasts up to 10{sup 6}, thanks in part to the use of geometric multigrid. Amore » test of accuracy when adaptive mesh refinement capabilities are utilized is also performed. An example problem relevant to nuclear reactor core simulation is presented, demonstrating the ability of the method to solve problems with realistic physical parameters.« less
  • We present an efficient method to solve numerically the equations of dissipative dynamics of the binary phase-field crystal model proposed by Elder et al. [K.R. Elder, M. Katakowski, M. Haataja, M. Grant, Phys. Rev. B 75 (2007) 064107] characterized by variable coefficients. Using the operator splitting method, the problem has been decomposed into sub-problems that can be solved more efficiently. A combination of non-trivial splitting with spectral semi-implicit solution leads to sets of algebraic equations of diagonal matrix form. Extensive testing of the method has been carried out to find the optimum balance among errors associated with time integration, spatialmore » discretization, and splitting. We show that our method speeds up the computations by orders of magnitude relative to the conventional explicit finite difference scheme, while the costs of the pointwise implicit solution per timestep remains low. Also we show that due to its numerical dissipation, finite differencing can not compete with spectral differencing in terms of accuracy. In addition, we demonstrate that our method can efficiently be parallelized for distributed memory systems, where an excellent scalability with the number of CPUs is observed.« less