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.

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.

Fourier analysis and boundary value problems

CERN Document Server

Gonzalez-Velasco, Enrique A

1996-01-01

Fourier Analysis and Boundary Value Problems provides a thorough examination of both the theory and applications of partial differential equations and the Fourier and Laplace methods for their solutions. Boundary value problems, including the heat and wave equations, are integrated throughout the book. Written from a historical perspective with extensive biographical coverage of pioneers in the field, the book emphasizes the important role played by partial differential equations in engineering and physics. In addition, the author demonstrates how efforts to deal with these problems have lead to wonderfully significant developments in mathematics.A clear and complete text with more than 500 exercises, Fourier Analysis and Boundary Value Problems is a good introduction and a valuable resource for those in the field.Key Features* Topics are covered from a historical perspective with biographical information on key contributors to the field* The text contains more than 500 exercises* Includes practical applicati...

A numerical solution of a singular boundary value problem arising in boundary layer theory.

Science.gov (United States)

Hu, Jiancheng

2016-01-01

In this paper, a second-order nonlinear singular boundary value problem is presented, which is equivalent to the well-known Falkner-Skan equation. And the one-dimensional third-order boundary value problem on interval [Formula: see text] is equivalently transformed into a second-order boundary value problem on finite interval [Formula: see text]. The finite difference method is utilized to solve the singular boundary value problem, in which the amount of computational effort is significantly less than the other numerical methods. The numerical solutions obtained by the finite difference method are in agreement with those obtained by previous authors.

Nonlinear Elliptic Boundary Value Problems at Resonance with Nonlinear Wentzell Boundary Conditions

Directory of Open Access Journals (Sweden)

Ciprian G. Gal

2017-01-01

