Method of interior boundaries in a mixed problem of acoustic scattering

Directory of Open Access Journals (Sweden)

P. A. Krutitskii

1999-01-01

Full Text Available The mixed problem for the Helmholtz equation in the exterior of several bodies (obstacles is studied in 2 and 3 dimensions. The Dirichlet boundary condition is given on some obstacles and the impedance boundary condition is specified on the rest. The problem is investigated by a special modification of the boundary integral equation method. This modification can be called ‘Method of interior boundaries’, because additional boundaries are introduced inside scattering bodies, where impedance boundary condition is given. The solution of the problem is obtained in the form of potentials on the whole boundary. The density in the potentials satisfies the uniquely solvable Fredholm equation of the second kind and can be computed by standard codes. In fact our method holds for any positive wave numbers. The Neumann, Dirichlet, impedance problems and mixed Dirichlet–Neumann problem are particular cases of our problem.

Adaptive boundary conditions for exterior flow problems

CERN Document Server

Boenisch, V; Wittwer, S

2003-01-01

We consider the problem of solving numerically the stationary incompressible Navier-Stokes equations in an exterior domain in two dimensions. This corresponds to studying the stationary fluid flow past a body. The necessity to truncate for numerical purposes the infinite exterior domain to a finite domain leads to the problem of finding appropriate boundary conditions on the surface of the truncated domain. We solve this problem by providing a vector field describing the leading asymptotic behavior of the solution. This vector field is given in the form of an explicit expression depending on a real parameter. We show that this parameter can be determined from the total drag exerted on the body. Using this fact we set up a self-consistent numerical scheme that determines the parameter, and hence the boundary conditions and the drag, as part of the solution process. We compare the values of the drag obtained with our adaptive scheme with the results from using traditional constant boundary conditions. Computati...

Analytical Study of Microwave Heating of the Coal Layer With Mixed Boundary Conditions of I and III Types

Directory of Open Access Journals (Sweden)

Salomatov Vladimir

2016-01-01

Full Text Available This work is dedicated to the search for the exact analytical dependences of microwave heating due to absorption of a plane electromagnetic wave by coal layer with asymmetric and non-uniform heat dissipation conditions I and III kind. Some of simplifications have been made, such as one-dimensional problem, uniformity and isotropic coal material, and the constancy of the electrical properties of thermal coal during heating of microwave radiation. This has led to the fact that the Maxwell’s task is solved separately from the Fourier’s task, and a heat source generated in the carbon layer is subject to Bouguer law. For the system of equations of heat transfer has been found a new dependent variable, which is to simplify the search for a final solution. All this has given the possibility of finding rigorous analytical solution of the problem of microwave heating of the coal layer in the presence of asymmetric and inhomogeneous boundary conditions I and III kind.

Solution of moving boundary problems with implicit boundary condition

International Nuclear Information System (INIS)

Moyano, E.A.

1990-01-01

An algorithm that solves numerically a model for studying one dimensional moving boundary problems, with implicit boundary condition, is described. Landau's transformation is used, in order to work with a fixed number of nodes at each instant. Then, it is necessary to deal with a parabolic partial differential equation, whose diffusive and convective terms have variable coefficients. The partial differential equation is implicitly discretized, using Laasonen's scheme, always stable, instead of employing Crank-Nicholson sheme, as it has been done by Ferris and Hill. Fixed time and space steps (Δt, Δξ) are used, and the iteration is made with variable positions of the interface, i.e. varying δs until a boundary condition is satisfied. The model has the same features of the oxygen diffusion in absorbing tissue. It would be capable of estimating time variant radiation treatments of cancerous tumors. (Author) [es

Trace expansions for mixed boundary problems

Energy Technology Data Exchange (ETDEWEB)

Seeley, Robert T

2002-01-01

We discuss the heat trace expansion for a mixed boundary problem for the Laplace operator acting on sections of some bundle V over a manifold M of dimension d. The boundary is divided in two parts N{sub D} and N{sub N}, intersecting in a smooth submanifold {sigma}. Dirichlet conditions are imposed on N{sub D} - {sigma}, and Neumann conditions on N{sub N} - {sigma}. It turns out that it is also necessary to impose a condition along {sigma}. We then obtain an expansion of the trace of the heat operator with these boundary conditions, containing integrals of the usual terms over the interior and the two parts of the boundary, together with integrals over {sigma} of terms that are 'global' in certain operators on a semicircle. The first nonzero such term is computed; it involves the zeta function of an operator on the semicircle, and depends on the boundary condition along {sigma}. We find that no logarithmic terms occur in the expansion.

Solving fuzzy two-point boundary value problem using fuzzy Laplace transform

OpenAIRE

Ahmad, Latif; Farooq, Muhammad; Ullah, Saif; Abdullah, Saleem

2014-01-01

A natural way to model dynamic systems under uncertainty is to use fuzzy boundary value problems (FBVPs) and related uncertain systems. In this paper we use fuzzy Laplace transform to find the solution of two-point boundary value under generalized Hukuhara differentiability. We illustrate the method for the solution of the well known two-point boundary value problem Schrodinger equation, and homogeneous boundary value problem. Consequently, we investigate the solutions of FBVPs under as a ne...

Antireflective Boundary Conditions for Deblurring Problems

Directory of Open Access Journals (Sweden)

Marco Donatelli

2010-01-01

Full Text Available This survey paper deals with the use of antireflective boundary conditions for deblurring problems where the issues that we consider are the precision of the reconstruction when the noise is not present, the linear algebra related to these boundary conditions, the iterative and noniterative regularization solvers when the noise is considered, both from the viewpoint of the computational cost and from the viewpoint of the quality of the reconstruction. In the latter case, we consider a reblurring approach that replaces the transposition operation with correlation. For many of the considered items, the anti-reflective algebra coming from the given boundary conditions is the optimal choice. Numerical experiments corroborating the previous statement and a conclusion section end the paper.

On a variational principle for shape optimization and elliptic free boundary problems

Directory of Open Access Journals (Sweden)

Raúl B. González De Paz

2009-02-01

Full Text Available A variational principle for several free boundary value problems using a relaxation approach is presented. The relaxed Energy functional is concave and it is defined on a convex set, so that the minimizing points are characteristic functions of sets. As a consequence of the first order optimality conditions, it is shown that the corresponding sets are domains bounded by free boundaries, so that the equivalence of the solution of the relaxed problem with the solution of several free boundary value problem is proved. Keywords: Calculus of variations, optimization, free boundary problems.

Mixed Boundary Value Problem on Hypersurfaces

Directory of Open Access Journals (Sweden)

R. DuDuchava

2014-01-01

Full Text Available The purpose of the present paper is to investigate the mixed Dirichlet-Neumann boundary value problems for the anisotropic Laplace-Beltrami equation divC(A∇Cφ=f on a smooth hypersurface C with the boundary Γ=∂C in Rn. A(x is an n×n bounded measurable positive definite matrix function. The boundary is decomposed into two nonintersecting connected parts Γ=ΓD∪ΓN and on ΓD the Dirichlet boundary conditions are prescribed, while on ΓN the Neumann conditions. The unique solvability of the mixed BVP is proved, based upon the Green formulae and Lax-Milgram Lemma. Further, the existence of the fundamental solution to divS(A∇S is proved, which is interpreted as the invertibility of this operator in the setting Hp,#s(S→Hp,#s-2(S, where Hp,#s(S is a subspace of the Bessel potential space and consists of functions with mean value zero.

New Boundary Constraints for Elliptic Systems used in Grid Generation Problems

Science.gov (United States)

Kaul, Upender K.; Clancy, Daniel (Technical Monitor)

2002-01-01

This paper discusses new boundary constraints for elliptic partial differential equations as used in grid generation problems in generalized curvilinear coordinate systems. These constraints, based on the principle of local conservation of thermal energy in the vicinity of the boundaries, are derived using the Green's Theorem. They uniquely determine the so called decay parameters in the source terms of these elliptic systems. These constraints' are designed for boundary clustered grids where large gradients in physical quantities need to be resolved adequately. It is observed that the present formulation also works satisfactorily for mild clustering. Therefore, a closure for the decay parameter specification for elliptic grid generation problems has been provided resulting in a fully automated elliptic grid generation technique. Thus, there is no need for a parametric study of these decay parameters since the new constraints fix them uniquely. It is also shown that for Neumann type boundary conditions, these boundary constraints uniquely determine the solution to the internal elliptic problem thus eliminating the non-uniqueness of the solution of an internal Neumann boundary value grid generation problem.

Solution of a Problem Linear Plane Elasticity with Mixed Boundary Conditions by the Method of Boundary Integrals

Directory of Open Access Journals (Sweden)

Nahed S. Hussein

2014-01-01

Full Text Available A numerical boundary integral scheme is proposed for the solution to the system of …eld equations of plane. The stresses are prescribed on one-half of the circle, while the displacements are given. The considered problem with mixed boundary conditions in the circle is replaced by two problems with homogeneous boundary conditions, one of each type, having a common solution. The equations are reduced to a system of boundary integral equations, which is then discretized in the usual way, and the problem at this stage is reduced to the solution to a rectangular linear system of algebraic equations. The unknowns in this system of equations are the boundary values of four harmonic functions which define the full elastic solution and the unknown boundary values of stresses or displacements on proper parts of the boundary. On the basis of the obtained results, it is inferred that a stress component has a singularity at each of the two separation points, thought to be of logarithmic type. The results are discussed and boundary plots are given. We have also calculated the unknown functions in the bulk directly from the given boundary conditions using the boundary collocation method. The obtained results in the bulk are discussed and three-dimensional plots are given. A tentative form for the singular solution is proposed and the corresponding singular stresses and displacements are plotted in the bulk. The form of the singular tangential stress is seen to be compatible with the boundary values obtained earlier. The efficiency of the used numerical schemes is discussed.

Existence of solutions to boundary value problem of fractional differential equations with impulsive

Directory of Open Access Journals (Sweden)

Weihua JIANG

2016-12-01

Full Text Available In order to solve the boundary value problem of fractional impulsive differential equations with countable impulses and integral boundary conditions on the half line, the existence of solutions to the boundary problem is specifically studied. By defining suitable Banach spaces, norms and operators, using the properties of fractional calculus and applying the contraction mapping principle and Krasnoselskii's fixed point theorem, the existence of solutions for the boundary value problem of fractional impulsive differential equations with countable impulses and integral boundary conditions on the half line is proved, and examples are given to illustrate the existence of solutions to this kind of equation boundary value problems.

Positive solutions of nonlinear fractional boundary value problems with Dirichlet boundary conditions

Directory of Open Access Journals (Sweden)

Qingkai Kong

2012-02-01

Full Text Available In this paper, we study the existence and multiplicity of positive solutions of a class of nonlinear fractional boundary value problems with  Dirichlet boundary conditions. By applying the fixed point theory on cones we establish a series of criteria for the existence of one, two, any arbitrary finite number, and an infinite number of positive solutions. A criterion for the nonexistence of positive solutions is also derived. Several examples are given for demonstration.

The unified method: III. Nonlinearizable problems on the interval

International Nuclear Information System (INIS)

Lenells, J; Fokas, A S

2012-01-01

Boundary value problems for integrable nonlinear evolution PDEs formulated on the finite interval can be analyzed by the unified method introduced by one of the authors and extensively used in the literature. The implementation of this general method to this particular class of problems yields the solution in terms of the unique solution of a matrix Riemann–Hilbert problem formulated in the complex k-plane (the Fourier plane), which has a jump matrix with explicit (x, t)-dependence involving six scalar functions of k, called the spectral functions. Two of these functions depend on the initial data, whereas the other four depend on all boundary values. The most difficult step of the new method is the characterization of the latter four spectral functions in terms of the given initial and boundary data, i.e. the elimination of the unknown boundary values. Here, we present an effective characterization of the spectral functions in terms of the given initial and boundary data. We present two different characterizations of this problem. One is based on the analysis of the so-called global relation, on the analysis of the equations obtained from the global relation via certain transformations leaving the dispersion relation of the associated linearized PDE invariant and on the computation of the large k asymptotics of the eigenfunctions defining the relevant spectral functions. The other is based on the analysis of the global relation and on the introduction of the so-called Gelfand–Levitan–Marchenko representations of the eigenfunctions defining the relevant spectral functions. We also show that these two different characterizations are equivalent and that in the limit when the length of the interval tends to infinity, the relevant formulas reduce to the analogous formulas obtained recently for the case of boundary value problems formulated on the half-line. (paper)

Free boundary problems in PDEs and particle systems

CERN Document Server

Carinci, Gioia; Giardinà, Cristian; Presutti, Errico

2016-01-01

In this volume a theory for models of transport in the presence of a free boundary is developed. Macroscopic laws of transport are described by PDE's. When the system is open, there are several mechanisms to couple the system with the external forces. Here a class of systems where the interaction with the exterior takes place in correspondence of a free boundary is considered. Both continuous and discrete models sharing the same structure are analysed. In Part I a free boundary problem related to the Stefan Problem is worked out in all details. For this model a new notion of relaxed solution is proposed for which global existence and uniqueness is proven. It is also shown that this is the hydrodynamic limit of the empirical mass density of the associated particle system. In Part II several other models are discussed. The expectation is that the results proved for the basic model extend to these other cases. All the models discussed in this volume have an interest in problems arising in several research fields...

Mixed problems for linear symmetric hyperbolic systems with characteristic boundary conditions

International Nuclear Information System (INIS)

Secchi, P.

1994-01-01

We consider the initial-boundary value problem for symmetric hyperbolic systems with characteristic boundary of constant multiplicity. In the linear case we give some results about the existence of regular solutions in suitable functions spaces which take in account the loss of regularity in the normal direction to the characteristic boundary. We also consider the equations of ideal magneto-hydrodynamics under perfectly conducting wall boundary conditions and give some results about the solvability of such mixed problem. (author). 16 refs

A Boundary Value Problem for Introductory Physics?

Science.gov (United States)

Grundberg, Johan

2008-01-01

The Laplace equation has applications in several fields of physics, and problems involving this equation serve as paradigms for boundary value problems. In the case of the Laplace equation in a disc there is a well-known explicit formula for the solution: Poisson's integral. We show how one can derive this formula, and in addition two equivalent…

Boundary value problems and dichotomic stability

NARCIS (Netherlands)

England, R.; Mattheij, R.M.M.

1988-01-01

Since the conditioning of a boundary value problem (BVP) is closely related to the existence of a dichotomic fundamental solution (i.e., where one set of modes is increasing and a complementary set is decreasing), it is important to have discretization methods that conserve this dichotomy property.

Heading off boundary problems: clinical supervision as risk management.

Science.gov (United States)

Walker, R; Clark, J J

1999-11-01

The effective management of risk in clinical practice includes steps to limit harm to clients resulting from ethical violations or professional misconduct. Boundary problems constitute some of the most damaging ethical violations. The authors propose an active use of clinical supervision to anticipate and head off possible ethical violations by intervening when signs of boundary problems appear. The authors encourage a facilitative, Socratic method, rather than directive approaches, to help supervisees maximize their learning about ethical complexities. Building on the idea of a slippery slope, in which seemingly insignificant acts can lead to unethical patterns of behavior, the authors discuss ten cues to potential boundary problems, including strong feelings about a client; extended sessions with clients; gift giving between clinician and client; loans, barter, and sale of goods; clinician self-disclosures; and touching and sex. The authors outline supervisory interventions to be made when the cues are detected.

Separable boundary-value problems in physics

CERN Document Server

Willatzen, Morten

2011-01-01

Innovative developments in science and technology require a thorough knowledge of applied mathematics, particularly in the field of differential equations and special functions. These are relevant in modeling and computing applications of electromagnetic theory and quantum theory, e.g. in photonics and nanotechnology. The problem of solving partial differential equations remains an important topic that is taught at both the undergraduate and graduate level. Separable Boundary-Value Problems in Physics is an accessible and comprehensive treatment of partial differential equations i

Vragov’s boundary value problem for an implicit equation of mixed type

Science.gov (United States)

Egorov, I. E.

2017-10-01

We study a Vragov boundary value problem for a third-order implicit equation of mixed type with an arbitrary manifold of type switch. These Sobolev-type equations arise in many important applied problems. Given certain constraints on the coefficients and the right-hand side of the equation, we demonstrate, using nonstationary Galerkin method and regularization method, the unique regular solvability of the boundary value problem. We also obtain an error estimate for approximate solutions of the boundary value problem in terms of the regularization parameter and the eigenvalues of the Dirichlet spectral problem for the Laplace operator.

Analytic Solution to Shell Boundary – Value Problems

Directory of Open Access Journals (Sweden)

Yu. I. Vinogradov

2015-01-01

Full Text Available Object of research is to find analytical solution to the shell boundary – value problems, i.e. to consider the solution for a class of problems concerning the mechanics of hoop closed shells strain.The objective of work is to create an analytical method to define a stress – strain state of shells under non-axisymmetric loading. Thus, a main goal is to derive the formulas – solutions of the linear ordinary differential equations with variable continuous coefficients.The partial derivative differential equations of mechanics of shells strain by Fourier's method of variables division are reduced to the system of the differential equations with ordinary derivatives. The paper presents the obtained formulas to define solutions of the uniform differential equations and received on their basis formulas to define a particular solution depending on a type of the right parts of the differential equations.The analytical algorithm of the solution of a boundary task uses an approach to transfer the boundary conditions to the randomly chosen point of an interval of changing independent variable through the solution of the canonical matrix ordinary differential equation with the subsequent solution of system of algebraic equations for compatibility of boundary conditions at this point. Efficiency of algorithm is based on the fact that the solution of the ordinary differential equations is defined as the values of Cauchy – Krylova functions, which meet initial arbitrary conditions.The results of researches presented in work are useful to experts in the field of calculus mathematics, dealing with solution of systems of linear ordinary differential equations and creation of effective analytical computing methods to solve shell boundary – value problems.

A fast direct solver for boundary value problems on locally perturbed geometries

Science.gov (United States)

Zhang, Yabin; Gillman, Adrianna

2018-03-01

Many applications including optimal design and adaptive discretization techniques involve solving several boundary value problems on geometries that are local perturbations of an original geometry. This manuscript presents a fast direct solver for boundary value problems that are recast as boundary integral equations. The idea is to write the discretized boundary integral equation on a new geometry as a low rank update to the discretized problem on the original geometry. Using the Sherman-Morrison formula, the inverse can be expressed in terms of the inverse of the original system applied to the low rank factors and the right hand side. Numerical results illustrate for problems where perturbation is localized the fast direct solver is three times faster than building a new solver from scratch.

A finite difference method for free boundary problems

KAUST Repository

Fornberg, Bengt

2010-04-01

Fornberg and Meyer-Spasche proposed some time ago a simple strategy to correct finite difference schemes in the presence of a free boundary that cuts across a Cartesian grid. We show here how this procedure can be combined with a minimax-based optimization procedure to rapidly solve a wide range of elliptic-type free boundary value problems. © 2009 Elsevier B.V. All rights reserved.

Multidimensional phase change problems by the dual-reciprocity boundary-element method

International Nuclear Information System (INIS)

Jo, J.C.; Shin, W.K.; Choi, C.Y.

1999-01-01

Transient heat transfer problems with phase changes (Stefan problems) occur in many engineering situations, including potential core melting and solidification during pressurized-water-reactor severe accidents, ablation of thermal shields, melting and solidification of alloys, and many others. This article addresses the numerical analysis of nonlinear transient heat transfer with melting or solidification. An effective and simple procedure is presented for the simulation of the motion of the boundary and the transient temperature field during the phase change process. To accomplish this purpose, an iterative implicit solution algorithm has been developed by employing the dual-reciprocity boundary-element method. The dual-reciprocity boundary-element approach provided in this article is much simpler than the usual boundary-element method in applying a reciprocity principle and an available technique for dealing with the domain integral of the boundary element formulation simultaneously. In this article, attention is focused on two-dimensional melting (ablation)/solidification problems for simplicity. The accuracy and effectiveness of the present analysis method have been illustrated through comparisons of the calculation results of some examples of one-phase ablation/solidification problems with their known semianalytical or numerical solutions where available

Group invariance in engineering boundary value problems

CERN Document Server

Seshadri, R

1985-01-01

REFEREN CES . 156 9 Transforma.tion of a Boundary Value Problem to an Initial Value Problem . 157 9.0 Introduction . 157 9.1 Blasius Equation in Boundary Layer Flow . 157 9.2 Longitudinal Impact of Nonlinear Viscoplastic Rods . 163 9.3 Summary . 168 REFERENCES . . . . . . . . . . . . . . . . . . 168 . 10 From Nonlinear to Linear Differential Equa.tions Using Transformation Groups. . . . . . . . . . . . . . 169 . 10.1 From Nonlinear to Linear Differential Equations . 170 10.2 Application to Ordinary Differential Equations -Bernoulli's Equation . . . . . . . . . . . 173 10.3 Application to Partial Differential Equations -A Nonlinear Chemical Exchange Process . 178 10.4 Limitations of the Inspectional Group Method . 187 10.5 Summary . 188 REFERENCES . . . . 188 11 Miscellaneous Topics . 190 11.1 Reduction of Differential Equations to Algebraic Equations 190 11.2 Reduction of Order of an Ordinary Differential Equation . 191 11.3 Transformat.ion From Ordinary to Partial Differential Equations-Search for First Inte...

A new Ellipsoidal Gravimetric-Satellite Altimetry Boundary Value Problem; Case study: High Resolution Geoid of Iran

Science.gov (United States)

Ardalan, A.; Safari, A.; Grafarend, E.

2003-04-01

A new ellipsoidal gravimetric-satellite altimetry boundary value problem has been developed and successfully tested. This boundary value problem has been constructed for gravity observables of the type (i) gravity potential (ii) gravity intensity (iii) deflection of vertical and (iv) satellite altimetry data. The developed boundary value problem is enjoying the ellipsoidal nature and as such can take advantage of high precision GPS observations in the set-up of the problem. The highlights of the solution are as follows: begin{itemize} Application of ellipsoidal harmonic expansion up to degree/order and ellipsoidal centrifugal field for the reduction of global gravity and isostasy effects from the gravity observable at the surface of the Earth. Application of ellipsoidal Newton integral on the equal area map projection surface for the reduction of residual mass effects within a radius of 55 km around the computational point. Ellipsoidal harmonic downward continuation of the residual observables from the surface of the earth down to the surface of reference ellipsoid using the ellipsoidal height of the observation points derived from GPS. Restore of the removed effects at the application points on the surface of reference ellipsoid. Conversion of the satellite altimetry derived heights of the water bodies into potential. Combination of the downward continued gravity information with the potential equivalent of the satellite altimetry derived heights of the water bodies. Application of ellipsoidal Bruns formula for converting the potential values on the surface of the reference ellipsoid into the geoidal heights (i.e. ellipsoidal heights of the geoid) with respect to the reference ellipsoid. Computation of the high-resolution geoid of Iran has successfully tested this new methodology!

Boundary element methods applied to two-dimensional neutron diffusion problems

International Nuclear Information System (INIS)

Itagaki, Masafumi

1985-01-01

The Boundary element method (BEM) has been applied to two-dimensional neutron diffusion problems. The boundary integral equation and its discretized form have been derived. Some numerical techniques have been developed, which can be applied to critical and fixed-source problems including multi-region ones. Two types of test programs have been developed according to whether the 'zero-determinant search' or the 'source iteration' technique is adopted for criticality search. Both programs require only the fluxes and currents on boundaries as the unknown variables. The former allows a reduction in computing time and memory in comparison with the finite element method (FEM). The latter is not always efficient in terms of computing time due to the domain integral related to the inhomogeneous source term; however, this domain integral can be replaced by the equivalent boundary integral for a region with a non-multiplying medium or with a uniform source, resulting in a significant reduction in computing time. The BEM, as well as the FEM, is well suited for solving irregular geometrical problems for which the finite difference method (FDM) is unsuited. The BEM also solves problems with infinite domains, which cannot be solved by the ordinary FEM and FDM. Some simple test calculations are made to compare the BEM with the FEM and FDM, and discussions are made concerning the relative merits of the BEM and problems requiring future solution. (author)

The vanishing discount problem and viscosity Mather measures. Part 2: boundary value problems

OpenAIRE

Ishii, Hitoshi; Mitake, Hiroyoshi; Tran, Hung V.

2016-01-01

In arXiv:1603.01051 (Part 1 of this series), we have introduced a variational approach to studying the vanishing discount problem for fully nonlinear, degenerate elliptic, partial differential equations in a torus. We develop this approach further here to handle boundary value problems. In particular, we establish new representation formulas for solutions of discount problems, critical values, and use them to prove convergence results for the vanishing discount problems.

Boundary value problems of holomorphic vector functions in 1D QCs

International Nuclear Information System (INIS)

Gao Yang; Zhao Yingtao; Zhao Baosheng

2007-01-01

By means of the generalized Stroh formalism, two-dimensional (2D) problems of one-dimensional (1D) quasicrystals (QCs) elasticity are turned into the boundary value problems of holomorphic vector functions in a given region. If the conformal mapping from an ellipse to a circle is known, a general method for solving the boundary value problems of holomorphic vector functions can be presented. To illustrate its utility, by using the necessary and sufficient condition of boundary value problems of holomorphic vector functions, we consider two basic 2D problems in 1D QCs, that is, an elliptic hole and a rigid line inclusion subjected to uniform loading at infinity. For the crack problem, the intensity factors of phonon and phason fields are determined, and the physical sense of the results relative to phason and the difference between mechanical behaviors of the crack problem in crystals and QCs are figured out. Moreover, the same procedure can be used to deal with the elastic problems for 2D and three-dimensional (3D) QCs

Laplace boundary-value problem in paraboloidal coordinates

International Nuclear Information System (INIS)

Duggen, L; Willatzen, M; Voon, L C Lew Yan

2012-01-01

This paper illustrates both a problem in mathematical physics, whereby the method of separation of variables, while applicable, leads to three ordinary differential equations that remain fully coupled via two separation constants and a five-term recurrence relation for series solutions, and an exactly solvable problem in electrostatics, as a boundary-value problem on a paraboloidal surface. In spite of the complex nature of the former, it is shown that the latter solution can be quite simple. Results are provided for the equipotential surfaces and electric field lines are given near a paraboloidal conductor. (paper)

A finite-volume method for convection problems with embedded moving boundaries

NARCIS (Netherlands)

Y.J. Hassen (Yunus); B. Koren (Barry)

2009-01-01

htmlabstractAn accurate method, using a novel immersed-boundary approach, is presented for numerically solving linear, scalar convection problems. Moving interior boundary conditions are embedded in the fixed-grid fluxes in the direct neighborhood of the moving boundaries. Tailor-made limiters are

Solving free-plasma-boundary problems with the SIESTA MHD code

Science.gov (United States)

Sanchez, R.; Peraza-Rodriguez, H.; Reynolds-Barredo, J. M.; Tribaldos, V.; Geiger, J.; Hirshman, S. P.; Cianciosa, M.

2017-10-01

SIESTA is a recently developed MHD equilibrium code designed to perform fast and accurate calculations of ideal MHD equilibria for 3D magnetic configurations. It is an iterative code that uses the solution obtained by the VMEC code to provide a background coordinate system and an initial guess of the solution. The final solution that SIESTA finds can exhibit magnetic islands and stochastic regions. In its original implementation, SIESTA addressed only fixed-boundary problems. This fixed boundary condition somewhat restricts its possible applications. In this contribution we describe a recent extension of SIESTA that enables it to address free-plasma-boundary situations, opening up the possibility of investigating problems with SIESTA in which the plasma boundary is perturbed either externally or internally. As an illustration, the extended version of SIESTA is applied to a configuration of the W7-X stellarator.

Simulating a singularity-free universe outside the problem boundary in poisson

International Nuclear Information System (INIS)

Halbach, K.; Schlueter, R.

1992-01-01

An exact analytical solution developed from the Dirichlet problem exterior to a circle is employed in the magnetostatics code POISSON to provide a boundary condition option which simulates a singularity-free universe external to the problem domain. Problems with domains of large unequal extents in perpendicular directions are treated by first conformally mapping the exterior of an ellipse onto the exterior of the unit circle. Problems exhibiting symmetry in one or two planes are modeled using a semi or quarter, respectively, in conjunction with the singularity-free rest-of-universe boundary condition

The boundary value problems of magnetotail equilibrium

International Nuclear Information System (INIS)

Birn, J.

1991-01-01

The equilibrium problem for the Earth's magnetotail is discussed under the assumption that the boundary of the tail can be prescribed or derived from the force balance with the solar wind. A general solution of this problem is presented for the two-dimensional case, where the dependence on the γ coordinate and the presence of Β gamma are neglected. These solutions are further generalized to include the γ dependence (but no Β gamma ) and an open magnetopause. In this formulation, a solution can be obtained by integration when the magnetopause boundary α(x,y), the total pressure function p(x), and the magnetic flux distribution A b (x,y) at the magnetopause are prescribed. Certain restrictions, however, may limit the free choice of these functions to yield physically reasonable, real solutions. When the interaction with the solar wind is included, the boundary location can no longer be chosen freely but follows from the force balance of the magnetotail with the solar wind. For a simplified description of this force balance a differential equation for the boundary location is derived, which generalizes an earlier result by Coroniti and Kennel (1972). It is shown that solutions of this differential equation are bounded by a maximum tail width if the plasma sheet thickness is limited. Several explicit solutions are presented, illustrating cases with and without tail flaring in the z direction, and including the restrictions of the force balance with the solar wind and of the conservation laws of adiabatic convection in a steady configuration

Positive solutions and eigenvalues of nonlocal boundary-value problems

Directory of Open Access Journals (Sweden)

Jifeng Chu

2005-07-01

Full Text Available We study the ordinary differential equation $x''+lambda a(tf(x=0$ with the boundary conditions $x(0=0$ and $x'(1=int_{eta}^{1}x'(sdg(s$. We characterize values of $lambda$ for which boundary-value problem has a positive solution. Also we find appropriate intervals for $lambda$ so that there are two positive solutions.

On a non-linear pseudodifferential boundary value problem

International Nuclear Information System (INIS)

Nguyen Minh Chuong.

1989-12-01

A pseudodifferential boundary value problem for operators with symbols taking values in Sobolev spaces and with non-linear right-hand side was studied. Existence and uniqueness theorems were proved. (author). 11 refs

Boundary integral equation methods in eigenvalue problems of elastodynamics and thin plates

CERN Document Server

Kitahara, M

1985-01-01

The boundary integral equation (BIE) method has been used more and more in the last 20 years for solving various engineering problems. It has important advantages over other techniques for numerical treatment of a wide class of boundary value problems and is now regarded as an indispensable tool for potential problems, electromagnetism problems, heat transfer, fluid flow, elastostatics, stress concentration and fracture problems, geomechanical problems, and steady-state and transient electrodynamics.In this book, the author gives a complete, thorough and detailed survey of the method. It pro

Efficient algorithms for analyzing the singularly perturbed boundary value problems of fractional order

Science.gov (United States)

Sayevand, K.; Pichaghchi, K.

2018-04-01

In this paper, we were concerned with the description of the singularly perturbed boundary value problems in the scope of fractional calculus. We should mention that, one of the main methods used to solve these problems in classical calculus is the so-called matched asymptotic expansion method. However we shall note that, this was not achievable via the existing classical definitions of fractional derivative, because they do not obey the chain rule which one of the key elements of the matched asymptotic expansion method. In order to accommodate this method to fractional derivative, we employ a relatively new derivative so-called the local fractional derivative. Using the properties of local fractional derivative, we extend the matched asymptotic expansion method to the scope of fractional calculus and introduce a reliable new algorithm to develop approximate solutions of the singularly perturbed boundary value problems of fractional order. In the new method, the original problem is partitioned into inner and outer solution equations. The reduced equation is solved with suitable boundary conditions which provide the terminal boundary conditions for the boundary layer correction. The inner solution problem is next solved as a solvable boundary value problem. The width of the boundary layer is approximated using appropriate resemblance function. Some theoretical results are established and proved. Some illustrating examples are solved and the results are compared with those of matched asymptotic expansion method and homotopy analysis method to demonstrate the accuracy and efficiency of the method. It can be observed that, the proposed method approximates the exact solution very well not only in the boundary layer, but also away from the layer.

Zero-energy eigenstates for the Dirac boundary problem

International Nuclear Information System (INIS)

Hortacsu, M.; Rothe, K.D.; Schroer, B.

1980-01-01

As an alternative to the method of spherical compactification for the Dirac operator in instanton background fields we study the correct method of 'box-quantization': the Atiyah-Patodi-Singer spectral boundary condition. This is the only self-adjoint boundary condition which respects the charge conjugation property and the γ 5 symmetry, apart form the usual breaking due to zero modes. We point out the relevance of this approach to the computation of instanton determinants and other problems involving Dirac spinors. (orig.)

Electromagnetic wave theory for boundary-value problems an advanced course on analytical methods

CERN Document Server

Eom, Hyo J

2004-01-01

Electromagnetic wave theory is based on Maxwell's equations, and electromagnetic boundary-value problems must be solved to understand electromagnetic scattering, propagation, and radiation. Electromagnetic theory finds practical applications in wireless telecommunications and microwave engineering. This book is written as a text for a two-semester graduate course on electromagnetic wave theory. As such, Electromagnetic Wave Theory for Boundary-Value Problems is intended to help students enhance analytic skills by solving pertinent boundary-value problems. In particular, the techniques of Fourier transform, mode matching, and residue calculus are utilized to solve some canonical scattering and radiation problems.

Existence results for anisotropic discrete boundary value problems

Directory of Open Access Journals (Sweden)

Avci Avci

2016-06-01

Full Text Available In this article, we prove the existence of nontrivial weak solutions for a class of discrete boundary value problems. The main tools used here are the variational principle and critical point theory.

Homology in Electromagnetic Boundary Value Problems

Directory of Open Access Journals (Sweden)

Pellikka Matti

2010-01-01

Full Text Available We discuss how homology computation can be exploited in computational electromagnetism. We represent various cellular mesh reduction techniques, which enable the computation of generators of homology spaces in an acceptable time. Furthermore, we show how the generators can be used for setting up and analysis of an electromagnetic boundary value problem. The aim is to provide a rationale for homology computation in electromagnetic modeling software.

Cost-effective computations with boundary interface operators in elliptic problems

International Nuclear Information System (INIS)

Khoromskij, B.N.; Mazurkevich, G.E.; Nikonov, E.G.

1993-01-01

The numerical algorithm for fast computations with interface operators associated with the elliptic boundary value problems (BVP) defined on step-type domains is presented. The algorithm is based on the asymptotically almost optimal technique developed for treatment of the discrete Poincare-Steklov (PS) operators associated with the finite-difference Laplacian on rectangles when using the uniform grid with a 'displacement by h/2'. The approach can be regarded as an extension of the method proposed for the partial solution of the finite-difference Laplace equation to the case of displaced grids and mixed boundary conditions. It is shown that the action of the PS operator for the Dirichlet problem and mixed BVP can be computed with expenses of the order of O(Nlog 2 N) both for arithmetical operations and computer memory needs, where N is the number of unknowns on the rectangle boundary. The single domain algorithm is applied to solving the multidomain elliptic interface problems with piecewise constant coefficients. The numerical experiments presented confirm almost linear growth of the computational costs and memory needs with respect to the dimension of the discrete interface problem. 14 refs., 3 figs., 4 tabs

IMPSOR, 3-D Boundary Problems Solution for Thermal Conductivity Calculation

International Nuclear Information System (INIS)

Wilson, D.G.; Williams, M.A.

1994-01-01

1 - Description of program or function: IMPSOR implements finite difference methods for multidimensional moving boundary problems with Dirichlet or Neumann boundary conditions. The geometry of the spatial domain is a rectangular parallelepiped with dimensions specified by the user. Dirichlet or Neumann boundary conditions may be specified on each face of the box independently. The user defines the initial and boundary conditions as well as the thermal and physical properties of the problem and several parameters for the numerical method, e.g. degree of implicitness, time-step size. 2 - Method of solution: The spatial domain is partitioned and the governing equation discretized, which yields a nonlinear system of equations at each time-step. This nonlinear system is solved using a successive over-relaxation (SOR) algorithm. For a given node, the previous iteration's temperature and thermal conductivity values are used for advanced points with current values at previous points. This constitutes a Gauss-Seidel iteration. Most of the computing time used by the numerical method is spent in the iterative solution of the nonlinear system. The SOR scheme employed is designed to accommodate vectorization on a Cray X-MP. 3 - Restrictions on the complexity of the problem: Maximum of 70,000 nodes

Initial-boundary value problems associated with the Ablowitz-Ladik system

Science.gov (United States)

Xia, Baoqiang; Fokas, A. S.

2018-02-01

We employ the Ablowitz-Ladik system as an illustrative example in order to demonstrate how to analyze initial-boundary value problems for integrable nonlinear differential-difference equations via the unified transform (Fokas method). In particular, we express the solutions of the integrable discrete nonlinear Schrödinger and integrable discrete modified Korteweg-de Vries equations in terms of the solutions of appropriate matrix Riemann-Hilbert problems. We also discuss in detail, for both the above discrete integrable equations, the associated global relations and the process of eliminating of the unknown boundary values.

Application of He's variational iteration method to the fifth-order boundary value problems

International Nuclear Information System (INIS)

Shen, S

2008-01-01

Variational iteration method is introduced to solve the fifth-order boundary value problems. This method provides an efficient approach to solve this type of problems without discretization and the computation of the Adomian polynomials. Numerical results demonstrate that this method is a promising and powerful tool for solving the fifth-order boundary value problems

A nonlinear free boundary problem with a self-driven Bernoulli condition

OpenAIRE

Dipierro, Serena; Karakhanyan, Aram; Valdinoci, Enrico

2017-01-01

We study a Bernoulli type free boundary problem with two phases J[u]=∫Ω|∇u(x)|2dx+Φ(M−(u),M+(u)),u−u¯∈W1,20(Ω), where u¯∈W1,2(Ω) is a given boundary datum. Here, M1 and M2 are weighted volumes of {u≤0}∩Ω and {u>0}∩Ω, respectively, and Φ is a nonnegative function of two real variables. We show that, for this problem, the Bernoulli constant, which determines the gradient jump condition across the free boundary, is of global type and it is indeed determined by the weighted volumes of the phas...

Analytic approximations to nonlinear boundary value problems modeling beam-type nano-electromechanical systems

Energy Technology Data Exchange (ETDEWEB)

Zou, Li [Dalian Univ. of Technology, Dalian City (China). State Key Lab. of Structural Analysis for Industrial Equipment; Liang, Songxin; Li, Yawei [Dalian Univ. of Technology, Dalian City (China). School of Mathematical Sciences; Jeffrey, David J. [Univ. of Western Ontario, London (Canada). Dept. of Applied Mathematics

2017-06-01

Nonlinear boundary value problems arise frequently in physical and mechanical sciences. An effective analytic approach with two parameters is first proposed for solving nonlinear boundary value problems. It is demonstrated that solutions given by the two-parameter method are more accurate than solutions given by the Adomian decomposition method (ADM). It is further demonstrated that solutions given by the ADM can also be recovered from the solutions given by the two-parameter method. The effectiveness of this method is demonstrated by solving some nonlinear boundary value problems modeling beam-type nano-electromechanical systems.

A Duality Approach for the Boundary Variation of Neumann Problems

DEFF Research Database (Denmark)

Bucur, Dorin; Varchon, Nicolas

2002-01-01

In two dimensions, we study the stability of the solution of an elliptic equation with Neumann boundary conditions for nonsmooth perturbations of the geometric domain. Using harmonic conjugates, we relate this problem to the shape stability of the solution of an elliptic equation with Dirichlet b...... boundary conditions. As a particular case, we prove the stability of the solution under a topological constraint ( uniform number of holes), which is analogous to Sverak's result for Dirichlet boundary conditions....

A duality approach or the boundary variation of Neumann problems

DEFF Research Database (Denmark)

Bucur, D.; Varchon, Nicolas

2002-01-01

In two dimensions, we study the stability of the solution of an elliptic equation with Neumann boundary conditions for nonsmooth perturbations of the geometric domain. Using harmonic conjugates, we relate this problem to the shape stability of the solution of an elliptic equation with Dirichlet b...... boundary conditions. As a particular case, we prove the stability of the solution under a topological constraint ( uniform number of holes), which is analogous to Sverak's result for Dirichlet boundary conditions....

Asymptotics of linear initial boundary value problems with periodic boundary data on the half-line and finite intervals

KAUST Repository

Dujardin, G. M.

2009-01-01

This paper deals with the asymptotic behaviour of the solutions of linear initial boundary value problems with constant coefficients on the half-line and on finite intervals. We assume that the boundary data are periodic in time and we investigate

A boundary-value problem in weighted Hölder spaces for elliptic equations which degenerate at the boundary of the domain

Energy Technology Data Exchange (ETDEWEB)

Bazalii, B V; Degtyarev, S P [Institute of Applied Mathematics and Mechanics, National Academy of Sciences of Ukraine, Donetsk (Ukraine)

2013-07-31

An elliptic boundary-value problem for second-order equations with nonnegative characteristic form is investigated in the situation when there is a weak degeneracy on the boundary of the domain. A priori estimates are obtained for solutions and the problem is proved to be solvable in some weighted Hölder spaces. Bibliography: 18 titles.

Positive Solutions of Two-Point Boundary Value Problems for Monge-Ampère Equations

Directory of Open Access Journals (Sweden)

Baoqiang Yan

2015-01-01

Full Text Available This paper considers the following boundary value problem: ((-u'(tn'=ntn-1f(u(t,  01 is odd. We establish the method of lower and upper solutions for some boundary value problems which generalizes the above equations and using this method we present a necessary and sufficient condition for the existence of positive solutions to the above boundary value problem and some sufficient conditions for the existence of positive solutions.

The Boundary Element Method Applied to the Two Dimensional Stefan Moving Boundary Problem

Science.gov (United States)

1991-03-15

Unc), - ( UGt )t - (UG,,),,] - (UG), If we integrate this equation with respect to r from 0 to t - c and with respect to and ij on the region 11(r...and others. "Moving Boundary Problems in Phase Change Mod- els," SIGNUM Newsletter, 20: 8-12 (1985). 21. Stefan, J. "Ober einige Probleme der Theorie ...ier Wirmelcitung," S.-B. \\Vein. Akad. Mat. Natur., 98: 173-484 (1889). 22.-. "flber (lie Theorie der Eisbildung insbesondere fiber die lisbildung im

Fourth-order discrete anisotropic boundary-value problems

Directory of Open Access Journals (Sweden)

Maciej Leszczynski

2015-09-01

Full Text Available In this article we consider the fourth-order discrete anisotropic boundary value problem with both advance and retardation. We apply the direct method of the calculus of variations and the mountain pass technique to prove the existence of at least one and at least two solutions. Non-existence of non-trivial solutions is also undertaken.

Boundary Element Solution of Geometrical Inverse Heat Conduction Problems for Development of IR CAT Scan

International Nuclear Information System (INIS)

Choi, C. Y.; Park, C. T.; Kim, T. H.; Han, K. N.; Choe, S. H.

1995-01-01

A geometrical inverse heat conduction problem is solved for the development of Infrared Computerized-Axial-Tomography (IR CAT) Scan by using a boundary element method in conjunction with regularization procedure. In this problem, an overspecified temperature condition by infrared scanning is provided on the surface, and is used together with other conditions to solve the position of an unknown boundary (cavity). An auxiliary problem is introduced in the solution of this problem. By defining a hypothetical inner boundary for the auxiliary problem domain, the cavity is located interior to the domain and its position is determined by solving a potential problem. Boundary element method with regularization procedure is used to solve this problem, and the effects of regularization on the inverse solution method are investigated by means of numerical analysis

Discontinuous Sturm-Liouville Problems with Eigenvalue Dependent Boundary Condition

Energy Technology Data Exchange (ETDEWEB)

Amirov, R. Kh., E-mail: emirov@cumhuriyet.edu.tr; Ozkan, A. S., E-mail: sozkan@cumhuriyet.edu.tr [Cumhuriyet University, Department of Mathematics Faculty of Art and Science (Turkey)

2014-12-15

In this study, an inverse problem for Sturm-Liouville differential operators with discontinuities is studied when an eigenparameter appears not only in the differential equation but it also appears in the boundary condition. Uniqueness theorems of inverse problems according to the Prüfer angle, the Weyl function and two different eigenvalues sets are proved.

Extended Subject Access to Hypertext Online Documentation. Part III: The Document-Boundaries Problem.

Science.gov (United States)

Girill, T. R.

1991-01-01

This article continues the description of DFT (Document, Find, Theseus), an online documentation system that provides computer-managed on-demand printing of software manuals as well as the interactive retrieval of reference passages. Document boundaries in the hypertext database are discussed, search vocabulary complexities are described, and text…

A combined analytic-numeric approach for some boundary-value problems

Directory of Open Access Journals (Sweden)

Mustafa Turkyilmazoglu

2016-02-01

Full Text Available A combined analytic-numeric approach is undertaken in the present work for the solution of boundary-value problems in the finite or semi-infinite domains. Equations to be treated arise specifically from the boundary layer analysis of some two and three-dimensional flows in fluid mechanics. The purpose is to find quick but accurate enough solutions. Taylor expansions at either boundary conditions are computed which are next matched to the other asymptotic or exact boundary conditions. The technique is applied to the well-known Blasius as well as Karman flows. Solutions obtained in terms of series compare favorably with the existing ones in the literature.

Application of Monte Carlo method to solving boundary value problem of differential equations

International Nuclear Information System (INIS)

Zuo Yinghong; Wang Jianguo

2012-01-01

This paper introduces the foundation of the Monte Carlo method and the way how to generate the random numbers. Based on the basic thought of the Monte Carlo method and finite differential method, the stochastic model for solving the boundary value problem of differential equations is built. To investigate the application of the Monte Carlo method to solving the boundary value problem of differential equations, the model is used to solve Laplace's equations with the first boundary condition and the unsteady heat transfer equation with initial values and boundary conditions. The results show that the boundary value problem of differential equations can be effectively solved with the Monte Carlo method, and the differential equations with initial condition can also be calculated by using a stochastic probability model which is based on the time-domain finite differential equations. Both the simulation results and theoretical analyses show that the errors of numerical results are lowered as the number of simulation particles is increased. (authors)

Continuum and Discrete Initial-Boundary Value Problems and Einstein's Field Equations

Directory of Open Access Journals (Sweden)

Olivier Sarbach

2012-08-01

Full Text Available Many evolution problems in physics are described by partial differential equations on an infinite domain; therefore, one is interested in the solutions to such problems for a given initial dataset. A prominent example is the binary black-hole problem within Einstein's theory of gravitation, in which one computes the gravitational radiation emitted from the inspiral of the two black holes, merger and ringdown. Powerful mathematical tools can be used to establish qualitative statements about the solutions, such as their existence, uniqueness, continuous dependence on the initial data, or their asymptotic behavior over large time scales. However, one is often interested in computing the solution itself, and unless the partial differential equation is very simple, or the initial data possesses a high degree of symmetry, this computation requires approximation by numerical discretization. When solving such discrete problems on a machine, one is faced with a finite limit to computational resources, which leads to the replacement of the infinite continuum domain with a finite computer grid. This, in turn, leads to a discrete initial-boundary value problem. The hope is to recover, with high accuracy, the exact solution in the limit where the grid spacing converges to zero with the boundary being pushed to infinity. The goal of this article is to review some of the theory necessary to understand the continuum and discrete initial boundary-value problems arising from hyperbolic partial differential equations and to discuss its applications to numerical relativity; in particular, we present well-posed initial and initial-boundary value formulations of Einstein's equations, and we discuss multi-domain high-order finite difference and spectral methods to solve them.

Continuum and Discrete Initial-Boundary Value Problems and Einstein's Field Equations.

Science.gov (United States)

Sarbach, Olivier; Tiglio, Manuel

2012-01-01

Many evolution problems in physics are described by partial differential equations on an infinite domain; therefore, one is interested in the solutions to such problems for a given initial dataset. A prominent example is the binary black-hole problem within Einstein's theory of gravitation, in which one computes the gravitational radiation emitted from the inspiral of the two black holes, merger and ringdown. Powerful mathematical tools can be used to establish qualitative statements about the solutions, such as their existence, uniqueness, continuous dependence on the initial data, or their asymptotic behavior over large time scales. However, one is often interested in computing the solution itself, and unless the partial differential equation is very simple, or the initial data possesses a high degree of symmetry, this computation requires approximation by numerical discretization. When solving such discrete problems on a machine, one is faced with a finite limit to computational resources, which leads to the replacement of the infinite continuum domain with a finite computer grid. This, in turn, leads to a discrete initial-boundary value problem. The hope is to recover, with high accuracy, the exact solution in the limit where the grid spacing converges to zero with the boundary being pushed to infinity. The goal of this article is to review some of the theory necessary to understand the continuum and discrete initial boundary-value problems arising from hyperbolic partial differential equations and to discuss its applications to numerical relativity; in particular, we present well-posed initial and initial-boundary value formulations of Einstein's equations, and we discuss multi-domain high-order finite difference and spectral methods to solve them.

Analytic Approximations to the Free Boundary and Multi-dimensional Problems in Financial Derivatives Pricing

Science.gov (United States)

Lau, Chun Sing

This thesis studies two types of problems in financial derivatives pricing. The first type is the free boundary problem, which can be formulated as a partial differential equation (PDE) subject to a set of free boundary condition. Although the functional form of the free boundary condition is given explicitly, the location of the free boundary is unknown and can only be determined implicitly by imposing continuity conditions on the solution. Two specific problems are studied in details, namely the valuation of fixed-rate mortgages and CEV American options. The second type is the multi-dimensional problem, which involves multiple correlated stochastic variables and their governing PDE. One typical problem we focus on is the valuation of basket-spread options, whose underlying asset prices are driven by correlated geometric Brownian motions (GBMs). Analytic approximate solutions are derived for each of these three problems. For each of the two free boundary problems, we propose a parametric moving boundary to approximate the unknown free boundary, so that the original problem transforms into a moving boundary problem which can be solved analytically. The governing parameter of the moving boundary is determined by imposing the first derivative continuity condition on the solution. The analytic form of the solution allows the price and the hedging parameters to be computed very efficiently. When compared against the benchmark finite-difference method, the computational time is significantly reduced without compromising the accuracy. The multi-stage scheme further allows the approximate results to systematically converge to the benchmark results as one recasts the moving boundary into a piecewise smooth continuous function. For the multi-dimensional problem, we generalize the Kirk (1995) approximate two-asset spread option formula to the case of multi-asset basket-spread option. Since the final formula is in closed form, all the hedging parameters can also be derived in

m-POINT BOUNDARY VALUE PROBLEM FOR SECOND ORDER IMPULSIVE DIFFERENTIAL EQUATION AT RESONANCE

Institute of Scientific and Technical Information of China (English)

无

2012-01-01

In his paper,we obtain a general theorem concerning the existence of solutions to an m-point boundary value problem for the second-order differential equation with impulses.Moreover,the result can also be applied to study the usual m-point boundary value problem at resonance without impulses.

Application of the perturbation iteration method to boundary layer type problems.

Science.gov (United States)

Pakdemirli, Mehmet

2016-01-01

The recently developed perturbation iteration method is applied to boundary layer type singular problems for the first time. As a preliminary work on the topic, the simplest algorithm of PIA(1,1) is employed in the calculations. Linear and nonlinear problems are solved to outline the basic ideas of the new solution technique. The inner and outer solutions are determined with the iteration algorithm and matched to construct a composite expansion valid within all parts of the domain. The solutions are contrasted with the available exact or numerical solutions. It is shown that the perturbation-iteration algorithm can be effectively used for solving boundary layer type problems.

APPLICATION OF BOUNDARY INTEGRAL EQUATION METHOD FOR THERMOELASTICITY PROBLEMS

Directory of Open Access Journals (Sweden)

Vorona Yu.V.

2015-12-01

Full Text Available Boundary Integral Equation Method is used for solving analytically the problems of coupled thermoelastic spherical wave propagation. The resulting mathematical expressions coincide with the solutions obtained in a conventional manner.

three solutions for a semilinear elliptic boundary value problem

Indian Academy of Sciences (India)

69

Keywords: The Laplacian operator, elliptic problem, Nehari man- ifold, three critical points, weak solution. 1. Introduction. Let Ω be a smooth bounded domain in RN , N ≥ 3 . In this work, we show the existence of at least three solutions for the semilinear elliptic boundary- value problem: (Pλ).. −∆u = f(x)|u(x)|p−2u(x) + ...

Parametrices and exact paralinearization of semi-linear boundary problems

DEFF Research Database (Denmark)

Johnsen, Jon

2008-01-01

The subject is parametrices for semi-linear problems, based on parametrices for linear boundary problems and on non-linearities that decompose into solution-dependent linear operators acting on the solutions. Non-linearities of product type are shown to admit this via exact paralinearization...... of homogeneous distributions, tensor products and halfspace extensions have been revised. Examples include the von Karman equation....

Numerical solution of fuzzy boundary value problems using Galerkin ...

Indian Academy of Sciences (India)

1 College of Mathematics and Statistics, Chongqing University, Chongqing 401331, China. 2 Department of ... exact solution of fuzzy first-order boundary value problems. (BVPs). ...... edge partial financial support by the Ministerio de Economıa.

Reconsidering the boundary conditions for a dynamic, transient mode I crack problem

KAUST Repository

Leise, Tanya

2008-11-01

A careful examination of a dynamic mode I crack problem leads to the conclusion that the commonly used boundary conditions do not always hold in the case of an applied crack face loading, so that a modification is required to satisfy the equations. In particular, a transient compressive stress wave travels along the crack faces, moving outward from the loading region on the crack face. This does not occur in the quasistatic or steady state problems, and is a special feature of the transient dynamic problem that is important during the time interval immediately following the application of crack face loading. We demonstrate why the usual boundary conditions lead to a prediction of crack face interpenetration, and then examine how to modify the boundary condition for a semi-infinite crack with a cohesive zone. Numerical simulations illustrate the resulting approach.

Description of internal flow problems by a boundary integral method with dipole panels

International Nuclear Information System (INIS)

Krieg, R.; Hailfinger, G.

1979-01-01

In reactor safety studies the failure of single components is postulated or sudden accident loadings are assumed and the consequences are investigated. Often as a first consequence highly transient three dimensional flow problems occur. In contrast to classical flow problems, in most of the above cases the fluid velocities are relatively small whereas the accelerations assume high values. As a consequence both, viscosity effects and dynamic pressures which are proportional to the square of the fluid velocities are usually negligible. For cases, where the excitation times are considerably longer than the times necessary for a wave to traverse characteristic regions of the fluid field, also the fluid compressibility is negligible. Under these conditions boundary integral methods are an appropriate tool to deal with the problem. Flow singularities are distributed over the fluid boundaries in such a way that pressure and velocity fields are obtained which satisfy the boundary conditions. In order to facilitate the numerical treatment the fluid boundaries are approximated by a finite number of panels with uniform singularity distributions on each of them. Consequently the pressure and velocity field of the given problem may be obtained by superposition of the corresponding fields due to these panels with their singularity intensities as unknown factors. Then satisfying the boundary conditions in so many boundary points as panels have been introduced, yields a system of linear equations which in general allows for a unique determination of the unknown intensities. (orig./RW)

An inverse boundary value problem for the Schroedinger operator with vector potentials in two dimensions

International Nuclear Information System (INIS)

Ziqi Sun

1993-01-01

During the past few years a considerable interest has been focused on the inverse boundary value problem for the Schroedinger operator with a scalar (electric) potential. The popularity gained by this subject seems to be due to its connection with the inverse scattering problem at fixed energy, the inverse conductivity problem and other important inverse problems. This paper deals with an inverse boundary value problem for the Schroedinger operator with vector (electric and magnetic) potentials. As in the case of the scalar potential, results of this study would have immediate consequences in the inverse scattering problem for magnetic field at fixed energy. On the other hand, inverse boundary value problems for elliptic operators are of independent interest. The study is partly devoted to the understanding of the inverse boundary value problem for a class of general elliptic operator of second order. Note that a self-adjoint elliptic operator of second order with Δ as its principal symbol can always be written as a Schroedinger operator with vector potentials

Numerical solution of system of boundary value problems using B-spline with free parameter

Science.gov (United States)

Gupta, Yogesh

2017-01-01

This paper deals with method of B-spline solution for a system of boundary value problems. The differential equations are useful in various fields of science and engineering. Some interesting real life problems involve more than one unknown function. These result in system of simultaneous differential equations. Such systems have been applied to many problems in mathematics, physics, engineering etc. In present paper, B-spline and B-spline with free parameter methods for the solution of a linear system of second-order boundary value problems are presented. The methods utilize the values of cubic B-spline and its derivatives at nodal points together with the equations of the given system and boundary conditions, ensuing into the linear matrix equation.

Numerical solutions of a three-point boundary value problem with an ...

African Journals Online (AJOL)

Numerical solutions of a three-point boundary value problem with an integral condition for a third-order partial differential equation by using Laplace transform method Solutions numeriques d'un probleme pour une classe d'equations differentielles d'ordr.

New formulations on the finite element method for boundary value problems with internal/external boundary layers

International Nuclear Information System (INIS)

Pereira, Luis Carlos Martins

1998-06-01

New Petrov-Galerkin formulations on the finite element methods for convection-diffusion problems with boundary layers are presented. Such formulations are based on a consistent new theory on discontinuous finite element methods. Existence and uniqueness of solutions for these problems in the new finite element spaces are demonstrated. Some numerical experiments shows how the new formulation operate and also their efficacy. (author)

METHOD OF GREEN FUNCTIONS IN MATHEMATICAL MODELLING FOR TWO-POINT BOUNDARY-VALUE PROBLEMS

Directory of Open Access Journals (Sweden)

E. V. Dikareva

2015-01-01

Full Text Available Summary. In many applied problems of control, optimization, system theory, theoretical and construction mechanics, for problems with strings and nods structures, oscillation theory, theory of elasticity and plasticity, mechanical problems connected with fracture dynamics and shock waves, the main instrument for study these problems is a theory of high order ordinary differential equations. This methodology is also applied for studying mathematical models in graph theory with different partitioning based on differential equations. Such equations are used for theoretical foundation of mathematical models but also for constructing numerical methods and computer algorithms. These models are studied with use of Green function method. In the paper first necessary theoretical information is included on Green function method for multi point boundary-value problems. The main equation is discussed, notions of multi-point boundary conditions, boundary functionals, degenerate and non-degenerate problems, fundamental matrix of solutions are introduced. In the main part the problem to study is formulated in terms of shocks and deformations in boundary conditions. After that the main results are formulated. In theorem 1 conditions for existence and uniqueness of solutions are proved. In theorem 2 conditions are proved for strict positivity and equal measureness for a pair of solutions. In theorem 3 existence and estimates are proved for the least eigenvalue, spectral properties and positivity of eigenfunctions. In theorem 4 the weighted positivity is proved for the Green function. Some possible applications are considered for a signal theory and transmutation operators.

Boundary conditions for free surface inlet and outlet problems

KAUST Repository

Taroni, M.; Breward, C. J. W.; Howell, P. D.; Oliver, J. M.

2012-01-01

We investigate and compare the boundary conditions that are to be applied to free-surface problems involving inlet and outlets of Newtonian fluid, typically found in coating processes. The flux of fluid is a priori known at an inlet, but unknown

On some boundary value problems in quantum statistical mechanics

International Nuclear Information System (INIS)

Angelescu, N.

1978-01-01

The following two topics of equilibrium quantum statistical mechanics are discussed in this thesis: (i) the independence of the thermodynamic limit of grand-canonical pressure on the boundary conditions; (ii) the magnetic properties of free quantum gases. Problem (i) is handled with a functional integration technique. Wiener-type conditional measures are constructed for a given domain and a general class of mixed conditions on its boundary, these measures are used to write down Feynman-Kac formulae for the kernels of exp(-βH), where H is the Hamiltonian of N interacting particles in the given domain. These measures share the property that they assign the same mass as the usual Wiener measure to any set of trajectories not intersecting the boundary. Local estimates on the kernels of exp(-βH) are derived, which imply independence of the pressure on the boundary conditions in the thermodynamic limit. Problem (ii) has a historical development: since Landau's work (1930), much discussion has been devoted to the influence of the finite size on the susceptibility. In finite volume, Dirichlet boundary conditions are imposed, on the ground that they ensure gauge invariance. The thermodynamic limit of the pressure is proved, using again functional integration. The functional measure is now complex but absolutely continuous with respect to Wiener measure, so the usual local estimates hold true. The controversy in the literature was concentrated on the commutativity of the operations of H-derivation and thermodynamic limit, so the existence of this limit for the zero-field susceptibility and its surface term are proved separately, demonstrating this commutativity. The proof relies on the following result of independent interest: the perturbation theory of self-adjoint trace-class semigroups is trace-class convergent and analytic. (author)

Boundary value problems on the half line in the theory of colloids

Directory of Open Access Journals (Sweden)

Ravi P. Agarwal

2002-01-01

Full Text Available We present existence results for some boundary value problems defined on infinite intervals. In particular our discussion includes a problem which arises in the theory of colloids.

About potential of double layer and boundary value problems for Laplace equation

International Nuclear Information System (INIS)

Aleshin, M.V.

1991-01-01

An integral operator raisen by a kernel of the double layer's potential is investigated. The kernel is defined on S (S - two-digit variety of C 2 class presented by a boundary of the finite domain in R 3 ). The operator is considered on C(S). Following results are received: the operator's spectrum belongs to [-1,1]; it's eigenvalues and eigenfunctions may be found by Kellog's method; knowledge of the operator's spectrum is enough to construct it's resolvent. These properties permit to point out the determined interation processes, solving boundary value problems for Laplace equation. One of such processes - solving of Roben problem - is generalized on electrostatic problems. 6 refs

Determination of free boundary problem of flow through porous media

International Nuclear Information System (INIS)

Tavares Junior, H.M.; Souza, A.J. de

1989-01-01

This paper deals with a free boundary problem of flow through porous media, which is solved by simplicial method conbined with mesh refinement. Variational method on fixed domain is utilized. (author)

A Boundary Element Solution to the Problem of Interacting AC Fields in Parallel Conductors

Directory of Open Access Journals (Sweden)

Einar M. Rønquist

1984-04-01

Full Text Available The ac fields in electrically insulated conductors will interact through the surrounding electromagnetic fields. The pertinent field equations reduce to the Helmholtz equation inside each conductor (interior problem, and to the Laplace equation outside the conductors (exterior problem. These equations are transformed to integral equations, with the magnetic vector potential and its normal derivative on the boundaries as unknowns. The integral equations are then approximated by sets of algebraic equations. The interior problem involves only unknowns on the boundary of each conductor, while the exterior problem couples unknowns from several conductors. The interior and the exterior problem are coupled through the field continuity conditions. The full set of equations is solved by standard Gaussian elimination. We also show how the total current and the dissipated power within each conductor can be expressed as boundary integrals. Finally, computational results for a sample problem are compared with a finite difference solution.

The homogeneous boundary value problem of the thick spherical shell

International Nuclear Information System (INIS)

Linder, F.

1975-01-01

With the aim to solve boundary value problems in the same manner as it is attained at thin shell theory (Superposition of Membrane solution to solution of boundary values), one has to search solutions of the equations of equilibrium of the three dimensional thick shell which produce tensions at the cut edge and are zero on the whole shell surface inside and outside. This problem was solved with the premissions of the linear theory of Elasticity. The gained solution is exact and contains the symmetric and non-symmetric behaviour and is described in relatively short analytical expressions for the deformations and tensions, after the problem of the coupled system had been solved. The static condition of the two surfaces (zero tension) leads to a homogeneous system of complex equations with the index of the Legendre spherical function as Eigenvalue. One symmetrical case is calculated numerically and is compared with the method of finite elements. This comparison results in good accordance. (Auth.)

A New Spectral Local Linearization Method for Nonlinear Boundary Layer Flow Problems

Directory of Open Access Journals (Sweden)

S. S. Motsa

2013-01-01

Full Text Available We propose a simple and efficient method for solving highly nonlinear systems of boundary layer flow problems with exponentially decaying profiles. The algorithm of the proposed method is based on an innovative idea of linearizing and decoupling the governing systems of equations and reducing them into a sequence of subsystems of differential equations which are solved using spectral collocation methods. The applicability of the proposed method, hereinafter referred to as the spectral local linearization method (SLLM, is tested on some well-known boundary layer flow equations. The numerical results presented in this investigation indicate that the proposed method, despite being easy to develop and numerically implement, is very robust in that it converges rapidly to yield accurate results and is more efficient in solving very large systems of nonlinear boundary value problems of the similarity variable boundary layer type. The accuracy and numerical stability of the SLLM can further be improved by using successive overrelaxation techniques.

Asymptotics for inhomogeneous Dirichlet initial-boundary value problem for the nonlinear Schrödinger equation

Energy Technology Data Exchange (ETDEWEB)

Kaikina, Elena I., E-mail: ekaikina@matmor.unam.mx [Centro de Ciencias Matemáticas, UNAM Campus Morelia, AP 61-3 (Xangari), Morelia CP 58089, Michoacán (Mexico)

2013-11-15

We consider the inhomogeneous Dirichlet initial-boundary value problem for the nonlinear Schrödinger equation, formulated on a half-line. We study traditionally important problems of the theory of nonlinear partial differential equations, such as global in time existence of solutions to the initial-boundary value problem and the asymptotic behavior of solutions for large time.

Asymptotics for inhomogeneous Dirichlet initial-boundary value problem for the nonlinear Schrödinger equation

International Nuclear Information System (INIS)

Kaikina, Elena I.

2013-01-01

We consider the inhomogeneous Dirichlet initial-boundary value problem for the nonlinear Schrödinger equation, formulated on a half-line. We study traditionally important problems of the theory of nonlinear partial differential equations, such as global in time existence of solutions to the initial-boundary value problem and the asymptotic behavior of solutions for large time

Order Reduction in High-Order Runge-Kutta Methods for Initial Boundary Value Problems

OpenAIRE

Rosales, Rodolfo Ruben; Seibold, Benjamin; Shirokoff, David; Zhou, Dong

2017-01-01

This paper studies the order reduction phenomenon for initial-boundary-value problems that occurs with many Runge-Kutta time-stepping schemes. First, a geometric explanation of the mechanics of the phenomenon is provided: the approximation error develops boundary layers, induced by a mismatch between the approximation error in the interior and at the boundaries. Second, an analysis of the modes of the numerical scheme is conducted, which explains under which circumstances boundary layers pers...

On nonseparated three-point boundary value problems for linear functional differential equations

Czech Academy of Sciences Publication Activity Database

Rontó, András; Rontó, M.

2011-01-01

Roč. 2011, - (2011), s. 326052 ISSN 1085-3375 Institutional research plan: CEZ:AV0Z10190503 Keywords : functional-differential equation * three-point boundary value problem * nonseparated boundary condition Subject RIV: BA - General Mathematics Impact factor: 1.318, year: 2011 http://www.hindawi.com/journals/ aaa /2011/326052/

Modified Differential Transform Method for Two Singular Boundary Values Problems

Directory of Open Access Journals (Sweden)

Yinwei Lin

2014-01-01

Full Text Available This paper deals with the two singular boundary values problems of second order. Two singular points are both boundary values points of the differential equation. The numerical solutions are developed by modified differential transform method (DTM for expanded point. Linear and nonlinear models are solved by this method to get more reliable and efficient numerical results. It can also solve ordinary differential equations where the traditional one fails. Besides, we give the convergence of this new method.

Algebraic structures in generalized Clifford analysis and applications to boundary value problems

Directory of Open Access Journals (Sweden)

José Játem

2015-12-01

Full Text Available The present article has a threefold purpose: First it is a survey of the algebraic structures of generalized Clifford-type algebras and shows the main results of the corresponding Clifford-type analysis and its application to boundary value problems known so far. Second it is aimed to implement algorithms to provide the fast and accurate computation of boundary value problems for inhomogeneous equations in the framework of the generalized Clifford analysis. Finally it is also aimed to encourage the development of a generalized discrete Clifford analysis.

Bibliography on moving boundary problems with key word index

International Nuclear Information System (INIS)

Wilson, D.G.; Solomon, A.D.; Trent, J.S.

1979-10-01

This bibliography concentrates mainly on time-dependent moving-boundary problems of heat and mass transfer. The bibliography is in two parts, a list of the references ordered by last name of the first author and a key word index to the titles. Few references from before 1965 are included

Bibliography on moving boundary problems with key word index

Energy Technology Data Exchange (ETDEWEB)

Wilson, D.G.; Solomon, A.D.; Trent, J.S.

1979-10-01

This bibliography concentrates mainly on time-dependent moving-boundary problems of heat and mass transfer. The bibliography is in two parts, a list of the references ordered by last name of the first author and a key word index to the titles. Few references from before 1965 are included. (RWR)

Asymptotic boundary value problems for evolution inclusions

Directory of Open Access Journals (Sweden)

Fürst Tomáš

2006-01-01

Full Text Available When solving boundary value problems on infinite intervals, it is possible to use continuation principles. Some of these principles take advantage of equipping the considered function spaces with topologies of uniform convergence on compact subintervals. This makes the representing solution operators compact (or condensing, but, on the other hand, spaces equipped with such topologies become more complicated. This paper shows interesting applications that use the strength of continuation principles and also presents a possible extension of such continuation principles to partial differential inclusions.

Asymptotic boundary value problems for evolution inclusions

Directory of Open Access Journals (Sweden)

Tomáš Fürst

2006-02-01

Full Text Available When solving boundary value problems on infinite intervals, it is possible to use continuation principles. Some of these principles take advantage of equipping the considered function spaces with topologies of uniform convergence on compact subintervals. This makes the representing solution operators compact (or condensing, but, on the other hand, spaces equipped with such topologies become more complicated. This paper shows interesting applications that use the strength of continuation principles and also presents a possible extension of such continuation principles to partial differential inclusions.

On Solutions of the Integrable Boundary Value Problem for KdV Equation on the Semi-Axis

International Nuclear Information System (INIS)

Ignatyev, M. Yu.

2013-01-01

This paper is concerned with the Korteweg–de Vries (KdV) equation on the semi-axis. The boundary value problem with inhomogeneous integrable boundary conditions is studied. We establish some characteristic properties of solutions of the problem. Also we construct a wide class of solutions of the problem using the inverse spectral method.

The extended-domain-eigenfunction method for solving elliptic boundary value problems with annular domains

Energy Technology Data Exchange (ETDEWEB)

Aarao, J; Bradshaw-Hajek, B H; Miklavcic, S J; Ward, D A, E-mail: Stan.Miklavcic@unisa.edu.a [School of Mathematics and Statistics, University of South Australia, Mawson Lakes, SA 5095 (Australia)

2010-05-07

Standard analytical solutions to elliptic boundary value problems on asymmetric domains are rarely, if ever, obtainable. In this paper, we propose a solution technique wherein we embed the original domain into one with simple boundaries where the classical eigenfunction solution approach can be used. The solution in the larger domain, when restricted to the original domain, is then the solution of the original boundary value problem. We call this the extended-domain-eigenfunction method. To illustrate the method's strength and scope, we apply it to Laplace's equation on an annular-like domain.

Numerical solution of large nonlinear boundary value problems by quadratic minimization techniques

International Nuclear Information System (INIS)

Glowinski, R.; Le Tallec, P.

1984-01-01

The objective of this paper is to describe the numerical treatment of large highly nonlinear two or three dimensional boundary value problems by quadratic minimization techniques. In all the different situations where these techniques were applied, the methodology remains the same and is organized as follows: 1) derive a variational formulation of the original boundary value problem, and approximate it by Galerkin methods; 2) transform this variational formulation into a quadratic minimization problem (least squares methods) or into a sequence of quadratic minimization problems (augmented lagrangian decomposition); 3) solve each quadratic minimization problem by a conjugate gradient method with preconditioning, the preconditioning matrix being sparse, positive definite, and fixed once for all in the iterative process. This paper will illustrate the methodology above on two different examples: the description of least squares solution methods and their application to the solution of the unsteady Navier-Stokes equations for incompressible viscous fluids; the description of augmented lagrangian decomposition techniques and their application to the solution of equilibrium problems in finite elasticity

The polynomial property of self-adjoint elliptic boundary-value problems and an algebraic description of their attributes

International Nuclear Information System (INIS)

Nazarov, S A

1999-01-01

We describe a wide class of boundary-value problems for which the application of elliptic theory can be reduced to elementary algebraic operations and which is characterized by the following polynomial property: the sesquilinear form corresponding to the problem degenerates only on some finite-dimensional linear space P of vector polynomials. Under this condition the boundary-value problem is elliptic, and its kernel and cokernel can be expressed in terms of P. For domains with piecewise-smooth boundary or infinite ends (conic, cylindrical, or periodic), we also present fragments of asymptotic formulae for the solutions, give specific versions of general conditional theorems on the Fredholm property (in particular, by modifying the ordinary weighted norms), and compute the index of the operator corresponding to the boundary-value problem. The polynomial property is also helpful for asymptotic analysis of boundary-value problems in thin domains and junctions of such domains. Namely, simple manipulations with P permit one to find the size of the system obtained by dimension reduction as well as the orders of the differential operators occurring in that system and provide complete information on the boundary layer structure. The results are illustrated by examples from elasticity and hydromechanics

Boundary-value problems in ODE

Science.gov (United States)

Tanriverdi, Tanfer

In this thesis we discuss two problems. The first problem is that of Fanno flow in a tube. In [10] the authors have discussed the mathematics of the Fanno model in much more detail than had been previously been done. The analysis in [10] indicates that the Fanno model becomes relevant, if t indicates the unscaled time and t=et , only when t is at least of order O(e- 1) . Indeed, two most important time scales are when t=O(e-1) and t=O(e- 2) . The authors, in the former case, set t=e- 1t1 (t1=t),x=e -11, and obtain the equation math> 62u6t 21- 62u 6x21=- 2u6 2u6t21 , ( 0.0.1) where u is the velocity of the gas, with p=1,6x1=0 (x1=0). One can follow the solution along the characteristic x1=t1 , and to match with the inviscid behaviour when t1-->0 , u=2+t1 (x1=t1). (0.0.2) In the region t=O(e2) , the authors set t=e2t2, x=e2x2,h= x2t2. For small e , the BC (0.0.02) now becomes u=t2 (x2=t 2), (0.0.3) so that (0.0.1) now has a similarity solution of the form u=t2g( h), u2=e- 1u, and (h2- 1)g'' +4hg'+2g=2g(g+hg' ),' =/ (0.0.4) with g(h)-->2 ash-->1- ,from(0.0.3) (0.0.5) g(h)-->∞ ash-->0- ,(fromthe pressure). ( 0.0.6) In a recent paper [11] the authors discuss the existence of a solution of (0.0.4)-(0.0.6) by using a two dimensional topological shooting method. We also discuss the existence of a solution of (0.0.4)-(0.0.6) by using a shooting method. We first turn the nonlinear ode (0.0.4) into an integral equation and then shoot from the singularity at ∞. The second problem arises when one considers eigenfunction expansions associated with second order ordinary differential equations, as Titchmarsh does in his book. One is concerned with the solutions of the equation - d2ydx2+ q(x)y=ly, (0.0.7) along with certain boundary conditions, where q(x)=-( n2- /)sech 2(x), n=n+/. The problem (0.0.7) has an application in the study of discrete reaction-diffusion equations. Our purpose in this problem is to look in some detail at the equation (0.0.7). We first use contour

The boundary element method applied to 3D magneto-electro-elastic dynamic problems

Science.gov (United States)

Igumnov, L. A.; Markov, I. P.; Kuznetsov, Iu A.

2017-11-01

Due to the coupling properties, the magneto-electro-elastic materials possess a wide number of applications. They exhibit general anisotropic behaviour. Three-dimensional transient analyses of magneto-electro-elastic solids can hardly be found in the literature. 3D direct boundary element formulation based on the weakly-singular boundary integral equations in Laplace domain is presented in this work for solving dynamic linear magneto-electro-elastic problems. Integral expressions of the three-dimensional fundamental solutions are employed. Spatial discretization is based on a collocation method with mixed boundary elements. Convolution quadrature method is used as a numerical inverse Laplace transform scheme to obtain time domain solutions. Numerical examples are provided to illustrate the capability of the proposed approach to treat highly dynamic problems.

On some examples of pollutant transport problems solved numerically using the boundary element method

Science.gov (United States)

Azis, Moh. Ivan; Kasbawati; Haddade, Amiruddin; Astuti Thamrin, Sri

2018-03-01

A boundary element method (BEM) is obtained for solving a boundary value problem of homogeneous anisotropic media governed by diffusion-convection equation. The application of the BEM is shown for two particular pollutant transport problems of Tello river and Unhas lake in Makassar Indonesia. For the two particular problems a variety of the coefficients of diffusion and the velocity components are taken. The results show that the solutions vary as the parameters change. And this suggests that one has to be careful in measuring or determining the values of the parameters.

Functional geometric method for solving free boundary problems for harmonic functions

Energy Technology Data Exchange (ETDEWEB)

Demidov, Aleksander S [M. V. Lomonosov Moscow State University, Moscow (Russian Federation)

2010-01-01

A survey is given of results and approaches for a broad spectrum of free boundary problems for harmonic functions of two variables. The main results are obtained by the functional geometric method. The core of these methods is an interrelated analysis of the functional and geometric characteristics of the problems under consideration and of the corresponding non-linear Riemann-Hilbert problems. An extensive list of open questions is presented. Bibliography: 124 titles.

L^p-continuity of solutions to parabolic free boundary problems

Directory of Open Access Journals (Sweden)

Abdeslem Lyaghfouri

2015-07-01

Full Text Available In this article, we consider a class of parabolic free boundary problems. We establish some properties of the solutions, including L^infinity-regularity in time and a monotonicity property, from which we deduce strong L^p-continuity in time.

Variational methods for boundary value problems for systems of elliptic equations

CERN Document Server

Lavrent'ev, M A

2012-01-01

Famous monograph by a distinguished mathematician presents an innovative approach to classical boundary value problems. The treatment employs the basic scheme first suggested by Hilbert and developed by Tonnelli. 1963 edition.

Parallel algorithms for boundary value problems

Science.gov (United States)

Lin, Avi

1991-01-01

A general approach to solve boundary value problems numerically in a parallel environment is discussed. The basic algorithm consists of two steps: the local step where all the P available processors work in parallel, and the global step where one processor solves a tridiagonal linear system of the order P. The main advantages of this approach are twofold. First, this suggested approach is very flexible, especially in the local step and thus the algorithm can be used with any number of processors and with any of the SIMD or MIMD machines. Secondly, the communication complexity is very small and thus can be used as easily with shared memory machines. Several examples for using this strategy are discussed.

Student Solutions Manual to Boundary Value Problems and Partial Differential Equations

CERN Document Server

Powers, David L

2005-01-01

This student solutions manual accompanies the text, Boundary Value Problems and Partial Differential Equations, 5e. The SSM is available in print via PDF or electronically, and provides the student with the detailed solutions of the odd-numbered problems contained throughout the book.Provides students with exercises that skillfully illustrate the techniques used in the text to solve science and engineering problemsNearly 900 exercises ranging in difficulty from basic drills to advanced problem-solving exercisesMany exercises based on current engineering applications

Boundary value problem for Caputo-Hadamard fractional differential equations

Directory of Open Access Journals (Sweden)

Yacine Arioua

2017-09-01

Full Text Available The aim of this work is to study the existence and uniqueness solutions for boundary value problem of nonlinear fractional differential equations with Caputo-Hadamard derivative in bounded domain. We used the standard and Krasnoselskii's fixed point theorems. Some new results of existence and uniqueness solutions for Caputo-Hadamard fractional equations are obtained.

Positive solutions for a nonlocal boundary-value problem with vector-valued response

Directory of Open Access Journals (Sweden)

Andrzej Nowakowski

2002-05-01

Full Text Available Using variational methods, we study the existence of positive solutions for a nonlocal boundary-value problem with vector-valued response. We develop duality and variational principles for this problem and present a numerical version which enables the approximation of solutions and gives a measure of a duality gap between primal and dual functional for approximate solutions for this problem.

POSITIVE SOLUTIONS OF A NONLINEAR THREE-POINT EIGENVALUE PROBLEM WITH INTEGRAL BOUNDARY CONDITIONS

Directory of Open Access Journals (Sweden)

FAOUZI HADDOUCHI

2015-11-01

Full Text Available In this paper, we study the existence of positive solutions of a three-point integral boundary value problem (BVP for the following second-order differential equation u''(t + \\lambda a(tf(u(t = 0; 0 0 is a parameter, 0 <\\eta < 1, 0 <\\alpha < 1/{\\eta}. . By using the properties of the Green's function and Krasnoselskii's fixed point theorem on cones, the eigenvalue intervals of the nonlinear boundary value problem are considered, some sufficient conditions for the existence of at least one positive solutions are established.

Solving inverse two-point boundary value problems using collage coding

Science.gov (United States)

Kunze, H.; Murdock, S.

2006-08-01

The method of collage coding, with its roots in fractal imaging, is the central tool in a recently established rigorous framework for solving inverse initial value problems for ordinary differential equations (Kunze and Vrscay 1999 Inverse Problems 15 745-70). We extend these ideas to solve the following inverse problem: given a function u(x) on [A, B] (which may be the interpolation of data points), determine a two-point boundary value problem on [A, B] which admits u(x) as a solution as closely as desired. The solution of such inverse problems may be useful in parameter estimation or determination of potential functional forms of the underlying differential equation. We discuss ways to improve results, including the development of a partitioning scheme. Several examples are considered.

The Dirichlet problem with L2-boundary data for elliptic linear equations

CERN Document Server

Chabrowski, Jan

1991-01-01

The Dirichlet problem has a very long history in mathematics and its importance in partial differential equations, harmonic analysis, potential theory and the applied sciences is well-known. In the last decade the Dirichlet problem with L2-boundary data has attracted the attention of several mathematicians. The significant features of this recent research are the use of weighted Sobolev spaces, existence results for elliptic equations under very weak regularity assumptions on coefficients, energy estimates involving L2-norm of a boundary data and the construction of a space larger than the usual Sobolev space W1,2 such that every L2-function on the boundary of a given set is the trace of a suitable element of this space. The book gives a concise account of main aspects of these recent developments and is intended for researchers and graduate students. Some basic knowledge of Sobolev spaces and measure theory is required.

Existence of 2m-1 Positive Solutions for Sturm-Liouville Boundary Value Problems with Linear Functional Boundary Conditions on the Half-Line

Directory of Open Access Journals (Sweden)

Yanmei Sun

2012-01-01

Full Text Available By using the Leggett-Williams fixed theorem, we establish the existence of multiple positive solutions for second-order nonhomogeneous Sturm-Liouville boundary value problems with linear functional boundary conditions. One explicit example with singularity is presented to demonstrate the application of our main results.

Problem of the Moving Boundary in Continuous Casting Solved by The Analytic-Numerical Method

Directory of Open Access Journals (Sweden)

Grzymkowski R.

2013-03-01

Full Text Available Mathematical modeling of thermal processes combined with the reversible phase transitions of type: solid phase - liquid phase leads to formulation of the parabolic or elliptic moving boundary problem. Solution of such defined problem requires, most often, to use some sophisticated numerical techniques and far advanced mathematical tools. The paper presents an analytic-numerical method, especially attractive from the engineer’s point of view, applied for finding the approximate solutions of the selected class of problems which can be reduced to the one-phase solidification problem of a plate with the unknown a priori, varying in time boundary of the region in which the solution is sought. Proposed method is based on the known formalism of initial expansion of a sought function, describing the field of temperature, into the power series, some coefficients of which are determined with the aid of boundary conditions, and on the approximation of a function defining the freezing front location with the broken line, parameters of which are determined numerically. The method represents a combination of the analytical and numerical techniques and seems to be an effective and relatively easy in using tool for solving problems of considered kind.

Problem of the Moving Boundary in Continuous Casting Solved by the Analytic-Numerical Method

Directory of Open Access Journals (Sweden)

R. Grzymkowski

2013-01-01

Full Text Available Mathematical modeling of thermal processes combined with the reversible phase transitions of type: solid phase – liquid phase leads to formulation of the parabolic or elliptic moving boundary problem. Solution of such defined problem requires, most often, to use some sophisticated numerical techniques and far advanced mathematical tools. The paper presents an analytic-numerical method, especially attractive from the engineer’s point of view, applied for finding the approximate solutions of the selected class of problems which can be reduced to the one-phase solidification problem of a plate with the unknown a priori, varying in time boundary of the region in which the solution is sought. Proposed method is based on the known formalism of initial expansion of a sought function, describing the field of temperature, into the power series, some coefficients of which are determined with the aid of boundary conditions, and on the approximation of a function defining the freezing front location with the broken line, parameters of which are determined numerically. The method represents a combination of the analytical and numerical techniques and seems to be an effective and relatively easy in using tool for solving problems of considered kind.

Construction of local boundary conditions for an eigenvalue problem using micro-local analysis: application to optical waveguide problems

International Nuclear Information System (INIS)

Barucq, Helene; Bekkey, Chokri; Djellouli, Rabia

2004-01-01

We present a general procedure based on the pseudo-differential calculus for deriving artificial boundary conditions for an eigenvalue problem that characterizes the propagation of guided modes in optical waveguides. This new approach allows the construction of local conditions that (a) are independent of the frequency regime, (b) preserve the sparsity pattern of the finite element discretization, and (c) are applicable to arbitrarily shaped convex artificial boundaries. The last feature has the potential for reducing the size of the computational domain. Numerical results are presented to highlight the potential of conditions of order 1/2 and 1, for improving significantly the computational efficiency of finite element methods for the solution of optical waveguide problems

On the radiative transfer problem in a spherical medium subject to Fresnel's reflective boundary conditions

International Nuclear Information System (INIS)

Mohammed, M.H.H.

2012-01-01

Radiation transfer problem for anisotropic scattering in a spherical homogeneous, turbid medium with angular dependent (specular) and diffuse reflecting boundary is considered. The angular dependent reflectivity of the boundary is considered as Fresnel's reflection probability function. The solution of the problem containing an energy source in a medium of specular and diffuse reflecting boundaries is given in terms of the solution of the source-free problem. The source-free problem for anisotropic scattering through a homogeneous solid sphere and two concentric spheres is solved by using the Pomraning- Eddington approximation method. This method transform the integro-differential equation into two differential equations for the radiance g (x) and net flux q (x) which has an analytical solution in terms of the modified Bessel function. Two different weight functions are used to verify the boundary conditions and so, find the solution constants. The partial heat fluxes at the boundaries of a solid sphere and spherical shell of transparent and reflecting boundaries are calculated. The media are taken with or without internal black-body radiation. The calculations are carried out for various values of refractive index and different radii. The results are compared with those of the Galerkin technique

The numerical solution of boundary value problems over an infinite domain

International Nuclear Information System (INIS)

Shepherd, M.; Skinner, R.

1976-01-01

A method is presented for the numerical solution of boundary value problems over infinite domains. An example that illustrates also the strength and accuracy of a numerical procedure for calculating Green's functions is described in detail

Multiple positive solutions for second order impulsive boundary value problems in Banach spaces

Directory of Open Access Journals (Sweden)

Zhi-Wei Lv

2010-06-01

Full Text Available By means of the fixed point index theory of strict set contraction operators, we establish new existence theorems on multiple positive solutions to a boundary value problem for second-order impulsive integro-differential equations with integral boundary conditions in a Banach space. Moreover, an application is given to illustrate the main result.

On the Approximate Controllability of Some Semilinear Parabolic Boundary-Value Problems

International Nuclear Information System (INIS)

Diaz, J. I.; Henry, J.; Ramos, A. M.

1998-01-01

We prove the approximate controllability of several nonlinear parabolic boundary-value problems by means of two different methods: the first one can be called a Cancellation method and the second one uses the Kakutani fixed-point theorem

A free-boundary value problem related to auto ignition of ...

African Journals Online (AJOL)

We examine a free boundary value problem related to auto ignition of combustible fluid in insulation materials. The criteria for the existence of similarity solution of the model equations are established. The conditions for the existence of unique solution are also stated. The numerical results which show the influence of ...

Moving-boundary problems for the time-fractional diffusion equation

Directory of Open Access Journals (Sweden)

Sabrina D. Roscani

2017-02-01

Full Text Available We consider a one-dimensional moving-boundary problem for the time-fractional diffusion equation. The time-fractional derivative of order $\\alpha\\in (0,1$ is taken in the sense of Caputo. We study the asymptotic behaivor, as t tends to infinity, of a general solution by using a fractional weak maximum principle. Also, we give some particular exact solutions in terms of Wright functions.

Use of Green's functions in the numerical solution of two-point boundary value problems

Science.gov (United States)

Gallaher, L. J.; Perlin, I. E.

1974-01-01

This study investigates the use of Green's functions in the numerical solution of the two-point boundary value problem. The first part deals with the role of the Green's function in solving both linear and nonlinear second order ordinary differential equations with boundary conditions and systems of such equations. The second part describes procedures for numerical construction of Green's functions and considers briefly the conditions for their existence. Finally, there is a description of some numerical experiments using nonlinear problems for which the known existence, uniqueness or convergence theorems do not apply. Examples here include some problems in finding rendezvous orbits of the restricted three body system.

How to approximate the heat equation with Neumann boundary conditions by nonlocal diffusion problems

OpenAIRE

Cortazar, C.; Elgueta, M.; Rossi, J. D.; Wolanski, N.

2006-01-01

We present a model for nonlocal diffusion with Neumann boundary conditions in a bounded smooth domain prescribing the flux through the boundary. We study the limit of this family of nonlocal diffusion operators when a rescaling parameter related to the kernel of the nonlocal operator goes to zero. We prove that the solutions of this family of problems converge to a solution of the heat equation with Neumann boundary conditions.

Existence and uniqueness for a two-point interface boundary value problem

Directory of Open Access Journals (Sweden)

Rakhim Aitbayev

2013-10-01

Full Text Available We obtain sufficient conditions, easily verifiable, for the existence and uniqueness of piecewise smooth solutions of a linear two-point boundary-value problem with general interface conditions. The coefficients of the differential equation may have jump discontinuities at the interface point. As an example, the conditions obtained are applied to a problem with typical interface such as perfect contact, non-perfect contact, and flux jump conditions.

Robust numerical methods for boundary-layer equations for a model problem of flow over a symmetric curved surface

NARCIS (Netherlands)

A.R. Ansari; B. Hossain; B. Koren (Barry); G.I. Shishkin (Gregori)

2007-01-01

textabstractWe investigate the model problem of flow of a viscous incompressible fluid past a symmetric curved surface when the flow is parallel to its axis. This problem is known to exhibit boundary layers. Also the problem does not have solutions in closed form, it is modelled by boundary-layer

Boundary-value problems for first and second order functional differential inclusions

Directory of Open Access Journals (Sweden)

Shihuang Hong

2003-03-01

Full Text Available This paper presents sufficient conditions for the existence of solutions to boundary-value problems of first and second order multi-valued differential equations in Banach spaces. Our results obtained using fixed point theorems, and lead to new existence principles.

Generalised functions method in the boundary value problems of elastodynamics by stationary running loads

International Nuclear Information System (INIS)

Alexeyeva, L.A.

2001-01-01

Investigation of diffraction processes of seismic waves on underground tunnels and pipelines with use of mathematical methods is related to solving boundary value problems (BVP) for hyperbolic system of differential equations in domains with cylindrical cavities when seismic disturbances propagate along boundaries with subsonic or transonic speeds. Also such classes of problems appear when it's necessary to study the behavior of underground constructions and Stress-strain State of environment. But in this case the velocities of running loads are less than velocities of wave propagation in surrounding medium. At present similar problems were solved only for constructions of circular cylindrical form with use of methods of full and not full dividing of variables. For cylindrical constructions of complex cross section strong mathematical theories for solving these problems were absent.(author)

Adaptive algorithm for solution of early exercise boundary problem of American put option implemented in Mathematica

Directory of Open Access Journals (Sweden)

Lukáš Ladislav

2017-01-01

Full Text Available The paper is focused on American option pricing problem. Assuming non-dividend paying American put option leads to two disjunctive regions, a continuation one and a stopping one, which are separated by an early exercise boundary. We present variational formulation of American option problem with special attention to early exercise action effect. Next, we discuss financially motivated additive decomposition of American option price into a European option price and another part due to the extra premium required by early exercising the option contract. As the optimal exercise boundary is a free boundary, its determination is coupled with the determination of the option price. Therefore, a closed-form expression of the free boundary is not attainable in general. We discuss in detail a derivation of an asymptotic expression of the early exercise boundary. Finally, we present some numerical results of determination of free boundary based upon this approach. All computations are performed by sw Mathematica, and suitable numerical procedure is discussed in detail, as well.

High order methods for incompressible fluid flow: Application to moving boundary problems

Energy Technology Data Exchange (ETDEWEB)

Bjoentegaard, Tormod

2008-04-15

Fluid flows with moving boundaries are encountered in a large number of real life situations, with two such types being fluid-structure interaction and free-surface flows. Fluid-structure phenomena are for instance apparent in many hydrodynamic applications; wave effects on offshore structures, sloshing and fluid induced vibrations, and aeroelasticity; flutter and dynamic response. Free-surface flows can be considered as a special case of a fluid-fluid interaction where one of the fluids are practically inviscid, such as air. This type of flows arise in many disciplines such as marine hydrodynamics, chemical engineering, material processing, and geophysics. The driving forces for free-surface flows may be of large scale such as gravity or inertial forces, or forces due to surface tension which operate on a much smaller scale. Free-surface flows with surface tension as a driving mechanism include the flow of bubbles and droplets, and the evolution of capillary waves. In this work we consider incompressible fluid flow, which are governed by the incompressible Navier-Stokes equations. There are several challenges when simulating moving boundary problems numerically, and these include - Spatial discretization - Temporal discretization - Imposition of boundary conditions - Solution strategy for the linear equations. These are some of the issues which will be addressed in this introduction. We will first formulate the problem in the arbitrary Lagrangian-Eulerian framework, and introduce the weak formulation of the problem. Next, we discuss the spatial and temporal discretization before we move to the imposition of surface tension boundary conditions. In the final section we discuss the solution of the resulting linear system of equations. (Author). refs., figs., tabs

An initial boundary value problem for modeling a piezoelectric dipolar body

Science.gov (United States)

Marin, Marin; Öchsner, Andreas

2018-03-01

This study deals with the first initial boundary value problem in elasticity of piezoelectric dipolar bodies. We consider the most general case of an anisotropic and inhomogeneous elastic body having a dipolar structure. For two different types of restrictions imposed on the problem data, we prove two results regarding the uniqueness of solution, by using a different but accessible method. Then, the mixed problem is transformed in a temporally evolutionary equation on a Hilbert space, conveniently constructed based on the problem data. With the help of a known result from the theory of semigroups of operators, the existence and uniqueness of the weak solution for this equation are proved.

Dual reciprocity boundary element analysis for the laminar forced heat convection problem in concentric annulus

International Nuclear Information System (INIS)

Choi, Chang Yong

1999-01-01

This paper presents a study of the Dual Reciprocity Boundary Element Method (DRBEM) for the laminar heat convection problem in a concentric annulus with constant heat flux boundary condition. DRBEM is one of the most successful technique used to transform the domain integrals arising from the nonhomogeneous term of the poisson equation into equivalent boundary only integrals. This recently developed and highly efficient numerical method is tested for the solution accuracy of the fluid flow and heat transfer study in a concentric annulus. Since their exact solutions are available, DRBEM solutions are verified with different number of boundary element discretization and internal points. The results obtained in this study are discussed with the relative error percentage of velocity and temperature solutions, and potential applicability of the method for the more complicated heat convection problems with arbitrary duct geometries

A Hamiltonian-based derivation of Scaled Boundary Finite Element Method for elasticity problems

International Nuclear Information System (INIS)

Hu Zhiqiang; Lin Gao; Wang Yi; Liu Jun

2010-01-01

The Scaled Boundary Finite Method (SBFEM) is a semi-analytical solution approach for solving partial differential equation. For problem in elasticity, the governing equations can be obtained by mechanically based formulation, Scaled-boundary-transformation-based formulation and principle of virtual work. The governing equations are described in the frame of Lagrange system and the unknowns are displacements. But in the solution procedure, the auxiliary variables are introduced and the equations are solved in the state space. Based on the observation that the duality system to solve elastic problem proposed by W.X. Zhong is similar to the above solution approach, the discretization of the SBFEM and the duality system are combined to derive the governing equations in the Hamilton system by introducing the dual variables in this paper. The Precise Integration Method (PIM) used in Duality system is also an efficient method for the solution of the governing equations of SBFEM in displacement and boundary stiffness matrix especially for the case which results some numerical difficulties in the usually uses the eigenvalue method. Numerical examples are used to demonstrate the validity and effectiveness of the PIM for solution of boundary static stiffness.

A New Numerical Algorithm for Two-Point Boundary Value Problems

OpenAIRE

Guo, Lihua; Wu, Boying; Zhang, Dazhi

2014-01-01

We present a new numerical algorithm for two-point boundary value problems. We first present the exact solution in the form of series and then prove that the n-term numerical solution converges uniformly to the exact solution. Furthermore, we establish the numerical stability and error analysis. The numerical results show the effectiveness of the proposed algorithm.

Numerical Analysis of Forth-Order Boundary Value Problems in Fluid Mechanics and Mathematics

DEFF Research Database (Denmark)

Hosseinzadeh, E.; Barari, Amin; Fouladi, F.

2011-01-01

In this paper He's variational iteration method is used to solve some examples of linear and non-linear forth-order boundary value problems. The first problem compared with homotopy analysis method solution and the other ones with the exact solution. The results show the high accuracy and speed o...

Numerical analysis of fourth-order boundary value problems in fluid mechanics and mathematics

DEFF Research Database (Denmark)

Hosseinzadeh, Elham; Barari, Amin; Fouladi, Fama

2010-01-01

In this paper He's variational iteration method is used to solve some examples of linear and non-linear forth-order boundary value problems. The first problem compared with homotopy analysis method solution and the other ones with the exact solution. The results show the high accuracy and speed o...

Mathematical and numerical study of nonlinear boundary problems related to plasma physics

International Nuclear Information System (INIS)

Sermange, M.

1982-06-01

After the study of some equations based on the Hodgkin-Huxley model, the work presented here is concerned with nonlinear boundary problems in MHD. They are gathered in two subjects: equilibrium equations and stability equations. The axisymmetric MHD equilibrium equations with free boundary have been studied by different authors, particularly the existence, regularity, unicity and non-unicity. Here, bifurcation, convergence of calculation methods existence of solutions in a discontinuous frame are studied. MHD stability can be determined by the principle of Bernstein et al; the mathematical work concerned here bears on the equivalence, in the case of two-dimensional or axisymmetric stability, between this model and a scalar eigenvalue problem which is introduced. At last, modules for computing MHD equilibrium for the simulation of plasma confinement in a tokamak are described [fr

Abraham Pais Prize Lecture: Shifting Problems and Boundaries in the History of Modern Physics

Science.gov (United States)

Nye, Mary-Jo

A long established category of study in the history of science is the ``history of physical sciences.'' It is a category that immediately begs the question of disciplinary boundaries for the problems and subjects addressed in historical inquiry. As a historian of the physical sciences, I often have puzzled over disciplinary boundaries and the means used to create or justify them. Scientists most often have been professionally identified with specific institutionalized fields since the late 19th century, but the questions they ask and the problems they solve are not neatly carved up by disciplinary perimeters. Like institutional departments or professorships, the Nobel Prizes in the 20th century often have delineated the scope of ``Physics'' or ``Chemistry'' (and ``Physiology or Medicine''), but the Prizes do not reflect disciplinary rigidity, despite some standard core subjects. In this paper I examine trends in Nobel Prize awards that indicate shifts in problem solving and in boundaries in twentieth century physics, tying those developments to changing themes in the history of physics and physical science in recent decades.

Existence of Three Positive Solutions to Some p-Laplacian Boundary Value Problems

Directory of Open Access Journals (Sweden)

Moulay Rchid Sidi Ammi

2013-01-01

Full Text Available We obtain, by using the Leggett-Williams fixed point theorem, sufficient conditions that ensure the existence of at least three positive solutions to some p-Laplacian boundary value problems on time scales.

Robust Wavelet Estimation to Eliminate Simultaneously the Effects of Boundary Problems, Outliers, and Correlated Noise

Directory of Open Access Journals (Sweden)

Alsaidi M. Altaher

2012-01-01

Full Text Available Classical wavelet thresholding methods suffer from boundary problems caused by the application of the wavelet transformations to a finite signal. As a result, large bias at the edges and artificial wiggles occur when the classical boundary assumptions are not satisfied. Although polynomial wavelet regression and local polynomial wavelet regression effectively reduce the risk of this problem, the estimates from these two methods can be easily affected by the presence of correlated noise and outliers, giving inaccurate estimates. This paper introduces two robust methods in which the effects of boundary problems, outliers, and correlated noise are simultaneously taken into account. The proposed methods combine thresholding estimator with either a local polynomial model or a polynomial model using the generalized least squares method instead of the ordinary one. A primary step that involves removing the outlying observations through a statistical function is considered as well. The practical performance of the proposed methods has been evaluated through simulation experiments and real data examples. The results are strong evidence that the proposed method is extremely effective in terms of correcting the boundary bias and eliminating the effects of outliers and correlated noise.

Integral methods of solving boundary-value problems of nonstationary heat conduction and their comparative analysis

Science.gov (United States)

Kot, V. A.

2017-11-01

The modern state of approximate integral methods used in applications, where the processes of heat conduction and heat and mass transfer are of first importance, is considered. Integral methods have found a wide utility in different fields of knowledge: problems of heat conduction with different heat-exchange conditions, simulation of thermal protection, Stefantype problems, microwave heating of a substance, problems on a boundary layer, simulation of a fluid flow in a channel, thermal explosion, laser and plasma treatment of materials, simulation of the formation and melting of ice, inverse heat problems, temperature and thermal definition of nanoparticles and nanoliquids, and others. Moreover, polynomial solutions are of interest because the determination of a temperature (concentration) field is an intermediate stage in the mathematical description of any other process. The following main methods were investigated on the basis of the error norms: the Tsoi and Postol’nik methods, the method of integral relations, the Gudman integral method of heat balance, the improved Volkov integral method, the matched integral method, the modified Hristov method, the Mayer integral method, the Kudinov method of additional boundary conditions, the Fedorov boundary method, the method of weighted temperature function, the integral method of boundary characteristics. It was established that the two last-mentioned methods are characterized by high convergence and frequently give solutions whose accuracy is not worse that the accuracy of numerical solutions.

An Existence Principle for Nonlocal Difference Boundary Value Problems with φ-Laplacian and Its Application to Singular Problems

Directory of Open Access Journals (Sweden)

Svatoslav Stanêk

2008-03-01

Full Text Available The paper presents an existence principle for solving a large class of nonlocal regular discrete boundary value problems with the φ-Laplacian. Applications of the existence principle to singular discrete problems are given.

Existence and uniqueness of entropy solution to initial boundary value problem for the inviscid Burgers equation

CERN Document Server

Zhu, C

2003-01-01

This paper is concerned with the existence and uniqueness of the entropy solution to the initial boundary value problem for the inviscid Burgers equation. To apply the method of vanishing viscosity to study the existence of the entropy solution, we first introduce the initial boundary value problem for the viscous Burgers equation, and as in Evans (1998 Partial Differential Equations (Providence, RI: American Mathematical Society) and Hopf (1950 Commun. Pure Appl. Math. 3 201-30), give the formula of the corresponding viscosity solutions by Hopf-Cole transformation. Secondly, we prove the convergence of the viscosity solution sequences and verify that the limiting function is an entropy solution. Finally, we give an example to show how our main result can be applied to solve the initial boundary value problem for the Burgers equation.

Existence and uniqueness of entropy solution to initial boundary value problem for the inviscid Burgers equation

International Nuclear Information System (INIS)

Zhu, Changjiang; Duan, Renjun

2003-01-01

This paper is concerned with the existence and uniqueness of the entropy solution to the initial boundary value problem for the inviscid Burgers equation. To apply the method of vanishing viscosity to study the existence of the entropy solution, we first introduce the initial boundary value problem for the viscous Burgers equation, and as in Evans (1998 Partial Differential Equations (Providence, RI: American Mathematical Society) and Hopf (1950 Commun. Pure Appl. Math. 3 201-30), give the formula of the corresponding viscosity solutions by Hopf-Cole transformation. Secondly, we prove the convergence of the viscosity solution sequences and verify that the limiting function is an entropy solution. Finally, we give an example to show how our main result can be applied to solve the initial boundary value problem for the Burgers equation

Multi-point boundary value problems for linear functional-differential equations

Czech Academy of Sciences Publication Activity Database

Domoshnitsky, A.; Hakl, Robert; Půža, Bedřich

2017-01-01

Roč. 24, č. 2 (2017), s. 193-206 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : boundary value problems * linear functional-differential equations * functional-differential inequalities Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.290, year: 2016 https://www.degruyter.com/view/j/gmj.2017.24.issue-2/gmj-2016-0076/gmj-2016-0076. xml

Multi-point boundary value problems for linear functional-differential equations

Czech Academy of Sciences Publication Activity Database

Domoshnitsky, A.; Hakl, Robert; Půža, Bedřich

2017-01-01

Roč. 24, č. 2 (2017), s. 193-206 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : boundary value problems * linear functional- differential equations * functional- differential inequalities Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.290, year: 2016 https://www.degruyter.com/view/j/gmj.2017.24.issue-2/gmj-2016-0076/gmj-2016-0076.xml

Infinitely many solutions for a fourth-order boundary-value problem

Directory of Open Access Journals (Sweden)

Seyyed Mohsen Khalkhali

2012-09-01

Full Text Available In this article we consider the existence of infinitely many solutions to the fourth-order boundary-value problem $$displaylines{ u^{iv}+alpha u''+eta(x u=lambda f(x,u+h(u,quad xin]0,1[cr u(0=u(1=0,cr u''(0=u''(1=0,. }$$ Our approach is based on variational methods and critical point theory.

Initial boundary value problems of nonlinear wave equations in an exterior domain

International Nuclear Information System (INIS)

Chen Yunmei.

1987-06-01

In this paper, we investigate the existence and uniqueness of the global solutions to the initial boundary value problems of nonlinear wave equations in an exterior domain. When the space dimension n >= 3, the unique global solution of the above problem is obtained for small initial data, even if the nonlinear term is fully nonlinear and contains the unknown function itself. (author). 10 refs

Asymptotics of linear initial boundary value problems with periodic boundary data on the half-line and finite intervals

KAUST Repository

Dujardin, G. M.

2009-08-12

This paper deals with the asymptotic behaviour of the solutions of linear initial boundary value problems with constant coefficients on the half-line and on finite intervals. We assume that the boundary data are periodic in time and we investigate whether the solution becomes time-periodic after sufficiently long time. Using Fokas\\' transformation method, we show that, for the linear Schrödinger equation, the linear heat equation and the linearized KdV equation on the half-line, the solutions indeed become periodic for large time. However, for the same linear Schrödinger equation on a finite interval, we show that the solution, in general, is not asymptotically periodic; actually, the asymptotic behaviour of the solution depends on the commensurability of the time period T of the boundary data with the square of the length of the interval over. © 2009 The Royal Society.

The Laplace equation boundary value problems on bounded and unbounded Lipschitz domains

CERN Document Server

Medková, Dagmar

2018-01-01

This book is devoted to boundary value problems of the Laplace equation on bounded and unbounded Lipschitz domains. It studies the Dirichlet problem, the Neumann problem, the Robin problem, the derivative oblique problem, the transmission problem, the skip problem and mixed problems. It also examines different solutions - classical, in Sobolev spaces, in Besov spaces, in homogeneous Sobolev spaces and in the sense of non-tangential limit. It also explains relations between different solutions. The book has been written in a way that makes it as readable as possible for a wide mathematical audience, and includes all the fundamental definitions and propositions from other fields of mathematics. This book is of interest to research students, as well as experts in partial differential equations and numerical analysis.

An analytical approximation scheme to two-point boundary value problems of ordinary differential equations

International Nuclear Information System (INIS)

Boisseau, Bruno; Forgacs, Peter; Giacomini, Hector

2007-01-01

A new (algebraic) approximation scheme to find global solutions of two-point boundary value problems of ordinary differential equations (ODEs) is presented. The method is applicable for both linear and nonlinear (coupled) ODEs whose solutions are analytic near one of the boundary points. It is based on replacing the original ODEs by a sequence of auxiliary first-order polynomial ODEs with constant coefficients. The coefficients in the auxiliary ODEs are uniquely determined from the local behaviour of the solution in the neighbourhood of one of the boundary points. The problem of obtaining the parameters of the global (connecting) solutions, analytic at one of the boundary points, reduces to find the appropriate zeros of algebraic equations. The power of the method is illustrated by computing the approximate values of the 'connecting parameters' for a number of nonlinear ODEs arising in various problems in field theory. We treat in particular the static and rotationally symmetric global vortex, the skyrmion, the Abrikosov-Nielsen-Olesen vortex, as well as the 't Hooft-Polyakov magnetic monopole. The total energy of the skyrmion and of the monopole is also computed by the new method. We also consider some ODEs coming from the exact renormalization group. The ground-state energy level of the anharmonic oscillator is also computed for arbitrary coupling strengths with good precision. (fast track communication)

Collisional plasma transport: two-dimensional scalar formulation of the initial boundary value problem and quasi one-dimensional models

International Nuclear Information System (INIS)

Mugge, J.W.

1979-10-01

The collisional plasma transport problem is formulated as an initial boundary value problem for general characteristic boundary conditions. Starting from the full set of hydrodynamic and electrodynamic equations an expansion in the electron-ion mass ratio together with a multiple timescale method yields simplified equations on each timescale. On timescales where many collisions have taken place for the simplified equations the initial boundary value problem is formulated. Through the introduction of potentials a two-dimensional scalar formulation in terms of quasi-linear integro-differential equations of second order for a domain consisting of plasma and vacuum sub-domains is obtained. (Auth.)

On a free boundary problem for a strongly degenerate quasilinear parabolic equation with an application to a model of pressure filtration

Energy Technology Data Exchange (ETDEWEB)

Buerger, R.; Frid, H.; Karlsen, K.H.

2002-07-01

We consider a free boundary problem of a quasilinear strongly degenerate parabolic equation arising from a model of pressure filtration of flocculated suspensions. We provide definitions of generalized solutions of the free boundary problem in the framework of L2 divergence-measure fields. The formulation of boundary conditions is based on a Gauss-Green theorem for divergence-measure fields on bounded domains with Lipschitz deformable boundaries and avoids referring to traces of the solution. This allows to consider generalized solutions from a larger class than BV. Thus it is not necessary to derive the usual uniform estimates on spatial and time derivatives of the solutions of the corresponding regularized problem requires in the BV approach. We first prove existence and uniqueness of the solution of the regularized parabolic free boundary problem and then apply the vanishing viscosity method to prove existence of a generalized solution to the degenerate free boundary problem. (author)

Solution matching for a three-point boundary-value problem on atime scale

Directory of Open Access Journals (Sweden)

Martin Eggensperger

2004-07-01

Full Text Available Let $mathbb{T}$ be a time scale such that $t_1, t_2, t_3 in mathbb{T}$. We show the existence of a unique solution for the three-point boundary value problem $$displaylines{ y^{DeltaDeltaDelta}(t = f(t, y(t, y^Delta(t, y^{DeltaDelta}(t, quad t in [t_1, t_3] cap mathbb{T},cr y(t_1 = y_1, quad y(t_2 = y_2, quad y(t_3 = y_3,. }$$ We do this by matching a solution to the first equation satisfying a two-point boundary conditions on $[t_1, t_2] cap mathbb{T}$ with a solution satisfying a two-point boundary conditions on $[t_2, t_3] cap mathbb{T}$.

On a boundary layer problem related to the gas flow in shales

KAUST Repository

Barenblatt, G. I.

2013-01-16

The development of gas deposits in shales has become a significant energy resource. Despite the already active exploitation of such deposits, a mathematical model for gas flow in shales does not exist. Such a model is crucial for optimizing the technology of gas recovery. In the present article, a boundary layer problem is formulated and investigated with respect to gas recovery from porous low-permeability inclusions in shales, which are the basic source of gas. Milton Van Dyke was a great master in the field of boundary layer problems. Dedicating this work to his memory, we want to express our belief that Van Dyke\\'s profound ideas and fundamental book Perturbation Methods in Fluid Mechanics (Parabolic Press, 1975) will live on-also in fields very far from the subjects for which they were originally invented. © 2013 US Government.

Appling Laplace Adomian decomposition method for delay differential equations with boundary value problems

Science.gov (United States)

Yousef, Hamood Mohammed; Ismail, Ahmad Izani

2017-11-01

In this paper, Laplace Adomian decomposition method (LADM) was applied to solve Delay differential equations with Boundary Value Problems. The solution is in the form of a convergent series which is easy to compute. This approach is tested on two test problem. The findings obtained exhibit the reliability and efficiency of the proposed method.

Chebyshev Finite Difference Method for Fractional Boundary Value Problems

Directory of Open Access Journals (Sweden)

Boundary

2015-09-01

Full Text Available This paper presents a numerical method for fractional differential equations using Chebyshev finite difference method. The fractional derivatives are described in the Caputo sense. Numerical results show that this method is of high accuracy and is more convenient and efficient for solving boundary value problems involving fractional ordinary differential equations. AMS Subject Classification: 34A08 Keywords and Phrases: Chebyshev polynomials, Gauss-Lobatto points, fractional differential equation, finite difference 1. Introduction The idea of a derivative which interpolates between the familiar integer order derivatives was introduced many years ago and has gained increasing importance only in recent years due to the development of mathematical models of a certain situations in engineering, materials science, control theory, polymer modelling etc. For example see [20, 22, 25, 26]. Most fractional order differential equations describing real life situations, in general do not have exact analytical solutions. Several numerical and approximate analytical methods for ordinary differential equation Received: December 2014; Accepted: March 2015 57 Journal of Mathematical Extension Vol. 9, No. 3, (2015, 57-71 ISSN: 1735-8299 URL: http://www.ijmex.com Chebyshev Finite Difference Method for Fractional Boundary Value Problems H. Azizi Taft Branch, Islamic Azad University Abstract. This paper presents a numerical method for fractional differential equations using Chebyshev finite difference method. The fractional derivative

The complex variable boundary element method: Applications in determining approximative boundaries

Science.gov (United States)

Hromadka, T.V.

1984-01-01

The complex variable boundary element method (CVBEM) is used to determine approximation functions for boundary value problems of the Laplace equation such as occurs in potential theory. By determining an approximative boundary upon which the CVBEM approximator matches the desired constant (level curves) boundary conditions, the CVBEM is found to provide the exact solution throughout the interior of the transformed problem domain. Thus, the acceptability of the CVBEM approximation is determined by the closeness-of-fit of the approximative boundary to the study problem boundary. ?? 1984.

B-spline solution of a singularly perturbed boundary value problem arising in biology

International Nuclear Information System (INIS)

Lin Bin; Li Kaitai; Cheng Zhengxing

2009-01-01

We use B-spline functions to develop a numerical method for solving a singularly perturbed boundary value problem associated with biology science. We use B-spline collocation method, which leads to a tridiagonal linear system. The accuracy of the proposed method is demonstrated by test problems. The numerical result is found in good agreement with exact solution.

On the asymptotic of solutions of elliptic boundary value problems in domains with edges

International Nuclear Information System (INIS)

Nkemzi, B.

2005-10-01

Solutions of elliptic boundary value problems in three-dimensional domains with edges may exhibit singularities. The usual procedure to study these singularities is by the application of the classical Mellin transformation or continuous Fourier transformation. In this paper, we show how the asymptotic behavior of solutions of elliptic boundary value problems in general three-dimensional domains with straight edges can be investigated by means of discrete Fourier transformation. We apply this approach to time-harmonic Maxwell's equations and prove that the singular solutions can fully be described in terms of Fourier series. The representation here can easily be used to approximate three-dimensional stress intensity factors associated with edge singularities. (author)

A simple finite element method for boundary value problems with a Riemann–Liouville derivative

KAUST Repository

Jin, Bangti; Lazarov, Raytcho; Lu, Xiliang; Zhou, Zhi

2016-01-01

© 2015 Elsevier B.V. All rights reserved. We consider a boundary value problem involving a Riemann-Liouville fractional derivative of order α∈(3/2,2) on the unit interval (0,1). The standard Galerkin finite element approximation converges slowly due to the presence of singularity term xα-¹ in the solution representation. In this work, we develop a simple technique, by transforming it into a second-order two-point boundary value problem with nonlocal low order terms, whose solution can reconstruct directly the solution to the original problem. The stability of the variational formulation, and the optimal regularity pickup of the solution are analyzed. A novel Galerkin finite element method with piecewise linear or quadratic finite elements is developed, and ^L2(D) error estimates are provided. The approach is then applied to the corresponding fractional Sturm-Liouville problem, and error estimates of the eigenvalue approximations are given. Extensive numerical results fully confirm our theoretical study.

A simple finite element method for boundary value problems with a Riemann–Liouville derivative

KAUST Repository

Jin, Bangti

2016-02-01

© 2015 Elsevier B.V. All rights reserved. We consider a boundary value problem involving a Riemann-Liouville fractional derivative of order α∈(3/2,2) on the unit interval (0,1). The standard Galerkin finite element approximation converges slowly due to the presence of singularity term xα-¹ in the solution representation. In this work, we develop a simple technique, by transforming it into a second-order two-point boundary value problem with nonlocal low order terms, whose solution can reconstruct directly the solution to the original problem. The stability of the variational formulation, and the optimal regularity pickup of the solution are analyzed. A novel Galerkin finite element method with piecewise linear or quadratic finite elements is developed, and ^L2(D) error estimates are provided. The approach is then applied to the corresponding fractional Sturm-Liouville problem, and error estimates of the eigenvalue approximations are given. Extensive numerical results fully confirm our theoretical study.

Triple solutions for multi-point boundary-value problem with p-Laplace operator

Directory of Open Access Journals (Sweden)

Yansheng Liu

2009-11-01

Full Text Available Using a fixed point theorem due to Avery and Peterson, this article shows the existence of solutions for multi-point boundary-value problem with p-Laplace operator and parameters. Also, we present an example to illustrate the results obtained.

The boundary value problem for discrete analytic functions

KAUST Repository

Skopenkov, Mikhail

2013-06-01

This paper is on further development of discrete complex analysis introduced by R.Isaacs, J.Ferrand, R.Duffin, and C.Mercat. We consider a graph lying in the complex plane and having quadrilateral faces. A function on the vertices is called discrete analytic, if for each face the difference quotients along the two diagonals are equal.We prove that the Dirichlet boundary value problem for the real part of a discrete analytic function has a unique solution. In the case when each face has orthogonal diagonals we prove that this solution uniformly converges to a harmonic function in the scaling limit. This solves a problem of S.Smirnov from 2010. This was proved earlier by R.Courant-K.Friedrichs-H.Lewy and L.Lusternik for square lattices, by D.Chelkak-S.Smirnov and implicitly by P.G.Ciarlet-P.-A.Raviart for rhombic lattices.In particular, our result implies uniform convergence of the finite element method on Delaunay triangulations. This solves a problem of A.Bobenko from 2011. The methodology is based on energy estimates inspired by alternating-current network theory. © 2013 Elsevier Ltd.

Modification of the Riemann problem and the application for the boundary conditions in computational fluid dynamics

Directory of Open Access Journals (Sweden)

Kyncl Martin

2017-01-01

Full Text Available We work with the system of partial differential equations describing the non-stationary compressible turbulent fluid flow. It is a characteristic feature of the hyperbolic equations, that there is a possible raise of discontinuities in solutions, even in the case when the initial conditions are smooth. The fundamental problem in this area is the solution of the so-called Riemann problem for the split Euler equations. It is the elementary problem of the one-dimensional conservation laws with the given initial conditions (LIC - left-hand side, and RIC - right-hand side. The solution of this problem is required in many numerical methods dealing with the 2D/3D fluid flow. The exact (entropy weak solution of this hyperbolical problem cannot be expressed in a closed form, and has to be computed by an iterative process (to given accuracy, therefore various approximations of this solution are being used. The complicated Riemann problem has to be further modified at the close vicinity of boundary, where the LIC is given, while the RIC is not known. Usually, this boundary problem is being linearized, or roughly approximated. The inaccuracies implied by these simplifications may be small, but these have a huge impact on the solution in the whole studied area, especially for the non-stationary flow. Using the thorough analysis of the Riemann problem we show, that the RIC for the local problem can be partially replaced by the suitable complementary conditions. We suggest such complementary conditions accordingly to the desired preference. This way it is possible to construct the boundary conditions by the preference of total values, by preference of pressure, velocity, mass flow, temperature. Further, using the suitable complementary conditions, it is possible to simulate the flow in the vicinity of the diffusible barrier. On the contrary to the initial-value Riemann problem, the solution of such modified problems can be written in the closed form for some

Partial differential equations and boundary-value problems with applications

CERN Document Server

Pinsky, Mark A

2011-01-01

Building on the basic techniques of separation of variables and Fourier series, the book presents the solution of boundary-value problems for basic partial differential equations: the heat equation, wave equation, and Laplace equation, considered in various standard coordinate systems-rectangular, cylindrical, and spherical. Each of the equations is derived in the three-dimensional context; the solutions are organized according to the geometry of the coordinate system, which makes the mathematics especially transparent. Bessel and Legendre functions are studied and used whenever appropriate th

Solutions Stability of Initial Boundary Problem, Modeling of Dynamics of Some Discrete Continuum Mechanical System

Directory of Open Access Journals (Sweden)

D. A. Eliseev

2015-01-01

Full Text Available The solution stability of an initial boundary problem for a linear hybrid system of differential equations, which models the rotation of a rigid body with two elastic rods located in the same plane is studied in the paper. To an axis passing through the mass center of the rigid body perpendicularly to the rods location plane is applied the stabilizing moment proportional to the angle of the system rotation, derivative of the angle, integral of the angle. The external moment provides a feedback. A method of studying the behavior of solutions of the initial boundary problem is proposed. This method allows to exclude from the hybrid system of differential equations partial differential equations, which describe the dynamics of distributed elements of a mechanical system. It allows us to build one equation for an angle of the system rotation. Its characteristic equation defines the stability of solutions of all the system. In the space of feedback-coefficients the areas that provide the asymptotic stability of solutions of the initial boundary problem are built up.

The boundary element method : errors and gridding for problems with hot spots

NARCIS (Netherlands)

Kakuba, G.

2011-01-01

Adaptive gridding methods are of fundamental importance both for industry and academia. As one of the computing methods, the Boundary Element Method (BEM) is used to simulate problems whose fundamental solutions are available. The method is usually characterised as constant elements BEM or linear

Simulation of Thermal Flow Problems via a Hybrid Immersed Boundary-Lattice Boltzmann Method

Directory of Open Access Journals (Sweden)

J. Wu

2012-01-01

Full Text Available A hybrid immersed boundary-lattice Boltzmann method (IB-LBM is presented in this work to simulate the thermal flow problems. In current approach, the flow field is resolved by using our recently developed boundary condition-enforced IB-LBM (Wu and Shu, (2009. The nonslip boundary condition on the solid boundary is enforced in simulation. At the same time, to capture the temperature development, the conventional energy equation is resolved. To model the effect of immersed boundary on temperature field, the heat source term is introduced. Different from previous studies, the heat source term is set as unknown rather than predetermined. Inspired by the idea in (Wu and Shu, (2009, the unknown is calculated in such a way that the temperature at the boundary interpolated from the corrected temperature field accurately satisfies the thermal boundary condition. In addition, based on the resolved temperature correction, an efficient way to compute the local and average Nusselt numbers is also proposed in this work. As compared with traditional implementation, no approximation for temperature gradients is required. To validate the present method, the numerical simulations of forced convection are carried out. The obtained results show good agreement with data in the literature.

Reflected forward-backward SDEs and obstacle problems with boundary conditions

Directory of Open Access Journals (Sweden)

Jin Ma

2001-01-01

Full Text Available In this paper we study a class of forward-backward stochastic differential equations with reflecting boundary conditions (FBSDER for short. More precisely, we consider the case in which the forward component of the FBSDER is restricted to a fixed, convex region, and the backward component will stay, at each fixed time, in a convex region that may depend on time and is possibly random. The solvability of such FBSDER is studied in a fairly general way. We also prove that if the coefficients are all deterministic and the backward equation is one-dimensional, then the adapted solution of such FBSDER will give the viscosity solution of a quasilinear variational inequality (obstacle problem with a Neumann boundary condition. As an application, we study how the solvability of FBSDERs is related to the solvability of an American game option.

A Monotone, Higher-Order Accurate, Fixed-Grid Finite-Volume Method for Advection Problems with Moving Boundaries

NARCIS (Netherlands)

Y.J. Hassen (Yunus); B. Koren (Barry)

2008-01-01

textabstractIn this paper, an accurate method, using a novel immersed-boundary approach, is presented for numerically solving linear, scalar convection problems. As is standard in immersed-boundary methods, moving bodies are embedded in a fixed Cartesian grid. The essence of the present method is

Homotopy Analysis Method for Boundary-Value Problem of Turbo Warrant Pricing under Stochastic Volatility

Directory of Open Access Journals (Sweden)

Hoi Ying Wong

2013-01-01

Full Text Available Turbo warrants are liquidly traded financial derivative securities in over-the-counter and exchange markets in Asia and Europe. The structure of turbo warrants is similar to barrier options, but a lookback rebate will be paid if the barrier is crossed by the underlying asset price. Therefore, the turbo warrant price satisfies a partial differential equation (PDE with a boundary condition that depends on another boundary-value problem (BVP of PDE. Due to the highly complicated structure of turbo warrants, their valuation presents a challenging problem in the field of financial mathematics. This paper applies the homotopy analysis method to construct an analytic pricing formula for turbo warrants under stochastic volatility in a PDE framework.

A uniformly valid approximation algorithm for nonlinear ordinary singular perturbation problems with boundary layer solutions.

Science.gov (United States)

Cengizci, Süleyman; Atay, Mehmet Tarık; Eryılmaz, Aytekin

2016-01-01

This paper is concerned with two-point boundary value problems for singularly perturbed nonlinear ordinary differential equations. The case when the solution only has one boundary layer is examined. An efficient method so called Successive Complementary Expansion Method (SCEM) is used to obtain uniformly valid approximations to this kind of solutions. Four test problems are considered to check the efficiency and accuracy of the proposed method. The numerical results are found in good agreement with exact and existing solutions in literature. The results confirm that SCEM has a superiority over other existing methods in terms of easy-applicability and effectiveness.

A Third-Order p-Laplacian Boundary Value Problem Solved by an SL(3,ℝ Lie-Group Shooting Method

Directory of Open Access Journals (Sweden)

Chein-Shan Liu

2013-01-01

Full Text Available The boundary layer problem for power-law fluid can be recast to a third-order p-Laplacian boundary value problem (BVP. In this paper, we transform the third-order p-Laplacian into a new system which exhibits a Lie-symmetry SL(3,ℝ. Then, the closure property of the Lie-group is used to derive a linear transformation between the boundary values at two ends of a spatial interval. Hence, we can iteratively solve the missing left boundary conditions, which are determined by matching the right boundary conditions through a finer tuning of r∈[0,1]. The present SL(3,ℝ Lie-group shooting method is easily implemented and is efficient to tackle the multiple solutions of the third-order p-Laplacian. When the missing left boundary values can be determined accurately, we can apply the fourth-order Runge-Kutta (RK4 method to obtain a quite accurate numerical solution of the p-Laplacian.

On Nonlinear Inverse Problems of Heat Transfer with Radiation Boundary Conditions: Application to Dehydration of Gypsum Plasterboards Exposed to Fire

OpenAIRE

Belmiloudi, A.; Mahé, F.

2014-01-01

International audience; The paper investigates boundary optimal controls and parameter estimates to the well-posedness nonlinear model of dehydration of thermic problems. We summarize the general formulations for the boundary control for initial-boundary value problem for nonlinear partial differential equations modeling the heat transfer and derive necessary optimality conditions, including the adjoint equation, for the optimal set of parameters minimizing objective functions J. Numerical si...

Sources of spurious force oscillations from an immersed boundary method for moving-body problems

Science.gov (United States)

Lee, Jongho; Kim, Jungwoo; Choi, Haecheon; Yang, Kyung-Soo

2011-04-01

When a discrete-forcing immersed boundary method is applied to moving-body problems, it produces spurious force oscillations on a solid body. In the present study, we identify two sources of these force oscillations. One source is from the spatial discontinuity in the pressure across the immersed boundary when a grid point located inside a solid body becomes that of fluid with a body motion. The addition of mass source/sink together with momentum forcing proposed by Kim et al. [J. Kim, D. Kim, H. Choi, An immersed-boundary finite volume method for simulations of flow in complex geometries, Journal of Computational Physics 171 (2001) 132-150] reduces the spurious force oscillations by alleviating this pressure discontinuity. The other source is from the temporal discontinuity in the velocity at the grid points where fluid becomes solid with a body motion. The magnitude of velocity discontinuity decreases with decreasing the grid spacing near the immersed boundary. Four moving-body problems are simulated by varying the grid spacing at a fixed computational time step and at a constant CFL number, respectively. It is found that the spurious force oscillations decrease with decreasing the grid spacing and increasing the computational time step size, but they depend more on the grid spacing than on the computational time step size.

Analytic solution of boundary-value problems for nonstationary model kinetic equations

International Nuclear Information System (INIS)

Latyshev, A.V.; Yushkanov, A.A.

1993-01-01

A theory for constructing the solutions of boundary-value problems for non-stationary model kinetic equations is constructed. This theory was incorrectly presented equation, separation of the variables is used, this leading to a characteristic equation. Eigenfunctions are found in the space of generalized functions, and the eigenvalue spectrum is investigated. An existence and uniqueness theorem for the expansion of the Laplace transform of the solution with respect to the eigenfunctions is proved. The proof is constructive and gives explicit expressions for the expansion coefficients. An application to the Rayleigh problem is obtained, and the corresponding result of Cercignani is corrected

Three-Field Modelling of Nonlinear Nonsmooth Boundary Value Problems and Stability of Differential Mixed Variational Inequalities

Directory of Open Access Journals (Sweden)

J. Gwinner

2013-01-01

Full Text Available The purpose of this paper is twofold. Firstly we consider nonlinear nonsmooth elliptic boundary value problems, and also related parabolic initial boundary value problems that model in a simplified way steady-state unilateral contact with Tresca friction in solid mechanics, respectively, stem from nonlinear transient heat conduction with unilateral boundary conditions. Here a recent duality approach, that augments the classical Babuška-Brezzi saddle point formulation for mixed variational problems to twofold saddle point formulations, is extended to the nonsmooth problems under consideration. This approach leads to variational inequalities of mixed form for three coupled fields as unknowns and to related differential mixed variational inequalities in the time-dependent case. Secondly we are concerned with the stability of the solution set of a general class of differential mixed variational inequalities. Here we present a novel upper set convergence result with respect to perturbations in the data, including perturbations of the associated nonlinear maps, the nonsmooth convex functionals, and the convex constraint set. We employ epiconvergence for the convergence of the functionals and Mosco convergence for set convergence. We impose weak convergence assumptions on the perturbed maps using the monotonicity method of Browder and Minty.

Existence of solutions to fractional boundary-value problems with a parameter

Directory of Open Access Journals (Sweden)

Ya-Ning Li

2013-06-01

Full Text Available This article concerns the existence of solutions to the fractional boundary-value problem $$displaylines{ -frac{d}{dt} ig(frac{1}{2} {}_0D_t^{-eta}+ frac{1}{2}{}_tD_{T}^{-eta}igu'(t=lambda u(t+abla F(t,u(t,quad hbox{a.e. } tin[0,T], cr u(0=0,quad u(T=0. }$$ First for the eigenvalue problem associated with it, we prove that there is a sequence of positive and increasing real eigenvalues; a characterization of the first eigenvalue is also given. Then under different assumptions on the nonlinearity F(t,u, we show the existence of weak solutions of the problem when $lambda$ lies in various intervals. Our main tools are variational methods and critical point theorems.

Conforming discretizations of boundary element solutions to the electroencephalography forward problem

Science.gov (United States)

Rahmouni, Lyes; Adrian, Simon B.; Cools, Kristof; Andriulli, Francesco P.

2018-01-01

In this paper, we present a new discretization strategy for the boundary element formulation of the Electroencephalography (EEG) forward problem. Boundary integral formulations, classically solved with the Boundary Element Method (BEM), are widely used in high resolution EEG imaging because of their recognized advantages, in several real case scenarios, in terms of numerical stability and effectiveness when compared with other differential equation based techniques. Unfortunately, however, it is widely reported in literature that the accuracy of standard BEM schemes for the forward EEG problem is often limited, especially when the current source density is dipolar and its location approaches one of the brain boundary surfaces. This is a particularly limiting problem given that during an high-resolution EEG imaging procedure, several EEG forward problem solutions are required, for which the source currents are near or on top of a boundary surface. This work will first present an analysis of standardly and classically discretized EEG forward problem operators, reporting on a theoretical issue of some of the formulations that have been used so far in the community. We report on the fact that several standardly used discretizations of these formulations are consistent only with an L2-framework, requiring the expansion term to be a square integrable function (i.e., in a Petrov-Galerkin scheme with expansion and testing functions). Instead, those techniques are not consistent when a more appropriate mapping in terms of fractional-order Sobolev spaces is considered. Such a mapping allows the expansion function term to be a less regular function, thus sensibly reducing the need for mesh refinements and low-precisions handling strategies that are currently required. These more favorable mappings, however, require a different and conforming discretization, which must be suitably adapted to them. In order to appropriately fulfill this requirement, we adopt a mixed

Efficient evaluation of influence coefficients in three-dimensional extended boundary-node method for potential problems

International Nuclear Information System (INIS)

Itoh, Taku; Saitoh, Ayumu; Kamitani, Atsushi; Nakamura, Hiroaki

2011-01-01

For the purpose of speed-up of the three-dimensional eXtended Boundary-Node Method (X-BNM), an efficient algorithm for evaluating influence coefficients has been developed. The algorithm can be easily implemented into the X-BNM without using any integration cells. By applying the resulting X-BNM to the Laplace problem, the performance of the algorithm is numerically investigated. The numerical experiments show that, by using the algorithm, computational costs for evaluating influence coefficients in the X-BNM are reduced considerably. Especially for a large-sized problem, the algorithm is efficiently performed, and the computational costs of the X-BNM are close to those of the Boundary-Element Method (BEM). In addition, for the problem, the X-BNM shows almost the same accuracy as that of the BEM. (author)

Three symmetric positive solutions of fourth-order singular nonlocal boundary value problems

Directory of Open Access Journals (Sweden)

Fuyi Xu

2011-12-01

Full Text Available In this paper, we study the existence of three positive solutions of fourth-order singular nonlocal boundary value problems. We show that there exist triple symmetric positive solutions by using Leggett-Williams fixed-point theorem. The conclusions in this paper essentially extend and improve some known results.

Elliptic boundary value problems with fractional regularity data the first order approach

CERN Document Server

Amenta, Alex

2018-01-01

In this monograph the authors study the well-posedness of boundary value problems of Dirichlet and Neumann type for elliptic systems on the upper half-space with coefficients independent of the transversal variable and with boundary data in fractional Hardy-Sobolev and Besov spaces. The authors use the so-called "first order approach" which uses minimal assumptions on the coefficients and thus allows for complex coefficients and for systems of equations. This self-contained exposition of the first order approach offers new results with detailed proofs in a clear and accessible way and will become a valuable reference for graduate students and researchers working in partial differential equations and harmonic analysis.

Symmetry analysis and exact solutions of one class of (1+3)-dimensional boundary-value problems of the Stefan type

OpenAIRE

Kovalenko, S. S.

2014-01-01

We present the group classification of one class of (1+3)-dimensional nonlinear boundary-value problems of the Stefan type that simulate the processes of melting and evaporation of metals. The results obtained are used for the construction of the exact solution of one boundary-value problem from the class under study.

Neurocognitive Problems in Children and Adolescents With Depression Using the CHC Theory and the WJ-III.

Science.gov (United States)

Basnet, Pravesh; Noggle, Chad A; Dean, Raymond S

2015-01-01

Depression has been commonly associated with both subjective complaints and objectively measured problems in cognition. Most commonly discussed in relation to the adult population, growing evidence has supported the idea that children and adolescents experience cognitive problems in relation to depression. The purpose of this study was to further examine the negative influence of depression on the cognitive functioning of children and adolescents. Additionally, the present study evaluated the sensitivity of the Woodcock-Johnson III Test of Cognitive Abilities (WJ-III-COG) and, in turn, the Cattell-Horn-Carroll (CHC) theory in measuring cognitive problems related to depression in children and adolescents. Participants included 420 children and adolescents aged 8 to 18 years old (M = 13.09, SD = 2.95) with a clinical diagnosis of depression. Comparisons were made against the normative mean. All participants completed 11 subtests of the WJ-III-COG including Visual-Auditory Learning, Spatial Relations, Sound Blending, Concept Formations, Visual Matching, Numbers Reversed, Auditory-Working Memory, Picture Recognition, Analysis Synthesis, Decision Speed, and Memory for Words. Children and adolescents with depression demonstrated significantly lower performance on subtests related to learning and memory (long-term retrieval), attentional capacity, working memory, reasoning, and processing speed. No problems were noted on subtests related to visual-spatial thinking and auditory processing. Findings suggested sensitivity of the WJ-III-COG and CHC theory in identifying cognitive problems associated with depression in children and adolescents.

Exact solution for a two-phase Stefan problem with variable latent heat and a convective boundary condition at the fixed face

Science.gov (United States)

Bollati, Julieta; Tarzia, Domingo A.

2018-04-01

Recently, in Tarzia (Thermal Sci 21A:1-11, 2017) for the classical two-phase Lamé-Clapeyron-Stefan problem an equivalence between the temperature and convective boundary conditions at the fixed face under a certain restriction was obtained. Motivated by this article we study the two-phase Stefan problem for a semi-infinite material with a latent heat defined as a power function of the position and a convective boundary condition at the fixed face. An exact solution is constructed using Kummer functions in case that an inequality for the convective transfer coefficient is satisfied generalizing recent works for the corresponding one-phase free boundary problem. We also consider the limit to our problem when that coefficient goes to infinity obtaining a new free boundary problem, which has been recently studied in Zhou et al. (J Eng Math 2017. https://doi.org/10.1007/s10665-017-9921-y).

Some problems in the magnetohydrodynamics of liquid metals

International Nuclear Information System (INIS)

Moffatt, H.K.

1978-01-01

When electric currents are caused to flow in an electrically conducting fluid, either by the external application of time-periodic magnetic fields or by the application of large electric potential gradients at the boundary, the associated Lorentz force is in general rotational and a fluid motion, which may be laminar or turbulent, is in general established. Three prototype problems, on which some progress has been made over the last decade, are reviewed: (i) the problem of the generation of rotation in a liquid metal by the application of a rotating magnetic field; (ii) the generation of cellular motion by the application of an alternating field of fixed direction; and (iii) the problem of the generation of fluid motion by the injection of steady current at a point electrode on the fluid boundary. All three problems are of importance in molten metal technology. (author)

An Adaptive Pseudospectral Method for Fractional Order Boundary Value Problems

Directory of Open Access Journals (Sweden)

Mohammad Maleki

2012-01-01

Full Text Available An adaptive pseudospectral method is presented for solving a class of multiterm fractional boundary value problems (FBVP which involve Caputo-type fractional derivatives. The multiterm FBVP is first converted into a singular Volterra integrodifferential equation (SVIDE. By dividing the interval of the problem to subintervals, the unknown function is approximated using a piecewise interpolation polynomial with unknown coefficients which is based on shifted Legendre-Gauss (ShLG collocation points. Then the problem is reduced to a system of algebraic equations, thus greatly simplifying the problem. Further, some additional conditions are considered to maintain the continuity of the approximate solution and its derivatives at the interface of subintervals. In order to convert the singular integrals of SVIDE into nonsingular ones, integration by parts is utilized. In the method developed in this paper, the accuracy can be improved either by increasing the number of subintervals or by increasing the degree of the polynomial on each subinterval. Using several examples including Bagley-Torvik equation the proposed method is shown to be efficient and accurate.

Asymptotic behaviour of solutions of the first boundary-value problem for strongly hyperbolic systems near a conical point at the boundary of the domain

International Nuclear Information System (INIS)

Hung, Nguyen M

1999-01-01

An existence and uniqueness theorem for generalized solutions of the first initial-boundary-value problem for strongly hyperbolic systems in bounded domains is established. The question of estimates in Sobolev spaces of the derivatives with respect to time of the generalized solution is discussed. It is shown that the smoothness of generalized solutions with respect to time is independent of the structure of the boundary of the domain but depends on the coefficients of the right-hand side. Results on the smoothness of the generalized solution and its asymptotic behaviour in a neighbourhood of a conical boundary point are also obtained

Variational multiscale enrichment method with mixed boundary conditions for elasto-viscoplastic problems

Science.gov (United States)

Zhang, Shuhai; Oskay, Caglar

2015-04-01

This manuscript presents the formulation and implementation of the variational multiscale enrichment (VME) method for the analysis of elasto-viscoplastic problems. VME is a global-local approach that allows accurate fine scale representation at small subdomains, where important physical phenomena are likely to occur. The response within far-fields is idealized using a coarse scale representation. The fine scale representation not only approximates the coarse grid residual, but also accounts for the material heterogeneity. A one-parameter family of mixed boundary conditions that range from Dirichlet to Neumann is employed to study the effect of the choice of the boundary conditions at the fine scale on accuracy. The inelastic material behavior is modeled using Perzyna type viscoplasticity coupled with flow stress evolution idealized by the Johnson-Cook model. Numerical verifications are performed to assess the performance of the proposed approach against the direct finite element simulations. The results of verification studies demonstrate that VME with proper boundary conditions accurately model the inelastic response accounting for material heterogeneity.

Analytic solutions to a family of boundary-value problems for Ginsburg-Landau type equations

Science.gov (United States)

Vassilev, V. M.; Dantchev, D. M.; Djondjorov, P. A.

2017-10-01

We consider a two-parameter family of nonlinear ordinary differential equations describing the behavior of a critical thermodynamic system, e.g., a binary liquid mixture, of film geometry within the framework of the Ginzburg-Landau theory by means of the order-parameter. We focus on the case in which the confining surfaces are strongly adsorbing but prefer different components of the mixture, i.e., the order-parameter tends to infinity at one of the boundaries and to minus infinity at the other one. We assume that the boundaries of the system are positioned at a finite distance from each other and give analytic solutions to the corresponding boundary-value problems in terms of Weierstrass and Jacobi elliptic functions.

Variational Homotopy Perturbation Method for Solving Higher Dimensional Initial Boundary Value Problems

Directory of Open Access Journals (Sweden)

Muhammad Aslam Noor

2008-01-01

Full Text Available We suggest and analyze a technique by combining the variational iteration method and the homotopy perturbation method. This method is called the variational homotopy perturbation method (VHPM. We use this method for solving higher dimensional initial boundary value problems with variable coefficients. The developed algorithm is quite efficient and is practically well suited for use in these problems. The proposed scheme finds the solution without any discritization, transformation, or restrictive assumptions and avoids the round-off errors. Several examples are given to check the reliability and efficiency of the proposed technique.

Layer potentials and boundary-value problems for second order elliptic operators with data in Besov spaces

CERN Document Server

Barton, Ariel

2016-01-01

This monograph presents a comprehensive treatment of second order divergence form elliptic operators with bounded measurable t-independent coefficients in spaces of fractional smoothness, in Besov and weighted L^p classes. The authors establish: (1) Mapping properties for the double and single layer potentials, as well as the Newton potential; (2) Extrapolation-type solvability results: the fact that solvability of the Dirichlet or Neumann boundary value problem at any given L^p space automatically assures their solvability in an extended range of Besov spaces; (3) Well-posedness for the non-homogeneous boundary value problems. In particular, the authors prove well-posedness of the non-homogeneous Dirichlet problem with data in Besov spaces for operators with real, not necessarily symmetric, coefficients.

Global-in-time smoothness of solutions to the vacuum free boundary problem for compressible isentropic Navier–Stokes equations

International Nuclear Information System (INIS)

Zeng, Huihui

2015-01-01

In this paper we establish the global existence of smooth solutions to vacuum free boundary problems of the one-dimensional compressible isentropic Navier–Stokes equations for which the smoothness extends all the way to the boundaries. The results obtained in this work include the physical vacuum for which the sound speed is C 1/2 -Hölder continuous near the vacuum boundaries when 1 < γ < 3. The novelty of this result is its global-in-time regularity which is in contrast to the previous main results of global weak solutions in the literature. Moreover, in previous studies of the one-dimensional free boundary problems of compressible Navier–Stokes equations, the Lagrangian mass coordinates method has often been used, but in the present work the particle path (flow trajectory) method is adopted, which has the advantage that the particle paths and, in particular, the free boundaries can be traced. (paper)

The boundary element method for the solution of the multidimensional inverse heat conduction problem

International Nuclear Information System (INIS)

Lagier, Guy-Laurent

1999-01-01

This work focuses on the solution of the inverse heat conduction problem (IHCP), which consists in the determination of boundary conditions from a given set of internal temperature measurements. This problem is difficult to solve due to its ill-posedness and high sensitivity to measurement error. As a consequence, numerical regularization procedures are required to solve this problem. However, most of these methods depend on the dimension and the nature, stationary or transient, of the problem. Furthermore, these methods introduce parameters, called hyper-parameters, which have to be chosen optimally, but can not be determined a priori. So, a new general method is proposed for solving the IHCP. This method is based on a Boundary Element Method formulation, and the use of the Singular Values Decomposition as a regularization procedure. Thanks to this method, it's possible to identify and eliminate the directions of the solution where the measurement error plays the major role. This algorithm is first validated on two-dimensional stationary and one-dimensional transient problems. Some criteria are presented in order to choose the hyper-parameters. Then, the methodology is applied to two-dimensional and three-dimensional, theoretical or experimental, problems. The results are compared with those obtained by a standard method and show the accuracy of the method, its generality, and the validity of the proposed criteria. (author) [fr

Discrete quintic spline for boundary value problem in plate deflation theory

Science.gov (United States)

Wong, Patricia J. Y.

2017-07-01

We propose a numerical scheme for a fourth-order boundary value problem arising from plate deflation theory. The scheme involves a discrete quintic spline, and it is of order 4 if a parameter takes a specific value, else it is of order 2. We also present a well known numerical example to illustrate the efficiency of our method as well as to compare with other numerical methods proposed in the literature.

The Boundary Function Method. Fundamentals

Science.gov (United States)

Kot, V. A.

2017-03-01

The boundary function method is proposed for solving applied problems of mathematical physics in the region defined by a partial differential equation of the general form involving constant or variable coefficients with a Dirichlet, Neumann, or Robin boundary condition. In this method, the desired function is defined by a power polynomial, and a boundary function represented in the form of the desired function or its derivative at one of the boundary points is introduced. Different sequences of boundary equations have been set up with the use of differential operators. Systems of linear algebraic equations constructed on the basis of these sequences allow one to determine the coefficients of a power polynomial. Constitutive equations have been derived for initial boundary-value problems of all the main types. With these equations, an initial boundary-value problem is transformed into the Cauchy problem for the boundary function. The determination of the boundary function by its derivative with respect to the time coordinate completes the solution of the problem.

On the Existence and Uniqueness of Rv-Generalized Solution for Dirichlet Problem with Singularity on All Boundary

Directory of Open Access Journals (Sweden)

V. Rukavishnikov

2014-01-01

Full Text Available The existence and uniqueness of the Rv-generalized solution for the first boundary value problem and a second order elliptic equation with coordinated and uncoordinated degeneracy of input data and with strong singularity solution on all boundary of a two-dimensional domain are established.

Variational Iteration Method for Fifth-Order Boundary Value Problems Using He's Polynomials

Directory of Open Access Journals (Sweden)

Muhammad Aslam Noor

2008-01-01

Full Text Available We apply the variational iteration method using He's polynomials (VIMHP for solving the fifth-order boundary value problems. The proposed method is an elegant combination of variational iteration and the homotopy perturbation methods and is mainly due to Ghorbani (2007. The suggested algorithm is quite efficient and is practically well suited for use in these problems. The proposed iterative scheme finds the solution without any discritization, linearization, or restrictive assumptions. Several examples are given to verify the reliability and efficiency of the method. The fact that the proposed technique solves nonlinear problems without using Adomian's polynomials can be considered as a clear advantage of this algorithm over the decomposition method.

Mixed problem with integral boundary condition for a high order mixed type partial differential equation

OpenAIRE

M. Denche; A. L. Marhoune

2003-01-01

In this paper, we study a mixed problem with integral boundary conditions for a high order partial differential equation of mixed type. We prove the existence and uniqueness of the solution. The proof is based on energy inequality, and on the density of the range of the operator generated by the considered problem.

The use of Fourier eigen transform to the boundary element method for transient elastodynamic problems

International Nuclear Information System (INIS)

Ji, X.; Chen, Y.M.

1989-01-01

The boundary element method (BEM) is developed from the boundary integral equation method and the discretization techniques. Compared with other numerical method, BEM has been shown to be a versatile and efficient method for a wide variety of engineering problems, including the wave propagation in elastic media. The first formulation and solution of the transient elastodynamic problem by combining BEM and Laplace transform is due to Cruse. Further improvement was achieved by introducing Durbin's method instead of Papoulis method of numerical Laplace inverse transform. However, a great deal of computer time is still needed for the inverse transformation. The alternative integral transform approach is BEM combining with Fourier transform. The numerical Fourier inverse transformation is also computer time consuming, even if the fast Fourier transform is used. In the present paper, the authors use BEM combining with Fourier transform and Fourier eigen transform (FET). The new approach is very attractive in saving on computer time. This paper illustrates the application of FET to BEM of 2-dimensional transient elastodynamic problem. The example of a half plane subjected to a discontinuous boundary load is solved on ELXSI 6400 computer. The CPU time is less than one minute. If Laplace or Fourier transform is adopted, the CPU time will be more than 10 minutes

The quantum-field renormalization group in the problem of a growing phase boundary

International Nuclear Information System (INIS)

Antonov, N.V.; Vasil'ev, A.N.

1995-01-01

Within the quantum-field renormalization-group approach we examine the stochastic equation discussed by S.I. Pavlik in describing a randomly growing phase boundary. We show that, in contrast to Pavlik's assertion, the model is not multiplicatively renormalizable and that its consistent renormalization-group analysis requires introducing an infinite number of counterterms and the respective coupling constants (open-quotes chargeclose quotes). An explicit calculation in the one-loop approximation shows that a two-dimensional surface of renormalization-group points exits in the infinite-dimensional charge space. If the surface contains an infrared stability region, the problem allows for scaling with the nonuniversal critical dimensionalities of the height of the phase boundary and time, δ h and δ t , which satisfy the exact relationship 2 δ h = δ t + d, where d is the dimensionality of the phase boundary. 23 refs., 1 tab

Boundary conditions for free surface inlet and outlet problems

KAUST Repository

Taroni, M.

2012-08-10

We investigate and compare the boundary conditions that are to be applied to free-surface problems involving inlet and outlets of Newtonian fluid, typically found in coating processes. The flux of fluid is a priori known at an inlet, but unknown at an outlet, where it is governed by the local behaviour near the film-forming meniscus. In the limit of vanishing capillary number Ca it is well known that the flux scales with Ca 2/3, but this classical result is non-uniform as the contact angle approaches π. By examining this limit we find a solution that is uniformly valid for all contact angles. Furthermore, by considering the far-field behaviour of the free surface we show that there exists a critical capillary number above which the problem at an inlet becomes over-determined. The implications of this result for the modelling of coating flows are discussed. © 2012 Cambridge University Press.

A classical Perron method for existence of smooth solutions to boundary value and obstacle problems for degenerate-elliptic operators via holomorphic maps

Science.gov (United States)

Feehan, Paul M. N.

2017-09-01

We prove existence of solutions to boundary value problems and obstacle problems for degenerate-elliptic, linear, second-order partial differential operators with partial Dirichlet boundary conditions using a new version of the Perron method. The elliptic operators considered have a degeneracy along a portion of the domain boundary which is similar to the degeneracy of a model linear operator identified by Daskalopoulos and Hamilton [9] in their study of the porous medium equation or the degeneracy of the Heston operator [21] in mathematical finance. Existence of a solution to the partial Dirichlet problem on a half-ball, where the operator becomes degenerate on the flat boundary and a Dirichlet condition is only imposed on the spherical boundary, provides the key additional ingredient required for our Perron method. Surprisingly, proving existence of a solution to this partial Dirichlet problem with ;mixed; boundary conditions on a half-ball is more challenging than one might expect. Due to the difficulty in developing a global Schauder estimate and due to compatibility conditions arising where the ;degenerate; and ;non-degenerate boundaries; touch, one cannot directly apply the continuity or approximate solution methods. However, in dimension two, there is a holomorphic map from the half-disk onto the infinite strip in the complex plane and one can extend this definition to higher dimensions to give a diffeomorphism from the half-ball onto the infinite ;slab;. The solution to the partial Dirichlet problem on the half-ball can thus be converted to a partial Dirichlet problem on the slab, albeit for an operator which now has exponentially growing coefficients. The required Schauder regularity theory and existence of a solution to the partial Dirichlet problem on the slab can nevertheless be obtained using previous work of the author and C. Pop [16]. Our Perron method relies on weak and strong maximum principles for degenerate-elliptic operators, concepts of

Unique solvability of some two-point boundary value problems for linear functional differential equations with singularities

Czech Academy of Sciences Publication Activity Database

Rontó, András; Samoilenko, A. M.

2007-01-01

Roč. 41, - (2007), s. 115-136 ISSN 1512-0015 R&D Projects: GA ČR(CZ) GA201/06/0254 Institutional research plan: CEZ:AV0Z10190503 Keywords : two-point problem * functional differential equation * singular boundary problem Subject RIV: BA - General Mathematics

Investigation of solutions of state-dependent multi-impulsive boundary value problems

Czech Academy of Sciences Publication Activity Database

Rontó, András; Rachůnková, I.; Rontó, M.; Rachůnek, L.

2017-01-01

Roč. 24, č. 2 (2017), s. 287-312 ISSN 1072-947X R&D Projects: GA ČR(CZ) GA14-06958S Institutional support: RVO:67985840 Keywords : state-dependent multi-impulsive systems * non-linear boundary value problem * parametrization technique Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.290, year: 2016 https://www.degruyter.com/view/j/gmj.2017.24.issue-2/gmj-2016-0084/gmj-2016-0084. xml

Investigation of solutions of state-dependent multi-impulsive boundary value problems

Czech Academy of Sciences Publication Activity Database

Rontó, András; Rachůnková, I.; Rontó, M.; Rachůnek, L.

2017-01-01

Roč. 24, č. 2 (2017), s. 287-312 ISSN 1072-947X R&D Projects: GA ČR(CZ) GA14-06958S Institutional support: RVO:67985840 Keywords : state-dependent multi-impulsive systems * non-linear boundary value problem * parametrization technique Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.290, year: 2016 https://www.degruyter.com/view/j/gmj.2017.24.issue-2/gmj-2016-0084/gmj-2016-0084.xml

Integral boundary-value problem for impulsive fractional functional differential equations with infinite delay

Directory of Open Access Journals (Sweden)

Archana Chauhan

2012-12-01

Full Text Available In this article, we establish a general framework for finding solutions for impulsive fractional integral boundary-value problems. Then, we prove the existence and uniqueness of solutions by applying well known fixed point theorems. The obtained results are illustrated with an example for their feasibility.

Geopotential coefficient determination and the gravimetric boundary value problem: A new approach

Science.gov (United States)

Sjoeberg, Lars E.

1989-01-01

New integral formulas to determine geopotential coefficients from terrestrial gravity and satellite altimetry data are given. The formulas are based on the integration of data over the non-spherical surface of the Earth. The effect of the topography to low degrees and orders of coefficients is estimated numerically. Formulas for the solution of the gravimetric boundary value problem are derived.

Electrodynamics as a problem of eigenvalues. III. Maxwell operator and perturbation theory

International Nuclear Information System (INIS)

Yudin, L.A.; Kapchinsky, M.I.; Korenev, I.L.; Efimov, S.P.

1996-01-01

This part of the work deals with perturbation theory. The standard technique permits us to solve problems of volume disturbances in free space, waveguides, and cavities, of disturbances of boundary conditions caused by finite conductivity of walls, and of disturbances of the boundary shape. Expressions for frequency and wave number shifts caused by the perturbation are obtained by the unified method. Several examples are presented: wave amplification in the waveguide filled with a resonant medium like an electron beam in a longitudinal magnetic field, the frequency shift (and increment/decrement) of a cavity loaded by an electron beam, and the field disturbance in the waveguide caused by finite conductivity of the walls. copyright 1996 American Institute of Physics

Application of the HN method to the critical slab problem for reflecting boundary conditions

International Nuclear Information System (INIS)

Tuereci, R.G.; Guelecyuez, M.C.; Kaskas, A.; Tezcan, C.

2004-01-01

The recently developed H N method is used to solve the critical slab problem for a slab which is surrounded by a reflector. In the special case for R=0 (the reflection coefficient) the problem reduces to the one under vacuum boundary conditions. It is shown that the method is concise and leads to fast converging numerical results. The presented numerical results are compared with the data available in literature

On Existence of Solutions to the Caputo Type Fractional Order Three-Point Boundary Value Problems

Directory of Open Access Journals (Sweden)

B.M.B. Krushna

2016-10-01

Full Text Available In this paper, we establish the existence of solutions to the fractional order three-point boundary value problems by utilizing Banach contraction principle and Schaefer's fixed point theorem.

The nonlocal boundary value problems for strongly singular higher-order nonlinear functional-differential equations

Czech Academy of Sciences Publication Activity Database

Mukhigulashvili, Sulkhan

-, č. 35 (2015), s. 23-50 ISSN 1126-8042 Institutional support: RVO:67985840 Keywords : higher order functional differential equations * Dirichlet boundary value problem * strong singularity Subject RIV: BA - General Mathematics http://ijpam.uniud.it/online_issue/201535/03-Mukhigulashvili.pdf

Polyharmonic boundary value problems positivity preserving and nonlinear higher order elliptic equations in bounded domains

CERN Document Server

Gazzola, Filippo; Sweers, Guido

2010-01-01

This monograph covers higher order linear and nonlinear elliptic boundary value problems in bounded domains, mainly with the biharmonic or poly-harmonic operator as leading principal part. Underlying models and, in particular, the role of different boundary conditions are explained in detail. As for linear problems, after a brief summary of the existence theory and Lp and Schauder estimates, the focus is on positivity or - since, in contrast to second order equations, a general form of a comparison principle does not exist - on “near positivity.” The required kernel estimates are also presented in detail. As for nonlinear problems, several techniques well-known from second order equations cannot be utilized and have to be replaced by new and different methods. Subcritical, critical and supercritical nonlinearities are discussed and various existence and nonexistence results are proved. The interplay with the positivity topic from the first part is emphasized and, moreover, a far-reaching Gidas-Ni-Nirenbe...

A boundary integral method for a dynamic, transient mode I crack problem with viscoelastic cohesive zone

KAUST Repository

Leise, Tanya L.

2009-08-19

We consider the problem of the dynamic, transient propagation of a semi-infinite, mode I crack in an infinite elastic body with a nonlinear, viscoelastic cohesize zone. Our problem formulation includes boundary conditions that preclude crack face interpenetration, in contrast to the usual mode I boundary conditions that assume all unloaded crack faces are stress-free. The nonlinear viscoelastic cohesive zone behavior is motivated by dynamic fracture in brittle polymers in which crack propagation is preceeded by significant crazing in a thin region surrounding the crack tip. We present a combined analytical/numerical solution method that involves reducing the problem to a Dirichlet-to-Neumann map along the crack face plane, resulting in a differo-integral equation relating the displacement and stress along the crack faces and within the cohesive zone. © 2009 Springer Science+Business Media B.V.

Existence of positive solutions for a multi-point four-order boundary-value problem

Directory of Open Access Journals (Sweden)

Le Xuan Truong

2011-10-01

Full Text Available The article shows sufficient conditions for the existence of positive solutions to a multi-point boundary-value problem for a fourth-order differential equation. Our main tools are the Guo-Krasnoselskii fixed point theorem and the monotone iterative technique. We also show that the set of positive solutions is compact.

Remark on periodic boundary-value problem for second-order linear ordinary differential equations

Czech Academy of Sciences Publication Activity Database

Dosoudilová, M.; Lomtatidze, Alexander

2018-01-01

Roč. 2018, č. 13 (2018), s. 1-7 ISSN 1072-6691 Institutional support: RVO:67985840 Keywords : second-order linear equation * periodic boundary value problem * unique solvability Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.954, year: 2016 https://ejde.math.txstate.edu/Volumes/2018/13/abstr.html

Initial-Boundary Value Problem Solution of the Nonlinear Shallow-water Wave Equations

Science.gov (United States)

Kanoglu, U.; Aydin, B.

2014-12-01

The hodograph transformation solutions of the one-dimensional nonlinear shallow-water wave (NSW) equations are usually obtained through integral transform techniques such as Fourier-Bessel transforms. However, the original formulation of Carrier and Greenspan (1958 J Fluid Mech) and its variant Carrier et al. (2003 J Fluid Mech) involve evaluation integrals. Since elliptic integrals are highly singular as discussed in Carrier et al. (2003), this solution methodology requires either approximation of the associated integrands by smooth functions or selection of regular initial/boundary data. It should be noted that Kanoglu (2004 J Fluid Mech) partly resolves this issue by simplifying the resulting integrals in closed form. Here, the hodograph transform approach is coupled with the classical eigenfunction expansion method rather than integral transform techniques and a new analytical model for nonlinear long wave propagation over a plane beach is derived. This approach is based on the solution methodology used in Aydın & Kanoglu (2007 CMES-Comp Model Eng) for wind set-down relaxation problem. In contrast to classical initial- or boundary-value problem solutions, here, the NSW equations are formulated to yield an initial-boundary value problem (IBVP) solution. In general, initial wave profile with nonzero initial velocity distribution is assumed and the flow variables are given in the form of Fourier-Bessel series. The results reveal that the developed method allows accurate estimation of the spatial and temporal variation of the flow quantities, i.e., free-surface height and depth-averaged velocity, with much less computational effort compared to the integral transform techniques such as Carrier et al. (2003), Kanoglu (2004), Tinti & Tonini (2005 J Fluid Mech), and Kanoglu & Synolakis (2006 Phys Rev Lett). Acknowledgments: This work is funded by project ASTARTE- Assessment, STrategy And Risk Reduction for Tsunamis in Europe. Grant 603839, 7th FP (ENV.2013.6.4-3 ENV

New implementation method for essential boundary condition to extended element-free Galerkin method. Application to nonlinear problem

International Nuclear Information System (INIS)

Saitoh, Ayumu; Matsui, Nobuyuki; Itoh, Taku; Kamitani, Atsushi; Nakamura, Hiroaki

2011-01-01

A new method has been proposed for implementing essential boundary conditions to the Element-Free Galerkin Method (EFGM) without using the Lagrange multiplier. Furthermore, the performance of the proposed method has been investigated for a nonlinear Poisson problem. The results of computations show that, as interpolation functions become closer to delta functions, the accuracy of the solution is improved on the boundary. In addition, the accuracy of the proposed method is higher than that of the conventional EFGM. Therefore, it might be concluded that the proposed method is useful for solving the nonlinear Poisson problem. (author)

Solving Singular Two-Point Boundary Value Problems Using Continuous Genetic Algorithm

Directory of Open Access Journals (Sweden)

Omar Abu Arqub

2012-01-01

Full Text Available In this paper, the continuous genetic algorithm is applied for the solution of singular two-point boundary value problems, where smooth solution curves are used throughout the evolution of the algorithm to obtain the required nodal values. The proposed technique might be considered as a variation of the finite difference method in the sense that each of the derivatives is replaced by an appropriate difference quotient approximation. This novel approach possesses main advantages; it can be applied without any limitation on the nature of the problem, the type of singularity, and the number of mesh points. Numerical examples are included to demonstrate the accuracy, applicability, and generality of the presented technique. The results reveal that the algorithm is very effective, straightforward, and simple.

Existence of the Optimal Control for Stochastic Boundary Control Problems Governed by Semilinear Parabolic Equations

Directory of Open Access Journals (Sweden)

Weifeng Wang

2014-01-01

Full Text Available We study an optimal control problem governed by a semilinear parabolic equation, whose control variable is contained only in the boundary condition. An existence theorem for the optimal control is obtained.

Problems and chances for probabilistic fracture mechanics in the analysis of steel pressure boundary reliability

Energy Technology Data Exchange (ETDEWEB)

Staat, M [Forschungszentrum Juelich GmbH (Germany). Inst. fuer Sicherheitsforschung und Reaktortechnik

1996-12-01

It is shown that the difficulty for probabilistic fracture mechanics (PFM) is the general problem of the high reliability of a small population. There is no way around the problem as yet. Therefore what PFM can contribute to the reliability of steel pressure boundaries is demonstrated with the example of a typical reactor pressure vessel and critically discussed. Although no method is distinguishable that could give exact failure probabilities, PFM has several additional chances. Upper limits for failure probability may be obtained together with trends for design and operating conditions. Further, PFM can identify the most sensitive parameters, improved control of which would increase reliability. Thus PFM should play a vital role in the analysis of steel pressure boundaries despite all shortcomings. (author). 19 refs, 7 figs, 1 tab.

Problems and chances for probabilistic fracture mechanics in the analysis of steel pressure boundary reliability

International Nuclear Information System (INIS)

Staat, M.

1996-01-01

It is shown that the difficulty for probabilistic fracture mechanics (PFM) is the general problem of the high reliability of a small population. There is no way around the problem as yet. Therefore what PFM can contribute to the reliability of steel pressure boundaries is demonstrated with the example of a typical reactor pressure vessel and critically discussed. Although no method is distinguishable that could give exact failure probabilities, PFM has several additional chances. Upper limits for failure probability may be obtained together with trends for design and operating conditions. Further, PFM can identify the most sensitive parameters, improved control of which would increase reliability. Thus PFM should play a vital role in the analysis of steel pressure boundaries despite all shortcomings. (author). 19 refs, 7 figs, 1 tab

Existence of Positive Solutions to a Boundary Value Problem for a Delayed Nonlinear Fractional Differential System

Directory of Open Access Journals (Sweden)

Chen Yuming

2011-01-01

Full Text Available Though boundary value problems for fractional differential equations have been extensively studied, most of the studies focus on scalar equations and the fractional order between 1 and 2. On the other hand, delay is natural in practical systems. However, not much has been done for fractional differential equations with delays. Therefore, in this paper, we consider a boundary value problem of a general delayed nonlinear fractional system. With the help of some fixed point theorems and the properties of the Green function, we establish several sets of sufficient conditions on the existence of positive solutions. The obtained results extend and include some existing ones and are illustrated with some examples for their feasibility.

Iterative method for solving a problem with mixed boundary conditions for biharmonic equation arising in fracture mechanics

Directory of Open Access Journals (Sweden)

Dang Quang A

2013-02-01

Full Text Available In this paper we consider a mixed boundary value problem for biharmonic equation of the Airy stress function which models a crack problem of a solid elastic plate. An iterative method for reducing the problem to a sequence of mixed problems for Poisson equations is proposed and investigated. The convergence of the method is established theoretically and illustrated on many numerical experiments.

The Method of Subsuper Solutions for Weighted p(r-Laplacian Equation Boundary Value Problems

Directory of Open Access Journals (Sweden)

Zhimei Qiu

2008-10-01

Full Text Available This paper investigates the existence of solutions for weighted p(r-Laplacian ordinary boundary value problems. Our method is based on Leray-Schauder degree. As an application, we give the existence of weak solutions for p(x-Laplacian partial differential equations.

Numerical solution of sixth-order boundary-value problems using Legendre wavelet collocation method

Science.gov (United States)

Sohaib, Muhammad; Haq, Sirajul; Mukhtar, Safyan; Khan, Imad

2018-03-01

An efficient method is proposed to approximate sixth order boundary value problems. The proposed method is based on Legendre wavelet in which Legendre polynomial is used. The mechanism of the method is to use collocation points that converts the differential equation into a system of algebraic equations. For validation two test problems are discussed. The results obtained from proposed method are quite accurate, also close to exact solution, and other different methods. The proposed method is computationally more effective and leads to more accurate results as compared to other methods from literature.

Mixed problem with nonlocal boundary conditions for a third-order partial differential equation of mixed type

OpenAIRE

Denche, M.; Marhoune, A. L.

2001-01-01

We study a mixed problem with integral boundary conditions for a third-order partial differential equation of mixed type. We prove the existence and uniqueness of the solution. The proof is based on two-sided a priori estimates and on the density of the range of the operator generated by the considered problem.

Numerical treatment of the unsteady hydromagnetic thermal boundary layer problem

International Nuclear Information System (INIS)

Drymonitou, M.A.; Geroyannis, V.S.; Goudas, C.L.

1980-01-01

This paper presents a suitable numerical method for the treatment of the unsteady hydromagnetic thermal boundary layer problem for flows past an infinite porous flat plate, the motion of which is governed by a general time-dependent law, under the influence of a transverse externally set magnetic field. The normal velocity of suction/injection at the plate is also assumed to be time-dependent. The results obtained on the basis of numerical approximations seem to compare favourably with earlier results (Pande et al., 1976; Tokis, 1978). Analytical approximations are given for the cases of a plate (i) generally accelerated and (ii) harmonically oscillating. The direct numerical treatment is obviously advantageous since it allows handling of cases where the known methods for analytical approximations are not applicable. This problem is closely related to the motions and heat transfer occurring locally on the surfaces of stars. (orig.)

Modified quasi-boundary value method for Cauchy problems of elliptic equations with variable coefficients

Directory of Open Access Journals (Sweden)

Hongwu Zhang

2011-08-01

Full Text Available In this article, we study a Cauchy problem for an elliptic equation with variable coefficients. It is well-known that such a problem is severely ill-posed; i.e., the solution does not depend continuously on the Cauchy data. We propose a modified quasi-boundary value regularization method to solve it. Convergence estimates are established under two a priori assumptions on the exact solution. A numerical example is given to illustrate our proposed method.

Uniqueness in some higher order elliptic boundary value problems in n dimensional domains

Directory of Open Access Journals (Sweden)

C.-P. Danet

2011-07-01

Full Text Available We develop maximum principles for several P functions which are defined on solutions to equations of fourth and sixth order (including a equation which arises in plate theory and bending of cylindrical shells. As a consequence, we obtain uniqueness results for fourth and sixth order boundary value problems in arbitrary n dimensional domains.

Positive solutions for a nonlinear periodic boundary-value problem with a parameter

Directory of Open Access Journals (Sweden)

Jingliang Qiu

2012-08-01

Full Text Available Using topological degree theory with a partially ordered structure of space, sufficient conditions for the existence and multiplicity of positive solutions for a second-order nonlinear periodic boundary-value problem are established. Inspired by ideas in Guo and Lakshmikantham [6], we study the dependence of positive periodic solutions as a parameter approaches infinity, $$ lim_{lambdao +infty}|x_{lambda}|=+infty,quadhbox{or}quad lim_{lambdao+infty}|x_{lambda}|=0. $$

Positive solutions of three-point boundary-value problems for p-Laplacian singular differential equations

Directory of Open Access Journals (Sweden)

George N. Galanis

2005-10-01

Full Text Available In this paper we prove the existence of positive solutions for the three-point singular boundary-value problem$$ -[phi _{p}(u']'=q(tf(t,u(t,quad 0boundary-value problem remains away from the origin for the case where the nonlinearity is sublinear and so we avoid its singularity at $u=0$.

Semi-classical approximation and the problem of boundary conditions in the theory of relativistic particle radiation

International Nuclear Information System (INIS)

Akhiezer, A.I.; Shul'ga, N.F.

1991-01-01

The process of relativistic particle radiation in an external field has been studied in the semi-classical approximation rather extensively. The main problem arising in the studies is in expressing the formula of the quantum theory of radiation in terms of classical quantities, for example of the classical trajectories. However, it still remains unclear how the particle trajectory is assigned, that is which particular initial or boundary conditions determine the trajectory in semi-classical approximation quantum theory of radiation. We shall try to solve this problem. Its importance comes from the fact that in some cases one and the same boundary conditions may give rise to two or more trajectories. We demonstrate that this fact must necessarily be taken into account on deriving the classical limit for the formulae of the quantum theory of radiation, since it leads to a specific interference effect in radiation. The method we used to deal with the problem is similar to the method employed by Fock to analyze the problem of a canonical transformation in classical and quantum mechanics. (author)

Fictitious domain methods for elliptic problems with general boundary conditions with an application to the numerical simulation of two phase flows

International Nuclear Information System (INIS)

Ramiere, I.

2006-09-01

This work is dedicated to the introduction of two original fictitious domain methods for the resolution of elliptic problems (mainly convection-diffusion problems) with general and eventually mixed boundary conditions: Dirichlet, Robin or Neumann. The originality lies in the approximation of the immersed boundary by an approximate interface derived from the fictitious domain Cartesian mesh, which is generally not boundary-fitted to the physical domain. The same generic numerical scheme is used to impose the embedded boundary conditions. Hence, these methods require neither a surface mesh of the immersed boundary nor the local modification of the numerical scheme. We study two modelling of the immersed boundary. In the first one, called spread interface, the approximate immersed boundary is the union of the cells crossed by the physical immersed boundary. In the second one, called thin interface, the approximate immersed boundary lies on sides of mesh cells. Additional algebraic transmission conditions linking both flux and solution jumps through the thin approximate interface are introduced. The fictitious problem to solve as well as the treatment of the embedded boundary conditions are detailed for the two methods. A Q1 finite element scheme is implemented for the numerical validation of the spread interface approach while a new cell-centered finite volume scheme is derived for the thin interface approach with immersed jumps. Each method is then combined to multilevel local mesh refinement algorithms (with solution or flux residual) to increase the precision of the solution in the vicinity of the immersed interface. A convergence analysis of a Q1 finite element method with non-boundary fitted meshes is also presented. This study proves the convergence rates of the present methods. Among the various industrial applications, the simulation on a model of heat exchanger in french nuclear power plants enables us to appreciate the performances of the fictitious domain

Comparison report for CSNI International Standard Problem 12 (ROSA-III Run 912)

International Nuclear Information System (INIS)

Tasaka, Kanji; Anoda, Yoshinari; Kumamaru, Hiroshige; Nakamura, Hideo; Iriko, Masanori; Yonomoto, Taisuke; Shiba, Masayoshi

1982-09-01

ROSA-III Run 912 was identified as International Standard Problem 12 by the Committee on the Safety of Nuclear Installations. Run 912 simulated a 5% split break LOCA condition in a BWR at the pump suction in the recirculation line with the HPCS failure. Comparisons between the test data and the calculations by eight international participants were made and discussed. (author)

A moving boundary problem for the Stokes equations involving osmosis : Variational modelling and short-time well-posedness

NARCIS (Netherlands)

Lippoth, F.; Peletier, M.A.; Prokert, G.

2016-01-01

Within the framework of variational modelling we derive a one-phase moving boundary problem describing the motion of a semipermeable membrane enclosing a viscous liquid, driven by osmotic pressure and surface tension of the membrane. For this problem we prove the existence of classical solutions for

Numerical solution of singularity-perturbed two-point boundary-value problems

International Nuclear Information System (INIS)

Masenge, R.W.P.

1993-07-01

Physical processes which involve transportation of slowly diffusing substances in a fast-flowing medium are mathematically modelled by so-called singularly-perturbed second order convection diffusion differential equations in which the convective first order terms dominate over the diffusive second order terms. In general, analytical solutions of such equations are characterized by having sharp solution fronts in some sections of the interior and/or the boundary of the domain of solution. The presence of these (usually very narrow) layer regions in the solution domain makes the task of globally approximating such solutions by standard numerical techniques very difficult. In this expository paper we use a simple one-dimensional prototype problem as a vehicle for analysing the nature of the numerical approximation difficulties involved. In the sequel we present, without detailed derivation, two practical numerical schemes which succeed in varying degrees in numerically resolving the layer of the solution to the prototype problem. (author). 3 refs, 1 fig., 1 tab

The mixed boundary value problem, Krein resolvent formulas and spectral asymptotic estimates

DEFF Research Database (Denmark)

Grubb, Gerd

2011-01-01

For a second-order symmetric strongly elliptic operator A on a smooth bounded open set in Rn, the mixed problem is defined by a Neumann-type condition on a part Σ+ of the boundary and a Dirichlet condition on the other part Σ−. We show a Kreĭn resolvent formula, where the difference between its...... to the area of Σ+, in the case where A is principally equal to the Laplacian...

Uniqueness in the inverse boundary value problem for piecewise homogeneous anisotropic elasticity

OpenAIRE

Cârstea, Cătălin I.; Honda, Naofumi; Nakamura, Gen

2016-01-01

Consider a three dimensional piecewise homogeneous anisotropic elastic medium $\\Omega$ which is a bounded domain consisting of a finite number of bounded subdomains $D_\\alpha$, with each $D_\\alpha$ a homogeneous elastic medium. One typical example is a finite element model with elements with curvilinear interfaces for an ansiotropic elastic medium. Assuming the $D_\\alpha$ are known and Lipschitz, we are concerned with the uniqueness in the inverse boundary value problem of identifying the ani...

Inverse eigenvalue problems for Sturm-Liouville equations with spectral parameter linearly contained in one of the boundary conditions

OpenAIRE

Guliyev, Namig J.

2008-01-01

International audience; Inverse problems of recovering the coefficients of Sturm–Liouville problems with the eigenvalue parameter linearly contained in one of the boundary conditions are studied: 1) from the sequences of eigenvalues and norming constants; 2) from two spectra. Necessary and sufficient conditions for the solvability of these inverse problems are obtained.

The Sturm-Liouville inverse spectral problem with boundary conditions depending on the spectral parameter

Directory of Open Access Journals (Sweden)

Cornelis van der Mee

2005-01-01

Full Text Available We present the complete version including proofs of the results announced in [van der Mee C., Pivovarchik V.: A Sturm-Liouville spectral problem with boundary conditions depending on the spectral parameter. Funct. Anal. Appl. 36 (2002, 315–317 [Funkts. Anal. Prilozh. 36 (2002, 74–77 (Russian

A Comprehensive Review of Boundary Integral Formulations of Acoustic Scattering Problems

Directory of Open Access Journals (Sweden)

S.I. Zaman

2000-12-01

Full Text Available This is a review presenting an overview of the developments in boundary integral formulations of the acoustic scattering problems. Generally, the problem is formulated in one of two ways viz. Green’s representation formula, and the Layer-theoretic formulation utilizing either a simple-layer or a double-layer potential. The review presents and expounds the major contributions in this area over the last four decades. The need for a robust and improved formulation of the exterior scattering problem (Neumann or Dirichlet arose due to the fact that the classical formulation failed to yield a unique solution at (acoustic wave-numbers which correspond to eigenvalues (eigenfrequencies of the corresponding interior scattering problem. Moreover, this correlation becomes more pronounced as the wave-numbers become larger i.e. as the (acoustic frequency increases. The robust integral formulations which are discussed here yield Fredholms integral equations of the second kind which are more amenable to computation than the first kind. However, the integral equation involves a hypersingular kernel which creates ill-conditioning in the final matrix representation. This is circumvented by a regularisation technique. An extensive useful list of references is also presented here for researchers in this area.

Smooth and robust solutions for Dirichlet boundary control of fluid-solid conjugate heat transfer problems

KAUST Repository

Yan, Yan

2015-01-01

We study a new optimization scheme that generates smooth and robust solutions for Dirichlet velocity boundary control (DVBC) of conjugate heat transfer (CHT) processes. The solutions to the DVBC of the incompressible Navier-Stokes equations are typically nonsmooth, due to the regularity degradation of the boundary stress in the adjoint Navier-Stokes equations. This nonsmoothness is inherited by the solutions to the DVBC of CHT processes, since the CHT process couples the Navier-Stokes equations of fluid motion with the convection-diffusion equations of fluid-solid thermal interaction. Our objective in the CHT boundary control problem is to select optimally the fluid inflow profile that minimizes an objective function that involves the sum of the mismatch between the temperature distribution in the fluid system and a prescribed temperature profile and the cost of the control.Our strategy to resolve the nonsmoothness of the boundary control solution is based on two features, namely, the objective function with a regularization term on the gradient of the control profile on both the continuous and the discrete levels, and the optimization scheme with either explicit or implicit smoothing effects, such as the smoothed Steepest Descent and the Limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) methods. Our strategy to achieve the robustness of the solution process is based on combining the smoothed optimization scheme with the numerical continuation technique on the regularization parameters in the objective function. In the section of numerical studies, we present two suites of experiments. In the first one, we demonstrate the feasibility and effectiveness of our numerical schemes in recovering the boundary control profile of the standard case of a Poiseuille flow. In the second one, we illustrate the robustness of our optimization schemes via solving more challenging DVBC problems for both the channel flow and the flow past a square cylinder, which use initial

Inverse Boundary Value Problem for Non-linear Hyperbolic Partial Differential Equations

OpenAIRE

Nakamura, Gen; Vashisth, Manmohan

2017-01-01

In this article we are concerned with an inverse boundary value problem for a non-linear wave equation of divergence form with space dimension $n\\geq 3$. This non-linear wave equation has a trivial solution, i.e. zero solution. By linearizing this equation at the trivial solution, we have the usual linear isotropic wave equation with the speed $\\sqrt{\\gamma(x)}$ at each point $x$ in a given spacial domain. For any small solution $u=u(t,x)$ of this non-linear equation, we have the linear isotr...

Singular integral equations boundary problems of function theory and their application to mathematical physics

CERN Document Server

Muskhelishvili, N I

2011-01-01

Singular integral equations play important roles in physics and theoretical mechanics, particularly in the areas of elasticity, aerodynamics, and unsteady aerofoil theory. They are highly effective in solving boundary problems occurring in the theory of functions of a complex variable, potential theory, the theory of elasticity, and the theory of fluid mechanics.This high-level treatment by a noted mathematician considers one-dimensional singular integral equations involving Cauchy principal values. Its coverage includes such topics as the Hölder condition, Hilbert and Riemann-Hilbert problem

Theorems on differential inequalities and periodic boundary value problem for second-order ordinary differential equations

Czech Academy of Sciences Publication Activity Database

Lomtatidze, Alexander

2016-01-01

Roč. 67, č. 1 (2016), s. 1-129 ISSN 1512-0015 Institutional support: RVO:67985840 Keywords : periodic boundary value problem * positive solution * singular equation Subject RIV: BA - General Mathematics http://rmi.tsu.ge/jeomj/memoirs/vol67/abs67-1.htm

Elasticity problems in domains with nonsmooth boundaries

International Nuclear Information System (INIS)

Esparza, David

2001-01-01

In the present work we study the behaviour of elastic stress fields in domains with non-regular boundaries. We consider three-dimensional problems in elastic media with thin conical defects (inclusions or cavities) and analyse the stress singularity at their vertices. To construct asymptotic expansions for the stress and displacement fields in terms of a small parameter ε related to the 'thickness' of the defect, we employ a technique based on the work by Kondrat'ev, Maz'ya, Nazarov and Plamenevskii. We first study the stress distribution in an elastic body with a thin conical notch. We derive an asymptotic representation for the stress singularity exponent by reducing the original problem to a spectral problem for a 9x9 matrix. The elements of this matrix are found to depend upon the geometry of the cross-section of the notch and the elastic properties of the medium. We specify the sets of eigenvalues and the corresponding eigenvectors for a circular, elliptical, 'triangular' and 'square' cross-section, and show that the strongest singularity is associated with the 'triangular' cross-section, and is generated by a non-axisymmetric load. We then analyse the stress distribution near a thin conical inclusion which is allowed to slide freely along its axis. We derive the representation for the stress singularity exponent for the case of a circular conical inclusion whose elastic properties differ from those of the medium. In the last chapter we study the stress distribution in the vicinity of a thin 'coated' conical inclusion. We show that a soft thin coating (perfectly bonded to the inclusion and the surrounding material) can be replaced by a so-called linear interface at which the normal displacement is discontinuous, and the stresses are proportional to the 'jump' in the normal displacement across the coating. We analyse the effect of the properties of the coating on the stress singularity exponent and compare the results with those for a perfectly bonded

A Hartman–Nagumo inequality for the vector ordinary -Laplacian and applications to nonlinear boundary value problems

Directory of Open Access Journals (Sweden)

Ureña Antonio J

2002-01-01

Full Text Available A generalization of the well-known Hartman–Nagumo inequality to the case of the vector ordinary -Laplacian and classical degree theory provide existence results for some associated nonlinear boundary value problems.

Existence of positive solutions for nonlocal second-order boundary value problem with variable parameter in Banach spaces

Directory of Open Access Journals (Sweden)

Zhang Peiguo

2011-01-01

Full Text Available Abstract By obtaining intervals of the parameter λ, this article investigates the existence of a positive solution for a class of nonlinear boundary value problems of second-order differential equations with integral boundary conditions in abstract spaces. The arguments are based upon a specially constructed cone and the fixed point theory in cone for a strict set contraction operator. MSC: 34B15; 34B16.

Partial differential equations & boundary value problems with Maple

CERN Document Server

Articolo, George A

2009-01-01

Partial Differential Equations and Boundary Value Problems with Maple presents all of the material normally covered in a standard course on partial differential equations, while focusing on the natural union between this material and the powerful computational software, Maple. The Maple commands are so intuitive and easy to learn, students can learn what they need to know about the software in a matter of hours- an investment that provides substantial returns. Maple''s animation capabilities allow students and practitioners to see real-time displays of the solutions of partial differential equations.  Maple files can be found on the books website. Ancillary list: Maple files- http://www.elsevierdirect.com/companion.jsp?ISBN=9780123747327  Provides a quick overview of the software w/simple commands needed to get startedIncludes review material on linear algebra and Ordinary Differential equations, and their contribution in solving partial differential equationsIncorporates an early introduction to Sturm-L...

The focal boundary value problem for strongly singular higher-order nonlinear functional-differential equations

Czech Academy of Sciences Publication Activity Database

Mukhigulashvili, Sulkhan; Půža, B.

2015-01-01

Roč. 2015, January (2015), s. 17 ISSN 1687-2770 Institutional support: RVO:67985840 Keywords : higher order nonlinear functional-differential equations * two-point right-focal boundary value problem * strong singularity Subject RIV: BA - General Mathematics Impact factor: 0.642, year: 2015 http://link.springer.com/article/10.1186%2Fs13661-014-0277-1

Solvability of fractional multi-point boundary-value problems with p-Laplacian operator at resonance

Directory of Open Access Journals (Sweden)

Tengfei Shen

2014-02-01

Full Text Available In this article, we consider the multi-point boundary-value problem for nonlinear fractional differential equations with $p$-Laplacian operator: $$\\displaylines{ D_{0^+}^\\beta \\varphi_p (D_{0^+}^\\alpha u(t = f(t,u(t,D_{0^+}^{\\alpha - 2} u(t,D_{0^+}^{\\alpha - 1} u(t, D_{0^+}^\\alpha u(t,\\quad t \\in (0,1, \\cr u(0 = u'(0=D_{0^+}^\\alpha u(0 = 0,\\quad D_{0^+}^{\\alpha - 1} u(1 = \\sum_{i = 1}^m {\\sigma_i D_{0^+}^{\\alpha - 1} u(\\eta_i } , }$$ where $2 < \\alpha \\le 3$, $0 < \\beta \\le 1$, $3 < \\alpha + \\beta \\le 4$, $\\sum_{i = 1}^m {\\sigma_i } = 1$, $D_{0^+}^\\alpha$ is the standard Riemann-Liouville fractional derivative. $\\varphi_{p}(s=|s|^{p-2}s$ is p-Laplacians operator. The existence of solutions for above fractional boundary value problem is obtained by using the extension of Mawhin's continuation theorem due to Ge, which enrich konwn results. An example is given to illustrate the main result.

Fermat collocation method for the solutions of nonlinear system of second order boundary value problems

Directory of Open Access Journals (Sweden)

Salih Yalcinbas

2016-01-01

Full Text Available In this study, a numerical approach is proposed to obtain approximate solutions of nonlinear system of second order boundary value problem. This technique is essentially based on the truncated Fermat series and its matrix representations with collocation points. Using the matrix method, we reduce the problem system of nonlinear algebraic equations. Numerical examples are also given to demonstrate the validity and applicability of the presented technique. The method is easy to implement and produces accurate results.

Positive Solutions of Three-Order Delayed Periodic Boundary Value Problems

Directory of Open Access Journals (Sweden)

Na Wang

2017-01-01

Full Text Available Our main purpose is to consider the existence of positive solutions for three-order two-point boundary value problem in the following form: u′′′(t+ρ3u(t=f(t,u(t-τ,  0≤t≤2π,  u(i(0=u(i(2π,  i=1,2,  u(t=σ,  -τ≤t≤0, where σ,ρ, and τ are given constants satisfying τ∈(0,π/2. Some inequality conditions on ρ3u-f(t,u guaranteeing the existence and nonexistence of positive solutions are presented. Our discussion is based on the fixed point theorem in cones.

Two-phase semilinear free boundary problem with a degenerate phase

KAUST Repository

Matevosyan, Norayr

2010-10-16

We study minimizers of the energy functional ∫D[{pipe}∇u{pipe}2 + λ(u+)p]dx for p ∈ (0, 1) without any sign restriction on the function u. The distinguished feature of the problem is the lack of nondegeneracy in the negative phase. The main result states that in dimension two the free boundaries Γ+ = ∂{u > 0} ∩ D andΓ- = ∂{u < 0} ∩ D are C1,α-regular, provided 1 - ∈0 < p < 1. The proof is obtained by a careful iteration of the Harnack inequality to obtain a nontrivial growth estimate in the negative phase, compensating for the apriori unknown nondegeneracy. © 2010 Springer-Verlag.

Oblique derivative problems for generalized Rassias equations of mixed type with several characteristic boundaries

Directory of Open Access Journals (Sweden)

Guo Chun Wen

2009-05-01

Full Text Available This article concerns the oblique derivative problems for second-order quasilinear degenerate equations of mixed type with several characteristic boundaries, which include the Tricomi problem as a special case. First we formulate the problem and obtain estimates of its solutions, then we show the existence of solutions by the successive iterations and the Leray-Schauder theorem. We use a complex analytic method: elliptic complex functions are used in the elliptic domain, and hyperbolic complex functions in the hyperbolic domain, such that second-order equations of mixed type with degenerate curve are reduced to the first order mixed complex equations with singular coefficients. An application of the complex analytic method, solves (1.1 below with $m=n=1$, $a=b=0$, which was posed as an open problem by Rassias.

Application of power series to the solution of the boundary value problem for a second order nonlinear differential equation

International Nuclear Information System (INIS)

Semenova, V.N.

2016-01-01

A boundary value problem for a nonlinear second order differential equation has been considered. A numerical method has been proposed to solve this problem using power series. Results of numerical experiments have been presented in the paper [ru

Monotone methods for solving a boundary value problem of second order discrete system

Directory of Open Access Journals (Sweden)

Wang Yuan-Ming

1999-01-01

Full Text Available A new concept of a pair of upper and lower solutions is introduced for a boundary value problem of second order discrete system. A comparison result is given. An existence theorem for a solution is established in terms of upper and lower solutions. A monotone iterative scheme is proposed, and the monotone convergence rate of the iteration is compared and analyzed. The numerical results are given.

Existence and Estimates of Positive Solutions for Some Singular Fractional Boundary Value Problems

Directory of Open Access Journals (Sweden)

Habib Mâagli

2014-01-01

fractional boundary value problem:Dαu(x=−a(xuσ(x, x∈(0,1 with the conditions limx→0+⁡x2−αu(x=0, u(1=0, where 1<α≤2, σ∈(−1,1, and a is a nonnegative continuous function on (0,1 that may be singular at x=0 or x=1. We also give the global behavior of such a solution.

Response of Non-Linear Shock Absorbers-Boundary Value Problem Analysis

Science.gov (United States)

Rahman, M. A.; Ahmed, U.; Uddin, M. S.

2013-08-01

A nonlinear boundary value problem of two degrees-of-freedom (DOF) untuned vibration damper systems using nonlinear springs and dampers has been numerically studied. As far as untuned damper is concerned, sixteen different combinations of linear and nonlinear springs and dampers have been comprehensively analyzed taking into account transient terms. For different cases, a comparative study is made for response versus time for different spring and damper types at three important frequency ratios: one at r = 1, one at r > 1 and one at r <1. The response of the system is changed because of the spring and damper nonlinearities; the change is different for different cases. Accordingly, an initially stable absorber may become unstable with time and vice versa. The analysis also shows that higher nonlinearity terms make the system more unstable. Numerical simulation includes transient vibrations. Although problems are much more complicated compared to those for a tuned absorber, a comparison of the results generated by the present numerical scheme with the exact one shows quite a reasonable agreement

Communication: modeling charge-sign asymmetric solvation free energies with nonlinear boundary conditions.

Science.gov (United States)

Bardhan, Jaydeep P; Knepley, Matthew G

2014-10-07

We show that charge-sign-dependent asymmetric hydration can be modeled accurately using linear Poisson theory after replacing the standard electric-displacement boundary condition with a simple nonlinear boundary condition. Using a single multiplicative scaling factor to determine atomic radii from molecular dynamics Lennard-Jones parameters, the new model accurately reproduces MD free-energy calculations of hydration asymmetries for: (i) monatomic ions, (ii) titratable amino acids in both their protonated and unprotonated states, and (iii) the Mobley "bracelet" and "rod" test problems [D. L. Mobley, A. E. Barber II, C. J. Fennell, and K. A. Dill, "Charge asymmetries in hydration of polar solutes," J. Phys. Chem. B 112, 2405-2414 (2008)]. Remarkably, the model also justifies the use of linear response expressions for charging free energies. Our boundary-element method implementation demonstrates the ease with which other continuum-electrostatic solvers can be extended to include asymmetry.

Communication: Modeling charge-sign asymmetric solvation free energies with nonlinear boundary conditions

International Nuclear Information System (INIS)

Bardhan, Jaydeep P.; Knepley, Matthew G.

2014-01-01

We show that charge-sign-dependent asymmetric hydration can be modeled accurately using linear Poisson theory after replacing the standard electric-displacement boundary condition with a simple nonlinear boundary condition. Using a single multiplicative scaling factor to determine atomic radii from molecular dynamics Lennard-Jones parameters, the new model accurately reproduces MD free-energy calculations of hydration asymmetries for: (i) monatomic ions, (ii) titratable amino acids in both their protonated and unprotonated states, and (iii) the Mobley “bracelet” and “rod” test problems [D. L. Mobley, A. E. Barber II, C. J. Fennell, and K. A. Dill, “Charge asymmetries in hydration of polar solutes,” J. Phys. Chem. B 112, 2405–2414 (2008)]. Remarkably, the model also justifies the use of linear response expressions for charging free energies. Our boundary-element method implementation demonstrates the ease with which other continuum-electrostatic solvers can be extended to include asymmetry

Communication: Modeling charge-sign asymmetric solvation free energies with nonlinear boundary conditions

Energy Technology Data Exchange (ETDEWEB)

Bardhan, Jaydeep P. [Department of Mechanical and Industrial Engineering, Northeastern University, Boston, Massachusetts 02115 (United States); Knepley, Matthew G. [Computation Institute, The University of Chicago, Chicago, Illinois 60637 (United States)

2014-10-07

We show that charge-sign-dependent asymmetric hydration can be modeled accurately using linear Poisson theory after replacing the standard electric-displacement boundary condition with a simple nonlinear boundary condition. Using a single multiplicative scaling factor to determine atomic radii from molecular dynamics Lennard-Jones parameters, the new model accurately reproduces MD free-energy calculations of hydration asymmetries for: (i) monatomic ions, (ii) titratable amino acids in both their protonated and unprotonated states, and (iii) the Mobley “bracelet” and “rod” test problems [D. L. Mobley, A. E. Barber II, C. J. Fennell, and K. A. Dill, “Charge asymmetries in hydration of polar solutes,” J. Phys. Chem. B 112, 2405–2414 (2008)]. Remarkably, the model also justifies the use of linear response expressions for charging free energies. Our boundary-element method implementation demonstrates the ease with which other continuum-electrostatic solvers can be extended to include asymmetry.

Closed form solution to a second order boundary value problem and its application in fluid mechanics

International Nuclear Information System (INIS)

Eldabe, N.T.; Elghazy, E.M.; Ebaid, A.

2007-01-01

The Adomian decomposition method is used by many researchers to investigate several scientific models. In this Letter, the modified Adomian decomposition method is applied to construct a closed form solution for a second order boundary value problem with singularity

On sign constant solutions of certain boundary value problems for second-order functional differential equations

Czech Academy of Sciences Publication Activity Database

Lomtatidze, Alexander; Vodstrčil, Petr

2005-01-01

Roč. 84, č. 2 (2005), s. 197-209 ISSN 0003-6811 Institutional research plan: CEZ:AV0Z10190503 Keywords : second order linear functional differential equations * nonnegative solution * two-point boundary value problem Subject RIV: BA - General Mathematics http://www.tandfonline.com/doi/full/10.1080/00036810410001724427

The Initial and Neumann Boundary Value Problem for a Class Parabolic Monge-Ampère Equation

Directory of Open Access Journals (Sweden)

Juan Wang

2013-01-01

Full Text Available We consider the existence, uniqueness, and asymptotic behavior of a classical solution to the initial and Neumann boundary value problem for a class nonlinear parabolic equation of Monge-Ampère type. We show that such solution exists for all times and is unique. It converges eventually to a solution that satisfies a Neumann type problem for nonlinear elliptic equation of Monge-Ampère type.

Positive solutions of boundary value problem for singular positone and semi-positone third-order difference equations

Directory of Open Access Journals (Sweden)

Gai Gongqi

2011-01-01

Full Text Available Abstract This article studies the boundary value problems for the third-order nonlinear singular difference equations Δ 3 u ( i - 2 + λ a ( i f ( i , u ( i = 0 , i ∈ [ 2 , T + 2 ] , satisfying five kinds of different boundary value conditions. This article shows the existence of positive solutions for positone and semi-positone type. The nonlinear term may be singular. Two examples are also given to illustrate the main results. The arguments are based upon fixed point theorems in a cone. MSC [2008]: 34B15; 39A10.

Quasisolutions of Inverse Boundary-Value Problem of Aerodynamics for Dense Airfoil Grids

Directory of Open Access Journals (Sweden)

A.M. Elizarov

2016-12-01

Full Text Available In the process of turbomachinery development, it is of great importance to accurately design impellers and select their blade shape. One of the promising approaches to solving this problem is based on the theory of inverse boundary-value problems in aerodynamics. It helps to develop methods for profiling airfoil grids with predetermined properties in the same way as it is done for isolated airfoils. In this paper, methods have been worked out to find quasisolutions of the inverse boundary-value problem in aerodynamics for a plane airfoil grid. Two methods of quasisolution have been described. The first “`formal” method is similar, in its essence, to the method used for construction of quasisolution for an isolated airfoil. It has been shown that such quasisolutions provide satisfactory results for grids having a sufficiently large relative airfoil pitch. If pitch values are low, this method is unacceptable, because “modified” velocity distribution in some areas is significantly different from the original one in this case. For this reason, areas with significant changes in the angle of the tangent line appear in the airfoil contour and the flow region becomes multivalent. To satisfy the conditions of solvability in the case of grids having a small airfoil pitch, a new quasisolution construction method taking into account the specifics of the problem has been suggested. The desired effect has been achieved due to changes in the weighting function of the minimized functional. The comparison of the results of construction of the new quasisolution with the results obtained by the “formal” method has demonstrated that the constructed airfoils are very similar when the pitch is large. In the case of dense grids, it is clear that preference should be given to the second method, as it brings less distortion to the initial velocity distribution and, thus, allows to physically find an actual airfoil contour.

Local non-similarity method through the Crocco's transformation in boundary layer problem

International Nuclear Information System (INIS)

Jardim, R.G.M.

1981-04-01

The coordinate transformation developed by L. Crocco to obtain the solution of the compressible fluid flows over isotermal flat plates is originally employed in the present work, with the purpose of adding its inherent advantage to the Non-Similarity Method idealized by E.M. Sparrow, in the solution of the incompressible non-similar boundary layers. The Crocco's transformation is applied to the conservation equation for forced convection, laminar, constant properties and two-dimensional flows over solids. Two non-similar problems arisen from freestream velocity distribution, the cylinder in crossflow and the Howarth's retarded flow, are solved with a view to illustrating the new procedure. In those solutions the effect of frictional heat is also considered. The results of hydrodynamic and thermal problems are compared with available published information and good agreement was observed. (Author) [pt

Multiple solutions of a free-boundary FRC equilibrium problem in a metal cylinder

International Nuclear Information System (INIS)

Spencer, R.L.; Hewett, D.W.

1981-01-01

A new approach to the computation of FRC equilibria that avoids previously encountered difficulties is presented. For arbitrary pressure profiles it is computationally expensive, but for one special pressure profile the problem is simple enough to require only minutes of Cray time; it is this problem that we have solved. We solve the Grad-Shafranov equation, Δ/sup */psi = r 2 p'(psi), in an infinitely long flux conserving cylinder of radius a with the boundary conditions that psi(a,z) = -psi/sub w/ and that delta psi/delta z = 0 as [z] approaches infinity. The pressure profile is p'(psi) = cH(psi) where c is a constant and where H(x) is the Heaviside function. We have found four solutions to this problem: There is a purely vacuum state, two z-independent plasma solutions, and an r-z-dependent plasma state

Classical boundary-value problem in Riemannian quantum gravity and self-dual Taub-NUT-(anti)de Sitter geometries

International Nuclear Information System (INIS)

Akbar, M.M.; D'Eath, P.D.

2003-01-01

The classical boundary-value problem of the Einstein field equations is studied with an arbitrary cosmological constant, in the case of a compact (S 3 ) boundary given a biaxial Bianchi-IX positive-definite three-metric, specified by two radii (a,b). For the simplest, four-ball, topology of the manifold with this boundary, the regular classical solutions are found within the family of Taub-NUT-(anti)de Sitter metrics with self-dual Weyl curvature. For arbitrary choice of positive radii (a,b), we find that there are three solutions for the infilling geometry of this type. We obtain exact solutions for them and for their Euclidean actions. The case of negative cosmological constant is investigated further. For reasonable squashing of the three-sphere, all three infilling solutions have real-valued actions which possess a 'cusp catastrophe' structure with a non-self-intersecting 'catastrophe manifold' implying that the dominant contribution comes from the unique real positive-definite solution on the ball. The positive-definite solution exists even for larger deformations of the three-sphere, as long as a certain inequality between a and b holds. The action of this solution is proportional to -a 3 for large a (∼b) and hence larger radii are favoured. The same boundary-value problem with more complicated interior topology containing a 'bolt' is investigated in a forthcoming paper

Optimal boundary control and boundary stabilization of hyperbolic systems

CERN Document Server

Gugat, Martin

2015-01-01

This brief considers recent results on optimal control and stabilization of systems governed by hyperbolic partial differential equations, specifically those in which the control action takes place at the boundary.  The wave equation is used as a typical example of a linear system, through which the author explores initial boundary value problems, concepts of exact controllability, optimal exact control, and boundary stabilization.  Nonlinear systems are also covered, with the Korteweg-de Vries and Burgers Equations serving as standard examples.  To keep the presentation as accessible as possible, the author uses the case of a system with a state that is defined on a finite space interval, so that there are only two boundary points where the system can be controlled.  Graduate and post-graduate students as well as researchers in the field will find this to be an accessible introduction to problems of optimal control and stabilization.

A boundary value problem for a third order hyperbolic equation with degeneration of order inside the domain

Directory of Open Access Journals (Sweden)

Ruzanna Kh. Makaova

2017-12-01

Full Text Available In this paper we study the boundary value problem for a degenerating third order equation of hyperbolic type in a mixed domain. The equation under consideration in the positive part of the domain coincides with the Hallaire equation, which is a pseudoparabolic type equation. Moreover, in the negative part of the domain it coincides with a degenerating hyperbolic equation of the first kind, the particular case of the Bitsadze–Lykov equation. The existence and uniqueness theorem for the solution is proved. The uniqueness of the solution to the problem is proved with the Tricomi method. Using the functional relationships of the positive and negative parts of the domain on the degeneration line, we arrive at the convolution type Volterra integral equation of the 2nd kind with respect to the desired solution by a derivative trace. With the Laplace transform method, we obtain the solution of the integral equation in its explicit form. At last, the solution to the problem under study is written out explicitly as the solution of the second boundary-value problem in the positive part of the domain for the Hallaire equation and as the solution to the Cauchy problem in the negative part of the domain for a degenerate hyperbolic equation of the first kind.

BOUNDARY VALUE PROBLEM FOR A LOADED EQUATION ELLIPTIC-HYPERBOLIC TYPE IN A DOUBLY CONNECTED DOMAIN

Directory of Open Access Journals (Sweden)

O.Kh. Abdullaev

2014-06-01

Full Text Available We study the existence and uniqueness of the solution of one boundary value problem for the loaded elliptic-hyperbolic equation of the second order with two lines of change of type in double-connected domain. Similar results have been received by D.M.Kuryhazov, when investigated domain is one-connected.

A Hierarchical Bayesian Setting for an Inverse Problem in Linear Parabolic PDEs with Noisy Boundary Conditions

KAUST Repository

Ruggeri, Fabrizio

2016-05-12

In this work we develop a Bayesian setting to infer unknown parameters in initial-boundary value problems related to linear parabolic partial differential equations. We realistically assume that the boundary data are noisy, for a given prescribed initial condition. We show how to derive the joint likelihood function for the forward problem, given some measurements of the solution field subject to Gaussian noise. Given Gaussian priors for the time-dependent Dirichlet boundary values, we analytically marginalize the joint likelihood using the linearity of the equation. Our hierarchical Bayesian approach is fully implemented in an example that involves the heat equation. In this example, the thermal diffusivity is the unknown parameter. We assume that the thermal diffusivity parameter can be modeled a priori through a lognormal random variable or by means of a space-dependent stationary lognormal random field. Synthetic data are used to test the inference. We exploit the behavior of the non-normalized log posterior distribution of the thermal diffusivity. Then, we use the Laplace method to obtain an approximated Gaussian posterior and therefore avoid costly Markov Chain Monte Carlo computations. Expected information gains and predictive posterior densities for observable quantities are numerically estimated using Laplace approximation for different experimental setups.

Extension Theory and Krein-type Resolvent Formulas for Nonsmooth Boundary Value Problems

DEFF Research Database (Denmark)

Abels, Helmut; Grubb, Gerd; Wood, Ian Geoffrey

2014-01-01

The theory of selfadjoint extensions of symmetric operators, and more generally the theory of extensions of dual pairs, was implemented some years ago for boundary value problems for elliptic operators on smooth bounded domains. Recently, the questions have been taken up again for nonsmooth domains....... In the present work we show that pseudodifferential methods can be used to obtain a full characterization, including Kreĭn resolvent formulas, of the realizations of nonselfadjoint second-order operators on

Hopf bifurcation of a free boundary problem modeling tumor growth with two time delays

International Nuclear Information System (INIS)

Xu Shihe

2009-01-01

In this paper, a free boundary problem modeling tumor growth with two discrete delays is studied. The delays respectively represents the time taken for cells to undergo mitosis and the time taken for the cell to modify the rate of cell loss due to apoptosis. We show the influence of time delays on the Hopf bifurcation when one of delays as a bifurcation parameter.

Multiple and sign-changing solutions for discrete Robin boundary value problem with parameter dependence

Directory of Open Access Journals (Sweden)

Long Yuhua

2017-12-01

Full Text Available In this paper, we study second-order nonlinear discrete Robin boundary value problem with parameter dependence. Applying invariant sets of descending flow and variational methods, we establish some new sufficient conditions on the existence of sign-changing solutions, positive solutions and negative solutions of the system when the parameter belongs to appropriate intervals. In addition, an example is given to illustrate our results.

Sufficient condition for existence of solutions for higher-order resonance boundary value problem with one-dimensional p-Laplacian

Directory of Open Access Journals (Sweden)

Liu Yang

2007-10-01

Full Text Available By using coincidence degree theory of Mawhin, existence results for some higher order resonance multipoint boundary value problems with one dimensional p-Laplacian operator are obtained.

Maximum Principles and Boundary Value Problems for First-Order Neutral Functional Differential Equations

Directory of Open Access Journals (Sweden)

Domoshnitsky Alexander

2009-01-01

Full Text Available We obtain the maximum principles for the first-order neutral functional differential equation where , and are linear continuous operators, and are positive operators, is the space of continuous functions, and is the space of essentially bounded functions defined on . New tests on positivity of the Cauchy function and its derivative are proposed. Results on existence and uniqueness of solutions for various boundary value problems are obtained on the basis of the maximum principles.

Existence of positive solutions for boundary value problems of fractional functional differential equations

Directory of Open Access Journals (Sweden)

Chuanzhi Bai

2010-06-01

Full Text Available This paper deals with the existence of positive solutions for a boundary value problem involving a nonlinear functional differential equation of fractional order $\\alpha$ given by $ D^{\\alpha} u(t + f(t, u_t = 0$, $t \\in (0, 1$, $2 < \\alpha \\le 3$, $ u^{\\prime}(0 = 0$, $u^{\\prime}(1 = b u^{\\prime}(\\eta$, $u_0 = \\phi$. Our results are based on the nonlinear alternative of Leray-Schauder type and Krasnosel'skii fixed point theorem.

Initial-boundary-value problem of the self-gravitating scalar field in the Bondi-Sachs gauge

International Nuclear Information System (INIS)

Frittelli, Simonetta; Gomez, Roberto

2007-01-01

It is shown that, in the Bondi-Sachs gauge that fixes the speed of incoming light rays to the value 1, the Einstein equations coupled to a scalar field in spherical symmetry are cast into a symmetric-hyperbolic system of equations for the scalar field, lapse and shift as fundamental variables. In this system of equations, the lapse and shift are incoming characteristic fields, and the scalar field has three components: incoming, outgoing and static. A constraint-preserving boundary condition is prescribed by imposing the projection of the Einstein equation normal to the boundary at the outer value of the radial coordinate. The boundary condition specifies one of the two incoming metric fields. The remaining incoming metric field and the incoming scalar field component need to be specified arbitrarily. Numerical simulations of the scattering of the scalar field by a black hole in the nonlinear regime are presented that illustrate interesting facts about black-hole physics and the behavior of the characteristic variables of the problem

Solvability conditions for non-local boundary value problems for two-dimensional half-linear differential systems

Czech Academy of Sciences Publication Activity Database

Kiguradze, I.; Šremr, Jiří

2011-01-01

Roč. 74, č. 17 (2011), s. 6537-6552 ISSN 0362-546X Institutional research plan: CEZ:AV0Z10190503 Keywords : half-linear differential system * non-local boundary value problem * solvability Subject RIV: BA - General Mathematics Impact factor: 1.536, year: 2011 http://www.sciencedirect.com/science/article/pii/S0362546X11004573

Communication: Modeling charge-sign asymmetric solvation free energies with nonlinear boundary conditions

Science.gov (United States)

Bardhan, Jaydeep P.; Knepley, Matthew G.

2014-01-01

We show that charge-sign-dependent asymmetric hydration can be modeled accurately using linear Poisson theory after replacing the standard electric-displacement boundary condition with a simple nonlinear boundary condition. Using a single multiplicative scaling factor to determine atomic radii from molecular dynamics Lennard-Jones parameters, the new model accurately reproduces MD free-energy calculations of hydration asymmetries for: (i) monatomic ions, (ii) titratable amino acids in both their protonated and unprotonated states, and (iii) the Mobley “bracelet” and “rod” test problems [D. L. Mobley, A. E. Barber II, C. J. Fennell, and K. A. Dill, “Charge asymmetries in hydration of polar solutes,” J. Phys. Chem. B 112, 2405–2414 (2008)]. Remarkably, the model also justifies the use of linear response expressions for charging free energies. Our boundary-element method implementation demonstrates the ease with which other continuum-electrostatic solvers can be extended to include asymmetry. PMID:25296776

Some free boundary problems in potential flow regime usinga based level set method

Energy Technology Data Exchange (ETDEWEB)

Garzon, M.; Bobillo-Ares, N.; Sethian, J.A.

2008-12-09

Recent advances in the field of fluid mechanics with moving fronts are linked to the use of Level Set Methods, a versatile mathematical technique to follow free boundaries which undergo topological changes. A challenging class of problems in this context are those related to the solution of a partial differential equation posed on a moving domain, in which the boundary condition for the PDE solver has to be obtained from a partial differential equation defined on the front. This is the case of potential flow models with moving boundaries. Moreover the fluid front will possibly be carrying some material substance which will diffuse in the front and be advected by the front velocity, as for example the use of surfactants to lower surface tension. We present a Level Set based methodology to embed this partial differential equations defined on the front in a complete Eulerian framework, fully avoiding the tracking of fluid particles and its known limitations. To show the advantages of this approach in the field of Fluid Mechanics we present in this work one particular application: the numerical approximation of a potential flow model to simulate the evolution and breaking of a solitary wave propagating over a slopping bottom and compare the level set based algorithm with previous front tracking models.

On inverse and direct free boundary problems in the theory of plasma equilibrium in a Tokamak

International Nuclear Information System (INIS)

Demidov, A.; Petrova, V.; Silantiev, V.

1996-01-01

Theorems of existence of simply connected 'plasma' domain for the cylindrical case of the Grad-Shafranov equation Δu = F(u) are given. For the inverse problem upper and lower estimates of normal derivative of u on the boundary of the 'plasma' domain are obtained. (author)

Triple solutions for a Dirichlet boundary value problem involving a perturbed discrete p(k-Laplacian operator

Directory of Open Access Journals (Sweden)

Khaleghi Moghadam Mohsen

2017-08-01

Full Text Available Triple solutions are obtained for a discrete problem involving a nonlinearly perturbed one-dimensional p(k-Laplacian operator and satisfying Dirichlet boundary conditions. The methods for existence rely on a Ricceri-local minimum theorem for differentiable functionals. Several examples are included to illustrate the main results.

Boundary fluxes for nonlocal diffusion

Science.gov (United States)

Cortazar, Carmen; Elgueta, Manuel; Rossi, Julio D.; Wolanski, Noemi

We study a nonlocal diffusion operator in a bounded smooth domain prescribing the flux through the boundary. This problem may be seen as a generalization of the usual Neumann problem for the heat equation. First, we prove existence, uniqueness and a comparison principle. Next, we study the behavior of solutions for some prescribed boundary data including blowing up ones. Finally, we look at a nonlinear flux boundary condition.

Existence of solutions for a fourth-order boundary value problem on the half-line via critical point theory

Directory of Open Access Journals (Sweden)

Mabrouk Briki

2016-05-01

Full Text Available In this paper, a fourth-order boundary value problem on the half-line is considered and existence of solutions is proved using a minimization principle and the mountain pass theorem.

Positive solutions for second-order boundary-value problems with phi-Laplacian

Directory of Open Access Journals (Sweden)

Diana-Raluca Herlea

2016-02-01

Full Text Available This article concerns the existence, localization and multiplicity of positive solutions for the boundary-value problem $$\\displaylines{ \\big(\\phi(u' \\big '+f(t,u =0, \\cr u(0 - a u'(0 = u'(1= 0, }$$ where $f:[0,1]\\times \\mathbb{R}_{+}\\to \\mathbb{R}_{+}$ is a continuous function and $\\phi :\\mathbb{R}\\to (-b,b$ is an increasing homeomorphism with $\\phi (0=0$. We obtain existence, localization and multiplicity results of positive solutions using Krasnosel'skii fixed point theorem in cones, and a weak Harnack type inequality. Concerning systems, the localization is established by the vector version of Krasnosel'skii theorem, where the compression-expansion conditions are expressed on components.

A symmetric solution of a multipoint boundary value problem at resonance

Directory of Open Access Journals (Sweden)

2006-01-01

Full Text Available We apply a coincidence degree theorem of Mawhin to show the existence of at least one symmetric solution of the nonlinear second-order multipoint boundary value problem u ″ ( t = f ( t , u ( t , | u ′ ( t | , t ∈ ( 0 , 1 , u ( 0 = ∑ i = 1 n μ i u ( ξ i , u ( 1 − t = u ( t , t ∈ ( 0 , 1 ] , where 0 < ξ 1 < ξ 2 < … ≤ ξ n 1 / 2 , ∑ i = 1 n μ i = 1 , f : [ 0 , 1 ] × ℝ 2 → ℝ with f ( t , x , y = f ( 1 − t , x , y , ( t , x , y ∈ [ 0 , 1 ] × ℝ 2 , satisfying the Carathéodory conditions.

Adaptive wavelet collocation methods for initial value boundary problems of nonlinear PDE's

Science.gov (United States)

Cai, Wei; Wang, Jian-Zhong

1993-01-01

We have designed a cubic spline wavelet decomposition for the Sobolev space H(sup 2)(sub 0)(I) where I is a bounded interval. Based on a special 'point-wise orthogonality' of the wavelet basis functions, a fast Discrete Wavelet Transform (DWT) is constructed. This DWT transform will map discrete samples of a function to its wavelet expansion coefficients in O(N log N) operations. Using this transform, we propose a collocation method for the initial value boundary problem of nonlinear PDE's. Then, we test the efficiency of the DWT transform and apply the collocation method to solve linear and nonlinear PDE's.

Existence of global solutions to free boundary value problems for bipolar Navier-Stokes-Possion systems

Directory of Open Access Journals (Sweden)

Jian Liu

2013-09-01

Full Text Available In this article, we consider the free boundary value problem for one-dimensional compressible bipolar Navier-Stokes-Possion (BNSP equations with density-dependent viscosities. For general initial data with finite energy and the density connecting with vacuum continuously, we prove the global existence of the weak solution. This extends the previous results for compressible NS [27] to NSP.

Nonlinear second-order multivalued boundary value problems

Indian Academy of Sciences (India)

R. Narasimhan (Krishtel eMaging) 1461 1996 Oct 15 13:05:22

Department of Mathematics, National Technical University, Zografou Campus,. Athens 15780 ... incorporates gradient systems, evolutionary variational inequalities and the classical boundary value ... We are led to an eventual application.

Integral Method of Boundary Characteristics: Neumann Condition

Science.gov (United States)

Kot, V. A.

2018-05-01

A new algorithm, based on systems of identical equalities with integral and differential boundary characteristics, is proposed for solving boundary-value problems on the heat conduction in bodies canonical in shape at a Neumann boundary condition. Results of a numerical analysis of the accuracy of solving heat-conduction problems with variable boundary conditions with the use of this algorithm are presented. The solutions obtained with it can be considered as exact because their errors comprise hundredths and ten-thousandths of a persent for a wide range of change in the parameters of a problem.

Performance improvement of extended boundary node method for solving elliptic boundary-value problems

International Nuclear Information System (INIS)

Saitoh, Ayumu; Kamitani, Atsushi; Takayama, Teruou; Nakamura, Hiroaki

2016-01-01

The extended boundary-node method (X-BNM) with the hierarchical-matrix (H-matrix) method has been developed and its performance has been investigated numerically. The results of computations show that the solver speed of the X-BNM with the H-matrix method is much faster than that of the standard X-BNM for the case where the number of boundary nodes exceeds a certain limit. Furthermore, the accuracy of the X-BNM with the H-matrix method is almost equal to that of the standard X-BNM. From these results, it is found that the H-matrix method is useful as the acceleration technique of the X-BNM. (author)

Indirect boundary element method for three dimensional problems. Analytical solution for contribution to wave field by triangular element; Sanjigen kansetsu kyokai yosoho. Sankakukei yoso no kiyo no kaisekikai

Energy Technology Data Exchange (ETDEWEB)

Yokoi, T [Building Research Institute, Tokyo (Japan); Sanchez-Sesma, F [Universidad National Autonoma de Mexico, (Mexico). Institute de Ingenieria

1997-05-27

Formulation is introduced for discretizing a boundary integral equation into an indirect boundary element method for the solution of 3-dimensional topographic problems. Yokoi and Takenaka propose an analytical solution-capable reference solution (solution for the half space elastic body with flat free surface) to problems of topographic response to seismic motion in a 2-dimensional in-plane field. That is to say, they propose a boundary integral equation capable of effectively suppressing the non-physical waves that emerge in the result of computation in the wake of the truncation of the discretized ground surface making use of the wave field in a semi-infinite elastic body with flat free surface. They apply the proposed boundary integral equation discretized into the indirect boundary element method to solve some examples, and succeed in proving its validity. In this report, the equation is expanded to deal with 3-dimensional topographic problems. A problem of a P-wave vertically landing on a flat and free surface is solved by the conventional boundary integral equation and the proposed boundary integral equation, and the solutions are compared with each other. It is found that the new method, different from the conventional one, can delete non-physical waves from the analytical result. 4 figs.

Two-point boundary value and Cauchy formulations in an axisymmetrical MHD equilibrium problem

International Nuclear Information System (INIS)

Atanasiu, C.V.; Subbotin, A.A.

1999-01-01

In this paper we present two equilibrium solvers for axisymmetrical toroidal configurations, both based on the expansion in poloidal angle method. The first one has been conceived as a two-point boundary value solver in a system of coordinates with straight field lines, while the second one uses a well-conditioned Cauchy formulation of the problem in a general curvilinear coordinate system. In order to check the capability of our moment methods to describe equilibrium accurately, a comparison of the moment solutions with analytical solutions obtained for a Solov'ev equilibrium has been performed. (author)

Existence and nonexistence results for a singular boundary value problem arising in the theory of epitaxial growth

Czech Academy of Sciences Publication Activity Database

Escudero, C.; Hakl, Robert; Peral, I.; Torres, P.J.

2014-01-01

Roč. 37, č. 6 (2014), s. 793-807 ISSN 0170-4214 Institutional support: RVO:67985840 Keywords : singular boundary value problem * epitaxial growth * radial solution Subject RIV: BA - General Mathematics Impact factor: 0.918, year: 2014 http://onlinelibrary.wiley.com/doi/10.1002/mma.2836/full

On Riemann boundary value problems for null solutions of the two dimensional Helmholtz equation

Science.gov (United States)

Bory Reyes, Juan; Abreu Blaya, Ricardo; Rodríguez Dagnino, Ramón Martin; Kats, Boris Aleksandrovich

2018-01-01

The Riemann boundary value problem (RBVP to shorten notation) in the complex plane, for different classes of functions and curves, is still widely used in mathematical physics and engineering. For instance, in elasticity theory, hydro and aerodynamics, shell theory, quantum mechanics, theory of orthogonal polynomials, and so on. In this paper, we present an appropriate hyperholomorphic approach to the RBVP associated to the two dimensional Helmholtz equation in R^2 . Our analysis is based on a suitable operator calculus.

Solution of problems with material nonlinearities with a coupled finite element/boundary element scheme using an iterative solver. Yucca Mountain Site Characterization Project

International Nuclear Information System (INIS)

Koteras, J.R.

1996-01-01

The prediction of stresses and displacements around tunnels buried deep within the earth is an important class of geomechanics problems. The material behavior immediately surrounding the tunnel is typically nonlinear. The surrounding mass, even if it is nonlinear, can usually be characterized by a simple linear elastic model. The finite element method is best suited for modeling nonlinear materials of limited volume, while the boundary element method is well suited for modeling large volumes of linear elastic material. A computational scheme that couples the finite element and boundary element methods would seem particularly useful for geomechanics problems. A variety of coupling schemes have been proposed, but they rely on direct solution methods. Direct solution techniques have large storage requirements that become cumbersome for large-scale three-dimensional problems. An alternative to direct solution methods is iterative solution techniques. A scheme has been developed for coupling the finite element and boundary element methods that uses an iterative solution method. This report shows that this coupling scheme is valid for problems where nonlinear material behavior occurs in the finite element region

Boundary conditions for the gravitational field

International Nuclear Information System (INIS)

Winicour, Jeffrey

2012-01-01

A review of the treatment of boundaries in general relativity is presented with the emphasis on application to the formulations of Einstein's equations used in numerical relativity. At present, it is known how to treat boundaries in the harmonic formulation of Einstein's equations and a tetrad formulation of the Einstein-Bianchi system. However, a universal approach valid for other formulations is not in hand. In particular, there is no satisfactory boundary theory for the 3+1 formulations which have been highly successful in binary black hole simulation. I discuss the underlying problems that make the initial-boundary-value problem much more complicated than the Cauchy problem. I review the progress that has been made and the important open questions that remain. Science is a differential equation. Religion is a boundary condition. (Alan Turing, quoted in J D Barrow, 'Theories of Everything') (topical review)

A phase change processor method for solving a one-dimensional phase change problem with convection boundary

Energy Technology Data Exchange (ETDEWEB)

Halawa, E.; Saman, W.; Bruno, F. [Institute for Sustainable Systems and Technologies, School of Advanced Manufacturing and Mechanical Engineering, University of South Australia, Mawson Lakes SA 5095 (Australia)

2010-08-15

A simple yet accurate iterative method for solving a one-dimensional phase change problem with convection boundary is described. The one-dimensional model takes into account the variation in the wall temperature along the direction of the flow as well as the sensible heat during preheating/pre-cooling of the phase change material (PCM). The mathematical derivation of convective boundary conditions has been integrated into a phase change processor (PCP) algorithm that solves the liquid fraction and temperature of the nodes. The algorithm is based on the heat balance at each node as it undergoes heating or cooling which inevitably involves phase change. The paper presents the model and its experimental validation. (author)

Existence of Positive Solutions to a Singular Semipositone Boundary Value Problem of Nonlinear Fractional Differential Systems

Directory of Open Access Journals (Sweden)

Xiaofeng Zhang

2017-12-01

Full Text Available In this paper, we consider the existence of positive solutions to a singular semipositone boundary value problem of nonlinear fractional differential equations. By applying the fixed point index theorem, some new results for the existence of positive solutions are obtained. In addition, an example is presented to demonstrate the application of our main results.

Accuracy Solution of Boundary Problems in Electrostatics for Systems "Conductors-Dielectrics" by Means of Auxiliary Charges

CERN Document Server

Topuriya, T P

2004-01-01

The analysis has been carried out on checking the influence of auxiliary charges on solution accuracy of boundary problems of electrostatics for systems "conductors-dielectrics". This accuracy depends on the number of charges and configuration of their allocation. The extended round dielectric in the electric field of a parallel-plate capacitor was taken as a physical model.

Optimal Wentzell Boundary Control of Parabolic Equations

International Nuclear Information System (INIS)

Luo, Yousong

2017-01-01

This paper deals with a class of optimal control problems governed by an initial-boundary value problem of a parabolic equation. The case of semi-linear boundary control is studied where the control is applied to the system via the Wentzell boundary condition. The differentiability of the state variable with respect to the control is established and hence a necessary condition is derived for the optimal solution in the case of both unconstrained and constrained problems. The condition is also sufficient for the unconstrained convex problems. A second order condition is also derived.

Optimal Wentzell Boundary Control of Parabolic Equations

Energy Technology Data Exchange (ETDEWEB)

Luo, Yousong, E-mail: yousong.luo@rmit.edu.au [RMIT University, School of Mathematical and Geospatial Sciences (Australia)

2017-04-15

This paper deals with a class of optimal control problems governed by an initial-boundary value problem of a parabolic equation. The case of semi-linear boundary control is studied where the control is applied to the system via the Wentzell boundary condition. The differentiability of the state variable with respect to the control is established and hence a necessary condition is derived for the optimal solution in the case of both unconstrained and constrained problems. The condition is also sufficient for the unconstrained convex problems. A second order condition is also derived.

Solution of free-boundary problems using finite-element/Newton methods and locally refined grids - Application to analysis of solidification microstructure

Science.gov (United States)

Tsiveriotis, K.; Brown, R. A.

1993-01-01

A new method is presented for the solution of free-boundary problems using Lagrangian finite element approximations defined on locally refined grids. The formulation allows for direct transition from coarse to fine grids without introducing non-conforming basis functions. The calculation of elemental stiffness matrices and residual vectors are unaffected by changes in the refinement level, which are accounted for in the loading of elemental data to the global stiffness matrix and residual vector. This technique for local mesh refinement is combined with recently developed mapping methods and Newton's method to form an efficient algorithm for the solution of free-boundary problems, as demonstrated here by sample calculations of cellular interfacial microstructure during directional solidification of a binary alloy.

Difference method for solving a nonlocal boundary value problem for a degenerating third-order pseudo-parabolic equation with variable coefficients

Science.gov (United States)

Beshtokov, M. Kh.

2016-10-01

A nonlocal boundary value problem for a degenerating third-order pseudo-parabolic equation with variable coefficients is considered. For solving this problem, a priori estimates in the differential and difference forms are obtained. The a priori estimates imply the uniqueness and stability of the solution on a layer with respect to the initial data and the right-hand side and the convergence of the solution of the difference problem to the solution of the differential problem.

Two-dimensional boundary-value problem for ion-ion diffusion

International Nuclear Information System (INIS)

Tuszewski, M.; Lichtenberg, A.J.

1977-01-01

Like-particle diffusion is usually negligible compared with unlike-particle diffusion because it is two orders higher in spatial derivatives. When the ratio of the ion gyroradius to the plasma transverse dimension is of the order of the fourth root of the mass ratio, previous one-dimensional analysis indicated that like-particle diffusion is significant. A two-dimensional boundary-value problem for ion-ion diffusion is investigated. Numerical solutions are found with models for which the nonlinear partial differential equation reduces to an ordinary fourth-order differential equation. These solutions indicate that the ion-ion losses are higher by a factor of six for a slab geometry, and by a factor of four for circular geometry, than estimated from dimensional analysis. The solutions are applied to a multiple mirror experiment stabilized with a quadrupole magnetic field which generates highly elliptical flux surfaces. It is found that the ion-ion losses dominate the electron-ion losses and that these classical radial losses contribute to a significant decrease of plasma lifetime, in qualitiative agreement with the experimental results

Boundary fluxes for non-local diffusion

OpenAIRE

Cortazar, C.; Elgueta, M.; Rossi, J. D.; Wolanski, N.

2006-01-01

We study a nonlocal diffusion operator in a bounded smooth domain prescribing the flux through the boundary. This problem may be seen as a generalization of the usual Neumann problem for the heat equation. First, we prove existence, uniqueness and a comparison principle. Next, we study the behavior of solutions for some prescribed boundary data including blowing up ones. Finally, we look at a nonlinear flux boundary condition.

Optimization of boundary controls of string vibrations

Energy Technology Data Exchange (ETDEWEB)

Il' in, V A; Moiseev, E I [Department of Computing Mathematics and Cybernetics, M.V. Lomonosov Moscow State University, Moscow (Russian Federation)

2005-12-31

For a large time interval T boundary controls of string vibrations are optimized in the following seven boundary-control problems: displacement control at one end (with the other end fixed or free); displacement control at both ends; elastic force control at one end (with the other end fixed or free); elastic force control at both ends; combined control (displacement control at one end and elastic force control at the other). Optimal boundary controls in each of these seven problems are sought as functions minimizing the corresponding boundary-energy integral under the constraints following from the initial and terminal conditions for the string at t=0 and t=T, respectively. For all seven problems, the optimal boundary controls are written out in closed analytic form.

Three-dimensional heat transfer analysis of the Doublet III beamline calorimeter

International Nuclear Information System (INIS)

Kamperschroer, J.H.; Pipkins, J.F.

1979-10-01

A general three-dimensional analysis has been formulated to study the flow of heat in a neutral beam calorimeter. The boundary value problem with an arbitrary incident heat flux has been solved using Fourier analysis and Laplace transform techniques. A general solution has been obtained and subsequently studied using numerical techniques as applied to the particular geometry and incident heat flux conditions of the Doublet III injection system. Negligible errors result in unfolding the incident heat flux through the use of thermocouples located near the rear surface, if data taking is initiated at the proper time and proceeds at a sufficiently rapid rate

A Boundary Property for Upper Domination

KAUST Repository

AbouEisha, Hassan M.

2016-08-08

An upper dominating set in a graph is a minimal (with respect to set inclusion) dominating set of maximum cardinality.The problem of finding an upper dominating set is generally NP-hard, but can be solved in polynomial time in some restricted graph classes, such as P4-free graphs or 2K2-free graphs.For classes defined by finitely many forbidden induced subgraphs, the boundary separating difficult instances of the problem from polynomially solvable ones consists of the so called boundary classes.However, none of such classes has been identified so far for the upper dominating set problem.In the present paper, we discover the first boundary class for this problem.

Extinction and focusing behaviour of spherical and annular flames described by a free boundary problem

Energy Technology Data Exchange (ETDEWEB)

Galaktionov, V A; Hulshof, J; Vazquez, J L

1997-09-01

We consider a free-boundary problem for the heat equation which arises in the description of premixed equi-diffusional flames in the limit of high activation energy. It consists of the heat equation u{sub t} = {Delta}u, u > 0, posed in an a priori unknown set {Omega} included in Q{sub T} = R{sup N} x (0,T) for some T >0 with boundary conditions on the free lateral boundary {Gamma}intersection between {partial_derivative}{Omega} et Q{sub T} (the flame front): u = 0 and {delta}u/{delta}{nu} = - 1. We impose initial condition u{sub 0}(x) {>=} 0 on the know initial domain {Omega}{sub 0} = interaction between {Omega}-bar et {l_brace} t = 0 {r_brace}. The paper establishes a theory of existence, uniqueness and regularity for radial symmetric solutions having bounded support. We remark that such solutions vanish in finite time (extinction phenomenon). in the paper we analyze the different types of possible extinction behaviour. We also investigate the focusing behaviour for solutions whose support expands in finite time to fill a hole. In all the cases the asymptotic behaviour is shown to be self-similar. (authors) 38 refs.

Recognition of boundary feedback systems

DEFF Research Database (Denmark)

Pedersen, Michael

1989-01-01

A system that has been the object of intense research is outlined. In view of that and recent progress of the theory of pseudodifferential boundary operator calculus, the author describes some features that could prove to be interesting in connection with the problems of boundary feedback stabili...... stabilizability. It is shown that it is possible to use the calculus to consider more general feedback systems in a variational setup.......A system that has been the object of intense research is outlined. In view of that and recent progress of the theory of pseudodifferential boundary operator calculus, the author describes some features that could prove to be interesting in connection with the problems of boundary feedback...

Polarity of translation boundaries in antiferroelectric PbZrO{sub 3}

Energy Technology Data Exchange (ETDEWEB)

Wei, Xian-Kui, E-mail: xiankui.wei@epfl.ch [Ceramics Laboratory, EPFL–Swiss Federal Institute of Technology, Lausanne 1015 (Switzerland); Peter Grünberg Institute and Ernst Ruska Center for Microscopy and Spectroscopy with Electrons, Research Center Jülich, 52425 Jülich (Germany); Jia, Chun-Lin [Peter Grünberg Institute and Ernst Ruska Center for Microscopy and Spectroscopy with Electrons, Research Center Jülich, 52425 Jülich (Germany); International Centre of Dielectric Research, The School of Electronic and Information Engineering, Xi’an Jiaotong University, Xi’an 710049 (China); Roleder, Krystian [Institute of Physics, University of Silesia, Katowice 40007 (Poland); Setter, Nava [Ceramics Laboratory, EPFL–Swiss Federal Institute of Technology, Lausanne 1015 (Switzerland)

2015-02-15

Graphical abstract: Strain-free rigid model and aberration-corrected transmission electron microscopes are used to investigate the polarity of translation boundaries in antiferroelectric PbZrO{sub 3}. - Highlights: • Domain boundaries in antiferroelectric PbZrO{sub 3} show polar and antipolar property. • The antiphase boundary can split into “sub-domains”. • Polarization reversal possibly exists inside the translation boundaries. • Thermal treatment can alter morphology and density of the translation boundaries. - Abstract: The polarity of translation boundaries (TBs) in antiferroelectric PbZrO{sub 3} is investigated. We show that previous experimentally reported polar property of R{sub III-1} type TB can be well approximated by a strain-free rigid model. Based on this, the modeling investigation suggests that there are two additional polar TBs, three antipolar-like TBs and one antipolar antiphase boundary. High-resolution scanning-transmission-electron-microscopy study reveals that the straight R{sub III-1} type TB can split into “sub-domains” with possible polarization reversal, suggesting the occurrence of ferroic orders at the TBs. In addition, dependence of morphology and density of the TBs on thermal treatments is discussed according to our results.

Existence and smoothness of solutions to second initial boundary value problems for Schrodinger systems in cylinders with non-smooth bases

Directory of Open Access Journals (Sweden)

Nguyen Manh Hung

2008-03-01

Full Text Available In this paper, we consider the second initial boundary value problem for strongly general Schrodinger systems in both the finite and the infinite cylinders $Q_T, 0problem are given.

Singular boundary value problem for the integrodifferential equation in an insurance model with stochastic premiums: Analysis and numerical solution

Science.gov (United States)

Belkina, T. A.; Konyukhova, N. B.; Kurochkin, S. V.

2012-10-01

A singular boundary value problem for a second-order linear integrodifferential equation with Volterra and non-Volterra integral operators is formulated and analyzed. The equation is defined on ℝ+, has a weak singularity at zero and a strong singularity at infinity, and depends on several positive parameters. Under natural constraints on the coefficients of the equation, existence and uniqueness theorems for this problem with given limit boundary conditions at singular points are proved, asymptotic representations of the solution are given, and an algorithm for its numerical determination is described. Numerical computations are performed and their interpretation is given. The problem arises in the study of the survival probability of an insurance company over infinite time (as a function of its initial surplus) in a dynamic insurance model that is a modification of the classical Cramer-Lundberg model with a stochastic process rate of premium under a certain investment strategy in the financial market. A comparative analysis of the results with those produced by the model with deterministic premiums is given.

Global solutions to the initial-boundary value problem for the quasilinear viscoelastic equation with a derivative nonlinearity

Directory of Open Access Journals (Sweden)

Mitsuhiro Nakao

2014-01-01

Full Text Available We prove the existence and uniqueness of a global decaying solution to the initial boundary value problem for the quasilinear wave equation with Kelvin-Voigt dissipation and a derivative nonlinearity. To derive the required estimates of the solutions we employ a 'loan' method and use a difference inequality on the energy.

Existence of Positive Solutions for a Coupled System of (p, q-Laplacian Fractional Higher Order Boundary Value Problems

Directory of Open Access Journals (Sweden)

K.R. Prasad

2015-11-01

Full Text Available In this paper, we establish the existence of at least three positive solutions for a system of (p,q-Laplacian fractional order two-point boundary value problems by applying five functionals fixed point theorem under suitable conditions on a cone in a Banach space.

Random walks in the quarter plane algebraic methods, boundary value problems, applications to queueing systems and analytic combinatorics

CERN Document Server

Fayolle, Guy; Malyshev, Vadim

2017-01-01

This monograph aims to promote original mathematical methods to determine the invariant measure of two-dimensional random walks in domains with boundaries. Such processes arise in numerous applications and are of interest in several areas of mathematical research, such as Stochastic Networks, Analytic Combinatorics, and Quantum Physics. This second edition consists of two parts. Part I is a revised upgrade of the first edition (1999), with additional recent results on the group of a random walk. The theoretical approach given therein has been developed by the authors since the early 1970s. By using Complex Function Theory, Boundary Value Problems, Riemann Surfaces, and Galois Theory, completely new methods are proposed for solving functional equations of two complex variables, which can also be applied to characterize the Transient Behavior of the walks, as well as to find explicit solutions to the one-dimensional Quantum Three-Body Problem, or to tackle a new class of Integrable Systems. Part II borrows spec...

Boundary-value problems with integral conditions for a system of Lame equations in the space of almost periodic functions

Directory of Open Access Journals (Sweden)

Volodymyr S. Il'kiv

2016-11-01

Full Text Available We study a problem with integral boundary conditions in the time coordinate for a system of Lame equations of dynamic elasticity theory of an arbitrary dimension. We find necessary and sufficient conditions for the existence and uniqueness of solution in the class of almost periodic functions in the spatial variables. To solve the problem of small denominators arising while constructing solutions, we use the metric approach.

Memory allocation and computations for Laplace’s equation of 3-D arbitrary boundary problems

Directory of Open Access Journals (Sweden)

Tsay Tswn-Syau

2017-01-01

Full Text Available Computation iteration schemes and memory allocation technique for finite difference method were presented in this paper. The transformed form of a groundwater flow problem in the generalized curvilinear coordinates was taken to be the illustrating example and a 3-dimensional second order accurate 19-point scheme was presented. Traditional element-by-element methods (e.g. SOR are preferred since it is simple and memory efficient but time consuming in computation. For efficient memory allocation, an index method was presented to store the sparse non-symmetric matrix of the problem. For computations, conjugate-gradient-like methods were reported to be computationally efficient. Among them, using incomplete Choleski decomposition as preconditioner was reported to be good method for iteration convergence. In general, the developed index method in this paper has the following advantages: (1 adaptable to various governing and boundary conditions, (2 flexible for higher order approximation, (3 independence of problem dimension, (4 efficient for complex problems when global matrix is not symmetric, (5 convenience for general sparse matrices, (6 computationally efficient in the most time consuming procedure of matrix multiplication, and (7 applicable to any developed matrix solver.

Differential and Difference Boundary Value Problem for Loaded Third-Order Pseudo-Parabolic Differential Equations and Difference Methods for Their Numerical Solution

Science.gov (United States)

Beshtokov, M. Kh.

2017-12-01

Boundary value problems for loaded third-order pseudo-parabolic equations with variable coefficients are considered. A priori estimates for the solutions of the problems in the differential and difference formulations are obtained. These a priori estimates imply the uniqueness and stability of the solution with respect to the initial data and the right-hand side on a layer, as well as the convergence of the solution of each difference problem to the solution of the corresponding differential problem.

On Impulsive Boundary Value Problems of Fractional Differential Equations with Irregular Boundary Conditions

Directory of Open Access Journals (Sweden)

Guotao Wang

2012-01-01

Full Text Available We study nonlinear impulsive differential equations of fractional order with irregular boundary conditions. Some existence and uniqueness results are obtained by applying standard fixed-point theorems. For illustration of the results, some examples are discussed.

EXISTENCE OF POSITIVE SOLUTION TO TWO-POINT BOUNDARY VALUE PROBLEM FOR A SYSTEM OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATIONS

Institute of Scientific and Technical Information of China (English)

无

2009-01-01

In this paper, we consider a two-point boundary value problem for a system of second order ordinary differential equations. Under some conditions, we show the existence of positive solution to the system of second order ordinary differential equa-tions.

Initial-boundary value problems for multi-term time-fractional diffusion equations with x-dependent coefficients

OpenAIRE

Li, Zhiyuan; Huang, Xinchi; Yamamoto, Masahiro

2018-01-01

In this paper, we discuss an initial-boundary value problem (IBVP) for the multi-term time-fractional diffusion equation with x-dependent coefficients. By means of the Mittag-Leffler functions and the eigenfunction expansion, we reduce the IBVP to an equivalent integral equation to show the unique existence and the analyticity of the solution for the equation. Especially, in the case where all the coefficients of the time-fractional derivatives are non-negative, by the Laplace and inversion L...

New formulations on the finite element method for boundary value problems with internal/external boundary layers; Novas formulacoes de elementos finitos para problemas de valor de contorno com camadas limite interna/externa

Energy Technology Data Exchange (ETDEWEB)

Pereira, Luis Carlos Martins

1998-06-15

New Petrov-Galerkin formulations on the finite element methods for convection-diffusion problems with boundary layers are presented. Such formulations are based on a consistent new theory on discontinuous finite element methods. Existence and uniqueness of solutions for these problems in the new finite element spaces are demonstrated. Some numerical experiments shows how the new formulation operate and also their efficacy. (author)

Lower and Upper Solutions Method for Positive Solutions of Fractional Boundary Value Problems

Directory of Open Access Journals (Sweden)

R. Darzi

2013-01-01

Full Text Available We apply the lower and upper solutions method and fixed-point theorems to prove the existence of positive solution to fractional boundary value problem D0+αut+ft,ut=0, 0

Einstein boundary conditions in relation to constraint propagation for the initial-boundary value problem of the Einstein equations

International Nuclear Information System (INIS)

Frittelli, Simonetta; Gomez, Roberto

2004-01-01

We show how the use of the normal projection of the Einstein tensor as a set of boundary conditions relates to the propagation of the constraints, for two representations of the Einstein equations with vanishing shift vector: the Arnowitt-Deser-Misner formulation, which is ill posed, and the Einstein-Christoffel formulation, which is symmetric hyperbolic. Essentially, the components of the normal projection of the Einstein tensor that act as nontrivial boundary conditions are linear combinations of the evolution equations with the constraints that are not preserved at the boundary, in both cases. In the process, the relationship of the normal projection of the Einstein tensor to the recently introduced 'constraint-preserving' boundary conditions becomes apparent

Triple Positive Solutions of a Nonlocal Boundary Value Problem for Singular Differential Equations with p-Laplacian

Directory of Open Access Journals (Sweden)

Jufang Wang

2013-01-01

Full Text Available We establish the existence of triple positive solutions of an m-point boundary value problem for the nonlinear singular second-order differential equations of mixed type with a p-Laplacian operator by Leggett-William fixed point theorem. At last, we give an example to demonstrate the use of the main result of this paper. The conclusions in this paper essentially extend and improve the known results.

Numeric treatment of nonlinear second order multi-point boundary value problems using ANN, GAs and sequential quadratic programming technique

Directory of Open Access Journals (Sweden)

Zulqurnain Sabir

2014-06-01

Full Text Available In this paper, computational intelligence technique are presented for solving multi-point nonlinear boundary value problems based on artificial neural networks, evolutionary computing approach, and active-set technique. The neural network is to provide convenient methods for obtaining useful model based on unsupervised error for the differential equations. The motivation for presenting this work comes actually from the aim of introducing a reliable framework that combines the powerful features of ANN optimized with soft computing frameworks to cope with such challenging system. The applicability and reliability of such methods have been monitored thoroughly for various boundary value problems arises in science, engineering and biotechnology as well. Comprehensive numerical experimentations have been performed to validate the accuracy, convergence, and robustness of the designed scheme. Comparative studies have also been made with available standard solution to analyze the correctness of the proposed scheme.

Moving-Boundary Problems Associated with Lyopreservation

Science.gov (United States)

Gruber, Christopher Andrew

The work presented in this Dissertation is motivated by research into the preservation of biological specimens by way of vitrification, a technique known as lyopreservation. The operative principle behind lyopreservation is that a glassy material forms as a solution of sugar and water is desiccated. The microstructure of this glass impedes transport within the material, thereby slowing metabolism and effectively halting the aging processes in a biospecimen. This Dissertation is divided into two segments. The first concerns the nature of diffusive transport within a glassy state. Experimental studies suggest that diffusion within a glass is anomalously slow. Scaled Brownian motion (SBM) is proposed as a mathematical model which captures the qualitative features of anomalously slow diffusion while minimizing computational expense. This model is applied to several moving-boundary problems and the results are compared to a more well-established model, fractional anomalous diffusion (FAD). The virtues of SBM are based on the model's relative mathematical simplicity: the governing equation under FAD dynamics involves a fractional derivative operator, which precludes the use of analytical methods in almost all circumstances and also entails great computational expense. In some geometries, SBM allows similarity solutions, though computational methods are generally required. The use of SBM as an approximation to FAD when a system is "nearly classical'' is also explored. The second portion of this Dissertation concerns spin-drying, which is an experimental approach to biopreservation in a laboratory setting. A biospecimen is adhered to a glass wafer and this substrate is covered with sugar solution and rapidly spun on a turntable while water is evaporated from the film surface. The mathematical model for the spin-drying process includes diffusion, viscous fluid flow, and evaporation, among other contributions to the dynamics. Lubrication theory is applied to the model and an

Grain boundary structure and properties

International Nuclear Information System (INIS)

Balluffi, R.W.

1979-01-01

An attempt is made to distinguish those fundamental aspects of grain boundaries which should be relevant to the problem of the time dependent fracture of high temperature structural materials. These include the basic phenomena which are thought to be associated with cavitation and cracking at grain boundaries during service and with the more general microstructural changes which occur during both processing and service. A very brief discussion of the current state of our knowledge of these fundamentals is given. Included are the following: (1) structure of ideal perfect boundaries; (2) defect structure of grain boundaries; (3) diffusion at grain boundaries; (4) grain boundaries as sources/sinks for point defects; (5) grain boundary migration; (6) dislocation phenomena at grain boundaries; (7) atomic bonding and cohesion at grain boundaries; (8) non-equilibrium properties of grain boundaries; and (9) techniques for studying grain boundaries

An analytical boundary element integral approach to track the boundary of a moving cavity using electrical impedance tomography

International Nuclear Information System (INIS)

Khambampati, Anil Kumar; Kim, Sin; Lee, Bo An; Kim, Kyung Youn

2012-01-01

This paper is about locating the boundary of a moving cavity within a homogeneous background from the voltage measurements recorded on the outer boundary. An inverse boundary problem of a moving cavity is formulated by considering a two-phase vapor–liquid flow in a pipe. The conductivity of the flow components (vapor and liquid) is assumed to be constant and known a priori while the location and shape of the inclusion (vapor) are the unknowns to be estimated. The forward problem is solved using the boundary element method (BEM) with the integral equations solved analytically. A special situation is considered such that the cavity changes its location and shape during the time taken to acquire a full set of independent measurement data. The boundary of a cavity is assumed to be elliptic and is parameterized with Fourier series. The inverse problem is treated as a state estimation problem with the Fourier coefficients that represent the center and radii of the cavity as the unknowns to be estimated. An extended Kalman filter (EKF) is used as an inverse algorithm to estimate the time varying Fourier coefficients. Numerical experiments are shown to evaluate the performance of the proposed method. Through the results, it can be noticed that the proposed BEM with EKF method is successful in estimating the boundary of a moving cavity. (paper)

An approximate moving boundary method for American option pricing

NARCIS (Netherlands)

Chockalingam, A.; Muthuraman, K.

2015-01-01

We present a method to solve the free-boundary problem that arises in the pricing of classical American options. Such free-boundary problems arise when one attempts to solve optimal-stopping problems set in continuous time. American option pricing is one of the most popular optimal-stopping problems

A problem of atomic diffusion in a moving boundary

International Nuclear Information System (INIS)

Bezerra, M.C.C.

1985-01-01

It is analysed the convergence of a numerical scheme for calculating approximate solutions of a model used for evaluating concentration of atoms in a diffusion process in the walls of nuclear reactors. The ion trapping process is admitted to be reversible and the wall corrosion process is also considered in the model so that must deal with a moving boundary. Some conditions for the motion of the boundary are established in such a way that convergence can be assured in more general settings than those of previous papers. (Author) [pt

Three Dimensional Energy Transmitting Boundary in the Time Domain

Directory of Open Access Journals (Sweden)

Naohiro eNakamura

2015-11-01

Full Text Available Although the energy transmitting boundary is accurate and efficient for the FEM earthquake response analysis, it could be applied in the frequency domain only. In the previous papers, the author proposed an earthquake response analysis method using the time domain energy transmitting boundary for two dimensional problems. In this paper, this technique is expanded for three dimensional problems. The inner field is supposed to be a hexahedron shape and the approximate time domain boundary is explained, first. Next, two dimensional anti-plane time domain boundary is studied for a part of the approximate three dimensional boundary method. Then, accuracy and efficiency of the proposed method are confirmed by example problems.

Non-surgical alternative in the treatment of skeletal Class III problems.

Science.gov (United States)

Jefferson, Y

1995-01-01

The dental profession is not static, but dynamic. New research findings, along with medical and technological advances, necessitate constant re-examination of treatment philosophies and techniques. What were acceptable treatment techniques in the past may not necessarily be the most effective and best techniques for our patients today. Currently, many practitioners feel that the only treatment for the correction of a skeletal Class III abnormality is via orthognathic surgery in older patients. In some cases it may be the only treatment option. But in most cases today, there are more conservative, non-surgical treatment alternatives in correcting Class III problems in younger aged children. In treating facial-skeletal problems, it must be emphasized that the human face is a biological masterpiece of form and function. Its importance has been documented in arts and sciences since the beginning of modern civilization. It is important enough so that individuals who are blessed with attractive features are afforded greater opportunities in our society. Attractive faces are associated with intelligence, honesty and good work ethics. With the advent of orthognathic surgery, functional appliance, functional regulator, and myofunctional therapy, the dental profession has the capability of leveling out the playing field for many individuals in our society. It does so by being able to correct problems closely associated with the human psyche--the human face. The ability to change facial features brings tremendous prestige to our profession. Along with this prestige comes greater responsibility. Our ability to change facial features entails greater understanding of facial balance and harmony. Ricketts states that the face must conform to stringent proportions known as the "divine proportion" in order for it to be esthetically pleasing. Also, our ability to move facial-skeletal structures entails greater understanding of the biomechanics of the human face. Without this

Regularity of the solutions to a nonlinear boundary problem with indefinite weight

Directory of Open Access Journals (Sweden)

Aomar Anane

2011-01-01

Full Text Available In this paper we study the regularity of the solutions to the problemDelta_p u = |u|^{p−2}u in the bounded smooth domainOmega ⊂ R^N,with|∇u|^{p−2} partial_{nu} u = lambda V (x|u|^{p−2}u + h as a nonlinear boundary condition, where partial Omega is C^{2,beta}, with beta ∈]0, 1[, and V is a weight in L^s(partial Omega and h ∈ L^s(partial Omega for some s ≥ 1. We prove that all solutions are in L^{infty}(Omega cap L^{infty}(Omega, and using the D.Debenedetto’s theorem of regularity in [1] we conclude that those solutions are in C^{1,alpha} overline{Omega} for some alpha ∈ ]0, 1[.

An asymptotic analytical solution to the problem of two moving boundaries with fractional diffusion in one-dimensional drug release devices

International Nuclear Information System (INIS)

Yin Chen; Xu Mingyu

2009-01-01

We set up a one-dimensional mathematical model with a Caputo fractional operator of a drug released from a polymeric matrix that can be dissolved into a solvent. A two moving boundaries problem in fractional anomalous diffusion (in time) with order α element of (0, 1] under the assumption that the dissolving boundary can be dissolved slowly is presented in this paper. The two-parameter regular perturbation technique and Fourier and Laplace transform methods are used. A dimensionless asymptotic analytical solution is given in terms of the Wright function

The theory of discrete barriers and its applications to linear boundary-value problems of the 'Dirichlet type'; Theorie des barrieres discretes et applications a des problemes lineaires elliptiques du ''type de dirichlet''

Energy Technology Data Exchange (ETDEWEB)

Jamet, P [Commissariat a l' Energie Atomique, Saclay (France). Centre d' Etudes Nucleaires

1967-07-01

This report gives a general presentation of barrier theory for finite difference operators, with its applications to some boundary value problems. (author) [French] Ce rapport est un expose synthetique de la theorie des barrieres pour les operateurs aux differences finies et ses applications a certaines classes de problemes lineaires elliptiques du 'type de Dirichlet'. (auteur)

A Matlab-based finite-difference solver for the Poisson problem with mixed Dirichlet-Neumann boundary conditions

Science.gov (United States)

Reimer, Ashton S.; Cheviakov, Alexei F.

2013-03-01

A Matlab-based finite-difference numerical solver for the Poisson equation for a rectangle and a disk in two dimensions, and a spherical domain in three dimensions, is presented. The solver is optimized for handling an arbitrary combination of Dirichlet and Neumann boundary conditions, and allows for full user control of mesh refinement. The solver routines utilize effective and parallelized sparse vector and matrix operations. Computations exhibit high speeds, numerical stability with respect to mesh size and mesh refinement, and acceptable error values even on desktop computers. Catalogue identifier: AENQ_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AENQ_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License v3.0 No. of lines in distributed program, including test data, etc.: 102793 No. of bytes in distributed program, including test data, etc.: 369378 Distribution format: tar.gz Programming language: Matlab 2010a. Computer: PC, Macintosh. Operating system: Windows, OSX, Linux. RAM: 8 GB (8, 589, 934, 592 bytes) Classification: 4.3. Nature of problem: To solve the Poisson problem in a standard domain with “patchy surface”-type (strongly heterogeneous) Neumann/Dirichlet boundary conditions. Solution method: Finite difference with mesh refinement. Restrictions: Spherical domain in 3D; rectangular domain or a disk in 2D. Unusual features: Choice between mldivide/iterative solver for the solution of large system of linear algebraic equations that arise. Full user control of Neumann/Dirichlet boundary conditions and mesh refinement. Running time: Depending on the number of points taken and the geometry of the domain, the routine may take from less than a second to several hours to execute.

Nature Inspired Computational Technique for the Numerical Solution of Nonlinear Singular Boundary Value Problems Arising in Physiology

Directory of Open Access Journals (Sweden)

Suheel Abdullah Malik

2014-01-01

Full Text Available We present a hybrid heuristic computing method for the numerical solution of nonlinear singular boundary value problems arising in physiology. The approximate solution is deduced as a linear combination of some log sigmoid basis functions. A fitness function representing the sum of the mean square error of the given nonlinear ordinary differential equation (ODE and its boundary conditions is formulated. The optimization of the unknown adjustable parameters contained in the fitness function is performed by the hybrid heuristic computation algorithm based on genetic algorithm (GA, interior point algorithm (IPA, and active set algorithm (ASA. The efficiency and the viability of the proposed method are confirmed by solving three examples from physiology. The obtained approximate solutions are found in excellent agreement with the exact solutions as well as some conventional numerical solutions.

A finite difference method for free boundary problems

KAUST Repository

Fornberg, Bengt

2010-01-01

Fornberg and Meyer-Spasche proposed some time ago a simple strategy to correct finite difference schemes in the presence of a free boundary that cuts across a Cartesian grid. We show here how this procedure can be combined with a minimax

Obstacle problems in mathematical physics

CERN Document Server

Rodrigues, J-F

1987-01-01

The aim of this research monograph is to present a general account of the applicability of elliptic variational inequalities to the important class of free boundary problems of obstacle type from a unifying point of view of classical Mathematical Physics.The first part of the volume introduces some obstacle type problems which can be reduced to variational inequalities. Part II presents some of the main aspects of the theory of elliptic variational inequalities, from the abstract hilbertian framework to the smoothness of the variational solution, discussing in general the properties of the free boundary and including some results on the obstacle Plateau problem. The last part examines the application to free boundary problems, namely the lubrication-cavitation problem, the elastoplastic problem, the Signorini (or the boundary obstacle) problem, the dam problem, the continuous casting problem, the electrochemical machining problem and the problem of the flow with wake in a channel past a profile.

Boundary value problems of the circular cylinders in the strain-gradient theory of linear elasticity

International Nuclear Information System (INIS)

Kao, B.G.

1979-11-01

Three boundary value problems in the strain-gradient theory of linear elasticity are solved for circular cylinders. They are the twisting of circular cylinder, uniformly pressuring of concentric circular cylinder, and pure-bending of simply connected cylinder. The comparisons of these solutions with the solutions in classical elasticity and in couple-stress theory reveal the differences in the stress fields as well as the apparent stress fields due to the influences of the strain-gradient. These aspects of the strain-gradient theory could be important in modeling the failure behavior of structural materials

About solution of multipoint boundary problem of static analysis of deep beam with the use of combined application of finite element method and discrete-continual finite element method. part 1: formulation of the problem and general principles of approximation

Directory of Open Access Journals (Sweden)

Lyakhovich Leonid

2017-01-01

Full Text Available This paper is devoted to formulation and general principles of approximation of multipoint boundary problem of static analysis of deep beam with the use of combined application of finite element method (FEM discrete-continual finite element method (DCFEM. The field of application of DCFEM comprises structures with regular physical and geometrical parameters in some dimension (“basic” dimension. DCFEM presupposes finite element approximation for non-basic dimension while in the basic dimension problem remains continual. DCFEM is based on analytical solutions of resulting multipoint boundary problems for systems of ordinary differential equations with piecewise-constant coefficients.

Asymptotic analysis and boundary layers

CERN Document Server

Cousteix, Jean

2007-01-01

This book presents a new method of asymptotic analysis of boundary-layer problems, the Successive Complementary Expansion Method (SCEM). The first part is devoted to a general comprehensive presentation of the tools of asymptotic analysis. It gives the keys to understand a boundary-layer problem and explains the methods to construct an approximation. The second part is devoted to SCEM and its applications in fluid mechanics, including external and internal flows. The advantages of SCEM are discussed in comparison with the standard Method of Matched Asymptotic Expansions. In particular, for the first time, the theory of Interactive Boundary Layer is fully justified. With its chapter summaries, detailed derivations of results, discussed examples and fully worked out problems and solutions, the book is self-contained. It is written on a mathematical level accessible to graduate and post-graduate students of engineering and physics with a good knowledge in fluid mechanics. Researchers and practitioners will estee...

Numerical analysis of melting/solidification phenomena using a moving boundary problem analysis method X-FEM

International Nuclear Information System (INIS)

Uchibori, Akihiro; Ohshima, Hiroyuki

2008-01-01

A numerical analysis method for melting/solidification phenomena has been developed to evaluate a feasibility of several candidate techniques in the nuclear fuel cycle. Our method is based on the eXtended Finite Element Method (X-FEM) which has been used for moving boundary problems. Key technique of the X-FEM is to incorporate signed distance function into finite element interpolation to represent a discontinuous gradient of the temperature at a moving solid-liquid interface. Construction of the finite element equation, the technique of quadrature and the method to solve the equation are reported here. The numerical solutions of the one-dimensional Stefan problem, solidification in a two-dimensional square corner and melting of pure gallium are compared to the exact solutions or to the experimental data. Through these analyses, validity of the newly developed numerical analysis method has been demonstrated. (author)

Using reciprocity in Boundary Element Calculations

DEFF Research Database (Denmark)

Juhl, Peter Møller; Cutanda Henriquez, Vicente

2010-01-01

The concept of reciprocity is widely used in both theoretical and experimental work. In Boundary Element calculations reciprocity is sometimes employed in the solution of computationally expensive scattering problems, which sometimes can be more efficiently dealt with when formulated...... as the reciprocal radiation problem. The present paper concerns the situation of having a point source (which is reciprocal to a point receiver) at or near a discretized boundary element surface. The accuracy of the original and the reciprocal problem is compared in a test case for which an analytical solution...

Exact multiplicity results for quasilinear boundary-value problems with cubic-like nonlinearities

Directory of Open Access Journals (Sweden)

Idris Addou

2000-01-01

Full Text Available We consider the boundary-value problem $$displaylines{ -(varphi_p (u'' =lambda f(u mbox{ in }(0,1 cr u(0 = u(1 =0,, }$$ where $p>1$, $lambda >0$ and $varphi_p (x =| x|^{p-2}x$. The nonlinearity $f$ is cubic-like with three distinct roots 0=a less than b less than c. By means of a quadrature method, we provide the exact number of solutions for all $lambda >0$. This way we extend a recent result, for $p=2$, by Korman et al. cite{KormanLiOuyang} to the general case $p>1$. We shall prove that when 1less than $pleq 2$ the structure of the solution set is exactly the same as that studied in the case $p=2$ by Korman et al. cite{KormanLiOuyang}, and strictly different in the case $p>2$.

A nonlinear boundary integral equations method for the solving of quasistatic elastic contact problem with Coulomb friction

Directory of Open Access Journals (Sweden)

Yurii M. Streliaiev

2016-06-01

Full Text Available Three-dimensional quasistatic contact problem of two linearly elastic bodies' interaction with Coulomb friction taken into account is considered. The boundary conditions of the problem have been simplified by the modification of the Coulomb's law of friction. This modification is based on the introducing of a delay in normal contact tractions that bound tangent contact tractions in the Coulomb's law of friction expressions. At this statement the problem is reduced to a sequence of similar systems of nonlinear integral equations describing bodies' interaction at each step of loading. A method for an approximate solution of the integral equations system corresponded to each step of loading is applied. This method consists of system regularization, discretization of regularized system and iterative process application for solving the discretized system. A numerical solution of a contact problem of an elastic sphere with an elastic half-space interaction under increasing and subsequently decreasing normal compressive force has been obtained.

Problems in mathematical analysis III integration

CERN Document Server

Kaczor, W J

2003-01-01

We learn by doing. We learn mathematics by doing problems. This is the third volume of Problems in Mathematical Analysis. The topic here is integration for real functions of one real variable. The first chapter is devoted to the Riemann and the Riemann-Stieltjes integrals. Chapter 2 deals with Lebesgue measure and integration. The authors include some famous, and some not so famous, integral inequalities related to Riemann integration. Many of the problems for Lebesgue integration concern convergence theorems and the interchange of limits and integrals. The book closes with a section on Fourier series, with a concentration on Fourier coefficients of functions from particular classes and on basic theorems for convergence of Fourier series. The book is primarily geared toward students in analysis, as a study aid, for problem-solving seminars, or for tutorials. It is also an excellent resource for instructors who wish to incorporate problems into their lectures. Solutions for the problems are provided in the boo...

Application of the Laplace transform method for the albedo boundary conditions in multigroup neutron diffusion eigenvalue problems in slab geometry

International Nuclear Information System (INIS)

Petersen, Claudio Zen; Vilhena, Marco T.; Barros, Ricardo C.

2009-01-01

In this paper the application of the Laplace transform method is described in order to determine the energy-dependent albedo matrix that is used in the boundary conditions multigroup neutron diffusion eigenvalue problems in slab geometry for nuclear reactor global calculations. In slab geometry, the diffusion albedo substitutes without approximation the baffle-reflector system around the active domain. Numerical results to typical test problems are shown to illustrate the accuracy and the efficiency of the Chebysheff acceleration scheme. (orig.)

Free-boundary models of a meltwater conduit

KAUST Repository

Dallaston, Michael C.

2014-08-01

© 2014 AIP Publishing LLC. We analyse the cross-sectional evolution of an englacial meltwater conduit that contracts due to inward creep of the surrounding ice and expands due to melting. Making use of theoretical methods from free-boundary problems in Stokes flow and Hele-Shaw squeeze flow we construct an exact solution to the coupled problem of external viscous creep and internal heating, in which we adopt a Newtonian approximation for ice flow and an idealized uniform heat source in the conduit. This problem provides an interesting variant on standard free-boundary problems, coupling different internal and external problems through the kinematic condition at the interface. The boundary in the exact solution takes the form of an ellipse that may contract or expand (depending on the magnitudes of effective pressure and heating rate) around fixed focal points. Linear stability analysis reveals that without the melting this solution is unstable to perturbations in the shape. Melting can stabilize the interface unless the aspect ratio is too small; in that case, instabilities grow largest at the thin ends of the ellipse. The predictions are corroborated with numerical solutions using boundary integral techniques. Finally, a number of extensions to the idealized model are considered, showing that a contracting circular conduit is unstable to all modes of perturbation if melting occurs at a uniform rate around the boundary, or if the ice is modelled as a shear-thinning fluid.

Axisymmetric MHD stability of sharp-boundary Tokamaks

International Nuclear Information System (INIS)

Rebhan, E.; Salat, A.

1976-09-01

For a sharp-boundary, constant pressure plasma model of axisymmetric equilibria the MHD stability problem of axisymmetric perturbations is solved by analytic reduction to a one-dimensional problem on the boundary and subsequent numerical treatment, using the energy principle. The stability boundaries are determined for arbitrary aspect ratio, arbitrary βsub(p) and elliptical, triangular and rectangular plasma cross-sections, wall stabilization not being taken into account. It is found that the axisymmetric stability strongly depends on the plasma shape and is almost independent of the safety factor q. (orig.) [de

The Use of Source-Sink and Doublet Distributions Extended to the Solution of Boundary-Value Problems in Supersonic Flow

Science.gov (United States)

Heaslet, Max A; Lomax, Harvard

1948-01-01

A direct analogy is established between the use of source-sink and doublet distributions in the solution of specific boundary-value problems in subsonic wing theory and the corresponding problems in supersonic theory. The correct concept of the "finite part" of an integral is introduced and used in the calculation of the improper integrals associated with supersonic doublet distributions. The general equations developed are shown to include several previously published results and particular examples are given for the loading on rolling and pitching triangular wings with supersonic leading edges.

Efficient Closed Form Cut-Off Planes and Propagation Planes Characteristics for Dielectric Slab Loaded Boundary Value Problems

OpenAIRE

Zafar, Junaid

2012-01-01

The geometrical relationship between the cut-off and propagating planes of any waveguide system is a prerequisite for any design process. The characterization of cut-off planes and optimisation are challenging for numerical methods, closed-form solutions are always preferred. In this paper Maxwells coupled field equations are used to characterise twin E-plane and H-plane slab loaded boundary value problems. The single mode bandwidths and dispersion characteristics of these structures are pres...

Fictitious domain methods for elliptic problems with general boundary conditions with an application to the numerical simulation of two phase flows; Methodes de domaine fictif pour des problemes elliptiques avec conditions aux limites generales en vue de la simulation numerique d'ecoulements diphasiques

Energy Technology Data Exchange (ETDEWEB)

Ramiere, I

2006-09-15

This work is dedicated to the introduction of two original fictitious domain methods for the resolution of elliptic problems (mainly convection-diffusion problems) with general and eventually mixed boundary conditions: Dirichlet, Robin or Neumann. The originality lies in the approximation of the immersed boundary by an approximate interface derived from the fictitious domain Cartesian mesh, which is generally not boundary-fitted to the physical domain. The same generic numerical scheme is used to impose the embedded boundary conditions. Hence, these methods require neither a surface mesh of the immersed boundary nor the local modification of the numerical scheme. We study two modelling of the immersed boundary. In the first one, called spread interface, the approximate immersed boundary is the union of the cells crossed by the physical immersed boundary. In the second one, called thin interface, the approximate immersed boundary lies on sides of mesh cells. Additional algebraic transmission conditions linking both flux and solution jumps through the thin approximate interface are introduced. The fictitious problem to solve as well as the treatment of the embedded boundary conditions are detailed for the two methods. A Q1 finite element scheme is implemented for the numerical validation of the spread interface approach while a new cell-centered finite volume scheme is derived for the thin interface approach with immersed jumps. Each method is then combined to multilevel local mesh refinement algorithms (with solution or flux residual) to increase the precision of the solution in the vicinity of the immersed interface. A convergence analysis of a Q1 finite element method with non-boundary fitted meshes is also presented. This study proves the convergence rates of the present methods. Among the various industrial applications, the simulation on a model of heat exchanger in french nuclear power plants enables us to appreciate the performances of the fictitious domain

The atmospheric boundary layer — advances in knowledge and application

Science.gov (United States)

Garratt, J. R.; Hess, G. D.; Physick, W. L.; Bougeault, P.

1996-02-01

We summarise major activities and advances in boundary-layer knowledge in the 25 years since 1970, with emphasis on the application of this knowledge to surface and boundary-layer parametrisation schemes in numerical models of the atmosphere. Progress in three areas is discussed: (i) the mesoscale modelling of selected phenomena; (ii) numerical weather prediction; and (iii) climate simulations. Future trends are identified, including the incorporation into models of advanced cloud schemes and interactive canopy schemes, and the nesting of high resolution boundary-layer schemes in global climate models.

Coherent structures in wave boundary layers. Part 2. Solitary motion

DEFF Research Database (Denmark)

Sumer, B. Mutlu; Jensen, Palle Martin; Sørensen, Lone B.

2010-01-01

This study continues the investigation of wave boundary layers reported by Carstensen, Sumer & Fredsøe (J. Fluid Mech., 2010, part 1 of this paper). The present paper summarizes the results of an experimental investigation of turbulent solitary wave boundary layers, simulated by solitary motion...... the boundary-layer flow experiences a regular array of vortex tubes near the bed over a short period of time during the deceleration stage; and (iii) transitional regime characterized with turbulent spots, revealed by single/multiple, or, sometimes, quite dense spikes in the bed shear stress traces...

Boundary as Bridge: An Analysis of the Educational Neuroscience Literature from a Boundary Perspective

Science.gov (United States)

Beauchamp, Catherine; Beauchamp, Miriam H.

2013-01-01

Within the emerging field of educational neuroscience, concerns exist that the impact of neuroscience research on education has been less effective than hoped. In seeking a way forward, it may be useful to consider the problems of integrating two complex fields in the context of disciplinary boundaries. Here, a boundary perspective is used as a…

Multiple Positive Solutions of a Nonlinear Four-Point Singular Boundary Value Problem with a p-Laplacian Operator on Time Scales

Directory of Open Access Journals (Sweden)

Shihuang Hong

2009-01-01

Full Text Available We present sufficient conditions for the existence of at least twin or triple positive solutions of a nonlinear four-point singular boundary value problem with a p-Laplacian dynamic equation on a time scale. Our results are obtained via some new multiple fixed point theorems.

Upper Domination: Towards a Dichotomy Through Boundary Properties

KAUST Repository

AbouEisha, Hassan M.

2017-07-14

An upper dominating set in a graph is a minimal dominating set of maximum cardinality. The problem of finding an upper dominating set is generally NP-hard. We study the complexity of this problem in finitely defined classes of graphs and conjecture that the problem admits a complexity dichotomy in this family. A helpful tool to study the complexity of an algorithmic problem is the notion of boundary classes. However, none of such classes has been identified so far for the upper dominating set problem. We discover the first boundary class for this problem and prove the dichotomy for monogenic classes.

Blurring Boundaries

DEFF Research Database (Denmark)

Neergaard, Ulla; Nielsen, Ruth

2010-01-01

of welfare functions into EU law both from an internal market law and a constitutional law perspective. The main problem areas covered by the Blurring Boundaries project were studied in sub-projects on: 1) Internal market law and welfare services; 2) Fundamental rights and non-discrimination law aspects......; and 3) Services of general interest. In the Blurring Boundaries project, three aspects of the European Social Model have been particularly highlighted: the constitutionalisation of the European Social Model, its multi-level legal character, and the clash between market access justice at EU level...... and distributive justice at national level....

Closed-form solution for piezoelectric layer with two collinear cracks parallel to the boundaries

Directory of Open Access Journals (Sweden)

B. M. Singh

2006-01-01

Full Text Available We consider the problem of determining the stress distribution in an infinitely long piezoelectric layer of finite width, with two collinear cracks of equal length and parallel to the layer boundaries. Within the framework of reigning piezoelectric theory under mode III, the cracked piezoelectric layer subjected to combined electromechanical loading is analyzed. The faces of the layers are subjected to electromechanical loading. The collinear cracks are located at the middle plane of the layer parallel to its face. By the use of Fourier transforms we reduce the problem to solving a set of triple integral equations with cosine kernel and a weight function. The triple integral equations are solved exactly. Closed form analytical expressions for stress intensity factors, electric displacement intensity factors, and shape of crack and energy release rate are derived. As the limiting case, the solution of the problem with one crack in the layer is derived. Some numerical results for the physical quantities are obtained and displayed graphically.

Convergence Analysis of the Preconditioned Group Splitting Methods in Boundary Value Problems

Directory of Open Access Journals (Sweden)

Norhashidah Hj. Mohd Ali

2012-01-01

Full Text Available The construction of a specific splitting-type preconditioner in block formulation applied to a class of group relaxation iterative methods derived from the centred and rotated (skewed finite difference approximations has been shown to improve the convergence rates of these methods. In this paper, we present some theoretical convergence analysis on this preconditioner specifically applied to the linear systems resulted from these group iterative schemes in solving an elliptic boundary value problem. We will theoretically show the relationship between the spectral radiuses of the iteration matrices of the preconditioned methods which affects the rate of convergence of these methods. We will also show that the spectral radius of the preconditioned matrices is smaller than that of their unpreconditioned counterparts if the relaxation parameter is in a certain optimum range. Numerical experiments will also be presented to confirm the agreement between the theoretical and the experimental results.

On Hydroelastic Body-Boundary Condition of Floating Structures

DEFF Research Database (Denmark)

Xia, Jinzhu

1996-01-01

A general linear body boundary condition of hydroelastic analysis of arbitrary shaped floating structures generalizes the classic kinematic rigid-body (Timman-Newman) boundary condition for seakeeping problems. The new boundary condition is consistent with the existing theories under certain...

Numerical experiments using CHIEF to treat the nonuniqueness in solving acoustic axisymmetric exterior problems via boundary integral equations

Directory of Open Access Journals (Sweden)

Adel A.K. Mohsen

2010-07-01

Full Text Available The problem of nonuniqueness (NU of the solution of exterior acoustic problems via boundary integral equations is discussed in this article. The efficient implementation of the CHIEF (Combined Helmholtz Integral Equations Formulation method to axisymmetric problems is studied. Interior axial fields are used to indicate the solution error and to select proper CHIEF points. The procedure makes full use of LU-decomposition as well as the forward solution derived in the solution. Implementations of the procedure for hard spheres are presented. Accurate results are obtained up to a normalised radius of ka = 20.983, using only one CHIEF point. The radiation from a uniformly vibrating sphere is also considered. Accurate results for ka up to 16.927 are obtained using two CHIEF points.

Comments on the comparison of global methods for linear two-point boundary value problems

International Nuclear Information System (INIS)

de Boor, C.; Swartz, B.

1977-01-01

A more careful count of the operations involved in solving the linear system associated with collocation of a two-point boundary value problem using a rough splines reverses results recently reported by others in this journal. In addition, it is observed that the use of the technique of ''condensation of parameters'' can decrease the computer storage required. Furthermore, the use of a particular highly localized basis can also reduce the setup time when the mesh is irregular. Finally, operation counts are roughly estimated for the solution of certain linear system associated with two competing collocation methods; namely, collocation with smooth splines and collocation of the equivalent first order system with continuous piecewise polynomials

Existence and Uniqueness of Solutions for a Discrete Fractional Mixed Type Sum-Difference Equation Boundary Value Problem

Directory of Open Access Journals (Sweden)

Weidong Lv

2015-01-01

Full Text Available By means of Schauder’s fixed point theorem and contraction mapping principle, we establish the existence and uniqueness of solutions to a boundary value problem for a discrete fractional mixed type sum-difference equation with the nonlinear term dependent on a fractional difference of lower order. Moreover, a suitable choice of a Banach space allows the solutions to be unbounded and two representative examples are presented to illustrate the effectiveness of the main results.

Adaptive Problem Solving

Science.gov (United States)

2017-03-01

Borrajo and Raquel Fuentetaja, Universidad Carlos III de Madrid on the meta-level search architecture for finding good combinations of representations and...heuristics on a problem-by-problem basis. The other is with Carlos Linares also from Universidad Carlos III de Madrid on developing effective

Can inequality be tamed through boundary work? A qualitative study of health promotion aimed at reducing health inequalities.

Science.gov (United States)

Pedersen, Pia Vivian; Hjelmar, Ulf; Høybye, Mette Terp; Rod, Morten Hulvej

2017-07-01

This paper examines the organisational dynamics that arise in health promotion aimed at reducing health inequalities. The paper draws on ethnographic fieldwork among public health officers in Danish municipalities and qualitative interviews from an evaluation of health promotion programmes targeting homeless and other marginalised citizens. Analytically, we focus on 'boundary work', i.e. the ways in which social and symbolic boundaries are established, maintained, transgressed and negotiated, both at the administrative level and among frontline professionals. The paper discusses three types of boundary work: (i) demarcating professional domains; (ii) setting the boundaries of the task itself; and (iii) managing administrative boundaries. The main argument is that the production, maintenance and transgression of these three types of boundaries constitute central and time-consuming aspects of the practices of public health professionals, and that boundary work constitutes an important element in professional practices seeking to 'tame a wicked problem', such as social inequalities in health. A cross-cutting feature of the three types of boundary work is the management of the divide between health and social issues, which the professionals seemingly seek to uphold and transgress at the same time. The paper thus contributes to ongoing discussions of intersectoral action to address health inequalities. Furthermore, it extends the scope and application of the concept of boundary work in the sociology of public health by suggesting that the focus in previous research on professional demarcation be broadened in order to capture other types of boundaries that shape, and are shaped by, professional practices. Copyright © 2017 Elsevier Ltd. All rights reserved.

A review on application of MHD theory to plasma boundary problems in tokamaks

International Nuclear Information System (INIS)

Itoh, Kimitaka.

1992-08-01

A survey is made on the problems of the edge plasmas, to which the analyses based on the MHD theory have been successfully applied. Also discussed are the efforts to extend the model equation to more general (and important as well) problems such as H-mode physics. An overview is first made on the advantages of the MHD picture, and the necessary supplementary physics are examined. Next, one- and two-dimensional models of the spatial structure of the edge plasma is discussed. The results on the stationary structure, both analytical and numerical, are reviewed: Typical example as well as the scaling law are shown. The instabilities associated with edge plasma is next reviewed. The surface kink mode, ballooning mode, interchange mode, resistive interchange mode and thermal instability are discussed. Role of the geometry such as the location of the X-point is studied. Influences of the atomic processes, and those of the radial electric field are also discussed. The analysis of the H-mode transition physics is finally discussed. The boundary plasma is a nonlinear media which possesses the possibility for bifurcation in which the radial electric field plays a key role. The model of the ion viscosity is also studied. Transition physics is developed. Analysis on the self-generating oscillation is shown and the relation with ELMs is discussed. After reviewing these problems, several comments are made to what directions the study can be deepened. (author) 53 refs

Solvability, regularity, and optimal control of boundary value problems for pdes in honour of Prof. Gianni Gilardi

CERN Document Server

Favini, Angelo; Rocca, Elisabetta; Schimperna, Giulio; Sprekels, Jürgen

2017-01-01

This volume gathers contributions in the field of partial differential equations, with a focus on mathematical models in phase transitions, complex fluids and thermomechanics. These contributions are dedicated to Professor Gianni Gilardi on the occasion of his 70th birthday. It particularly develops the following thematic areas: nonlinear dynamic and stationary equations; well-posedness of initial and boundary value problems for systems of PDEs; regularity properties for the solutions; optimal control problems and optimality conditions; feedback stabilization and stability results. Most of the articles are presented in a self-contained manner, and describe new achievements and/or the state of the art in their line of research, providing interested readers with an overview of recent advances and future research directions in PDEs.

The stability boundary of group-III transition metal diboride ScB 2 (0 0 0 1) surfaces

Science.gov (United States)

Zhao, Hui; Qin, Na

2012-01-01

Experimental observations and theoretical investigations exhibit that a group-IV(V) transition metal diboride (0 0 0 1) surface is terminated with a 1 × 1 TM(B) layer. As to a group-III transition metal diboride, we have investigated the stability boundary of ScB2 (0 0 0 1) surfaces using first principles total energy plane-wave pseudopotential method based on density functional theory. The Mulliken charge population analysis shows that Sc atoms in the second layer cannot provide B atoms in the first layer with sufficient electrons to form a complete graphene-like boron layer. We also found that the charge transfer between the first and the second layer for the B-terminated surface is more than that for Sc-terminated surface. It elucidates the reason that the outermost interlayer spacing contract more strongly in the B-terminated surface than in the Sc-terminated surface. The surface energies of both terminated ScB2 (0 0 0 1) surfaces as a function of the chemical potential of B are also calculated to check the relative stability of the two surface structures.

Quantum Gravitational Effects on the Boundary

Science.gov (United States)

James, F.; Park, I. Y.

2018-04-01

Quantum gravitational effects might hold the key to some of the outstanding problems in theoretical physics. We analyze the perturbative quantum effects on the boundary of a gravitational system and the Dirichlet boundary condition imposed at the classical level. Our analysis reveals that for a black hole solution, there is a contradiction between the quantum effects and the Dirichlet boundary condition: the black hole solution of the one-particle-irreducible action no longer satisfies the Dirichlet boundary condition as would be expected without going into details. The analysis also suggests that the tension between the Dirichlet boundary condition and loop effects is connected with a certain mechanism of information storage on the boundary.

Finite-volume discretizations and immersed boundaries

NARCIS (Netherlands)

Y.J. Hassen (Yunus); B. Koren (Barry)

2009-01-01

htmlabstractIn this chapter, an accurate method, using a novel immersed-boundary approach, is presented for numerically solving linear, scalar convection problems. As is standard in immersed-boundary methods, moving bodies are embedded in a fixed `Cartesian' grid. The essence of the present method

Finite-volume discretizations and immersed boundaries

NARCIS (Netherlands)

Y.J. Hassen (Yunus); B. Koren (Barry)

2010-01-01

textabstractIn this chapter, an accurate method, using a novel immersed-boundary approach, is presented for numerically solving linear, scalar convection problems. As is standard in immersed-boundary methods, moving bodies are embedded in a fixed Cartesian grid. The essence of the present method is

Subgroup report on grain boundary and interphase boundary structure and properties

International Nuclear Information System (INIS)

Balluffi, R.W.; Cannon, R.M.; Clarke, D.R.; Heuer, A.H.; Ho, P.S.; Kear, B.H.; Vitek, V.; Weertman, J.R.; White, C.L.

1979-01-01

In many high temperature structural applications, the performance characteristics of a materials system are largely controlled by the properties of its grain and interphase boundaries. Failure in creep and fatigue frequently occurs by cavitation, or cracking along grain boundaries. In a few special cases, this failure problem has been overcome by directional alignment of grain and interphase boundaries by various types of metallurgical processing such as directional solidification and directional recrystallization. A good example is to be found in the application of directionally aligned structures in high performance gas-turbine airfoils. However, where fine, equiaxed grain structures are desirable, other methods of controlling grain boundary properties have been developed. Important among these has been the introduction of improvements in primary melting practices, designed to control important impurities. This is of decisive importance because even traces of certain impurity elements present in grain boundaries in high temperature materials can seriously affect properties. Impurities are deleterious and need to be removed. However, in certain cases, (e.g., creep fracture) controlled impurity additions can be beneficial and result in improved properties

Communication Reducing Algorithms for Distributed Hierarchical N-Body Problems with Boundary Distributions

KAUST Repository

AbdulJabbar, Mustafa Abdulmajeed

2017-05-11

Reduction of communication and efficient partitioning are key issues for achieving scalability in hierarchical N-Body algorithms like Fast Multipole Method (FMM). In the present work, we propose three independent strategies to improve partitioning and reduce communication. First, we show that the conventional wisdom of using space-filling curve partitioning may not work well for boundary integral problems, which constitute a significant portion of FMM’s application user base. We propose an alternative method that modifies orthogonal recursive bisection to relieve the cell-partition misalignment that has kept it from scaling previously. Secondly, we optimize the granularity of communication to find the optimal balance between a bulk-synchronous collective communication of the local essential tree and an RDMA per task per cell. Finally, we take the dynamic sparse data exchange proposed by Hoefler et al. [1] and extend it to a hierarchical sparse data exchange, which is demonstrated at scale to be faster than the MPI library’s MPI_Alltoallv that is commonly used.

Numerical continuation methods for dynamical systems path following and boundary value problems

CERN Document Server

Krauskopf, Bernd; Galan-Vioque, Jorge

2007-01-01

Path following in combination with boundary value problem solvers has emerged as a continuing and strong influence in the development of dynamical systems theory and its application. It is widely acknowledged that the software package AUTO - developed by Eusebius J. Doedel about thirty years ago and further expanded and developed ever since - plays a central role in the brief history of numerical continuation. This book has been compiled on the occasion of Sebius Doedel''s 60th birthday. Bringing together for the first time a large amount of material in a single, accessible source, it is hoped that the book will become the natural entry point for researchers in diverse disciplines who wish to learn what numerical continuation techniques can achieve. The book opens with a foreword by Herbert B. Keller and lecture notes by Sebius Doedel himself that introduce the basic concepts of numerical bifurcation analysis. The other chapters by leading experts discuss continuation for various types of systems and objects ...

Necessary and Sufficient Conditions for the Existence of Positive Solution for Singular Boundary Value Problems on Time Scales

Directory of Open Access Journals (Sweden)

Zhang Xuemei

2009-01-01

Full Text Available By constructing available upper and lower solutions and combining the Schauder's fixed point theorem with maximum principle, this paper establishes sufficient and necessary conditions to guarantee the existence of as well as positive solutions for a class of singular boundary value problems on time scales. The results significantly extend and improve many known results for both the continuous case and more general time scales. We illustrate our results by one example.

The steady Navier–Stokes problem with the inhomogeneous Navier-type boundary conditions in a 2D multiply-connected bounded domain

Czech Academy of Sciences Publication Activity Database

Neustupa, Jiří

2015-01-01

Roč. 35, č. 3 (2015), s. 201-212 ISSN 0174-4747 R&D Projects: GA ČR GA13-00522S Institutional support: RVO:67985840 Keywords : steady Navier-Stokes problem * slip boundary conditions Subject RIV: BA - General Mathematics http://www.degruyter.com/view/j/anly.2015.35.issue-3/anly-2014-1304/anly-2014-1304. xml

An algorithm for full parametric solution of problems on the statics of orthotropic plates by the method of boundary states with perturbations

Science.gov (United States)

Penkov, V. B.; Ivanychev, D. A.; Novikova, O. S.; Levina, L. V.

2018-03-01

The article substantiates the possibility of building full parametric analytical solutions of mathematical physics problems in arbitrary regions by means of computer systems. The suggested effective means for such solutions is the method of boundary states with perturbations, which aptly incorporates all parameters of an orthotropic medium in a general solution. We performed check calculations of elastic fields of an anisotropic rectangular region (test and calculation problems) for a generalized plane stress state.

Vortex sheet approximation of boundary layers

International Nuclear Information System (INIS)

Chorin, A.J.

1978-01-01

a grid free method for approximating incomprssible boundary layers is introduced. The computational elements are segments of vortex sheets. The method is related to the earlier vortex method; simplicity is achieved at the cost of replacing the Navier-Stokes equations by the Prandtl boundary layer equations. A new method for generating vorticity at boundaries is also presented; it can be used with the earlier voartex method. The applications presented include (i) flat plate problems, and (ii) a flow problem in a model cylinder- piston assembly, where the new method is used near walls and an improved version of the random choice method is used in the interior. One of the attractive features of the new method is the ease with which it can be incorporated into hybrid algorithms

Unique solvability of a non-linear non-local boundary-value problem for systems of non-linear functional differential equations

Czech Academy of Sciences Publication Activity Database

Dilna, N.; Rontó, András

2010-01-01

Roč. 60, č. 3 (2010), s. 327-338 ISSN 0139-9918 R&D Projects: GA ČR(CZ) GA201/06/0254 Institutional research plan: CEZ:AV0Z10190503 Keywords : non-linear boundary value-problem * functional differential equation * non-local condition * unique solvability * differential inequality Subject RIV: BA - General Mathematics Impact factor: 0.316, year: 2010 http://link.springer.com/article/10.2478%2Fs12175-010-0015-9

A One-Dimensional Wave Equation with White Noise Boundary Condition

International Nuclear Information System (INIS)

Kim, Jong Uhn

2006-01-01

We discuss the Cauchy problem for a one-dimensional wave equation with white noise boundary condition. We also establish the existence of an invariant measure when the noise is additive. Similar problems for parabolic equations were discussed by several authors. To our knowledge, there is only one work which investigated the initial-boundary value problem for a wave equation with random noise at the boundary. We handle a more general case by a different method. Our result on the existence of an invariant measure relies on the author's recent work on a certain class of stochastic evolution equations

Multiple positive solutions to nonlinear boundary value problems of a system for fractional differential equations.

Science.gov (United States)

Zhai, Chengbo; Hao, Mengru

2014-01-01

By using Krasnoselskii's fixed point theorem, we study the existence of at least one or two positive solutions to a system of fractional boundary value problems given by -D(0+)(ν1)y1(t) = λ1a1(t)f(y1(t), y2(t)), - D(0+)(ν2)y2(t) = λ2a2(t)g(y1(t), y2(t)), where D(0+)(ν) is the standard Riemann-Liouville fractional derivative, ν1, ν2 ∈ (n - 1, n] for n > 3 and n ∈ N, subject to the boundary conditions y1((i))(0) = 0 = y ((i))(0), for 0 ≤ i ≤ n - 2, and [D(0+)(α)y1(t)] t=1 = 0 = [D(0+ (α)y2(t)] t=1, for 1 ≤ α ≤ n - 2, or y1((i))(0) = 0 = y ((i))(0), for 0 ≤ i ≤ n - 2, and [D(0+)(α)y1(t)] t=1 = ϕ1(y1), [D(0+)(α)y2(t)] t=1 = ϕ2(y2), for 1 ≤ α ≤ n - 2, ϕ1, ϕ2 ∈ C([0,1], R). Our results are new and complement previously known results. As an application, we also give an example to demonstrate our result.

Numerical analysis of the asymptotic behavior of solutions of a boundary problem for a nonlinear parabolic equation

International Nuclear Information System (INIS)

Vasileva, D.P.

1993-01-01

Blow-up and global time self-similar solutions of a boundary problem for a nonlinear equation u t = Δ u σ+1 + u β are found in the case β = σ + 1. It is shown that they describe the asymptotic behavior of a wide class of initial perturbations. A numerical investigation of the solutions in the case β>σ + 1 is also made. A hypothesis is done that the behavior for large times of global time solutions is described by the self-similar solutions of the equation without source.(author). 20 refs.; 9 figs

An advanced boundary element method (BEM) implementation for the forward problem of electromagnetic source imaging

International Nuclear Information System (INIS)

Akalin-Acar, Zeynep; Gencer, Nevzat G

2004-01-01

The forward problem of electromagnetic source imaging has two components: a numerical model to solve the related integral equations and a model of the head geometry. This study is on the boundary element method (BEM) implementation for numerical solutions and realistic head modelling. The use of second-order (quadratic) isoparametric elements and the recursive integration technique increase the accuracy in the solutions. Two new formulations are developed for the calculation of the transfer matrices to obtain the potential and magnetic field patterns using realistic head models. The formulations incorporate the use of the isolated problem approach for increased accuracy in solutions. If a personal computer is used for computations, each transfer matrix is calculated in 2.2 h. After this pre-computation period, solutions for arbitrary source configurations can be obtained in milliseconds for a realistic head model. A hybrid algorithm that uses snakes, morphological operations, region growing and thresholding is used for segmentation. The scalp, skull, grey matter, white matter and eyes are segmented from the multimodal magnetic resonance images and meshes for the corresponding surfaces are created. A mesh generation algorithm is developed for modelling the intersecting tissue compartments, such as eyes. To obtain more accurate results quadratic elements are used in the realistic meshes. The resultant BEM implementation provides more accurate forward problem solutions and more efficient calculations. Thus it can be the firm basis of the future inverse problem solutions

Causal boundary for stably causal space-times

International Nuclear Information System (INIS)

Racz, I.

1987-12-01

The usual boundary constructions for space-times often yield an unsatisfactory boundary set. This problem is reviewed and a new solution is proposed. An explicit identification rule is given on the set of the ideal points of the space-time. This construction leads to a satisfactory boundary point set structure for stably causal space-times. The topological properties of the resulting causal boundary construction are examined. For the stably causal space-times each causal curve has a unique endpoint on the boundary set according to the extended Alexandrov topology. The extension of the space-time through the boundary is discussed. To describe the singularities the defined boundary sets have to be separated into two disjoint sets. (D.Gy.) 8 refs

Polyakov's quantized string with boundary terms

International Nuclear Information System (INIS)

Durhuus, B.; Olesen, P.; Petersen, J.L.

1981-11-01

The authors compute the boundary terms needed in Polyakov's method for calculating averages of functionals defined on surfaces. The method used is due to Seeley, who found recursive relations yielding the boundary terms. These relations are solved for a general second order elliptic differential operator. This solution is then applied to Polyakov's problem. (Auth.)

An iterative kernel based method for fourth order nonlinear equation with nonlinear boundary condition

Science.gov (United States)

Azarnavid, Babak; Parand, Kourosh; Abbasbandy, Saeid

2018-06-01

This article discusses an iterative reproducing kernel method with respect to its effectiveness and capability of solving a fourth-order boundary value problem with nonlinear boundary conditions modeling beams on elastic foundations. Since there is no method of obtaining reproducing kernel which satisfies nonlinear boundary conditions, the standard reproducing kernel methods cannot be used directly to solve boundary value problems with nonlinear boundary conditions as there is no knowledge about the existence and uniqueness of the solution. The aim of this paper is, therefore, to construct an iterative method by the use of a combination of reproducing kernel Hilbert space method and a shooting-like technique to solve the mentioned problems. Error estimation for reproducing kernel Hilbert space methods for nonlinear boundary value problems have yet to be discussed in the literature. In this paper, we present error estimation for the reproducing kernel method to solve nonlinear boundary value problems probably for the first time. Some numerical results are given out to demonstrate the applicability of the method.

Boundary value problems of finite elasticity local theorems on existence, uniqueness, and analytic dependence on data

CERN Document Server

Valent, Tullio

1988-01-01

In this book I present, in a systematic form, some local theorems on existence, uniqueness, and analytic dependence on the load, which I have recently obtained for some types of boundary value problems of finite elasticity. Actually, these results concern an n-dimensional (n ~ 1) formal generalization of three-dimensional elasticity. Such a generalization, be­ sides being quite spontaneous, allows us to consider a great many inter­ esting mathematical situations, and sometimes allows us to clarify certain aspects of the three-dimensional case. Part of the matter presented is unpublished; other arguments have been only partially published and in lesser generality. Note that I concentrate on simultaneous local existence and uniqueness; thus, I do not deal with the more general theory of exis­ tence. Moreover, I restrict my discussion to compressible elastic bodies and I do not treat unilateral problems. The clever use of the inverse function theorem in finite elasticity made by STOPPELLI [1954, 1957a, 1957b]...

Variational principle for gravity with null and non-null boundaries: a unified boundary counter-term

Energy Technology Data Exchange (ETDEWEB)

Parattu, Krishnamohan; Chakraborty, Sumanta; Padmanabhan, T. [IUCAA, Post Bag 4, Pune (India)

2016-03-15

It is common knowledge that the Einstein-Hilbert action does not furnish a well-posed variational principle. The usual solution to this problem is to add an extra boundary term to the action, called a counter-term, so that the variational principle becomes well-posed. When the boundary is spacelike or timelike, the Gibbons-Hawking-York counter-term is the most widely used. For null boundaries, we had proposed a counter-term in a previous paper. In this paper, we extend the previous analysis and propose a counter-term that can be used to eliminate variations of the ''off-the-surface'' derivatives of the metric on any boundary, regardless of its spacelike, timelike or null nature. (orig.)

Regularity of pointwise boundary control systems

DEFF Research Database (Denmark)

Pedersen, Michael

1992-01-01

We will in these notes address some problems arising in "real-life" control application, namely problems concerning distributional control inputs on the boundary of the spatial domain. We extend the classical variational approach and give easily checkable sufficient conditions for the solutions...

Singular boundary perturbations of distributed systems

DEFF Research Database (Denmark)

Pedersen, Michael

1990-01-01

Some problems arising in real-life control applications are addressed--namely, problems concerning non-smooth control inputs on the boundary of the spatial domain. The classical variational approach is extended, and sufficient conditions are given for the solutions to continuous functions of time...

Boundary conditions for the numerical solution of elliptic equations in exterior regions

International Nuclear Information System (INIS)

Bayliss, A.; Gunzburger, M.; Turkel, E.

1982-01-01

Elliptic equations in exterior regions frequently require a boundary condition at infinity to ensure the well-posedness of the problem. Examples of practical applications include the Helmholtz equation and Laplace's equation. Computational procedures based on a direct discretization of the elliptic problem require the replacement of the condition at infinity by a boundary condition on a finite artificial surface. Direct imposition of the condition at infinity along the finite boundary results in large errors. A sequence of boundary conditions is developed which provides increasingly accurate approximations to the problem in the infinite domain. Estimates of the error due to the finite boundary are obtained for several cases. Computations are presented which demonstrate the increased accuracy that can be obtained by the use of the higher order boundary conditions. The examples are based on a finite element formulation but finite difference methods can also be used

Open boundary condition, Wilson flow and the scalar glueball mass

International Nuclear Information System (INIS)

Chowdhury, Abhishek; Harindranath, A.; Maiti, Jyotirmoy

2014-01-01

A major problem with periodic boundary condition on the gauge fields used in current lattice gauge theory simulations is the trapping of topological charge in a particular sector as the continuum limit is approached. To overcome this problem open boundary condition in the temporal direction has been proposed recently. One may ask whether open boundary condition can reproduce the observables calculated with periodic boundary condition. In this work we find that the extracted lowest glueball mass using open and periodic boundary conditions at the same lattice volume and lattice spacing agree for the range of lattice scales explored in the range 3 GeV≤(1/a)≤5 GeV. The problem of trapping is overcome to a large extent with open boundary and we are able to extract the glueball mass at even larger lattice scale ≈ 5.7 GeV. To smoothen the gauge fields we have used recently proposed Wilson flow which, compared to HYP smearing, exhibits better systematics in the extraction of glueball mass. The extracted glueball mass shows remarkable insensitivity to the lattice spacings in the range explored in this work, 3 GeV≤(1/a)≤5.7 GeV.

Internal and external 2-d boundary layer flows

Science.gov (United States)

Crawford, M. E.; Kays, W. M.

1978-01-01

Computer program computes general two dimensional turbulent boundary-layer flow using finite-difference techniques. Structure allows for user modification to accommodate unique problems. Program should prove useful in many applications where accurate boundary-layer flow calculations are required.

Multiple Solutions of Nonlinear Boundary Value Problems of Fractional Order: A New Analytic Iterative Technique

Directory of Open Access Journals (Sweden)

Omar Abu Arqub

2014-01-01

Full Text Available The purpose of this paper is to present a new kind of analytical method, the so-called residual power series, to predict and represent the multiplicity of solutions to nonlinear boundary value problems of fractional order. The present method is capable of calculating all branches of solutions simultaneously, even if these multiple solutions are very close and thus rather difficult to distinguish even by numerical techniques. To verify the computational efficiency of the designed proposed technique, two nonlinear models are performed, one of them arises in mixed convection flows and the other one arises in heat transfer, which both admit multiple solutions. The results reveal that the method is very effective, straightforward, and powerful for formulating these multiple solutions.

Boundary value problems for multi-term fractional differential equations

Science.gov (United States)

Daftardar-Gejji, Varsha; Bhalekar, Sachin

2008-09-01

Multi-term fractional diffusion-wave equation along with the homogeneous/non-homogeneous boundary conditions has been solved using the method of separation of variables. It is observed that, unlike in the one term case, solution of multi-term fractional diffusion-wave equation is not necessarily non-negative, and hence does not represent anomalous diffusion of any kind.

E-coil: an inverse boundary element method for a quasi-static problem

International Nuclear Information System (INIS)

Sanchez, Clemente Cobos; Garcia, Salvador Gonzalez; Power, Henry

2010-01-01

Boundary element methods represent a valuable approach for designing gradient coils; these methods are based on meshing the current carrying surface into an array of boundary elements. The temporally varying magnetic fields produced by gradient coils induce electric currents in conducting tissues and so the exposure of human subjects to these magnetic fields has become a safety concern, especially with the increase in the strength of the field gradients used in magnetic resonance imaging. Here we present a boundary element method for the design of coils that minimize the electric field induced in prescribed conducting systems. This work also details some numerical examples of the application of this coil design method. The reduction of the electric field induced in a prescribed region inside the coils is also evaluated.

E-coil: an inverse boundary element method for a quasi-static problem

Energy Technology Data Exchange (ETDEWEB)

Sanchez, Clemente Cobos; Garcia, Salvador Gonzalez [Depto. Electromagnetismo y F. de la Materia Facultad de Ciencias University of Granada Avda. Fuentenueva E-18071 (Spain); Power, Henry, E-mail: ccobos@ugr.e [School of Mechanical, Materials and Manufacturing Engineering, The University of Nottingham, Nottingham Park, Nottingham NG7 2RD (United Kingdom)

2010-06-07

Boundary element methods represent a valuable approach for designing gradient coils; these methods are based on meshing the current carrying surface into an array of boundary elements. The temporally varying magnetic fields produced by gradient coils induce electric currents in conducting tissues and so the exposure of human subjects to these magnetic fields has become a safety concern, especially with the increase in the strength of the field gradients used in magnetic resonance imaging. Here we present a boundary element method for the design of coils that minimize the electric field induced in prescribed conducting systems. This work also details some numerical examples of the application of this coil design method. The reduction of the electric field induced in a prescribed region inside the coils is also evaluated.

Accurate Solution of Multi-Region Continuum Biomolecule Electrostatic Problems Using the Linearized Poisson-Boltzmann Equation with Curved Boundary Elements

Science.gov (United States)

Altman, Michael D.; Bardhan, Jaydeep P.; White, Jacob K.; Tidor, Bruce

2009-01-01

We present a boundary-element method (BEM) implementation for accurately solving problems in biomolecular electrostatics using the linearized Poisson–Boltzmann equation. Motivating this implementation is the desire to create a solver capable of precisely describing the geometries and topologies prevalent in continuum models of biological molecules. This implementation is enabled by the synthesis of four technologies developed or implemented specifically for this work. First, molecular and accessible surfaces used to describe dielectric and ion-exclusion boundaries were discretized with curved boundary elements that faithfully reproduce molecular geometries. Second, we avoided explicitly forming the dense BEM matrices and instead solved the linear systems with a preconditioned iterative method (GMRES), using a matrix compression algorithm (FFTSVD) to accelerate matrix-vector multiplication. Third, robust numerical integration methods were employed to accurately evaluate singular and near-singular integrals over the curved boundary elements. Finally, we present a general boundary-integral approach capable of modeling an arbitrary number of embedded homogeneous dielectric regions with differing dielectric constants, possible salt treatment, and point charges. A comparison of the presented BEM implementation and standard finite-difference techniques demonstrates that for certain classes of electrostatic calculations, such as determining absolute electrostatic solvation and rigid-binding free energies, the improved convergence properties of the BEM approach can have a significant impact on computed energetics. We also demonstrate that the improved accuracy offered by the curved-element BEM is important when more sophisticated techniques, such as non-rigid-binding models, are used to compute the relative electrostatic effects of molecular modifications. In addition, we show that electrostatic calculations requiring multiple solves using the same molecular geometry

Polyakov's quantized string with boundary terms

International Nuclear Information System (INIS)

Durhuus, B.; Olesen, P.; Petersen, J.L.

1982-01-01

We compute the boundary terms due to the conformal anomaly which are needed in Polyakov's method of calculating averages of functionals defined on surfaces. The method we use is due to Seeley, who found recursive relations yielding the boundary terms. We solve these relations for a general second-order elliptic differential operator. This solution is then applied to Polyakov's problem. (orig.)

Absorbing boundary conditions for Einstein's field equations

Energy Technology Data Exchange (ETDEWEB)

Sarbach, Olivier [Instituto de Fisica y Matematicas, Universidad Michoacana de San Nicolas de Hidalgo, Edificio C-3, Cd. Universitaria. C. P. 58040 Morelia, Michoacan (Mexico)

2007-11-15

A common approach for the numerical simulation of wave propagation on a spatially unbounded domain is to truncate the domain via an artificial boundary, thus forming a finite computational domain with an outer boundary. Absorbing boundary conditions must then be specified at the boundary such that the resulting initial-boundary value problem is well posed and such that the amount of spurious reflection is minimized. In this article, we review recent results on the construction of absorbing boundary conditions in General Relativity and their application to numerical relativity.

Beyond the numbers: expanding the boundaries of neuropsychology.

Science.gov (United States)

Perry, William

2009-02-01

Beyond the Numbers: Expanding the Boundaries of Neuropsychology was Dr Perry's 2007 presidential address in the annual conference of the National Academy of Neuropsychology. In his address he discussed the achievements of the science of neuropsychology and highlighted some areas that exemplified the expansion of the boundaries of neuropsychology. These areas are: (i) the study of neuropsychological functioning in new or non-traditional populations, particularly seemingly healthy people and people with non-brain diseases; (ii) the interface of cognition and genetics; (iii) the use of the process approach as a means of understanding brain functioning; and (iv) a translational application to the science of neuropsychology.

Beyond the Numbers: Expanding the Boundaries of Neuropsychology†

Science.gov (United States)

Perry, William

2009-01-01

Beyond the Numbers: Expanding the Boundaries of Neuropsychology was Dr Perry's 2007 presidential address in the annual conference of the National Academy of Neuropsychology. In his address he discussed the achievements of the science of neuropsychology and highlighted some areas that exemplified the expansion of the boundaries of neuropsychology. These areas are: (i) the study of neuropsychological functioning in new or non-traditional populations, particularly seemingly healthy people and people with non-brain diseases; (ii) the interface of cognition and genetics; (iii) the use of the process approach as a means of understanding brain functioning; and (iv) a translational application to the science of neuropsychology. PMID:19395354

Topography and instability of monolayers near domain boundaries

International Nuclear Information System (INIS)

Diamant, H.; Witten, T. A.; Ege, C.; Gopal, A.; Lee, K. Y. C.

2001-01-01

We theoretically study the topography of a biphasic surfactant monolayer in the vicinity of domain boundaries. The differing elastic properties of the two phases generally lead to a nonflat topography of 'mesas,' where domains of one phase are elevated with respect to the other phase. The mesas are steep but low, having heights of up to 10 nm. As the monolayer is laterally compressed, the mesas develop overhangs and eventually become unstable at a surface tension of about K(δc 0 ) 2 (δc 0 being the difference in spontaneous curvature and K a bending modulus). In addition, the boundary is found to undergo a topography-induced rippling instability upon compression, if its line tension is smaller than about Kδc 0 . The effect of diffuse boundaries on these features and the topographic behavior near a critical point are also examined. We discuss the relevance of our findings to several experimental observations related to surfactant monolayers: (i) small topographic features recently found near domain boundaries; (ii) folding behavior observed in mixed phospholipid monolayers and model lung surfactants; (iii) roughening of domain boundaries seen under lateral compression; (iv) the absence of biphasic structures in tensionless surfactant films

Boundary Shape Control of the Navier-Stokes Equations and Applications

Institute of Scientific and Technical Information of China (English)

Kaitai LI; Jian SU; Aixiang HUANG

2010-01-01

In this paper,the geometrical design for the blade's surface(s)in an impeller or for the profile of an aircraft,is modeled from the mathematical point of view by a boundary shape control problem for the Navier-Stokes equations.The objective function is the sum of a global dissipative function and the power of the fluid.The control variables are the geometry of the boundary and the state equations are the Navier-Stokes equations.The Euler-Lagrange equations of the optimal control problem are derived,which are an elliptic boundary value system of fourth order,coupled with the Navier-Stokes equations.The authors also prove the existence of the solution of the optimal control problem,the existence of the solution of the Navier-Stokes equations with mixed boundary conditions,the weak continuity of the solution of the Navier-Stokes equations with respect to the geometry shape of the blade's surface and the existence of solutions of the equations for the G(a)teaux derivative of the solution of the Navier-Stokes equations with respect to the geometry of the boundary.

III-V microelectronics

CERN Document Server

Nougier, JP

1991-01-01

As is well known, Silicon widely dominates the market of semiconductor devices and circuits, and in particular is well suited for Ultra Large Scale Integration processes. However, a number of III-V compound semiconductor devices and circuits have recently been built, and the contributions in this volume are devoted to those types of materials, which offer a number of interesting properties. Taking into account the great variety of problems encountered and of their mutual correlations when fabricating a circuit or even a device, most of the aspects of III-V microelectronics, from fundamental p

Stabilizing local boundary conditions for two-dimensional shallow water equations

KAUST Repository

Dia, Ben Mansour

2018-03-27

In this article, we present a sub-critical two-dimensional shallow water flow regulation. From the energy estimate of a set of one-dimensional boundary stabilization problems, we obtain a set of polynomial equations with respect to the boundary values as a requirement for the energy decrease. Using the Riemann invariant analysis, we build stabilizing local boundary conditions that guarantee the stability of the hydrodynamical state around a given steady state. Numerical results for the controller applied to the nonlinear problem demonstrate the performance of the method.

Laplace Boundary-Value Problem in Paraboloidal Coordinates

Science.gov (United States)

Duggen, L.; Willatzen, M.; Voon, L. C. Lew Yan

2012-01-01

This paper illustrates both a problem in mathematical physics, whereby the method of separation of variables, while applicable, leads to three ordinary differential equations that remain fully coupled via two separation constants and a five-term recurrence relation for series solutions, and an exactly solvable problem in electrostatics, as a…

A free boundary problem describing the saturated-unsaturated flow in a porous medium

Directory of Open Access Journals (Sweden)

Gabriela Marinoschi

2004-01-01

Full Text Available This paper presents a functional approach to a nonlinear model describing the complete physical process of water infiltration into an unsaturated soil, including the saturation occurrence and the advance of the wetting front. The model introduced in this paper involves a multivalued operator covering the simultaneous saturated and unsaturated flow behaviors and enhances the study of the displacement of the free boundary between these two flow regimes. The model resides in Richards' equation written in pressure form with an initial condition and boundary conditions which in this work express the inflow due to the rain on the soil surface on the one hand, and characterize a certain permeability corresponding to the underground boundary, on the other hand. Existence, uniqueness, and regularity results for the transformed model in diffusive form, that is, for the moisture of the soil, and the existence of the weak solution for the pressure form are proved in the 3D case. The main part of the paper focuses on the existence of the free boundary between the saturated and unsaturated parts of the soil, and this is proved, in the 1D case, for certain stronger assumptions on the initial data and boundary conditions.

CFD simulation of the atmospheric boundary layer: wall function problems

NARCIS (Netherlands)

Blocken, B.J.E.; Stathopoulos, T.; Carmeliet, J.

2007-01-01

Accurate Computational Fluid Dynamics (CFD) simulations of atmospheric boundary layer (ABL) flow are essential for a wide variety of atmospheric studies including pollutant dispersion and deposition. The accuracy of such simulations can be seriously compromised when wall-function roughness

Necessary and Sufficient Conditions for the Existence of Positive Solution for Singular Boundary Value Problems on Time Scales

Directory of Open Access Journals (Sweden)

Meiqiang Feng

2009-01-01

Full Text Available By constructing available upper and lower solutions and combining the Schauder's fixed point theorem with maximum principle, this paper establishes sufficient and necessary conditions to guarantee the existence of Cld[0,1]𝕋 as well as CldΔ[0,1]𝕋 positive solutions for a class of singular boundary value problems on time scales. The results significantly extend and improve many known results for both the continuous case and more general time scales. We illustrate our results by one example.

A free boundary problem for a reaction-diffusion system with nonlinear memory

DEFF Research Database (Denmark)

Lin, Zhigui; Ling, Zhi; Pedersen, Michael

2013-01-01

We consider a integro-partial differential equation with a free boundary which appears in the theory of the nuclear dynamics. First, local existence and uniqueness are obtained by using the contraction mapping theorem. Then, the behavior of the free boundary and the blow-up criteria are obtained........ Finally, we examine the long-time behavior of the global solution. We show that the solution is global and fast if the initial data are small....

Grain boundary structure and properties

International Nuclear Information System (INIS)

Balluffi, R.W.

1979-05-01

An attempt is made to distinguish those fundamental aspects of grain boundaries which should be relevant to the problem of the time dependent fracture of high temperature structural materials. These include the basic phenomena which are thought to be associated with cavitation and cracking at grain boundaries during service and with the more general microstructural changes which occur during both processing and service. A very brief discussion of the current state of knowledge of these fundamentals is given

Multiplicity results for classes of one-dimensional p-Laplacian boundary-value problems with cubic-like nonlinearities

Directory of Open Access Journals (Sweden)

Idris Addou

2000-07-01

Full Text Available We study boundary-value problems of the type $$displaylines{ -(varphi_{p}( u' ' =lambda f( u ,hbox{ in }(0,1 cr u( 0 =u( 1 =0, }$$ where $p>1$, $varphi_{p}( x =left| x ight| ^{p-2}x$, and $lambda >0$. We provide multiplicity results when $f$ behaves like a cubic with three distinct roots, at which it satisfies Lipschitz-type conditions involving a parameter $q>1$. We shall show how changes in the position of $q$ with respect to $p$ lead to different behavior of the solution set. When dealing with sign-changing solutions, we assume that $f$ is {it half-odd}; a condition generalizing the usual oddness. We use a quadrature method.

The laminar boundary layer equations

CERN Document Server

Curle, N

2017-01-01

Thorough introduction to boundary layer problems offers an ordered, logical presentation accessible to undergraduates. The text's careful expositions of the limitations and accuracy of various methods will also benefit professionals. 1962 edition.

Boundary layer and fundamental problems of hydrodynamics (compatibility of a logarithmic velocity profile in a turbulent boundary layer with the experience values)

Science.gov (United States)

Zaryankin, A. E.

2017-11-01

The compatibility of the semiempirical turbulence theory of L. Prandtl with the actual flow pattern in a turbulent boundary layer is considered in this article, and the final calculation results of the boundary layer is analyzed based on the mentioned theory. It shows that accepted additional conditions and relationships, which integrate the differential equation of L. Prandtl, associating the turbulent stresses in the boundary layer with the transverse velocity gradient, are fulfilled only in the near-wall region where the mentioned equation loses meaning and are inconsistent with the physical meaning on the main part of integration. It is noted that an introduced concept about the presence of a laminar sublayer between the wall and the turbulent boundary layer is the way of making of a physical meaning to the logarithmic velocity profile, and can be defined as adjustment of the actual flow to the formula that is inconsistent with the actual boundary conditions. It shows that coincidence of the experimental data with the actual logarithmic profile is obtained as a result of the use of not particular physical value, as an argument, but function of this value.

The zonal satellite problem. III Symmetries

Directory of Open Access Journals (Sweden)

Mioc V.

2002-01-01

Full Text Available The two-body problem associated with a force field described by a potential of the form U =Sum(k=1,n ak/rk (r = distance between particles, ak = real parameters is resumed from the only standpoint of symmetries. Such symmetries, expressed in Hamiltonian coordinates, or in standard polar coordinates, are recovered for McGehee-type coordinates of both collision-blow-up and infinity-blow-up kind. They form diffeomorphic commutative groups endowed with a Boolean structure. Expressed in Levi-Civita’s coordinates, the problem exhibits a larger group of symmetries, also commutative and presenting a Boolean structure.

Review: the atmospheric boundary layer

Science.gov (United States)

Garratt, J. R.

1994-10-01

An overview is given of the atmospheric boundary layer (ABL) over both continental and ocean surfaces, mainly from observational and modelling perspectives. Much is known about ABL structure over homogeneous land surfaces, but relatively little so far as the following are concerned, (i) the cloud-topped ABL (over the sea predominantly); (ii) the strongly nonhomogeneous and nonstationary ABL; (iii) the ABL over complex terrain. These three categories present exciting challenges so far as improved understanding of ABL behaviour and improved representation of the ABL in numerical models of the atmosphere are concerned.

Asymptotic stability of shear-flow solutions to incompressible viscous free boundary problems with and without surface tension

Science.gov (United States)

Tice, Ian

2018-04-01

This paper concerns the dynamics of a layer of incompressible viscous fluid lying above a rigid plane and with an upper boundary given by a free surface. The fluid is subject to a constant external force with a horizontal component, which arises in modeling the motion of such a fluid down an inclined plane, after a coordinate change. We consider the problem both with and without surface tension for horizontally periodic flows. This problem gives rise to shear-flow equilibrium solutions, and the main thrust of this paper is to study the asymptotic stability of the equilibria in certain parameter regimes. We prove that there exists a parameter regime in which sufficiently small perturbations of the equilibrium at time t=0 give rise to global-in-time solutions that return to equilibrium exponentially in the case with surface tension and almost exponentially in the case without surface tension. We also establish a vanishing surface tension limit, which connects the solutions with and without surface tension.

Numerical solutions of multi-dimensional solidification/melting problems by the dual reciprocity boundary element method

International Nuclear Information System (INIS)

Jo, Jong Chull; Shin, Won Ky

1997-01-01

This paper presents an effective and simple procedure for the simulation of the motion of the solid-liquid interfacial boundary and the transient temperature field during phase change process. To accomplish this purpose, an iterative implicit solution algorithm has been developed by employing the dual reciprocity boundary element method. The dual reciprocity boundary element approach provided in this paper is much simpler than the usual boundary element method applying a reciprocity principle and an available technique for dealing with domain integral of boundary element formulation simultaneously. The effectiveness of the present analysis method have been illustrated through comparisons of the calculation results of an example with its semi-analytical or other numerical solutions where available

Boundary Value Problems and Approximate Solutions

African Journals Online (AJOL)

Tadesse

Department of Mathematics, College of Natural and Computational Scineces, Mekelle ..... In this section, the Variational Iteration Method is applied to different forms of .... Some problems in non-Newtonian fluid mechanics, Ph.D. thesis, Wales.

Hydrogeological boundary settings in SR 97. Uncertainties in regional boundary settings and transfer of boundary conditions to site-scale models

International Nuclear Information System (INIS)

Follin, S.

1999-06-01

The SR 97 project presents a performance assessment (PA) of the overall safety of a hypothetical deep repository at three sites in Sweden arbitrarily named Aberg, Beberg and Ceberg. One component of this PA assesses the uncertainties in the hydrogeological modelling. This study focuses on uncertainties in boundary settings (size of model domain and boundary conditions) in the regional and site-scale hydrogeological modelling of the three sites used to simulating the possible transport of radionuclides from the emplacement waste packages through the host rock to the accessible environment. Model uncertainties associated with, for instance, parameter heterogeneity and structural interpretations are addressed in other studies. This study concludes that the regional modelling of the SR 97 project addresses uncertainties in the choice of boundary conditions and size of model domain differently at each site, although the overall handling is acceptable and in accordance with common modelling practice. For example, the treatment of uncertainties with regard to the ongoing post-glacial flushing of the Baltic Shield is creditably addressed although not exhaustive from a modelling point of view. A significant contribution of the performed modelling is the study of nested numerical models, i.e., the numerical interplay between regional and site-scale numerical models. In the site-scale modelling great efforts are made to address problems associated with (i) the telescopic mesh refinement (TMR) technique with regard to the stochastic continuum approach, and (ii) the transfer of boundary conditions between variable-density flow systems and flow systems that are constrained to treat uniform density flow. This study concludes that the efforts made to handle these problems are acceptable with regards to the objectives of the SR 97 project

Boundary layer flow past a circular cylinder in axial flow

International Nuclear Information System (INIS)

Sawchuk, S.P.; Zamir, M.; Camiletti, S.E.

1985-01-01

This paper discusses a study of the laminar boundary layer on a semi-infinite circular cylinder in axial incompressible flow. Unlike previous studies, the present study investigates a full range of this boundary layer problem to determine skin friction, heat transfer and other integral properties of the boundary layer

Heat conduction in a plate-type fuel element with time-dependent boundary conditions

International Nuclear Information System (INIS)

Faya, A.J.G.; Maiorino, J.R.

1981-01-01

A method for the solution of boundary-value problems with variable boundary conditions is applied to solve a heat conduction problem in a plate-type fuel element with time dependent film coefficient. The numerical results show the feasibility of the method in the solution of this class of problems. (Author) [pt

Second-order wave diffraction by a circular cylinder using scaled boundary finite element method

International Nuclear Information System (INIS)

Song, H; Tao, L

2010-01-01

The scaled boundary finite element method (SBFEM) has achieved remarkable success in structural mechanics and fluid mechanics, combing the advantage of both FEM and BEM. Most of the previous works focus on linear problems, in which superposition principle is applicable. However, many physical problems in the real world are nonlinear and are described by nonlinear equations, challenging the application of the existing SBFEM model. A popular idea to solve a nonlinear problem is decomposing the nonlinear equation to a number of linear equations, and then solves them individually. In this paper, second-order wave diffraction by a circular cylinder is solved by SBFEM. By splitting the forcing term into two parts, the physical problem is described as two second-order boundary-value problems with different asymptotic behaviour at infinity. Expressing the velocity potentials as a series of depth-eigenfunctions, both of the 3D boundary-value problems are decomposed to a number of 2D boundary-value sub-problems, which are solved semi-analytically by SBFEM. Only the cylinder boundary is discretised with 1D curved finite-elements on the circumference of the cylinder, while the radial differential equation is solved completely analytically. The method can be extended to solve more complex wave-structure interaction problems resulting in direct engineering applications.

Numerical solutions of multi-dimensional solidification/melting problems by the dual reciprocity boundary element method

Energy Technology Data Exchange (ETDEWEB)

Jo, Jong Chull; Shin, Won Ky [Korea Institute of Nuclear Safety, Taejon (Korea, Republic of)

1998-12-31

This paper presents an effective and simple procedure for the simulation of the motion of the solid-liquid interfacial boundary and the transient temperature field during phase change process. To accomplish this purpose, an iterative implicit solution algorithm has been developed by employing the dual reciprocity boundary element method. The dual reciprocity boundary element approach provided in this paper is much simpler than the usual boundary element method applying a reciprocity principle and an available technique for dealing with domain integral of boundary element formulation simultaneously. The effectiveness of the present analysis method have been illustrated through comparisons of the calculation results of an example with its semi-analytical or other numerical solutions where available. 22 refs., 3 figs. (Author)