19,673 research outputs found
Numerical solving unsteady space-fractional problems with the square root of an elliptic operator
An unsteady problem is considered for a space-fractional equation in a
bounded domain. A first-order evolutionary equation involves the square root of
an elliptic operator of second order. Finite element approximation in space is
employed. To construct approximation in time, regularized two-level schemes are
used. The numerical implementation is based on solving the equation with the
square root of the elliptic operator using an auxiliary Cauchy problem for a
pseudo-parabolic equation. The scheme of the second-order accuracy in time is
based on a regularization of the three-level explicit Adams scheme. More
general problems for the equation with convective terms are considered, too.
The results of numerical experiments are presented for a model two-dimensional
problem.Comment: 21 pages, 7 figures. arXiv admin note: substantial text overlap with
arXiv:1412.570
A Method for Geometry Optimization in a Simple Model of Two-Dimensional Heat Transfer
This investigation is motivated by the problem of optimal design of cooling
elements in modern battery systems. We consider a simple model of
two-dimensional steady-state heat conduction described by elliptic partial
differential equations and involving a one-dimensional cooling element
represented by a contour on which interface boundary conditions are specified.
The problem consists in finding an optimal shape of the cooling element which
will ensure that the solution in a given region is close (in the least squares
sense) to some prescribed target distribution. We formulate this problem as
PDE-constrained optimization and the locally optimal contour shapes are found
using a gradient-based descent algorithm in which the Sobolev shape gradients
are obtained using methods of the shape-differential calculus. The main novelty
of this work is an accurate and efficient approach to the evaluation of the
shape gradients based on a boundary-integral formulation which exploits certain
analytical properties of the solution and does not require grids adapted to the
contour. This approach is thoroughly validated and optimization results
obtained in different test problems exhibit nontrivial shapes of the computed
optimal contours.Comment: Accepted for publication in "SIAM Journal on Scientific Computing"
(31 pages, 9 figures
q-deformed KZB heat equation: completeness, modular properties and SL(3,Z)
We study the properties of one-dimensional hypergeometric integral solutions
of the q-difference ("quantum") analogue of the Knizhnik-Zamolodchikov-Bernard
equations on tori. We show that they also obey a difference KZB heat equation
in the modular parameter, give formulae for modular transformations, and prove
a completeness result, by showing that the associated Fourier transform is
invertible. These results are based on SL(3,Z) transformation properties
parallel to those of elliptic gamma functions.Comment: 39 page
Immersed Boundary Smooth Extension: A high-order method for solving PDE on arbitrary smooth domains using Fourier spectral methods
The Immersed Boundary method is a simple, efficient, and robust numerical
scheme for solving PDE in general domains, yet it only achieves first-order
spatial accuracy near embedded boundaries. In this paper, we introduce a new
high-order numerical method which we call the Immersed Boundary Smooth
Extension (IBSE) method. The IBSE method achieves high-order accuracy by
smoothly extending the unknown solution of the PDE from a given smooth domain
to a larger computational domain, enabling the use of simple Cartesian-grid
discretizations (e.g. Fourier spectral methods). The method preserves much of
the flexibility and robustness of the original IB method. In particular, it
requires minimal geometric information to describe the boundary and relies only
on convolution with regularized delta-functions to communicate information
between the computational grid and the boundary. We present a fast algorithm
for solving elliptic equations, which forms the basis for simple, high-order
implicit-time methods for parabolic PDE and implicit-explicit methods for
related nonlinear PDE. We apply the IBSE method to solve the Poisson, heat,
Burgers', and Fitzhugh-Nagumo equations, and demonstrate fourth-order pointwise
convergence for Dirichlet problems and third-order pointwise convergence for
Neumann problems
- …