Dynamical Casimir effect on a cavity with mixed boundary conditions

International Nuclear Information System (INIS)

Alves, Danilo T.; Farina, Carlos; Maia Neto, Paulo Americo

2002-01-01

The most well-known mechanical effect related to the quantum vacuum is the Casimir force between two mirrors at rest. A new effect appears when the mirrors are set to move. In this case, the vacuum field may exert a dissipative force, damping the motion. As a consequence of energy conservation, there will be creation of real particles. If the motion is non-relativistic and has a small amplitude, the dynamical Casimir force can be found via a perturbative method proposed by Ford and Vilenkin. Using their technique, the electromagnetic dynamical Casimir problem, considered when the oscillating cavity is formed by two parallel plates of the same nature (perfectly conducting or perfectly permeable), can be divided into two separated boundary condition problems, namely: one involving Dirichlet BC, related to the transverse electric polarization and the other involving a Neumann BC, related to the transverse magnetic mode. The case of conducting plates can be found in the literature. However, another interesting case, the mixed oscillating cavity where the plates are of different nature, namely, a perfectly conducting plate and a perfectly permeable one (Boyer plates), has not been studied yet. We show that,for this case, the transverse electric models will be related to mixed boundary conditions: Dirichlet-like BC at the conducting plate and Neumann-like BC at the permeable plate. Analogously, the magnetic modes are related to a Neumann BC at the conducting plate and to a Dirichlet BC at the permeable one. As a first step before attacking the three-dimensional electromagnetic problem with mixed BC, we present here a simpler model: a one-dimensional cavity, where a massless scalar field is submitted to mixed (Dirichlet-Neumann) BC. For simplicity, we consider a non-relativistic motion for the conducting wall (Dirichlet BC) and suppose that the perfectly permeable wall (Neumann BC) is at rest. From this model we can extract insights about the dynamical Casimir

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.

Mixed global anomalies and boundary conformal field theories

OpenAIRE

Numasawa, Tokiro; Yamaguchi, Satoshi

2017-01-01

We consider the relation of mixed global gauge gravitational anomalies and boundary conformal field theory in WZW models for simple Lie groups. The discrete symmetries of consideration are the centers of the simple Lie groups. These mixed anomalies prevent to gauge them i.e, take the orbifold by the center. The absence of anomalies impose conditions on the levels of WZW models. Next, we study the conformal boundary conditions for the original theories. We consider the existence of a conformal...

Mixed basin boundary structures of chaotic systems

International Nuclear Information System (INIS)

Rosa, E. Jr.; Ott, E.

1999-01-01

Motivated by recent numerical observations on a four-dimensional continuous-time dynamical system, we consider different types of basin boundary structures for chaotic systems. These general structures are essentially mixtures of the previously known types of basin boundaries where the character of the boundary assumes features of the previously known boundary types at different points arbitrarily finely interspersed in the boundary. For example, we discuss situations where an everywhere continuous boundary that is otherwise smooth and differentiable at almost every point has an embedded uncountable, zero Lebesgue measure set of points at which the boundary curve is nondifferentiable. Although the nondifferentiable set is only of zero Lebesgue measure, the curve close-quote s fractal dimension may (depending on parameters) still be greater than one. In addition, we discuss bifurcations from such a mixed boundary to a 'pure' boundary that is a fractal nowhere differentiable curve or surface and to a pure nonfractal boundary that is everywhere smooth. copyright 1999 The American Physical Society

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.

Heat Kernel Asymptotics of Zaremba Boundary Value Problem

Energy Technology Data Exchange (ETDEWEB)

Avramidi, Ivan G. [Department of Mathematics, New Mexico Institute of Mining and Technology (United States)], E-mail: iavramid@nmt.edu

2004-03-15

The Zaremba boundary-value problem is a boundary value problem for Laplace-type second-order partial differential operators acting on smooth sections of a vector bundle over a smooth compact Riemannian manifold with smooth boundary but with discontinuous boundary conditions, which include Dirichlet boundary conditions on one part of the boundary and Neumann boundary conditions on another part of the boundary. We study the heat kernel asymptotics of Zaremba boundary value problem. The construction of the asymptotic solution of the heat equation is described in detail and the heat kernel is computed explicitly in the leading approximation. Some of the first nontrivial coefficients of the heat kernel asymptotic expansion are computed explicitly.

On a Mixed Nonlinear One Point Boundary Value Problem for an Integrodifferential Equation

Directory of Open Access Journals (Sweden)

Mesloub Said

2008-01-01

Full Text Available This paper is devoted to the study of a mixed problem for a nonlinear parabolic integro-differential equation which mainly arise from a one dimensional quasistatic contact problem. We prove the existence and uniqueness of solutions in a weighted Sobolev space. Proofs are based on some a priori estimates and on the Schauder fixed point theorem. we also give a result which helps to establish the regularity of a solution.

Rethinking workforce boundaries: roles, responsibilities and skill mix and readiness for change in general practice

OpenAIRE

Topping, Annie; Nkosana Nyawata, Idah; Stephenson, John; Featherstone, Valerie A.

2012-01-01

Title: Rethinking workforce boundaries: roles, responsibilities and skill mix and readiness for change in general practice \\ud The Problem \\ud The last 10 years has seen major changes in the way services are delivered in primary care. Skill mix, has offered many practices real opportunities for doing things differently. As the introduction of advanced nurse practitioners (ANPs) and health care assistants (HCAs) into the primary care workforce demonstrate. While workforce redesign has its crit...

Boundary value problems and partial differential equations

CERN Document Server

Powers, David L

2005-01-01

Boundary Value Problems is the leading text on boundary value problems and Fourier series. The author, David Powers, (Clarkson) has written a thorough, theoretical overview of solving boundary value problems involving partial differential equations by the methods of separation of variables. Professors and students agree that the author is a master at creating linear problems that adroitly illustrate the techniques of separation of variables used to solve science and engineering.* CD with animations and graphics of solutions, additional exercises and chapter review questions* Nearly 900 exercises ranging in difficulty* Many fully worked examples

Topology optimization of acoustic-structure interaction problems using a mixed finite element formulation

DEFF Research Database (Denmark)

Yoon, Gil Ho; Jensen, Jens Stissing; Sigmund, Ole

2007-01-01

given during the optimization process. In this paper we circumvent the explicit boundary representation by using a mixed finite element formulation with displacements and pressure as primary variables (a u/p-formulation). The Helmholtz equation is obtained as a special case of the mixed formulation...... for the elastic shear modulus equating to zero. Hence, by spatial variation of the mass density, shear and bulk moduli we are able to solve the coupled problem by the mixed formulation. Using this modelling approach, the topology optimization procedure is simply implemented as a standard density approach. Several...... two-dimensional acoustic-structure problems are optimized in order to verify the proposed method....

Mixed FEM for Second Order Elliptic Problems on Polygonal Meshes with BEM-Based Spaces

KAUST Repository

Efendiev, Yalchin

2014-01-01

We present a Boundary Element Method (BEM)-based FEM for mixed formulations of second order elliptic problems in two dimensions. The challenge, we would like to address, is a proper construction of H(div)-conforming vector valued trial functions on arbitrary polygonal partitions of the domain. The proposed construction generates trial functions on polygonal elements which inherit some of the properties of the unknown solution. In the numerical realization, the relevant local problems are treated by means of boundary integral formulations. We test the accuracy of the method on two model problems. © 2014 Springer-Verlag.

The determination of the mixing height. Current progress and problems

Energy Technology Data Exchange (ETDEWEB)

Gryning, S.E.; Beyrich, F.; Batchvarova, E. [eds.

1997-10-01

This report contains extended abstracts of presentations given at a EURASAP Workshop on The Determination of the Mixing Height - Current Progress and Problems. The Workshop, initiated from discussions with Peter Builtjes, was held at Risoe National Laboratory 1-3 October 1997 within the framework of EURASAP (European Association for the Sciences of Air Pollution). The specific topics and chairpersons of the Workshop were: Theoretical Considerations (Sven-Erik Gryning), Mixing Height Estimation from Turbulence Measurements and In-Situ Soundings (Douw Steyn), Mixing Height Determination from NWP-Models (Han van Dop), Climatology and Global Aspects (Werner Klug), Mixing Height Determination from Remote Systems (Werner Klug), Verification of Mixing Height Parameterizations and Models (Frank Beyrich), Mixing Height over Complex Terrain (Ekaterina Batchvarova), Internal Boundary Layers: Mixing Height in Coastal Areas and Over Cities (Allen White). The discussion at the end of the Workshop was chaired by Robert Bornstein. (au)

State-dependent impulses boundary value problems on compact interval

CERN Document Server

Rachůnková, Irena

2015-01-01

This book offers the reader a new approach to the solvability of boundary value problems with state-dependent impulses and provides recently obtained existence results for state dependent impulsive problems with general linear boundary conditions. It covers fixed-time impulsive boundary value problems both regular and singular and deals with higher order differential equations or with systems that are subject to general linear boundary conditions. We treat state-dependent impulsive boundary value problems, including a new approach giving effective conditions for the solvability of the Dirichlet problem with one state-dependent impulse condition and we show that the depicted approach can be extended to problems with a finite number of state-dependent impulses. We investigate the Sturm–Liouville boundary value problem for a more general right-hand side of a differential equation. Finally, we offer generalizations to higher order differential equations or differential systems subject to general linear boundary...

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.

Numerical stability of finite difference algorithms for electrochemical kinetic simulations: Matrix stability analysis of the classic explicit, fully implicit and Crank-Nicolson methods and typical problems involving mixed boundary conditions

DEFF Research Database (Denmark)

Bieniasz, Leslaw K.; Østerby, Ole; Britz, Dieter

1995-01-01

The stepwise numerical stability of the classic explicit, fully implicit and Crank-Nicolson finite difference discretizations of example diffusional initial boundary value problems from electrochemical kinetics has been investigated using the matrix method of stability analysis. Special attention...... has been paid to the effect of the discretization of the mixed, linear boundary condition with time-dependent coefficients on stability, assuming the two-point forward-difference approximations for the gradient at the left boundary (electrode). Under accepted assumptions one obtains the usual...... stability criteria for the classic explicit and fully implicit methods. The Crank-Nicolson method turns out to be only conditionally stable in contrast to the current thought regarding this method....

A general mixed boundary model reduction method for component mode synthesis

NARCIS (Netherlands)

Voormeeren, S.N.; Van der Valk, P.L.C.; Rixen, D.J.

2010-01-01

A classic issue in component mode synthesis (CMS) methods is the choice for fixed or free boundary conditions at the interface degrees of freedom (DoF) and the associated vibration modes in the components reduction base. In this paper, a novel mixed boundary CMS method called the “Mixed

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.

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)

Mixed FEM for Second Order Elliptic Problems on Polygonal Meshes with BEM-Based Spaces

KAUST Repository

Efendiev, Yalchin; Galvis, Juan; Lazarov, Raytcho; Weiß er, Steffen

2014-01-01

We present a Boundary Element Method (BEM)-based FEM for mixed formulations of second order elliptic problems in two dimensions. The challenge, we would like to address, is a proper construction of H(div)-conforming vector valued trial functions

Analytical solution for the convectively-mixed atmospheric boundary layer

NARCIS (Netherlands)

Ouwersloot, H.G.; Vilà-Guerau de Arellano, J.

2013-01-01

Based on the prognostic equations of mixed-layer theory assuming a zeroth order jump at the entrainment zone, analytical solutions for the boundary-layer height evolution are derived with different degrees of accuracy. First, an exact implicit expression for the boundary-layer height for a situation

Modified Taylor series method for solving nonlinear differential equations with mixed boundary conditions defined on finite intervals.

Science.gov (United States)

Vazquez-Leal, Hector; Benhammouda, Brahim; Filobello-Nino, Uriel Antonio; Sarmiento-Reyes, Arturo; Jimenez-Fernandez, Victor Manuel; Marin-Hernandez, Antonio; Herrera-May, Agustin Leobardo; Diaz-Sanchez, Alejandro; Huerta-Chua, Jesus

2014-01-01

In this article, we propose the application of a modified Taylor series method (MTSM) for the approximation of nonlinear problems described on finite intervals. The issue of Taylor series method with mixed boundary conditions is circumvented using shooting constants and extra derivatives of the problem. In order to show the benefits of this proposal, three different kinds of problems are solved: three-point boundary valued problem (BVP) of third-order with a hyperbolic sine nonlinearity, two-point BVP for a second-order nonlinear differential equation with an exponential nonlinearity, and a two-point BVP for a third-order nonlinear differential equation with a radical nonlinearity. The result shows that the MTSM method is capable to generate easily computable and highly accurate approximations for nonlinear equations. 34L30.

Numerical Methods for Free Boundary Problems

CERN Document Server

1991-01-01

About 80 participants from 16 countries attended the Conference on Numerical Methods for Free Boundary Problems, held at the University of Jyviiskylii, Finland, July 23-27, 1990. The main purpose of this conference was to provide up-to-date information on important directions of research in the field of free boundary problems and their numerical solutions. The contributions contained in this volume cover the lectures given in the conference. The invited lectures were given by H.W. Alt, V. Barbu, K-H. Hoffmann, H. Mittelmann and V. Rivkind. In his lecture H.W. Alt considered a mathematical model and existence theory for non-isothermal phase separations in binary systems. The lecture of V. Barbu was on the approximate solvability of the inverse one phase Stefan problem. K-H. Hoff­ mann gave an up-to-date survey of several directions in free boundary problems and listed several applications, but the material of his lecture is not included in this proceedings. H.D. Mittelmann handled the stability of thermo capi...

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

A general mixed boundary model reduction method for component mode synthesis

International Nuclear Information System (INIS)

Voormeeren, S N; Van der Valk, P L C; Rixen, D J

2010-01-01

A classic issue in component mode synthesis (CMS) methods is the choice for fixed or free boundary conditions at the interface degrees of freedom (DoF) and the associated vibration modes in the components reduction base. In this paper, a novel mixed boundary CMS method called the 'Mixed Craig-Bampton' method is proposed. The method is derived by dividing the substructure DoF into a set of internal DoF, free interface DoF and fixed interface DoF. To this end a simple but effective scheme is introduced that, for every pair of interface DoF, selects a free or fixed boundary condition for each DoF individually. Based on this selection a reduction basis is computed consisting of vibration modes, static constraint modes and static residual flexibility modes. In order to assemble the reduced substructures a novel mixed assembly procedure is developed. It is shown that this approach leads to relatively sparse reduced matrices, whereas other mixed boundary methods often lead to full matrices. As such, the Mixed Craig-Bampton method forms a natural generalization of the classic Craig-Bampton and more recent Dual Craig-Bampton methods. Finally, the method is applied to a finite element test model. Analysis reveals that the proposed method has comparable or better accuracy and superior versatility with respect to the existing methods.

Free Convection over a Permeable Horizontal Flat Plate Embedded in a Porous Medium with Radiation Effects and Mixed Thermal Boundary Conditions

OpenAIRE

Najiyah S. Khasi'ie; Roziena Khairuddin; Najihah Mohamed; Mohd Zuki Salleh; Roslinda Nazar; Ioan Pop

2012-01-01

Problem statement: In this study, the mathematical modeling of free convection boundary layer flow over a permeable horizontal flat plate embedded in a porous medium under mixed thermal boundary conditions and radiation effects is considered. Approach: The transformed boundary layer equations are solved numerically using the shooting method. Results: Numerical solutions are obtained for the wall temperature, the heat transfer coefficient, as well as the velocity and temperature profiles. The ...

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.

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

Stochastic Theory of Turbulence Mixing by Finite Eddies in the Turbulent Boundary Layer

NARCIS (Netherlands)

Dekker, H.; Leeuw, G. de; Maassen van den Brink, A.

1995-01-01

Turbulence mixing is treated by means of a novel formulation of nonlocal K-theory, involving sample paths and a stochastic hypothesis. The theory simplifies for mixing by exchange (strong-eddies) and is then applied to the boundary layer (involving scaling). This maps boundary layer turbulence onto

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.

Thin-film superconducting rings in the critical state: the mixed boundary value approach

Science.gov (United States)

Brambilla, Roberto; Grilli, Francesco

2015-02-01

In this paper, we describe the critical state of a thin superconducting ring (and of a perfectly conducting ring as a limiting case) as a mixed boundary value problem. The disc is characterized by a three-part boundary division of the positive real axis, so this work is an extension of the procedure used in a previous work of ours for describing superconducting discs and strips, which are characterized by a two-part boundary division of the real axis. Here, we present the mathematical tools to solve this kind of problems—the Erdélyi-Kober operators—in a frame that can be immediately used. Contrary to the two-part problems considered in our previous work, three-part problems do not generally have analytical solutions and the numerical work takes on a significant heaviness. Nevertheless, this work is remunerated by three clear advantages: firstly, all the cases are afforded in the same way, without the necessity of any brilliant invention or ability; secondly, in the case of superconducting rings, the penetration of the magnetic field in the internal/external rims is a result of the method itself and does not have to be imposed, as it is commonly done with other methods presented in the literature; thirdly, the method can be extended to investigate even more complex cases (four-part problems). In this paper, we consider the cases of rings in uniform field and with transport current, with or without flux trapping in the hole and the case without net current, corresponding to a cut ring (washer), as used in some SQUID applications.

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

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

