Systems of grid equations that approximate elliptic boundary value problems on locally modified grids are considered. The triangulation, which approximates the boundary with second order of accuracy, is generated from an initial uniform triangulation by shifting nodes near the boundary according to special rules. This "locally modified" grid possesses several significant features: this triangulation has a regular structure, the generation of the triangulation is rather fast, this construction allows to use multilevel preconditioning (BPX-like) methods. The proposed iterative methods for solving elliptic boundary value problems approximately are based on two approaches: The fictitious space method, i.e. the reduction of the original problem to a problem in an auxiliary (fictitious) space, and the multilevel decomposition method, i.e. the construction of preconditioners by decomposing functions on hierarchical grids. The convergence rate of the corresponding iterative process with the preconditioner obtained is independent of the mesh size. The construction of the grid and the preconditioning operator for the three dimensional problem can be done in the same way.

京都大学 / 0048 / 新制・課程博士 / 博士(情報学) / 甲第21212号 / 情博第665号 / 新制||情||115(附属図書館) / 京都大学大学院情報学研究科先端数理科学専攻 / (主査)教授 磯 祐介, 教授 木上 淳, 助手 藤原 宏志 / 学位規則第4条第1項該当 / Doctor of Informatics / Kyoto University / DFAM

No description available.

The inadequacy of Fick’s law to incorporate causality can be overcome by replacing it with the Green–Naghdi type II (GNII) flux relation. Combining the GNII assumption and conservation of mass leads to [see document for equation] where r (x, t) is the density function, S(p) is a source term and c¥ is a positive constant which carries (SI) units of m/sec. A general source term given by [see document for equation] is proposed. Here, the constants y and ps are the rate coefficient and saturation density respectively. The travelling wave solutions and numerical analysis of four special cases of equation (2), namely: Pearl-Verhulst Growth law, Zel’dovich Law, Newmann Law and Stefan- Boltzmann Law are investigated. For both analysis, results are compared with the available literature and extended for other cases. The numerical analysis is carried out by imposing well-studied Initial Boundary Value Problem and implementing a built-in method in the software package Mathematica 9. For Pearl-Verhulst source type, the results are compared to those found in literature [1]. Confirming the validity of built-in method for Pearl-Verhulst law, the generic built-in method is extended to study the transient signal response for similar initial boundary value problems when the source terms are Zel’dovich law, Newmann law and Stefan-Boltzmann law.

The analytical solution of a given inhomogeneous boundary value problem of a linear, elliptic, inhomogeneous partial differential equation and a set of inhomogeneous mixed Dirichlet- and Neumann-type boundary conditions is derived in the present paper. In the context of elasticity theory, the problem arises for a non-conservative symmetric ansatz and an extended constitutive law shown earlier. For convenient user application, the scalar function expressed in cylindrical coordinates is primarily obtained for the general case before being expatiated on a special case of linear boundary conditions.

In addition to previous publications, the paper presents the analytical solution of a special boundary value problem which arises in the context of elasticity theory for an extended constitutive law and a non-conservative symmetric ansatz. Besides deriving the general analytical solution, a specific form for linear boundary conditions is given for user convenience.

No description available.

This thesis consists of three self-contained parts, each with its own abstract, body, references and page numbering. Part I, 'Potential theory, path integrals and the Laplacian of the indicator', finds the transition density of absorbed or reflected Brownian motion in a d-dimensional domain as a Feynman-Kac functional involving the Laplacian of the indicator, thereby relating the hitherto unrelated fields of classical potential theory and path integrals. Part II, 'The problem of alternatives', considers parallel investment in alternative technologies or drugs developed over time, where there can be only one winner. Parallel investment accelerates the search for the winner, and increases the winner's expected performance, but is also costly. To determine which candidates show sufficient performance and/or promise, we find an integral equation for the boundary of the optimal continuation region. Part III, 'Optimal support for renewable deployment', considers the role of government subsidies for renewable technologies. Rapidly diminishing subsidies are cheaper for taxpayers, but could prematurely kill otherwise successful technologies. By contrast, high subsidies are not only expensive but can also prop up uneconomical technologies. To analyse this trade-off we present a new model for technology learning that makes capacity expansion endogenous. There are two reasons for this standalone structure. First, the target readership is divergent. Part I concerns mathematical physics, Part II operations research, and Part III policy. Readers interested in specific parts can thus read these in isolation. Those interested in the thesis as a whole may prefer to read the three introductions first. Second, the separate parts are only partially interconnected. Each uses some theory from the preceding part, but not all of it; e.g. Part II uses only a subset of the theory from Part I. The quickest route to Part III is therefore not through the entirety of the preceding parts. Furthermore, those instances where results from previous parts are used are clearly indicated.

No description available.

One of the main problems of current physical geodesy is the unification of local height systems and the creation of the unified global vertical reference frame, whose primary definition constant will be the geopotential value of W0 on the level surface, such as geoid. This problem encounters several pitfalls, such as different types of physical heights used in the world, ways of eliminating the effects of tides on the shape of the Earth's body. The first parts of the thesis describe the theoretical foundations concerning the description of the Earth's gravitational field, the basics of height theory and the solution of the boundary value problems used in geodesy, together with the solution of the Hotine's and Stokes' integral by convolution. Due to the fact that the data of directly measured gravity acceleration for the territory of the Czechia and the Slovakia are not freely available with sufficient coverage, these data were replaced by gravity disturbances calculated from the Global Gravity Model of the Earth. To improve the data obtained from geopotential models, so-called residual terrain modelling was used. It is a spectral combination of gravity field models with relevant Earth's gravitational field quantities, which are obtained by modeling from a digital terrain model and a height model that represents the mean value of topography. The combination of these data consists of calculating those frequencies of the gravitational signal from the digital terrain model that are not a part of the signal obtained from geopotential models that in general have a lower spatial resolution. Two methods were used to connect the local height systems of the Czechia and the Slovakia. The first method was to use the solution of Molodensky's problem. The second method used is originally a method developed to test geopotential models. Both methods were tested on a set of measuring points (GNSS/levelling points).