<abstract xmlns="http://www.w3.org/1999/xhtml">

<p>This paper presents a numerical method for solving the two-dimensional Laplace equation in curvilinear coordinates using a finite difference scheme. Laplace equation is a differential equation used to describe various physical phenomena studied within different branches of engineering. The two-dimensional Laplace equation has analytic solutions only for simple shaped domains and simple boundary conditions (first and second kind). Even so, the results are usually complicated to interpret and use. The present paper proposes a method to solve the equation for domains bounded by complex inner and outer curves. The boundary curves are defined by a set of points from which the equations of curves are approximated by using Fourier series. The domain is generated by blending the two curves using a blending function. The results are a set of continuous parametric equations which are defining the domain in angle and radius variables. Following a method presented in the literature, the Laplace equation is transformed from Cartesian coordinates to curvilinear coordinates, as the Fourier series is easily differentiable. Next, the new Laplace equation is approximated using the finite difference method and a MATLAB code is developed to solve the system of linear equations using successive over-relaxation method (SOR). The presented method is used to simulate steady two-dimensional heat conduction in a polygonal shape domain, with uniform and non-uniform Dirichlet boundary conditions. Results are analyzed in order to determine the model’s precision, speed and convergence criteria.</p>
</abstract>

This paper presents a numerical method for solving the two-dimensional Laplace equation in curvilinear coordinates using a finite difference scheme. Laplace equation is a differential equation used to describe various physical phenomena studied within different branches of engineering. The two-dimensional Laplace equation has analytic solutions only for simple shaped domains and simple boundary conditions (first and second kind). Even so, the results are usually complicated to interpret and use. The present paper proposes a method to solve the equation for domains bounded by complex inner and outer curves. The boundary curves are defined by a set of points from which the equations of curves are approximated by using Fourier series. The domain is generated by blending the two curves using a blending function. The results are a set of continuous parametric equations which are defining the domain in angle and radius variables. Following a method presented in the literature, the Laplace equation is transformed from Cartesian coordinates to curvilinear coordinates, as the Fourier series is easily differentiable. Next, the new Laplace equation is approximated using the finite difference method and a MATLAB code is developed to solve the system of linear equations using successive over-relaxation method (SOR). The presented method is used to simulate steady two-dimensional heat conduction in a polygonal shape domain, with uniform and non-uniform Dirichlet boundary conditions. Results are analyzed in order to determine the model’s precision, speed and convergence criteria.

Finite Difference Method for Solving the Two-Dimensional Laplace Equation in Curvilinear Coordinates

gabriel-dumitru.tcaciuc@academic.tuiasi.ro

“Gheorghe Asachi” Technical University of Iaşi

Bulletin of the Polytechnic Institute of Iași. Machine constructions Section

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

This paper presents a numerical method for solving the two-dimensional Laplace equation in curvilinear coordinates using a finite difference scheme. Laplace...

Finite Difference Method for Solving the Two-Dimensional Laplace Equation in Curvilinear Coordinates

Published Online: Jan 22, 2024

Page range: 109 - 121

Received: Dec 15, 2023

Accepted: Dec 30, 2023

DOI: https://doi.org/10.2478/bipcm-2023-0028

Keywords
finite difference method, curvilinear coordinates, heat conduction, Fourier series, Laplace equation

© 2023 Gabriel-Dumitru Tcaciuc et al., published by Sciendo

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

Finite Difference Method for Solving the Two-Dimensional Laplace Equation in Curvilinear Coordinates

Published Online: Jan 22, 2024

Page range: 109 - 121

Received: Dec 15, 2023

Accepted: Dec 30, 2023

DOI: https://doi.org/10.2478/bipcm-2023-0028

Keywordsfinite difference method, curvilinear coordinates, heat conduction, Fourier series, Laplace equation

© 2023 Gabriel-Dumitru Tcaciuc et al., published by Sciendo

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

Keywords
finite difference method, curvilinear coordinates, heat conduction, Fourier series, Laplace equation