9,831 research outputs found
Inverse heat conduction problems by using particular solutions
Based on the method of fundamental solutions, we develop in this paper a new computational method to solve two-dimensional transient heat conduction inverse problems. The main idea is to use particular solutions as radial basis functions (PSRBF) for approximation of the solutions to the inverse heat conduction problems. The heat conduction equations are first analyzed in the Laplace transformed domain and the Durbin inversion method is then used to determine the solutions in the time domain. Least-square and singular value decomposition (SVD) techniques are adopted to solve the ill-conditioned linear system of algebraic equations obtained from the proposed PSRBF method. To demonstrate the effectiveness and simplicity of this approach, several numerical examples are given with satisfactory accuracy and stability.Peer reviewe
Transparent boundary conditions based on the Pole Condition for time-dependent, two-dimensional problems
The pole condition approach for deriving transparent boundary conditions is
extended to the time-dependent, two-dimensional case. Non-physical modes of the
solution are identified by the position of poles of the solution's spatial
Laplace transform in the complex plane. By requiring the Laplace transform to
be analytic on some problem dependent complex half-plane, these modes can be
suppressed. The resulting algorithm computes a finite number of coefficients of
a series expansion of the Laplace transform, thereby providing an approximation
to the exact boundary condition. The resulting error decays super-algebraically
with the number of coefficients, so relatively few additional degrees of
freedom are sufficient to reduce the error to the level of the discretization
error in the interior of the computational domain. The approach shows good
results for the Schr\"odinger and the drift-diffusion equation but, in contrast
to the one-dimensional case, exhibits instabilities for the wave and
Klein-Gordon equation. Numerical examples are shown that demonstrate the good
performance in the former and the instabilities in the latter case
Boundary knot method: A meshless, exponential convergence, integration-free, and boundary-only RBF technique
Based on the radial basis function (RBF), non-singular general solution and
dual reciprocity principle (DRM), this paper presents an inheretnly meshless,
exponential convergence, integration-free, boundary-only collocation techniques
for numerical solution of general partial differential equation systems. The
basic ideas behind this methodology are very mathematically simple and
generally effective. The RBFs are used in this study to approximate the
inhomogeneous terms of system equations in terms of the DRM, while non-singular
general solution leads to a boundary-only RBF formulation. The present method
is named as the boundary knot method (BKM) to differentiate it from the other
numerical techniques. In particular, due to the use of non-singular general
solutions rather than singular fundamental solutions, the BKM is different from
the method of fundamental solution in that the former does no need to introduce
the artificial boundary and results in the symmetric system equations under
certain conditions. It is also found that the BKM can solve nonlinear partial
differential equations one-step without iteration if only boundary knots are
used. The efficiency and utility of this new technique are validated through
some typical numerical examples. Some promising developments of the BKM are
also discussed.Comment: 36 pages, 2 figures, Welcome to contact me on this paper: Email:
[email protected] or [email protected]
The energetics of melting fertile heterogeneities within the depleted mantle
To explore the consequences of mantle heterogeneity for primary melt production, we develop a
mathematical model of energy conservation for an upwelling, melting body of recycled oceanic crust
embedded in the depleted upper mantle. We consider the endâmember geometric cases of spherical blobs
and tabular veins. The model predicts that thermal diffusion into the heterogeneity can cause a factorâofâ
two increase in the degree of melting for bodies with minimum dimension smaller than âŒ1 km, yielding
melt fractions between 50 and 80%. The role of diffusion is quantified by an appropriately defined Peclet
number, which represents the balance of diffusionâdriven and adiabatic melting. At intermediate Peclet
number, we show that melting a heterogeneity can cool the ambient mantle by up to âŒ20 K (spherical)
or âŒ60 K (tabular) within a distance of two times the characteristic size of the body. At small Peclet
number, where heterogeneities are expected to be in thermal equilibrium with the ambient mantle, we
calculate the energetic effect of pyroxenite melting on the surrounding peridotite; we find that each 5%
of recycled oceanic crust diminishes the peridotite degree of melting by 1â2%. Injection of the magma from
highly molten bodies of recycled oceanic crust into a melting region of depleted upper mantle may nucleate
reactiveâdissolution channels that remain chemically isolated from the surrounding peridotite
Efficient Large Scale Transient Heat Conduction Analysis Using A Parallelized Boundary Element Method
A parallel domain decomposition Laplace transform Boundary Element Method, BEM, algorithm for the solution of large-scale transient heat conduction problems will be developed. This is accomplished by building on previous work by the author and including several new additions (most note-worthy is the extension to 3-D) aimed at extending the scope and improving the efficiency of this technique for large-scale problems. A Laplace transform method is utilized to avoid time marching and a Proper Orthogonal Decomposition, POD, interpolation scheme is used to improve the efficiency of the numerical Laplace inversion process. A detailed analysis of the Stehfest Transform (numerical Laplace inversion) is performed to help optimize the procedure for heat transfer problems. A domain decomposition process is described in detail and is used to significantly reduce the size of any single problem for the BEM, which greatly reduces the storage and computational burden of the BEM. The procedure is readily implemented in parallel and renders the BEM applicable to large-scale transient conduction problems on even modest computational platforms. A major benefit of the Laplace space approach described herein, is that it readily allows adaptation and integration of traditional BEM codes, as the resulting governing equations are time independent. This work includes the adaptation of two such traditional BEM codes for steady-state heat conduction, in both two and three dimensions. Verification and validation example problems are presented which show the accuracy and efficiency of the techniques. Additionally, comparisons to commercial Finite Volume Method results are shown to further prove the effectiveness
- âŠ