An A Posteriori Error Analysis of Mixed Finite Element Galerkin Approximations to Second Order Linear Parabolic Problems

KAUST Repository

Memon, Sajid; Nataraj, Neela; Pani, Amiya Kumar

2012-01-01

In this article, a posteriori error estimates are derived for mixed finite element Galerkin approximations to second order linear parabolic initial and boundary value problems. Using mixed elliptic reconstructions, a posteriori error estimates in L∞(L2)- and L2(L2)-norms for the solution as well as its flux are proved for the semidiscrete scheme. Finally, based on a backward Euler method, a completely discrete scheme is analyzed and a posteriori error bounds are derived, which improves upon earlier results on a posteriori estimates of mixed finite element approximations to parabolic problems. Results of numerical experiments verifying the efficiency of the estimators have also been provided. © 2012 Society for Industrial and Applied Mathematics.

Solution of moving boundary problems with implicit boundary condition

International Nuclear Information System (INIS)

Moyano, E.A.

1990-01-01

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

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

OpenAIRE

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

2014-01-01

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

Mixed Element Formulation for the Finite Element-Boundary Integral Method

National Research Council Canada - National Science Library

Meese, J; Kempel, L. C; Schneider, S. W

2006-01-01

A mixed element approach using right hexahedral elements and right prism elements for the finite element-boundary integral method is presented and discussed for the study of planar cavity-backed antennas...

Casimir energy in d-dimensional rectangular geometries, under mixed boundary conditions

International Nuclear Information System (INIS)

Silva, J.C. da; Placido, Hebe Q.; Santana, A.E.; M Neto, Arthur

1997-01-01

The Casimir energy and its temperature corrections are presented for the electromagnetic field confined in a d-dimensional hypercavity. The expressions are derived considering Dirichlet boundary conditions for each pair of hyperplanes defining a confined direction (the homogeneous case); or yet, by choosing different boundary conditions (Dirichlet or Neumann) at each hyperplane of the pair (the mixed case). (author)

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

Science.gov (United States)

Reimer, Ashton S.; Cheviakov, Alexei F.

2013-03-01

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

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.

Numerical simulation of small-scale mixing processes in the upper ocean and atmospheric boundary layer

International Nuclear Information System (INIS)

Druzhinin, O; Troitskaya, Yu; Zilitinkevich, S

2016-01-01

The processes of turbulent mixing and momentum and heat exchange occur in the upper ocean at depths up to several dozens of meters and in the atmospheric boundary layer within interval of millimeters to dozens of meters and can not be resolved by known large- scale climate models. Thus small-scale processes need to be parameterized with respect to large scale fields. This parameterization involves the so-called bulk coefficients which relate turbulent fluxes with large-scale fields gradients. The bulk coefficients are dependent on the properties of the small-scale mixing processes which are affected by the upper-ocean stratification and characteristics of surface and internal waves. These dependencies are not well understood at present and need to be clarified. We employ Direct Numerical Simulation (DNS) as a research tool which resolves all relevant flow scales and does not require closure assumptions typical of Large-Eddy and Reynolds Averaged Navier-Stokes simulations (LES and RANS). Thus DNS provides a solid ground for correct parameterization of small-scale mixing processes and also can be used for improving LES and RANS closure models. In particular, we discuss the problems of the interaction between small-scale turbulence and internal gravity waves propagating in the pycnocline in the upper ocean as well as the impact of surface waves on the properties of atmospheric boundary layer over wavy water surface. (paper)

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

Directory of Open Access Journals (Sweden)

Raúl B. González De Paz

2009-02-01

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

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

Directory of Open Access Journals (Sweden)

Jufang Wang

2013-01-01

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

Problems of matter-antimatter boundary layers

International Nuclear Information System (INIS)

Lehnert, B.

1975-01-01

This paper outlines the problems of the quasi-steady matter-antimatter boundary layers discussed in Klein-Alfven's cosmological theory, and a crude model of the corresponding ambiplasma balance is presented: (i) at interstellar particle densities, no well-defined boundary layer can exist in presence of neutral gas, nor can such a layer be sustained in an unmagnetized fully ionized ambiplasma. (ii) Within the limits of applicability of the present model, sharply defined boundary layers are under certain conditions found to exist in a magnetized ambiplasma. Thus, at beta values less than unity, a steep pressure drop of the low-energy components of matter and antimatter can be balanced by a magnetic field and the electric currents in the ambiplasma. (iii) The boundary layer thickness is of the order of 2x 0 approximately 10/BT 0 sup(1/4) meters, where B is the magnetic field strength in MKS units and T 0 the characteristic temperature of the low-energy components in the layer. (Auth.)

Spectral results for mixed problems and fractional elliptic operators,

DEFF Research Database (Denmark)

Grubb, Gerd

2015-01-01

In the first part of the paper we show Weyl type spectral asymptotic formulas for pseudodifferential operators P a of order 2a, with type and factorization index a  ∈ R +, restricted to compact sets with boundary; this includes fractional powers of the Laplace operator. The domain...... and the regularity of eigenfunctions is described. In the second part, we apply this in a study of realizations A χ,Σ+ in L 2( Ω ) of mixed problems for a second-order strongly elliptic symmetric differential operator A on a bounded smooth set Ω ⊂ R n; here the boundary ∂Ω=Σ is partioned smoothly into Σ......=Σ_∪Σ+, the Dirichlet condition γ0u=0 is imposed on Σ_, and a Neumann or Robin condition χu=0 is imposed on Σ+. It is shown that the Dirichlet-to-Neumann operator Pγ,χ is principally of type 1/2 with factorization index 1/2, relative to Σ+. The above theory allows a detailed description of D (Aχ,Σ_+) with singular...

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.

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.

Weak Thermocline Mixing in the North Pacific Low-Latitude Western Boundary Current System

Science.gov (United States)

Liu, Zhiyu; Lian, Qiang; Zhang, Fangtao; Wang, Lei; Li, Mingming; Bai, Xiaolin; Wang, Jianing; Wang, Fan

2017-10-01

Despite its potential importance in the global climate system, mixing properties of the North Pacific low-latitude western boundary current system (LLWBC) remained unsampled until very recently. We report here on the first measurements of turbulence microstructure associated with these currents, made in the western boundary region of the tropical North Pacific east of the Philippines. The results suggest that thermocline mixing in the North Pacific LLWBC is generally weak with the diapycnal diffusivity κρ˜O(10-6) m2 s-1. This is consistent with predictions from internal wave-wave interaction theory that mixing due to internal wave breaking is significantly reduced at low latitudes. Enhanced mixing is found to be associated with a permanent cyclonic eddy, the Mindanao Eddy, but mainly at its south and north flanks. There, κρ is elevated by an order of magnitude due to eddy-induced geostrophic shear. Mixing in the eddy core is at the background level with no indication of enhancement.

Performance of mixed formulations for the particle finite element method in soil mechanics problems

Science.gov (United States)

Monforte, Lluís; Carbonell, Josep Maria; Arroyo, Marcos; Gens, Antonio

2017-07-01