Full Text Available Given a bounded domain Ω⊂RN with a Lipschitz boundary ∂Ω and p,q∈(1,+∞, we consider the quasilinear elliptic equation -Δpu+α1u=f in Ω complemented with the generalized Wentzell-Robin type boundary conditions of the form bx∇up-2∂nu-ρbxΔq,Γu+α2u=g on ∂Ω. In the first part of the article, we give necessary and sufficient conditions in terms of the given functions f, g and the nonlinearities α1, α2, for the solvability of the above nonlinear elliptic boundary value problems with the nonlinear boundary conditions. In other words, we establish a sort of “nonlinear Fredholm alternative” for our problem which extends the corresponding Landesman and Lazer result for elliptic problems with linear homogeneous boundary conditions. In the second part, we give some additional results on existence and uniqueness and we study the regularity of the weak solutions for these classes of nonlinear problems. More precisely, we show some global a priori estimates for these weak solutions in an L∞-setting.

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.

A two-dimensional embedded-boundary method for convection problems with moving boundaries

NARCIS (Netherlands)

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

2010-01-01

htmlabstractIn this work, a two-dimensional embedded-boundary algorithm for convection problems is presented. A moving body of arbitrary boundary shape is immersed in a Cartesian finite-volume grid, which is fixed in space. The boundary surface is reconstructed in such a way that only certain fluxes

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.

Boundary Value Problems Arising in Kalman Filtering

Directory of Open Access Journals (Sweden)

Sinem Ertürk

2009-01-01

Full Text Available The classic Kalman filtering equations for independent and correlated white noises are ordinary differential equations (deterministic or stochastic with the respective initial conditions. Changing the noise processes by taking them to be more realistic wide band noises or delayed white noises creates challenging partial differential equations with initial and boundary conditions. In this paper, we are aimed to give a survey of this connection between Kalman filtering and boundary value problems, bringing them into the attention of mathematicians as well as engineers dealing with Kalman filtering and boundary value problems.

Boundary Value Problems Arising in Kalman Filtering

Directory of Open Access Journals (Sweden)

Bashirov Agamirza

2008-01-01

Full Text Available The classic Kalman filtering equations for independent and correlated white noises are ordinary differential equations (deterministic or stochastic with the respective initial conditions. Changing the noise processes by taking them to be more realistic wide band noises or delayed white noises creates challenging partial differential equations with initial and boundary conditions. In this paper, we are aimed to give a survey of this connection between Kalman filtering and boundary value problems, bringing them into the attention of mathematicians as well as engineers dealing with Kalman filtering and boundary value problems.

Positive solutions for a fourth order boundary value problem

Directory of Open Access Journals (Sweden)

Bo Yang

2005-02-01

Full Text Available We consider a boundary value problem for the beam equation, in which the boundary conditions mean that the beam is embedded at one end and free at the other end. Some new estimates to the positive solutions to the boundary value problem are obtained. Some sufficient conditions for the existence of at least one positive solution for the boundary value problem are established. An example is given at the end of the paper to illustrate the main results.

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.

Modeling the summertime Arctic cloudy boundary layer

Energy Technology Data Exchange (ETDEWEB)

Curry, J.A.; Pinto, J.O. [Univ. of Colorado, Boulder, CO (United States); McInnes, K.L. [CSIRO Division of Atmospheric Research, Mordialloc (Australia)

1996-04-01

Global climate models have particular difficulty in simulating the low-level clouds during the Arctic summer. Model problems are exacerbated in the polar regions by the complicated vertical structure of the Arctic boundary layer. The presence of multiple cloud layers, a humidity inversion above cloud top, and vertical fluxes in the cloud that are decoupled from the surface fluxes, identified in Curry et al. (1988), suggest that models containing sophisticated physical parameterizations would be required to accurately model this region. Accurate modeling of the vertical structure of multiple cloud layers in climate models is important for determination of the surface radiative fluxes. This study focuses on the problem of modeling the layered structure of the Arctic summertime boundary-layer clouds and in particular, the representation of the more complex boundary layer type consisting of a stable foggy surface layer surmounted by a cloud-topped mixed layer. A hierarchical modeling/diagnosis approach is used. A case study from the summertime Arctic Stratus Experiment is examined. A high-resolution, one-dimensional model of turbulence and radiation is tested against the observations and is then used in sensitivity studies to infer the optimal conditions for maintaining two separate layers in the Arctic summertime boundary layer. A three-dimensional mesoscale atmospheric model is then used to simulate the interaction of this cloud deck with the large-scale atmospheric dynamics. An assessment of the improvements needed to the parameterizations of the boundary layer, cloud microphysics, and radiation in the 3-D model is made.

Uniqueness of solution to a stationary boundary kinetic problem

International Nuclear Information System (INIS)

Zhykharsky, A.V.

1992-01-01

The paper treats the question of uniqueness of solution to the boundary kinetic problem. This analysis is based on the accurate solutions to the stationary one-dimensional boundary kinetic problem for the limited plasma system. In the paper a simplified problem statement is used (no account is taken of the external magnetic field, a simplest form of boundary conditions is accepted) which, however, covers all features of the problem considered. Omitting the details of the conclusion we will write a set of Vlasov stationary kinetic equations for the cases of plane, cylindrical and spherical geometry of the problem. (author) 1 ref

The use of MACSYMA for solving elliptic boundary value problems

Science.gov (United States)

Thejll, Peter; Gilbert, Robert P.

1990-01-01

A boundary method is presented for the solution of elliptic boundary value problems. An approach based on the use of complete systems of solutions is emphasized. The discussion is limited to the Dirichlet problem, even though the present method can possibly be adapted to treat other boundary value problems.

On the solvability of initial boundary value problems for nonlinear ...

African Journals Online (AJOL)

In this paper, we study the initial boundary value problems for a non-linear time dependent Schrödinger equation with Dirichlet and Neumann boundary conditions, respectively. We prove the existence and uniqueness of solutions of the initial boundary value problems by using Galerkin's method. Keywords: Initial boundary ...

State space approach to mixed boundary value problems.

Science.gov (United States)

Chen, C. F.; Chen, M. M.

1973-01-01

A state-space procedure for the formulation and solution of mixed boundary value problems is established. This procedure is a natural extension of the method used in initial value problems; however, certain special theorems and rules must be developed. The scope of the applications of the approach includes beam, arch, and axisymmetric shell problems in structural analysis, boundary layer problems in fluid mechanics, and eigenvalue problems for deformable bodies. Many classical methods in these fields developed by Holzer, Prohl, Myklestad, Thomson, Love-Meissner, and others can be either simplified or unified under new light shed by the state-variable approach. A beam problem is included as an illustration.

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.

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...

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.

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 ...

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.)

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

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

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.

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...

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 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.

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.

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.

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.

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…

Bifurcation of solutions to Hamiltonian boundary value problems

Science.gov (United States)

McLachlan, R. I.; Offen, C.

2018-06-01

A bifurcation is a qualitative change in a family of solutions to an equation produced by varying parameters. In contrast to the local bifurcations of dynamical systems that are often related to a change in the number or stability of equilibria, bifurcations of boundary value problems are global in nature and may not be related to any obvious change in dynamical behaviour. Catastrophe theory is a well-developed framework which studies the bifurcations of critical points of functions. In this paper we study the bifurcations of solutions of boundary-value problems for symplectic maps, using the language of (finite-dimensional) singularity theory. We associate certain such problems with a geometric picture involving the intersection of Lagrangian submanifolds, and hence with the critical points of a suitable generating function. Within this framework, we then study the effect of three special cases: (i) some common boundary conditions, such as Dirichlet boundary conditions for second-order systems, restrict the possible types of bifurcations (for example, in generic planar systems only the A-series beginning with folds and cusps can occur); (ii) integrable systems, such as planar Hamiltonian systems, can exhibit a novel periodic pitchfork bifurcation; and (iii) systems with Hamiltonian symmetries or reversing symmetries can exhibit restricted bifurcations associated with the symmetry. This approach offers an alternative to the analysis of critical points in function spaces, typically used in the study of bifurcation of variational problems, and opens the way to the detection of more exotic bifurcations than the simple folds and cusps that are often found in examples.

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.

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.

The SMM Model as a Boundary Value Problem Using the Discrete Diffusion Equation

Science.gov (United States)

Campbell, Joel

2007-01-01

A generalized single step stepwise mutation model (SMM) is developed that takes into account an arbitrary initial state to a certain partial difference equation. This is solved in both the approximate continuum limit and the more exact discrete form. A time evolution model is developed for Y DNA or mtDNA that takes into account the reflective boundary modeling minimum microsatellite length and the original difference equation. A comparison is made between the more widely known continuum Gaussian model and a discrete model, which is based on modified Bessel functions of the first kind. A correction is made to the SMM model for the probability that two individuals are related that takes into account a reflecting boundary modeling minimum microsatellite length. This method is generalized to take into account the general n-step model and exact solutions are found. A new model is proposed for the step distribution.

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.

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

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

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.

An efficient realization of frequency dependent boundary conditions in an acoustic finite-difference time-domain model

DEFF Research Database (Denmark)

Escolano-Carrasco, José; Jacobsen, Finn; López, J.J.

2008-01-01

The finite-difference time-domain (FDTD) method provides a simple and accurate way of solving initial boundary value problems. However, most acoustic problems involve frequency dependent boundary conditions, and it is not easy to include such boundary conditions in an FDTD model. Although solutions...... to this problem exist, most of them have high computational costs, and stability cannot always be ensured. In this work, a solution is proposed based on "mixing modelling strategies"; this involves separating the FDTD mesh and the boundary conditions (a digital filter representation of the impedance...

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.

On a price formation free boundary model by Lasry and Lions

KAUST Repository

Caffarelli, Luis A.

2011-06-01

We discuss global existence and asymptotic behaviour of a price formation free boundary model introduced by Lasry and Lions in 2007. Our results are based on a construction which transforms the problem into the heat equation with specially prepared initial datum. The key point is that the free boundary present in the original problem becomes the zero level set of this solution. Using the properties of the heat operator we can show global existence, regularity and asymptotic results of the free boundary. 2011 Académie des sciences.

On a price formation free boundary model by Lasry and Lions

KAUST Repository

Caffarelli, Luis A.; Markowich, Peter A.; Pietschmann, Jan-F.

2011-01-01

We discuss global existence and asymptotic behaviour of a price formation free boundary model introduced by Lasry and Lions in 2007. Our results are based on a construction which transforms the problem into the heat equation with specially prepared initial datum. The key point is that the free boundary present in the original problem becomes the zero level set of this solution. Using the properties of the heat operator we can show global existence, regularity and asymptotic results of the free boundary. 2011 Académie des sciences.

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...

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

A minimally-resolved immersed boundary model for reaction-diffusion problems

OpenAIRE

Pal Singh Bhalla, A; Griffith, BE; Patankar, NA; Donev, A

2013-01-01

We develop an immersed boundary approach to modeling reaction-diffusion processes in dispersions of reactive spherical particles, from the diffusion-limited to the reaction-limited setting. We represent each reactive particle with a minimally-resolved "blob" using many fewer degrees of freedom per particle than standard discretization approaches. More complicated or more highly resolved particle shapes can be built out of a collection of reactive blobs. We demonstrate numerically that the blo...

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.

Micro worlds versus boundary objects in group model building; evidence from the literature on problem definition and model conceptulization

Energy Technology Data Exchange (ETDEWEB)

Zagonel, Aldo A. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States). Systems Engineering & Analysis; Andersen, David F. [University in Albany, NY (United States). The Rockefeller College of Public Affairs & Policy

2007-03-01

Based upon participant observation in group model building and content analysis of the system dynamics literature, we postulate that modeling efforts have a dual nature. On one hand, the modeling process aims to create a useful representation of a real-world system. This must be done, however, while aligning the clients’ mental models around a shared view of the system. There is significant overlap and confusion between these two goals and how they play out on a practical level. This research clarifies these distinctions by establishing an ideal-type dichotomy. To highlight the differences, we created two straw men: “micro world” characterizes a model that represents reality and “boundary object” represents a socially negotiated model. Using this framework, the literature was examined, revealing evidence for several competing views on problem definition and model conceptualization. The results are summarized in the text of this article, substantiated with strikingly polarized citations, often from the same authors. We also introduce hypotheses for the duality across the remaining phases of the modeling process. Finally, understanding and appreciation of the differences between these ideal types can promote constructive debate on their balance in system dynamics theory and practice.

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.

Coastal boundary layers in ocean modelling: an application to the Adriatic Sea

International Nuclear Information System (INIS)

Malanotte Rizzoli, P.; Dell'Orto, F.

1981-01-01

Boundary layers play an important role in modelling geophysical fluid-dynamical flows, in as much as they constitute regions of ageostrophic dynamics in which the physical balances characterizing the main interior of the water mass break down. A short synopsis is given of important boundary layers in ocean circulation modelling with specific emphasis drawn upon side wall boundary layers, namely those adjacent to the coastlines of the considered basin. Application of boundary layer analysis is thereafter made for one specific phenomenological situation, namely the Northern Adriatic Sea and the problem posed by its wintertime seasonal circulation. The analysis furnishes a mathematical model fo the coastal strip adjacent to the Italian shoreline, treated as a boundary layer in the density field, starting from general model equations valid throughout the interior of the Northern Adriatic. The boundary layer model is consequently used to modify the side wall boundary condition for the interior density field. Related numerical experiments are shown and compared with previous standard experiments in which the boundary layer contribution to the density field has not been considered. (author)

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.

Existence and uniqueness of solution for a model problem of transonic flow

International Nuclear Information System (INIS)

Tangmanee, S.

1985-11-01

A model problem of transonic flow ''the Tricomi equation'' bounded by the rectangular-curve boundary is studied. We transform the model problem into a symmetric positive system and an admissible boundary condition is posed. We show that with some conditions the existence and uniqueness of the solution are guaranteed. (author)

A moving boundary problem and orthogonal collocation in solving a dynamic liquid surfactant membrane model including osmosis and breakage

Directory of Open Access Journals (Sweden)

E.C. Biscaia Junior

2001-06-01

Full Text Available A dynamic kinetic-diffusive model for the extraction of metallic ions from aqueous liquors using liquid surfactant membranes is proposed. The model incorporates undesirable intrinsic phenomena such as swelling and breakage of the emulsion globules that have to be controlled during process operation. These phenomena change the spatial location of the chemical reaction during the course of extraction, resulting in a transient moving boundary problem. The orthogonal collocation method was used to transform the partial differential equations into an ordinary differential equation set that was solved by an implicit numerical routine. The model was found to be numerically stable and reliable in predicting the behaviour of zinc extraction with acidic extractant for long residence times.

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

Boundary regularity of Nevanlinna domains and univalent functions in model subspaces

International Nuclear Information System (INIS)

Baranov, Anton D; Fedorovskiy, Konstantin Yu

2011-01-01

In the paper we study boundary regularity of Nevanlinna domains, which have appeared in problems of uniform approximation by polyanalytic polynomials. A new method for constructing Nevanlinna domains with essentially irregular nonanalytic boundaries is suggested; this method is based on finding appropriate univalent functions in model subspaces, that is, in subspaces of the form K Θ =H 2 ominus ΘH 2 , where Θ is an inner function. To describe the irregularity of the boundaries of the domains obtained, recent results by Dolzhenko about boundary regularity of conformal mappings are used. Bibliography: 18 titles.

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...

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.

Oblique radiation lateral open boundary conditions for a regional climate atmospheric model

Science.gov (United States)

Cabos Narvaez, William; De Frutos Redondo, Jose Antonio; Perez Sanz, Juan Ignacio; Sein, Dmitry

2013-04-01

The prescription of lateral boundary conditions in regional atmospheric models represent a very important issue for limited area models. The ill-posed nature of the open boundary conditions makes it necessary to devise schemes in order to filter spurious wave reflections at boundaries, being desirable to have one boundary condition per variable. On the other side, due to the essentially hyperbolic nature of the equations solved in state of the art atmospheric models, external data is required only for inward boundary fluxes. These circumstances make radiation lateral boundary conditions a good choice for the filtering of spurious wave reflections. Here we apply the adaptive oblique radiation modification proposed by Mikoyada and Roseti to each of the prognostic variables of the REMO regional atmospheric model and compare it to the more common normal radiation condition used in REMO. In the proposed scheme, special attention is paid to the estimation of the radiation phase speed, essential to detecting the direction of boundary fluxes. One of the differences with the classical scheme is that in case of outward propagation, the adaptive nudging imposed in the boundaries allows to minimize under and over specifications problems, adequately incorporating the external information.

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.

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.

Modeling of uncertainties in statistical inverse problems

International Nuclear Information System (INIS)

Kaipio, Jari

2008-01-01

In all real world problems, the models that tie the measurements to the unknowns of interest, are at best only approximations for reality. While moderate modeling and approximation errors can be tolerated with stable problems, inverse problems are a notorious exception. Typical modeling errors include inaccurate geometry, unknown boundary and initial data, properties of noise and other disturbances, and simply the numerical approximations of the physical models. In principle, the Bayesian approach to inverse problems, in which all uncertainties are modeled as random variables, is capable of handling these uncertainties. Depending on the type of uncertainties, however, different strategies may be adopted. In this paper we give an overview of typical modeling errors and related strategies within the Bayesian framework.

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.

On mathematical modelling and numerical simulation of transient compressible flow across open boundaries

Energy Technology Data Exchange (ETDEWEB)

Rian, Kjell Erik

2003-07-01

In numerical simulations of turbulent reacting compressible flows, artificial boundaries are needed to obtain a finite computational domain when an unbounded physical domain is given. Artificial boundaries which fluids are free to cross are called open boundaries. When calculating such flows, non-physical reflections at the open boundaries may occur. These reflections can pollute the solution severely, leading to inaccurate results, and the generation of spurious fluctuations may even cause the numerical simulation to diverge. Thus, a proper treatment of the open boundaries in numerical simulations of turbulent reacting compressible flows is required to obtain a reliable solution for realistic conditions. A local quasi-one-dimensional characteristic-based open-boundary treatment for the Favre-averaged governing equations for time-dependent three-dimensional multi-component turbulent reacting compressible flow is presented. A k-{epsilon} model for turbulent compressible flow and Magnussen's EDC model for turbulent combustion is included in the analysis. The notion of physical boundary conditions is incorporated in the method, and the conservation equations themselves are applied on the boundaries to complement the set of physical boundary conditions. A two-dimensional finite-difference-based computational fluid dynamics code featuring high-order accurate numerical schemes was developed for the numerical simulations. Transient numerical simulations of the well-known, one-dimensional shock-tube problem, a two-dimensional pressure-tower problem in a decaying turbulence field, and a two-dimensional turbulent reacting compressible flow problem have been performed. Flow- and combustion-generated pressure waves seem to be well treated by the non-reflecting subsonic open-boundary conditions. Limitations of the present open-boundary treatment are demonstrated and discussed. The simple and solid physical basis of the method makes it both favourable and relatively easy to

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.

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.

Exact solution of the one-dimensional Hubbard model with arbitrary boundary magnetic fields

Energy Technology Data Exchange (ETDEWEB)

Li, Yuan-Yuan; Cao, Junpeng [Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190 (China); Yang, Wen-Li [Institute of Modern Physics, Northwest University, Xian 710069 (China); Beijing Center for Mathematics and Information Interdisciplinary Sciences, Beijing, 100048 (China); Shi, Kangjie [Institute of Modern Physics, Northwest University, Xian 710069 (China); Wang, Yupeng, E-mail: yupeng@iphy.ac.cn [Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190 (China)

2014-02-15

The one-dimensional Hubbard model with arbitrary boundary magnetic fields is solved exactly via the Bethe ansatz methods. With the coordinate Bethe ansatz in the charge sector, the second eigenvalue problem associated with the spin sector is constructed. It is shown that the second eigenvalue problem can be transformed into that of the inhomogeneous XXX spin chain with arbitrary boundary fields which can be solved via the off-diagonal Bethe ansatz method.

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.

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

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....

Model-based estimation with boundary side information or boundary regularization

International Nuclear Information System (INIS)

Chiao, P.C.; Rogers, W.L.; Fessler, J.A.; Clinthorne, N.H.; Hero, A.O.

1994-01-01

The authors have previously developed a model-based strategy for joint estimation of myocardial perfusion and boundaries using ECT (Emission Computed Tomography). The authors have also reported difficulties with boundary estimation in low contrast and low count rate situations. In this paper, the authors propose using boundary side information (obtainable from high resolution MRI and CT images) or boundary regularization to improve both perfusion and boundary estimation in these situations. To fuse boundary side information into the emission measurements, the authors formulate a joint log-likelihood function to include auxiliary boundary measurements as well as ECT projection measurements. In addition, the authors introduce registration parameters to align auxiliary boundary measurements with ECT measurements and jointly estimate these parameters with other parameters of interest from the composite measurements. In simulated PET O-15 water myocardial perfusion studies using a simplified model, the authors show that the joint estimation improves perfusion estimation performance and gives boundary alignment accuracy of <0.5 mm even at 0.2 million counts. The authors implement boundary regularization through formulating a penalized log-likelihood function. The authors also demonstrate in simulations that simultaneous regularization of the epicardial boundary and myocardial thickness gives comparable perfusion estimation accuracy with the use of boundary side information

Oracle estimation of parametric models under boundary constraints.

Science.gov (United States)

Wong, Kin Yau; Goldberg, Yair; Fine, Jason P

2016-12-01

In many classical estimation problems, the parameter space has a boundary. In most cases, the standard asymptotic properties of the estimator do not hold when some of the underlying true parameters lie on the boundary. However, without knowledge of the true parameter values, confidence intervals constructed assuming that the parameters lie in the interior are generally over-conservative. A penalized estimation method is proposed in this article to address this issue. An adaptive lasso procedure is employed to shrink the parameters to the boundary, yielding oracle inference which adapt to whether or not the true parameters are on the boundary. When the true parameters are on the boundary, the inference is equivalent to that which would be achieved with a priori knowledge of the boundary, while if the converse is true, the inference is equivalent to that which is obtained in the interior of the parameter space. The method is demonstrated under two practical scenarios, namely the frailty survival model and linear regression with order-restricted parameters. Simulation studies and real data analyses show that the method performs well with realistic sample sizes and exhibits certain advantages over standard methods. © 2016, The International Biometric Society.

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...

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

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.

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.

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)

A Numerical Model of Anisotropic Mass Transport Through Grain Boundary Networks

Science.gov (United States)

Wang, Yibo

Tin (Sn) thin films are commonly used in electronic circuit applications as coatings on contacts and solders for joining components. It is widely observed, for some such system, that whiskers---long, thin crystalline structures---emerge and grow from the film. The Sn whisker phenomenon has become a highly active research area since Sn whiskers have caused a large amount of damage and loss in manufacturing, military, medical and power industries. Though lead (Pb) addition to Sn has been used to solve this problem for over five decades, the adverse environmental and health effects of Pb have motivated legislation to severely constrain Pb use in society. People are researching and seeking the reasons which cause whiskers and corresponding methods to solve the problem. The contributing factors to cause a Sn whisker are potentially many and much still remains unknown. Better understanding of fundamental driving forces should point toward strategies to improve (a) the accuracy with which we can predict whisker formation, and (b) our ability to mitigate the phenomenon. This thesis summarizes recent important research achievements in understanding Sn whisker formation and growth, both experimentally and theoretically. Focus is then placed on examining the role that anisotropy in grain boundary diffusivity plays in determining whisker characteristics (specifically, whether they form and, if so, where on a surface). To study this aspect of the problem and to enable future studies on stress driven grain boundary diffusion, this thesis presents a numerical anisotropic mass transport model. In addition to presenting details of the model and implementation, model predictions for a set of increasingly complex grain boundary networks are discussed. Preliminary results from the model provide evidence that anisotropic grain boundary diffusion may be a primary driving mechanism in whisker formation.

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.

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

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.

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)

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...

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

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

Conformal boundary loop models

International Nuclear Information System (INIS)

Jacobsen, Jesper Lykke; Saleur, Hubert

2008-01-01

We study a model of densely packed self-avoiding loops on the annulus, related to the Temperley-Lieb algebra with an extra idempotent boundary generator. Four different weights are given to the loops, depending on their homotopy class and whether they touch the outer rim of the annulus. When the weight of a contractible bulk loop x≡q+q -1 element of (-2,2], this model is conformally invariant for any real weight of the remaining three parameters. We classify the conformal boundary conditions and give exact expressions for the corresponding boundary scaling dimensions. The amplitudes with which the sectors with any prescribed number and types of non-contractible loops appear in the full partition function Z are computed rigorously. Based on this, we write a number of identities involving Z which hold true for any finite size. When the weight of a contractible boundary loop y takes certain discrete values, y r ≡([r+1] q )/([r] q ) with r integer, other identities involving the standard characters K r,s of the Virasoro algebra are established. The connection with Dirichlet and Neumann boundary conditions in the O(n) model is discussed in detail, and new scaling dimensions are derived. When q is a root of unity and y=y r , exact connections with the A m type RSOS model are made. These involve precise relations between the spectra of the loop and RSOS model transfer matrices, valid in finite size. Finally, the results where y=y r are related to the theory of Temperley-Lieb cabling

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.

Boundary correlators in supergroup WZNW models

Energy Technology Data Exchange (ETDEWEB)

Creutzig, T.; Schomerus, V.

2008-04-15

We investigate correlation functions for maximally symmetric boundary conditions in the WZNW model on GL(11). Special attention is payed to volume filling branes. Generalizing earlier ideas for the bulk sector, we set up a Kac-Wakimotolike formalism for the boundary model. This first order formalism is then used to calculate bulk-boundary 2-point functions and the boundary 3-point functions of the model. The note ends with a few comments on correlation functions of atypical fields, point-like branes and generalizations to other supergroups. (orig.)

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

A Source-Term Based Boundary Layer Bleed/Effusion Model for Passive Shock Control

Science.gov (United States)

Baurle, Robert A.; Norris, Andrew T.

2011-01-01

A modeling framework for boundary layer effusion has been developed based on the use of source (or sink) terms instead of the usual practice of specifying bleed directly as a boundary condition. This framework allows the surface boundary condition (i.e. isothermal wall, adiabatic wall, slip wall, etc.) to remain unaltered in the presence of bleed. This approach also lends itself to easily permit the addition of empirical models for second order effects that are not easily accounted for by simply defining effective transpiration values. Two effusion models formulated for supersonic flows have been implemented into this framework; the Doerffer/Bohning law and the Slater formulation. These models were applied to unit problems that contain key aspects of the flow physics applicable to bleed systems designed for hypersonic air-breathing propulsion systems. The ability of each model to predict bulk bleed properties was assessed, as well as the response of the boundary layer as it passes through and downstream of a porous bleed system. The model assessment was performed with and without the presence of shock waves. Three-dimensional CFD simulations that included the geometric details of the porous plate bleed systems were also carried out to supplement the experimental data, and provide additional insights into the bleed flow physics. Overall, both bleed formulations fared well for the tests performed in this study. However, the sample of test problems considered in this effort was not large enough to permit a comprehensive validation of the models.

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.)

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

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

Integro-differential equations of fractional order with nonlocal fractional boundary conditions associated with financial asset model

Directory of Open Access Journals (Sweden)

Bashir Ahmad

2013-02-01

Full Text Available In this article, we discuss the existence of solutions for a boundary-value problem of integro-differential equations of fractional order with nonlocal fractional boundary conditions by means of some standard tools of fixed point theory. Our problem describes a more general form of fractional stochastic dynamic model for financial asset. An illustrative example is also presented.

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/

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)

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

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

New boundary conditions for 3D RF modelling

International Nuclear Information System (INIS)

Ko, K.; Nelson, E.; Fitze, H.

1990-01-01

The new capabilities are being implemented into the 3D particle-in-cell code, ARGUS, which will reduce substantially both problem size and computing time when modeling realistic geometries with high accuracies. In the time domain, a cylindrical radiative boundary condition will enable traveling wave propagation to be simulated in accelerator structures. An application of interest is the input coupler in the SLAC x-band high-gradient structure where local field gradients and impedance matching are important issues. In the frequency domain, a quasi-periodic boundary condition will facilitate the cold-test analysis of 3D periodic structures where many calculations are required to generate an ω β diagram. Present applications include the crossed-field amplifier cavity and the cluster klystron cavity

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.

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

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.

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

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

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.

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.

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

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.

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.

Model-based estimation with boundary side information or boundary regularization [cardiac emission CT].

Science.gov (United States)

Chiao, P C; Rogers, W L; Fessler, J A; Clinthorne, N H; Hero, A O

1994-01-01

The authors have previously developed a model-based strategy for joint estimation of myocardial perfusion and boundaries using ECT (emission computed tomography). They have also reported difficulties with boundary estimation in low contrast and low count rate situations. Here they propose using boundary side information (obtainable from high resolution MRI and CT images) or boundary regularization to improve both perfusion and boundary estimation in these situations. To fuse boundary side information into the emission measurements, the authors formulate a joint log-likelihood function to include auxiliary boundary measurements as well as ECT projection measurements. In addition, they introduce registration parameters to align auxiliary boundary measurements with ECT measurements and jointly estimate these parameters with other parameters of interest from the composite measurements. In simulated PET O-15 water myocardial perfusion studies using a simplified model, the authors show that the joint estimation improves perfusion estimation performance and gives boundary alignment accuracy of <0.5 mm even at 0.2 million counts. They implement boundary regularization through formulating a penalized log-likelihood function. They also demonstrate in simulations that simultaneous regularization of the epicardial boundary and myocardial thickness gives comparable perfusion estimation accuracy with the use of boundary side information.

Explicit treatment for Dirichlet, Neumann and Cauchy boundary conditions in POD-based reduction of groundwater models

Science.gov (United States)

Gosses, Moritz; Nowak, Wolfgang; Wöhling, Thomas

2018-05-01

In recent years, proper orthogonal decomposition (POD) has become a popular model reduction method in the field of groundwater modeling. It is used to mitigate the problem of long run times that are often associated with physically-based modeling of natural systems, especially for parameter estimation and uncertainty analysis. POD-based techniques reproduce groundwater head fields sufficiently accurate for a variety of applications. However, no study has investigated how POD techniques affect the accuracy of different boundary conditions found in groundwater models. We show that the current treatment of boundary conditions in POD causes inaccuracies for these boundaries in the reduced models. We provide an improved method that splits the POD projection space into a subspace orthogonal to the boundary conditions and a separate subspace that enforces the boundary conditions. To test the method for Dirichlet, Neumann and Cauchy boundary conditions, four simple transient 1D-groundwater models, as well as a more complex 3D model, are set up and reduced both by standard POD and POD with the new extension. We show that, in contrast to standard POD, the new method satisfies both Dirichlet and Neumann boundary conditions. It can also be applied to Cauchy boundaries, where the flux error of standard POD is reduced by its head-independent contribution. The extension essentially shifts the focus of the projection towards the boundary conditions. Therefore, we see a slight trade-off between errors at model boundaries and overall accuracy of the reduced model. The proposed POD extension is recommended where exact treatment of boundary conditions is required.

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.

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

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.

Three dimensional grain boundary modeling in polycrystalline plasticity

Science.gov (United States)

Yalçinkaya, Tuncay; Özdemir, Izzet; Fırat, Ali Osman

2018-05-01

At grain scale, polycrystalline materials develop heterogeneous plastic deformation fields, localizations and stress concentrations due to variation of grain orientations, geometries and defects. Development of inter-granular stresses due to misorientation are crucial for a range of grain boundary (GB) related failure mechanisms, such as stress corrosion cracking (SCC) and fatigue cracking. Local crystal plasticity finite element modelling of polycrystalline metals at micron scale results in stress jumps at the grain boundaries. Moreover, the concepts such as the transmission of dislocations between grains and strength of the grain boundaries are not included in the modelling. The higher order strain gradient crystal plasticity modelling approaches offer the possibility of defining grain boundary conditions. However, these conditions are mostly not dependent on misorientation of grains and can define only extreme cases. For a proper definition of grain boundary behavior in plasticity, a model for grain boundary behavior should be incorporated into the plasticity framework. In this context, a particular grain boundary model ([l]) is incorporated into a strain gradient crystal plasticity framework ([2]). In a 3-D setting, both bulk and grain boundary models are implemented as user-defined elements in Abaqus. The strain gradient crystal plasticity model works in the bulk elements and considers displacements and plastic slips as degree of freedoms. Interface elements model the plastic slip behavior, yet they do not possess any kind of mechanical cohesive behavior. The physical aspects of grain boundaries and the performance of the model are addressed through numerical examples.

Mathematical problems in modeling artificial heart

Directory of Open Access Journals (Sweden)

Ahmed N. U.

1995-01-01

Full Text Available In this paper we discuss some problems arising in mathematical modeling of artificial hearts. The hydrodynamics of blood flow in an artificial heart chamber is governed by the Navier-Stokes equation, coupled with an equation of hyperbolic type subject to moving boundary conditions. The flow is induced by the motion of a diaphragm (membrane inside the heart chamber attached to a part of the boundary and driven by a compressor (pusher plate. On one side of the diaphragm is the blood and on the other side is the compressor fluid. For a complete mathematical model it is necessary to write the equation of motion of the diaphragm and all the dynamic couplings that exist between its position, velocity and the blood flow in the heart chamber. This gives rise to a system of coupled nonlinear partial differential equations; the Navier-Stokes equation being of parabolic type and the equation for the membrane being of hyperbolic type. The system is completed by introducing all the necessary static and dynamic boundary conditions. The ultimate objective is to control the flow pattern so as to minimize hemolysis (damage to red blood cells by optimal choice of geometry, and by optimal control of the membrane for a given geometry. The other clinical problems, such as compatibility of the material used in the construction of the heart chamber, and the membrane, are not considered in this paper. Also the dynamics of the valve is not considered here, though it is also an important element in the overall design of an artificial heart. We hope to model the valve dynamics in later paper.

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.

A priori bounds for solutions of two-point boundary value problems using differential inequalities

International Nuclear Information System (INIS)

Vidossich, G.

1979-01-01

Two point boundary value problems for systems of differential equations are studied with a new approach based on differential inequalities of first order. This leads to the following results: (i) one-sided conditions are enough, in the sense that the inner product is substituted to the norm; (ii) the upper bound exists for practically any kind of equations and boundary value problem if the interval is sufficiently small since it depends on the Peano existence theorem; (iii) the bound seems convenient when the equation has some singularity in t as well as when sigular problems are considered. (author)

Modifications of the κ-ε model in GOTHIC near physical boundaries

International Nuclear Information System (INIS)

Analytis, G.Th.; Andreani, M.

2000-01-01

In CFD-type codes like the containment analysis code GOTHIC, one of the options that can be used for modelling of turbulence is the k - ε model. Though, in contrast to other CFD codes which are tailored for performing detailed CFD calculations with a large number of spatial meshes, in codes like GOTHIC which are primarily aiming at calculating transients in reactor containments, one generally uses coarse meshes. The solution of the two parabolic k - ε model equations requires the definition of boundary conditions at physical boundaries and this, in rum, requires very small spatial meshes near these boundaries. Hence, while in codes like CFX this is properly done, in codes like GOTHIC, this is done in an indirect and non-rigorous fashion, exactly due to the fact that the spatial meshes are usually large; this can have catastrophic consequences during the calculation of a transient and in this work, we shall give some examples of this and outline a method by which this problem can be by-passed. (author)

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.

DYNAMIC SURFACE BOUNDARY-CONDITIONS - A SIMPLE BOUNDARY MODEL FOR MOLECULAR-DYNAMICS SIMULATIONS

NARCIS (Netherlands)

JUFFER, AH; BERENDSEN, HJC

1993-01-01

A simple model for the treatment of boundaries in molecular dynamics simulations is presented. The method involves the positioning of boundary atoms on a surface that surrounds a system of interest. The boundary atoms interact with the inner region and represent the effect of atoms outside the

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

discretization based on dual boundary elements residing on a suitably defined dual mesh. We devote also a particular attention to implementation-oriented details of our new technique that will allow the rapid incorporation of our finding in one's own EEG forward solution technology. We conclude by showing how the resulting forward EEG problems show favorable properties with respect to previously proposed schemes, and we show their applicability to real-case modeling scenarios obtained from Magnetic Resonance Imaging (MRI) data. xml:lang="fr" Malheureusement, il est également reconnu dans la littérature que leur précision diminue particulièrement lorsque la source de courant est dipolaire et se situe près de la surface du cerveau. Ce défaut constitue une importante limitation, étant donné qu'au cours d'une session d'imagerie EEG à haute résolution, plusieurs solutions du problème direct EEG sont requises, pour lesquelles les sources de courant sont proches ou sur la surface de cerveau. Ce travail présente d'abord une analyse des opérateurs intervenant dans le problème direct et leur discrétisation. Nous montrons que plusieurs discrétisations couramment utilisées ne conviennent que dans un cadre L2, nécessitant que le terme d'expansion soit une fonction de carré intégrable. Dès lors, ces techniques ne sont pas cohérentes avec les propriétés spectrales des opérateurs en termes d'espaces de Sobolev d'ordre fractionnaire. Nous développons ensuite une nouvelle stratégie de discrétisation conforme aux espaces de Sobolev avec des fonctions d'expansion moins régulières, donnant lieu à une nouvelle formulation intégrale. Le solveur résultant présente des propriétés favorables par rapport aux méthodes existantes et réduit sensiblement le recours à un maillage adaptatif et autres stratégies actuellement requises pour améliorer la précision du calcul. Les résultats numériques présentés corroborent les développements théoriques et mettent en �

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

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

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.

On one model problem for the reaction-diffusion-advection equation

Science.gov (United States)

Davydova, M. A.; Zakharova, S. A.; Levashova, N. T.

2017-09-01

The asymptotic behavior of the solution with boundary layers in the time-independent mathematical model of reaction-diffusion-advection arising when describing the distribution of greenhouse gases in the surface atmospheric layer is studied. On the basis of the asymptotic method of differential inequalities, the existence of a boundary-layer solution and its asymptotic Lyapunov stability as a steady-state solution of the corresponding parabolic problem is proven. One of the results of this work is the determination of the local domain of the attraction of a boundary-layer solution.

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.

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.

Research on a Small Signal Stability Region Boundary Model of the Interconnected Power System with Large-Scale Wind Power

Directory of Open Access Journals (Sweden)

Wenying Liu

2015-03-01

Full Text Available For the interconnected power system with large-scale wind power, the problem of the small signal stability has become the bottleneck of restricting the sending-out of wind power as well as the security and stability of the whole power system. Around this issue, this paper establishes a small signal stability region boundary model of the interconnected power system with large-scale wind power based on catastrophe theory, providing a new method for analyzing the small signal stability. Firstly, we analyzed the typical characteristics and the mathematic model of the interconnected power system with wind power and pointed out that conventional methods can’t directly identify the topological properties of small signal stability region boundaries. For this problem, adopting catastrophe theory, we established a small signal stability region boundary model of the interconnected power system with large-scale wind power in two-dimensional power injection space and extended it to multiple dimensions to obtain the boundary model in multidimensional power injection space. Thirdly, we analyzed qualitatively the topological property’s changes of the small signal stability region boundary caused by large-scale wind power integration. Finally, we built simulation models by DIgSILENT/PowerFactory software and the final simulation results verified the correctness and effectiveness of the proposed model.

Tokamak plasma boundary layer model

International Nuclear Information System (INIS)

Volkov, T.F.; Kirillov, V.D.

1983-01-01

A model has been developed for the limiter layer and for the boundary region of the plasma column in a tokamak to facilitate analytic calculations of the thickness of the limiter layers, the profiles and boundary values of the temperature and the density under various conditions, and the difference between the electron and ion temperatures. This model can also be used to analyze the recycling of neutrals, the energy and particle losses to the wall and the limiter, and other characteristics

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.

Tricritical Ising model with a boundary

International Nuclear Information System (INIS)

De Martino, A.; Moriconi, M.

1998-03-01

We study the integrable and supersymmetric massive φ (1,3) deformation of the tricritical Ising model in the presence of a boundary. We use constraints from supersymmetry in order to compute the exact boundary S-matrices, which turn out to depend explicitly on the topological charge of the supersymmetry algebra. We also solve the general boundary Yang-Baxter equation and show that in appropriate limits the general reflection matrices go over the supersymmetry preserving solutions. Finally, we briefly discuss the possible connection between our reflection matrices and boundary perturbations within the framework of perturbed boundary conformal field theory. (author)

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

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...

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....

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.

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

Slarti: A boundary condition editor for a coupled climate model

Science.gov (United States)

Mickelson, S. A.; Jacob, R. L.; Pierrehumbert, R.

2006-12-01

One of the largest barriers to making climate models more flexible is the difficulty in creating new boundary conditions, especially for "deep time" paleoclimate cases where continents are in different positions. Climate models consist of several mutually-interacting component models and the boundary conditions must be consistent between them. We have developed a program called Slarti which uses a Graphical User Interface and a set of consistency rules to aid researchers in creating new, consistent, boundary condition files for the Fast Ocean Atmosphere Model (FOAM). Users can start from existing mask, topography, or bathymetry data or can build a "world" entirely from scratch (e.g. a single island continent). Once a case has been started, users can modify mask, vegetation, bathymetry, topography, and river flow fields by drawing new data through a "paint" interface. Users activate a synchronization button which goes through the fields to eliminate inconsistencies. When the changes are complete and save is selected, Slarti creates all the necessary files for an initial run of FOAM. The data is edited at the highest resolution (the ocean-land surface in FOAM) and then interpolated to the atmosphere resolution. Slarti was implemented in Java to maintain portability across platforms. We also relied heavily on Java Swing components to create the interface. This allowed us to create an object-oriented interface that could be used on many different systems. Since Slarti allows users to visualize their changes, they are able to see areas that may cause problems when the model is ran. Some examples would be lakes from the river flow field and narrow trenches within the bathymetry. Through different checks and options available through its interface, Slarti makes the process of creating new boundary conditions for FOAM easier and faster while reducing the chance for user errors.

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.

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)

Blurring Boundaries: From the Danish Welfare State to the European Social Model?

DEFF Research Database (Denmark)

Neergaard, Ulla; Nielsen, Ruth

Abstract: This paper builds on the results obtained in the so-called Blurring Boundaries project which was undertaken at the Law Department, Copenhagen Business School, in the period from 2007 to 2009. It looks at the sustainability of the Danish welfare state in an EU law context and on the inte......Abstract: This paper builds on the results obtained in the so-called Blurring Boundaries project which was undertaken at the Law Department, Copenhagen Business School, in the period from 2007 to 2009. It looks at the sustainability of the Danish welfare state in an EU law context...... and on the integration 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...

A high-latitude, low-latitude boundary layer model of the convection current system

International Nuclear Information System (INIS)

Siscoe, G.L.; Lotko, W.; Sonnerup, B.U.O.

1991-01-01

Observations suggest that both the high- and low-latitude boundary layers contribute to magnetospheric convection, and that their contributions are linked. In the interpretation pursued here, the high-latitude boundary layer (HBL) generates the voltage while the low-latitude boundary layer (LBL) generates the current for the part of the convection electric circuit that closes through the ionosphere. This paper gives a model that joins the high- and low-latitude boundary layers consistently with the ionospheric Ohm's law. It describes an electric circuit linking both boundary layers, the region 1 Birkeland currents, and the ionospheric Pedersen closure currents. The model works by using the convection electric field that the ionosphere receives from the HBL to determine two boundary conditions to the equations that govern viscous LBL-ionosphere coupling. The result provides the needed self-consistent coupling between the two boundary layers and fully specifies the solution for the viscous LBL-ionosphere coupling equations. The solution shows that in providing the current required by the ionospheric Ohm's law, the LBL needs only a tenth of the voltage that spans the HBL. The solution also gives the latitude profiles of the ionospheric electric field, parallel currents, and parallel potential. It predicts that the plasma in the inner part of the LBL moves sunward instead of antisunward and that, as the transpolar potential decreases below about 40 kV, reverse polarity (region 0) currents appear at the poleward border of the region 1 currents. A possible problem with the model is its prediction of a thin boundary layer (∼1000 km), whereas thicknesses inferred from satellite data tend to be greater

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}$.

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.

Multi-scale model analysis of boundary layer ozone over East Asia

Directory of Open Access Journals (Sweden)

M. Lin

2009-05-01

Full Text Available This study employs the regional Community Multiscale Air Quality (CMAQ model to examine seasonal and diurnal variations of boundary layer ozone (O₃ over East Asia. We evaluate the response of model simulations of boundary layer O₃ to the choice of chemical mechanisms, meteorological fields, boundary conditions, and model resolutions. Data obtained from surface stations, aircraft measurements, and satellites are used to advance understanding of O₃ chemistry and mechanisms over East Asia and evaluate how well the model represents the observed features. Satellite measurements and model simulations of summertime rainfall are used to assess the impact of the Asian monsoon on O₃ production. Our results suggest that summertime O₃ over Central Eastern China is highly sensitive to cloud cover and monsoonal rainfall over this region. Thus, accurate simulation of the East Asia summer monsoon is critical to model analysis of atmospheric chemistry over China. Examination of hourly summertime O₃ mixing ratios from sites in Japan confirms the important role of diurnal boundary layer fluctuations in controlling ground-level O₃. By comparing five different model configurations with observations at six sites, the specific mechanisms responsible for model behavior are identified and discussed. In particular, vertical mixing, urban chemistry, and dry deposition depending on boundary layer height strongly affect model ability to capture observed behavior. Central Eastern China appears to be the most sensitive region in our study to the choice of chemical mechanisms. Evaluation with TRACE-P aircraft measurements reveals that neither the CB4 nor the SAPRC99 mechanisms consistently capture observed behavior of key photochemical oxidants in springtime. However, our analysis finds that SAPRC99 performs somewhat better in simulating mixing ratios of H₂O₂ (hydrogen peroxide

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.

Boundary representation modelling techniques

CERN Document Server

2006-01-01

Provides the most complete presentation of boundary representation solid modelling yet publishedOffers basic reference information for software developers, application developers and users Includes a historical perspective as well as giving a background for modern research.

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.

International Conference on Computer Modelling of Seas and Coastal Regions and Boundary Elements and Fluid Dynamics

CERN Document Server

Partridge, P; Boundary Elements in Fluid Dynamics

1992-01-01

This book Boundary Elements in Fluid Dynamics is the second volume of the two volume proceedings of the International Conference on Computer Modelling of Seas and Coastal Regions and Boundary Elements and Fluid Dynamics, held in Southampton, U.K., in April 1992. The Boundary Element Method (BEM) is now fully established as an ac­ curate and successful technique for solving engineering problems in a wide range of fields. The success of the method is due to its advantages in data reduction, as only the boundary of the region is modelled. Thus moving boundaries may be more easily handled, which is not the case if domain methods are used. In addition, the method is easily able to model regions to extending to infinity. Fluid mechanics is traditionally one of the most challenging areas of engi­ neering, the simulation of fluid motion, particularly in three dimensions, is always a serious test for any numerical method, and is an area in which BEM analysis may be used taking full advantage of its special character...

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

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)

Boundary layer heights and surface fluxes of momentum and heat derived from ECMWF data for use in pollutant dispersion models - problems with data accuracy

Energy Technology Data Exchange (ETDEWEB)

Wotawa, G. [Univ. of Agricultural Sciences, Inst. of Meteorology and Physics, Vienna (Austria); Stohl, A. [Ludwig-Maximilians-Univ. Muenchen, Munich (Germany)

1997-10-01

Certain boundary layer parameters, especially boundary layer heights, are very important for pollutant dispersion modelling. On the regional scale (>- 100 km), data of the numerical weather prediction model of the European Centre for Medium-Range Weather Forecasts are often used for that purpose. Based on ECMWF data, the meteorological preprocessor FLEXTRA for Lagrangian air quality simulation models and the Lagrangian particle diffusion model FLEXPART have been developed. Using analyses and short term forecasts, a temporal resolution of three hours can be achieved. Some alternative methods to obtain boundary layer parameters can be applied, producing different results which affect all subsequent calculations, for instance the calculation of boundary layer trajectories and the dispersion of air pollutants. (au)

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

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.

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

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

Free-boundary models of a meltwater conduit

KAUST Repository

Dallaston, Michael C.; Hewitt, Ian J.

2014-01-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

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

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

Diffuse-interface polycrystal plasticity: expressing grain boundaries as geometrically necessary dislocations

Science.gov (United States)

Admal, Nikhil Chandra; Po, Giacomo; Marian, Jaime

2017-12-01

The standard way of modeling plasticity in polycrystals is by using the crystal plasticity model for single crystals in each grain, and imposing suitable traction and slip boundary conditions across grain boundaries. In this fashion, the system is modeled as a collection of boundary-value problems with matching boundary conditions. In this paper, we develop a diffuse-interface crystal plasticity model for polycrystalline materials that results in a single boundary-value problem with a single crystal as the reference configuration. Using a multiplicative decomposition of the deformation gradient into lattice and plastic parts, i.e. F( X,t)= F L( X,t) F P( X,t), an initial stress-free polycrystal is constructed by imposing F L to be a piecewise constant rotation field R 0( X), and F P= R 0( X)T, thereby having F( X,0)= I, and zero elastic strain. This model serves as a precursor to higher order crystal plasticity models with grain boundary energy and evolution.

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

A boundary-value inverse model and its application to the calculation of tidal oscillation systems in the Western South Atlantic Ocean

International Nuclear Information System (INIS)

Miranda-Alonso, S.

1991-01-01

A Cauchy-Riemann problem is solved for the case of the linearized equations for long waves. The initial-values are amplitudes and phases measured at the coast. No boundary values are made use of. This inverse-problem is solved by starting the calculations at the coast and continuing outwards to the open ocean in a rectangular areas with one side at the coast and the other three at the open ocean. The initial values were expanded into the complex plane to get a platform to perform with the calculations. This non-well-posed problem was solved by means of two different mathematical techniques for comparison. The results produced with the inverse model were compared with those produced with a 'classical' model initialized at the three open boundaries with the results of the inverse model. The oscillating systems produced by both models were quite similar, giving validity to this invese modeling approach which should be a useful technique to solve problems when only initial values are known. (orig.)

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.

On two-point boundary correlations in the six-vertex model with domain wall boundary conditions

Science.gov (United States)

Colomo, F.; Pronko, A. G.

2005-05-01

The six-vertex model with domain wall boundary conditions on an N × N square lattice is considered. The two-point correlation function describing the probability of having two vertices in a given state at opposite (top and bottom) boundaries of the lattice is calculated. It is shown that this two-point boundary correlator is expressible in a very simple way in terms of the one-point boundary correlators of the model on N × N and (N - 1) × (N - 1) lattices. In alternating sign matrix (ASM) language this result implies that the doubly refined x-enumerations of ASMs are just appropriate combinations of the singly refined ones.

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.

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.

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.

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.

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

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.

A Probabilistic Approach for Breast Boundary Extraction in Mammograms

Directory of Open Access Journals (Sweden)

Hamed Habibi Aghdam

2013-01-01

Full Text Available The extraction of the breast boundary is crucial to perform further analysis of mammogram. Methods to extract the breast boundary can be classified into two categories: methods based on image processing techniques and those based on models. The former use image transformation techniques such as thresholding, morphological operations, and region growing. In the second category, the boundary is extracted using more advanced techniques, such as the active contour model. The problem with thresholding methods is that it is a hard to automatically find the optimal threshold value by using histogram information. On the other hand, active contour models require defining a starting point close to the actual boundary to be able to successfully extract the boundary. In this paper, we propose a probabilistic approach to address the aforementioned problems. In our approach we use local binary patterns to describe the texture around each pixel. In addition, the smoothness of the boundary is handled by using a new probability model. Experimental results show that the proposed method reaches 38% and 50% improvement with respect to the results obtained by the active contour model and threshold-based methods respectively, and it increases the stability of the boundary extraction process up to 86%.

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)

The Inverse Problem of Identification of Hydrogen Permeability Model

Directory of Open Access Journals (Sweden)

Yury V. Zaika

2018-01-01

Full Text Available One of the technological challenges for hydrogen materials science is the currently active search for structural materials with important applications (including the ITER project and gas-separation plants. One had to estimate the parameters of diffusion and sorption to numerically model the different scenarios and experimental conditions of the material usage (including extreme ones. The article presents boundary value problems of hydrogen permeability and thermal desorption with dynamical boundary conditions. A numerical method is developed for TDS spectrum simulation, where only integration of a nonlinear system of low order ordinary differential equations is required. The main final output of the article is a noise-resistant algorithm for solving the inverse problem of parametric identification for the aggregated experiment where desorption and diffusion are dynamically interrelated (without the artificial division of studies into the diffusion limited regime (DLR and the surface limited regime (SLR.

On relevant boundary perturbations of unitary minimal models

International Nuclear Information System (INIS)

Recknagel, A.; Roggenkamp, D.; Schomerus, V.

2000-01-01

We consider unitary Virasoro minimal models on the disk with Cardy boundary conditions and discuss deformations by certain relevant boundary operators, analogous to tachyon condensation in string theory. Concentrating on the least relevant boundary field, we can perform a perturbative analysis of renormalization group fixed points. We find that the systems always flow towards stable fixed points which admit no further (non-trivial) relevant perturbations. The new conformal boundary conditions are in general given by superpositions of 'pure' Cardy boundary conditions

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.

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).

NLIE of Dirichlet sine-Gordon model for boundary bound states

International Nuclear Information System (INIS)

Ahn, Changrim; Bajnok, Zoltan; Palla, Laszlo; Ravanini, Francesco

2008-01-01

We investigate boundary bound states of sine-Gordon model on the finite-size strip with Dirichlet boundary conditions. For the purpose we derive the nonlinear integral equation (NLIE) for the boundary excited states from the Bethe ansatz equation of the inhomogeneous XXZ spin 1/2 chain with boundary imaginary roots discovered by Saleur and Skorik. Taking a large volume (IR) limit we calculate boundary energies, boundary reflection factors and boundary Luescher corrections and compare with the excited boundary states of the Dirichlet sine-Gordon model first considered by Dorey and Mattsson. We also consider the short distance limit and relate the IR scattering data with that of the UV conformal field theory

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

On non-equilibrium states in QFT model with boundary interaction

International Nuclear Information System (INIS)

Bazhanov, Vladimir V.; Lukyanov, Sergei L.; Zamolodchikov, Alexander B.

1999-01-01

We prove that certain non-equilibrium expectation values in the boundary sine-Gordon model coincide with associated equilibrium-state expectation values in the systems which differ from the boundary sine-Gordon in that certain extra boundary degrees of freedom (q-oscillators) are added. Applications of this result to actual calculation of non-equilibrium characteristics of the boundary sine-Gordon model are also discussed

Classically integrable boundary conditions for symmetric-space sigma models

International Nuclear Information System (INIS)

MacKay, N.J.; Young, C.A.S.

2004-01-01

We investigate boundary conditions for the non-linear sigma model on the compact symmetric space G/H. The Poisson brackets and the classical local conserved charges necessary for integrability are preserved by boundary conditions which correspond to involutions which commute with the involution defining H. Applied to SO(3)/SO(2), the non-linear sigma model on S 2 , these yield the great circles as boundary submanifolds. Applied to GxG/G, they reproduce known results for the principal chiral model

Coupled wake boundary layer model of windfarms

Science.gov (United States)

Stevens, Richard; Gayme, Dennice; Meneveau, Charles

2014-11-01

We present a coupled wake boundary layer (CWBL) model that describes the distribution of the power output in a windfarm. The model couples the traditional, industry-standard wake expansion/superposition approach with a top-down model for the overall windfarm boundary layer structure. Wake models capture the effect of turbine positioning, while the top-down approach represents the interaction between the windturbine wakes and the atmospheric boundary layer. Each portion of the CWBL model requires specification of a parameter that is unknown a-priori. The wake model requires the wake expansion rate, whereas the top-down model requires the effective spanwise turbine spacing within which the model's momentum balance is relevant. The wake expansion rate is obtained by matching the mean velocity at the turbine from both approaches, while the effective spanwise turbine spacing is determined from the wake model. Coupling of the constitutive components of the CWBL model is achieved by iterating these parameters until convergence is reached. We show that the CWBL model predictions compare more favorably with large eddy simulation results than those made with either the wake or top-down model in isolation and that the model can be applied successfully to the Horns Rev and Nysted windfarms. The `Fellowships for Young Energy Scientists' (YES!) of the Foundation for Fundamental Research on Matter supported by NWO, and NSF Grant #1243482.

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

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)

Explicitly represented polygon wall boundary model for the explicit MPS method

Science.gov (United States)

Mitsume, Naoto; Yoshimura, Shinobu; Murotani, Kohei; Yamada, Tomonori

2015-05-01

This study presents an accurate and robust boundary model, the explicitly represented polygon (ERP) wall boundary model, to treat arbitrarily shaped wall boundaries in the explicit moving particle simulation (E-MPS) method, which is a mesh-free particle method for strong form partial differential equations. The ERP model expresses wall boundaries as polygons, which are explicitly represented without using the distance function. These are derived so that for viscous fluids, and with less computational cost, they satisfy the Neumann boundary condition for the pressure and the slip/no-slip condition on the wall surface. The proposed model is verified and validated by comparing computed results with the theoretical solution, results obtained by other models, and experimental results. Two simulations with complex boundary movements are conducted to demonstrate the applicability of the E-MPS method to the ERP model.

Measurement and modeling of radiation-induced grain boundary grain boundary segregation in stainless steels

International Nuclear Information System (INIS)

Bruemmer, S.M.; Charlot, L.A.; Simonen, E.P.

1995-08-01

Grain boundary radiation-induced segregation (RIS) in Fe-Ni-Cr stainless alloys has been measured and modelled as a function of irradiation temperature and dose. Heavy-ion irradiation was used to produce damage levels from 1 to 20 displacements per atom (dpa) at temperatures from 175 to 550 degrees C. Measured Fe, Ni, and Cr segregation increased sharply with irradiation dose (from 0 to 5 dpa) and temperature (from 175 to about 350 degrees C). However, grain boundary concentrations did not change significantly as dose or temperatures were further increased. Impurity segregation (Si and P) was also measured, but only Si enrichment appeared to be radiation-induced. Grain boundary Si levels peaked at an intermediate temperature of ∼325 degrees C reaching levels of ∼8 at. %. Equilibrium segregation of P was measured in the high-P alloys, but interfacial concentration did not increase with irradiation exposure. Examination of reported RIS in neutron-irradiated stainless steels revealed similar effects of irradiation dose on grain boundary compositional changes for both major alloying and impurity element's. The Inverse Kirkendall model accurately predicted major alloying element RIS in ion- and neutron-irradiated alloys over the wide range of temperature and dose conditions. In addition, preliminary calculations indicate that the Johnson-Lam model can reasonably estimate grain boundary Si enrichment if back diffusion is enhanced

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.

Thermodynamic Bethe ansatz for boundary sine-Gordon model

International Nuclear Information System (INIS)

Lee, Taejun; Rim, Chaiho

2003-01-01

(R-channel) TBA is elaborated to find the effective central charge dependence on the boundary parameters for the massless boundary sine-Gordon model with the coupling constant (8π)/β 2 =1+λ with λ a positive integer. Numerical analysis of the massless boundary TBA demonstrates that at an appropriate boundary parameter range (cusp point) there exists a singularity crossing phenomena and this effect should be included in TBA to have the right behavior of the effective central charge

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.

Boundary element method for modelling creep behaviour

International Nuclear Information System (INIS)

Zarina Masood; Shah Nor Basri; Abdel Majid Hamouda; Prithvi Raj Arora

2002-01-01

A two dimensional initial strain direct boundary element method is proposed to numerically model the creep behaviour. The boundary of the body is discretized into quadratic element and the domain into quadratic quadrilaterals. The variables are also assumed to have a quadratic variation over the elements. The boundary integral equation is solved for each boundary node and assembled into a matrix. This matrix is solved by Gauss elimination with partial pivoting to obtain the variables on the boundary and in the interior. Due to the time-dependent nature of creep, the solution has to be derived over increments of time. Automatic time incrementation technique and backward Euler method for updating the variables are implemented to assure stability and accuracy of results. A flowchart of the solution strategy is also presented. (Author)

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

Benthic boundary layer modelling studies

International Nuclear Information System (INIS)

Richards, K.J.

1984-01-01

A numerical model has been developed to study the factors which control the height of the benthic boundary layer in the deep ocean and the dispersion of a tracer within and directly above the layer. This report covers tracer clouds of horizontal scales of 10 to 100 km. The dispersion of a tracer has been studied in two ways. Firstly, a number of particles have been introduced into the flow. The trajectories of these particles provide information on dispersion rates. For flow conditions similar to those observed in the abyssal N.E. Atlantic the diffusivity of a tracer was found to be 5 x 10 6 cm 2 s -1 for a tracer within the boundary layer and 8 x 10 6 cm 2 s -1 for a tracer above the boundary layer. The results are in accord with estimates made from current meter measurements. The second method of studying dispersion was to calculate the evolution of individual tracer clouds. Clouds within and above the benthic boundary layer often show quite different behaviour from each other although the general structure of the clouds in the two regions were found to have no significant differences. (author)

Modelling stable atmospheric boundary layers over snow

NARCIS (Netherlands)

Sterk, H.A.M.

2015-01-01

Thesis entitled:

Modelling Stable Atmospheric Boundary Layers over Snow

H.A.M. Sterk

Wageningen, 29th of April, 2015

Summary

The emphasis of this thesis is on the understanding and forecasting of the Stable Boundary Layer (SBL) over snow-covered surfaces. SBLs

Twisted quantum double model of topological order with boundaries

Science.gov (United States)

Bullivant, Alex; Hu, Yuting; Wan, Yidun

2017-10-01

We generalize the twisted quantum double model of topological orders in two dimensions to the case with boundaries by systematically constructing the boundary Hamiltonians. Given the bulk Hamiltonian defined by a gauge group G and a 3-cocycle in the third cohomology group of G over U (1 ) , a boundary Hamiltonian can be defined by a subgroup K of G and a 2-cochain in the second cochain group of K over U (1 ) . The consistency between the bulk and boundary Hamiltonians is dictated by what we call the Frobenius condition that constrains the 2-cochain given the 3-cocyle. We offer a closed-form formula computing the ground-state degeneracy of the model on a cylinder in terms of the input data only, which can be naturally generalized to surfaces with more boundaries. We also explicitly write down the ground-state wave function of the model on a disk also in terms of the input data only.

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

Boundary rings and N=2 coset models

International Nuclear Information System (INIS)

Lerche, W.; Walcher, J.

2002-01-01

We investigate boundary states of N=2 coset models based on Grassmannians Gr(n,n+k), and find that the underlying intersection geometry is given by the fusion ring of U(n). This is isomorphic to the quantum cohomology ring of Gr(n,n+k+1), which in turn can be encoded in a 'boundary' superpotential whose critical points correspond to the boundary states. In this way the intersection properties can be represented in terms of a soliton graph that forms a generalized, Z n+k+1 symmetric McKay quiver. We investigate the spectrum of bound states and find that the rational boundary CFT produces only a small subset of the possible quiver representations

Modelling Velocity Spectra in the Lower Part of the Planetary Boundary Layer

DEFF Research Database (Denmark)

Olesen, H.R.; Larsen, Søren Ejling; Højstrup, Jørgen

1984-01-01

of the planetary boundary layer. Knowledge of the variation with stability of the (reduced) frequency f, for the spectral maximum is utilized in this modelling. Stable spectra may be normalized so that they adhere to one curve only, irrespective of stability, and unstable w-spectra may also be normalized to fit...... one curve. The problem of using filtered velocity variances when modelling spectra is discussed. A simplified procedure to provide a first estimate of the filter effect is given. In stable, horizontal velocity spectra, there is often a ‘gap’ at low frequencies. Using dimensional considerations...... and the spectral model previously derived, an expression for the gap frequency is found....

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.

On modeling the sound propagation through a lined duct with a modified Ingard-Myers boundary condition

Science.gov (United States)

Yang, Cheng; Fang, Yi; Zhao, Chao; Zhang, Xin

2018-06-01

A duct acoustics model is an essential component of an impedance eduction technique and its computation cost determines the impedance measurement efficiency. In this paper, a model is developed for the sound propagation through a lined duct carrying a uniform mean flow. In contrast to many existing models, the interface between the liner and the duct field is defined with a modified Ingard-Myers boundary condition that takes account of the effect of the boundary layer above the liner. A mode-matching method is used to couple the unlined and lined duct segments for the model development. For the lined duct segment, the eigenvalue problem resulted from the modified boundary condition is solved by an integration scheme which, on the one hand, allows the lined duct modes to be computed in an efficient manner, and on the other hand, orders the modes automatically. The duct acoustics model developed from the solved lined duct modes is shown to converge more rapidly than the one developed from the rigid-walled duct modes. Validation against the experiment data in the literature shows that the proposed model is able to predict more accurately the liner performance measured by the two-source method. This, however, cannot be made by a duct acoustics model associated with the conventional Ingard-Myers boundary condition. The proposed model has the potential to be integrated into an impedance eduction technique for more reliable liner measurement.

Boundary Value Problems for a Super-Sublinear Asymmetric Oscillator: The Exact Number of Solutions

Directory of Open Access Journals (Sweden)

Armands Gritsans

2013-01-01

Full Text Available Properties of asymmetric oscillator described by the equation (i, where and , are studied. A set of such that the problem (i, (ii, and (iii have a nontrivial solution, is called α-spectrum. We give full description of α-spectra in terms of solution sets and solution surfaces. The exact number of nontrivial solutions of the two-parameter Dirichlet boundary value problem (i, and (ii is given.

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

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.

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

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

Boundary scattering in the ϕ{sup 4} model

Energy Technology Data Exchange (ETDEWEB)

Dorey, Patrick [Department of Mathematical Sciences, Durham University,Durham DH1 3LE (United Kingdom); Halavanau, Aliaksei [Department of Theoretical Physics and Astrophysics,BSU, Minsk Independence Avenue 4 (Belarus); Fermi National Laboratory,Pine St. and Kirk Rd., ZIP 60511, Mail Station 221, Batavia, Illinois (United States); Mercer, James [Department of Mathematical Sciences, Durham University,Durham DH1 3LE (United Kingdom); Deloitte MCS Limited,Hill House, 1 Little New Street, London, EC4A 3TR (United Kingdom); Romanczukiewicz, Tomasz [Institute of Physics, Jagiellonian University,Lojasiewicza 11, 30-348 Krakow (Poland); Shnir, Yasha [Department of Theoretical Physics and Astrophysics,BSU, Minsk Independence Avenue 4 (Belarus); BLTP, JINR,141980 Dubna (Russian Federation); Institute of Physics, Oldenburg University,Postfach 2503 D-26111 Oldenburg (Germany)

2017-05-19

We study boundary scattering in the ϕ{sup 4} model on a half-line with a one-parameter family of Neumann-type boundary conditions. A rich variety of phenomena is observed, which extends previously-studied behaviour on the full line to include regimes of near-elastic scattering, the restoration of a missing scattering window, and the creation of a kink or oscillon through the collision-induced decay of a metastable boundary state. We also study the decay of the vibrational boundary mode, and explore different scenarios for its relaxation and for the creation of kinks.

Beyond the climate envelope: using trait filtering models to predict biome boundaries from plant physiology.

Science.gov (United States)

Fisher, R.; Hoffmann, W. A.; Muszala, S.

2014-12-01

The introduction of second-generation dynamic vegetation models - which simulate the distribution of light resources between plant types along the vertical canopy profile, and therefore facilitate the representation of plant competition explicitly - is a large increase in the complexity and fidelity with which the terrestrial biosphere is abstracted into Earth System Models. In this new class of model, biome boundaries are predicted as the emergent properties of plant physiology, and are therefore sensitive to the high-dimensional parameterizations of plant functional traits. These new approaches offer the facility to quantitatively test ecophysiological hypotheses of plant distribution at large scales, a field which remains surprisingly under-developed. Here we describe experiments conducted with the Community Land Model Ecosystem Demography component, CLM(ED), in which we reduce the complexity of the problem by testing how individual plant functional trait changes to control the location of biome boundaries between functional types. Specifically, we investigate which physiological trade-offs determine the boundary between frequently burned savanna and forest biomes, and attempt to distinguish how each strategic life-history trade-off (carbon storage, bark investment, re-sprouting strategy) contributes towards the maintenance of sharp geographical gradients between fire adapted and typically inflammable closed canopy ecosystems. This study forms part of the planning for a model-inspired fire manipulation experiment at the cerrado-forest boundary in South-Eastern Brazil, and the results will be used to guide future data-collection and analysis strategies.

Nested Bethe Ansatz for Spin Ladder Model with Open Boundary Conditions

International Nuclear Information System (INIS)

Wu Junfang; Zhang Chunmin; Yue Ruihong; Li Runling

2005-01-01

The nested Bethe ansatz (BA) method is applied to find the eigenvalues and the eigenvectors of the transfer matrix for spin-ladder model with open boundary conditions. Based on the reflection equation, we find the general diagonal solution, which determines the general boundary interaction in the Hamiltonian. We introduce the spin-ladder model with open boundary conditions. By finding the solution K ± of the reflection equation which determines the nontrivial boundary terms in the Hamiltonian, we diagonalize the transfer matrix of the spin-ladder model with open boundary conditions in the framework of nested BA.

Modelling classroom conditions with different boundary conditions

DEFF Research Database (Denmark)

Marbjerg, Gerd Høy; Jeong, Cheol-Ho; Brunskog, Jonas

2014-01-01

A model that combines image source modelling and acoustical radiosity with complex boundary condition, thus including phase shifts on reflection has been developed. The model is called PARISM (Phased Acoustical Radiosity and Image Source Model). It has been developed in order to be able to model...

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

Reflection of equatorial Kelvin waves at eastern ocean boundaries Part I: hypothetical boundaries

Directory of Open Access Journals (Sweden)

J. Soares

1999-06-01

Full Text Available A baroclinic shallow-water model is developed to investigate the effect of the orientation of the eastern ocean boundary on the behavior of equatorial Kelvin waves. The model is formulated in a spherical polar coordinate system and includes dissipation and non-linear terms, effects which have not been previously included in analytical approaches to the problem. Both equatorial and middle latitude response are considered given the large latitudinal extent used in the model. Baroclinic equatorial Kelvin waves of intraseasonal, seasonal and annual periods are introduced into the domain as pulses of finite width. Their subsequent reflection, transmission and dissipation are investigated. It is found that dissipation is very important for the transmission of wave energy along the boundary and for reflections from the boundary. The dissipation was found to be dependent not only on the presence of the coastal Kelvin waves in the domain, but also on the period of these coastal waves. In particular the dissipation increases with wave period. It is also shown that the equatorial β-plane approximation can allow an anomalous generation of Rossby waves at higher latitudes. Nonlinearities generally have a small effect on the solutions, within the confines of this model.Key words. Oceanography: general (equatorial oceanography; numerical modeling · Oceanography: physical (eastern boundary currents

Fused integrable lattice models with quantum impurities and open boundaries

International Nuclear Information System (INIS)

Doikou, Anastasia

2003-01-01

The alternating integrable spin chain and the RSOS(q 1 ,q 2 ;p) model in the presence of a quantum impurity are investigated. The boundary free energy due to the impurity is derived, the ratios of the corresponding g functions at low and high temperature are specified and their relevance to boundary flows in unitary minimal and generalized coset models is discussed. Finally, the alternating spin chain with diagonal and non-diagonal integrable boundaries is studied, and the corresponding boundary free energy and g functions are derived

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.

Development of stress boundary conditions in smoothed particle hydrodynamics (SPH) for the modeling of solids deformation

Science.gov (United States)

Douillet-Grellier, Thomas; Pramanik, Ranjan; Pan, Kai; Albaiz, Abdulaziz; Jones, Bruce D.; Williams, John R.

2017-10-01

This paper develops a method for imposing stress boundary conditions in smoothed particle hydrodynamics (SPH) with and without the need for dummy particles. SPH has been used for simulating phenomena in a number of fields, such as astrophysics and fluid mechanics. More recently, the method has gained traction as a technique for simulation of deformation and fracture in solids, where the meshless property of SPH can be leveraged to represent arbitrary crack paths. Despite this interest, application of boundary conditions within the SPH framework is typically limited to imposed velocity or displacement using fictitious dummy particles to compensate for the lack of particles beyond the boundary interface. While this is enough for a large variety of problems, especially in the case of fluid flow, for problems in solid mechanics there is a clear need to impose stresses upon boundaries. In addition to this, the use of dummy particles to impose a boundary condition is not always suitable or even feasibly, especially for those problems which include internal boundaries. In order to overcome these difficulties, this paper first presents an improved method for applying stress boundary conditions in SPH with dummy particles. This is then followed by a proposal of a formulation which does not require dummy particles. These techniques are then validated against analytical solutions to two common problems in rock mechanics, the Brazilian test and the penny-shaped crack problem both in 2D and 3D. This study highlights the fact that SPH offers a good level of accuracy to solve these problems and that results are reliable. This validation work serves as a foundation for addressing more complex problems involving plasticity and fracture propagation.

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.

ON A PARABOLIC FREE BOUNDARY EQUATION MODELING PRICE FORMATION

KAUST Repository

MARKOWICH, P. A.

2009-10-01

We discuss existence and uniqueness of solutions for a one-dimensional parabolic evolution equation with a free boundary. This problem was introduced by Lasry and Lions as description of the dynamical formation of the price of a trading good. Short time existence and uniqueness is established by a contraction argument. Then we discuss the issue of global-in-time-extension of the local solution which is closely related to the regularity of the free boundary. We also present numerical results. © 2009 World Scientific Publishing Company.

ON A PARABOLIC FREE BOUNDARY EQUATION MODELING PRICE FORMATION

KAUST Repository

MARKOWICH, P. A.; MATEVOSYAN, N.; PIETSCHMANN, J.-F.; WOLFRAM, M.-T.

2009-01-01

We discuss existence and uniqueness of solutions for a one-dimensional parabolic evolution equation with a free boundary. This problem was introduced by Lasry and Lions as description of the dynamical formation of the price of a trading good. Short time existence and uniqueness is established by a contraction argument. Then we discuss the issue of global-in-time-extension of the local solution which is closely related to the regularity of the free boundary. We also present numerical results. © 2009 World Scientific Publishing Company.

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.

Finite element approximation to a model problem of transonic flow

International Nuclear Information System (INIS)

Tangmanee, S.

1986-12-01

A model problem of transonic flow ''the Tricomi equation'' in Ω is contained in IR 2 bounded by the rectangular-curve boundary is posed in the form of symmetric positive differential equations. The finite element method is then applied. When the triangulation of Ω-bar is made of quadrilaterals and the approximation space is the Lagrange polynomial, we get the error estimates. 14 refs, 1 fig

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. $$

Stability of stationary solutions for inflow problem on the micropolar fluid model

Science.gov (United States)

Yin, Haiyan

2017-04-01

In this paper, we study the asymptotic behavior of solutions to the initial boundary value problem for the micropolar fluid model in a half-line R+:=(0,∞). We prove that the corresponding stationary solutions of the small amplitude to the inflow problem for the micropolar fluid model are time asymptotically stable under small H1 perturbations in both the subsonic and degenerate cases. The microrotation velocity brings us some additional troubles compared with Navier-Stokes equations in the absence of the microrotation velocity. The proof of asymptotic stability is based on the basic energy method.

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

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$.

Convergence of a continuous BGK model for initial boundary-value problems for conservation laws

Directory of Open Access Journals (Sweden)

Driss Seghir

2001-11-01

Full Text Available We consider a scalar conservation law in the quarter plane. This equation is approximated in a continuous kinetic Bhatnagar-Gross-Krook (BGK model. The convergence of the model towards the unique entropy solution is established in the space of functions of bounded variation, using kinetic entropy inequalities, without special restriction on the flux nor on the equilibrium problem's data. As an application, we establish the hydrodynamic limit for a $2imes2$ relaxation system with general data. Also we construct a new family of convergent continuous BGK models with simple maxwellians different from the $chi$ models.

A model of anelastic relaxation associated with polygonization boundary

International Nuclear Information System (INIS)

Yan, S.C.

1990-01-01

A model of anelastic relaxation associated with polygonization boundary is proposed in order to explain internal friction peaks and other experimental phenomena observed recently. The model, which is referred to as vacancy-thermal jog model, shows that under conditions of high temperature and low applied stress with lower frequencies of vibration, thermal jog pairs are generated on dislocation segments of the boundaries. These jogs are in saturation with vacancies in the vicinity of them, and the vacancy current due to the concentration gradient of vacancy drifts among the boundaries. As a result, a diffusional creep is produced and a part of energy is dissipated. For vacancy drift, it is required that the thermal jogs emit (absorb) vacancies, which brings climbing bow of segments into operation, and another part of energy is dissipated so that there are two parts of energy dissipated in the strain process connected with polygonization boundary. Based on this point of view, the mathematical expressions of internal friction and modulus defect associated with polygonization boundary were subsequently derived and found to be in satisfactory agreement with experiments. (author). 13 refs, 6 figs

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)

Crossing boundaries in a collaborative modeling workspace

Science.gov (United States)

Morisette, Jeffrey T.; Cravens, Amanda; Miller, Brian W.; Talbert, Marian; Talbert, Colin; Jarnevich, Catherine S.; Fink, Michelle; Decker, Karin; Odell, Eric

2017-01-01

There is substantial literature on the importance of bridging across disciplinary and science–management boundaries. One of the ways commonly suggested to cross boundaries is for participants from both sides of the boundary to jointly produce information (i.e., knowledge co-production). But simply providing tools or bringing people together in the same room is not sufficient. Here we present a case study documenting the mechanisms by which managers and scientists collaborated to incorporate climate change projections into Colorado’s State Wildlife Action Plan. A critical component of the project was the use of a collaborative modeling and visualization workspace: the U.S. Geological Survey’s Resource for Advanced Modeling (RAM). Using video analysis and pre/post surveys from this case study, we examine how the RAM facilitated cognitive and social processes that co-produced a more salient and credible end product. This case provides practical suggestions to scientists and practitioners who want to implement actionable science.

A finite element model for the quench front evolution problem

International Nuclear Information System (INIS)

Folescu, J.; Galeao, A.C.N.R.; Carmo, E.G.D. do.

1985-01-01

A model for the rewetting problem associated with the loss of coolant accident in a PWR reactor is proposed. A variational formulation for the time-dependent heat conduction problem on fuel rod cladding is used, and appropriate boundary conditions are assumed in order to simulate the thermal interaction between the fuel rod cladding and the fluid. A numerical procedure which uses the finite element method for the spatial discretization and a Crank-Nicolson-like method for the step-by-step integration is developed. Some numerical results are presented showing the quench front evolution and its stationary profile. (Author) [pt

Switching moving boundary models for two-phase flow evaporators and condensers

Science.gov (United States)

Bonilla, Javier; Dormido, Sebastián; Cellier, François E.

2015-03-01

The moving boundary method is an appealing approach for the design, testing and validation of advanced control schemes for evaporators and condensers. When it comes to advanced control strategies, not only accurate but fast dynamic models are required. Moving boundary models are fast low-order dynamic models, and they can describe the dynamic behavior with high accuracy. This paper presents a mathematical formulation based on physical principles for two-phase flow moving boundary evaporator and condenser models which support dynamic switching between all possible flow configurations. The models were implemented in a library using the equation-based object-oriented Modelica language. Several integrity tests in steady-state and transient predictions together with stability tests verified the models. Experimental data from a direct steam generation parabolic-trough solar thermal power plant is used to validate and compare the developed moving boundary models against finite volume models.

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...

Initializing a Mesoscale Boundary-Layer Model with Radiosonde Observations

Science.gov (United States)

Berri, Guillermo J.; Bertossa, Germán

2018-01-01

A mesoscale boundary-layer model is used to simulate low-level regional wind fields over the La Plata River of South America, a region characterized by a strong daily cycle of land-river surface-temperature contrast and low-level circulations of sea-land breeze type. The initial and boundary conditions are defined from a limited number of local observations and the upper boundary condition is taken from the only radiosonde observations available in the region. The study considers 14 different upper boundary conditions defined from the radiosonde data at standard levels, significant levels, level of the inversion base and interpolated levels at fixed heights, all of them within the first 1500 m. The period of analysis is 1994-2008 during which eight daily observations from 13 weather stations of the region are used to validate the 24-h surface-wind forecast. The model errors are defined as the root-mean-square of relative error in wind-direction frequency distribution and mean wind speed per wind sector. Wind-direction errors are greater than wind-speed errors and show significant dispersion among the different upper boundary conditions, not present in wind speed, revealing a sensitivity to the initialization method. The wind-direction errors show a well-defined daily cycle, not evident in wind speed, with the minimum at noon and the maximum at dusk, but no systematic deterioration with time. The errors grow with the height of the upper boundary condition level, in particular wind direction, and double the errors obtained when the upper boundary condition is defined from the lower levels. The conclusion is that defining the model upper boundary condition from radiosonde data closer to the ground minimizes the low-level wind-field errors throughout the region.

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.

An outgoing energy flux boundary condition for finite difference ICRP antenna models

International Nuclear Information System (INIS)

Batchelor, D.B.; Carter, M.D.

1992-11-01

For antennas at the ion cyclotron range of frequencies (ICRF) modeling in vacuum can now be carried out to a high level of detail such that shaping of the current straps, isolating septa, and discrete Faraday shield structures can be included. An efficient approach would be to solve for the fields in the vacuum region near the antenna in three dimensions by finite methods and to match this solution at the plasma-vacuum interface to a solution obtained in the plasma region in one dimension by Fourier methods. This approach has been difficult to carry out because boundary conditions must be imposed at the edge of the finite difference grid on a point-by-point basis, whereas the condition for outgoing energy flux into the plasma is known only in terms of the Fourier transform of the plasma fields. A technique is presented by which a boundary condition can be imposed on the computational grid of a three-dimensional finite difference, or finite element, code by constraining the discrete Fourier transform of the fields at the boundary points to satisfy an outgoing energy flux condition appropriate for the plasma. The boundary condition at a specific grid point appears as a coupling to other grid points on the boundary, with weighting determined by a kemel calctdated from the plasma surface impedance matrix for the various plasma Fourier modes. This boundary condition has been implemented in a finite difference solution of a simple problem in two dimensions, which can also be solved directly by Fourier transformation. Results are presented, and it is shown that the proposed boundary condition does enforce outgoing energy flux and yields the same solution as is obtained by Fourier methods

Diagnostic analysis of turbulent boundary layer data by a trivariate Lagrangian partitioning method

Energy Technology Data Exchange (ETDEWEB)

Welsh, P.T. [Florida State Univ., Tallahassee, FL (United States)

1994-12-31

The rapid scientific and technological advances in meteorological theory and modeling predominantly have occurred on the large (or synoptic) scale flow characterized by the extratropical cyclone. Turbulent boundary layer flows, in contrast, have been slower in developing both theoretically and in accuracy for several reasons. There are many existing problems in boundary layer models, among them are limits to computational power available, the inability to handle countergradient fluxes, poor growth matching to real boundary layers, and inaccuracy in calculating the diffusion of scalar concentrations. Such transport errors exist within the boundary layer as well as into the free atmosphere above. This research uses a new method, which can provide insight into these problems, and ultimately improve boundary layer models. There are several potential applications of the insights provided by this approach, among them are estimation of cloud contamination of satellite remotely sensed surface parameters, improved flux and vertical transport calculations, and better understanding of the diurnal boundary layer growth process and its hysteresis cycle.

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

Micromechanical Modeling of Grain Boundaries Damage in a Copper Alloy Under Creep

International Nuclear Information System (INIS)

Voese, Markus

2015-01-01

In order to include the processes on the scale of the grain structure into the description of the creep behaviour of polycrystalline materials, the damage development of a single grain boundary has been initially investigated in the present work. For this purpose, a special simulationmethod has been used, whose resolution procedure based on holomorphic functions. The mechanisms taken into account for the simulations include nucleation, growth by grain boundary diffusion, coalescence and shrinkage until complete sintering of grain boundary cavities. These studies have then been used to develop a simplified cavitation model, which describes the grain boundary damage by two state variables and the time-dependent development by a mechanism-oriented rate formulation. To include the influence of grain boundaries within continuum mechanical considerations of polycrystals, an interface model has been developed, that incorporates both damage according to the simplified cavitation model and grain boundary sliding in dependence of a phenomenological grain boundary viscosity. Furthermore a micromechanical model of a polycrystal has been developed that allows to include a material's grain structure into the simulation of the creep behaviour by means of finite element simulations. Thereby, the deformations of individual grains are expressed by a viscoplastic single crystal model and the grain boundaries are described by the proposed interface model. The grain structure is represented by a finite element model, in which the grain boundaries are modelled by cohesive elements. From the evaluation of experimental creep data, the micromechanical model of a polycrystal has been calibrated for a copper-antimony alloy at a temperature of 823 K. Thereby, the adjustment of the single crystal model has been carried out on the basis of creep rates of pure copper single crystal specimens. The experimental determination of grain boundary sliding and grain boundary porosity for coarse

Modeling Simple Driving Tasks with a One-Boundary Diffusion Model

Science.gov (United States)

Ratcliff, Roger; Strayer, David

2014-01-01

A one-boundary diffusion model was applied to the data from two experiments in which subjects were performing a simple simulated driving task. In the first experiment, the same subjects were tested on two driving tasks using a PC-based driving simulator and the psychomotor vigilance test (PVT). The diffusion model fit the response time (RT) distributions for each task and individual subject well. Model parameters were found to correlate across tasks which suggests common component processes were being tapped in the three tasks. The model was also fit to a distracted driving experiment of Cooper and Strayer (2008). Results showed that distraction altered performance by affecting the rate of evidence accumulation (drift rate) and/or increasing the boundary settings. This provides an interpretation of cognitive distraction whereby conversing on a cell phone diverts attention from the normal accumulation of information in the driving environment. PMID:24297620

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

The viscous lee wave problem and its implications for ocean modelling

Science.gov (United States)

Shakespeare, Callum J.; Hogg, Andrew McC.

2017-05-01

Ocean circulation models employ 'turbulent' viscosity and diffusivity to represent unresolved sub-gridscale processes such as breaking internal waves. Computational power has now advanced sufficiently to permit regional ocean circulation models to be run at sufficiently high (100 m-1 km) horizontal resolution to resolve a significant part of the internal wave spectrum. Here we develop theory for boundary generated internal waves in such models, and in particular, where the waves dissipate their energy. We focus specifically on the steady lee wave problem where stationary waves are generated by a large-scale flow acting across ocean bottom topography. We generalise the energy flux expressions of [Bell, T., 1975. Topographically generated internal waves in the open ocean. J. Geophys. Res. 80, 320-327] to include the effect of arbitrary viscosity and diffusivity. Applying these results for realistic parameter choices we show that in the present generation of models with O(1) m2s-1 horizontal viscosity/diffusivity boundary-generated waves will inevitably dissipate the majority of their energy within a few hundred metres of the boundary. This dissipation is a direct consequence of the artificially high viscosity/diffusivity, which is not always physically justified in numerical models. Hence, caution is necessary in comparing model results to ocean observations. Our theory further predicts that O(10-2) m2s-1 horizontal and O(10-4) m2s-1 vertical viscosity/diffusivity is required to achieve a qualitatively inviscid representation of internal wave dynamics in ocean models.

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.

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.

Exact solutions to plaquette Ising models with free and periodic boundaries

International Nuclear Information System (INIS)

Mueller, Marco; Johnston, Desmond A.; Janke, Wolfhard

2017-01-01

An anisotropic limit of the 3d plaquette Ising model, in which the plaquette couplings in one direction were set to zero, was solved for free boundary conditions by Suzuki (1972) , who later dubbed it the fuki-nuke, or “no-ceiling”, model. Defining new spin variables as the product of nearest-neighbour spins transforms the Hamiltonian into that of a stack of (standard) 2d Ising models and reveals the planar nature of the magnetic order, which is also present in the fully isotropic 3d plaquette model. More recently, the solution of the fuki-nuke model was discussed for periodic boundary conditions, which require a different approach to defining the product spin transformation, by Castelnovo et al. (2010) . We clarify the exact relation between partition functions with free and periodic boundary conditions expressed in terms of original and product spin variables for the 2d plaquette and 3d fuki-nuke models, noting that the differences are already present in the 1d Ising model. In addition, we solve the 2d plaquette Ising model with helical boundary conditions. The various exactly solved examples illustrate how correlations can be induced in finite systems as a consequence of the choice of boundary conditions.

Investigation of turbulence models with compressibility corrections for hypersonic boundary flows

Directory of Open Access Journals (Sweden)

Han Tang

2015-12-01

Full Text Available The applications of pressure work, pressure-dilatation, and dilatation-dissipation (Sarkar, Zeman, and Wilcox models to hypersonic boundary flows are investigated. The flat plate boundary layer flows of Mach number 5–11 and shock wave/boundary layer interactions of compression corners are simulated numerically. For the flat plate boundary layer flows, original turbulence models overestimate the heat flux with Mach number high up to 10, and compressibility corrections applied to turbulence models lead to a decrease in friction coefficients and heating rates. The pressure work and pressure-dilatation models yield the better results. Among the three dilatation-dissipation models, Sarkar and Wilcox corrections present larger deviations from the experiment measurement, while Zeman correction can achieve acceptable results. For hypersonic compression corner flows, due to the evident increase of turbulence Mach number in separation zone, compressibility corrections make the separation areas larger, thus cannot improve the accuracy of calculated results. It is unreasonable that compressibility corrections take effect in separation zone. Density-corrected model by Catris and Aupoix is suitable for shock wave/boundary layer interaction flows which can improve the simulation accuracy of the peak heating and have a little influence on separation zone.

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

Rayleigh's, Stoneley's, and Scholte's Interface Waves in Elastic Models Using a Boundary Element Method

Directory of Open Access Journals (Sweden)

Esteban Flores-Mendez

2012-01-01

Full Text Available This work is focused on studying interface waves for three canonical models, that is, interfaces formed by vacuum-solid, solid-solid, and liquid-solid. These interfaces excited by dynamic loads cause the emergence of Rayleigh's, Stoneley's, and Scholte's waves, respectively. To perform the study, the indirect boundary element method is used, which has proved to be a powerful tool for numerical modeling of problems in elastodynamics. In essence, the method expresses the diffracted wave field of stresses, pressures, and displacements by a boundary integral, also known as single-layer representation, whose shape can be regarded as a Fredholm's integral representation of second kind and zero order. This representation can be considered as an exemplification of Huygens' principle, which is equivalent to Somigliana's representation theorem. Results in frequency domain for the three types of interfaces are presented; then, using the fourier discrete transform, we derive the results in time domain, where the emergence of interface waves is highlighted.

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.

Self-consistent model of the low-latitude boundary layer

International Nuclear Information System (INIS)

Phan, T.D.; Sonnerup, B.U.Oe.; Lotko, W.

1989-01-01

A simple two-dimensional, steady state, viscous model of the dawnside and duskside low-latitude boundary layer (LLBL) has been developed. It incorporates coupling to the ionosphere via field-aligned currents and associated field-aligned potential drops, governed by a simple conductance law, and it describes boundary layer currents, magnetic fields, and plasma flow in a self-consistent manner. The magnetic field induced by these currents leads to two effects: (1) a diamagnetic depression of the magnetic field in the equatorial region and (2) bending of the field lines into parabolas in the xz plane with their vertices in the equatorial plane, at z = 0, and pointing in the flow direction, i.e., tailward. Both effects are strongest at the magnetopause edge of the boundary layer and vanish at the magnetospheric edge. The diamagnetic depression corresponds to an excess of plasma pressure in the equatorial boundary layer near the magnetopause. The boundary layer structure is governed by a fourth-order, nonlinear, ordinary differential equation in which one nondimensional parameter, the Hartmann number M, appears. A second parameter, introduced via the boundary conditions, is a nondimensional flow velocity v 0 * at the magnetopause. Numerical results from the model are presented and the possible use of observations to determine the model parameters is discussed. The main new contribution of the study is to provide a better description of the field and plasma configuration in the LLBL itself and to clarify in quantitative terms the circumstances in which induced magnetic fields become important

RANS Modeling of Benchmark Shockwave / Boundary Layer Interaction Experiments

Science.gov (United States)

Georgiadis, Nick; Vyas, Manan; Yoder, Dennis

2010-01-01

This presentation summarizes the computations of a set of shock wave / turbulent boundary layer interaction (SWTBLI) test cases using the Wind-US code, as part of the 2010 American Institute of Aeronautics and Astronautics (AIAA) shock / boundary layer interaction workshop. The experiments involve supersonic flows in wind tunnels with a shock generator that directs an oblique shock wave toward the boundary layer along one of the walls of the wind tunnel. The Wind-US calculations utilized structured grid computations performed in Reynolds-averaged Navier-Stokes mode. Three turbulence models were investigated: the Spalart-Allmaras one-equation model, the Menter Shear Stress Transport wavenumber-angular frequency two-equation model, and an explicit algebraic stress wavenumber-angular frequency formulation. Effects of grid resolution and upwinding scheme were also considered. The results from the CFD calculations are compared to particle image velocimetry (PIV) data from the experiments. As expected, turbulence model effects dominated the accuracy of the solutions with upwinding scheme selection indicating minimal effects.!

A new magneto-cardiogram study using a vector model with a virtual heart and the boundary element method

International Nuclear Information System (INIS)

Zhang Chen; Lu Hong; Hua Ning; Tang Xue-Zheng; Tang Fa-Kuan; Shou Guo-Fa; Xia Ling; Ma Ping

2013-01-01

A cardiac vector model is presented and verified, and then the forward problem for cardiac magnetic fields and electric potential are discussed based on this model and the realistic human torso volume conductor model, including lungs. A torso—cardiac vector model is used for a 12-lead electrocardiographic (ECG) and magneto-cardiogram (MCG) simulation study by using the boundary element method (BEM). Also, we obtain the MCG wave picture using a compound four-channel HT c ·SQUID system in a magnetically shielded room. By comparing the simulated results and experimental results, we verify the cardiac vector model and then do a preliminary study of the forward problem of MCG and ECG. Therefore, the results show that the vector model is reasonable in cardiac electrophysiology. (general)

A non-reflecting boundary for use in a finite element beam model of a railway track

Science.gov (United States)

Yang, Jiannan; Thompson, David J.

2015-02-01

Some beam-like structures such as a railway track are effectively infinite in nature. Analytical solutions exist for simple structures but numerical methods like the finite element (FE) method are often employed to study more complicated problems. However, when the FE method is used for structures of infinite extent it is essential to introduce artificial boundaries to limit the area of computation. Here, a non-reflecting boundary is developed using a damped tapered tip for application in a finite element model representing an infinite supported beam. The FE model of the tapered tip is validated against an analytical model based on Bessel functions. The reflection characteristics of the FE tapered tip are quantified using a wave/FE superposition method. It is shown that the damped tapered tip is much more effective than its constant counterpart and achieves reduction of the model size. The damped tapered tip is applied to a simple FE railway track model and good agreement is found when its point mobility is compared with an analytical infinite track model.

Statistical surrogate model based sampling criterion for stochastic global optimization of problems with constraints

Energy Technology Data Exchange (ETDEWEB)

Cho, Su Gil; Jang, Jun Yong; Kim, Ji Hoon; Lee, Tae Hee [Hanyang University, Seoul (Korea, Republic of); Lee, Min Uk [Romax Technology Ltd., Seoul (Korea, Republic of); Choi, Jong Su; Hong, Sup [Korea Research Institute of Ships and Ocean Engineering, Daejeon (Korea, Republic of)

2015-04-15

Sequential surrogate model-based global optimization algorithms, such as super-EGO, have been developed to increase the efficiency of commonly used global optimization technique as well as to ensure the accuracy of optimization. However, earlier studies have drawbacks because there are three phases in the optimization loop and empirical parameters. We propose a united sampling criterion to simplify the algorithm and to achieve the global optimum of problems with constraints without any empirical parameters. It is able to select the points located in a feasible region with high model uncertainty as well as the points along the boundary of constraint at the lowest objective value. The mean squared error determines which criterion is more dominant among the infill sampling criterion and boundary sampling criterion. Also, the method guarantees the accuracy of the surrogate model because the sample points are not located within extremely small regions like super-EGO. The performance of the proposed method, such as the solvability of a problem, convergence properties, and efficiency, are validated through nonlinear numerical examples with disconnected feasible regions.

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!

Exact Solution of the Six-Vertex Model with Domain Wall Boundary Conditions. Disordered Phase

CERN Document Server

Bleher, P M

2005-01-01

The six-vertex model, or the square ice model, with domain wall boundary conditions (DWBC) has been introduced and solved for finite $N$ by Korepin and Izergin. The solution is based on the Yang-Baxter equations and it represents the free energy in terms of an $N\\times N$ Hankel determinant. Paul Zinn-Justin observed that the Izergin-Korepin formula can be re-expressed in terms of the partition function of a random matrix model with a nonpolynomial interaction. We use this observation to obtain the large $N$ asymptotics of the six-vertex model with DWBC in the disordered phase. The solution is based on the Riemann-Hilbert approach and the Deift-Zhou nonlinear steepest descent method. As was noticed by Kuperberg, the problem of enumeration of alternating sign matrices (the ASM problem) is a special case of the the six-vertex model. We compare the obtained exact solution of the six-vertex model with known exact results for the 1, 2, and 3 enumerations of ASMs, and also with the exact solution on the so-called f...

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.

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.

Spectral combination of spherical gravitational curvature boundary-value problems

Science.gov (United States)

PitoÅák, Martin; Eshagh, Mehdi; Šprlák, Michal; Tenzer, Robert; Novák, Pavel

2018-04-01

Four solutions of the spherical gravitational curvature boundary-value problems can be exploited for the determination of the Earth's gravitational potential. In this article we discuss the combination of simulated satellite gravitational curvatures, i.e., components of the third-order gravitational tensor, by merging these solutions using the spectral combination method. For this purpose, integral estimators of biased- and unbiased-types are derived. In numerical studies, we investigate the performance of the developed mathematical models for the gravitational field modelling in the area of Central Europe based on simulated satellite measurements. Firstly, we verify the correctness of the integral estimators for the spectral downward continuation by a closed-loop test. Estimated errors of the combined solution are about eight orders smaller than those from the individual solutions. Secondly, we perform a numerical experiment by considering the Gaussian noise with the standard deviation of 6.5× 10-17 m-1s-2 in the input data at the satellite altitude of 250 km above the mean Earth sphere. This value of standard deviation is equivalent to a signal-to-noise ratio of 10. Superior results with respect to the global geopotential model TIM-r5 are obtained by the spectral downward continuation of the vertical-vertical-vertical component with the standard deviation of 2.104 m2s-2, but the root mean square error is the largest and reaches 9.734 m2s-2. Using the spectral combination of all gravitational curvatures the root mean square error is more than 400 times smaller but the standard deviation reaches 17.234 m2s-2. The combination of more components decreases the root mean square error of the corresponding solutions while the standard deviations of the combined solutions do not improve as compared to the solution from the vertical-vertical-vertical component. The presented method represents a weight mean in the spectral domain that minimizes the root mean square error

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.

An Efficient Constraint Boundary Sampling Method for Sequential RBDO Using Kriging Surrogate Model

Energy Technology Data Exchange (ETDEWEB)

Kim, Jihoon; Jang, Junyong; Kim, Shinyu; Lee, Tae Hee [Hanyang Univ., Seoul (Korea, Republic of); Cho, Sugil; Kim, Hyung Woo; Hong, Sup [Korea Research Institute of Ships and Ocean Engineering, Busan (Korea, Republic of)

2016-06-15

Reliability-based design optimization (RBDO) requires a high computational cost owing to its reliability analysis. A surrogate model is introduced to reduce the computational cost in RBDO. The accuracy of the reliability depends on the accuracy of the surrogate model of constraint boundaries in the surrogated-model-based RBDO. In earlier researches, constraint boundary sampling (CBS) was proposed to approximate accurately the boundaries of constraints by locating sample points on the boundaries of constraints. However, because CBS uses sample points on all constraint boundaries, it creates superfluous sample points. In this paper, efficient constraint boundary sampling (ECBS) is proposed to enhance the efficiency of CBS. ECBS uses the statistical information of a kriging surrogate model to locate sample points on or near the RBDO solution. The efficiency of ECBS is verified by mathematical examples.

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

Converting boundary representation solid models to half-space representation models for Monte Carlo analysis

International Nuclear Information System (INIS)

Davis, J. E.; Eddy, M. J.; Sutton, T. M.; Altomari, T. J.

2007-01-01

Solid modeling computer software systems provide for the design of three-dimensional solid models used in the design and analysis of physical components. The current state-of-the-art in solid modeling representation uses a boundary representation format in which geometry and topology are used to form three-dimensional boundaries of the solid. The geometry representation used in these systems is cubic B-spline curves and surfaces - a network of cubic B-spline functions in three-dimensional Cartesian coordinate space. Many Monte Carlo codes, however, use a geometry representation in which geometry units are specified by intersections and unions of half-spaces. This paper describes an algorithm for converting from a boundary representation to a half-space representation. (authors)

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.

A "1"3"7Cs erosion model with moving boundary

International Nuclear Information System (INIS)

Yin, Chuan; Ji, Hongbing

2015-01-01

A novel quantitative model of the relationship between diffused concentration changes and erosion rates using assessment of soil losses was developed. It derived from the analysis of surface soil "1"3"7Cs flux variation under persistent erosion effect and based on the principle of geochemistry kinetics moving boundary. The new moving boundary model improves the basic simplified transport model (Zhang et al., 2008), and mainly applies to uniform rainfall areas which show a long-time soil erosion. The simulation results for this kind of erosion show under a long-time soil erosion, the influence of "1"3"7Cs concentration will decrease exponentially with increasing depth. Using the new model fit to the measured "1"3"7Cs depth distribution data in Zunyi site, Guizhou Province, China which has typical uniform rainfall provided a good fit with R"2 = 0.92. To compare the soil erosion rates calculated by the simple transport model and the new model, we take the Kaixian reference profile as example. The soil losses estimated by the previous simplified transport model are greater than those estimated by the new moving boundary model, which is consistent with our expectations. - Highlights: • The diffused moving boundary principle analysing "1"3"7Cs flux variation. • The new erosion model applies to uniform rainfall areas. • The erosion effect on "1"3"7Cs will decrease exponentially with increasing depth. • The new model provides two methods of calculating erosion rate.

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.

OPERATOR-RELATED FORMULATION OF THE EIGENVALUE PROBLEM FOR THE BOUNDARY PROBLEM OF ANALYSIS OF A THREE-DIMENSIONAL STRUCTURE WITH PIECEWISE-CONSTANT PHYSICAL AND GEOMETRICAL PARAMETERS ALONGSIDE THE BASIC DIRECTION WITHIN THE FRAMEWORK OF THE DISCRETE-CON

Directory of Open Access Journals (Sweden)

Akimov Pavel Alekseevich

2012-10-01

Full Text Available The proposed paper covers the operator-related formulation of the eigenvalue problem of analysis of a three-dimensional structure that has piecewise-constant physical and geometrical parameters alongside the so-called basic direction within the framework of a discrete-continual approach (a discrete-continual finite element method, a discrete-continual variation method. Generally, discrete-continual formulations represent contemporary mathematical models that become available for computer implementation. They make it possible for a researcher to consider the boundary effects whenever particular components of the solution represent rapidly varying functions. Another feature of discrete-continual methods is the absence of any limitations imposed on lengths of structures. The three-dimensional problem of elasticity is used as the design model of a structure. In accordance with the so-called method of extended domain, the domain in question is embordered by an extended one of an arbitrary shape. At the stage of numerical implementation, relative key features of discrete-continual methods include convenient mathematical formulas, effective computational patterns and algorithms, simple data processing, etc. The authors present their formulation of the problem in question for an isotropic medium with allowance for supports restrained by elastic elements while standard boundary conditions are also taken into consideration.

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.

A simplified treatment of the boundary conditions of the k- ε model in coarse-mesh CFD-type codes

International Nuclear Information System (INIS)

Analytis, G.Th.; Andreani, M.

1999-01-01

In coarse-mesh, CFD-type codes such as the containment analysis code GOTHIC, one of the options that can be used for modelling of turbulence is the k - ε model. However, in contrast to most other CFD codes which are designed to perform detailed CFD calculations with a large number of spatial meshes, codes such as GOTHIC are primarily aimed at simplified calculation of transients in large spaces (e.g., reactor containments), and generally use coarse meshes. The solution of the two parabolic equations for the k - ε model requires the definition of boundary conditions at physical boundaries and this, in turn, requires very small spatial meshes near these boundaries. Hence, while in codes like CFX this is done in a rigorous and consistent manner, codes like GOTHIC adopt an indirect and heuristic approach, due to the fact that the spatial meshes are usually large. This can have adverse consequences during the calculation of a transient and in this work, we shall give some examples of this and outline a method by which this problem can be avoided. (author)

A model problem for estimation of moving-film time relaxation at sudden change of boundary conditions

Science.gov (United States)

Smirnovsky, Alexander A.; Eliseeva, Viktoria O.

2018-05-01

The study of the film flow occurred under the influence of a gas slug flow is of definite interest in heat and mass transfer during the motion of a coolant in the second circuit of a nuclear water-water reactor. Thermohydraulic codes are usually used for analysis of the such problems in which the motion of the liquid film and the vapor is modeled on the basis of a one-dimensional balance equations. Due to a greater inertia of the liquid film motion, film flow parameters changes with a relaxation compared with gas flow. We consider a model problem of film flow under the influence of friction from gas slug flow neglecting such effects as wave formation, droplet breakage and deposition on the film surface, evaporation and condensation. Such a problem is analogous to the well-known problems of Couette and Stokes flows. An analytical solution has been obtained for laminar flow. Numerical RANS-based simulation of turbulent flow was performed using OpenFOAM. It is established that the relaxation process is almost self-similar. This fact opens a possibility of obtaining valuable correlations for the relaxation time.

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

Spline-Interpolation Solution of One Elasticity Theory Problem

CERN Document Server

Shirakova, Elena A

2011-01-01

The book presents methods of approximate solution of the basic problem of elasticity for special types of solids. Engineers can apply the approximate methods (Finite Element Method, Boundary Element Method) to solve the problems but the application of these methods may not be correct for solids with the certain singularities or asymmetrical boundary conditions. The book is recommended for researchers and professionals working on elasticity modeling. It explains methods of solving elasticity problems for special solids. Approximate methods (Finite Element Method, Boundary Element Method) have b

Symbolic boundary calculus and feedback operators

DEFF Research Database (Denmark)

Pedersen, Michael

1991-01-01

The recent developments in microlocal analysis and pseudodifferential boundary calculus are well suited tools in the investigation of a large number of problems occurring in control theory for partial differential equations. We explain some of the basic ideas of a pseudodifferential model...

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

Representing the atmospheric boundary layer in climate models of intermediate compexity

NARCIS (Netherlands)

Ronda, R.J.; Haarsma, R.J.; Holtslag, A.A.M.

2003-01-01

In this study the role of atmospheric boundary layer schemes in climate models is investigated. Including a boundary layer scheme in an Earth system model of intermediate complexity (EMIC) produces only minor differences in the estimated global distribution of sensible and latent heat fluxes over

Boundary condition histograms for modulated phases

International Nuclear Information System (INIS)

Benakli, M.; Gabay, M.; Saslow, W.M.

1997-11-01

Boundary conditions strongly affect the results of numerical computations for finite size inhomogeneous or incommensurate structures. We present a method which allows to deal with this problem, both for ground state and for critical properties: it combines fluctuating boundary conditions and specific histogram techniques. Our approach concerns classical as well as quantum systems. In particular, current-current correlation functions, which probe large scale coherence of the states, can be accurately evaluated. We illustrate our method on a frustrated two dimensional XY model. (author)

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.

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.

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

Directory of Open Access Journals (Sweden)