This paper presents a computational framework for the numerical analysis of fluid-saturated porous media at large strains. The proposal relies, on one hand, on the particle finite element method (PFEM), known for its capability to tackle large deformations and rapid changing boundaries, and, on the other hand, on constitutive descriptions well established in current geotechnical analyses (Darcy's law; Modified Cam Clay; Houlsby hyperelasticity). An important feature of this kind of problem is that incompressibility may arise either from undrained conditions or as a consequence of material behaviour; incompressibility may lead to volumetric locking of the low-order elements that are typically used in PFEM. In this work, two different three-field mixed formulations for the coupled hydromechanical problem are presented, in which either the effective pressure or the Jacobian are considered as nodal variables, in addition to the solid skeleton displacement and water pressure. Additionally, several mixed formulations are described for the simplified single-phase problem due to its formal similitude to the poromechanical case and its relevance in geotechnics, since it may approximate the saturated soil behaviour under undrained conditions. In order to use equal-order interpolants in displacements and scalar fields, stabilization techniques are used in the mass conservation equation of the biphasic medium and in the rest of scalar equations. Finally, all mixed formulations are assessed in some benchmark problems and their performances are compared. It is found that mixed formulations that have the Jacobian as a nodal variable perform better.

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.

Free boundary problems in PDEs and particle systems

CERN Document Server

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

2016-01-01

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

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.

Regularization of moving boundaries in a Laplacian field by a mixed Dirichlet-Neumann boundary condition : exact results

NARCIS (Netherlands)

B.J. Meulenbroek (Bernard); U. M. Ebert (Ute); L. Schäfer

2005-01-01

textabstractThe dynamics of ionization fronts that generate a conducting body, are in simplest approximation equivalent to viscous fingering without regularization. Going beyond this approximation, we suggest that ionization fronts can be modeled by a mixed Dirichlet-Neumann boundary condition. We

Heading off boundary problems: clinical supervision as risk management.

Science.gov (United States)

Walker, R; Clark, J J

1999-11-01

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

Separable boundary-value problems in physics

CERN Document Server

Willatzen, Morten

2011-01-01

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

Regularization of moving boundaries in a Laplacian field by a mixed dirichlet-neumann boundary condition: Exact results

NARCIS (Netherlands)

Meulenbroek, B.; Ebert, U.; Schäfer, L.

2005-01-01

The dynamics of ionization fronts that generate a conducting body are in the simplest approximation equivalent to viscous fingering without regularization. Going beyond this approximation, we suggest that ionization fronts can be modeled by a mixed Dirichlet-Neumann boundary condition. We derive

Analytic Solution to Shell Boundary – Value Problems

Directory of Open Access Journals (Sweden)

Yu. I. Vinogradov

2015-01-01

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

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

Science.gov (United States)

Zhang, Yabin; Gillman, Adrianna

2018-03-01

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

A finite difference method for free boundary problems

KAUST Repository

Fornberg, Bengt

2010-04-01

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

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

International Nuclear Information System (INIS)

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

1999-01-01

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

Group invariance in engineering boundary value problems

CERN Document Server

Seshadri, R

1985-01-01

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

A 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

Boundary element methods applied to two-dimensional neutron diffusion problems

International Nuclear Information System (INIS)

Itagaki, Masafumi

1985-01-01

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

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

OpenAIRE

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

2016-01-01

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

Boundary value problems of holomorphic vector functions in 1D QCs

International Nuclear Information System (INIS)

Gao Yang; Zhao Yingtao; Zhao Baosheng

2007-01-01

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

Laplace boundary-value problem in paraboloidal coordinates

International Nuclear Information System (INIS)

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

2012-01-01

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

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

Directory of Open Access Journals (Sweden)

Ruzanna Kh. Makaova

2017-12-01

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

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

NARCIS (Netherlands)

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

2009-01-01

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

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

Science.gov (United States)

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

2017-10-01

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

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

International Nuclear Information System (INIS)

Halbach, K.; Schlueter, R.

1992-01-01

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

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

Modal and Mixed Specifications: Key Decision Problems and their Complexities

DEFF Research Database (Denmark)

Antonik, Adam; Huth, Michael; Larsen, Kim Guldstrand

2010-01-01

, and whether all implementations of one specification are implementations of another one. For each of these decision problems we investigate the worst-case computational complexity for the modal and mixed case. We show that the first decision problem is EXPTIME-complete for modal as well as for mixed......Modal and mixed transition systems are specification formalisms that allow mixing of over- and under-approximation. We discuss three fundamental decision problems for such specifications: whether a set of specifications has a common implementation, whether a sole specification has an implementation...... specifications. We prove that the second decision problem is EXPTIME-complete for mixed specifications (while it is known to be trivial for modal ones). The third decision problem is furthermore demonstrated to be EXPTIME-complete for mixed specifications....

The boundary value problems of magnetotail equilibrium

International Nuclear Information System (INIS)

Birn, J.

1991-01-01

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

Positive solutions and eigenvalues of nonlocal boundary-value problems

Directory of Open Access Journals (Sweden)

Jifeng Chu

2005-07-01

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

On a non-linear pseudodifferential boundary value problem

International Nuclear Information System (INIS)

Nguyen Minh Chuong.

1989-12-01

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

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

CERN Document Server

Kitahara, M

1985-01-01

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

Eigenvalue inequalities for the Laplacian with mixed boundary conditions

Czech Academy of Sciences Publication Activity Database

Lotoreichik, Vladimir; Rohleder, J.

2017-01-01

Roč. 263, č. 1 (2017), s. 491-508 ISSN 0022-0396 R&D Projects: GA ČR(CZ) GA14-06818S Institutional support: RVO:61389005 Keywords : Laplace operator * mixed boundary conditions * eigenvalue inequality * polyhedral domain * Lipschitz domain Subject RIV: BE - Theoretical Physics OBOR OECD: Atomic, molecular and chemical physics (physics of atoms and molecules including collision, interaction with radiation, magnetic resonances, Mössbauer effect) Impact factor: 1.988, year: 2016

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

Science.gov (United States)

Sayevand, K.; Pichaghchi, K.

2018-04-01

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

Zero-energy eigenstates for the Dirac boundary problem

International Nuclear Information System (INIS)

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

1980-01-01

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

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

CERN Document Server

Eom, Hyo J

2004-01-01

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

Existence results for anisotropic discrete boundary value problems

Directory of Open Access Journals (Sweden)

Avci Avci

2016-06-01

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

Homology in Electromagnetic Boundary Value Problems

Directory of Open Access Journals (Sweden)

Pellikka Matti

2010-01-01

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

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

Boundary conditions for the use of personal ventilation over mixing ventilation in open plan offices

DEFF Research Database (Denmark)

Petersen, Steffen; Hviid, Christian Anker

2013-01-01

This paper investigates the boundary conditions for choosing a combined Personal Ventilation (PV) and Mixing Ventilation (MV) over conventional mixing ventilation in an office with multiple workers. A simplified procedure for annual performance assessment of PV/MV systems in terms of air quality...

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

Science.gov (United States)

Xia, Baoqiang; Fokas, A. S.

2018-02-01

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

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

International Nuclear Information System (INIS)

Shen, S

2008-01-01

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

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

OpenAIRE

Dipierro, Serena; Karakhanyan, Aram; Valdinoci, Enrico

2017-01-01

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

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

Energy Technology Data Exchange (ETDEWEB)

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

2017-06-01

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

A Duality Approach for the Boundary Variation of Neumann Problems

DEFF Research Database (Denmark)

Bucur, Dorin; Varchon, Nicolas

2002-01-01

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

A duality approach or the boundary variation of Neumann problems

DEFF Research Database (Denmark)

Bucur, D.; Varchon, Nicolas

2002-01-01

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

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

KAUST Repository

Dujardin, G. M.

2009-01-01

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

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

Energy Technology Data Exchange (ETDEWEB)

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

2013-07-31

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

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

Directory of Open Access Journals (Sweden)

Baoqiang Yan

2015-01-01

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

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

Science.gov (United States)

1991-03-15

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

Fourth-order discrete anisotropic boundary-value problems

Directory of Open Access Journals (Sweden)

Maciej Leszczynski

2015-09-01

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

Mixed-hybrid finite element method for the transport equation and diffusion approximation of transport problems

International Nuclear Information System (INIS)

Cartier, J.

2006-04-01

This thesis focuses on mathematical analysis, numerical resolution and modelling of the transport equations. First of all, we deal with numerical approximation of the solution of the transport equations by using a mixed-hybrid scheme. We derive and study a mixed formulation of the transport equation, then we analyse the related variational problem and present the discretization and the main properties of the scheme. We particularly pay attention to the behavior of the scheme and we show its efficiency in the diffusion limit (when the mean free path is small in comparison with the characteristic length of the physical domain). We present academical benchmarks in order to compare our scheme with other methods in many physical configurations and validate our method on analytical test cases. Unstructured and very distorted meshes are used to validate our scheme. The second part of this thesis deals with two transport problems. The first one is devoted to the study of diffusion due to boundary conditions in a transport problem between two plane plates. The second one consists in modelling and simulating radiative transfer phenomenon in case of the industrial context of inertial confinement fusion. (author)

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.

Mixed convection in inclined lid driven cavity by Lattice Boltzmann Method and heat flux boundary condition

International Nuclear Information System (INIS)

D'Orazio, A; Karimipour, A; Nezhad, A H; Shirani, E

2014-01-01

Laminar mixed convective heat transfer in two-dimensional rectangular inclined driven cavity is studied numerically by means of a double population thermal Lattice Boltzmann method. Through the top moving lid the heat flux enters the cavity whereas it leaves the system through the bottom wall; side walls are adiabatic. The counter-slip internal energy density boundary condition, able to simulate an imposed non zero heat flux at the wall, is applied, in order to demonstrate that it can be effectively used to simulate heat transfer phenomena also in case of moving walls. Results are analyzed over a range of the Richardson numbers and tilting angles of the enclosure, encompassing the dominating forced convection, mixed convection, and dominating natural convection flow regimes. As expected, heat transfer rate increases as increases the inclination angle, but this effect is significant for higher Richardson numbers, when buoyancy forces dominate the problem; for horizontal cavity, average Nusselt number decreases with the increase of Richardson number because of the stratified field configuration

Bivium as a Mixed Integer Programming Problem

DEFF Research Database (Denmark)

Borghoff, Julia; Knudsen, Lars Ramkilde; Stolpe, Mathias

2009-01-01

over $GF(2)$ into a combinatorial optimization problem. We convert the Boolean equation system into an equation system over $\\mathbb{R}$ and formulate the problem of finding a $0$-$1$-valued solution for the system as a mixed-integer programming problem. This enables us to make use of several...

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

Directory of Open Access Journals (Sweden)

Mustafa Turkyilmazoglu

2016-02-01

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

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

International Nuclear Information System (INIS)

Zuo Yinghong; Wang Jianguo

2012-01-01

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

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

Directory of Open Access Journals (Sweden)

Olivier Sarbach

2012-08-01

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

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

Science.gov (United States)

Sarbach, Olivier; Tiglio, Manuel

2012-01-01

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

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

Science.gov (United States)

Lau, Chun Sing

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

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

Institute of Scientific and Technical Information of China (English)

无

2012-01-01

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

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

Science.gov (United States)

Pakdemirli, Mehmet

2016-01-01

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

APPLICATION OF BOUNDARY INTEGRAL EQUATION METHOD FOR THERMOELASTICITY PROBLEMS

Directory of Open Access Journals (Sweden)

Vorona Yu.V.

2015-12-01

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

three solutions for a semilinear elliptic boundary value problem

Indian Academy of Sciences (India)

69

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

Parametrices and exact paralinearization of semi-linear boundary problems

DEFF Research Database (Denmark)

Johnsen, Jon

2008-01-01

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

Conforming discretizations of boundary element solutions to the electroencephalography forward problem

Science.gov (United States)

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

2018-01-01

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

Influence of convective conditions on three dimensional mixed convective hydromagnetic boundary layer flow of Casson nanofluid

Energy Technology Data Exchange (ETDEWEB)

Rauf, A., E-mail: raufamar@ciitsahiwal.edu.pk [Department of Mathematics, Comsats Institute of Information Technology, Sahiwal 57000 (Pakistan); Siddiq, M.K. [Centre for Advanced Studies in Pure and Applied Mathematics, Department of Mathematics, Bahauddin Zakariya University, Multan 63000 (Pakistan); Abbasi, F.M. [Department of Mathematics, Comsats Institute of Information Technology, Islamabad 44000 (Pakistan); Meraj, M.A. [Department of Mathematics, Comsats Institute of Information Technology, Sahiwal 57000 (Pakistan); Ashraf, M. [Centre for Advanced Studies in Pure and Applied Mathematics, Department of Mathematics, Bahauddin Zakariya University, Multan 63000 (Pakistan); Shehzad, S.A. [Department of Mathematics, Comsats Institute of Information Technology, Sahiwal 57000 (Pakistan)

2016-10-15

The present work deals with the steady laminar three-dimensional mixed convective magnetohydrodynamic (MHD) boundary layer flow of Casson nanofluid over a bidirectional stretching surface. A uniform magnetic field is applied normal to the flow direction. Similarity variables are implemented to convert the non-linear partial differential equations into ordinary ones. Convective boundary conditions are utilized at surface of the sheet. A numerical technique of Runge–Kutta–Fehlberg (RFK45) is used to obtain the results of velocity, temperature and concentration fields. The physical dimensionless parameters are discussed through tables and graphs. - Highlights: • Mixed convective boundary layer flow of Casson nanofluid is taken into account. • Impact of magnetic field is examined. • Convective heat and mass conditions are imposed. • Numerical solutions are presented and discussed.

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

Indian Academy of Sciences (India)

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

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

KAUST Repository

Leise, Tanya

2008-11-01

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

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

International Nuclear Information System (INIS)

Krieg, R.; Hailfinger, G.

1979-01-01

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

Accounting rigid support at the border in a mixed model the finite element method in problems of ice cover destruction

Directory of Open Access Journals (Sweden)

V. V. Knyazkov

2014-01-01

Full Text Available To evaluate the force to damage the ice covers is necessary for estimation of icebreaking capability of vessels, as well as of hull strength of icebreakers, and navigation of ships in ice conditions. On the other hand, the use of ice cover support to arrange construction works from the ice is also of practical interest.By the present moment a great deal of investigations of ice cover deformation have been carried out to result, usually, in approximate calculations formula which was obtained after making a variety of assumptions. Nevertheless, we believe that it is possible to make further improvement in calculations. Application numerical methods, and, for example, FEM, makes possible to avoid numerous drawbacks of analytical methods dealing with both complex boundaries and load application areas and other problem peculiarities.The article considers an application of mixed models of FEM for investigating ice cover deformation. A simple flexible triangle element of mixed type was taken to solve this problem. Vector of generalized coordinates of the element contains apices flexures and normal bending moments in the middle of its sides. Compared to other elements mixed models easily satisfy compatibility requirements on the boundary of adjacent elements and do not require numerical displacement differentiation to define bending moments, because bending moments are included in vector of element generalized coordinates.The method of account of rigid support plate is proposed. The resulting ratio, taking into account the "stiffening", reduces the number of resolving systems of equations by the number of elements on the plate contour.To evaluate further the results the numerical realization of ice cover stress-strained problem it becomes necessary and correct to check whether calculation results correspond to accurate solution. Using an example of circular plate the convergence of numerical solutions to analytical solutions is showed.The article

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

A three-point Taylor algorithm for three-point boundary value problems

NARCIS (Netherlands)

J.L. López; E. Pérez Sinusía; N.M. Temme (Nico)

2011-01-01

textabstractWe consider second-order linear differential equations $\\varphi(x)y''+f(x)y'+g(x)y=h(x)$ in the interval $(-1,1)$ with Dirichlet, Neumann or mixed Dirichlet-Neumann boundary conditions given at three points of the interval: the two extreme points $x=\\pm 1$ and an interior point

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

Science.gov (United States)

Gupta, Yogesh

2017-01-01

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

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

African Journals Online (AJOL)

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

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

International Nuclear Information System (INIS)

Pereira, Luis Carlos Martins

1998-06-01

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

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

Directory of Open Access Journals (Sweden)

E. V. Dikareva

2015-01-01

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

Boundary conditions for free surface inlet and outlet problems

KAUST Repository

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

2012-01-01

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

A conceptual framework to quantify the influence of convective boundary layer development on carbon dioxide mixing ratios

NARCIS (Netherlands)

Pino, D.; Vilà-Guerau de Arellano, J.; Peters, W.; Schröter, J.; van Heerwaarden, C. C.; Krol, M. C.

2012-01-01

Interpretation of observed diurnal carbon dioxide (CO2) mixing ratios near the surface requires knowledge of the local dynamics of the planetary boundary layer. In this paper, we study the relationship between the boundary layer dynamics and the CO2 budget in convective conditions through a newly

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.

Mixed convection boundary-layer flow from a horizontal circular cylinder with a constant surface heat flux

Energy Technology Data Exchange (ETDEWEB)

Nazar, R.; Amin, N. [Department of Mathematics, Universiti Teknologi Malaysia, 81310 Johor Bahru, Johor (Malaysia); Pop, I. [Faculty of Mathematics, University of Cluj, R-3400 Cluj, CP 253 (Romania)

2004-02-01

The laminar mixed convection boundary-layer flow of a viscous and incompressible fluid past a horizontal circular cylinder, which is maintained at a constant heat flux and is placed in a stream flowing vertically upward has been theoretically studied in this paper. The solutions for the flow and heat transfer characteristics are evaluated numerically for different values of the mixed convection parameter {lambda} with the Prandtl number Pr = 1 and 7, respectively. It is found, as for the case of a heated or cooled cylinder, considered by Merkin [5], that assisting flow delays separation of the boundary-layer and can, if the assisting flow is strong enough, suppress it completely. The opposing flow, on the other side, brings the separation point nearer to the lower stagnation point and for sufficiently strong opposing flows there will not be a boundary-layer on the cylinder. (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

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

An analytical solution for the Marangoni mixed convection boundary layer flow

DEFF Research Database (Denmark)

Moghimi, M. A.; Kimiaeifar, Amin; Rahimpour, M.

2010-01-01

In this article, an analytical solution for a Marangoni mixed convection boundary layer flow is presented. A similarity transform reduces the Navier-Stokes equations to a set of nonlinear ordinary differential equations, which are solved analytically by means of the homotopy analysis method (HAM...... the convergence of the solution. The numerical solution of the similarity equations is developed and the results are in good agreement with the analytical results based on the HAM....

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.

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.

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

Energy Technology Data Exchange (ETDEWEB)

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

2013-11-15

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

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

International Nuclear Information System (INIS)

Kaikina, Elena I.

2013-01-01

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

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

OpenAIRE

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

2017-01-01

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

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

Czech Academy of Sciences Publication Activity Database

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

2011-01-01

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

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

A theoretical study of mixing downstream of transverse injection into a supersonic boundary layer

Science.gov (United States)

Baker, A. J.; Zelazny, S. W.

1972-01-01

A theoretical and analytical study was made of mixing downstream of transverse hydrogen injection, from single and multiple orifices, into a Mach 4 air boundary layer over a flat plate. Numerical solutions to the governing three-dimensional, elliptic boundary layer equations were obtained using a general purpose computer program. Founded upon a finite element solution algorithm. A prototype three-dimensional turbulent transport model was developed using mixing length theory in the wall region and the mass defect concept in the outer region. Excellent agreement between the computed flow field and experimental data for a jet/freestream dynamic pressure ratio of unity was obtained in the centerplane region of the single-jet configuration. Poorer agreement off centerplane suggests an inadequacy of the extrapolated two-dimensional turbulence model. Considerable improvement in off-centerplane computational agreement occured for a multi-jet configuration, using the same turbulent transport model.

Modified Differential Transform Method for Two Singular Boundary Values Problems

Directory of Open Access Journals (Sweden)

Yinwei Lin

2014-01-01

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

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

Directory of Open Access Journals (Sweden)

José Játem

2015-12-01

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

Bibliography on moving boundary problems with key word index

International Nuclear Information System (INIS)

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

1979-10-01

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

Bibliography on moving boundary problems with key word index

Energy Technology Data Exchange (ETDEWEB)

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

1979-10-01

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

A new model for the redundancy allocation problem with component mixing and mixed redundancy strategy

International Nuclear Information System (INIS)

Gholinezhad, Hadi; Zeinal Hamadani, Ali

2017-01-01

This paper develops a new model for redundancy allocation problem. In this paper, like many recent papers, the choice of the redundancy strategy is considered as a decision variable. But, in our model each subsystem can exploit both active and cold-standby strategies simultaneously. Moreover, the model allows for component mixing such that components of different types may be used in each subsystem. The problem, therefore, boils down to determining the types of components, redundancy levels, and number of active and cold-standby units of each type for each subsystem to maximize system reliability by considering such constraints as available budget, weight, and space. Since RAP belongs to the NP-hard class of optimization problems, a genetic algorithm (GA) is developed for solving the problem. Finally, the performance of the proposed algorithm is evaluated by applying it to a well-known test problem from the literature with relatively satisfactory results. - Highlights: • A new model for the redundancy allocation problem in series–parallel systems is proposed. • The redundancy strategy of each subsystem is considered as a decision variable and can be active, cold-standby or mixed. • Component mixing is allowed, in other words components of any subsystem can be non-identical. • A genetic algorithm is developed for solving the problem. • Computational experiments demonstrate that the new model leads to interesting results.

Asymptotic boundary value problems for evolution inclusions

Directory of Open Access Journals (Sweden)

Fürst Tomáš

2006-01-01

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

Asymptotic boundary value problems for evolution inclusions

Directory of Open Access Journals (Sweden)

Tomáš Fürst

2006-02-01

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

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

International Nuclear Information System (INIS)

Ignatyev, M. Yu.

2013-01-01

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

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

Energy Technology Data Exchange (ETDEWEB)

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

2010-05-07

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

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

International Nuclear Information System (INIS)

Glowinski, R.; Le Tallec, P.

1984-01-01

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

Local Similarity in the Stable Boundary Layer and Mixing-Length Approaches : Consistency of Concepts

NARCIS (Netherlands)

Van de Wiel, B.J.H.; Moene, A.F.; De Ronde, W.H.; Jonker, H.J.J.

2008-01-01

In stably stratified flows vertical movement of eddies is limited by the fact that kinetic energy is converted into potential energy, leading to a buoyancy displacement scale z B . Our new mixing-length concept for turbulent transport in the stable boundary layer follows a rigid-wall analogy, in the

Local similarity in the stable boundary layer and mixing-length approaches: consistency of concepts

NARCIS (Netherlands)

Wiel, van de B.J.H.; Moene, A.F.; Ronde, W.H.; Jonker, H.J.J.

2008-01-01

In stably stratified flows vertical movement of eddies is limited by the fact that kinetic energy is converted into potential energy, leading to a buoyancy displacement scale z B . Our new mixing-length concept for turbulent transport in the stable boundary layer follows a rigid-wall analogy, in the

Local similarity in the stable boundary layer and mixing-length approaches : consistency of concepts

NARCIS (Netherlands)

Wiel, van de B.J.H.; Moene, A.F.; Ronde, de W.H.; Jonker, H.J.J.

2008-01-01

In stably stratified flows vertical movement of eddies is limited by the fact that kinetic energy is converted into potential energy, leading to a buoyancy displacement scale zB. Our new mixing-length concept for turbulent transport in the stable boundary layer follows a rigid-wall analogy, in the

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

International Nuclear Information System (INIS)

Nazarov, S A

1999-01-01

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

Boundary-value problems in ODE

Science.gov (United States)

Tanriverdi, Tanfer

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

The unbounded solution of a periodic mixed Sturm–Liouville problem in an infinite strip for the Laplacian

Directory of Open Access Journals (Sweden)

M.G. Elsheikh

2013-10-01

Full Text Available The unbounded solution, at the points where the boundary conditions change, for a mixed Sturm–Liouville problem of the Dirichlet–Neumann type can be obtained using the method of the integral equation formulation. Since this formulation is usually reduced to an infinite algebraic system in which the unknowns are the Fourier coefficients of the unknown unbounded entity, a study of ℓp-solutions imposes itself concerning the influence of the truncation on such systems. This study is achieved and the well-known theorem on the ℓ2-solutions of the infinite algebraic systems is generalized.

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

Science.gov (United States)

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

2018-03-01

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

Functional geometric method for solving free boundary problems for harmonic functions

Energy Technology Data Exchange (ETDEWEB)

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

2010-01-01

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

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

Directory of Open Access Journals (Sweden)

Abdeslem Lyaghfouri

2015-07-01

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

Variational methods for boundary value problems for systems of elliptic equations

CERN Document Server

Lavrent'ev, M A

2012-01-01

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

Parallel algorithms for boundary value problems

Science.gov (United States)

Lin, Avi

1991-01-01

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

Student Solutions Manual to Boundary Value Problems and Partial Differential Equations

CERN Document Server

Powers, David L

2005-01-01

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

Boundary value problem for Caputo-Hadamard fractional differential equations

Directory of Open Access Journals (Sweden)

Yacine Arioua

2017-09-01

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

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

Directory of Open Access Journals (Sweden)

Andrzej Nowakowski

2002-05-01

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

Numerical Solution of Mixed Problems of the Theory of Elasticity with One-Sided Constraints

Directory of Open Access Journals (Sweden)

I. V. Stankevich

2017-01-01

enable us to ensure continuous approximation of not only displacements, but also stresses and strains. Mixed schemes to solve the boundary value problems lead to the saddle-point problems. Their solutions use various iterative techniques. One of the most effective techniques is a modified SSOR (MSSOR method, based on the SOR (Successive Over Relaxation one.The paper considers one of the options of the finite element method in the framework of mixed scheme that uses a Reissner functional. The procedures of the algorithm proposed in the paper are used to solve the problem of contact interaction when an elastic body of the finite dimensions, being under a load of the external forces, relies on the absolutely rigid half-space. The contact occurs with the distinguished contact surface, which in the general case can change its size during thermo-mechanical loading. The algorithm is implemented as an application software complex. The numerical study of the one-sided contact interaction between the elastic plate and the perfectly rigid half-space has shown a fairly high efficiency of the developed algorithm and the code that implements it.

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

Directory of Open Access Journals (Sweden)

FAOUZI HADDOUCHI

2015-11-01

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

Solving inverse two-point boundary value problems using collage coding

Science.gov (United States)

Kunze, H.; Murdock, S.

2006-08-01

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

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

CERN Document Server

Chabrowski, Jan

1991-01-01

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

Mixed isogeometric finite cell methods for the stokes problem

NARCIS (Netherlands)

Hoang, T.; Verhoosel, C.V.; Auricchio, F.; van Brummelen, E.H.; Reali, A.

2017-01-01

We study the application of the Isogeometric Finite Cell Method (IGA-FCM) to mixed formulations in the context of the Stokes problem. We investigate the performance of the IGA-FCM when utilizing some isogeometric mixed finite elements, namely: Taylor-Hood, Sub-grid, Raviart-Thomas, and Nédélec

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)

Stability of mixing layers

Science.gov (United States)

Tam, Christopher; Krothapalli, A

1993-01-01

The research program for the first year of this project (see the original research proposal) consists of developing an explicit marching scheme for solving the parabolized stability equations (PSE). Performing mathematical analysis of the computational algorithm including numerical stability analysis and the determination of the proper boundary conditions needed at the boundary of the computation domain are implicit in the task. Before one can solve the parabolized stability equations for high-speed mixing layers, the mean flow must first be found. In the past, instability analysis of high-speed mixing layer has mostly been performed on mean flow profiles calculated by the boundary layer equations. In carrying out this project, it is believed that the boundary layer equations might not give an accurate enough nonparallel, nonlinear mean flow needed for parabolized stability analysis. A more accurate mean flow can, however, be found by solving the parabolized Navier-Stokes equations. The advantage of the parabolized Navier-Stokes equations is that its accuracy is consistent with the PSE method. Furthermore, the method of solution is similar. Hence, the major part of the effort of the work of this year has been devoted to the development of an explicit numerical marching scheme for the solution of the Parabolized Navier-Stokes equation as applied to the high-seed mixing layer problem.

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

Directory of Open Access Journals (Sweden)

Yanmei Sun

2012-01-01

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

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

Directory of Open Access Journals (Sweden)

Grzymkowski R.

2013-03-01

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

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

Directory of Open Access Journals (Sweden)

R. Grzymkowski

2013-01-01

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

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

International Nuclear Information System (INIS)

Barucq, Helene; Bekkey, Chokri; Djellouli, Rabia

2004-01-01

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

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

International Nuclear Information System (INIS)

Mohammed, M.H.H.

2012-01-01

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

The numerical solution of boundary value problems over an infinite domain

International Nuclear Information System (INIS)

Shepherd, M.; Skinner, R.

1976-01-01

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

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.

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

Directory of Open Access Journals (Sweden)

Zhi-Wei Lv

2010-06-01

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

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

International Nuclear Information System (INIS)

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

1998-01-01

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

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

African Journals Online (AJOL)

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

Moving-boundary problems for the time-fractional diffusion equation

Directory of Open Access Journals (Sweden)

Sabrina D. Roscani

2017-02-01

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

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

Science.gov (United States)

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

1974-01-01

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

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

OpenAIRE

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

2006-01-01

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

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

Directory of Open Access Journals (Sweden)

Rakhim Aitbayev

2013-10-01

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

Incremental localized boundary-domain integro-differential equations of elastic damage mechanics for inhomogeneous body

OpenAIRE

Mikhailov, SE

2006-01-01

Copyright @ 2006 Tech Science Press A quasi-static mixed boundary value problem of elastic damage mechanics for a continuously inhomogeneous body is considered. Using the two-operator Green-Betti formula and the fundamental solution of an auxiliary homogeneous linear elasticity with frozen initial, secant or tangent elastic coe±cients, a boundary-domain integro-differential formulation of the elasto-plastic problem with respect to the displacement rates and their gradients is derived. Usin...

Mixed convection boundary layer flow over a vertical surface embedded in a thermally stratified porous medium

International Nuclear Information System (INIS)

Ishak, Anuar; Nazar, Roslinda; Pop, Ioan

2008-01-01

The mixed convection boundary layer flow through a stable stratified porous medium bounded by a vertical surface is investigated. The external velocity and the surface temperature are assumed to vary as x m , where x is measured from the leading edge of the vertical surface and m is a constant. Numerical solutions for the governing Darcy and energy equations are obtained. The results indicate that the thermal stratification significantly affects the surface shear stress as well as the surface heat transfer, besides delays the boundary layer separation

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

Discrete maximum principle for Poisson equation with mixed boundary conditions solved by hp-FEM

Czech Academy of Sciences Publication Activity Database

Vejchodský, Tomáš; Šolín, P.

2009-01-01

Roč. 1, č. 2 (2009), s. 201-214 ISSN 2070-0733 R&D Projects: GA AV ČR IAA100760702; GA ČR(CZ) GA102/07/0496; GA ČR GA102/05/0629 Institutional research plan: CEZ:AV0Z10190503 Keywords : discrete maximum principle * hp-FEM * Poisson equation * mixed boundary conditions Subject RIV: BA - General Mathematics

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

NARCIS (Netherlands)

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

2007-01-01

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

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

Directory of Open Access Journals (Sweden)

Shihuang Hong

2003-03-01

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

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

International Nuclear Information System (INIS)

Alexeyeva, L.A.

2001-01-01

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

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

Directory of Open Access Journals (Sweden)

Lukáš Ladislav

2017-01-01

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

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

Energy Technology Data Exchange (ETDEWEB)

Bjoentegaard, Tormod

2008-04-15

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

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.

Topological Design for Acoustic-Structure Interaction Problems with a Mixed Finite Element Method

DEFF Research Database (Denmark)

Yoon, Gil Ho; Jensen, Jakob Søndergaard; Sigmund, Ole

2006-01-01

to subdomain interfaces evolving during the optimization process. In this paper, we propose to use a mixed finite element formulation with displacements and pressure as primary variables (u/p formulation) which eliminates the need for explicit boundary representation. In order to describe the Helmholtz...... equation and the linear elasticity equation, the mass density as well as the shear and bulk moduli are interpolated with the design variables. In this formulation, the coupled interface boundary conditions are automatically satisfied without having to compute surface coupling integrals. Two dimensional...

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

Topology optimization for acoustic-structure interaction problems

DEFF Research Database (Denmark)

Yoon, Gil Ho; Jensen, Jakob Søndergaard; Sigmund, Ole

2006-01-01

We propose a gradient based topology optimization algorithm for acoustic-structure (vibro-acoustic) interaction problems without an explicit interfacing boundary representation. In acoustic-structure interaction problems, the pressure field and the displacement field are governed by the Helmholtz...... to subdomain interfaces evolving during the optimization process. In this paper, we propose to use a mixed finite element formulation with displacements and pressure as primary variables (u/p formulation) which eliminates the need for explicit boundary representation. In order to describe the Helmholtz......-dimensional acoustic-structure interaction problems are optimized to show the validity of the proposed method....

Numerical solution of the problems for plates on partial internal supports of complicated configurations

International Nuclear Information System (INIS)

Quang A, Dang; Hai, Truong Ha

2014-01-01

Very recently in the work S imple Iterative Method for Solving Problems for Plates with Partial Internal Supports, Journal of Engineering Mathematics, DOI: 10.1007/s10665-013-9652-7 (in press) , we proposed a numerical method for solving some problems of plates on one and two line partial internal supports (LPIS). In the essence they are problems with strongly mixed boundary conditions for biharmonic equation. Using this method we reduced the problems to a sequence of boundary value problems for the Poisson equation with weakly mixed boundary conditions, which are easily solved numerically. The advantages of the method over other ones were shown. In this paper we apply the method to plates on internal supports of more complicated configurations. Namely, we consider the case of three LPIS and the case of the cross support. The convergence of the method is established theoretically and its efficiency is confirmed on numerical experiments

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

International Nuclear Information System (INIS)

Sermange, M.

1982-06-01

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

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

Science.gov (United States)

Nye, Mary-Jo

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

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

Directory of Open Access Journals (Sweden)

Moulay Rchid Sidi Ammi

2013-01-01

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

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

Directory of Open Access Journals (Sweden)

Alsaidi M. Altaher

2012-01-01

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

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

Science.gov (United States)

Kot, V. A.

2017-11-01

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

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

Directory of Open Access Journals (Sweden)

Svatoslav Stanêk

2008-03-01

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

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

CERN Document Server

Zhu, C

2003-01-01

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

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

International Nuclear Information System (INIS)

Zhu, Changjiang; Duan, Renjun

2003-01-01

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

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

Czech Academy of Sciences Publication Activity Database

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

2017-01-01

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

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

Czech Academy of Sciences Publication Activity Database

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

2017-01-01

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

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

Directory of Open Access Journals (Sweden)

Seyyed Mohsen Khalkhali

2012-09-01

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

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

International Nuclear Information System (INIS)

Chen Yunmei.

1987-06-01

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

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

KAUST Repository

Dujardin, G. M.

2009-08-12

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

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

International Nuclear Information System (INIS)

Boisseau, Bruno; Forgacs, Peter; Giacomini, Hector

2007-01-01

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

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

International Nuclear Information System (INIS)

Mugge, J.W.

1979-10-01

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

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

Energy Technology Data Exchange (ETDEWEB)

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

2002-07-01

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

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

Directory of Open Access Journals (Sweden)

Martin Eggensperger

2004-07-01

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

A Viscous Fluid Flow through a Thin Channel with Mixed Rigid-Elastic Boundary: Variational and Asymptotic Analysis

Directory of Open Access Journals (Sweden)

R. Fares

2012-01-01

Full Text Available We study the nonsteady Stokes flow in a thin tube structure composed by two thin rectangles with lateral elastic boundaries which are connected by a domain with rigid boundaries. After a variational approach of the problem which gives us existence, uniqueness, regularity results, and some a priori estimates, we construct an asymptotic solution. The existence of a junction region between the two rectangles imposes to consider, as part of the asymptotic solution, some boundary layer correctors that correspond to this region. We present and solve the problems for all the terms of the asymptotic expansion. For two different cases, we describe the order of steps of the algorithm of solving the problem and we construct the main term of the asymptotic expansion. By means of the a priori estimates, we justify our asymptotic construction, by obtaining a small error between the exact and the asymptotic solutions.

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

KAUST Repository

Barenblatt, G. I.

2013-01-16

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

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

Science.gov (United States)

Yousef, Hamood Mohammed; Ismail, Ahmad Izani

2017-11-01

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

Chebyshev Finite Difference Method for Fractional Boundary Value Problems

Directory of Open Access Journals (Sweden)

Boundary

2015-09-01

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

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

Science.gov (United States)

Hromadka, T.V.

1984-01-01

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

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

International Nuclear Information System (INIS)

Lin Bin; Li Kaitai; Cheng Zhengxing

2009-01-01

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

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

International Nuclear Information System (INIS)

Nkemzi, B.

2005-10-01

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

Multiparametric programming based algorithms for pure integer and mixed-integer bilevel programming problems

KAUST Repository

Domínguez, Luis F.

2010-12-01

This work introduces two algorithms for the solution of pure integer and mixed-integer bilevel programming problems by multiparametric programming techniques. The first algorithm addresses the integer case of the bilevel programming problem where integer variables of the outer optimization problem appear in linear or polynomial form in the inner problem. The algorithm employs global optimization techniques to convexify nonlinear terms generated by a reformulation linearization technique (RLT). A continuous multiparametric programming algorithm is then used to solve the reformulated convex inner problem. The second algorithm addresses the mixed-integer case of the bilevel programming problem where integer and continuous variables of the outer problem appear in linear or polynomial forms in the inner problem. The algorithm relies on the use of global multiparametric mixed-integer programming techniques at the inner optimization level. In both algorithms, the multiparametric solutions obtained are embedded in the outer problem to form a set of single-level (M)(I)(N)LP problems - which are then solved to global optimality using standard fixed-point (global) optimization methods. Numerical examples drawn from the open literature are presented to illustrate the proposed algorithms. © 2010 Elsevier Ltd.

Transient reaction of an elastic half-plane on a source of a concentrated boundary disturbance

Science.gov (United States)

Okonechnikov, A. S.; Tarlakovski, D. V.; Ul'yashina, A. N.; Fedotenkov, G. V.

2016-11-01

One of the key problems in studying the non-stationary processes of solid mechanics is obtaining of influence functions. These functions serve as solutions for the problems of effect of sudden concentrated loads on a body with linear elastic properties. Knowledge of the influence functions allows us to obtain the solutions for the problems with non-mixed boundary and initial conditions in the form of quadrature formulae with the help of superposition principle, as well as get the integral governing equations for the problems with mixed boundary and initial conditions. This paper offers explicit derivations for all nonstationary surface influence functions of an elastic half-plane in a plane strain condition. It is achieved with the help of combined inverse transform of a Fourier-Laplace integral transformation. The external disturbance is both dynamic and kinematic. The derived functions in xτ-domain are studied to find and describe singularities and are supplemented with graphs.

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

KAUST Repository

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

2016-01-01

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

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

KAUST Repository

Jin, Bangti

2016-02-01

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

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

Directory of Open Access Journals (Sweden)

Yansheng Liu

2009-11-01

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

The boundary value problem for discrete analytic functions

KAUST Repository

Skopenkov, Mikhail

2013-06-01

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

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

Directory of Open Access Journals (Sweden)

Kyncl Martin

2017-01-01

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

Partial differential equations and boundary-value problems with applications

CERN Document Server

Pinsky, Mark A

2011-01-01

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

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

Directory of Open Access Journals (Sweden)

D. A. Eliseev

2015-01-01

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

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

NARCIS (Netherlands)

Kakuba, G.

2011-01-01

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

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

Directory of Open Access Journals (Sweden)

J. Wu

2012-01-01

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

Effect of Variable Viscosity on Vortex Instability of Non-Darcy Mixed Convection Boundary Layer Flow Adjacent to a Nonisothermal Horizontal Surface in a Porous Medium

Directory of Open Access Journals (Sweden)

A. M. Elaiw

2012-01-01

Full Text Available We study the effect of variable viscosity on the flow and vortex instability for non-Darcy mixed convection boundary layer flow on a nonisothermal horizontal plat surface in a saturated porous medium. The variation of viscosity is expressed as an exponential function of temperature. The analysis of the disturbance flow is based on linear stability theory. The base flow equations and the resulting eigenvalue problem are solved using finite difference schemes. It is found that the variable viscosity effect enhances the heat transfer rate and destabilizes the flow for liquid heating, while the opposite trend is true for gas heating.

Application of the Least Squares Method in Axisymmetric Biharmonic Problems

Directory of Open Access Journals (Sweden)

Vasyl Chekurin

2016-01-01

Full Text Available An approach for solving of the axisymmetric biharmonic boundary value problems for semi-infinite cylindrical domain was developed in the paper. On the lateral surface of the domain homogeneous Neumann boundary conditions are prescribed. On the remaining part of the domain’s boundary four different biharmonic boundary pieces of data are considered. To solve the formulated biharmonic problems the method of least squares on the boundary combined with the method of homogeneous solutions was used. That enabled reducing the problems to infinite systems of linear algebraic equations which can be solved with the use of reduction method. Convergence of the solution obtained with developed approach was studied numerically on some characteristic examples. The developed approach can be used particularly to solve axisymmetric elasticity problems for cylindrical bodies, the heights of which are equal to or exceed their diameters, when on their lateral surface normal and tangential tractions are prescribed and on the cylinder’s end faces various types of boundary conditions in stresses in displacements or mixed ones are given.

hp Spectral element methods for three dimensional elliptic problems

Indian Academy of Sciences (India)

elliptic boundary value problems on non-smooth domains in R3. For Dirichlet problems, ... of variable degree bounded by W. Let N denote the number of layers in the geomet- ric mesh ... We prove a stability theorem for mixed problems when the spectral element functions vanish ..... Applying Theorem 3.1,. ∫ r l. |Mu|2dx −.

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

Modelling the development of mixing height in near equatorial region

Energy Technology Data Exchange (ETDEWEB)

Samah, A.A. [Univ. of Malaya, Air Pollution Research Unit, Kuala Lumpur (Malaysia)

1997-10-01

Most current air pollution models were developed for mid-latitude conditions and as such many of the empirical parameters used were based on observations taken in the mid-latitude boundary layer which is physically different from that of the equatorial boundary layer. In the equatorial boundary layer the Coriolis parameter f is small or zero and moisture plays a more important role in the control of stability and the surface energy balance. Therefore air pollution models such as the OMLMULTI or the ADMS which were basically developed for mid-latitude conditions must be applied with some caution and would need some adaptation to properly simulate the properties of equatorial boundary layer. This work elucidates some of the problems of modelling the evolution of mixing height in the equatorial region. The mixing height estimates were compared with routine observations taken during a severe air pollution episodes in Malaysia. (au)

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.

The Jump Problem for Mixed-Type Equations with Defects on the Type Change Line

Directory of Open Access Journals (Sweden)

Ahmed Maher

2010-01-01

Full Text Available The jump problem and problems with defects on the type change line for model mixed-type equations in the mixed domains are investigated. The explicit solutions of the jump problem are obtained by the method of integral equations and by the Fourier transformation method. The problems with defects are reduced to singular integral equations. Some results for the solution of the equation under consideration are discussed concerning the existence and uniqueness for the solution of the suggested problem.

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

Directory of Open Access Journals (Sweden)

Chein-Shan Liu

2013-01-01

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

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

OpenAIRE

Belmiloudi, A.; Mahé, F.

2014-01-01

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

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

Science.gov (United States)

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

2011-04-01

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

Experiences of ZAMG on mixing height determination

Energy Technology Data Exchange (ETDEWEB)

Piringer, M. [Zentralanstalt fuer Meteorologie und Geodynamik, ZAMG, Vienna (Austria)

1997-10-01

Temperature inversions in the boundary layer occur quite often, esp. in mountainous terrain by which Austria is covered to a large extent, and can lead to enhanced pollution at the surface because the air volume available for dilution is then vertically limited. The Department of Environmental Meteorology of ZAMG therefore set up several field programs in the past to study such conditions at a variety of sites in Austria, using tethersondes and Sodars. Early investigations aimed at comparing Sodar echo profiles to the tethersonde temperature profiles to derive mixing heights from the Sodar echo structure. More recently, evolving from KONGEX, the `convective boundary layer experiment`, mixing heights calculated for Vienna by the OML model were compared to those derived from radiosonde and tethersonde potential temperature profiles. Results of these investigations will be presented, focussing on the problems when using the different methods. New efforts to derive mixing heights from data were also undertaken and are discussed separately. (au)

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

International Nuclear Information System (INIS)

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

1993-01-01

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

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.

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

International Nuclear Information System (INIS)

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

2011-01-01

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

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

Directory of Open Access Journals (Sweden)

Fuyi Xu

2011-12-01

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

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

CERN Document Server

Amenta, Alex

2018-01-01

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

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

OpenAIRE

Kovalenko, S. S.

2014-01-01

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

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

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.

Engineering Mixed Ionic Electronic Conduction in La 0.8 Sr 0.2 MnO 3+ δ Nanostructures through Fast Grain Boundary Oxygen Diffusivity

KAUST Repository

Saranya, Aruppukottai M.; Pla, Dolors; Morata, Alex; Cavallaro, Andrea; Canales-Vá zquez, Jesú s; Kilner, John A.; Burriel, Mó nica; Tarancó n, Albert

2015-01-01

to implement in nanostructures. Here, an artificial mixed ionic electronic conducting oxide is fabricated by grain boundary (GB) engineering thin films of La0.8Sr0.2MnO3+δ. This electronic conductor is converted into a good mixed ionic electronic conductor

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

Representing Sheared Convective Boundary Layer by Zeroth- and First-Order-Jump Mixed-Layer Models: Large-Eddy Simulation Verification

NARCIS (Netherlands)

Pino, D.; Vilà-Guerau de Arellano, J.; Kim, S.W.

2006-01-01

Dry convective boundary layers characterized by a significant wind shear on the surface and at the inversion are studied by means of the mixed-layer theory. Two different representations of the entrainment zone, each of which has a different closure of the entrainment heat flux, are considered. The

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.

Mixing height measurements from UHF wind profiling radar

Energy Technology Data Exchange (ETDEWEB)

Angevine, W.M.; Grimsdell, A.W. [CIRES, Univ. of Colorado, and NOAA Aeronomy Lab., Boulder, Colorado (United States)

1997-10-01

Mixing height in convective boundary layers can be detected by wind profiling radars (profilers) operating at or near 915 MHZ. We have made such measurements in a variety of settings including Alabama in 1992; Nova Scotia, Canada, during the North Atlantic Regional Experiment (NARE) 1993; Tennessee during the Southern Oxidant Study (SOS) 1994; near a 450 m tower in Wisconsin in 1995; and extensively in Illinois during the Flatland95, `96, and `97 experiments, as well as continuous operations at the Flatland Atmospheric Observatory. Profiler mixing height measurements, like all measurements, are subject to some limitations. The most important of these are due to rainfall, minimum height, and height resolution. Profilers are very sensitive to rain, which dominates the reflectivity and prevents the mixing height from being detected. Because the best height resolution is currently 60 m and the minimum height is 120-150 m AGL, the profiler is not suited for detecting mixing height in stable or nocturnal boundary layers. Problems may also arise in very dry or cold environments. (au) 12 refs.

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.

The effect of adiabatic and conducting wall boundary conditions on LES of a thermal mixing tee

International Nuclear Information System (INIS)

Howard, Richard J.A.; Pasutto, Thomas

2009-01-01

In this paper preliminary LES simulations are carried out of the FATHERINO mixing T junction experiment. In this experiment 80degC hot water enters a lateral steel pipe which has a diameter of D=0.054m, at a speed of 1.04m/s and meets 5degC cold water which enters a perpendicular steel pipe branch that also has a diameter D=0.054m but this time at a lower speed of 0.26m/s. The modelling of the steel pipe walls is tested by comparing adiabatic and 1D conducting wall boundary conditions. The numerical grid used contains approximately 440,000 hexahedral elements. The near wall refinement is not sufficient to resolve the near wall boundary layer (y + approx. = 32) and a standard logarithmic boundary condition is used. A method known as the synthetic eddy method is used to generate the turbulent flow at the pipe inlets. Three different LES models are used (Smagorinsky, dynamic Smagorinsky and wale) to resolve the subgrid turbulent motion beyond the wall grid. An additional test is carried out where no subgrid model is used with only the wall modelling being applied. The results show that the wale model generates much less resolved turbulence than the other cases and this model shows virtually no difference between the two methods of wall thermal modelling. The dynamic Smagorinsky model shows that, downstream of the mixing T, the lower wall remains at a lower temperature for longer when the adiabatic boundary condition is applied. The Smagorinsky model is found to produce the highest level of resolved temperature fluctuation. For this model the 1D thermal modelling approach increases the unsteadiness of both the velocity and temperature fields at the onset of the mixing and in the middle of the pipe downstream of the T junction. However near the lower wall the 1D thermal modelling approach tends to reduce the unsteadiness. The case with no subgrid modelling shows higher levels of turbulence kinetic energy but lower levels of temperature fluctuation than the cases with

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)

Variational P1 approximations of general-geometry multigroup transport problems

International Nuclear Information System (INIS)

Rulko, R.P.; Tomasevic, D.; Larsen, E.W.

1995-01-01

A variational approximation is developed for general-geometry multigroup transport problems with arbitrary anisotropic scattering. The variational principle is based on a functional that approximates a reaction rate in a subdomain of the system. In principle, approximations that result from this functional ''optimally'' determine such reaction rates. The functional contains an arbitrary parameter α and requires the approximate solutions of a forward and an adjoint transport problem. If the basis functions for the forward and adjoint solutions are chosen to be linear functions of the angular variable Ω, the functional yields the familiar multigroup P 1 equations for all values of α. However, the boundary conditions that result from the functional depend on α. In particular, for problems with vacuum boundaries, one obtains the conventional mixed boundary condition, but with an extrapolation distance that depends continuously on α. The choice α = 0 yields a generalization of boundary conditions derived earlier by Federighi and Pomraning for a more limited class of problems. The choice α = 1 yields a generalization of boundary conditions derived previously by Davis for monoenergetic problems. Other boundary conditions are obtained by choosing different values of α. The authors discuss this indeterminancy of α in conjunction with numerical experiments

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

Multiparametric programming based algorithms for pure integer and mixed-integer bilevel programming problems

KAUST Repository

Domí nguez, Luis F.; Pistikopoulos, Efstratios N.

2010-01-01

continuous multiparametric programming algorithm is then used to solve the reformulated convex inner problem. The second algorithm addresses the mixed-integer case of the bilevel programming problem where integer and continuous variables of the outer problem

Discrete quintic spline for boundary value problem in plate deflation theory

Science.gov (United States)

Wong, Patricia J. Y.

2017-07-01

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

The Boundary Function Method. Fundamentals

Science.gov (United States)

Kot, V. A.

2017-03-01

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

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

Directory of Open Access Journals (Sweden)

V. Rukavishnikov

2014-01-01

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

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

Directory of Open Access Journals (Sweden)

Muhammad Aslam Noor

2008-01-01

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

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

International Nuclear Information System (INIS)

Ji, X.; Chen, Y.M.

1989-01-01

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

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

International Nuclear Information System (INIS)

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

1995-01-01

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

The nonlocal problem for a hyperbolic equation with Bessel operator in a rectangular domain

Directory of Open Access Journals (Sweden)

Natalya V. Zaitseva

2016-12-01

Full Text Available We consider a boundary value problem for a hyperbolic equation with Bessel differential operator in a rectangular domain with integral nonlocal boundary value condition of the first kind. The equivalence between boundary value problem with integral nonlocal condition of the first kind and a local boundary value problem with mixed boundary conditions of the first and third kinds is proved. The existence and uniqueness of solution of the equivalent problem are established by means of the spectral method. At the uniqueness proof the completeness of the eigenfunction system of the spectral problem is used . At the existence proof the assessment of coefficients of series, the asymptotic formula for Bessel function of the first kind and asymptotic formula for eigenvalues are used. Sufficient conditions on the functions defining initial data of the problem are received. The solution of the problem is obtained in explicit form. The solution is obtained in the form of the Fourier–Bessel series. Its convergence is proved in the class of regular solutions.

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.

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

Boundary control of fluid flow through porous media

DEFF Research Database (Denmark)

Hasan, Agus; Foss, Bjarne; Sagatun, Svein Ivar

2010-01-01

The flow of fluids through porous media can be described by the Boussinesq’s equation with mixed boundary conditions; a Neumann’s boundary condition and a nonlinear boundary condition. The nonlinear boundary condition provides a means to control the fluid flow through porous media. In this paper,......, some stabilizing controllers are constructed for various cases using Lyapunov design.......The flow of fluids through porous media can be described by the Boussinesq’s equation with mixed boundary conditions; a Neumann’s boundary condition and a nonlinear boundary condition. The nonlinear boundary condition provides a means to control the fluid flow through porous media. In this paper...

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

Czech Academy of Sciences Publication Activity Database

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

2017-01-01

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

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

Czech Academy of Sciences Publication Activity Database

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

2017-01-01

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

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

Directory of Open Access Journals (Sweden)

Archana Chauhan

2012-12-01

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

Fluid Structure Interaction Techniques For Extrusion And Mixing Processes

Science.gov (United States)

Valette, Rudy; Vergnes, Bruno; Coupez, Thierry

2007-05-01

This work focuses on the development of numerical techniques devoted to the simulation of mixing processes of complex fluids such as twin-screw extrusion or batch mixing. In mixing process simulation, the absence of symmetry of the moving boundaries (the screws or the rotors) implies that their rigid body motion has to be taken into account by using a special treatment We therefore use a mesh immersion technique (MIT), which consists in using a P1+/P1-based (MINI-element) mixed finite element method for solving the velocity-pressure problem and then solving the problem in the whole barrel cavity by imposing a rigid motion (rotation) to nodes found located inside the so called immersed domain, each sub-domain (screw, rotor) being represented by a surface CAD mesh (or its mathematical equation in simple cases). The independent meshes are immersed into a unique background computational mesh by computing the distance function to their boundaries. Intersections of meshes are accounted for, allowing to compute a fill factor usable as for the VOF methodology. This technique, combined with the use of parallel computing, allows to compute the time-dependent flow of generalized Newtonian fluids including yield stress fluids in a complex system such as a twin screw extruder, including moving free surfaces, which are treated by a "level set" and Hamilton-Jacobi method.

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.

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

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

Mixed-hybrid finite element method for the transport equation and diffusion approximation of transport problems; Resolution de l'equation du transport par une methode d'elements finis mixtes-hybrides et approximation par la diffusion de problemes de transport

Energy Technology Data Exchange (ETDEWEB)

Cartier, J

2006-04-15

This thesis focuses on mathematical analysis, numerical resolution and modelling of the transport equations. First of all, we deal with numerical approximation of the solution of the transport equations by using a mixed-hybrid scheme. We derive and study a mixed formulation of the transport equation, then we analyse the related variational problem and present the discretization and the main properties of the scheme. We particularly pay attention to the behavior of the scheme and we show its efficiency in the diffusion limit (when the mean free path is small in comparison with the characteristic length of the physical domain). We present academical benchmarks in order to compare our scheme with other methods in many physical configurations and validate our method on analytical test cases. Unstructured and very distorted meshes are used to validate our scheme. The second part of this thesis deals with two transport problems. The first one is devoted to the study of diffusion due to boundary conditions in a transport problem between two plane plates. The second one consists in modelling and simulating radiative transfer phenomenon in case of the industrial context of inertial confinement fusion. (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...

Aviation NOx-induced CH4 effect: Fixed mixing ratio boundary conditions versus flux boundary conditions

Science.gov (United States)

Khodayari, Arezoo; Olsen, Seth C.; Wuebbles, Donald J.; Phoenix, Daniel B.

2015-07-01

Atmospheric chemistry-climate models are often used to calculate the effect of aviation NOx emissions on atmospheric ozone (O3) and methane (CH4). Due to the long (∼10 yr) atmospheric lifetime of methane, model simulations must be run for long time periods, typically for more than 40 simulation years, to reach steady-state if using CH4 emission fluxes. Because of the computational expense of such long runs, studies have traditionally used specified CH4 mixing ratio lower boundary conditions (BCs) and then applied a simple parameterization based on the change in CH4 lifetime between the control and NOx-perturbed simulations to estimate the change in CH4 concentration induced by NOx emissions. In this parameterization a feedback factor (typically a value of 1.4) is used to account for the feedback of CH4 concentrations on its lifetime. Modeling studies comparing simulations using CH4 surface fluxes and fixed mixing ratio BCs are used to examine the validity of this parameterization. The latest version of the Community Earth System Model (CESM), with the CAM5 atmospheric model, was used for this study. Aviation NOx emissions for 2006 were obtained from the AEDT (Aviation Environmental Design Tool) global commercial aircraft emissions. Results show a 31.4 ppb change in CH4 concentration when estimated using the parameterization and a 1.4 feedback factor, and a 28.9 ppb change when the concentration was directly calculated in the CH4 flux simulations. The model calculated value for CH4 feedback on its own lifetime agrees well with the 1.4 feedback factor. Systematic comparisons between the separate runs indicated that the parameterization technique overestimates the CH4 concentration by 8.6%. Therefore, it is concluded that the estimation technique is good to within ∼10% and decreases the computational requirements in our simulations by nearly a factor of 8.

Classical solutions of mixed problems for quasilinear first order PFDEs on a cylindrical domain

Directory of Open Access Journals (Sweden)

Wojciech Czernous

2014-01-01

Full Text Available We abandon the setting of the domain as a Cartesian product of real intervals, customary for first order PFDEs (partial functional differential equations with initial boundary conditions. We give a new set of conditions on the possibly unbounded domain \\(\\Omega\\ with Lipschitz differentiable boundary. Well-posedness is then reliant on a variant of the normal vector condition. There is a neighbourhood of \\(\\partial\\Omega\\ with the property that if a characteristic trajectory has a point therein, then its every earlier point lies there as well. With local assumptions on coefficients and on the free term, we prove existence and Lipschitz dependence on data of classical solutions on \\((0,c\\times\\Omega\\ to the initial boundary value problem, for small \\(c\\. Regularity of solutions matches this domain, and the proof uses the Banach fixed-point theorem. Our general model of functional dependence covers problems with deviating arguments and integro-differential equations.

A Duality Theory for Non-convex Problems in the Calculus of Variations

Science.gov (United States)

Bouchitté, Guy; Fragalà, Ilaria

2018-02-01

We present a new duality theory for non-convex variational problems, under possibly mixed Dirichlet and Neumann boundary conditions. The dual problem reads nicely as a linear programming problem, and our main result states that there is no duality gap. Further, we provide necessary and sufficient optimality conditions, and we show that our duality principle can be reformulated as a min-max result which is quite useful for numerical implementations. As an example, we illustrate the application of our method to a celebrated free boundary problem. The results were announced in Bouchitté and Fragalà (C R Math Acad Sci Paris 353(4):375-379, 2015).

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

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.

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

Science.gov (United States)

Kanoglu, U.; Aydin, B.

2014-12-01

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

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

International Nuclear Information System (INIS)

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

2011-01-01

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

A branch-and-cut-and-price algorithm for the mixed capacitated general routing problem

DEFF Research Database (Denmark)

Bach, Lukas; Wøhlk, Sanne; Lysgaard, Jens

2016-01-01

In this paper, we consider the Mixed Capacitated General Routing Problem which is a combination of the Capacitated Vehicle Routing Problem and the Capacitated Arc Routing Problem. The problem is also known as the Node, Edge, and Arc Routing Problem. We propose a Branch-and-Cut-and-Price algorithm...

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.

Systems of general nonlinear set-valued mixed variational inequalities problems in Hilbert spaces

Directory of Open Access Journals (Sweden)

Cho Yeol

2011-01-01

Full Text Available Abstract In this paper, the existing theorems and methods for finding solutions of systems of general nonlinear set-valued mixed variational inequalities problems in Hilbert spaces are studied. To overcome the difficulties, due to the presence of a proper convex lower semi-continuous function, φ and a mapping g, which appeared in the considered problem, we have used some applications of the resolvent operator technique. We would like to point out that although many authors have proved results for finding solutions of the systems of nonlinear set-valued (mixed variational inequalities problems, it is clear that it cannot be directly applied to the problems that we have considered in this paper because of φ and g. 2000 AMS Subject Classification: 47H05; 47H09; 47J25; 65J15.

Implementing Mixed Method of Peer Teaching and Problem Solving on Undergraduate Students

Directory of Open Access Journals (Sweden)

A. Firli

2017-02-01

Full Text Available This study examined the application of problem solving method combined with student centered learning (peer teaching method as a mixed method to improve student’s passing level of financial management course. The object of this study was the 84 students of financial management course separated within two classes during the odd semester period 2014/2015, July until December 2015 with fourteen meeting courses. Data used to measure the results of the application is mid and final exam scores of both classes. Researcher used observation, interview and documentation as data collect technique also triangulation technique as data validity check. This study used problem solving method combined with student centered learning (peer teaching method as a mixed method which included into the Classroom Action Research. The final results show the increase in class A passing level is 17%. Class B passing level increased 3%. From the research we also know that in practical use of mixed method learning, leader’s quality and conducive learning environment are influencing factors in improving student’s learning performance. While the result confirms that mixed method improving learning performance, this study also founds additional factors that might be considerably affecting the results of learning performance when implementing the mixed method.

Turbulent fluxes in stably stratified boundary layers

International Nuclear Information System (INIS)

L'vov, Victor S; Procaccia, Itamar; Rudenko, Oleksii

2008-01-01

We present here an extended version of an invited talk we gave at the international conference 'Turbulent Mixing and Beyond'. The dynamical and statistical description of stably stratified turbulent boundary layers with the important example of the stable atmospheric boundary layer in mind is addressed. Traditional approaches to this problem, based on the profiles of mean quantities, velocity second-order correlations and dimensional estimates of the turbulent thermal flux, run into a well-known difficulty, predicting the suppression of turbulence at a small critical value of the Richardson number, in contradiction to observations. Phenomenological attempts to overcome this problem suffer from various theoretical inconsistencies. Here, we present an approach taking into full account all the second-order statistics, which allows us to respect the conservation of total mechanical energy. The analysis culminates in an analytic solution of the profiles of all mean quantities and all second-order correlations, removing the unphysical predictions of previous theories. We propose that the approach taken here is sufficient to describe the lower parts of the atmospheric boundary layer, as long as the Richardson number does not exceed an order of unity. For much higher Richardson numbers, the physics may change qualitatively, requiring careful consideration of the potential Kelvin-Helmoholtz waves and their interaction with the vortical turbulence.

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.

Effect of mixing rule boundary conditions on high pressure (liquid + liquid) equilibrium prediction

International Nuclear Information System (INIS)

Hsieh, Min-Kang; Lin, Shiang-Tai

2012-01-01

Highlights: ► Prediction of LLE from the combined use of EOS and liquid model are examined. ► The mixing rule used affects the predicted pressure dependence of LLE. ► MHV1 mixing rule predicts decent LLE at low pressures. ► WS mixing rule predicts more accurate excess volume and LLE at high pressures. ► The hybrid of MHV1 and WS mixing rule gives overall the best predictions. - Abstract: We examine the prediction of high pressure (liquid + liquid) equilibrium (LLE) from the Peng–Robinson equation with three excess Gibbs free energy (G ex )-based mixing rules (MR): the first order modified Huron–Vidal (MHV1), the Wong–Sandler (WS), and a hybrid of these two (referred to as G ex B 2 ). These mixing rules differ by the boundary conditions used for determination of the temperature and composition dependence of parameters a and b in the PR EOS. The condition of matching the excess Gibbs free energy from the EOS at zero pressure to that from the G ex model, used in MHV1 and G ex B 2 MR, leads to a similar miscibility gap from PR EOS and the G ex model used. On the other hand, the condition of matching excess Helmholtz energy from the EOS at infinite pressure to that from the G ex model, used in the WS MR, shows remarkable deviations. The condition of quadratic composition dependence in the second virial coefficient (B 2 ), used in WS and G ex B 2 MR, allows for both positive and negative values in the molar excess volume. Depending on the mixture, either the increase or decrease of the miscibility gap with pressure can be observed when the WS or the G ex B 2 MR is used. The condition of linear combination of molecular sizes of each component used in the MHV1 MR, however, often leads to small, positive molar excess volumes. As a consequence, the predicted LLE from using the MHV1 MR are insensitive to pressure. Therefore, we find that the G ex B 2 mixing rule provides the best predictive power for the LLE over a wide range of temperature and pressure.

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.

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

Is radioactive mixed waste packaging and transportation really a problem

International Nuclear Information System (INIS)

McCall, D.L.; Calihan, T.W. III.

1992-01-01

Recently, there has been significant concern expressed in the nuclear community over the packaging and transportation of radioactive mixed waste under US Department of Transportation regulation. This concern has grown more intense over the last 5 to 10 years. Generators and regulators have realized that much of the waste shipped as ''low-level radioactive waste'' was in fact ''radioactive mixed waste'' and that these wastes pose unique transportation and disposal problems. Radioactive mixed wastes must, therefore, be correctly identified and classed for shipment. If must also be packaged, marked, labeled, and otherwise prepared to ensure safe transportation and meet applicable storage and disposal requirements, when established. This paper discusses regulations applicable to the packaging and transportation of radioactive mixed waste and identifies effective methods that waste shippers can adopt to meet the current transportation requirements. This paper will include a characterization and description of the waste, authorized packaging, and hazard communication requirements during transportation. Case studies will be sued to assist generators in understanding mixed waste shipment requirements and clarify the requirements necessary to establish a waste shipment program. Although management and disposal of radioactive mixed waste is clearly a critical issue, packaging and transportation of these waste materials is well defined in existing US Department of Transportation hazardous material regulations

The problem with time in mixed continuous/discrete time modelling

NARCIS (Netherlands)

Rovers, K.C.; Kuper, Jan; Smit, Gerardus Johannes Maria

The design of cyber-physical systems requires the use of mixed continuous time and discrete time models. Current modelling tools have problems with time transformations (such as a time delay) or multi-rate systems. We will present a novel approach that implements signals as functions of time,

On a mixed problem for a coupled nonlinear system

Directory of Open Access Journals (Sweden)

Marcondes R. Clark

1997-03-01

Full Text Available In this article we prove the existence and uniqueness of solutions to the mixed problem associated with the nonlinear system $$ u_{tt}-M(int_Omega |abla u|^2dxDelta u+|u|^ ho u+heta =f $$ $$ heta _t -Delta heta +u_{t}=g $$ where $M$ is a positive real function, and $f$ and $g$ are known real functions.

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

Directory of Open Access Journals (Sweden)

George N. Galanis

2005-10-01

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

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

International Nuclear Information System (INIS)

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

1991-01-01

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

Observation of interior and boundary-layer mixing processes due to near-inertial waves in a stratified basin without tides

Science.gov (United States)

van der Lee, Eefke; Umlauf, Lars

2010-05-01

Near-inertial waves form an important contribution to oceanic energy and shear spectra, and thus play a major role in mixing the ocean's interior. Here, we compare internal-wave mixing processes in the interior of a stratified basin to those occurring on the sloping boundaries. We use the virtually tideless Baltic Sea as a natural laboratory, allowing us to isolate the effect of near-inertial waves that is otherwise (often) overshadowed by internal tides. The measurements presented here consist of moored ADCPs and CTD loggers in the center of the basin and on the slopes, combined with densely spaced shear-microstructure and ADCP cross-slope transects. During summer stratification, a three-layer density structure, with a thermocline and a deeper halocline, was observed with clear signals of downward near-inertial energy propagation after a short wind event. These motions are interpreted as near-inertial wave modes interacting with the sloping topography. Dissipation rates observed in the center of the basin scale with shear and stratification parameters in the way suggested by MacKinnon and Gregg (2003) for the shelf. On the slopes, microstructure transects reveal a periodic near-bed dissipation rate signal and a growing and decaying bottom boundary layer (BBL) thickness; both signals are triggered by near-bottom currents oscillating with a near-inertial frequency. Near-bottom dissipation rates are greatly enhanced compared to the interior, and, due to the straining of lateral density gradients by the cross-slope velocity, mixing is rather efficient, and contributes significantly to the basin-scale mixing.

Mixed-hybrid finite element method for the transport equation and diffusion approximation of transport problems; Resolution de l'equation du transport par une methode d'elements finis mixtes-hybrides et approximation par la diffusion de problemes de transport

Energy Technology Data Exchange (ETDEWEB)

Cartier, J

2006-04-15

This thesis focuses on mathematical analysis, numerical resolution and modelling of the transport equations. First of all, we deal with numerical approximation of the solution of the transport equations by using a mixed-hybrid scheme. We derive and study a mixed formulation of the transport equation, then we analyse the related variational problem and present the discretization and the main properties of the scheme. We particularly pay attention to the behavior of the scheme and we show its efficiency in the diffusion limit (when the mean free path is small in comparison with the characteristic length of the physical domain). We present academical benchmarks in order to compare our scheme with other methods in many physical configurations and validate our method on analytical test cases. Unstructured and very distorted meshes are used to validate our scheme. The second part of this thesis deals with two transport problems. The first one is devoted to the study of diffusion due to boundary conditions in a transport problem between two plane plates. The second one consists in modelling and simulating radiative transfer phenomenon in case of the industrial context of inertial confinement fusion. (author)

Mixed hybrid finite elements and streamline computation for the potential flow problem

NARCIS (Netherlands)

Kaasschieter, E.F.; Huijben, A.J.M.

1992-01-01

An important class of problems in mathematical physics involves equations of the form -¿ · (A¿¿) = f. In a variety of problems it is desirable to obtain an accurate approximation of the flow quantity u = -A¿¿. Such an accurate approximation can be determined by the mixed finite element method. In

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

NARCIS (Netherlands)

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

2016-01-01

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

A property of assignment type mixed integer linear programming problems

NARCIS (Netherlands)

Benders, J.F.; van Nunen, J.A.E.E.

1982-01-01

In this paper we will proof that rather tight upper bounds can be given for the number of non-unique assignments that are achieved after solving the linear programming relaxation of some types of mixed integer linear assignment problems. Since in these cases the number of splitted assignments is

Modified Mixed Lagrangian-Eulerian Method Based on Numerical Framework of MT3DMS on Cauchy Boundary.

Science.gov (United States)

Suk, Heejun

2016-07-01

MT3DMS, a modular three-dimensional multispecies transport model, has long been a popular model in the groundwater field for simulating solute transport in the saturated zone. However, the method of characteristics (MOC), modified MOC (MMOC), and hybrid MOC (HMOC) included in MT3DMS did not treat Cauchy boundary conditions in a straightforward or rigorous manner, from a mathematical point of view. The MOC, MMOC, and HMOC regard the Cauchy boundary as a source condition. For the source, MOC, MMOC, and HMOC calculate the Lagrangian concentration by setting it equal to the cell concentration at an old time level. However, the above calculation is an approximate method because it does not involve backward tracking in MMOC and HMOC or allow performing forward tracking at the source cell in MOC. To circumvent this problem, a new scheme is proposed that avoids direct calculation of the Lagrangian concentration on the Cauchy boundary. The proposed method combines the numerical formulations of two different schemes, the finite element method (FEM) and the Eulerian-Lagrangian method (ELM), into one global matrix equation. This study demonstrates the limitation of all MT3DMS schemes, including MOC, MMOC, HMOC, and a third-order total-variation-diminishing (TVD) scheme under Cauchy boundary conditions. By contrast, the proposed method always shows good agreement with the exact solution, regardless of the flow conditions. Finally, the successful application of the proposed method sheds light on the possible flexibility and capability of the MT3DMS to deal with the mass transport problems of all flow regimes. © 2016, National Ground Water Association.

A Mixed Integer Linear Programming Approach to Electrical Stimulation Optimization Problems.

Science.gov (United States)

Abouelseoud, Gehan; Abouelseoud, Yasmine; Shoukry, Amin; Ismail, Nour; Mekky, Jaidaa

2018-02-01

Electrical stimulation optimization is a challenging problem. Even when a single region is targeted for excitation, the problem remains a constrained multi-objective optimization problem. The constrained nature of the problem results from safety concerns while its multi-objectives originate from the requirement that non-targeted regions should remain unaffected. In this paper, we propose a mixed integer linear programming formulation that can successfully address the challenges facing this problem. Moreover, the proposed framework can conclusively check the feasibility of the stimulation goals. This helps researchers to avoid wasting time trying to achieve goals that are impossible under a chosen stimulation setup. The superiority of the proposed framework over alternative methods is demonstrated through simulation examples.

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

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

Analysis of the type II robotic mixed-model assembly line balancing problem

Science.gov (United States)

Çil, Zeynel Abidin; Mete, Süleyman; Ağpak, Kürşad

2017-06-01

In recent years, there has been an increasing trend towards using robots in production systems. Robots are used in different areas such as packaging, transportation, loading/unloading and especially assembly lines. One important step in taking advantage of robots on the assembly line is considering them while balancing the line. On the other hand, market conditions have increased the importance of mixed-model assembly lines. Therefore, in this article, the robotic mixed-model assembly line balancing problem is studied. The aim of this study is to develop a new efficient heuristic algorithm based on beam search in order to minimize the sum of cycle times over all models. In addition, mathematical models of the problem are presented for comparison. The proposed heuristic is tested on benchmark problems and compared with the optimal solutions. The results show that the algorithm is very competitive and is a promising tool for further research.

A matheuristic for the Cargo Mix Problem with Block Stowage

DEFF Research Database (Denmark)

Christensen, Jonas Mark; Pacino, Dario

2017-01-01

The cargo-mix problem aims at selecting the amount of containers of a given type to load on a vessel. In this article we present an extended definition that includes the analysis of a circular route with draft restrictions, limitations on expected cargo and the use of a block stowage strategy...

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.

Recent Development of Multigrid Algorithms for Mixed and Noncomforming Methods for Second Order Elliptical Problems

Science.gov (United States)

Chen, Zhangxin; Ewing, Richard E.

1996-01-01

Multigrid algorithms for nonconforming and mixed finite element methods for second order elliptic problems on triangular and rectangular finite elements are considered. The construction of several coarse-to-fine intergrid transfer operators for nonconforming multigrid algorithms is discussed. The equivalence between the nonconforming and mixed finite element methods with and without projection of the coefficient of the differential problems into finite element spaces is described.

Mixed-norm estimates for the M/EEG inverse problem using accelerated gradient methods.

Science.gov (United States)

Gramfort, Alexandre; Kowalski, Matthieu; Hämäläinen, Matti

2012-04-07

Magneto- and electroencephalography (M/EEG) measure the electromagnetic fields produced by the neural electrical currents. Given a conductor model for the head, and the distribution of source currents in the brain, Maxwell's equations allow one to compute the ensuing M/EEG signals. Given the actual M/EEG measurements and the solution of this forward problem, one can localize, in space and in time, the brain regions that have produced the recorded data. However, due to the physics of the problem, the limited number of sensors compared to the number of possible source locations, and measurement noise, this inverse problem is ill-posed. Consequently, additional constraints are needed. Classical inverse solvers, often called minimum norm estimates (MNE), promote source estimates with a small ℓ₂ norm. Here, we consider a more general class of priors based on mixed norms. Such norms have the ability to structure the prior in order to incorporate some additional assumptions about the sources. We refer to such solvers as mixed-norm estimates (MxNE). In the context of M/EEG, MxNE can promote spatially focal sources with smooth temporal estimates with a two-level ℓ₁/ℓ₂ mixed-norm, while a three-level mixed-norm can be used to promote spatially non-overlapping sources between different experimental conditions. In order to efficiently solve the optimization problems of MxNE, we introduce fast first-order iterative schemes that for the ℓ₁/ℓ₂ norm give solutions in a few seconds making such a prior as convenient as the simple MNE. Furthermore, thanks to the convexity of the optimization problem, we can provide optimality conditions that guarantee global convergence. The utility of the methods is demonstrated both with simulations and experimental MEG data.

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

Directory of Open Access Journals (Sweden)

Cornelis van der Mee

2005-01-01

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

A Comprehensive Review of Boundary Integral Formulations of Acoustic Scattering Problems

Directory of Open Access Journals (Sweden)

S.I. Zaman

2000-12-01

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

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

KAUST Repository

Yan, Yan

2015-01-01

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

Mixed convective heat transfer to Sisko fluid over a radially stretching sheet in the presence of convective boundary conditions

Energy Technology Data Exchange (ETDEWEB)

Khan, Masood; Malik, Rabia, E-mail: rabiamalik.qau@gmail.com; Munir, Asif [Department of Mathematics, Quaid-i-Azam University, Islamabad 44000 (Pakistan)

2015-08-15

In this article, the mixed convective heat transfer to Sisko fluid over a radially stretching surface in the presence of convective boundary conditions is investigated. The viscous dissipation and thermal radiation effects are also taken into account. The suitable transformations are applied to convert the governing partial differential equations into a set of nonlinear coupled ordinary differential equations. The analytical solution of the governing problem is obtained by using the homotopy analysis method (HAM). Additionally, these analytical results are compared with the numerical results obtained by the shooting technique. The obtained results for the velocity and temperature are analyzed graphically for several physical parameters for the assisting and opposing flows. It is found that the effect of buoyancy parameter is more prominent in case of the assisting flow as compared to the opposing flow. Further, in tabular form the numerical values are given for the local skin friction coefficient and local Nusselt number. A remarkable agreement is noticed by comparing the present results with the results reported in the literature as a special case.

Mixed convective heat transfer to Sisko fluid over a radially stretching sheet in the presence of convective boundary conditions

International Nuclear Information System (INIS)

Khan, Masood; Malik, Rabia; Munir, Asif

2015-01-01

In this article, the mixed convective heat transfer to Sisko fluid over a radially stretching surface in the presence of convective boundary conditions is investigated. The viscous dissipation and thermal radiation effects are also taken into account. The suitable transformations are applied to convert the governing partial differential equations into a set of nonlinear coupled ordinary differential equations. The analytical solution of the governing problem is obtained by using the homotopy analysis method (HAM). Additionally, these analytical results are compared with the numerical results obtained by the shooting technique. The obtained results for the velocity and temperature are analyzed graphically for several physical parameters for the assisting and opposing flows. It is found that the effect of buoyancy parameter is more prominent in case of the assisting flow as compared to the opposing flow. Further, in tabular form the numerical values are given for the local skin friction coefficient and local Nusselt number. A remarkable agreement is noticed by comparing the present results with the results reported in the literature as a special case

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

Mixed-symmetry fields in AdS(5), conformal fields, and AdS/CFT

Energy Technology Data Exchange (ETDEWEB)

Metsaev, R.R. [Department of Theoretical Physics, P.N. Lebedev Physical Institute,Leninsky prospect 53, Moscow 119991 (Russian Federation)

2015-01-15

Mixed-symmetry arbitrary spin massive, massless, and self-dual massive fields in AdS(5) are studied. Light-cone gauge actions for such fields leading to decoupled equations of motion are constructed. Light-cone gauge formulation of mixed-symmetry anomalous conformal currents and shadows in 4d flat space is also developed. AdS/CFT correspondence for normalizable and non-normalizable modes of mixed-symmetry AdS fields and the respective boundary mixed-symmetry anomalous conformal currents and shadows is studied. We demonstrate that the light-cone gauge action for massive mixed-symmetry AdS field evaluated on solution of the Dirichlet problem amounts to the light-cone gauge 2-point vertex of mixed-symmetry anomalous shadow. Also we show that UV divergence of the action for mixed-symmetry massive AdS field with some particular value of mass parameter evaluated on the Dirichlet problem amounts to the action of long mixed-symmetry conformal field, while UV divergence of the action for mixed-symmetry massless AdS field evaluated on the Dirichlet problem amounts to the action of short mixed-symmetry conformal field. We speculate on string theory interpretation of a model which involves short low-spin conformal fields and long higher-spin conformal fields.

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

CERN Document Server

Muskhelishvili, N I

2011-01-01

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

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

Czech Academy of Sciences Publication Activity Database

Lomtatidze, Alexander

2016-01-01

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

A model study of mixing and entrainment in the horizontally evolving atmospheric convective boundary layer

Energy Technology Data Exchange (ETDEWEB)

Fedorovich, E.; Kaiser, R. [Univ. Karlsruhe, Inst. fuer Hydrologie und Wasserwirtschaft (Germany)

1997-10-01

We present results from a parallel wind-tunnel/large-eddy simulation (LES) model study of mixing and entrainment in the atmospheric convective boundary layer (CBL) longitudinally developing over a heated surface. The advection-type entrainment of warmer air from upper turbulence-free layers into the growing CBL has been investigated. Most of numerical and laboratory model studies of the CBL carried out so far dealt with another type of entrainment, namely the non-steady one, regarding the CBL growth as a non-stationary process. In the atmosphere, both types of the CBL development can take place, often being superimposed. (au)

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.

The height of the atmospheric boundary layer during unstable conditions

Energy Technology Data Exchange (ETDEWEB)

Gryning, S.E.

2005-11-01

over land, especially extensive data sets from the Lower Fraser Valley in Vancouver, Canada were useful. The performance of the model over water was examined on measurements of marine boundary layer developments over Christiansoe in the Baltic Sea. The model is found to perform relative well for both small and large distance downwind from the shoreline. The validation covers a scale of 100 km. In some practical applications meteorological data are taken as the output from operational Numerical Weather Prediction models. The mixing height is not a part of the model output but has to be determined from the available data. This aspect is dealt with in chapter 4. Methods that are used to determine the mixing height from the available data are examined and improvements suggested. Numerical Weather Prediction models have a rather coarse resolution. The issue of grid resolution and its effect on the estimated mixing height is discussed. It is found that the mixing height in the coastal zone is not resolved. In chapter 5 a new method for determining the regional fluxes of momentum and sensible heat over chessboard type patchy landscapes is presented. It is based on inversion of the model for the growth of the mixed-layer in combination with a parameterized form of the blending height concept. An example from sub-arctic Lapland illustrates that the forest dominates both the aggregated (area averaged) momentum and heat fluxes, but in quite different ways. Chapter 6 provides a description of present problems and gives an outlook for future research in the field. (au)

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

Directory of Open Access Journals (Sweden)

Zhang Peiguo

2011-01-01

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

Partial differential equations & boundary value problems with Maple

CERN Document Server

Articolo, George A

2009-01-01

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

Approximate solutions of some problems of scattering of surface ...

Indian Academy of Sciences (India)

A Choudhary

Abstract. A class of mixed boundary value problems (bvps), occurring in the study of scattering of surface water waves by thin vertical rigid barriers placed in water of finite depth, is examined for their approximate solutions. Two different placings of vertical barriers are analyzed, namely, (i) a partially immersed barrier and.

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

Czech Academy of Sciences Publication Activity Database

Mukhigulashvili, Sulkhan; Půža, B.

2015-01-01

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

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

Directory of Open Access Journals (Sweden)

Tengfei Shen

2014-02-01

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

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

Directory of Open Access Journals (Sweden)

Salih Yalcinbas

2016-01-01

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

Surface mixing and biological activity in the four Eastern Boundary Upwelling Systems

Directory of Open Access Journals (Sweden)

V. Rossi

2009-08-01

Full Text Available Eastern Boundary Upwelling Systems (EBUS are characterized by a high productivity of plankton associated with large commercial fisheries, thus playing key biological and socio-economical roles. Since they are populated by several physical oceanic structures such as filaments and eddies, which interact with the biological processes, it is a major challenge to study this sub- and mesoscale activity in connection with the chlorophyll distribution. The aim of this work is to make a comparative study of these four upwelling systems focussing on their surface stirring, using the Finite Size Lyapunov Exponents (FSLEs, and their biological activity, based on satellite data. First, the spatial distribution of horizontal mixing is analysed from time averages and from probability density functions of FSLEs, which allow us to divide each areas in two different subsystems. Then we studied the temporal variability of surface stirring focussing on the annual and seasonal cycle. We also proposed a ranking of the four EBUS based on the averaged mixing intensity. When investigating the links with chlorophyll concentration, the previous subsystems reveal distinct biological signatures. There is a global negative correlation between surface horizontal mixing and chlorophyll standing stocks over the four areas. To try to better understand this inverse relationship, we consider the vertical dimension by looking at the Ekman-transport and vertical velocities. We suggest the possibility of a changing response of the phytoplankton to sub/mesoscale turbulence, from a negative effect in the very productive coastal areas to a positive one in the open ocean. This study provides new insights for the understanding of the variable biological productivity in the ocean, which results from both dynamics of the marine ecosystem and of the 3-D turbulent medium.

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!

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.

On discrete maximum principles for nonlinear elliptic problems

Czech Academy of Sciences Publication Activity Database

Karátson, J.; Korotov, S.; Křížek, Michal

2007-01-01

Roč. 76, č. 1 (2007), s. 99-108 ISSN 0378-4754 R&D Projects: GA MŠk 1P05ME749; GA AV ČR IAA1019201 Institutional research plan: CEZ:AV0Z10190503 Keywords : nonlinear elliptic problem * mixed boundary conditions * finite element method Subject RIV: BA - General Mathematics Impact factor: 0.738, year: 2007

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

A new and self-contained presentation of the theory of boundary operators for slit diffraction and their logarithmic approximations

Energy Technology Data Exchange (ETDEWEB)

Gorenflo, Norbert [Beuth Hochschule fuer Technik Berlin (Germany). Fachbereich II; Kunik, Matthias [Magdeburg Univ. (Germany). Inst. fuer Analysis und Numerik

2009-07-01

We present a new and self-contained theory for mapping properties of the boundary operators for slit diffraction occurring in Sommerfeld's diffraction theory, covering two different cases of the polarisation of the light. This theory is entirely developed in the context of the boundary operators with a Hankel kernel and not based on the corresponding mixed boundary value problem for the Helmholtz equation. For a logarithmic approximation of the Hankel kernel we also study the corresponding mapping properties and derive explicit solutions together with certain regularity results. (orig.)

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

Directory of Open Access Journals (Sweden)

Wang Yuan-Ming

1999-01-01

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

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

Directory of Open Access Journals (Sweden)

Habib Mâagli

2014-01-01

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

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

Science.gov (United States)

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

2013-08-01

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

An Algorithm for the Mixed Transportation Network Design Problem.

Science.gov (United States)

Liu, Xinyu; Chen, Qun

2016-01-01

This paper proposes an optimization algorithm, the dimension-down iterative algorithm (DDIA), for solving a mixed transportation network design problem (MNDP), which is generally expressed as a mathematical programming with equilibrium constraint (MPEC). The upper level of the MNDP aims to optimize the network performance via both the expansion of the existing links and the addition of new candidate links, whereas the lower level is a traditional Wardrop user equilibrium (UE) problem. The idea of the proposed solution algorithm (DDIA) is to reduce the dimensions of the problem. A group of variables (discrete/continuous) is fixed to optimize another group of variables (continuous/discrete) alternately; then, the problem is transformed into solving a series of CNDPs (continuous network design problems) and DNDPs (discrete network design problems) repeatedly until the problem converges to the optimal solution. The advantage of the proposed algorithm is that its solution process is very simple and easy to apply. Numerical examples show that for the MNDP without budget constraint, the optimal solution can be found within a few iterations with DDIA. For the MNDP with budget constraint, however, the result depends on the selection of initial values, which leads to different optimal solutions (i.e., different local optimal solutions). Some thoughts are given on how to derive meaningful initial values, such as by considering the budgets of new and reconstruction projects separately.

An Algorithm for the Mixed Transportation Network Design Problem.

Directory of Open Access Journals (Sweden)

Xinyu Liu

Full Text Available This paper proposes an optimization algorithm, the dimension-down iterative algorithm (DDIA, for solving a mixed transportation network design problem (MNDP, which is generally expressed as a mathematical programming with equilibrium constraint (MPEC. The upper level of the MNDP aims to optimize the network performance via both the expansion of the existing links and the addition of new candidate links, whereas the lower level is a traditional Wardrop user equilibrium (UE problem. The idea of the proposed solution algorithm (DDIA is to reduce the dimensions of the problem. A group of variables (discrete/continuous is fixed to optimize another group of variables (continuous/discrete alternately; then, the problem is transformed into solving a series of CNDPs (continuous network design problems and DNDPs (discrete network design problems repeatedly until the problem converges to the optimal solution. The advantage of the proposed algorithm is that its solution process is very simple and easy to apply. Numerical examples show that for the MNDP without budget constraint, the optimal solution can be found within a few iterations with DDIA. For the MNDP with budget constraint, however, the result depends on the selection of initial values, which leads to different optimal solutions (i.e., different local optimal solutions. Some thoughts are given on how to derive meaningful initial values, such as by considering the budgets of new and reconstruction projects separately.

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

A mixed nonoverlapping covolume method on quadrilateral grids for elliptic problems

NARCIS (Netherlands)

Zhao, X.; Chen, Y.; Lv, J.

2016-01-01

A covolume method is proposed for the mixed formulation of second-order elliptic problems. The solution domain is divided by a quadrilateral grid, corresponding to which a nonoverlapping dual grid is constructed. The velocity and pressure are approximated by the lowest-order Raviart–Thomas space on

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 Kind of Nonlinear Programming Problem Based on Mixed Fuzzy Relation Equations Constraints

Science.gov (United States)

Li, Jinquan; Feng, Shuang; Mi, Honghai

In this work, a kind of nonlinear programming problem with non-differential objective function and under the constraints expressed by a system of mixed fuzzy relation equations is investigated. First, some properties of this kind of optimization problem are obtained. Then, a polynomial-time algorithm for this kind of optimization problem is proposed based on these properties. Furthermore, we show that this algorithm is optimal for the considered optimization problem in this paper. Finally, numerical examples are provided to illustrate our algorithms.

A quadratic approximation-based algorithm for the solution of multiparametric mixed-integer nonlinear programming problems

KAUST Repository

Domí nguez, Luis F.; Pistikopoulos, Efstratios N.

2012-01-01

An algorithm for the solution of convex multiparametric mixed-integer nonlinear programming problems arising in process engineering problems under uncertainty is introduced. The proposed algorithm iterates between a multiparametric nonlinear

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

Directory of Open Access Journals (Sweden)

Gai Gongqi

2011-01-01

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

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

Directory of Open Access Journals (Sweden)

A.M. Elizarov

2016-12-01

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

Boundary induced nonlinearities at small Reynolds numbers

NARCIS (Netherlands)

Sbragaglia, M.; Sugiyama, K.

2007-01-01

We investigate the importance of boundary slip at finite Reynolds numbers for mixed boundary conditions. Nonlinear effects are induced by the non-homogeneity of the boundary condition and change the symmetry properties of the flow with an overall mean flow reduction. To explain the observed drag

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

A hybrid Constraint Programming/Mixed Integer Programming framework for the preventive signaling maintenance crew scheduling problem

DEFF Research Database (Denmark)

Pour, Shahrzad M.; Drake, John H.; Ejlertsen, Lena Secher

2017-01-01

A railway signaling system is a complex and interdependent system which should ensure the safe operation of trains. We introduce and address a mixed integer optimisation model for the preventive signal maintenance crew scheduling problem in the Danish railway system. The problem contains many...... to feed as ‘warm start’ solutions to a Mixed Integer Programming (MIP) solver for further optimisation. We apply the CP/MIP framework to a section of the Danish rail network and benchmark our results against both direct application of a MIP solver and modelling the problem as a Constraint Optimisation...

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

Dual and primal mixed Petrov-Galerkin finite element methods in heat transfer problems

International Nuclear Information System (INIS)

Loula, A.F.D.; Toledo, E.M.

1988-12-01

New mixed finite element formulations for the steady state heat transfer problem are presented with no limitation in the choice of conforming finite element spaces. Adding least square residual forms of the governing equations of the classical Galerkin formulation the original saddle point problem is transformed into a minimization problem. Stability analysis, error estimates and numerical results are presented, confirming the error estimates and the good performance of this new formulation. (author) [pt

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.

Visual, Algebraic and Mixed Strategies in Visually Presented Linear Programming Problems.

Science.gov (United States)

Shama, Gilli; Dreyfus, Tommy

1994-01-01

Identified and classified solution strategies of (n=49) 10th-grade students who were presented with linear programming problems in a predominantly visual setting in the form of a computerized game. Visual strategies were developed more frequently than either algebraic or mixed strategies. Appendix includes questionnaires. (Contains 11 references.)…

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

Directory of Open Access Journals (Sweden)

O.Kh. Abdullaev

2014-06-01

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

Periodic Contact and Crack Problems in Plane Elasticity

DEFF Research Database (Denmark)

Krenk, Steen

1976-01-01

By use of singular integral equations it is demonstrated how some periodic contact and crack problems can be solved in closed form. The integral equation in question is the same which is encountered when dealing with mixed boundary conditions on a circle. As analytical evaluation of the solution ...... may be quite complicated attention is drawn to a numerical quadrature method...

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.

Domain Immersion Technique And Free Surface Computations Applied To Extrusion And Mixing Processes

Science.gov (United States)

Valette, Rudy; Vergnes, Bruno; Basset, Olivier; Coupez, Thierry

2007-04-01

This work focuses on the development of numerical techniques devoted to the simulation of mixing processes of complex fluids such as twin-screw extrusion or batch mixing. In mixing process simulation, the absence of symmetry of the moving boundaries (the screws or the rotors) implies that their rigid body motion has to be taken into account by using a special treatment. We therefore use a mesh immersion technique (MIT), which consists in using a P1+/P1-based (MINI-element) mixed finite element method for solving the velocity-pressure problem and then solving the problem in the whole barrel cavity by imposing a rigid motion (rotation) to nodes found located inside the so called immersed domain, each subdomain (screw, rotor) being represented by a surface CAD mesh (or its mathematical equation in simple cases). The independent meshes are immersed into a unique backgound computational mesh by computing the distance function to their boundaries. Intersections of meshes are accounted for, allowing to compute a fill factor usable as for the VOF methodology. This technique, combined with the use of parallel computing, allows to compute the time-dependent flow of generalized Newtonian fluids including yield stress fluids in a complex system such as a twin screw extruder, including moving free surfaces, which are treated by a "level set" and Hamilton-Jacobi method.

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

DEFF Research Database (Denmark)

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

2014-01-01

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

A multigrid solution method for mixed hybrid finite elements

Energy Technology Data Exchange (ETDEWEB)

Schmid, W. [Universitaet Augsburg (Germany)

1996-12-31

We consider the multigrid solution of linear equations arising within the discretization of elliptic second order boundary value problems of the form by mixed hybrid finite elements. Using the equivalence of mixed hybrid finite elements and non-conforming nodal finite elements, we construct a multigrid scheme for the corresponding non-conforming finite elements, and, by this equivalence, for the mixed hybrid finite elements, following guidelines from Arbogast/Chen. For a rectangular triangulation of the computational domain, this non-conforming schemes are the so-called nodal finite elements. We explicitly construct prolongation and restriction operators for this type of non-conforming finite elements. We discuss the use of plain multigrid and the multilevel-preconditioned cg-method and compare their efficiency in numerical tests.

A quadratic approximation-based algorithm for the solution of multiparametric mixed-integer nonlinear programming problems

KAUST Repository

Domínguez, Luis F.

2012-06-25

An algorithm for the solution of convex multiparametric mixed-integer nonlinear programming problems arising in process engineering problems under uncertainty is introduced. The proposed algorithm iterates between a multiparametric nonlinear programming subproblem and a mixed-integer nonlinear programming subproblem to provide a series of parametric upper and lower bounds. The primal subproblem is formulated by fixing the integer variables and solved through a series of multiparametric quadratic programming (mp-QP) problems based on quadratic approximations of the objective function, while the deterministic master subproblem is formulated so as to provide feasible integer solutions for the next primal subproblem. To reduce the computational effort when infeasibilities are encountered at the vertices of the critical regions (CRs) generated by the primal subproblem, a simplicial approximation approach is used to obtain CRs that are feasible at each of their vertices. The algorithm terminates when there does not exist an integer solution that is better than the one previously used by the primal problem. Through a series of examples, the proposed algorithm is compared with a multiparametric mixed-integer outer approximation (mp-MIOA) algorithm to demonstrate its computational advantages. © 2012 American Institute of Chemical Engineers (AIChE).

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

International Nuclear Information System (INIS)

Xu Shihe

2009-01-01

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

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

Directory of Open Access Journals (Sweden)

Long Yuhua

2017-12-01

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

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

Directory of Open Access Journals (Sweden)

Liu Yang

2007-10-01

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

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

Directory of Open Access Journals (Sweden)

Domoshnitsky Alexander

2009-01-01

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

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

Directory of Open Access Journals (Sweden)

Chuanzhi Bai

2010-06-01

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

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

International Nuclear Information System (INIS)

Frittelli, Simonetta; Gomez, Roberto

2007-01-01

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

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

Czech Academy of Sciences Publication Activity Database

Kiguradze, I.; Šremr, Jiří

2011-01-01

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

Second-order polynomial model to solve the least-cost lumber grade mix problem

Science.gov (United States)

Urs Buehlmann; Xiaoqiu Zuo; R. Edward. Thomas

2010-01-01

Material costs when cutting solid wood parts from hardwood lumber for secondary wood products manufacturing account for 20 to 50 percent of final product cost. These costs can be minimized by proper selection of the lumber quality used. The lumber quality selection problem is referred to as the least-cost lumber grade mix problem in the industry. The objective of this...

A General Iterative Method of Fixed Points for Mixed Equilibrium Problems and Variational Inclusion Problems

Directory of Open Access Journals (Sweden)

Phayap Katchang

2010-01-01

Full Text Available The purpose of this paper is to investigate the problem of finding a common element of the set of solutions for mixed equilibrium problems, the set of solutions of the variational inclusions with set-valued maximal monotone mappings and inverse-strongly monotone mappings, and the set of fixed points of a family of finitely nonexpansive mappings in the setting of Hilbert spaces. We propose a new iterative scheme for finding the common element of the above three sets. Our results improve and extend the corresponding results of the works by Zhang et al. (2008, Peng et al. (2008, Peng and Yao (2009, as well as Plubtieng and Sriprad (2009 and some well-known results in the literature.

Problems of mixed convection flow regime map in a vertical cylinder

International Nuclear Information System (INIS)

Kang, Gyeong Uk; Chung, Bum Jin

2012-01-01

One of the technical issues by the development of the VHTR is the mixed convection, which is the regime of heat transfer that occurs when the driving forces of both forced and natural convection are of comparable orders of magnitude. In vertical internal flows, the buoyancy force acts upward only, but forced flows can move either upward or downward. Thus, there are two types of mixed convection flows, depending on the direction of the forced flow. When the directions of the forced flow and buoyancy are the same, the flow is a buoyancy aided flow; when they are opposite, the flow is a buoyancy opposed flow. In laminar flows, buoyancy aided flow shows enhanced heat transfer compared to the pure forced convection and buoyancy opposed flow shows impaired heat transfer due to the flow velocity affected by the buoyancy forces. In turbulent flows, however, buoyancy opposed flows shows enhanced heat transfer due to increased turbulence production and buoyancy aided flow shows impaired heat transfer at low buoyancy forces and as the buoyancy increases, the heat transfer restores and at further increases of the buoyancy forces, the heat transfer is enhanced. It is of primary interests to classify which convection regime is mainly dominant. The methods most used to classify between forced, mixed and natural convection have been to refer to the classical flow regime map suggested by Meta is and Eckert. During the course of fundamental literature studies on this topic, it is found that there are some problems on the flow regime map in a vertical cylinder. This paper is to discuss problems identified through reviewing the papers composed in the classical flow regime map. We have tried to reproduce the flow regime map independently using the data obtained from the literatures and compared with the classical flow regime map and finally, the problems on this topic were discussed

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

Energy Technology Data Exchange (ETDEWEB)

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

2008-12-09

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

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

International Nuclear Information System (INIS)

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

1996-01-01

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

Mixed integer linear programming model for dynamic supplier selection problem considering discounts

Directory of Open Access Journals (Sweden)

Adi Wicaksono Purnawan

2018-01-01

Full Text Available Supplier selection is one of the most important elements in supply chain management. This function involves evaluation of many factors such as, material costs, transportation costs, quality, delays, supplier capacity, storage capacity and others. Each of these factors varies with time, therefore, supplier identified for one period is not necessarily be same for the next period to supply the same product. So, mixed integer linear programming (MILP was developed to overcome the dynamic supplier selection problem (DSSP. In this paper, a mixed integer linear programming model is built to solve the lot-sizing problem with multiple suppliers, multiple periods, multiple products and quantity discounts. The buyer has to make a decision for some products which will be supplied by some suppliers for some periods cosidering by discount. To validate the MILP model with randomly generated data. The model is solved by Lingo 16.

Diffusive boundary layers over varying topography

KAUST Repository

Dell, R.  W.; Pratt, L.  J.

2015-01-01

Diffusive bottom boundary layers can produce upslope flows in a stratified fluid. Accumulating observations suggest that these boundary layers may drive upwelling and mixing in mid-ocean ridge flank canyons. However, most studies of diffusive bottom

Half-Space Temperature Field with a Movable Thermally Thin-Coated Boundary Under External Heat Flux

Directory of Open Access Journals (Sweden)

P. A. Vlasov

2014-01-01

Full Text Available In engineering practice analytical methods of the mathematical theory of heat conduction hold a special place. This is due to many reasons, in particular, because of the fact that the solutions of the relevant problems represented in analytically closed form, can be used not only for a parametric analysis of the studied temperature field and to explore the specific features of its formation, but also to test the developed computational algorithms, which are aimed at solving real-world application heat and mass transfer problems. Difficulties arising when using the analytical mathematical theory methods of heat conduction in practice are well known. Also they are significantly exacerbated if the boundaries of the system under study are movable, even in the simplest case, when the law of motion is known.The main goal of the conducted research is to have an analytically closed-form problem solution for finding the orthotropic half-space temperature field, a boundary of which has thermally thin coating exposed to extremely concentrated stationary external heat flux and uniformly moves parallel to itself.The assumption that the covering of the boundary is thermally thin, allowed to realize the idea of \\concentrated capacity", that is to accept the hypothesis that the mean-thickness coating temperature is equal to the temperature of its boundaries. This assumption allowed us to reduce the problem under consideration to a mixed problem for a parabolic equation with a specific boundary condition.The Hankel integral transform of zero order with respect to the radial variable and the Laplace transform with respect to the temporal variable were used to solve the reduced problem. These techniques have allowed us to submit the required solution as an iterated integral.

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

Directory of Open Access Journals (Sweden)

Khaleghi Moghadam Mohsen

2017-08-01

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

Boundary fluxes for nonlocal diffusion

Science.gov (United States)

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

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

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

Directory of Open Access Journals (Sweden)

Mabrouk Briki

2016-05-01

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

A Bee Colony Optimization Approach for Mixed Blocking Constraints Flow Shop Scheduling Problems

Directory of Open Access Journals (Sweden)

Mostafa Khorramizadeh

2015-01-01

Full Text Available The flow shop scheduling problems with mixed blocking constraints with minimization of makespan are investigated. The Taguchi orthogonal arrays and path relinking along with some efficient local search methods are used to develop a metaheuristic algorithm based on bee colony optimization. In order to compare the performance of the proposed algorithm, two well-known test problems are considered. Computational results show that the presented algorithm has comparative performance with well-known algorithms of the literature, especially for the large sized problems.

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

Directory of Open Access Journals (Sweden)

Diana-Raluca Herlea

2016-02-01

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

A symmetric solution of a multipoint boundary value problem at resonance

Directory of Open Access Journals (Sweden)

2006-01-01

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

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

Science.gov (United States)

Cai, Wei; Wang, Jian-Zhong

1993-01-01

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

How can we describe the entrainment processes in sheared convective boundary layers?: a large-eddy simulation and mixed-layer theory/model comparison study

NARCIS (Netherlands)

Pino, D.; Vilà-Guerau de Arellano, J.; Kim, S.W.

2006-01-01

Dry convective boundary layers characterized by a significant wind shear on the surface and at the inversion zone are studied by means of the mixed layer theory. Two different representations of the entrainment zone, each of which has a different closure of the entrainment heat flux, are considered.

The Solution Set Characterization and Error Bound for the Extended Mixed Linear Complementarity Problem

Directory of Open Access Journals (Sweden)

Hongchun Sun

2012-01-01

Full Text Available For the extended mixed linear complementarity problem (EML CP, we first present the characterization of the solution set for the EMLCP. Based on this, its global error bound is also established under milder conditions. The results obtained in this paper can be taken as an extension for the classical linear complementarity problems.

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

Directory of Open Access Journals (Sweden)

Jian Liu

2013-09-01

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

Nonlinear second-order multivalued boundary value problems

Indian Academy of Sciences (India)

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

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

Computational model for turbulent flow around a grid spacer with mixing vane

International Nuclear Information System (INIS)

Tsutomu Ikeno; Takeo Kajishima

2005-01-01

Turbulent mixing coefficient and pressure drop are important factors in subchannel analysis to predict onset of DNB. However, universal correlations are difficult since these factors are significantly affected by the geometry of subchannel and a grid spacer with mixing vane. Therefore, we propose a computational model to estimate these factors. Computational model: To represent the effect of geometry of grid spacer in computational model, we applied a large eddy simulation (LES) technique in couple with an improved immersed-boundary method. In our previous work (Ikeno, et al., NURETH-10), detailed properties of turbulence in subchannel were successfully investigated by developing the immersed boundary method in LES. In this study, additional improvements are given: new one-equation dynamic sub-grid scale (SGS) model is introduced to account for the complex geometry without any artificial modification; the higher order accuracy is maintained by consistent treatment for boundary conditions for velocity and pressure. NUMERICAL TEST AND DISCUSSION: Turbulent mixing coefficient and pressure drop are affected strongly by the arrangement and inclination of mixing vane. Therefore, computations are carried out for each of convolute and periodic arrangements, and for each of 30 degree and 20 degree inclinations. The difference in turbulent mixing coefficient due to these factors is reasonably predicted by our method. (An example of this numerical test is shown in Fig. 1.) Turbulent flow of the problem includes unsteady separation behind the mixing vane and vortex shedding in downstream. Anisotropic distribution of turbulent stress is also appeared in rod gap. Therefore, our computational model has advantage for assessing the influence of arrangement and inclination of mixing vane. By coarser computational mesh, one can screen several candidates for spacer design. Then, by finer mesh, more quantitative analysis is possible. By such a scheme, we believe this method is useful

A Mixed Traffic Assignment Problem and Its Upper Bound of Efficiency Loss

Directory of Open Access Journals (Sweden)

Xiao-Jun Yu

2014-01-01

Full Text Available A variational inequality (VI formulation of a mixed traffic assignment problem associated with user equilibrium (UE player and altruistic players (UE-AP is developed. The users governed by the UE player use the paths with minimum actual travel cost, while the users governed by altruistic players use the paths with minimum perceived travel cost, which is a linear combination of the actual and marginal travel costs. A nonlinear programming method is adopted to derive the efficiency loss of the UE-AP mixed equilibrium with polynomial cost function. Finally, numerical examples are given to verify our analytical results.

Domain decomposition methods for the mixed dual formulation of the critical neutron diffusion problem

International Nuclear Information System (INIS)

Guerin, P.

2007-12-01

The neutronic simulation of a nuclear reactor core is performed using the neutron transport equation, and leads to an eigenvalue problem in the steady-state case. Among the deterministic resolution methods, diffusion approximation is often used. For this problem, the MINOS solver based on a mixed dual finite element method has shown his efficiency. In order to take advantage of parallel computers, and to reduce the computing time and the local memory requirement, we propose in this dissertation two domain decomposition methods for the resolution of the mixed dual form of the eigenvalue neutron diffusion problem. The first approach is a component mode synthesis method on overlapping sub-domains. Several Eigenmodes solutions of a local problem solved by MINOS on each sub-domain are taken as basis functions used for the resolution of the global problem on the whole domain. The second approach is a modified iterative Schwarz algorithm based on non-overlapping domain decomposition with Robin interface conditions. At each iteration, the problem is solved on each sub domain by MINOS with the interface conditions deduced from the solutions on the adjacent sub-domains at the previous iteration. The iterations allow the simultaneous convergence of the domain decomposition and the eigenvalue problem. We demonstrate the accuracy and the efficiency in parallel of these two methods with numerical results for the diffusion model on realistic 2- and 3-dimensional cores. (author)