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)

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

Classically integrable boundary conditions for symmetric-space sigma models

International Nuclear Information System (INIS)

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

2004-01-01

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

Existence and Analytic Approximation of Solutions of Duffing Type Nonlinear Integro-Differential Equation with Integral Boundary Conditions

Directory of Open Access Journals (Sweden)

Alsaedi Ahmed

2009-01-01

Full Text Available A generalized quasilinearization technique is developed to obtain a sequence of approximate solutions converging monotonically and quadratically to a unique solution of a boundary value problem involving Duffing type nonlinear integro-differential equation with integral boundary conditions. The convergence of order for the sequence of iterates is also established. It is found that the work presented in this paper not only produces new results but also yields several old results in certain limits.

Boundary Controllability of Nonlinear Fractional Integrodifferential Systems

Directory of Open Access Journals (Sweden)

Ahmed HamdyM

2010-01-01

Full Text Available Sufficient conditions for boundary controllability of nonlinear fractional integrodifferential systems in Banach space are established. The results are obtained by using fixed point theorems. We also give an application for integropartial differential equations of fractional order.

Breather type solutions of the vector nonlinear Schroedinger equation with quasi-constant boundary conditions

International Nuclear Information System (INIS)

Makhan'kov, V.G.; Slavov, S.I.

1989-01-01

Vector nonlinear Schroedinger equations (VS3) is investigated under quasi-constant boundary conditions. New two-soliton solutions are obtained with such non-trivial dynamics that they may be called the breather solutions. A version of the basic Novikov-Dubrovin-Krichever algebro-geometrical approach is applied to obtain breather like solutions existing for all types of internal symmetry is specified are formulated in terms of the soliton velocity expressed via the parameters of the problem. 4 refs

Electromagnetic pulses at the boundary of a nonlinear plasma

International Nuclear Information System (INIS)

Satorius, E.H.

1975-01-01

An investigation was made of the behavior of strong electromagnetic pulses at the boundary of a nonlinear, cold, collisionless, and uniform plasma. The nonlinearity considered here is due to the nonlinear terms in the fluid equation which is used to describe the plasma. Two cases are studied. First, the case where there is a voltage pulse applied across the plane boundary of a semi-infinite, nonlinear plasma. Two different voltage pulses are considered, i.e., a delta function pulse and a suddenly turned-on sinusoidal pulse. The resulting electromagnetic fields propagating in the nonlinear plasma are found in this case. In the second case, the reflection of incident E-polarized and H-polarized, electromagnetic pulses at various angles of incidence from a nonlinear, semi-infinite plasma are considered. Again, two forms of incident pulses are considered: a delta function pulse and a suddenly turned-on sinusoidal pulse. In case two, the reflected electromagnetic fields are found. In both cases, the method used for finding the fields is to first solve the fluid equation (which describes the plasma) for the nonlinear conduction current in terms of the electric field using a perturbation method (since the nonlinear effects are assumed to be small). Next, this current is substituted into Maxwell's equations, and finally the electromagnetic fields which satisfy the boundary conditions are found. (U.S.)

Entropy Stable Wall Boundary Conditions for the Compressible Navier-Stokes Equations

Science.gov (United States)

Parsani, Matteo; Carpenter, Mark H.; Nielsen, Eric J.

2014-01-01

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite volume, finite difference, discontinuous Galerkin, and flux reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

Discretely Conservative Finite-Difference Formulations for Nonlinear Conservation Laws in Split Form: Theory and Boundary Conditions

Science.gov (United States)

Fisher, Travis C.; Carpenter, Mark H.; Nordstroem, Jan; Yamaleev, Nail K.; Swanson, R. Charles

2011-01-01

Simulations of nonlinear conservation laws that admit discontinuous solutions are typically restricted to discretizations of equations that are explicitly written in divergence form. This restriction is, however, unnecessary. Herein, linear combinations of divergence and product rule forms that have been discretized using diagonal-norm skew-symmetric summation-by-parts (SBP) operators, are shown to satisfy the sufficient conditions of the Lax-Wendroff theorem and thus are appropriate for simulations of discontinuous physical phenomena. Furthermore, special treatments are not required at the points that are near physical boundaries (i.e., discrete conservation is achieved throughout the entire computational domain, including the boundaries). Examples are presented of a fourth-order, SBP finite-difference operator with second-order boundary closures. Sixth- and eighth-order constructions are derived, and included in E. Narrow-stencil difference operators for linear viscous terms are also derived; these guarantee the conservative form of the combined operator.

Effect of plate permeability on nonlinear stability of the asymptotic suction boundary layer.

Science.gov (United States)

Wedin, Håkan; Cherubini, Stefania; Bottaro, Alessandro

2015-07-01

The nonlinear stability of the asymptotic suction boundary layer is studied numerically, searching for finite-amplitude solutions that bifurcate from the laminar flow state. By changing the boundary conditions for disturbances at the plate from the classical no-slip condition to more physically sound ones, the stability characteristics of the flow may change radically, both for the linearized as well as the nonlinear problem. The wall boundary condition takes into account the permeability K̂ of the plate; for very low permeability, it is acceptable to impose the classical boundary condition (K̂=0). This leads to a Reynolds number of approximately Re(c)=54400 for the onset of linearly unstable waves, and close to Re(g)=3200 for the emergence of nonlinear solutions [F. A. Milinazzo and P. G. Saffman, J. Fluid Mech. 160, 281 (1985); J. H. M. Fransson, Ph.D. thesis, Royal Institute of Technology, KTH, Sweden, 2003]. However, for larger values of the plate's permeability, the lower limit for the existence of linear and nonlinear solutions shifts to significantly lower Reynolds numbers. For the largest permeability studied here, the limit values of the Reynolds numbers reduce down to Re(c)=796 and Re(g)=294. For all cases studied, the solutions bifurcate subcritically toward lower Re, and this leads to the conjecture that they may be involved in the very first stages of a transition scenario similar to the classical route of the Blasius boundary layer initiated by Tollmien-Schlichting (TS) waves. The stability of these nonlinear solutions is also investigated, showing a low-frequency main unstable mode whose growth rate decreases with increasing permeability and with the Reynolds number, following a power law Re(-ρ), where the value of ρ depends on the permeability coefficient K̂. The nonlinear dynamics of the flow in the vicinity of the computed finite-amplitude solutions is finally investigated by direct numerical simulations, providing a viable scenario for

Nonlinear Seismic Behavior of Different Boundary Conditions of Transmission Line Systems under Earthquake Loading

Directory of Open Access Journals (Sweden)

Li Tian

2016-01-01

Full Text Available Nonlinear seismic behaviors of different boundary conditions of transmission line system under earthquake loading are investigated in this paper. The transmission lines are modeled by cable element accounting for the nonlinearity of the cable. For the suspension type, three towers and two span lines with spring model (Model 1 and three towers and four span lines’ model (Model 2 are established, respectively. For the tension type, three towers and two span lines’ model (Model 3 and three towers and four span lines’ model (Model 4 are created, respectively. The frequencies of the transmission towers and transmission lines of the suspension type and tension type are calculated, respectively. The responses of the suspension type and tension type are investigated using nonlinear time history analysis method, respectively. The results show that the responses of the transmission tower and transmission line of the two models of the suspension type are slightly different. However, the responses of transmission tower and transmission line of the two models of the tension type are significantly different. Therefore, in order to obtain accurate results, a reasonable model should be considered. The results could provide a reference for the seismic analysis of the transmission tower-line system.

Dynamic relaxation method for nonlinear buckling analysis of moderately thick FG cylindrical panels with various boundary conditions

Energy Technology Data Exchange (ETDEWEB)

Golmakani, M. E.; Far, M. N. Sadraee; Moravej, M. [Dept. of Mechanical Engineering, Mashhad Branch, Islamic Azad University, Mashhad (Iran, Islamic Republic of)

2016-12-15

Using new approach proposed by Dynamic relaxation (DR) method, buckling analysis of moderately thick Functionally graded (FG) cylindrical panels subjected to axial compression is investigated for various boundary conditions. The mechanical properties of FG panel are assumed to vary continuously along the thickness direction by the simple rule of mixture and Mori-Tanaka model. The incremental form of nonlinear formulations are derived based on First-order shear deformation theory (FSDT) and large deflection von Karman equations. The DR method combined with the finite difference discretization technique is employed to solve the incremental form of equilibrium equations. The critical mechanical buckling load is determined based on compressive load-displacement curve by adding the incremental displacements in each load step to the displacements obtained from the previous ones. A detailed parametric study is carried out to investigate the influences of the boundary conditions, rule of mixture, grading index, radius-to-thickness ratio, length-to-radius ratio and panel angle on the mechanical buckling load. The results reveal that with increase of grading index the effect of radius-to-thickness ratio on the buckling load decreases. It is also observed that effect of distribution rules on the buckling load is dependent to the type of boundary conditions.

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

Directory of Open Access Journals (Sweden)

Guotao Wang

2012-01-01

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

On the wave equation with semilinear porous acoustic boundary conditions

KAUST Repository

Graber, Philip Jameson; Said-Houari, Belkacem

2012-01-01

The goal of this work is to study a model of the wave equation with semilinear porous acoustic boundary conditions with nonlinear boundary/interior sources and a nonlinear boundary/interior damping. First, applying the nonlinear semigroup theory, we show the existence and uniqueness of local in time solutions. The main difficulty in proving the local existence result is that the Neumann boundary conditions experience loss of regularity due to boundary sources. Using an approximation method involving truncated sources and adapting the ideas in Lasiecka and Tataru (1993) [28], we show that the existence of solutions can still be obtained. Second, we prove that under some restrictions on the source terms, then the local solution can be extended to be global in time. In addition, it has been shown that the decay rates of the solution are given implicitly as solutions to a first order ODE and depends on the behavior of the damping terms. In several situations, the obtained ODE can be easily solved and the decay rates can be given explicitly. Third, we show that under some restrictions on the initial data and if the interior source dominates the interior damping term and if the boundary source dominates the boundary damping, then the solution ceases to exists and blows up in finite time. Moreover, in either the absence of the interior source or the boundary source, then we prove that the solution is unbounded and grows as an exponential function. © 2012 Elsevier Inc.

On the wave equation with semilinear porous acoustic boundary conditions

KAUST Repository

Graber, Philip Jameson

2012-05-01

The goal of this work is to study a model of the wave equation with semilinear porous acoustic boundary conditions with nonlinear boundary/interior sources and a nonlinear boundary/interior damping. First, applying the nonlinear semigroup theory, we show the existence and uniqueness of local in time solutions. The main difficulty in proving the local existence result is that the Neumann boundary conditions experience loss of regularity due to boundary sources. Using an approximation method involving truncated sources and adapting the ideas in Lasiecka and Tataru (1993) [28], we show that the existence of solutions can still be obtained. Second, we prove that under some restrictions on the source terms, then the local solution can be extended to be global in time. In addition, it has been shown that the decay rates of the solution are given implicitly as solutions to a first order ODE and depends on the behavior of the damping terms. In several situations, the obtained ODE can be easily solved and the decay rates can be given explicitly. Third, we show that under some restrictions on the initial data and if the interior source dominates the interior damping term and if the boundary source dominates the boundary damping, then the solution ceases to exists and blows up in finite time. Moreover, in either the absence of the interior source or the boundary source, then we prove that the solution is unbounded and grows as an exponential function. © 2012 Elsevier Inc.

Entropy Stable Wall Boundary Conditions for the Three-Dimensional Compressible Navier-Stokes Equations

Science.gov (United States)

Parsani, Matteo; Carpenter, Mark H.; Nielsen, Eric J.

2015-01-01

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the three-dimensional compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators on unstructured grids are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite difference, finite volume, discontinuous Galerkin, and flux reconstruction/correction procedure via reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

Existence and asymptotic behavior of the wave equation with dynamic boundary conditions

KAUST Repository

Graber, Philip Jameson; Said-Houari, Belkacem

2012-01-01

The goal of this work is to study a model of the strongly damped wave equation with dynamic boundary conditions and nonlinear boundary/interior sources and nonlinear boundary/interior damping. First, applying the nonlinear semigroup theory, we show the existence and uniqueness of local in time solutions. In addition, we show that in the strongly damped case solutions gain additional regularity for positive times t>0. Second, we show that under some restrictions on the initial data and if the interior source dominates the interior damping term and if the boundary source dominates the boundary damping, then the solution grows as an exponential function. Moreover, in the absence of the strong damping term, we prove that the solution ceases to exists and blows up in finite time. © 2012 Springer Science+Business Media, LLC.

Existence and asymptotic behavior of the wave equation with dynamic boundary conditions

KAUST Repository

Graber, Philip Jameson

2012-03-07

The goal of this work is to study a model of the strongly damped wave equation with dynamic boundary conditions and nonlinear boundary/interior sources and nonlinear boundary/interior damping. First, applying the nonlinear semigroup theory, we show the existence and uniqueness of local in time solutions. In addition, we show that in the strongly damped case solutions gain additional regularity for positive times t>0. Second, we show that under some restrictions on the initial data and if the interior source dominates the interior damping term and if the boundary source dominates the boundary damping, then the solution grows as an exponential function. Moreover, in the absence of the strong damping term, we prove that the solution ceases to exists and blows up in finite time. © 2012 Springer Science+Business Media, LLC.

Sum-frequency nonlinear Cherenkov radiation generated on the boundary of bulk medium crystal.

Science.gov (United States)

Wang, Xiaojing; Cao, Jianjun; Zhao, Xiaohui; Zheng, Yuanlin; Ren, Huaijin; Deng, Xuewei; Chen, Xianfeng

2015-12-14

We demonstrated experimentally a method to generate the sum-frequency Nonlinear Cherenkov radiation (NCR) on the boundary of bulk medium by using two synchronized laser beam with wavelength of 1300 nm and 800 nm. It is also an evidence that the polarization wave is always confined to the boundary. Critical conditions of surface sum-frequency NCR under normal and anomalous dispersion condition is discussed.

Nonlinear streak computation using boundary region equations

Energy Technology Data Exchange (ETDEWEB)

Martin, J A; Martel, C, E-mail: juanangel.martin@upm.es, E-mail: carlos.martel@upm.es [Depto. de Fundamentos Matematicos, E.T.S.I Aeronauticos, Universidad Politecnica de Madrid, Plaza Cardenal Cisneros 3, 28040 Madrid (Spain)

2012-08-01

The boundary region equations (BREs) are applied for the simulation of the nonlinear evolution of a spanwise periodic array of streaks in a flat plate boundary layer. The well-known BRE formulation is obtained from the complete Navier-Stokes equations in the high Reynolds number limit, and provides the correct asymptotic description of three-dimensional boundary layer streaks. In this paper, a fast and robust streamwise marching scheme is introduced to perform their numerical integration. Typical streak computations present in the literature correspond to linear streaks or to small-amplitude nonlinear streaks computed using direct numerical simulation (DNS) or the nonlinear parabolized stability equations (PSEs). We use the BREs to numerically compute high-amplitude streaks, a method which requires much lower computational effort than DNS and does not have the consistency and convergence problems of the PSE. It is found that the flow configuration changes substantially as the amplitude of the streaks grows and the nonlinear effects come into play. The transversal motion (in the wall normal-streamwise plane) becomes more important and strongly distorts the streamwise velocity profiles, which end up being quite different from those of the linear case. We analyze in detail the resulting flow patterns for the nonlinearly saturated streaks and compare them with available experimental results. (paper)

Nonlinear Transient Growth and Boundary Layer Transition

Science.gov (United States)

Paredes, Pedro; Choudhari, Meelan M.; Li, Fei

2016-01-01

Parabolized stability equations (PSE) are used in a variational approach to study the optimal, non-modal disturbance growth in a Mach 3 at plate boundary layer and a Mach 6 circular cone boundary layer. As noted in previous works, the optimal initial disturbances correspond to steady counter-rotating streamwise vortices, which subsequently lead to the formation of streamwise-elongated structures, i.e., streaks, via a lift-up effect. The nonlinear evolution of the linearly optimal stationary perturbations is computed using the nonlinear plane-marching PSE for stationary perturbations. A fully implicit marching technique is used to facilitate the computation of nonlinear streaks with large amplitudes. To assess the effect of the finite-amplitude streaks on transition, the linear form of plane- marching PSE is used to investigate the instability of the boundary layer flow modified by spanwise periodic streaks. The onset of bypass transition is estimated by using an N- factor criterion based on the amplification of the streak instabilities. Results show that, for both flow configurations of interest, streaks of sufficiently large amplitude can lead to significantly earlier onset of transition than that in an unperturbed boundary layer without any streaks.

Analysis of boundary layer flow over a porous nonlinearly stretching sheet with partial slip at

Directory of Open Access Journals (Sweden)

Swati Mukhopadhyay

2013-12-01

Full Text Available The boundary layer flow of a viscous incompressible fluid toward a porous nonlinearly stretching sheet is considered in this analysis. Velocity slip is considered instead of no-slip condition at the boundary. Similarity transformations are used to convert the partial differential equation corresponding to the momentum equation into nonlinear ordinary differential equation. Numerical solution of this equation is obtained by shooting method. It is found that the horizontal velocity decreases with increasing slip parameter.

Nonlinear radiative heat transfer in magnetohydrodynamic (MHD stagnation point flow of nanofluid past a stretching sheet with convective boundary condition

Directory of Open Access Journals (Sweden)

Wubshet Ibrahim

2015-12-01

Full Text Available Two-dimensional boundary layer flow of nanofluid fluid past a stretching sheet is examined. The paper reveals the effect of non-linear radiative heat transfer on magnetohydrodynamic (MHD stagnation point flow past a stretching sheet with convective heating. Condition of zero normal flux of nanoparticles at the wall for the stretched flow is considered. The nanoparticle fractions on the boundary are considered to be passively controlled. The solution for the velocity, temperature and nanoparticle concentration depends on parameters viz. Prandtl number Pr, velocity ratio parameter A, magnetic parameter M, Lewis number Le, Brownian motion Nb, and the thermophoresis parameter Nt. Moreover, the problem is governed by temperature ratio parameter (Nr=TfT∞ and radiation parameter Rd. Similarity transformation is used to reduce the governing non-linear boundary-value problems into coupled higher order non-linear ordinary differential equation. These equations were numerically solved using the function bvp4c from the matlab software for different values of governing parameters. Numerical results are obtained for velocity, temperature and concentration, as well as the skin friction coefficient and local Nusselt number. The results indicate that the skin friction coefficient Cf increases as the values of magnetic parameter M increase and decreases as the values of velocity ratio parameter A increase. The local Nusselt number −θ′(0 decreases as the values of thermophoresis parameter Nt and radiation parameter Nr increase and it increases as the values of both Biot number Bi and Prandtl number Pr increase. Furthermore, radiation has a positive effect on temperature and concentration profiles.

Lyapunov Based Estimation of Flight Stability Boundary under Icing Conditions

Directory of Open Access Journals (Sweden)

Binbin Pei

2017-01-01

Full Text Available Current fight boundary of the envelope protection in icing conditions is usually defined by the critical values of state parameters; however, such method does not take the interrelationship of each parameter and the effect of the external disturbance into consideration. This paper proposes constructing the stability boundary of the aircraft in icing conditions through analyzing the region of attraction (ROA around the equilibrium point. Nonlinear icing effect model is proposed according to existing wind tunnel test results. On this basis, the iced polynomial short period model can be deduced further to obtain the stability boundary under icing conditions using ROA analysis. Simulation results for a series of icing severity demonstrate that, regardless of the icing severity, the boundary of the calculated ROA can be treated as an estimation of the stability boundary around an equilibrium point. The proposed methodology is believed to be a promising way for ROA analysis and stability boundary construction of the aircraft in icing conditions, and it will provide theoretical support for multiple boundary protection of icing tolerant flight.

Boundary condition effect on response modification factor of X-braced steel frames

Directory of Open Access Journals (Sweden)

Walid A. Attia

2018-04-01

Full Text Available Design of the structures to resist seismic force depends on the theory of dissipation in elastic energy that already exists in response modification factor “R-factor”. The main problem in codes gives a constant value for R-factor, since change in boundary conditions of building change in behavior of braced steel frame structures and that effects on R-factor. This study is an attempt to assess overstrength, ductility and response modification factor of X-braced steel frame under change in boundary conditions, as change in the direction of strong axis of column and connection support type of column besides variation in storey and bays numbers to be 21 frames and each frame has 8 different boundary conditions as sum of 168 cases for analysis. These frames were analyzed by using nonlinear static “pushover” analysis. As results of this study change in support type and direction of strong axis of column give large change in value of R-factor; the minimum value was 4.37 and maximum value 10.97. Minimum value is close to code value that’s mean the code is more conservative in suggesting of R-factor and gives a large factor of safety. Change in the location of bracing gives change in value of R-factor for all boundary conditions. Change in direction of strong axis of columns and support type didn’t give change in value of fundamental period, all boundary conditions. Keywords: Response modification factor, Ductility reduction factor, Overstrength factor, Boundary conditions, Brace frame, Nonlinear static analysis “Pushover”

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

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

Stabilizing local boundary conditions for two-dimensional shallow water equations

KAUST Repository

Dia, Ben Mansour

2018-03-27

In this article, we present a sub-critical two-dimensional shallow water flow regulation. From the energy estimate of a set of one-dimensional boundary stabilization problems, we obtain a set of polynomial equations with respect to the boundary values as a requirement for the energy decrease. Using the Riemann invariant analysis, we build stabilizing local boundary conditions that guarantee the stability of the hydrodynamical state around a given steady state. Numerical results for the controller applied to the nonlinear problem demonstrate the performance of the method.

An Efficient Numerical Approach for Solving Nonlinear Coupled Hyperbolic Partial Differential Equations with Nonlocal Conditions

Directory of Open Access Journals (Sweden)

A. H. Bhrawy

2014-01-01

Full Text Available One of the most important advantages of collocation method is the possibility of dealing with nonlinear partial differential equations (PDEs as well as PDEs with variable coefficients. A numerical solution based on a Jacobi collocation method is extended to solve nonlinear coupled hyperbolic PDEs with variable coefficients subject to initial-boundary nonlocal conservation conditions. This approach, based on Jacobi polynomials and Gauss-Lobatto quadrature integration, reduces solving the nonlinear coupled hyperbolic PDEs with variable coefficients to a system of nonlinear ordinary differential equation which is far easier to solve. In fact, we deal with initial-boundary coupled hyperbolic PDEs with variable coefficients as well as initial-nonlocal conditions. Using triangular, soliton, and exponential-triangular solutions as exact solutions, the obtained results show that the proposed numerical algorithm is efficient and very accurate.

Artificial boundary conditions for certain evolution PDEs with cubic nonlinearity for non-compactly supported initial data

Science.gov (United States)

Vaibhav, V.

2011-04-01

The paper addresses the problem of constructing non-reflecting boundary conditions for two types of one dimensional evolution equations, namely, the cubic nonlinear Schrödinger (NLS) equation, ∂tu+Lu-iχ|u|2u=0 with L≡-i∂x2, and the equation obtained by letting L≡∂x3. The usual restriction of compact support of the initial data is relaxed by allowing it to have a constant amplitude along with a linear phase variation outside a compact domain. We adapt the pseudo-differential approach developed by Antoine et al. (2006) [5] for the NLS equation to the second type of evolution equation, and further, extend the scheme to the aforementioned class of initial data for both of the equations. In addition, we discuss efficient numerical implementation of our scheme and produce the results of several numerical experiments demonstrating its effectiveness.

A new approach to implement absorbing boundary condition in biomolecular electrostatics.

Science.gov (United States)

Goni, Md Osman

2013-01-01

This paper discusses a novel approach to employ the absorbing boundary condition in conjunction with the finite-element method (FEM) in biomolecular electrostatics. The introduction of Bayliss-Turkel absorbing boundary operators in electromagnetic scattering problem has been incorporated by few researchers. However, in the area of biomolecular electrostatics, this boundary condition has not been investigated yet. The objective of this paper is twofold. First, to solve nonlinear Poisson-Boltzmann equation using Newton's method and second, to find an efficient and acceptable solution with minimum number of unknowns. In this work, a Galerkin finite-element formulation is used along with a Bayliss-Turkel absorbing boundary operator that explicitly accounts for the open field problem by mapping the Sommerfeld radiation condition from the far field to near field. While the Bayliss-Turkel condition works well when the artificial boundary is far from the scatterer, an acceptable tolerance of error can be achieved with the second order operator. Numerical results on test case with simple sphere show that the treatment is able to reach the same level of accuracy achieved by the analytical method while using a lower grid density. Bayliss-Turkel absorbing boundary condition (BTABC) combined with the FEM converges to the exact solution of scattering problems to within discretization error.

MHD axisymmetric flow of power-law fluid over an unsteady stretching sheet with convective boundary conditions

Directory of Open Access Journals (Sweden)

Jawad Ahmed

Full Text Available This paper examines the boundary layer flow and heat transfer characteristic in power law fluid model over unsteady radially stretching sheet under the influence of convective boundary conditions. A uniform magnetic field is applied transversely to the direction of the flow. The governing time dependent nonlinear boundary layer equations are reduced into nonlinear ordinary differential equations with the help of similarity transformations. The transformed coupled ordinary differential equations are then solved analytically by homotopy analysis method (HAM and numerically by shooting procedure. Effects of various governing parameters like, power law index n, magnetic parameter M, unsteadiness A, suction/injection S, Biot number γ and generalized Prandtl number Pr on velocity, temperature, local skin friction and the local Nusselt number are studied and discussed. It is found from the analysis that the magnetic parameter diminishes the velocity profile and the corresponding thermal boundary layer thickness. Keywords: Axisymmetric flow, Power law fluid, Unsteady stretching, Convective boundary conditions

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.

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

On Perturbative Cubic Nonlinear Schrodinger Equations under Complex Nonhomogeneities and Complex Initial Conditions

Directory of Open Access Journals (Sweden)

Magdy A. El-Tawil

2009-01-01

Full Text Available A perturbing nonlinear Schrodinger equation is studied under general complex nonhomogeneities and complex initial conditions for zero boundary conditions. The perturbation method together with the eigenfunction expansion and variational parameters methods are used to introduce an approximate solution for the perturbative nonlinear case for which a power series solution is proved to exist. Using Mathematica, the symbolic solution algorithm is tested through computing the possible approximations under truncation procedures. The method of solution is illustrated through case studies and figures.

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

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

Czech Academy of Sciences Publication Activity Database

Dilna, N.; Rontó, András

2010-01-01

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

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

Nonlinear parabolic problems with Neumann-type boundary conditions and L^1-data

Directory of Open Access Journals (Sweden)

Abderrahmane El Hachimi

2007-11-01

$$ \\frac{\\partial u}{\\partial t}-\\triangle_{p}u+\\alpha(u=f \\quad \\text{in } ]0,\\ T[\\times\\Omega, $$ with Neumann-type boundary conditions and initial data in $L^1$. Our approach is based essentially on the time discretization technique by Euler forward scheme.

Boundary control of nonlinear coupled heat systems using backstepping

KAUST Repository

Bendevis, Paul

2016-10-20

A state feedback boundary controller is designed for a 2D coupled PDE system modelling heat transfer in a membrane distillation system for water desalination. Fluid is separated into two compartments with nonlinear coupling at a membrane boundary. The controller sets the temperature on one boundary in order to track a temperature difference across the membrane boundary. The control objective is achieved by an extension of backstepping methods to these coupled equations. Stability of the target system via Lyapunov like methods, and the invertibility of the integral transformation are used to show the stability of the tracking error.

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

Directory of Open Access Journals (Sweden)

J. Gwinner

2013-01-01

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

Nonlinear interaction of the surface waves at a plasma boundary

International Nuclear Information System (INIS)

Dolgopolov, V.V.; El-Naggar, I.A.; Hussein, A.M.; Khalil, Sh.M.

1976-01-01

Amplitudes of electromagnetic waves with combination frequencies, radiating from the plasma boundary due to nonlinear interaction of the surface waves, have been found. Previous papers on this subject did not take into account that the tangential components of the electric field of waves with combination frequencies were discontinuous at the plasma boundary. (Auth.)

Buckling of Monopod Bucket Foundations – Influence of Boundary Conditions and Soil-structure Interaction

DEFF Research Database (Denmark)

Madsen, Søren; Pinna, Rodney; Randolph, M. F.

2015-01-01

of large-diameter bucket foundations. Since shell structures are generally sensitive to initially imperfect geometries, eigenmode-affine imperfections are introduced in a nonlinear finite-element analysis. The influence of modelling the real lid structure compared to classic boundary conditions...

Consistent boundary conditions for open strings

International Nuclear Information System (INIS)

Lindstroem, Ulf; Rocek, Martin; Nieuwenhuizen, Peter van

2003-01-01

We study boundary conditions for the bosonic, spinning (NSR) and Green-Schwarz open string, as well as for (1+1)-dimensional supergravity. We consider boundary conditions that arise from (1) extremizing the action, (2) BRST, rigid or local supersymmetry, or κ(Siegel)-symmetry of the action, (3) closure of the set of boundary conditions under the symmetry transformations, and (4) the boundary limits of bulk Euler-Lagrange equations that are 'conjugate' to other boundary conditions. We find corrections to Neumann boundary conditions in the presence of a bulk tachyon field. We discuss a boundary superspace formalism. We also find that path integral quantization of the open string requires an infinite tower of boundary conditions that can be interpreted as a smoothness condition on the doubled interval; we interpret this to mean that for a path-integral formulation of open strings with only Neuman boundary conditions, the description in terms of orientifolds is not just natural, but is actually fundamental

Lattice Boltzmann simulations of pressure-driven flows in microchannels using Navier–Maxwell slip boundary conditions

KAUST Repository

Reis, Tim

2012-01-01

We present lattice Boltzmann simulations of rarefied flows driven by pressure drops along two-dimensional microchannels. Rarefied effects lead to non-zero cross-channel velocities, nonlinear variations in the pressure along the channel. Both effects are absent in flows driven by uniform body forces. We obtain second-order accuracy for the two components of velocity the pressure relative to asymptotic solutions of the compressible Navier-Stokes equations with slip boundary conditions. Since the common lattice Boltzmann formulations cannot capture Knudsen boundary layers, we replace the usual discrete analogs of the specular diffuse reflection conditions from continuous kinetic theory with a moment-based implementation of the first-order Navier-Maxwell slip boundary conditions that relate the tangential velocity to the strain rate at the boundary. We use these conditions to solve for the unknown distribution functions that propagate into the domain across the boundary. We achieve second-order accuracy by reformulating these conditions for the second set of distribution functions that arise in the derivation of the lattice Boltzmann method by an integration along characteristics. Our moment formalism is also valuable for analysing the existing boundary conditions. It reveals the origin of numerical slip in the bounce-back other common boundary conditions that impose conditions on the higher moments, not on the local tangential velocity itself. © 2012 American Institute of Physics.

Hydromagnetic nonlinear thermally radiative nanoliquid flow with Newtonian heat and mass conditions

Directory of Open Access Journals (Sweden)

Muhammad Ijaz Khan

Full Text Available This paper communicates the analysis of MHD three-dimensional flow of Jeffrey nanoliquid over a stretchable surface. Flow due to a bidirectional surface is considered. Heat and mass transfer subject to volume fraction of nanoparticles, heat generation and nonlinear solar radiation are examined. Newtonian heat and mass transportation conditions are employed at surface. Concept of boundary layer is utilized to developed the mathematical problem. The boundary value problem is dictated by ten physical parameters: Deborah number, Hartman number, ratio of stretching rates, thermophoretic parameter, Brownian motion parameter, Prandtl number, temperature ratio parameter, conjugate heat and mass parameters and Lewis number. Convergent solutions are obtained using homotopic procedure. Convergence zone for obtained results is explicitly identified. The obtained solutions are interpreted physically. Keywords: Hydromagnetic flow, Viscoelastic nanofluid, Thermophoretic and Brownian moment, Nonlinear thermal radiation, Heat generation

Dressing method and quadratic bundles related to symmetric spaces. Vanishing boundary conditions

Science.gov (United States)

Valchev, T. I.

2016-02-01

We consider quadratic bundles related to Hermitian symmetric spaces of the type SU(m + n)/S(U(m) × U(n)). The simplest representative of the corresponding integrable hierarchy is given by a multi-component Kaup-Newell derivative nonlinear Schrödinger equation which serves as a motivational example for our general considerations. We extensively discuss how one can apply Zakharov-Shabat's dressing procedure to derive reflectionless potentials obeying zero boundary conditions. Those could be used for one to construct fast decaying solutions to any nonlinear equation belonging to the same hierarchy. One can distinguish between generic soliton type solutions and rational solutions.

Outcome of homogeneous and heterogeneous reactions in Darcy-Forchheimer flow with nonlinear thermal radiation and convective condition

Directory of Open Access Journals (Sweden)

T. Hayat

Full Text Available The present analysis aims to report the consequences of nonlinear radiation, convective condition and heterogeneous-homogeneous reactions in Darcy-Forchheimer flow over a non-linear stretching sheet with variable thickness. Non-uniform magnetic field and nonuniform heat generation/absorption are accounted. The governing boundary layer partial differential equations are converted into a system of nonlinear ordinary differential equations. The computations are organized and the effects of physical variables such as thickness parameter, power index, Hartman number, inertia and porous parameters, radiation parameter, Biot number, Prandtl number, ratio parameter, heat generation parameter and homogeneous-heterogeneous reaction parameter are investigated. The variations of skin friction coefficient and Nusselt number for different interesting variables are plotted and discussed. It is noticed that Biot number and heat generation variable lead to enhance the temperature distribution. The solutal boundary layer thickness decreases for larger homogeneous variable while reverse trend is seen for heterogeneous reaction. Keywords: Variable sheet thickness, Darcy-Forchheimer flow, Homogeneous-heterogeneous reactions, Power-law surface velocity, Convective condition, Heat generation/absorption, Nonlinear radiation

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

Directory of Open Access Journals (Sweden)

Suheel Abdullah Malik

2014-01-01

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

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

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

Directory of Open Access Journals (Sweden)

S. S. Motsa

2013-01-01

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

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

Energy Technology Data Exchange (ETDEWEB)

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

2013-11-15

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

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

International Nuclear Information System (INIS)

Kaikina, Elena I.

2013-01-01

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

On the removal of boundary errors caused by Runge-Kutta integration of non-linear partial differential equations

Science.gov (United States)

Abarbanel, Saul; Gottlieb, David; Carpenter, Mark H.

1994-01-01

It has been previously shown that the temporal integration of hyperbolic partial differential equations (PDE's) may, because of boundary conditions, lead to deterioration of accuracy of the solution. A procedure for removal of this error in the linear case has been established previously. In the present paper we consider hyperbolic (PDE's) (linear and non-linear) whose boundary treatment is done via the SAT-procedure. A methodology is present for recovery of the full order of accuracy, and has been applied to the case of a 4th order explicit finite difference scheme.

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

Spline Collocation Method for Nonlinear Multi-Term Fractional Differential Equation

OpenAIRE

Choe, Hui-Chol; Kang, Yong-Suk

2013-01-01

We study an approximation method to solve nonlinear multi-term fractional differential equations with initial conditions or boundary conditions. First, we transform the nonlinear multi-term fractional differential equations with initial conditions and boundary conditions to nonlinear fractional integral equations and consider the relations between them. We present a Spline Collocation Method and prove the existence, uniqueness and convergence of approximate solution as well as error estimatio...

Reconstruction of boundary conditions from internal conditions using viability theory

KAUST Repository

Hofleitner, Aude; Claudel, Christian G.; Bayen, Alexandre M.

2012-01-01

This article presents a method for reconstructing downstream boundary conditions to a HamiltonJacobi partial differential equation for which initial and upstream boundary conditions are prescribed as piecewise affine functions and an internal condition is prescribed as an affine function. Based on viability theory, we reconstruct the downstream boundary condition such that the solution of the Hamilton-Jacobi equation with the prescribed initial and upstream conditions and reconstructed downstream boundary condition satisfies the internal value condition. This work has important applications for estimation in flow networks with unknown capacity reductions. It is applied to urban traffic, to reconstruct signal timings and temporary capacity reductions at intersections, using Lagrangian sensing such as GPS devices onboard vehicles.

Reconstruction of boundary conditions from internal conditions using viability theory

KAUST Repository

Hofleitner, Aude

2012-06-01

This article presents a method for reconstructing downstream boundary conditions to a HamiltonJacobi partial differential equation for which initial and upstream boundary conditions are prescribed as piecewise affine functions and an internal condition is prescribed as an affine function. Based on viability theory, we reconstruct the downstream boundary condition such that the solution of the Hamilton-Jacobi equation with the prescribed initial and upstream conditions and reconstructed downstream boundary condition satisfies the internal value condition. This work has important applications for estimation in flow networks with unknown capacity reductions. It is applied to urban traffic, to reconstruct signal timings and temporary capacity reductions at intersections, using Lagrangian sensing such as GPS devices onboard vehicles.

Numerical study of magnetohydrodynamics (MHD boundary layer slip flow of a Maxwell nanofluid over an exponentially stretching surface with convective boundary condition

Directory of Open Access Journals (Sweden)

P.BalaAnki Reddy

2017-12-01

Full Text Available This paper focuses on a theoretical analysis of a steady two-dimensional magnetohydrodynamic boundary layer flow of a Maxwell fluid over an exponentially stretching surface in the presence of velocity slip and convective boundary condition. This model is used for a nanofluid, which incorporates the effects of Brownian motion and thermophoresis. The resulting non-linear partial differential equations of the governing flow field are converted into a system of coupled non-linear ordinary differential equations by using suitable similarity transformations, and the resultant equations are then solved numerically by using Runge-Kutta fourth order method along with shooting technique. A parametric study is conducted to illustrate the behavior of the velocity, temperature and concentration. The influence of significant parameters on velocity, temperature, concentration, skin friction coefficient and Nusselt number has been studied and numerical results are presented graphically and in tabular form. The reported numerical results are compared with previously published works on various special cases and are found to be an in excellent agreement. It is found that momentum boundary layer thickness decreases with the increase of magnetic parameter. It can also be found that the thermal boundary layer thickness increases with Brownian motion and thermophoresis parameters.

Boundary conditions in random sequential adsorption

Science.gov (United States)

Cieśla, Michał; Ziff, Robert M.

2018-04-01

The influence of different boundary conditions on the density of random packings of disks is studied. Packings are generated using the random sequential adsorption algorithm with three different types of boundary conditions: periodic, open, and wall. It is found that the finite size effects are smallest for periodic boundary conditions, as expected. On the other hand, in the case of open and wall boundaries it is possible to introduce an effective packing size and a constant correction term to significantly improve the packing densities.

Quantum “violation” of Dirichlet boundary condition

Directory of Open Access Journals (Sweden)

I.Y. Park

2017-02-01

Full Text Available Dirichlet boundary conditions have been widely used in general relativity. They seem at odds with the holographic property of gravity simply because a boundary configuration can be varying and dynamic instead of dying out as required by the conditions. In this work we report what should be a tension between the Dirichlet boundary conditions and quantum gravitational effects, and show that a quantum-corrected black hole solution of the 1PI action no longer obeys, in the naive manner one may expect, the Dirichlet boundary conditions imposed at the classical level. We attribute the ‘violation’ of the Dirichlet boundary conditions to a certain mechanism of the information storage on the boundary.

Quantum “violation” of Dirichlet boundary condition

Energy Technology Data Exchange (ETDEWEB)

Park, I.Y., E-mail: inyongpark05@gmail.com

2017-02-10

Dirichlet boundary conditions have been widely used in general relativity. They seem at odds with the holographic property of gravity simply because a boundary configuration can be varying and dynamic instead of dying out as required by the conditions. In this work we report what should be a tension between the Dirichlet boundary conditions and quantum gravitational effects, and show that a quantum-corrected black hole solution of the 1PI action no longer obeys, in the naive manner one may expect, the Dirichlet boundary conditions imposed at the classical level. We attribute the ‘violation’ of the Dirichlet boundary conditions to a certain mechanism of the information storage on the boundary.

Boundary Observability and Stabilization for Westervelt Type Wave Equations without Interior Damping

International Nuclear Information System (INIS)

Kaltenbacher, Barbara

2010-01-01

In this paper we show boundary observability and boundary stabilizability by linear feedbacks for a class of nonlinear wave equations including the undamped Westervelt model used in nonlinear acoustics. We prove local existence for undamped generalized Westervelt equations with homogeneous Dirichlet boundary conditions as well as global existence and exponential decay with absorbing type boundary conditions.

Investigating effects of boundary conditions on the evaluation of R-factor of un-braced steel frames

Directory of Open Access Journals (Sweden)

Masood M.M. Irheem

2017-08-01

Full Text Available Design of Structures to resist seismic load depends on the theory of dissipation in elastic of energy that already exists in response modification factor “R-factor”. The main problem in codes gives a constant value for R-factor, since change in boundary conditions of building change in behavior of steel frame structures and that effect on R-factor. This study is an attempt to assess overstrength, ductility and response modification factor of un-braced steel frames under change in boundary conditions as change in the direction of strong axis of column and support type beside to variation in story and bay number to be 9 frame and each frame has 8 different boundary conditions as sum of 72 case for analysis. These frames were analyzed by using nonlinear static “pushover” analysis using SAP2000 program. As a result of this study R-factor does not has a constant value, when change in boundary conditions R-factor directly changes, minimum value of 8 boundary conditions is close to the code value that is mean the code is more conservative and give a large factor of safety. Ductility reduction factor increases with increasing number of story for all boundary conditions, but overstrength has different rule. Response modification factor, overstrength factor and ductility reduction factor decrease when fundamentals period increasing for the studied frames.

Investigation of nonlinear 2D bottom transportation dynamics in coastal zone on optimal curvilinear boundary adaptive grids

Directory of Open Access Journals (Sweden)

Sukhinov Alexander

2017-01-01

Full Text Available One of the practically important tasks of hydrophysics for sea coastal systems is the problem of modeling and forecasting bottom sediment transportation. A number of problems connected to ship safety traffic, water medium condition near the coastal line etc. depends on forecasting bottom deposit transportation under natural and technogenic influences. Coastal systems are characterized by a complicated form of coastline - the presence of long, narrow and curvilinear peninsulas and bays. Water currents and waves near the beach are strongly depend on complicated coastal line and in turn, exert on the bottom sediment transportation near the shore. The use of rectangular grids in the construction of discrete models leads to significant errors in both the specification of boundary conditions and in the modeling of hydrophysical processes in the coastal zone. In this paper, we consider the construction of a finite-element approximation of the initial-boundary value problem for the spatially two-dimensional linearized equation of sediment transportation using optimal boundary-adaptive grid. First, the linearization of a spatially two-dimensional nonlinear parabolic equation on the time grid is performed-when the coefficients of the equation that are nonlinearly dependent on the bottom relief function are set on the previous time layer, and the corresponding initial conditions are used on the first time layer. The algorithm for constructing the grid is based on the procedure for minimizing the generalized Dirichlet functional. On the constructed grid, finite element approximation using bilinear basis functions is performed, which completes the construction of a discrete model for the given problem. The using of curvilinear boundary adaptive grids leads to decreasing of total grid number in 5-20 times and respectively the total modeling time and/or it allows to improve modeling accuracy.

Nonlinear Lyapunov-based boundary control of distributed heat transfer mechanisms in membrane distillation plant

KAUST Repository

Eleiwi, Fadi; Laleg-Kirati, Taous-Meriem

2015-01-01

This paper presents a nonlinear Lyapunov-based boundary control for the temperature difference of a membrane distillation boundary layers. The heat transfer mechanisms inside the process are modeled with a 2D advection-diffusion equation. The model

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

Directory of Open Access Journals (Sweden)

Shihuang Hong

2009-01-01

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

Microlocal approach towards construction of nonreflecting boundary conditions

Science.gov (United States)

Vaibhav, V.

2014-09-01

This paper addresses the problem of construction of non-reflecting boundary condition for certain second-order nonlinear dispersive equations. It is shown that using the concept of microlocality it is possible to relax the requirement of compact support of the initial data. The method is demonstrated for a class of initial data such that outside the computational domain it behaves like a continuous-wave. The generalization is detailed for two existing schemes in the framework of pseudo-differential calculus, namely, Szeftel's method (Szeftel (2006) [1]) and gauge transformation strategy (Antoine et al. (2006) [2]). Efficient numerical implementation is discussed and a comparative performance analysis is presented. The paper also briefly surveys the possibility of extension of the method to higher-dimensional PDEs.

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

The influence of compressibility on nonlinear spectral energy transfer - Part 2: Effect on hypersonic boundary layer transition

Science.gov (United States)

Mittal, Ankita; Girimaji, Sharath

2017-11-01

We examine the effect of compressible spectral energy transfer in the nonlinear regime of transition to turbulence of hypersonic boundary layers. The nature of spectral energy transfer between perturbation modes is profoundly influenced by two compressibility mechanisms. First and foremost, the emergence of nonlinear pressure-dilatation mechanism leads to kinetic-internal energy exchange within the perturbation field. Such interchange is absent in incompressible flow as pressure merely reorients the perturbation amplitude vector while conserving kinetic energy. Secondly, the nature of triadic interactions also changes due to variability in density. In this work, we demonstrate that the efficiency of nonlinear spectral energy transfer is diminished in compressible boundary layers. Emergence of new perturbation modes or `broad-banding' of the perturbation field is significantly delayed in comparison to incompressible boundary layer undergoing transition. A significant amount of perturbation energy is transformed to internal energy and thus unavailable for `tripping' the flow into turbulent state. These factors profoundly change the nature of the nonlinear stage of transition in compressible boundary layer leading to delayed onset of full-fledged turbulence.

The topology of non-linear global carbon dynamics: from tipping points to planetary boundaries

International Nuclear Information System (INIS)

Anderies, J M; Carpenter, S R; Steffen, Will; Rockström, Johan

2013-01-01

We present a minimal model of land use and carbon cycle dynamics and use it to explore the relationship between non-linear dynamics and planetary boundaries. Only the most basic interactions between land cover and terrestrial, atmospheric, and marine carbon stocks are considered in the model. Our goal is not to predict global carbon dynamics as it occurs in the actual Earth System. Rather, we construct a conceptually reasonable heuristic model of a feedback system between different carbon stocks that captures the qualitative features of the actual Earth System and use it to explore the topology of the boundaries of what can be called a ‘safe operating space’ for humans. The model analysis illustrates the existence of dynamic, non-linear tipping points in carbon cycle dynamics and the potential complexity of planetary boundaries. Finally, we use the model to illustrate some challenges associated with navigating planetary boundaries. (letter)

The topology of non-linear global carbon dynamics: from tipping points to planetary boundaries

Science.gov (United States)

Anderies, J. M.; Carpenter, S. R.; Steffen, Will; Rockström, Johan

2013-12-01

We present a minimal model of land use and carbon cycle dynamics and use it to explore the relationship between non-linear dynamics and planetary boundaries. Only the most basic interactions between land cover and terrestrial, atmospheric, and marine carbon stocks are considered in the model. Our goal is not to predict global carbon dynamics as it occurs in the actual Earth System. Rather, we construct a conceptually reasonable heuristic model of a feedback system between different carbon stocks that captures the qualitative features of the actual Earth System and use it to explore the topology of the boundaries of what can be called a ‘safe operating space’ for humans. The model analysis illustrates the existence of dynamic, non-linear tipping points in carbon cycle dynamics and the potential complexity of planetary boundaries. Finally, we use the model to illustrate some challenges associated with navigating planetary boundaries.

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.

Outcome of homogeneous and heterogeneous reactions in Darcy-Forchheimer flow with nonlinear thermal radiation and convective condition

Science.gov (United States)

Hayat, T.; Shah, Faisal; Alsaedi, A.; Hussain, Zakir

The present analysis aims to report the consequences of nonlinear radiation, convective condition and heterogeneous-homogeneous reactions in Darcy-Forchheimer flow over a non-linear stretching sheet with variable thickness. Non-uniform magnetic field and nonuniform heat generation/absorption are accounted. The governing boundary layer partial differential equations are converted into a system of nonlinear ordinary differential equations. The computations are organized and the effects of physical variables such as thickness parameter, power index, Hartman number, inertia and porous parameters, radiation parameter, Biot number, Prandtl number, ratio parameter, heat generation parameter and homogeneous-heterogeneous reaction parameter are investigated. The variations of skin friction coefficient and Nusselt number for different interesting variables are plotted and discussed. It is noticed that Biot number and heat generation variable lead to enhance the temperature distribution. The solutal boundary layer thickness decreases for larger homogeneous variable while reverse trend is seen for heterogeneous reaction.

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.

Molecular dynamics simulation of UO2 nanocrystals melting under isolated and periodic boundary conditions

International Nuclear Information System (INIS)

Boyarchenkov, A.S.; Potashnikov, S.I.; Nekrasov, K.A.; Kupryazhkin, A.Ya.

2012-01-01

Highlights: ► We perform MD simulation of UO 2 nanocrystals melting (in range of 768–49 152 ions). ► T(P) melting curves intersect zero near −20 GPa and saturate near 25 GPa. ► Reciprocal size dependences of nanocrystal melting point decrease nonlinearly. ► Linear and parabolic extrapolations to macroscopic values are considered. ► Melting point and density jump are reproduced, but heat of fusion is underestimated. - Abstract: Melting of uranium dioxide (UO 2 ) nanocrystals has been studied by molecular dynamics (MD) simulation. Ten recent and widely used sets of pair potentials were assessed in the rigid ion approximation. Both isolated (in vacuum) and periodic boundary conditions (PBC) were explored. Using barostat under PBC the pressure dependences of melting point were obtained. These curves intersected zero near −20 GPa, saturated near 25 GPa and increased nonlinearly in between. Using simulation of surface under isolated boundary conditions (IBC) recommended melting temperature and density jump were successfully reproduced. However, the heat of fusion is still underestimated. These melting characteristics were calculated for nanocrystals of cubic shape in the range of 768–49 152 particles (volume range of 10–1000 nm 3 ). The obtained reciprocal size dependences decreased nonlinearly. Linear and parabolic extrapolations to macroscopic values are considered. The parabolic one is found to be better suited for analysis of the data on temperature and heat of melting.

Spectral boundary conditions and solitonic solutions in a classical Sellmeier dielectric

Energy Technology Data Exchange (ETDEWEB)

Belgiorno, F. [Politecnico di Milano, Dipartimento di Matematica, Milan (Italy); INdAM-GNFM, Rome (Italy); INFN, Milan (Italy); Cacciatori, S.L. [Universita dell' Insubria, Department of Science and High Technology, Como (Italy); INFN, Milan (Italy); Vigano, A. [Universita degli Studi di Milano, Dipartimento di Fisica, Milan (Italy)

2017-06-15

Electromagnetic field interactions in a dielectric medium represent a longstanding field of investigation, both at the classical level and at the quantum one. We propose a 1+1 dimensional toy-model which consists of an half-line filling dielectric medium, with the aim to set up a simplified situation where technicalities related to gauge invariance and, as a consequence, physics of constrained systems are avoided, and still interesting features appear. In particular, we simulate the electromagnetic field and the polarization field by means of two coupled scalar fields φ, ψ, respectively, in a Hopfield-like model. We find that, in order to obtain a physically meaningful behavior for the model, one has to introduce spectral boundary conditions depending on the particle spectrum one is dealing with. This is the first interesting achievement of our analysis. The second relevant achievement is that, by introducing a nonlinear contribution in the polarization field ψ, with the aim of mimicking a third order nonlinearity in a nonlinear dielectric, we obtain solitonic solutions in the Hopfield model framework, whose classical behavior is analyzed too. (orig.)

On the Stability of Three-Dimensional Boundary Layers. Part 1; Linear and Nonlinear Stability

Science.gov (United States)

Janke, Erik; Balakumar, Ponnampalam

1999-01-01

The primary stability of incompressible three-dimensional boundary layers is investigated using the Parabolized Stability Equations (PSE). We compute the evolution of stationary and traveling disturbances in the linear and nonlinear region prior to transition. As model problems, we choose Swept Hiemenz Flow and the DLR Transition Experiment. The primary stability results for Swept Hiemenz Flow agree very well with computations by Malik et al. For the DLR Experiment, the mean flow profiles are obtained by solving the boundary layer equations for the measured pressure distribution. Both linear and nonlinear results show very good agreement with the experiment.

Nonlinear evolution of tearing and coalescence instability with free boundary conditions

International Nuclear Information System (INIS)

Malara, F.; Veltri, P.; Carbone, V.

1990-01-01

The nonlinear evolution of a reconnection instability in a plane current sheet is described. In particular, the appearance of coalescence instability was studied, which follows the formation of a chain of magnetic islands due to the tearing instability. In order to describe realistically this phonemenon, the time evolution of all the unstable modes which are present in the spectrum at the same time is considered. Moreover, this study allows to investigate the turbulent energy cascade which forms owing to the nonlinear coupling between such modes. (R.P.) 9 refs.; 6 figs

Stabilization of Hypersonic Boundary Layers by Linear and Nonlinear Optimal Perturbations

Science.gov (United States)

Paredes, Pedro; Choudhari, Meelan M.; Li, Fei

2017-01-01

The effect of stationary, finite-amplitude, linear and nonlinear optimal perturbations on the modal disturbance growth in a Mach 6 axisymmetric flow over a 7 deg. half-angle cone with 0:126 mm nose radius and 0:305 m length is investigated. The freestream parameters (M = 6, Re(exp 1) = 18 x 10(exp. 6) /m) are selected to match the flow conditions of a previous experiment in the VKI H3 hypersonic tunnel. Plane-marching parabolized stability equations are used in conjunction with a partial-differential equation based planar eigenvalue analysis to characterize the boundary layer instability in the presence of azimuthally periodic streaks. The streaks are observed to stabilize nominally planar Mack mode instabilities, although oblique Mack mode and first-mode disturbances are destabilized. Experimentally measured transition onset in the absence of any streaks correlates with an amplification factor of N = 6 for the planar Mack modes. For high enough streak amplitudes, the transition threshold of N = 6 is not reached by the Mack mode instabilities within the length of the cone; however, subharmonic first-mode instabilities, which are destabilized by the presence of the streaks, do reach N = 6 near the end of the cone. The highest stabilization is observed at streak amplitudes of approximately 20 percent of the freestream velocity. Because the use of initial disturbance profiles based on linear optimal growth theory may yield suboptimal control in the context of nonlinear streaks, the computational predictions are extended to nonlinear optimal growth theory. Results show that by using nonlinearly optimal perturbation leads to slightly enhanced stabilization of plane Mack mode disturbances as well as reduced destabilization of subharmonic first-mode disturbances.

Waveguides with Absorbing Boundaries: Nonlinearity Controlled by an Exceptional Point and Solitons

Science.gov (United States)

Midya, Bikashkali; Konotop, Vladimir V.

2017-07-01

We reveal the existence of continuous families of guided single-mode solitons in planar waveguides with weakly nonlinear active core and absorbing boundaries. Stable propagation of TE and TM-polarized solitons is accompanied by attenuation of all other modes, i.e., the waveguide features properties of conservative and dissipative systems. If the linear spectrum of the waveguide possesses exceptional points, which occurs in the case of TM polarization, an originally focusing (defocusing) material nonlinearity may become effectively defocusing (focusing). This occurs due to the geometric phase of the carried eigenmode when the surface impedance encircles the exceptional point. In its turn, the change of the effective nonlinearity ensures the existence of dark (bright) solitons in spite of focusing (defocusing) Kerr nonlinearity of the core. The existence of an exceptional point can also result in anomalous enhancement of the effective nonlinearity. In terms of practical applications, the nonlinearity of the reported waveguide can be manipulated by controlling the properties of the absorbing cladding.

Nonlocal beam models for buckling of nanobeams using state-space method regarding different boundary conditions

International Nuclear Information System (INIS)

Sahmani, S.; Ansari, R.

2011-01-01

Buckling analysis of nanobeams is investigated using nonlocal continuum beam models of the different classical beam theories namely as Euler-Bernoulli beam theory (EBT), Timoshenko beam theory (TBT), and Levinson beam theory (LBT). To this end, Eringen's equations of nonlocal elasticity are incorporated into the classical beam theories for buckling of nanobeams with rectangular cross-section. In contrast to the classical theories, the nonlocal elastic beam models developed here have the capability to predict critical buckling loads that allowing for the inclusion of size effects. The values of critical buckling loads corresponding to four commonly used boundary conditions are obtained using state-space method. The results are presented for different geometric parameters, boundary conditions, and values of nonlocal parameter to show the effects of each of them in detail. Then the results are fitted with those of molecular dynamics simulations through a nonlinear least square fitting procedure to find the appropriate values of nonlocal parameter for the buckling analysis of nanobeams relevant to each type of nonlocal beam model and boundary conditions analysis

Nonlocal beam models for buckling of nanobeams using state-space method regarding different boundary conditions

Energy Technology Data Exchange (ETDEWEB)

Sahmani, S.; Ansari, R. [University of Guilan, Rasht (Iran, Islamic Republic of)

2011-09-15

Buckling analysis of nanobeams is investigated using nonlocal continuum beam models of the different classical beam theories namely as Euler-Bernoulli beam theory (EBT), Timoshenko beam theory (TBT), and Levinson beam theory (LBT). To this end, Eringen's equations of nonlocal elasticity are incorporated into the classical beam theories for buckling of nanobeams with rectangular cross-section. In contrast to the classical theories, the nonlocal elastic beam models developed here have the capability to predict critical buckling loads that allowing for the inclusion of size effects. The values of critical buckling loads corresponding to four commonly used boundary conditions are obtained using state-space method. The results are presented for different geometric parameters, boundary conditions, and values of nonlocal parameter to show the effects of each of them in detail. Then the results are fitted with those of molecular dynamics simulations through a nonlinear least square fitting procedure to find the appropriate values of nonlocal parameter for the buckling analysis of nanobeams relevant to each type of nonlocal beam model and boundary conditions analysis.

On Hydroelastic Body-Boundary Condition of Floating Structures

DEFF Research Database (Denmark)

Xia, Jinzhu

1996-01-01

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

Nonlinear evolution of Mack modes in a hypersonic boundary layer

Science.gov (United States)

Chokani, Ndaona

2005-01-01

In hypersonic boundary layer flows the nonlinear disturbance evolution occurs relatively slowly over a very long length scale and has a profound effect on boundary layer transition. In the case of low-level freestream disturbances and negligible surface roughness, the transition is due to the modal growth of exponentially growing Mack modes that are destabilized by wall cooling. Cross-bicoherence measurements, derived from hot-wire data acquired in a quiet hypersonic tunnel, are used to identify and quantify phase-locked, quadratic sum and difference interactions involving the Mack modes. In the early stages of the nonlinear disturbance evolution, cross-bicoherence measurements indicate that the energy exchange between the Mack mode and the mean flow first occurs to broaden the sidebands; this is immediately followed by a sum interaction of the Mack mode to generate the first harmonic. In the next stages of the nonlinear disturbance evolution, there is a difference interaction of the first harmonic, which is also thought to contribute to the mean flow distortion. This difference interaction, in the latter stages, is also accompanied by a difference interaction between Mack mode and first harmonic, and a sum interaction, which forces the second harmonic. Analysis using the digital complex demodulation technique, shows that the low-frequency, phase-locked interaction that is identified in the cross bicoherence when the Mack mode and first harmonic have large amplitudes, arises due to the amplitude modulation of Mack mode and first harmonic.

Application of the comparison principle to analysis of nonlinear systems. [using Lipschitz condition and differential equations

Science.gov (United States)

Gunderson, R. W.

1975-01-01

A comparison principle based on a Kamke theorem and Lipschitz conditions is presented along with its possible applications and modifications. It is shown that the comparison lemma can be used in the study of such areas as classical stability theory, higher order trajectory derivatives, Liapunov functions, boundary value problems, approximate dynamic systems, linear and nonlinear systems, and bifurcation analysis.

A non-local computational boundary condition for duct acoustics

Science.gov (United States)

Zorumski, William E.; Watson, Willie R.; Hodge, Steve L.

1994-01-01

A non-local boundary condition is formulated for acoustic waves in ducts without flow. The ducts are two dimensional with constant area, but with variable impedance wall lining. Extension of the formulation to three dimensional and variable area ducts is straightforward in principle, but requires significantly more computation. The boundary condition simulates a nonreflecting wave field in an infinite duct. It is implemented by a constant matrix operator which is applied at the boundary of the computational domain. An efficient computational solution scheme is developed which allows calculations for high frequencies and long duct lengths. This computational solution utilizes the boundary condition to limit the computational space while preserving the radiation boundary condition. The boundary condition is tested for several sources. It is demonstrated that the boundary condition can be applied close to the sound sources, rendering the computational domain small. Computational solutions with the new non-local boundary condition are shown to be consistent with the known solutions for nonreflecting wavefields in an infinite uniform duct.

Inverse scattering transform for the vector nonlinear Schroedinger equation with nonvanishing boundary conditions

International Nuclear Information System (INIS)

Prinari, Barbara; Ablowitz, Mark J.; Biondini, Gino

2006-01-01

The inverse scattering transform for the vector defocusing nonlinear Schroedinger (NLS) equation with nonvanishing boundary values at infinity is constructed. The direct scattering problem is formulated on a two-sheeted covering of the complex plane. Two out of the six Jost eigenfunctions, however, do not admit an analytic extension on either sheet of the Riemann surface. Therefore, a suitable modification of both the direct and the inverse problem formulations is necessary. On the direct side, this is accomplished by constructing two additional analytic eigenfunctions which are expressed in terms of the adjoint eigenfunctions. The discrete spectrum, bound states and symmetries of the direct problem are then discussed. In the most general situation, a discrete eigenvalue corresponds to a quartet of zeros (poles) of certain scattering data. The inverse scattering problem is formulated in terms of a generalized Riemann-Hilbert (RH) problem in the upper/lower half planes of a suitable uniformization variable. Special soliton solutions are constructed from the poles in the RH problem, and include dark-dark soliton solutions, which have dark solitonic behavior in both components, as well as dark-bright soliton solutions, which have one dark and one bright component. The linear limit is obtained from the RH problem and is shown to correspond to the Fourier transform solution obtained from the linearized vector NLS system

On nonlinear statistical thermodynamics of boundary plasma with postactions

International Nuclear Information System (INIS)

Temko, S.W.; Temko, K.W.; Kuz'min, S.K.

1992-01-01

The authors use the statistical thermodynamics of small systems proposed before their publications for boundary weakly ionized plasma with postaction. Boundary properties of the plasma is taken into account by two ways: (1) suppose that only small number of very quick particles are able to leave the cloud having done entrance into outer medium work; (2) take into account the interaction between particles and inner surface of the cloud. Interactions in the boundary plasma are described by corresponding potential functions. The potential functions are mathematical models of real interactions in boundary plasma. Choosing of potential functions, their numerical parameters, geometrical form and dimensions of the cloud is made by using the methods of optimal experiment planning, maximum likelihood and computer experiment. Free energy of the cloud is a likelihood function. State of boundary plasma with admixtures is described by vector-density of particles distribution. Term ''distribution'' is used here in Sobolev-Schwartc sense. The authors obtain the vector-density of particles distribution in cloud which gives the condition minimum of free energy for every time moment under quasistatistical equilibrium. The system of conditions for free energy conditional minimizing for every time moment includes integral equilibrium equations, ''non-hard normalization'' and additional conditions taken as a result of analyzing physical and physical-chemical nature of boundary plasma. To obtain conditional minimum of free energy it is necessary to solve the system of conditions. First of all they solve equilibrium problem by the authors method. They obtain vector-density of particles distribution in the cloud. Then using method of random walk with postaction between sets of random walk process they build distribution function of random vector-density

Nonlinear models for autoregressive conditional heteroskedasticity

DEFF Research Database (Denmark)

Teräsvirta, Timo

This paper contains a brief survey of nonlinear models of autore- gressive conditional heteroskedasticity. The models in question are parametric nonlinear extensions of the original model by Engle (1982). After presenting the individual models, linearity testing and parameter estimation are discu...

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.

Integral Method of Boundary Characteristics: Neumann Condition

Science.gov (United States)

Kot, V. A.

2018-05-01

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

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

Inverse scattering transform and soliton solutions for square matrix nonlinear Schrödinger equations with non-zero boundary conditions

Science.gov (United States)

Prinari, Barbara; Demontis, Francesco; Li, Sitai; Horikis, Theodoros P.

2018-04-01

The inverse scattering transform (IST) with non-zero boundary conditions at infinity is developed for an m × m matrix nonlinear Schrödinger-type equation which, in the case m = 2, has been proposed as a model to describe hyperfine spin F = 1 spinor Bose-Einstein condensates with either repulsive interatomic interactions and anti-ferromagnetic spin-exchange interactions (self-defocusing case), or attractive interatomic interactions and ferromagnetic spin-exchange interactions (self-focusing case). The IST for this system was first presented by Ieda et al. (2007) , using a different approach. In our formulation, both the direct and the inverse problems are posed in terms of a suitable uniformization variable which allows to develop the IST on the standard complex plane, instead of a two-sheeted Riemann surface or the cut plane with discontinuities along the cuts. Analyticity of the scattering eigenfunctions and scattering data, symmetries, properties of the discrete spectrum, and asymptotics are derived. The inverse problem is posed as a Riemann-Hilbert problem for the eigenfunctions, and the reconstruction formula of the potential in terms of eigenfunctions and scattering data is provided. In addition, the general behavior of the soliton solutions is analyzed in detail in the 2 × 2 self-focusing case, including some special solutions not previously discussed in the literature.

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

Directory of Open Access Journals (Sweden)

Aomar Anane

2011-01-01

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

Fluid-solid boundary conditions for multiparticle collision dynamics

International Nuclear Information System (INIS)

Whitmer, Jonathan K; Luijten, Erik

2010-01-01

The simulation of colloidal particles suspended in solvent requires an accurate representation of the interactions between the colloids and the solvent molecules. Using the multiparticle collision dynamics method, we examine several proposals for stick boundary conditions, studying their properties in both plane Poiseuille flow (where fluid interacts with the boundary of a stationary macroscopic solid) and particle-based colloid simulations (where the boundaries are thermally affected and in motion). In addition, our simulations compare various collision rules designed to remove spurious slip near solid surfaces, and the effects of these rules on the thermal motion of colloidal particles. Furthermore, we demonstrate that stochastic reflection of the fluid at solid boundaries fails to faithfully represent stick boundary conditions, and conclude that bounce-back conditions should be applied at both mobile and stationary surfaces. Finally, we generalize these ideas to create partial slip boundary conditions at both stationary and mobile surfaces.

Surface free energy for systems with integrable boundary conditions

International Nuclear Information System (INIS)

Goehmann, Frank; Bortz, Michael; Frahm, Holger

2005-01-01

The surface free energy is the difference between the free energies for a system with open boundary conditions and the same system with periodic boundary conditions. We use the quantum transfer matrix formalism to express the surface free energy in the thermodynamic limit of systems with integrable boundary conditions as a matrix element of certain projection operators. Specializing to the XXZ spin-1/2 chain we introduce a novel 'finite temperature boundary operator' which characterizes the thermodynamical properties of surfaces related to integrable boundary conditions

Nonlinear Lyapunov-based boundary control of distributed heat transfer mechanisms in membrane distillation plant

KAUST Repository

Eleiwi, Fadi

2015-07-01

This paper presents a nonlinear Lyapunov-based boundary control for the temperature difference of a membrane distillation boundary layers. The heat transfer mechanisms inside the process are modeled with a 2D advection-diffusion equation. The model is semi-descretized in space, and a nonlinear state-space representation is provided. The control is designed to force the temperature difference along the membrane sides to track a desired reference asymptotically, and hence a desired flux would be generated. Certain constraints are put on the control law inputs to be within an economic range of energy supplies. The effect of the controller gain is discussed. Simulations with real process parameters for the model, and the controller are provided. © 2015 American Automatic Control Council.

A two-dimensional vibration analysis of piezoelectrically actuated microbeam with nonideal boundary conditions

Science.gov (United States)

Rezaei, M. P.; Zamanian, M.

2017-01-01

In this paper, the influences of nonideal boundary conditions (due to flexibility) on the primary resonant behavior of a piezoelectrically actuated microbeam have been studied, for the first time. The structure has been assumed to treat as an Euler-Bernoulli beam, considering the effects of geometric nonlinearity. In this work, the general nonideal supports have been modeled as a the combination of horizontal, vertical and rotational springs, simultaneously. Allocating particular values to the stiffness of these springs provides the mathematical models for the majority of boundary conditions. This consideration leads to use a two-dimensional analysis of the multiple scales method instead of previous works' method (one-dimensional analysis). If one neglects the nonideal effects, then this paper would be an effort to solve the two-dimensional equations of motion without a need of a combination of these equations using the shortening or stretching effect. Letting the nonideal effects equal to zero and comparing their results with the results of previous approaches have been demonstrated the accuracy of the two-dimensional solutions. The results have been identified the unique effects of constraining and stiffening of boundaries in horizontal, vertical and rotational directions. This means that it is inaccurate to suppose the nonideality of supports only in one or two of these directions like as previous works. The findings are of vital importance as a better prediction of the frequency response for the nonideal supports. Furthermore, the main findings of this effort can help to choose appropriate boundary conditions for desired systems.

Existence theory for sequential fractional differential equations with anti-periodic type boundary conditions

Directory of Open Access Journals (Sweden)

Aqlan Mohammed H.

2016-01-01

Full Text Available We develop the existence theory for sequential fractional differential equations involving Liouville-Caputo fractional derivative equipped with anti-periodic type (non-separated and nonlocal integral boundary conditions. Several existence criteria depending on the nonlinearity involved in the problems are presented by means of a variety of tools of the fixed point theory. The applicability of the results is shown with the aid of examples. Our results are not only new in the given configuration but also yield some new special cases for specific choices of parameters involved in the problems.

Discrete transparent boundary conditions for Schroedinger-type equations

International Nuclear Information System (INIS)

Schmidt, F.; Yevick, D.

1997-01-01

We present a general technique for constructing nonlocal transparent boundary conditions for one-dimensional Schroedinger-type equations. Our method supplies boundary conditions for the θ-family of implicit one-step discretizations of Schroedinger's equation in time. The use of Mikusinski's operator approach in time avoids direct and inverse transforms between time and frequency domains and thus implements the boundary conditions in a direct manner. 14 refs., 9 figs

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.

A Galleria Boundary Element Method for two-dimensional nonlinear magnetostatics

Science.gov (United States)

Brovont, Aaron D.

The Boundary Element Method (BEM) is a numerical technique for solving partial differential equations that is used broadly among the engineering disciplines. The main advantage of this method is that one needs only to mesh the boundary of a solution domain. A key drawback is the myriad of integrals that must be evaluated to populate the full system matrix. To this day these integrals have been evaluated using numerical quadrature. In this research, a Galerkin formulation of the BEM is derived and implemented to solve two-dimensional magnetostatic problems with a focus on accurate, rapid computation. To this end, exact, closed-form solutions have been derived for all the integrals comprising the system matrix as well as those required to compute fields in post-processing; the need for numerical integration has been eliminated. It is shown that calculation of the system matrix elements using analytical solutions is 15-20 times faster than with numerical integration of similar accuracy. Furthermore, through the example analysis of a c-core inductor, it is demonstrated that the present BEM formulation is a competitive alternative to the Finite Element Method (FEM) for linear magnetostatic analysis. Finally, the BEM formulation is extended to analyze nonlinear magnetostatic problems via the Dual Reciprocity Method (DRBEM). It is shown that a coarse, meshless analysis using the DRBEM is able to achieve RMS error of 3-6% compared to a commercial FEM package in lightly saturated conditions.

A purely nonlinear route to transition approaching the edge of chaos in a boundary layer

International Nuclear Information System (INIS)

Cherubini, S; De Palma, P; Robinet, J-Ch; Bottaro, A

2012-01-01

The understanding of transition in shear flows has recently progressed along new paradigms based on the central role of coherent flow structures and their nonlinear interactions. We follow such paradigms to identify, by means of a nonlinear optimization of the energy growth at short time, the initial perturbation which most easily induces transition in a boundary layer. Moreover, a bisection procedure has been used to identify localized flow structures living on the edge of chaos, found to be populated by hairpin vortices and streaks. Such an edge structure appears to act as a relative attractor for the trajectory of the laminar base state perturbed by the initial finite-amplitude disturbances, mediating the route to turbulence of the flow, via the triggering of a regeneration cycle of Λ and hairpin structures at different space and time scales. These findings introduce a new, purely nonlinear scenario of transition in a boundary-layer flow. (paper)

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.

Graph Theory-Based Technique for Isolating Corrupted Boundary Conditions in Continental-Scale River Network Hydrodynamic Simulation

Science.gov (United States)

Yu, C. W.; Hodges, B. R.; Liu, F.

2017-12-01

Development of continental-scale river network models creates challenges where the massive amount of boundary condition data encounters the sensitivity of a dynamic nu- merical model. The topographic data sets used to define the river channel characteristics may include either corrupt data or complex configurations that cause instabilities in a numerical solution of the Saint-Venant equations. For local-scale river models (e.g. HEC- RAS), modelers typically rely on past experience to make ad hoc boundary condition adjustments that ensure a stable solution - the proof of the adjustment is merely the sta- bility of the solution. To date, there do not exist any formal methodologies or automated procedures for a priori detecting/fixing boundary conditions that cause instabilities in a dynamic model. Formal methodologies for data screening and adjustment are a critical need for simulations with a large number of river reaches that draw their boundary con- dition data from a wide variety of sources. At the continental scale, we simply cannot assume that we will have access to river-channel cross-section data that has been ade- quately analyzed and processed. Herein, we argue that problematic boundary condition data for unsteady dynamic modeling can be identified through numerical modeling with the steady-state Saint-Venant equations. The fragility of numerical stability increases with the complexity of branching in river network system and instabilities (even in an unsteady solution) are typically triggered by the nonlinear advection term in Saint-Venant equations. It follows that the behavior of the simpler steady-state equations (which retain the nonlin- ear term) can be used to screen the boundary condition data for problematic regions. In this research, we propose a graph-theory based method to isolate the location of corrupted boundary condition data in a continental-scale river network and demonstrate its utility with a network of O(10^4) elements. Acknowledgement

Temperature jump boundary conditions in radiation diffusion

International Nuclear Information System (INIS)

Alonso, C.T.

1976-12-01

The radiation diffusion approximation greatly simplifies radiation transport problems. Yet the application of this method has often been unnecessarily restricted to optically thick regions, or has been extended through the use of such ad hoc devices as flux limiters. The purpose of this paper is to review and draw attention to the use of the more physically appropriate temperature jump boundary conditions for extending the range of validity of the diffusion approximation. Pioneering work has shown that temperature jump boundary conditions remove the singularity in flux that occurs in ordinary diffusion at small optical thicknesses. In this review paper Deissler's equations for frequency-dependent jump boundary conditions are presented and specific geometric examples are calculated analytically for steady state radiation transfer. When jump boundary conditions are applied to radiation diffusion, they yield exact solutions which are naturally flux- limited and geometry-corrected. We believe that the presence of temperature jumps on source boundaries is probably responsible in some cases for the past need for imposing ad hoc flux-limiting constraints on pure diffusion solutions. The solution for transfer between plane slabs, which is exact to all orders of optical thickness, also provides a useful tool for studying the accuracy of computer codes

Boundary condition histograms for modulated phases

International Nuclear Information System (INIS)

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

1997-11-01

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

Casimir pistons with general boundary conditions

Directory of Open Access Journals (Sweden)

Guglielmo Fucci

2015-02-01

Full Text Available In this work we analyze the Casimir energy and force for a scalar field endowed with general self-adjoint boundary conditions propagating in a higher dimensional piston configuration. The piston is constructed as a direct product I×N, with I=[0,L]⊂R and N a smooth, compact Riemannian manifold with or without boundary. The study of the Casimir energy and force for this configuration is performed by employing the spectral zeta function regularization technique. The obtained analytic results depend explicitly on the spectral zeta function associated with the manifold N and the parameters describing the general boundary conditions imposed. These results are then specialized to the case in which the manifold N is a d-dimensional sphere.

Modelling classroom conditions with different boundary conditions

DEFF Research Database (Denmark)

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

2014-01-01

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

Magnetohydrodynamics (MHD flow of a tangent hyperbolic fluid with nanoparticles past a stretching sheet with second order slip and convective boundary condition

Directory of Open Access Journals (Sweden)

Wubshet Ibrahim

Full Text Available This article presents the effect of thermal radiation on magnetohydrodynamic flow of tangent hyperbolic fluid with nanoparticle past an enlarging sheet with second order slip and convective boundary condition. Condition of zero normal flux of nanoparticles at the wall is used for the concentration boundary condition, which is the current topic that have yet to be studied extensively. The solution for the velocity, temperature and nanoparticle concentration is governed by parameters viz. power-law index (n, Weissenberg number We, Biot number Bi, Prandtl number Pr, velocity slip parameters δ and γ, Lewis number Le, Brownian motion parameter Nb and the thermophoresis parameter Nt. Similarity transformation is used to metamorphosed the governing non-linear boundary-value problem into coupled higher order non-linear ordinary differential equation. The succeeding equations were numerically solved using the function bvp4c from the matlab for different values of emerging parameters. Numerical results are deliberated through graphs and tables for velocity, temperature, concentration, the skin friction coefficient and local Nusselt number. The results designate that the skin friction coefficient Cf deplete as the values of Weissenberg number We, slip parameters γ and δ upturn and it rises as the values of power-law index n increase. The local Nusselt number -θ′(0 decreases as slip parameters γ and δ, radiation parameter Nr, Weissenberg number We, thermophoresis parameter Nt and power-law index n increase. However, the local Nusselt number increases as the Biot number Bi increase. Keywords: Tangent hyperbolic fluid, Second order slip flow, MHD, Convective boundary condition, Radiation effect, Passive control of nanoparticles

Hydromagnetic natural convection flow between vertical parallel plates with time-periodic boundary conditions

International Nuclear Information System (INIS)

Adesanya, S.O.; Oluwadare, E.O.; Falade, J.A.; Makinde, O.D.

2015-01-01

In this paper, the free convective flow of magnetohydrodynamic fluid through a channel with time periodic boundary condition is investigated by taking the effects of Joule dissipation into consideration. Based on simplifying assumptions, the coupled governing equations are reduced to a set of nonlinear boundary valued problem. Approximate solutions are obtained by using semi-analytical Adomian decomposition method. The effect of pertinent parameters on the fluid velocity, temperature distribution, Nusselt number and skin friction are presented graphically and discussed. The result of the computation shows that an increase in the magnetic field intensity has significant influence on the fluid flow. - Highlights: • The influence of magnetic field on the free convective fluid flow is considered. • The coupled equations are solved by using Adomian decomposition method. • The Adomian series solution agreed with previously obtained result. • Magnetic field decreases the velocity maximum but enhances temperature field

Absorbing boundary conditions for Einstein's field equations

Energy Technology Data Exchange (ETDEWEB)

Sarbach, Olivier [Instituto de Fisica y Matematicas, Universidad Michoacana de San Nicolas de Hidalgo, Edificio C-3, Cd. Universitaria. C. P. 58040 Morelia, Michoacan (Mexico)

2007-11-15

A common approach for the numerical simulation of wave propagation on a spatially unbounded domain is to truncate the domain via an artificial boundary, thus forming a finite computational domain with an outer boundary. Absorbing boundary conditions must then be specified at the boundary such that the resulting initial-boundary value problem is well posed and such that the amount of spurious reflection is minimized. In this article, we review recent results on the construction of absorbing boundary conditions in General Relativity and their application to numerical relativity.

Twin Positive Solutions of a Nonlinear m-Point Boundary Value Problem for Third-Order p-Laplacian Dynamic Equations on Time Scales

Directory of Open Access Journals (Sweden)

Wei Han

2008-01-01

Full Text Available Several existence theorems of twin positive solutions are established for a nonlinear m-point boundary value problem of third-order p-Laplacian dynamic equations on time scales by using a fixed point theorem. We present two theorems and four corollaries which generalize the results of related literature. As an application, an example to demonstrate our results is given. The obtained conditions are different from some known results.

Trickle-down boundary conditions in aeolian dune-field pattern formation

Science.gov (United States)

Ewing, R. C.; Kocurek, G.

2015-12-01

One the one hand, wind-blown dune-field patterns emerge within the overarching boundary conditions of climate, tectonics and eustasy implying the presence of these signals in the aeolian geomorphic and stratigraphic record. On the other hand, dune-field patterns are a poster-child of self-organization, in which autogenic processes give rise to patterned landscapes despite remarkable differences in the geologic setting (i.e., Earth, Mars and Titan). How important are climate, tectonics and eustasy in aeolian dune field pattern formation? Here we develop the hypothesis that, in terms of pattern development, dune fields evolve largely independent of the direct influence of 'system-scale' boundary conditions, such as climate, tectonics and eustasy. Rather, these boundary conditions set the stage for smaller-scale, faster-evolving 'event-scale' boundary conditions. This 'trickle-down' effect, in which system-scale boundary conditions indirectly influence the event scale boundary conditions provides the uniqueness and richness of dune-field patterned landscapes. The trickle-down effect means that the architecture of the stratigraphic record of dune-field pattern formation archives boundary conditions, which are spatially and temporally removed from the overarching geologic setting. In contrast, the presence of an aeolian stratigraphic record itself, reflects changes in system-scale boundary conditions that drive accumulation and preservation of aeolian strata.

Constraint-preserving boundary treatment for a harmonic formulation of the Einstein equations

Energy Technology Data Exchange (ETDEWEB)

Seiler, Jennifer; Szilagyi, Bela; Pollney, Denis; Rezzolla, Luciano [Max-Planck-Institut fuer Gravitationsphysik, Albert-Einstein-Institut, Golm (Germany)

2008-09-07

We present a set of well-posed constraint-preserving boundary conditions for a first-order in time, second-order in space, harmonic formulation of the Einstein equations. The boundary conditions are tested using robust stability, linear and nonlinear waves, and are found to be both less reflective and constraint preserving than standard Sommerfeld-type boundary conditions.

Constraint-preserving boundary treatment for a harmonic formulation of the Einstein equations

International Nuclear Information System (INIS)

Seiler, Jennifer; Szilagyi, Bela; Pollney, Denis; Rezzolla, Luciano

2008-01-01

We present a set of well-posed constraint-preserving boundary conditions for a first-order in time, second-order in space, harmonic formulation of the Einstein equations. The boundary conditions are tested using robust stability, linear and nonlinear waves, and are found to be both less reflective and constraint preserving than standard Sommerfeld-type boundary conditions

Integrable boundary conditions and modified Lax equations

International Nuclear Information System (INIS)

Avan, Jean; Doikou, Anastasia

2008-01-01

We consider integrable boundary conditions for both discrete and continuum classical integrable models. Local integrals of motion generated by the corresponding 'transfer' matrices give rise to time evolution equations for the initial Lax operator. We systematically identify the modified Lax pairs for both discrete and continuum boundary integrable models, depending on the classical r-matrix and the boundary matrix

Energy principle with included boundary conditions

International Nuclear Information System (INIS)

Lehnert, B.

1994-01-01

Earlier comments by the author on the limitations of the classical form of the extended energy principle are supported by a complementary analysis on the potential energy change arising from free-boundary displacements of a magnetically confined plasma. In the final formulation of the extended principle, restricted displacements, satisfying pressure continuity by means of plasma volume currents in a thin boundary layer, are replaced by unrestricted (arbitrary) displacements which can give rise to induced surface currents. It is found that these currents contribute to the change in potential energy, and that their contribution is not taken into account by such a formulation. A general expression is further given for surface currents induced by arbitrary displacements. The expression is used to reformulate the energy principle for the class of displacements which satisfy all necessary boundary conditions, including that of the pressure balance. This makes a minimization procedure of the potential energy possible, for the class of all physically relevant test functions which include the constraints imposed by the boundary conditions. Such a procedure is also consistent with a corresponding variational calculus. (Author)

Nonlinear Dynamics of Vortices in Different Types of Grain Boundaries

Science.gov (United States)

Sheikhzada, Ahmad K.

As a major component of linear particle accelerators, superconducting radio-frequency (SRF) resonator cavities are required to operate with lowest energy dissipation and highest accelerating gradient. SRF cavities are made of polycrystalline materials in which grain boundaries can limit maximum RF currents and produce additional power dissipation sources due to local penetration of Josephson vortices. The essential physics of vortex penetration and mechanisms of dissipation of vortices driven by strong RF currents along networks of grain boundaries and their contribution to the residual surface resistance have not been well understood. To evaluate how GBs can limit the performance of SRF materials, particularly Nb and Nb3Sn, we performed extensive numerical simulations of nonlinear dynamics of Josephson vortices in grain boundaries under strong dc and RF fields. The RF power due to penetration of vortices both in weakly-coupled and strongly-coupled grain boundaries was calculated as functions of the RF field and frequency. The result of this calculation manifested a quadratic dependence of power to field amplitude at strong RF currents, an illustration of resistive behavior of grain boundaries. Our calculations also showed that the surface resistance is a complicated function of field controlled by penetration and annihilation of vortices and antivortices in strong RF fields which ultimately saturates to normal resistivity of grain boundary. We found that Cherenkov radiation of rapidly moving vortices in grain boundaries can produce a new instability causing generation of expanding vortex-antivortex pair which ultimately drives the entire GB in a resistive state. This effect is more pronounced in polycrystalline thin film and multilayer coating structures in which it can cause significant increase in power dissipation and results in hysteresis effects in I-V characteristics, particularly at low temperatures.

Nonlinear Dynamics of Vortices in Different Types of Grain Boundaries

Energy Technology Data Exchange (ETDEWEB)

Sheikhzada, Ahmad [Old Dominion Univ., Norfolk, VA (United States)

2017-05-01

As a major component of linear particle accelerators, superconducting radio-frequency (SRF) resonator cavities are required to operate with lowest energy dissipation and highest accelerating gradient. SRF cavities are made of polycrystalline materials in which grain boundaries can limit maximum RF currents and produce additional power dissipation sources due to local penetration of Josephson vortices. The essential physics of vortex penetration and mechanisms of dissipation of vortices driven by strong RF currents along networks of grain boundaries and their contribution to the residual surface resistance have not been well understood. To evaluate how GBs can limit the performance of SRF materials, particularly Nb and Nb3Sn, we performed extensive numerical simulations of nonlinear dynamics of Josephson vortices in grain boundaries under strong dc and RF fields. The RF power due to penetration of vortices both in weakly-coupled and strongly-coupled grain boundaries was calculated as functions of the RF field and frequency. The result of this calculation manifested a quadratic dependence of power to field amplitude at strong RF currents, an illustration of resistive behavior of grain boundaries. Our calculations also showed that the surface resistance is a complicated function of field controlled by penetration and annihilation of vortices and antivortices in strong RF fields which ultimately saturates to normal resistivity of grain boundary. We found that Cherenkov radiation of rapidly moving vortices in grain boundaries can produce a new instability causing generation of expanding vortex-antivortex pair which ultimately drives the entire GB in a resistive state. This effect is more pronounced in polycrystalline thin film and multilayer coating structures in which it can cause significant increase in power dissipation and results in hysteresis effects in I-V characteristics, particularly at low temperatures.

Integrability and boundary conditions of supersymmetric systems

International Nuclear Information System (INIS)

Yue Ruihong; Liang Hong

1996-01-01

By studying the solutions of the reflection equations, we find out a series of integrable supersymmetric systems with different boundary conditions. The Hamiltonian contains four free parameters which describe the contribution of the boundary terms

Simulation of the steady-state energy transfer in rigid bodies, with convective-radiative boundary conditions, employing a minimum principle

International Nuclear Information System (INIS)

Gama, R.M.S. da.

1992-08-01

The energy transfer phenomenon in a rigid and opaque body that exchanges energy, with the environment, by convection and by diffuse thermal radiation is studied. The considered phenomenon is described by a partial differential equation, subjected to (nonlinear) boundary conditions. A minimum principle, suitable for a large class of energy transfer problems is presented. Some particular cases are simulated. (author)

Boundary fluxes for non-local diffusion

OpenAIRE

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

2006-01-01

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

A Parameter Estimation Method for Nonlinear Systems Based on Improved Boundary Chicken Swarm Optimization

Directory of Open Access Journals (Sweden)

Shaolong Chen

2016-01-01

Full Text Available Parameter estimation is an important problem in nonlinear system modeling and control. Through constructing an appropriate fitness function, parameter estimation of system could be converted to a multidimensional parameter optimization problem. As a novel swarm intelligence algorithm, chicken swarm optimization (CSO has attracted much attention owing to its good global convergence and robustness. In this paper, a method based on improved boundary chicken swarm optimization (IBCSO is proposed for parameter estimation of nonlinear systems, demonstrated and tested by Lorenz system and a coupling motor system. Furthermore, we have analyzed the influence of time series on the estimation accuracy. Computer simulation results show it is feasible and with desirable performance for parameter estimation of nonlinear systems.

Boundary conditions in rational conformal field theories

International Nuclear Information System (INIS)

Behrend, Roger E.; Pearce, Paul A.; Petkova, Valentina B.; Zuber, Jean-Bernard

2000-01-01

We develop further the theory of Rational Conformal Field Theories (RCFTs) on a cylinder with specified boundary conditions emphasizing the role of a triplet of algebras: the Verlinde, graph fusion and Pasquier algebras. We show that solving Cardy's equation, expressing consistency of a RCFT on a cylinder, is equivalent to finding integer valued matrix representations of the Verlinde algebra. These matrices allow us to naturally associate a graph G to each RCFT such that the conformal boundary conditions are labelled by the nodes of G. This approach is carried to completion for sl(2) theories leading to complete sets of conformal boundary conditions, their associated cylinder partition functions and the A-D-E classification. We also review the current status for WZW sl(3) theories. Finally, a systematic generalisation of the formalism of Cardy-Lewellen is developed to allow for multiplicities arising from more general representations of the Verlinde algebra. We obtain information on the bulk-boundary coefficients and reproduce the relevant algebraic structures from the sewing constraints

Einstein boundary conditions for the 3+1 Einstein equations

International Nuclear Information System (INIS)

Frittelli, Simonetta; Gomez, Roberto

2003-01-01

In the 3+1 framework of the Einstein equations for the case of a vanishing shift vector and arbitrary lapse, we calculate explicitly the four boundary equations arising from the vanishing of the projection of the Einstein tensor along the normal to the boundary surface of the initial-boundary value problem. Such conditions take the form of evolution equations along (as opposed to across) the boundary for certain components of the extrinsic curvature and for certain space derivatives of the three-metric. We argue that, in general, such boundary conditions do not follow necessarily from the evolution equations and the initial data, but need to be imposed on the boundary values of the fundamental variables. Using the Einstein-Christoffel formulation, which is strongly hyperbolic, we show how three of the boundary equations up to linear combinations should be used to prescribe the values of some incoming characteristic fields. Additionally, we show that the fourth one imposes conditions on some outgoing fields

On the Boussinesq-Burgers equations driven by dynamic boundary conditions

Science.gov (United States)

Zhu, Neng; Liu, Zhengrong; Zhao, Kun

2018-02-01

We study the qualitative behavior of the Boussinesq-Burgers equations on a finite interval subject to the Dirichlet type dynamic boundary conditions. Assuming H1 ×H2 initial data which are compatible with boundary conditions and utilizing energy methods, we show that under appropriate conditions on the dynamic boundary data, there exist unique global-in-time solutions to the initial-boundary value problem, and the solutions converge to the boundary data as time goes to infinity, regardless of the magnitude of the initial data.

Quadratic Functionals with General Boundary Conditions

International Nuclear Information System (INIS)

Dosla, Z.; Dosly, O.

1997-01-01

The purpose of this paper is to give the Reid 'Roundabout Theorem' for quadratic functionals with general boundary conditions. In particular, we describe the so-called coupled point and regularity condition introduced in terms of Riccati equation solutions

On filter boundary conditions in topology optimization

DEFF Research Database (Denmark)

Clausen, Anders; Andreassen, Erik

2017-01-01

Most research papers on topology optimization involve filters for regularization. Typically, boundary effects from the filters are ignored. Despite significant drawbacks the inappropriate homogeneous Neumann boundary conditions are used, probably because they are trivial to implement. In this paper...

On the instability of a 3-dimensional attachment line boundary layer: Weakly nonlinear theory and a numerical approach

Science.gov (United States)

Hall, P.; Malik, M. R.

1984-01-01

The instability of a three dimensional attachment line boundary layer is considered in the nonlinear regime. Using weakly nonlinear theory, it is found that, apart from a small interval near the (linear) critical Reynolds number, finite amplitude solutions bifurcate subcritically from the upper branch of the neutral curve. The time dependent Navier-Stokes equations for the attachment line flow have been solved using a Fourier-Chebyshev spectral method and the subcritical instability is found at wavenumbers that correspond to the upper branch. Both the theory and the numerical calculations show the existence of supercritical finite amplitude (equilibrium) states near the lower branch which explains why the observed flow exhibits a preference for the lower branch modes. The effect of blowing and suction on nonlinear stability of the attachment line boundary layer is also investigated.

On the instability of a three-dimensional attachment-line boundary layer - Weakly nonlinear theory and a numerical approach

Science.gov (United States)

Hall, P.; Malik, M. R.

1986-01-01

The instability of a three-dimensional attachment-line boundary layer is considered in the nonlinear regime. Using weakly nonlinear theory, it is found that, apart from a small interval near the (linear) critical Reynolds number, finite-amplitude solutions bifurcate subcritically from the upper branch of the neutral curve. The time-dependent Navier-Stokes equations for the attachment-line flow have been solved using a Fourier-Chebyshev spectral method and the subcritical instability is found at wavenumbers that correspond to the upper branch. Both the theory and the numerical calculations show the existence of supercritical finite-amplitude (equilibrium) states near the lower branch which explains why the observed flow exhibits a preference for the lower branch modes. The effect of blowing and suction on nonlinear stability of the attachment-line boundary layer is also investigated.

Boundary conditions for the gravitational field

International Nuclear Information System (INIS)

Winicour, Jeffrey

2012-01-01

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

Stability and boundary stabilization of 1-D hyperbolic systems

CERN Document Server

Bastin, Georges

2016-01-01

This monograph explores the modeling of conservation and balance laws of one-dimensional hyperbolic systems using partial differential equations. It presents typical examples of hyperbolic systems for a wide range of physical engineering applications, allowing readers to understand the concepts in whichever setting is most familiar to them. With these examples, it also illustrates how control boundary conditions may be defined for the most commonly used control devices. The authors begin with the simple case of systems of two linear conservation laws and then consider the stability of systems under more general boundary conditions that may be differential, nonlinear, or switching. They then extend their discussion to the case of nonlinear conservation laws and demonstrate the use of Lyapunov functions in this type of analysis. Systems of balance laws are considered next, starting with the linear variety before they move on to more general cases of nonlinear ones. They go on to show how the problem of boundary...

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.

Nonlinear hydromagnetic Rayleigh-Taylor instability for strong viscous fluids in porous media

CERN Document Server

El-Dib, Y O

2003-01-01

In the present work a weakly nonlinear stability for magnetic fluid is discussed. The research of an interface between two strong viscous homogeneous incompressible fluids through porous medium is investigated theoretically and graphically. The effect of the vertical magnetic field has been demonstrated in this study. The linear form of equation of motion is solved in the light of the nonlinear boundary conditions. The boundary value problem leads to construct nonlinear characteristic equation having complex coefficients in elevation function. The nonlinearity is kept to third-order expansion. The nonlinear characteristic equation leads to derive the well-known nonlinear Schroedinger equation. This equation having complex coefficients of the disturbance amplitude varies in both space and time. Stability criteria have been performed for nonlinear Chanderasekhar dispersion relation including the porous effects. Stability conditions are discussed through the assumption of equal kinematic viscosity. The calculati...

Advances in nonlinear vibration analysis of structures. Part-I. Beams

Indian Academy of Sciences (India)

Unknown

element analysis of nonlinear beams under static and dynamic loads. ... linearization, substitution of inplane boundary conditions at element level rather .... Modelling the nonlinear vibration problems using finite elements, albeit with a couple.

Convection and reaction in a diffusive boundary layer in a porous medium: nonlinear dynamics.

Science.gov (United States)

Andres, Jeanne Therese H; Cardoso, Silvana S S

2012-09-01

We study numerically the nonlinear interactions between chemical reaction and convective fingering in a diffusive boundary layer in a porous medium. The reaction enhances stability by consuming a solute that is unstably distributed in a gravitational field. We show that chemical reaction profoundly changes the dynamics of the system, by introducing a steady state, shortening the evolution time, and altering the spatial patterns of velocity and concentration of solute. In the presence of weak reaction, finger growth and merger occur effectively, driving strong convective currents in a thick layer of solute. However, as the reaction becomes stronger, finger growth is inhibited, tip-splitting is enhanced and the layer of solute becomes much thinner. Convection enhances the mass flux of solute consumed by reaction in the boundary layer but has a diminishing effect as reaction strength increases. This nonlinear behavior has striking differences to the density fingering of traveling reaction fronts, for which stronger chemical kinetics result in more effective finger merger owing to an increase in the speed of the front. In a boundary layer, a strong stabilizing effect of reaction can maintain a long-term state of convection in isolated fingers of wavelength comparable to that at onset of instability.

Asymptotic boundary conditions for dissipative waves: General theory

Science.gov (United States)

Hagstrom, Thomas

1990-01-01

An outstanding issue in the computational analysis of time dependent problems is the imposition of appropriate radiation boundary conditions at artificial boundaries. Accurate conditions are developed which are based on the asymptotic analysis of wave propagation over long ranges. Employing the method of steepest descents, dominant wave groups are identified and simple approximations to the dispersion relation are considered in order to derive local boundary operators. The existence of a small number of dominant wave groups may be expected for systems with dissipation. Estimates of the error as a function of domain size are derived under general hypotheses, leading to convergence results. Some practical aspects of the numerical construction of the asymptotic boundary operators are also discussed.

Asymptotic boundary conditions for dissipative waves - General theory

Science.gov (United States)

Hagstrom, Thomas

1991-01-01

An outstanding issue in computational analysis of time dependent problems is the imposition of appropriate radiation boundary conditions at artificial boundaries. Accurate conditions are developed which are based on the asymptotic analysis of wave propagation over long ranges. Employing the method of steepest descents, dominant wave groups are identified and simple approximations to the dispersion relation are considered in order to derive local boundary operators. The existence of a small number of dominant wave groups may be expected for systems with dissipation. Estimates of the error as a function of domain size are derived under general hypotheses, leading to convergence results. Some practical aspects of the numerical construction of the asymptotic boundary operators are also discussed.

Nonlinear model of a rotating hub-beams structure: Equations of motion

Science.gov (United States)

Warminski, Jerzy

2018-01-01

Dynamics of a rotating structure composed of a rigid hub and flexible beams is presented in the paper. A nonlinear model of a beam takes into account bending, extension and nonlinear curvature. The influence of geometric nonlinearity and nonconstant angular velocity on dynamics of the rotating structure is presented. The exact equations of motion and associated boundary conditions are derived on the basis of the Hamilton's principle. The simplification of the exact nonlinear mathematical model is proposed taking into account the second order approximation. The reduced partial differential equations of motion together with associated boundary conditions can be used to study natural or forced vibrations of a rotating structure considering constant or nonconstant angular speed of a rigid hub and an arbitrary number of flexible blades.

A numerical study of magnetohydrodynamics flow in Casson nanofluid combined with Joule heating and slip boundary conditions

Directory of Open Access Journals (Sweden)

A. Kamran

Full Text Available A numerical study of Casson nanofluid past horizontal stretching surface with magnetic effect and Joule heating are presented. Slip and thermal convective boundary conditions are considered in the study. A numerical technique of Keller box is applied to the nonlinear ODEs which are obtained by applying the similarity transformation to the nonlinear partial differential equations. The magnetic field and Joule heating effects are observed graphically. Also the strength of convective heat exchange (Nusselt number and the strength of mass exchange (Sherwood number are analyzed. It is noted that Nusselt number declines whereas Sherwood number rises by increasing Eckert number. The impact of increasing Hartman number resulted in the decrease of both Sherwood and Nusselt number. Keywords: Casson nanofluid, Magnetohydrodynamic, Joule heating, Keller box method

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.

Exponential convergence for nonlinear diffusion problems with positive lateral boundary conditions

International Nuclear Information System (INIS)

Holland, C.J.; Berryman, J.G.

1985-01-01

It is established that the solution u of u/sub t/ = Δ(u/sup m/)>0, with positive initial data, positive lateral boundary data, and positive exponent m, converges exponentially to the solution v of the corresponding stationary equation Δ(v/sup m/) = 0. The analysis also provides the form of the leading contribution to the difference

Nonlinear core deflection in injection molding

Science.gov (United States)

Poungthong, P.; Giacomin, A. J.; Saengow, C.; Kolitawong, C.; Liao, H.-C.; Tseng, S.-C.

2018-05-01

Injection molding of thin slender parts is often complicated by core deflection. This deflection is caused by molten plastics race tracking through the slit between the core and the rigid cavity wall. The pressure of this liquid exerts a lateral force of the slender core causing the core to bend, and this bending is governed by a nonlinear fifth order ordinary differential equation for the deflection that is not directly in the position along the core. Here we subject this differential equation to 6 sets of boundary conditions, corresponding to 6 commercial core constraints. For each such set of boundary conditions, we develop an explicit approximate analytical solution, including both a linear term and a nonlinear term. By comparison with finite difference solutions, we find our new analytical solutions to be accurate. We then use these solutions to derive explicit analytical approximations for maximum deflections and for the core position of these maximum deflections. Our experiments on the base-gated free-tip boundary condition agree closely with our new explicit approximate analytical solution.

Boundary conditions and subelliptic estimates for geometric Kramers-Fokker-Planck operators on manifolds with boundaries

CERN Document Server

Nier, Francis

2018-01-01

This article is concerned with the maximal accretive realizations of geometric Kramers-Fokker-Planck operators on manifolds with boundaries. A general class of boundary conditions is introduced which ensures the maximal accretivity and some global subelliptic estimates. Those estimates imply nice spectral properties as well as exponential decay properties for the associated semigroup. Admissible boundary conditions cover a wide range of applications for the usual scalar Kramer-Fokker-Planck equation or Bismut's hypoelliptic laplacian.

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

Minimization of heat slab nodes with higher order boundary conditions

International Nuclear Information System (INIS)

Solbrig, C.W.

1992-01-01

The accuracy of a numerical solution can be limited by the numerical approximation to the boundary conditions rather than the accuracy of the equations which describe the interior. The study presented in this paper compares the results from two different numerical formulations of the convective boundary condition on the face of a heat transfer slab. The standard representation of the boundary condition in a test problem yielded an unacceptable error even when the heat transfer slab was partitioned into over 300 nodes. A higher order boundary condition representation was obtained by using a second order approximation for the first derivative at the boundary and combining it with the general equation used for inner nodes. This latter formulation produced reasonable results when as few as ten nodes were used

Transport synthetic acceleration with opposing reflecting boundary conditions

Energy Technology Data Exchange (ETDEWEB)

Zika, M R; Adams, M L

2000-02-01

The transport synthetic acceleration (TSA) scheme is extended to problems with opposing reflecting boundary conditions. This synthetic method employs a simplified transport operator as its low-order approximation. A procedure is developed that allows the use of the conjugate gradient (CG) method to solve the resulting low-order system of equations. Several well-known transport iteration algorithms are cast in a linear algebraic form to show their equivalence to standard iterative techniques. Source iteration in the presence of opposing reflecting boundary conditions is shown to be equivalent to a (poorly) preconditioned stationary Richardson iteration, with the preconditioner defined by the method of iterating on the incident fluxes on the reflecting boundaries. The TSA method (and any synthetic method) amounts to a further preconditioning of the Richardson iteration. The presence of opposing reflecting boundary conditions requires special consideration when developing a procedure to realize the CG method for the proposed system of equations. The CG iteration may be applied only to symmetric positive definite matrices; this condition requires the algebraic elimination of the boundary angular corrections from the low-order equations. As a consequence of this elimination, evaluating the action of the resulting matrix on an arbitrary vector involves two transport sweeps and a transmission iteration. Results of applying the acceleration scheme to a simple test problem are presented.

Revisit boundary conditions for the self-adjoint angular flux formulation

Energy Technology Data Exchange (ETDEWEB)

Wang, Yaqi [Idaho National Lab. (INL), Idaho Falls, ID (United States); Gleicher, Frederick N. [Idaho National Lab. (INL), Idaho Falls, ID (United States)

2015-03-01

We revisit the boundary conditions for SAAF. We derived the equivalent parity variational form ready for coding up. The more rigorous approach of evaluating odd parity should be solving the odd parity equation coupled with the even parity. We proposed a symmetric reflecting boundary condition although neither positive definiteness nor even-odd decoupling is achieved. A simple numerical test verifies the validity of these boundary conditions.

Nonlinear Cointegration Approach for Condition Monitoring of Wind Turbines

Directory of Open Access Journals (Sweden)

Konrad Zolna

2015-01-01

Full Text Available Monitoring of trends and removal of undesired trends from operational/process parameters in wind turbines is important for their condition monitoring. This paper presents the homoscedastic nonlinear cointegration for the solution to this problem. The cointegration approach used leads to stable variances in cointegration residuals. The adapted Breusch-Pagan test procedure is developed to test for the presence of heteroscedasticity in cointegration residuals obtained from the nonlinear cointegration analysis. Examples using three different time series data sets—that is, one with a nonlinear quadratic deterministic trend, another with a nonlinear exponential deterministic trend, and experimental data from a wind turbine drivetrain—are used to illustrate the method and demonstrate possible practical applications. The results show that the proposed approach can be used for effective removal of nonlinear trends form various types of data, allowing for possible condition monitoring applications.

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

Evaluation of general non-reflecting boundary conditions for industrial CFD applications

Science.gov (United States)

Basara, Branislav; Frolov, Sergei; Lidskii, Boris; Posvyanskii, Vladimir

2007-11-01

The importance of having proper boundary conditions for the calculation domain is a known issue in Computational Fluid Dynamics (CFD). In many situations, it is very difficult to define a correct boundary condition. The flow may enter and leave the computational domain at the same time and at the same boundary. In such circumstances, it is important that numerical implementation of boundary conditions enforces certain physical constraints leading to correct results which then ensures a better convergence rate. The aim of this paper is to evaluate recently proposed non-reflecting boundary conditions (Frolov et al., 2001, Advances in Chemical Propulsion) on industrial CFD applications. Derivation of the local non-reflecting boundary conditions at the open boundary is based on finding the solution of linearized Euler equations vanishing at infinity for both incompressible and compressible formulations. This is implemented into the in-house CFD package AVL FIRE and some numerical details will be presented as well. The key applications in this paper are from automotive industry, e.g. an external car aerodynamics, an intake port, etc. The results will show benefits of using effective non-reflecting boundary conditions.

Nonlinear Displacement Discontinuity Model for Generalized Rayleigh Wave in Contact Interface

Energy Technology Data Exchange (ETDEWEB)

Kim, No Hyu; Yang, Seung Yong [Korea University of Technology and Education, Cheonan (Korea, Republic of)

2007-12-15

Imperfectly jointed interface serves as mechanical waveguide for elastic waves and gives rise to two distinct kinds of guided wave propagating along the interface. Contact acoustic nonlinearity (CAN) is known to plays major role in the generation of these interface waves called generalized Rayleigh waves in non-welded interface. Closed crack is modeled as non-welded interface that has nonlinear discontinuity condition in displacement across its boundary. Mathematical analysis of boundary conditions and wave equation is conducted to investigate the dispersive characteristics of the interface waves. Existence of the generalized Rayleigh wave(interface wave) in nonlinear contact interface is verified in theory where the dispersion equation for the interface wave is formulated and analyzed. It reveals that the interface waves have two distinct modes and that the phase velocity of anti-symmetric wave mode is highly dependent on contact conditions represented by linear and nonlinear dimensionless specific stiffness

Nonlinear Displacement Discontinuity Model for Generalized Rayleigh Wave in Contact Interface

International Nuclear Information System (INIS)

Kim, No Hyu; Yang, Seung Yong

2007-01-01

Imperfectly jointed interface serves as mechanical waveguide for elastic waves and gives rise to two distinct kinds of guided wave propagating along the interface. Contact acoustic nonlinearity (CAN) is known to plays major role in the generation of these interface waves called generalized Rayleigh waves in non-welded interface. Closed crack is modeled as non-welded interface that has nonlinear discontinuity condition in displacement across its boundary. Mathematical analysis of boundary conditions and wave equation is conducted to investigate the dispersive characteristics of the interface waves. Existence of the generalized Rayleigh wave(interface wave) in nonlinear contact interface is verified in theory where the dispersion equation for the interface wave is formulated and analyzed. It reveals that the interface waves have two distinct modes and that the phase velocity of anti-symmetric wave mode is highly dependent on contact conditions represented by linear and nonlinear dimensionless specific stiffness

Thermal Simulations, Open Boundary Conditions and Switches

Science.gov (United States)

Burnier, Yannis; Florio, Adrien; Kaczmarek, Olaf; Mazur, Lukas

2018-03-01

SU(N) gauge theories on compact spaces have a non-trivial vacuum structure characterized by a countable set of topological sectors and their topological charge. In lattice simulations, every topological sector needs to be explored a number of times which reflects its weight in the path integral. Current lattice simulations are impeded by the so-called freezing of the topological charge problem. As the continuum is approached, energy barriers between topological sectors become well defined and the simulations get trapped in a given sector. A possible way out was introduced by Lüscher and Schaefer using open boundary condition in the time extent. However, this solution cannot be used for thermal simulations, where the time direction is required to be periodic. In this proceedings, we present results obtained using open boundary conditions in space, at non-zero temperature. With these conditions, the topological charge is not quantized and the topological barriers are lifted. A downside of this method are the strong finite-size effects introduced by the boundary conditions. We also present some exploratory results which show how these conditions could be used on an algorithmic level to reshuffle the system and generate periodic configurations with non-zero topological charge.

Thermal Simulations, Open Boundary Conditions and Switches

Directory of Open Access Journals (Sweden)

Burnier Yannis

2018-01-01

Full Text Available SU(N gauge theories on compact spaces have a non-trivial vacuum structure characterized by a countable set of topological sectors and their topological charge. In lattice simulations, every topological sector needs to be explored a number of times which reflects its weight in the path integral. Current lattice simulations are impeded by the so-called freezing of the topological charge problem. As the continuum is approached, energy barriers between topological sectors become well defined and the simulations get trapped in a given sector. A possible way out was introduced by Lüscher and Schaefer using open boundary condition in the time extent. However, this solution cannot be used for thermal simulations, where the time direction is required to be periodic. In this proceedings, we present results obtained using open boundary conditions in space, at non-zero temperature. With these conditions, the topological charge is not quantized and the topological barriers are lifted. A downside of this method are the strong finite-size effects introduced by the boundary conditions. We also present some exploratory results which show how these conditions could be used on an algorithmic level to reshuffle the system and generate periodic configurations with non-zero topological charge.

Unified treatment of microscopic boundary conditions and efficient algorithms for estimating tangent operators of the homogenized behavior in the computational homogenization method

Science.gov (United States)

Nguyen, Van-Dung; Wu, Ling; Noels, Ludovic

2017-03-01

This work provides a unified treatment of arbitrary kinds of microscopic boundary conditions usually considered in the multi-scale computational homogenization method for nonlinear multi-physics problems. An efficient procedure is developed to enforce the multi-point linear constraints arising from the microscopic boundary condition either by the direct constraint elimination or by the Lagrange multiplier elimination methods. The macroscopic tangent operators are computed in an efficient way from a multiple right hand sides linear system whose left hand side matrix is the stiffness matrix of the microscopic linearized system at the converged solution. The number of vectors at the right hand side is equal to the number of the macroscopic kinematic variables used to formulate the microscopic boundary condition. As the resolution of the microscopic linearized system often follows a direct factorization procedure, the computation of the macroscopic tangent operators is then performed using this factorized matrix at a reduced computational time.

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

Open boundary condition, Wilson flow and the scalar glueball mass

International Nuclear Information System (INIS)

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

2014-01-01

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

Nonlinear observer-based Lyapunov boundary control of distributed heat transfer mechanisms for membrane distillation plant

KAUST Repository

Eleiwi, Fadi

2016-09-19

This paper presents a nonlinear observer-based Lyapunov control for a membrane distillation (MD) process. The control considers the inlet temperatures of the feed and the permeate solutions as inputs, transforming it to boundary control process, and seeks to maintain the temperature difference along the membrane boundaries around a sufficient level to promote water production. MD process is modeled with advection diffusion equation model in two dimensions, where the diffusion and convection heat transfer mechanisms are best described. Model analysis, effective order reduction and parameters physical interpretation, are provided. Moreover, a nonlinear observer has been designed to provide the control with estimates of the temperature evolution at each time instant. In addition, physical constraints are imposed on the control to have an acceptable range of feasible inputs, and consequently, better energy consumption. Numerical simulations for the complete process with real membrane parameter values are provided, in addition to detailed explanations for the role of the controller and the observer. (C) 2016 Elsevier Ltd. All rights reserved.

Entropy Stability and the No-Slip Wall Boundary Condition

KAUST Repository

Svä rd, Magnus; Carpenter, Mark H.; Parsani, Matteo

2018-01-01

We present an entropy stable numerical scheme subject to no-slip wall boundary conditions. To enforce entropy stability only the no-penetration boundary condition and a temperature condition are needed at a wall, and this leads to an L bound on the conservative variables. In this article, we take the next step and design a finite difference scheme that also bounds the velocity gradients. This necessitates the use of the full no-slip conditions.

Entropy Stability and the No-Slip Wall Boundary Condition

KAUST Repository

Svärd, Magnus

2018-01-18

We present an entropy stable numerical scheme subject to no-slip wall boundary conditions. To enforce entropy stability only the no-penetration boundary condition and a temperature condition are needed at a wall, and this leads to an L bound on the conservative variables. In this article, we take the next step and design a finite difference scheme that also bounds the velocity gradients. This necessitates the use of the full no-slip conditions.

Exploring exotic states with twisted boundary conditions

International Nuclear Information System (INIS)

Agadjanov, Dimitri

2017-01-01

he goal of this thesis is to develop methods to study the nature and properties of exotic hadrons from lattice simulations. The main focus lies in the application of twisted boundary conditions. The thesis consists of a general introduction and the collection of three papers, represented respectively in three chapters. The introduction of the thesis reviews the theoretical background, which is further used in the rest of the thesis. Further implementing partially twisted boundary conditions in the scalar sector of lattice QCD is studied. Then we develop a method to study the content of the exotic hadrons by determining the wave function renormalization constant from lattice simulations, exploiting the dependence of the spectrum on the twisted boundary conditions. The final chapter deals with a novel method to study the multi-channel scattering problem in a finite volume, which is relevant for exotic states. Its key idea is to extract the complex hadron-hadron optical potential, avoiding the difficulties, associated with the solution of the multi-channel Luescher equation.

Exploring exotic states with twisted boundary conditions

Energy Technology Data Exchange (ETDEWEB)

Agadjanov, Dimitri

2017-09-11

he goal of this thesis is to develop methods to study the nature and properties of exotic hadrons from lattice simulations. The main focus lies in the application of twisted boundary conditions. The thesis consists of a general introduction and the collection of three papers, represented respectively in three chapters. The introduction of the thesis reviews the theoretical background, which is further used in the rest of the thesis. Further implementing partially twisted boundary conditions in the scalar sector of lattice QCD is studied. Then we develop a method to study the content of the exotic hadrons by determining the wave function renormalization constant from lattice simulations, exploiting the dependence of the spectrum on the twisted boundary conditions. The final chapter deals with a novel method to study the multi-channel scattering problem in a finite volume, which is relevant for exotic states. Its key idea is to extract the complex hadron-hadron optical potential, avoiding the difficulties, associated with the solution of the multi-channel Luescher equation.

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

Divergence-Free Wavelets on the Hypercube : General Boundary Conditions

NARCIS (Netherlands)

Stevenson, R.

2016-01-01

On the n-dimensional hypercube, for given k∈N, wavelet Riesz bases are constructed for the subspace of divergence-free vector fields of the Sobolev space Hk((0,1)n)n with general homogeneous Dirichlet boundary conditions, including slip or no-slip boundary conditions. Both primal and suitable dual

Spectral asymmetry for bag boundary conditions

International Nuclear Information System (INIS)

Beneventano, C G; Santangelo, E M; Wipf, A

2002-01-01

We give an expression, in terms of boundary spectral functions, for the spectral asymmetry of the Euclidean Dirac operator in two dimensions, when its domain is determined by local boundary conditions and the manifold is of product type. As an application, we explicitly evaluate the asymmetry in the case of a finite-length cylinder and check that the outcome is consistent with our general result. Finally, we study the asymmetry in a disc, which is a non-product case, and propose an interpretation

Nonlinear interaction of s-polarized surface waves at the boundary of a semibounded magnetized plasma

International Nuclear Information System (INIS)

Amein, W.H.; El-Siragy, N.M.; Nagy, O.Z.; Sayed, Y.A.

1981-01-01

Nonlinear interaction of S-Polarized surface waves at the boundary of a semibounded magnetized plasma is investigated. The expressions of the amplitudes of the generated waves are found. It is shown that, the generated waves with combined frequencies are equally radiated from the transient layer into plasma and vacuum

Diffusive growth of a single droplet with three different boundary conditions

Science.gov (United States)

Tavassoli, Z.; Rodgers, G. J.

2000-02-01

We study a single, motionless three-dimensional droplet growing by adsorption of diffusing monomers on a 2D substrate. The diffusing monomers are adsorbed at the aggregate perimeter of the droplet with different boundary conditions. Models with both an adsorption boundary condition and a radiation boundary condition, as well as a phenomenological model, are considered and solved in a quasistatic approximation. The latter two models allow particle detachment. In the short time limit, the droplet radius grows as a power of the time with exponents of 1/4, 1/2 and 3/4 for the models with adsorption, radiation and phenomenological boundary conditions, respectively. In the long time limit a universal growth rate as $[t/\\ln(t)]^{1/3}$ is observed for the radius of the droplet for all models independent of the boundary conditions. This asymptotic behaviour was obtained by Krapivsky \\cite{krapquasi} where a similarity variable approach was used to treat the growth of a droplet with an adsorption boundary condition based on a quasistatic approximation. Another boundary condition with a constant flux of monomers at the aggregate perimeter is also examined. The results exhibit a power law growth rate with an exponent of 1/3 for all times.

Nonlinear Elliptic Differential Equations with Multivalued Nonlinearities

Indian Academy of Sciences (India)

In this paper we study nonlinear elliptic boundary value problems with monotone and nonmonotone multivalued nonlinearities. First we consider the case of monotone nonlinearities. In the first result we assume that the multivalued nonlinearity is defined on all R R . Assuming the existence of an upper and of a lower ...

A suitable boundary condition for bounded plasma simulation without sheath resolution

International Nuclear Information System (INIS)

Parker, S.E.; Procassini, R.J.; Birdsall, C.K.; Cohen, B.I.

1993-01-01

We have developed a technique that allows for a sheath boundary layer without having to resolve the inherently small space and time scales of the sheath region. We refer to this technique as the logical sheath boundary condition. This boundary condition, when incorporated into a direct-implicit particle code, permits large space- and time-scale simulations of bounded systems, which would otherwise be impractical on current supercomputers. The lack of resolution of the collector sheath potential drop obtained from conventional implicit simulations at moderate values of ω pe Δt and Δz/λ De provides the motivation for the development of the logical sheath boundary condition. The algorithm for use of the logical sheath boundary condition in a particle simulation is presented. Results from simulations which use the logical sheath boundary condition are shown to compare reasonably well with those from an analytic theory and simulations in which the sheath is resolved

Nested Bethe Ansatz for Spin Ladder Model with Open Boundary Conditions

International Nuclear Information System (INIS)

Wu Junfang; Zhang Chunmin; Yue Ruihong; Li Runling

2005-01-01

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

Quasi-periodic solutions of nonlinear beam equations with quintic quasi-periodic nonlinearities

Directory of Open Access Journals (Sweden)

Qiuju Tuo

2015-01-01

Full Text Available In this article, we consider the one-dimensional nonlinear beam equations with quasi-periodic quintic nonlinearities $$ u_{tt}+u_{xxxx}+(B+ \\varepsilon\\phi(tu^5=0 $$ under periodic boundary conditions, where B is a positive constant, $\\varepsilon$ is a small positive parameter, $\\phi(t$ is a real analytic quasi-periodic function in t with frequency vector $\\omega=(\\omega_1,\\omega_2,\\dots,\\omega_m$. It is proved that the above equation admits many quasi-periodic solutions by KAM theory and partial Birkhoff normal form.

Effective Stress Law in Unconventional Reservoirs under Different Boundary Conditions

Science.gov (United States)

Saurabh, S.; Harpalani, S.

2017-12-01

Unconventional reservoirs have attracted a great deal of research interest worldwide during the past two decades. Low permeability and specialized techniques required to exploit these resources present opportunities for improvement in both production rates and ultimate recovery. Understanding subsurface stress modifications and permeability evolution are valuable when evaluating the prospects of unconventional reservoirs. These reservoir properties are functions of effective stress. As a part of this study, effective stress law, specifically the variation of anisotropic Biot's coefficient under various boundary conditions believed to exist in gas reservoirs by different researchers, has been established. Pressure-dependent-permeability (PdK) experiments were carried out on San Juan coal under different boundary conditions, that is, uniaxial strain condition and constant volume condition. Stress and strain in the vertical and horizontal directions were monitored throughout the experiment. Data collected during the experiments was used to determine the Biot's coefficient in vertical and horizontal directions under these two boundary conditions, treating coal as transversely isotropic. The variation of Biot's coefficient was found to be well correlated with the variation in coal permeability. Based on the estimated values of Biot's coefficients, a theory of variation in its value is presented for other boundary conditions. The findings of the study shed light on the inherent behavior of Biot's coefficient under different reservoir boundary conditions. This knowledge can improve the modeling work requiring estimation of effective stress in reservoirs, such as, pressure-/stress- dependent permeability. At the same time, if the effective stresses are known with more certainty by other methods, it enables assessment of the unknown reservoir boundary conditions.

On a strong solution of the non-stationary Navier-Stokes equations under slip or leak boundary conditions of friction type

Science.gov (United States)

Kashiwabara, Takahito

Strong solutions of the non-stationary Navier-Stokes equations under non-linearized slip or leak boundary conditions are investigated. We show that the problems are formulated by a variational inequality of parabolic type, to which uniqueness is established. Using Galerkin's method and deriving a priori estimates, we prove global and local existence for 2D and 3D slip problems respectively. For leak problems, under no-leak assumption at t=0 we prove local existence in 2D and 3D cases. Compatibility conditions for initial states play a significant role in the estimates.

An(1) affine Toda field theories with integrable boundary conditions revisited

International Nuclear Information System (INIS)

Doikou, Anastasia

2008-01-01

Generic classically integrable boundary conditions for the A n (1) affine Toda field theories (ATFT) are investigated. The present analysis rests primarily on the underlying algebra, defined by the classical version of the reflection equation. We use as a prototype example the first non-trivial model of the hierarchy i.e. the A 2 (1) ATFT, however our results may be generalized for any A n (1) (n > 1). We assume here two distinct types of boundary conditions called some times soliton preserving (SP), and soliton non-preserving (SNP) associated to two distinct algebras, i.e. the reflection algebra and the (q) twisted Yangian respectively. The boundary local integrals of motion are then systematically extracted from the asymptotic expansion of the associated transfer matrix. In the case of SNP boundary conditions we recover previously known results. The other type of boundary conditions (SP), associated to the reflection algebra, are novel in this context and lead to a different set of conserved quantities that depend on free boundary parameters. It also turns out that the number of local integrals of motion for SP boundary conditions is 'double' compared to those of the SNP case.

Periodic oscillations in linear continuous media coupled with nonlinear discrete systems

International Nuclear Information System (INIS)

Lupini, R.

1998-01-01

A general derivation of partial differential equations with boundary conditions in the form of ordinary differential equations is obtained using the principle of stationary action for a Lagrangian function composed of continuous plus discrete parts in interaction across the boundaries of a 1-dimensional medium. This approach leads directly to the theorem of energy conservation. For linear continuous medium, homogeneous Dirichlet condition at one boundary, and nonlinear oscillator at the other boundary, the entire differential problem reduces to a nonlinear differential-difference equation of neutral type and of the second order. The lag parameter is τ = l/c, where c is the phase speed, l the length of the continuum. The Author investigate the problem of the occurrence of periodic solutions of period integer multiple of the lag (super harmonic solutions) in the case of zero inertia of the boundary system. The problem for such oscillations is shown to reduce to systems of ordinary differential equations with matching conditions in a phase space of lower dimensionality: Phase-plane techniques are used to determine solutions of period 4τ, 8τ and 6τ

On boundary conditions in three-dimensional AdS gravity

Energy Technology Data Exchange (ETDEWEB)

Miskovic, Olivera [Instituto de Fisica, P. Universidad Catolica de Valparaiso, Casilla 4059, Valparaiso (Chile) and Departamento de Fisica, P. Universidad Catolica de Chile, Casilla 306, Santiago 22 (Chile)]. E-mail: olivera.miskovic@ucv.cl; Olea, Rodrigo [Departamento de Fisica, P. Universidad Catolica de Chile, Casilla 306, Santiago 22 (Chile) and Centro Multidisciplinar de Astrofisica, CENTRA, Departamento de Fisica, Instituto Superior Tecnico, Universidade Tecnica de Lisboa, Av. Rovisco Pais 1, 1049-001 Lisbon (Portugal)]. E-mail: rolea@fisica.ist.utl.pt

2006-09-07

A finite action principle for three-dimensional gravity with negative cosmological constant, based on a boundary condition for the asymptotic extrinsic curvature, is considered. The bulk action appears naturally supplemented by a boundary term that is one half the Gibbons-Hawking term, that makes the Euclidean action and the Noether charges finite without additional Dirichlet counterterms. The consistency of this boundary condition with the Dirichlet problem in AdS gravity and the Chern-Simons formulation in three dimensions, and its suitability for the higher odd-dimensional case, are also discussed.

Homogenization of the stochastic Navier–Stokes equation with a stochastic slip boundary condition

KAUST Repository

Bessaih, Hakima

2015-11-02

The two-dimensional Navier–Stokes equation in a perforated domain with a dynamical slip boundary condition is considered. We assume that the dynamic is driven by a stochastic perturbation on the interior of the domain and another stochastic perturbation on the boundaries of the holes. We consider a scaling (ᵋ for the viscosity and 1 for the density) that will lead to a time-dependent limit problem. However, the noncritical scaling (ᵋ, β > 1) is considered in front of the nonlinear term. The homogenized system in the limit is obtained as a Darcy’s law with memory with two permeabilities and an extra term that is due to the stochastic perturbation on the boundary of the holes. The nonhomogeneity on the boundary contains a stochastic part that yields in the limit an additional term in the Darcy’s law. We use the two-scale convergence method after extending the solution with 0 inside the holes to pass to the limit. By Itô stochastic calculus, we get uniform estimates on the solution in appropriate spaces. Due to the stochastic integral, the pressure that appears in the variational formulation does not have enough regularity in time. This fact made us rely only on the variational formulation for the passage to the limit on the solution. We obtain a variational formulation for the limit that is solution of a Stokes system with two pressures. This two-scale limit gives rise to three cell problems, two of them give the permeabilities while the third one gives an extra term in the Darcy’s law due to the stochastic perturbation on the boundary of the holes.

Automated Boundary Conditions for Wind Tunnel Simulations

Science.gov (United States)

Carlson, Jan-Renee

2018-01-01

Computational fluid dynamic (CFD) simulations of models tested in wind tunnels require a high level of fidelity and accuracy particularly for the purposes of CFD validation efforts. Considerable effort is required to ensure the proper characterization of both the physical geometry of the wind tunnel and recreating the correct flow conditions inside the wind tunnel. The typical trial-and-error effort used for determining the boundary condition values for a particular tunnel configuration are time and computer resource intensive. This paper describes a method for calculating and updating the back pressure boundary condition in wind tunnel simulations by using a proportional-integral-derivative controller. The controller methodology and equations are discussed, and simulations using the controller to set a tunnel Mach number in the NASA Langley 14- by 22-Foot Subsonic Tunnel are demonstrated.

Nonlinear Kirchhoff-Carrier wave equation in a unit membrane with mixed homogeneous boundary conditions

Directory of Open Access Journals (Sweden)

Nguyen Thanh Long

2005-12-01

Full Text Available In this paper we consider the nonlinear wave equation problem $$displaylines{ u_{tt}-Big(|u|_0^2,|u_{r}|_0^2ig(u_{rr}+frac{1}{r}u_{r} =f(r,t,u,u_{r},quad 0less than r less than 1,; 0 less than t less than T, ig|lim_{ro 0^+}sqrt{r}u_{r}(r,tig| less than infty, u_{r}(1,t+hu(1,t=0, u(r,0=widetilde{u}_0(r, u_{t}(r,0=widetilde{u}_1(r. }$$ To this problem, we associate a linear recursive scheme for which the existence of a local and unique weak solution is proved, in weighted Sobolev using standard compactness arguments. In the latter part, we give sufficient conditions for quadratic convergence to the solution of the original problem, for an autonomous right-hand side independent on $u_{r}$ and a coefficient function $B$ of the form $B=B(|u|_0^2=b_0+|u|_0^2$ with $b_0$ greater than 0.

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

International Nuclear Information System (INIS)

Koteras, J.R.

1996-01-01

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

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

International Nuclear Information System (INIS)

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

1982-01-01

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

On domain wall boundary conditions for the XXZ spin Hamiltonian

DEFF Research Database (Denmark)

Orlando, Domenico; Reffert, Susanne; Reshetikhin, Nicolai

In this note, we derive the spectrum of the infinite quantum XXZ spin chain with domain wall boundary conditions. The eigenstates are constructed as limits of Bethe states for the finite XXZ spin chain with quantum sl(2) invariant boundary conditions....

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.

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

Directory of Open Access Journals (Sweden)

Nahed S. Hussein

2014-01-01

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

Simulations of QCD and QED with C* boundary conditions

Science.gov (United States)

Hansen, Martin; Lucini, Biagio; Patella, Agostino; Tantalo, Nazario

2018-03-01

We present exploratory results from dynamical simulations of QCD in isolation, as well as QCD coupled to QED, with C* boundary conditions. In finite volume, the use of C* boundary conditions allows for a gauge invariant and local formulation of QED without zero modes. In particular we show that the simulations reproduce known results and that masses of charged mesons can be extracted in a completely gauge invariant way. For the simulations we use a modified version of the HiRep code. The primary features of the simulation code are presented and we discuss some details regarding the implementation of C* boundary conditions and the simulated lattice action. Preprint: CP3-Origins-2017-046 DNRF90, CERN-TH-2017-214

Critical effects of downstream boundary conditions on vortex breakdown

Science.gov (United States)

Kandil, Osama; Kandil, Hamdy A.; Liu, C. H.

1992-01-01

The unsteady, compressible, full Navier-Stokes (NS) equations are used to study the critical effects of the downstream boundary conditions on the supersonic vortex breakdown. The present study is applied to two supersonic vortex breakdown cases. In the first case, quasi-axisymmetric supersonic swirling flow is considered in a configured circular duct, and in the second case, quasi-axisymmetric supersonic swirling jet, that is issued from a nozzle into a supersonic jet of lower Mach number, is considered. For the configured duct flow, four different types of downstream boundary conditions are used, and for the swirling jet flow from the nozzle, two types of downstream boundary conditions are used. The solutions are time accurate which are obtained using an implicit, upwind, flux-difference splitting, finite-volume scheme.

Boundary conditions for plasma fluid models at the magnetic presheath entrance

International Nuclear Information System (INIS)

Loizu, J.; Ricci, P.; Halpern, F. D.; Jolliet, S.

2012-01-01

The proper boundary conditions at the magnetic presheath entrance for plasma fluid turbulence models based on the drift approximation are derived, focusing on a weakly collisional plasma sheath with T i ≪T e and a magnetic field oblique to a totally absorbing wall. First, the location of the magnetic presheath entrance is rigorously derived. Then boundary conditions at the magnetic presheath entrance are analytically deduced for v ||i , v ||e , n, φ, T e , and for the vorticity ω=∇ ⊥ 2 φ. The effects of E × B and diamagnetic drifts on the boundary conditions are also investigated. Kinetic simulations are performed that confirm the analytical results. Finally, the new set of boundary conditions is implemented in a three-dimensional global fluid code for the simulation of plasma turbulence and, as an example, the results of a tokamak scrape-off layer simulation are discussed. The framework presented can be generalized to obtain boundary conditions at the magnetic presheath entrance in more complex scenarios.

A device adaptive inflow boundary condition for Wigner equations of quantum transport

International Nuclear Information System (INIS)

Jiang, Haiyan; Lu, Tiao; Cai, Wei

2014-01-01

In this paper, an improved inflow boundary condition is proposed for Wigner equations in simulating a resonant tunneling diode (RTD), which takes into consideration the band structure of the device. The original Frensley inflow boundary condition prescribes the Wigner distribution function at the device boundary to be the semi-classical Fermi–Dirac distribution for free electrons in the device contacts without considering the effect of the quantum interaction inside the quantum device. The proposed device adaptive inflow boundary condition includes this effect by assigning the Wigner distribution to the value obtained from the Wigner transform of wave functions inside the device at zero external bias voltage, thus including the dominant effect on the electron distribution in the contacts due to the device internal band energy profile. Numerical results on computing the electron density inside the RTD under various incident waves and non-zero bias conditions show much improvement by the new boundary condition over the traditional Frensley inflow boundary condition

On the symmetry of the boundary conditions of the volume potential

Science.gov (United States)

Kal'menov, Tynysbek Sh.; Arepova, Gaukhar; Suragan, Durvudkhan

2017-09-01

It is well known that the volume potential determines the mass or the charge distributed over the domain with density f. The volume potential is extensively used in function theory and embedding theorems. It is also well known that the volume potential gives a solution to an inhomogeneous equation. And it generates a linear self-adjoint operator. It is known that self-adjoint differential operators are generated by boundary conditions. In our previous papers for an arbitrary domain a boundary condition on the volume potential is given. In the past, it was not possible to prove the self-adjointness of these obtained boundary conditions. In the present paper, we prove the symmetry of boundary condition for the volume potential.

Neutron transport assembly calculation with non-zero net current boundary condition

International Nuclear Information System (INIS)

Jo, Chang Keun

1993-02-01

Fuel assembly calculation for the homogenized group constants is one of the most important parts in the reactor core analysis. The homogenized group constants of one a quarter assembly are usually generated for the nodal calculation of the reactor core. In the current nodal calculation, one or a quarter of the fuel assembly corresponds to a unit node. The homogenized group constant calculation for a fuel assembly proceeds through cell spectrum calculations, group condensation and cell homogenization calculations, two dimensional fuel assembly calculation, and then depletion calculations of fuel rods. To obtain the assembly wise homogenized group constants, the two dimensional transport calculation is usually performed. Most codes for the assembly wise homogenized group constants employ a zero net current boundary condition. CASMO-3 is such a code that is in wide use. The zero net current boundary condition is plausible and valid in an infinite reactor composed of the same kind of assemblies. However, the reactor is finite and the core is constructed by different kinds of assemblies. Hence, the assumption of the zero net current boundary condition is not valid in the actual reactor. The objective of this study is to develop a homogenization methodology that can treat any actual boundary condition, i.e. non-zero net current boundary condition. In order to treat the non-zero net current boundary condition, we modify CASMO-3. For the two-dimensional treatment in CASMO-3, a multigroup integral transport routine based on the method of transmission probability is used. The code performs assembly calculation with zero net current boundary condition. CASMO-3 is modified to consider the inhomogeneous source at the assembly boundary surface due to the non-zero net current. The modified version of CASMO-3 is called CASMO-3M. CASMO-3M is applied to several benchmark problems. In order to obtain the inhomogeneous source, the global calculation is performed. The local calculation

The analytical solution for drug delivery system with nonhomogeneous moving boundary condition

Science.gov (United States)

Saudi, Muhamad Hakimi; Mahali, Shalela Mohd; Harun, Fatimah Noor

2017-08-01

This paper discusses the development and the analytical solution of a mathematical model based on drug release system from a swelling delivery device. The mathematical model is represented by a one-dimensional advection-diffusion equation with nonhomogeneous moving boundary condition. The solution procedures consist of three major steps. Firstly, the application of steady state solution method, which is used to transform the nonhomogeneous moving boundary condition to homogeneous boundary condition. Secondly, the application of the Landau transformation technique that gives a significant impact in removing the advection term in the system of equation and transforming the moving boundary condition to a fixed boundary condition. Thirdly, the used of separation of variables method to find the analytical solution for the resulted initial boundary value problem. The results show that the swelling rate of delivery device and drug release rate is influenced by value of growth factor r.

Classically integrable boundary conditions for affine Toda field theories

International Nuclear Information System (INIS)

Bowcock, P.; Corrigan, E.; Dorey, P.E.; Rietdijk, R.H.

1995-01-01

Boundary conditions compatible with classical integrability are studied both directly, using an approach based on the explicit construction of conserved quantities, and indirectly by first developing a generalisation of the Lax pair idea. The latter approach is closer to the spirit of earlier work by Sklyanin and yields a complete set of conjectures for permissible boundary conditions for any affine Toda field theory. (orig.)

Facilitating conditions for boundary-spanning behavior in governance networks

OpenAIRE

Meerkerk, Ingmar; Edelenbos, Jurian

2017-01-01

textabstractThis article examines the impact of two facilitating conditions for boundary-spanning behaviour in urban governance networks. While research on boundary spanning is growing, there is little attention for antecedents. Combining governance network literature on project management and organizational literature on facilitative and servant leadership, we examine two potential conditions: a facilitative project management style and executive support. We conducted survey research among p...

A simple and efficient outflow boundary condition for the incompressible Navier–Stokes equations

Directory of Open Access Journals (Sweden)

Yibao Li

2017-01-01

Full Text Available Many researchers have proposed special treatments for outlet boundary conditions owing to lack of information at the outlet. Among them, the simplest method requires a large enough computational domain to prevent or reduce numerical errors at the boundaries. However, an efficient method generally requires special treatment to overcome the problems raised by the outlet boundary condition used. For example, mass flux is not conserved and the fluid field is not divergence-free at the outlet boundary. Overcoming these problems requires additional computational cost. In this paper, we present a simple and efficient outflow boundary condition for the incompressible Navier–Stokes equations, aiming to reduce the computational domain for simulating flow inside a long channel in the streamwise direction. The proposed outflow boundary condition is based on the transparent equation, where a weak formulation is used. The pressure boundary condition is derived by using the Navier–Stokes equations and the outlet flow boundary condition. In the numerical algorithm, a staggered marker-and-cell grid is used and temporal discretization is based on a projection method. The intermediate velocity boundary condition is consistently adopted to handle the velocity–pressure coupling. Characteristic numerical experiments are presented to demonstrate the robustness and accuracy of the proposed numerical scheme. Furthermore, the agreement of computational results from small and large domains suggests that our proposed outflow boundary condition can significantly reduce computational domain sizes.

Propagation of Quasi-plane Nonlinear Waves in Tubes

Directory of Open Access Journals (Sweden)

P. Koníček

2002-01-01

Full Text Available This paper deals with possibilities of using the generalized Burgers equation and the KZK equation to describe nonlinear waves in circular ducts. A new method for calculating of diffraction effects taking into account boundary layer effects is described. The results of numerical solutions of the model equations are compared. Finally, the limits of validity of the used model equations are discussed with respect to boundary conditions and the radius of the circular duct. The limits of applicability of the KZK equation and the GBE equation for describing nonlinear waves in tubes are discussed.

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

Science.gov (United States)

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

2013-04-01

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

Spectral distribution of scalar particles created by a moving boundary with Robin boundary condition

International Nuclear Information System (INIS)

Mintz, B.; Farina, C; Maia Neto, P.A.; Rodrigues, R.B.

2006-01-01

We consider a massless scalar field in 1+1 dimensions satisfying a Robin boundary condition (BC) at a non-relativistic boundary. By deriving a Bogoliubov transformation between the input and output bosonic field operators, we calculate the spectral distribution of created particles. The particular cases of Dirichlet and Neumann BC may be obtained from our result as limiting cases, yielding equal spectra (this result is valid only in this space-time dimensionality). The creation effect for the field under Dirichlet BC turns out to be an upper bound for the spectra derived for Robin BC. Also, we show that the particle creation phenomenon with Robin conditions can be considerably reduced (with respect to the Dirichlet or Neumann cases) by selecting a particular mechanical oscillation frequency of the moving boundary. (author)

Scattering through a straight quantum waveguide with combined boundary conditions

Czech Academy of Sciences Publication Activity Database

Briet, Ph.; Dittrich, Jaroslav; Soccorsi, E.

2014-01-01

Roč. 55, č. 11 (2014), s. 112104 ISSN 0022-2488 R&D Projects: GA ČR(CZ) GA14-06818S Institutional support: RVO:61389005 Keywords : quantum waveguide * scattering * Dirichlet boundary condition * Neumann boundary condition Subject RIV: BE - Theoretical Physics Impact factor: 1.243, year: 2014

'Duality twisted'boundary conditions in n-state Potts Models

International Nuclear Information System (INIS)

Schuetz, G.

1992-11-01

We discuss a new class of toroidal boundary conditions for one-dimensional quantum Hamiltonian with S n symmetry which are related to two-dimensional n-state Potts models in the extreme anisotropic Hamiltonian limit. At their self-dual point (a point were a second-order phase transition occurs for n=2,3,4) the duality transformation is shown to be an additional symmetry giving rise to a new class of 'duality twisted' toroidal boundary conditions. This corresponding Hamiltonians are given in terms of generators of the periodic Temprely-Lieb algebra with an odd number of generators. We discuss as an example the critical Ising model. Here the complete spectrum for the new boundary conditions can be obtained from a projection mechanism in the spin-1/2 XXZ Heisenberg chain. (author)

Differential quadrature method of nonlinear bending of functionally graded beam

Science.gov (United States)

Gangnian, Xu; Liansheng, Ma; Wang, Youzhi; Quan, Yuan; Weijie, You

2018-02-01

Using the third-order shear deflection beam theory (TBT), nonlinear bending of functionally graded (FG) beams composed with various amounts of ceramic and metal is analyzed utilizing the differential quadrature method (DQM). The properties of beam material are supposed to accord with the power law index along to thickness. First, according to the principle of stationary potential energy, the partial differential control formulae of the FG beams subjected to a distributed lateral force are derived. To obtain numerical results of the nonlinear bending, non-dimensional boundary conditions and control formulae are dispersed by applying the DQM. To verify the present solution, several examples are analyzed for nonlinear bending of homogeneous beams with various edges. A minute parametric research is in progress about the effect of the law index, transverse shear deformation, distributed lateral force and boundary conditions.

Temperature field conduction solution by incomplete boundary condition

Energy Technology Data Exchange (ETDEWEB)

Novakovic, M; Petrasinovic, Lj; Djuric, M [Tehnoloski fakultet, Novi Sad (Yugoslavia); Perovic, N [Institut za Nuklearne Nauke Boris Kidric, Belgrade (Yugoslavia)

1977-01-01

The problem of determination of one part boundary conditions temperatures for Fourier partial differential equation when the other part of boundary condition and derivates (heat fluxes) are known is a practical interest as it enables one to determine and accessible temperature by measuring temperatures on other side, of the wall. Method developed and applied here consist of transforming the Fourier partial differential equation by time discretisation in sets of pairs of ordinary differential equations for temperature and heat flux. Such pair of differential equations of first order was solved by Runge-Kutta method. The integration proceeds along space interval simultaneosly for all time intervals. It is interesting to note that this procedure does not require the initial condition.

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

International Nuclear Information System (INIS)

Frittelli, Simonetta; Gomez, Roberto

2004-01-01

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

Downstream-Conditioned Maximum Entropy Method for Exit Boundary Conditions in the Lattice Boltzmann Method

Directory of Open Access Journals (Sweden)

Javier A. Dottori

2015-01-01

Full Text Available A method for modeling outflow boundary conditions in the lattice Boltzmann method (LBM based on the maximization of the local entropy is presented. The maximization procedure is constrained by macroscopic values and downstream components. The method is applied to fully developed boundary conditions of the Navier-Stokes equations in rectangular channels. Comparisons are made with other alternative methods. In addition, the new downstream-conditioned entropy is studied and it was found that there is a correlation with the velocity gradient during the flow development.

Experimental verification of free-space singular boundary conditions in an invisibility cloak

International Nuclear Information System (INIS)

Wu, Qiannan; Gao, Fei; Song, Zhengyong; Lin, Xiao; Zhang, Youming; Zhang, Baile; Chen, Huanyang

2016-01-01

A major issue in invisibility cloaking, which caused intense mathematical discussions in the past few years but still remains physically elusive, is the plausible singular boundary conditions associated with the singular metamaterials at the inner boundary of an invisibility cloak. The perfect cloaking phenomenon, as originally proposed by Pendry et al for electromagnetic waves, cannot be treated as physical before a realistic inner boundary of a cloak is demonstrated. Although a recent demonstration has been done in a waveguide environment, the exotic singular boundary conditions should apply to a general environment as in free space. Here we fabricate a metamaterial surface that exhibits the singular boundary conditions and demonstrate its performance in free space. Particularly, the phase information of waves reflected from this metamaterial surface is explicitly measured, confirming the singular responses of boundary conditions for an invisibility cloak. (paper)

Experimental verification of free-space singular boundary conditions in an invisibility cloak

Science.gov (United States)

Wu, Qiannan; Gao, Fei; Song, Zhengyong; Lin, Xiao; Zhang, Youming; Chen, Huanyang; Zhang, Baile

2016-04-01

A major issue in invisibility cloaking, which caused intense mathematical discussions in the past few years but still remains physically elusive, is the plausible singular boundary conditions associated with the singular metamaterials at the inner boundary of an invisibility cloak. The perfect cloaking phenomenon, as originally proposed by Pendry et al for electromagnetic waves, cannot be treated as physical before a realistic inner boundary of a cloak is demonstrated. Although a recent demonstration has been done in a waveguide environment, the exotic singular boundary conditions should apply to a general environment as in free space. Here we fabricate a metamaterial surface that exhibits the singular boundary conditions and demonstrate its performance in free space. Particularly, the phase information of waves reflected from this metamaterial surface is explicitly measured, confirming the singular responses of boundary conditions for an invisibility cloak.

Boundary conditions for the diffusion equation in radiative transfer

International Nuclear Information System (INIS)

Haskell, R.C.; Svaasand, L.O.; Tsay, T.; Feng, T.; McAdams, M.S.; Tromberg, B.J.

1994-01-01

Using the method of images, we examine the three boundary conditions commonly applied to the surface of a semi-infinite turbid medium. We find that the image-charge configurations of the partial-current and extrapolated-boundary conditions have the same dipole and quadrupole moments and that the two corresponding solutions to the diffusion equation are approximately equal. In the application of diffusion theory to frequency-domain photon-migration (FDPM) data, these two approaches yield values for the scattering and absorption coefficients that are equal to within 3%. Moreover, the two boundary conditions can be combined to yield a remarkably simple, accurate, and computationally fast method for extracting values for optical parameters from FDPM data. FDPM data were taken both at the surface and deep inside tissue phantoms, and the difference in data between the two geometries is striking. If one analyzes the surface data without accounting for the boundary, values deduced for the optical coefficients are in error by 50% or more. As expected, when aluminum foil was placed on the surface of a tissue phantom, phase and modulation data were closer to the results for an infinite-medium geometry. Raising the reflectivity of a tissue surface can, in principle, eliminate the effect of the boundary. However, we find that phase and modulation data are highly sensitive to the reflectivity in the range of 80--100%, and a minimum value of 98% is needed to mimic an infinite-medium geometry reliably. We conclude that noninvasive measurements of optically thick tissue require a rigorous treatment of the tissue boundary, and we suggest a unified partial-current--extrapolated boundary approach

Nonlinear to Linear Elastic Code Coupling in 2-D Axisymmetric Media.

Energy Technology Data Exchange (ETDEWEB)

Preston, Leiph [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

2017-08-01

Explosions within the earth nonlinearly deform the local media, but at typical seismological observation distances, the seismic waves can be considered linear. Although nonlinear algorithms can simulate explosions in the very near field well, these codes are computationally expensive and inaccurate at propagating these signals to great distances. A linearized wave propagation code, coupled to a nonlinear code, provides an efficient mechanism to both accurately simulate the explosion itself and to propagate these signals to distant receivers. To this end we have coupled Sandia's nonlinear simulation algorithm CTH to a linearized elastic wave propagation code for 2-D axisymmetric media (axiElasti) by passing information from the nonlinear to the linear code via time-varying boundary conditions. In this report, we first develop the 2-D axisymmetric elastic wave equations in cylindrical coordinates. Next we show how we design the time-varying boundary conditions passing information from CTH to axiElasti, and finally we demonstrate the coupling code via a simple study of the elastic radius.

On Stable Wall Boundary Conditions for the Hermite Discretization of the Linearised Boltzmann Equation

Science.gov (United States)

Sarna, Neeraj; Torrilhon, Manuel

2018-01-01

We define certain criteria, using the characteristic decomposition of the boundary conditions and energy estimates, which a set of stable boundary conditions for a linear initial boundary value problem, involving a symmetric hyperbolic system, must satisfy. We first use these stability criteria to show the instability of the Maxwell boundary conditions proposed by Grad (Commun Pure Appl Math 2(4):331-407, 1949). We then recognise a special block structure of the moment equations which arises due to the recursion relations and the orthogonality of the Hermite polynomials; the block structure will help us in formulating stable boundary conditions for an arbitrary order Hermite discretization of the Boltzmann equation. The formulation of stable boundary conditions relies upon an Onsager matrix which will be constructed such that the newly proposed boundary conditions stay close to the Maxwell boundary conditions at least in the lower order moments.

Finite element solution of quasistationary nonlinear magnetic field

International Nuclear Information System (INIS)

Zlamal, Milos

1982-01-01

The computation of quasistationary nonlinear two-dimensional magnetic field leads to the following problem. There is given a bounded domain OMEGA and an open nonempty set R included in OMEGA. We are looking for the magnetic vector potential u(x 1 , x 2 , t) which satisifies: 1) a certain nonlinear parabolic equation and an initial condition in R: 2) a nonlinear elliptic equation in S = OMEGA - R which is the stationary case of the above mentioned parabolic equation; 3) a boundary condition on delta OMEGA; 4) u as well as its conormal derivative are continuous accross the common boundary of R and S. This problem is formulated in two equivalent abstract ways. There is constructed an approximate solution completely discretized in space by a generalized Galerkin method (straight finite elements are a special case) and by backward A-stable differentiation methods in time. Existence and uniqueness of a weak solution is proved as well as a weak and strong convergence of the approximate solution to this solution. There are also derived error bounds for the solution of the two-dimensional nonlinear magnetic field equations under the assumption that the exact solution is sufficiently smooth

Device Applications of Nonlinear Dynamics

CERN Document Server

Baglio, Salvatore

2006-01-01

This edited book is devoted specifically to the applications of complex nonlinear dynamic phenomena to real systems and device applications. While in the past decades there has been significant progress in the theory of nonlinear phenomena under an assortment of system boundary conditions and preparations, there exist comparatively few devices that actually take this rich behavior into account. "Device Applications of Nonlinear Dynamics" applies and exploits this knowledge to make devices which operate more efficiently and cheaply, while affording the promise of much better performance. Given the current explosion of ideas in areas as diverse as molecular motors, nonlinear filtering theory, noise-enhanced propagation, stochastic resonance and networked systems, the time is right to integrate the progress of complex systems research into real devices.

Essential Boundary Conditions with Straight C1 Finite Elements in Curved Domains

International Nuclear Information System (INIS)

Ferraro, N.M.; Jardin, S.C.; Luo, X.

2010-01-01

The implementation of essential boundary conditions in C1 finite element analysis requires proper treatment of both the boundary conditions on second-order differentials of the solution and the curvature of the domain boundary. A method for the imposition of essential boundary conditions using straight elements (where the elements are not deformed to approximate a curved domain) is described. It is shown that pre-multiplication of the matrix equation by the local rotation matrix at each boundary node is not the optimal transformation. The uniquely optimal transformation is found, which does not take the form of a similarity transformation due to the non-orthogonality of the transformation to curved coordinates.

The boundary conditions for point transformed electromagnetic invisibility cloaks

International Nuclear Information System (INIS)

Weder, Ricardo

2008-01-01

In this paper we study point transformed electromagnetic invisibility cloaks in transformation media that are obtained by transformation from general anisotropic media. We assume that there are several point transformed electromagnetic cloaks located in different points in space. Our results apply in particular to the first-order invisibility cloaks introduced by Pendry et al and to the high-order invisibility cloaks introduced by Hendi et al and by Cai et al. We identify the appropriate cloaking boundary conditions that the solutions of Maxwell equations have to satisfy at the outside, ∂K + , and at the inside, ∂K - , of the boundary of the cloaked object K in the case where the permittivity and the permeability are bounded below and above in K. Namely, that the tangential components of the electric and the magnetic fields have to vanish at ∂K + -which is always true-and that the normal components of the curl of the electric and the magnetic fields have to vanish at ∂K - . These results are proven requiring that energy be conserved. In the case of one spherical cloak with a spherically stratified K and a radial current at ∂K we verify by an explicit calculation that our cloaking boundary conditions are satisfied and that cloaking of active devices holds, even if the current is at the boundary of the cloaked object. As we prove our results for media that are obtained by transformation from general anisotropic media, our results apply to the cloaking of objects with passive and active devices contained in general anisotropic media, in particular to objects with passive and active devices contained inside general crystals. Our results suggest a method to enhance cloaking in the approximate transformation media that are used in practice. Namely, to coat the boundary of the cloaked object (the inner boundary of the cloak) with a material that imposes the boundary conditions above. As these boundary conditions have to be satisfied for exact transformation

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

Canonical group quantization and boundary conditions

International Nuclear Information System (INIS)

Jung, Florian

2012-01-01

In the present thesis, we study quantization of classical systems with non-trivial phase spaces using the group-theoretical quantization technique proposed by Isham. Our main goal is a better understanding of global and topological aspects of quantum theory. In practice, the group-theoretical approach enables direct quantization of systems subject to constraints and boundary conditions in a natural and physically transparent manner -- cases for which the canonical quantization method of Dirac fails. First, we provide a clarification of the quantization formalism. In contrast to prior treatments, we introduce a sharp distinction between the two group structures that are involved and explain their physical meaning. The benefit is a consistent and conceptually much clearer construction of the Canonical Group. In particular, we shed light upon the 'pathological' case for which the Canonical Group must be defined via a central Lie algebra extension and emphasise the role of the central extension in general. In addition, we study direct quantization of a particle restricted to a half-line with 'hard wall' boundary condition. Despite the apparent simplicity of this example, we show that a naive quantization attempt based on the cotangent bundle over the half-line as classical phase space leads to an incomplete quantum theory; the reflection which is a characteristic aspect of the 'hard wall' is not reproduced. Instead, we propose a different phase space that realises the necessary boundary condition as a topological feature and demonstrate that quantization yields a suitable quantum theory for the half-line model. The insights gained in the present special case improve our understanding of the relation between classical and quantum theory and illustrate how contact interactions may be incorporated.

The determination of an unknown boundary condition in a fractional diffusion equation

KAUST Repository

Rundell, William

2013-07-01

In this article we consider an inverse boundary problem, in which the unknown boundary function ∂u/∂v = f(u) is to be determined from overposed data in a time-fractional diffusion equation. Based upon the free space fundamental solution, we derive a representation for the solution f as a nonlinear Volterra integral equation of second kind with a weakly singular kernel. Uniqueness and reconstructibility by iteration is an immediate result of a priori assumption on f and applying the fixed point theorem. Numerical examples are presented to illustrate the validity and effectiveness of the proposed method. © 2013 Copyright Taylor and Francis Group, LLC.

Periodic Boundary Conditions in the ALEGRA Finite Element Code

International Nuclear Information System (INIS)

Aidun, John B.; Robinson, Allen C.; Weatherby, Joe R.

1999-01-01

This document describes the implementation of periodic boundary conditions in the ALEGRA finite element code. ALEGRA is an arbitrary Lagrangian-Eulerian multi-physics code with both explicit and implicit numerical algorithms. The periodic boundary implementation requires a consistent set of boundary input sets which are used to describe virtual periodic regions. The implementation is noninvasive to the majority of the ALEGRA coding and is based on the distributed memory parallel framework in ALEGRA. The technique involves extending the ghost element concept for interprocessor boundary communications in ALEGRA to additionally support on- and off-processor periodic boundary communications. The user interface, algorithmic details and sample computations are given

Computation of Nonlinear Backscattering Using a High-Order Numerical Method

Science.gov (United States)

Fibich, G.; Ilan, B.; Tsynkov, S.

2001-01-01

The nonlinear Schrodinger equation (NLS) is the standard model for propagation of intense laser beams in Kerr media. The NLS is derived from the nonlinear Helmholtz equation (NLH) by employing the paraxial approximation and neglecting the backscattered waves. In this study we use a fourth-order finite-difference method supplemented by special two-way artificial boundary conditions (ABCs) to solve the NLH as a boundary value problem. Our numerical methodology allows for a direct comparison of the NLH and NLS models and for an accurate quantitative assessment of the backscattered signal.

Discrete Green's Theorem, Green's Functions and Stable Radiative FDTD Boundary Conditions

NARCIS (Netherlands)

Arnold, J.M.; Hon, de B.P.

2007-01-01

We propose a radiative boundary condition for the discrete-grid formulation of Helmholtz’ equation, based on rational approximation in the frequency domain of a Green’s function for the discretised system. This boundary condition is free from instabilities.

Nonlinear electrorheological instability of two Rivlin-Ericksen elastico-viscous fluids

CERN Document Server

El-Dib, Y O

2003-01-01

The behaviour of surface waves propagating between two Rivlin-Ericksen elastico-viscous fluids is examined. The investigation is made in the presence of a vertical electric field and a relative horizontal constant velocity. The influence of both surface tension and gravity force is taken into account. Due to the inclusion of streaming flow a mathematical simplification is considered. The viscoelastic contribution is demonstrated in the boundary conditions. From this point of view the approximation equations of motion are solved in the absence of viscoelastic effects. The solutions of the linearized equations of motion under nonlinear boundary conditions lead to derivation of a nonlinear equation governing the interfacial displacement and having damping terms with complex coefficients. This equation is accomplished by utilizing the cubic nonlinearity. The use of the Gardner-Morikawa transformation yields a simplified linear dispersion relation so that the periodic solution for the linear form is utilized. The ...

Repulsive Casimir force from fractional Neumann boundary conditions

International Nuclear Information System (INIS)

Lim, S.C.; Teo, L.P.

2009-01-01

This Letter studies the finite temperature Casimir force acting on a rectangular piston associated with a massless fractional Klein-Gordon field at finite temperature. Dirichlet boundary conditions are imposed on the walls of a d-dimensional rectangular cavity, and a fractional Neumann condition is imposed on the piston that moves freely inside the cavity. The fractional Neumann condition gives an interpolation between the Dirichlet and Neumann conditions, where the Casimir force is known to be always attractive and always repulsive respectively. For the fractional Neumann boundary condition, the attractive or repulsive nature of the Casimir force is governed by the fractional order which takes values from zero (Dirichlet) to one (Neumann). When the fractional order is larger than 1/2, the Casimir force is always repulsive. For some fractional orders that are less than but close to 1/2, it is shown that the Casimir force can be either attractive or repulsive depending on the aspect ratio of the cavity and the temperature.

A Modified Conjugacy Condition and Related Nonlinear Conjugate Gradient Method

Directory of Open Access Journals (Sweden)

Shengwei Yao

2014-01-01

Full Text Available The conjugate gradient (CG method has played a special role in solving large-scale nonlinear optimization problems due to the simplicity of their very low memory requirements. In this paper, we propose a new conjugacy condition which is similar to Dai-Liao (2001. Based on this condition, the related nonlinear conjugate gradient method is given. With some mild conditions, the given method is globally convergent under the strong Wolfe-Powell line search for general functions. The numerical experiments show that the proposed method is very robust and efficient.

Nonlinear traveling waves in rotating Rayleigh-Bacute enard convection: Stability boundaries and phase diffusion

International Nuclear Information System (INIS)

Liu, Y.; Ecke, R.E.

1999-01-01

We present experimental measurements of a sidewall traveling wave in rotating Rayleigh-Bacute enard convection. The fluid, water with Prandtl number about 6.3, was confined in a 1-cm-high cylindrical cell with radius-to-height ratio Γ=5. We used simultaneous optical-shadowgraph, heat-transport, and local temperature measurements to determine the stability and characteristics of the traveling-wave state for dimensionless rotation rates 60<Ω<420. The state is well described by the one-dimensional complex Ginzburg-Landau (CGL) equation for which the linear and nonlinear coefficients were determined for Ω=274. The Eckhaus-Benjamin-Feir-stability boundary was established and the phase-diffusion coefficient and nonlinear group velocity were determined in the stable regime. Higher-order corrections to the CGL equation were also investigated. copyright 1999 The American Physical Society

Non-linear processes in the Earth atmosphere boundary layer

Science.gov (United States)

Grunskaya, Lubov; Valery, Isakevich; Dmitry, Rubay

2013-04-01

The work is connected with studying electromagnetic fields in the resonator Earth-Ionosphere. There is studied the interconnection of tide processes of geophysical and astrophysical origin with the Earth electromagnetic fields. On account of non-linear property of the resonator Earth-Ionosphere the tides (moon and astrophysical tides) in the electromagnetic Earth fields are kinds of polyharmonic nature. It is impossible to detect such non-linear processes with the help of the classical spectral analysis. Therefore to extract tide processes in the electromagnetic fields, the method of covariance matrix eigen vectors is used. Experimental investigations of electromagnetic fields in the atmosphere boundary layer are done at the distance spaced stations, situated on Vladimir State University test ground, at Main Geophysical Observatory (St. Petersburg), on Kamchatka pen., on Lake Baikal. In 2012 there was continued to operate the multichannel synchronic monitoring system of electrical and geomagnetic fields at the spaced apart stations: VSU physical experimental proving ground; the station of the Institute of Solar and Terrestrial Physics of Russian Academy of Science (RAS) at Lake Baikal; the station of the Institute of volcanology and seismology of RAS in Paratunka; the station in Obninsk on the base of the scientific and production society "Typhoon". Such investigations turned out to be possible after developing the method of scanning experimental signal of electromagnetic field into non- correlated components. There was used a method of the analysis of the eigen vectors ofthe time series covariance matrix for exposing influence of the moon tides on Ez. The method allows to distribute an experimental signal into non-correlated periodicities. The present method is effective just in the situation when energetical deposit because of possible influence of moon tides upon the electromagnetic fields is little. There have been developed and realized in program components

Linear and nonlinear dynamic analysis by boundary element method. Ph.D. Thesis, 1986 Final Report

Science.gov (United States)

Ahmad, Shahid

1991-01-01

An advanced implementation of the direct boundary element method (BEM) applicable to free-vibration, periodic (steady-state) vibration and linear and nonlinear transient dynamic problems involving two and three-dimensional isotropic solids of arbitrary shape is presented. Interior, exterior, and half-space problems can all be solved by the present formulation. For the free-vibration analysis, a new real variable BEM formulation is presented which solves the free-vibration problem in the form of algebraic equations (formed from the static kernels) and needs only surface discretization. In the area of time-domain transient analysis, the BEM is well suited because it gives an implicit formulation. Although the integral formulations are elegant, because of the complexity of the formulation it has never been implemented in exact form. In the present work, linear and nonlinear time domain transient analysis for three-dimensional solids has been implemented in a general and complete manner. The formulation and implementation of the nonlinear, transient, dynamic analysis presented here is the first ever in the field of boundary element analysis. Almost all the existing formulation of BEM in dynamics use the constant variation of the variables in space and time which is very unrealistic for engineering problems and, in some cases, it leads to unacceptably inaccurate results. In the present work, linear and quadratic isoparametric boundary elements are used for discretization of geometry and functional variations in space. In addition, higher order variations in time are used. These methods of analysis are applicable to piecewise-homogeneous materials, such that not only problems of the layered media and the soil-structure interaction can be analyzed but also a large problem can be solved by the usual sub-structuring technique. The analyses have been incorporated in a versatile, general-purpose computer program. Some numerical problems are solved and, through comparisons

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

Directory of Open Access Journals (Sweden)

Xiaofeng Zhang

2017-12-01

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

A novel investigation of a micropolar fluid characterized by nonlinear constitutive diffusion model in boundary layer flow and heat transfer

Science.gov (United States)

Sui, Jize; Zhao, Peng; Cheng, Zhengdong; Zheng, Liancun; Zhang, Xinxin

2017-02-01

The rheological and heat-conduction constitutive models of micropolar fluids (MFs), which are important non-Newtonian fluids, have been, until now, characterized by simple linear expressions, and as a consequence, the non-Newtonian performance of such fluids could not be effectively captured. Here, we establish the novel nonlinear constitutive models of a micropolar fluid and apply them to boundary layer flow and heat transfer problems. The nonlinear power law function of angular velocity is represented in the new models by employing generalized "n-diffusion theory," which has successfully described the characteristics of non-Newtonian fluids, such as shear-thinning and shear-thickening fluids. These novel models may offer a new approach to the theoretical understanding of shear-thinning behavior and anomalous heat transfer caused by the collective micro-rotation effects in a MF with shear flow according to recent experiments. The nonlinear similarity equations with a power law form are derived and the approximate analytical solutions are obtained by the homotopy analysis method, which is in good agreement with the numerical solutions. The results indicate that non-Newtonian behaviors involving a MF depend substantially on the power exponent n and the modified material parameter K 0 introduced by us. Furthermore, the relations of the engineering interest parameters, including local boundary layer thickness, local skin friction, and Nusselt number are found to be fitted by a quadratic polynomial to n with high precision, which enables the extraction of the rapid predictions from a complex nonlinear boundary-layer transport system.

Stabilization of exact nonlinear Timoshenko beams in space by boundary feedback

Science.gov (United States)

Do, K. D.

2018-05-01

Boundary feedback controllers are designed to stabilize Timoshenko beams with large translational and rotational motions in space under external disturbances. The exact nonlinear partial differential equations governing motion of the beams are derived and used in the control design. The designed controllers guarantee globally practically asymptotically (and locally practically exponentially) stability of the beam motions at the reference state. The control design, well-posedness and stability analysis are based on various relationships between the earth-fixed and body-fixed coordinates, Sobolev embeddings, and a Lyapunov-type theorem developed to study well-posedness and stability for a class of evolution systems in Hilbert space. Simulation results are included to illustrate the effectiveness of the proposed control design.

The influence of the magnetic boundary conditions on the nature of astrophysical convection

International Nuclear Information System (INIS)

Lopez, J.M.; Murphy, J.O.

1983-01-01

The effects of employing two forms of the boundary conditions for the magnetic field disturbance, H, are demonstrated. The appropriate conditions on H for current-free boundaries can be written as DHt-aH=0. The second case uses the conditions DH=0 at the lower and upper boundaries

SEACAS Theory Manuals: Part 1. Problem Formulation in Nonlinear Solid Mechancis

Energy Technology Data Exchange (ETDEWEB)

Attaway, S.W.; Laursen, T.A.; Zadoks, R.I.

1998-08-01

This report gives an introduction to the basic concepts and principles involved in the formulation of nonlinear problems in solid mechanics. By way of motivation, the discussion begins with a survey of some of the important sources of nonlinearity in solid mechanics applications, using wherever possible simple one dimensional idealizations to demonstrate the physical concepts. This discussion is then generalized by presenting generic statements of initial/boundary value problems in solid mechanics, using linear elasticity as a template and encompassing such ideas as strong and weak forms of boundary value problems, boundary and initial conditions, and dynamic and quasistatic idealizations. The notational framework used for the linearized problem is then extended to account for finite deformation of possibly inelastic solids, providing the context for the descriptions of nonlinear continuum mechanics, constitutive modeling, and finite element technology given in three companion reports.

Canonical group quantization and boundary conditions

Energy Technology Data Exchange (ETDEWEB)

Jung, Florian

2012-07-16

In the present thesis, we study quantization of classical systems with non-trivial phase spaces using the group-theoretical quantization technique proposed by Isham. Our main goal is a better understanding of global and topological aspects of quantum theory. In practice, the group-theoretical approach enables direct quantization of systems subject to constraints and boundary conditions in a natural and physically transparent manner -- cases for which the canonical quantization method of Dirac fails. First, we provide a clarification of the quantization formalism. In contrast to prior treatments, we introduce a sharp distinction between the two group structures that are involved and explain their physical meaning. The benefit is a consistent and conceptually much clearer construction of the Canonical Group. In particular, we shed light upon the 'pathological' case for which the Canonical Group must be defined via a central Lie algebra extension and emphasise the role of the central extension in general. In addition, we study direct quantization of a particle restricted to a half-line with 'hard wall' boundary condition. Despite the apparent simplicity of this example, we show that a naive quantization attempt based on the cotangent bundle over the half-line as classical phase space leads to an incomplete quantum theory; the reflection which is a characteristic aspect of the 'hard wall' is not reproduced. Instead, we propose a different phase space that realises the necessary boundary condition as a topological feature and demonstrate that quantization yields a suitable quantum theory for the half-line model. The insights gained in the present special case improve our understanding of the relation between classical and quantum theory and illustrate how contact interactions may be incorporated.

Lie and Q-Conditional Symmetries of Reaction-Diffusion-Convection Equations with Exponential Nonlinearities and Their Application for Finding Exact Solutions

Directory of Open Access Journals (Sweden)

Roman Cherniha

2018-04-01

Full Text Available This review is devoted to search for Lie and Q-conditional (nonclassical symmetries and exact solutions of a class of reaction-diffusion-convection equations with exponential nonlinearities. A complete Lie symmetry classification of the class is derived via two different algorithms in order to show that the result depends essentially on the type of equivalence transformations used for the classification. Moreover, a complete description of Q-conditional symmetries for PDEs from the class in question is also presented. It is shown that all the well-known results for reaction-diffusion equations with exponential nonlinearities follow as particular cases from the results derived for this class of reaction-diffusion-convection equations. The symmetries obtained for constructing exact solutions of the relevant equations are successfully applied. The exact solutions are compared with those found by means of different techniques. Finally, an application of the exact solutions for solving boundary-value problems arising in population dynamics is presented.

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

Science.gov (United States)

Colomo, F.; Pronko, A. G.

2005-05-01

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

An effective absorbing layer for the boundary condition in acoustic seismic wave simulation

Science.gov (United States)

Yao, Gang; da Silva, Nuno V.; Wu, Di

2018-04-01

Efficient numerical simulation of seismic wavefields generally involves truncating the Earth model in order to keep computing time and memory requirements down. Absorbing boundary conditions, therefore, are applied to remove the boundary reflections caused by this truncation, thereby allowing for accurate modeling of wavefields. In this paper, we derive an effective absorbing boundary condition for both acoustic and elastic wave simulation, through the simplification of the damping term of the split perfectly matched layer (SPML) boundary condition. This new boundary condition is accurate, cost-effective, and easily implemented, especially for high-performance computing. Stability analysis shows that this boundary condition is effectively as stable as normal (non-absorbing) wave equations for explicit time-stepping finite differences. We found that for full-waveform inversion (FWI), the strengths of the effective absorbing layer—a reduction of the computational and memory cost coupled with a simplistic implementation—significantly outweighs the limitation of incomplete absorption of outgoing waves relative to the SPML. More importantly, we demonstrate that this limitation can easily be overcome through the use of two strategies in FWI, namely variable cell size and model extension thereby fully compensating for the imperfectness of the proposed absorbing boundary condition.

Magnetohydrodynamic three-dimensional flow of viscoelastic nanofluid in the presence of nonlinear thermal radiation

Energy Technology Data Exchange (ETDEWEB)

Hayat, T. [Department of Mathematics, Quaid-I-Azam University 45320, Islamabad 44000 (Pakistan); Nonlinear Analysis and Applied Mathematics (NAAM) Research Group, Faculty of Science, King Abdulaziz University, P. O. Box 80203, Jeddah 21589 (Saudi Arabia); Muhammad, Taseer, E-mail: taseer_qau@yahoo.com [Department of Mathematics, Quaid-I-Azam University 45320, Islamabad 44000 (Pakistan); Alsaedi, A.; Alhuthali, M.S. [Nonlinear Analysis and Applied Mathematics (NAAM) Research Group, Faculty of Science, King Abdulaziz University, P. O. Box 80203, Jeddah 21589 (Saudi Arabia)

2015-07-01

Magnetohydrodynamic (MHD) three-dimensional flow of couple stress nanofluid in the presence of thermophoresis and Brownian motion effects is analyzed. Energy equation subject to nonlinear thermal radiation is taken into account. The flow is generated by a bidirectional stretching surface. Fluid is electrically conducting in the presence of a constant applied magnetic field. The induced magnetic field is neglected for a small magnetic Reynolds number. Mathematical formulation is performed using boundary layer analysis. Newly proposed boundary condition requiring zero nanoparticle mass flux is employed. The governing nonlinear mathematical problems are first converted into dimensionless expressions and then solved for the series solutions of velocities, temperature and nanoparticles concentration. Convergence of the constructed solutions is verified. Effects of emerging parameters on the temperature and nanoparticles concentration are plotted and discussed. Skin friction coefficients and Nusselt number are also computed and analyzed. It is found that the thermal boundary layer thickness is an increasing function of radiative effect. - Highlights: • Three-dimensional boundary layer flow of viscoelastic nanofluid is examined. • Nonlinear thermal radiation is analyzed. • Brownian motion and thermophoresis effects are present. • Recently developed condition requiring zero nanoparticle mass flux is implemented. • Construction of convergent solutions of nonlinear flow is possible.

Nonlinear soil-structure interaction analysis based on the boundary-element method in time domain with application to embedded foundation

International Nuclear Information System (INIS)

Wolf, J.P.; Darbre, G.R.

1985-01-01

The computational procedure of the so-called truncated indirect boundary-element method is derived. The latter, which is non-local in space and time, represents a rigorous generally applicable procedure for taking into account a layered halfspace in a non-linear soil-structure interaction analysis. As an example, the non-linear soil-structure interaction analysis of a structure embedded in a halfspace with partial uplift of the basement and separation of the side wall is investigated. (orig.)

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.

Single particle nonlocality, geometric phases and time-dependent boundary conditions

Science.gov (United States)

Matzkin, A.

2018-03-01

We investigate the issue of single particle nonlocality in a quantum system subjected to time-dependent boundary conditions. We discuss earlier claims according to which the quantum state of a particle remaining localized at the center of an infinite well with moving walls would be specifically modified by the change in boundary conditions due to the wall’s motion. We first prove that the evolution of an initially localized Gaussian state is not affected nonlocally by a linearly moving wall: as long as the quantum state has negligible amplitude near the wall, the boundary motion has no effect. This result is further extended to related confined time-dependent oscillators in which the boundary’s motion is known to give rise to geometric phases: for a Gaussian state remaining localized far from the boundaries, the effect of the geometric phases is washed out and the particle dynamics shows no traces of a nonlocal influence that would be induced by the moving boundaries.

Preconditioned characteristic boundary conditions based on artificial compressibility method for solution of incompressible flows

Science.gov (United States)

Hejranfar, Kazem; Parseh, Kaveh

2017-09-01

The preconditioned characteristic boundary conditions based on the artificial compressibility (AC) method are implemented at artificial boundaries for the solution of two- and three-dimensional incompressible viscous flows in the generalized curvilinear coordinates. The compatibility equations and the corresponding characteristic variables (or the Riemann invariants) are mathematically derived and then applied as suitable boundary conditions in a high-order accurate incompressible flow solver. The spatial discretization of the resulting system of equations is carried out by the fourth-order compact finite-difference (FD) scheme. In the preconditioning applied here, the value of AC parameter in the flow field and also at the far-field boundary is automatically calculated based on the local flow conditions to enhance the robustness and performance of the solution algorithm. The code is fully parallelized using the Concurrency Runtime standard and Parallel Patterns Library (PPL) and its performance on a multi-core CPU is analyzed. The incompressible viscous flows around a 2-D circular cylinder, a 2-D NACA0012 airfoil and also a 3-D wavy cylinder are simulated and the accuracy and performance of the preconditioned characteristic boundary conditions applied at the far-field boundaries are evaluated in comparison to the simplified boundary conditions and the non-preconditioned characteristic boundary conditions. It is indicated that the preconditioned characteristic boundary conditions considerably improve the convergence rate of the solution of incompressible flows compared to the other boundary conditions and the computational costs are significantly decreased.

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 Wavelet Bicoherence-Based Quadratic Nonlinearity Feature for Translational Axis Condition Monitoring

Directory of Open Access Journals (Sweden)

Yong Li

2014-01-01

Full Text Available The translational axis is one of the most important subsystems in modern machine tools, as its degradation may result in the loss of the product qualification and lower the control precision. Condition-based maintenance (CBM has been considered as one of the advanced maintenance schemes to achieve effective, reliable and cost-effective operation of machine systems, however, current vibration-based maintenance schemes cannot be employed directly in the translational axis system, due to its complex structure and the inefficiency of commonly used condition monitoring features. In this paper, a wavelet bicoherence-based quadratic nonlinearity feature is proposed for translational axis condition monitoring by using the torque signature of the drive servomotor. Firstly, the quadratic nonlinearity of the servomotor torque signature is discussed, and then, a biphase randomization wavelet bicoherence is introduced for its quadratic nonlinear detection. On this basis, a quadratic nonlinearity feature is proposed for condition monitoring of the translational axis. The properties of the proposed quadratic nonlinearity feature are investigated by simulations. Subsequently, this feature is applied to the real-world servomotor torque data collected from the X-axis on a high precision vertical machining centre. All the results show that the performance of the proposed feature is much better than that of original condition monitoring features.

Second order bounce back boundary condition for the lattice Boltzmann fluid simulation

International Nuclear Information System (INIS)

Kim, In Chan

2000-01-01

A new bounce back boundary method of the second order in error is proposed for the lattice Boltzmann fluid simulation. This new method can be used for the arbitrarily irregular lattice geometry of a non-slip boundary. The traditional bounce back boundary condition for the lattice Boltzmann simulation is of the first order in error. Since the lattice Boltzmann method is the second order scheme by itself, a boundary technique of the second order has been desired to replace the first order bounce back method. This study shows that, contrary to the common belief that the bounce back boundary condition is unilaterally of the first order, the second order bounce back boundary condition can be realized. This study also shows that there exists a generalized bounce back technique that can be characterized by a single interpolation parameter. The second order bounce back method can be obtained by proper selection of this parameter in accordance with the detailed lattice geometry of the boundary. For an illustrative purpose, the transient Couette and the plane Poiseuille flows are solved by the lattice Boltzmann simulation with various boundary conditions. The results show that the generalized bounce back method yields the second order behavior in the error of the solution, provided that the interpolation parameter is properly selected. Coupled with its intuitive nature and the ease of implementation, the bounce back method can be as good as any second order boundary method

Vibration Analysis of Annular Sector Plates under Different Boundary Conditions

Directory of Open Access Journals (Sweden)

Dongyan Shi

2014-01-01

Full Text Available An analytical framework is developed for the vibration analysis of annular sector plates with general elastic restraints along each edge of plates. Regardless of boundary conditions, the displacement solution is invariably expressed as a new form of trigonometric expansion with accelerated convergence. The expansion coefficients are treated as the generalized coordinates and determined using the Rayleigh-Ritz technique. This work allows a capability of modeling annular sector plates under a variety of boundary conditions and changing the boundary conditions as easily as modifying the material properties or dimensions of the plates. Of equal importance, the proposed approach is universally applicable to annular sector plates of any inclusion angles up to 2π. The reliability and accuracy of the current method are adequately validated through numerical examples.

Self-Similar Nanocavity Design with Ultrasmall Mode Volume for Single-Photon Nonlinearities

DEFF Research Database (Denmark)

Choi, Hyeongrak; Heuck, Mikkel; Englund, Dirk R.

2017-01-01

We propose a photonic crystal nanocavity design with self-similar electromagnetic boundary conditions, achieving ultrasmall mode volume (V-eff). The electric energy density of a cavity mode can be maximized in the air or dielectric region, depending on the choice of boundary conditions. We illust...... at the single-photon level. These features open new directions in cavity quantum electrodynamics, spectroscopy, and quantum nonlinear optics....

External boundary effects on simultaneous diffusion and reaction processes

International Nuclear Information System (INIS)

Le Roux, M.N.; Wilhelmsson, H.

1989-01-01

External boundaries influence the spatial and temporal structure of evolution of dynamic systems governed by reaction-diffusion equations. Critical limits, i.e. thresholds for explosive growth or onset of diffusion dominated decay, are found to be caused by the presence of the boundary and to depend on: the position of the boundary, where the density is assumed to be zero at any instant of time: the mutual weights (coefficients) and powers of the nonlinear reaction and diffusion processes; and the initial spatial distribution. However, for particular relations between the nonlinear powers of the reaction and diffusion terms the critical limits do not depend on the initial conditions. The results are obtained by simulation experiment for one, two and three dimensions. Trends in the dynamic evolution of the system with an external boundary imposed are compared with the corresponding analytic results obtained for free boundary. Interesting applications are found in various areas, e.g. in the field of high temperature fusion plasma where the evolution of the temperature profile for the so-called H-mode (constant plasma density) is described

Decoupling in an expanding universe boundary RG-flow affects initial conditions for inflation

CERN Document Server

Schalm, K; Van der Schaar, J P; Schalm, Koenraad; Shiu, Gary; Schaar, Jan Pieter van der

2004-01-01

We study decoupling in FRW spacetimes, emphasizing a Lagrangian description throughout. To account for the vacuum choice ambiguity in cosmological settings, we introduce an arbitrary boundary action representing the initial conditions. RG flow in these spacetimes naturally affects the boundary interactions. As a consequence the boundary conditions are sensitive to high-energy physics through irrelevant terms in the boundary action. Using scalar field theory as an example, we derive the leading dimension four irrelevant boundary operators. We discuss how the known vacuum choices, e.g. the Bunch-Davies vacuum, appear in the Lagrangian description and square with decoupling. For all choices of boundary conditions encoded by relevant boundary operators, of which the known ones are a subset, backreaction is under control. All, moreover, will generically feel the influence of high-energy physics through irrelevant (dimension four) boundary corrections. Having established a coherent effective field theory framework ...

Approximate source conditions for nonlinear ill-posed problems—chances and limitations

International Nuclear Information System (INIS)

Hein, Torsten; Hofmann, Bernd

2009-01-01

In the recent past the authors, with collaborators, have published convergence rate results for regularized solutions of linear ill-posed operator equations by avoiding the usual assumption that the solutions satisfy prescribed source conditions. Instead the degree of violation of such source conditions is expressed by distance functions d(R) depending on a radius R ≥ 0 which is an upper bound of the norm of source elements under consideration. If d(R) tends to zero as R → ∞ an appropriate balancing of occurring regularization error terms yields convergence rates results. This approach was called the method of approximate source conditions, originally developed in a Hilbert space setting. The goal of this paper is to formulate chances and limitations of an application of this method to nonlinear ill-posed problems in reflexive Banach spaces and to complement the field of low order convergence rates results in nonlinear regularization theory. In particular, we are going to establish convergence rates for a variant of Tikhonov regularization. To keep structural nonlinearity conditions simple, we update the concept of degree of nonlinearity in Hilbert spaces to a Bregman distance setting in Banach spaces

Galerkin methods for Boltzmann-Poisson transport with reflection conditions on rough boundaries

Science.gov (United States)

Morales Escalante, José A.; Gamba, Irene M.

2018-06-01

We consider in this paper the mathematical and numerical modeling of reflective boundary conditions (BC) associated to Boltzmann-Poisson systems, including diffusive reflection in addition to specularity, in the context of electron transport in semiconductor device modeling at nano scales, and their implementation in Discontinuous Galerkin (DG) schemes. We study these BC on the physical boundaries of the device and develop a numerical approximation to model an insulating boundary condition, or equivalently, a pointwise zero flux mathematical condition for the electron transport equation. Such condition balances the incident and reflective momentum flux at the microscopic level, pointwise at the boundary, in the case of a more general mixed reflection with momentum dependant specularity probability p (k →). We compare the computational prediction of physical observables given by the numerical implementation of these different reflection conditions in our DG scheme for BP models, and observe that the diffusive condition influences the kinetic moments over the whole domain in position space.

Nonlinear saturation of the Rayleigh Taylor instability

International Nuclear Information System (INIS)

Das, A.; Mahajan, S.; Kaw, P.; Sen, A.; Benkadda, S.; Verga, A.

1997-01-01

The problem of the nonlinear saturation of the 2 dimensional Rayleigh Taylor instability is re-examined to put various earlier results in a proper perspective. The existence of a variety of final states can be attributed to the differences in the choice of boundary conditions and initial conditions in earlier numerical modeling studies. Our own numerical simulations indicate that the RT instability saturates by the self consistent generation of shear flow even in situations (with periodic boundaries) where, in principle, an infinite amount of gravitational energy can be tapped. Such final states can be achieved for suitable values of the Prandtl number. (author)

δ'-function perturbations and Neumann boundary-conditions by path integration

International Nuclear Information System (INIS)

Grosche, C.

1994-02-01

δ'-function perturbations and Neumann boundary conditions are incorporated into the path integral formalism. The starting point is the consideration of the path integral representation for the one dimensional Dirac particle together with a relativistic point interaction. The non-relativistic limit yields either a usual δ-function or a δ'-function perturbation; making their strengths infinitely repulsive one obtains Dirichlet, respectively Neumann boundary conditions in the path integral. (orig.)

NLIE of Dirichlet sine-Gordon model for boundary bound states

International Nuclear Information System (INIS)

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

2008-01-01

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

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

Matrix albedo for discrete ordinates infinite-medium boundary condition

International Nuclear Information System (INIS)

Mathews, K.; Dishaw, J.

2007-01-01

Discrete ordinates problems with an infinite exterior medium (reflector) can be more efficiently computed by eliminating grid cells in the exterior medium and applying a matrix albedo boundary condition. The albedo matrix is a discretized bidirectional reflection distribution function (BRDF) that accounts for the angular quadrature set, spatial quadrature method, and spatial grid that would have been used to model a portion of the exterior medium. The method is exact in slab geometry, and could be used as an approximation in multiple dimensions or curvilinear coordinates. We present an adequate method for computing albedo matrices and demonstrate their use in verifying a discrete ordinates code in slab geometry by comparison with Ganapol's infinite medium semi-analytic TIEL benchmark. With sufficient resolution in the spatial and angular grids and iteration tolerance to yield solutions converged to 6 digits, the conventional (scalar) albedo boundary condition yielded 2-digit accuracy at the boundary, but the matrix albedo solution reproduced the benchmark scalar flux at the boundary to all 6 digits. (authors)

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.

The influence of boundary conditions on domain structure stability in spin wave approximation

International Nuclear Information System (INIS)

Wachinewski, A.

1974-01-01

Instead of the usually used Born-Karman cyclic conditions, boundary conditions which take into account the situation of the boundary lattice sites lying on the crystal's surface are assumed. It is shown that the particular choice of the boundary conditions secures the stability of domain structure in ferromagnet (positive spin wave energies), without including the Winter term in Hamiltonian. (author)

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

Directory of Open Access Journals (Sweden)

Mitsuhiro Nakao

2014-01-01

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

Magnetohydrodynamic boundary layer flow past a porous substrate with Beavers-Joseph boundary condition

International Nuclear Information System (INIS)

Jat, R.N.; Chaudhary, Santosh

2009-01-01

The flow of an electrically conducting fluid past a porous substrate attached to the flat plate with Beavers-Joseph boundary condition under the influence of a uniform transverse magnetic field has been studied. Taking suitable similar variables, the momentum equation is transformed to ordinary differential equation and solved by standard techniques. The energy equation is solved by considering two boundary layers, one in the porous substrate and the other above the porous substrate. The velocity and temperature distributions along with Nusselt number are discussed numerically and presented through graphs. (author)

Unconditionally stable methods for simulating multi-component two-phase interface models with Peng-Robinson equation of state and various boundary conditions

KAUST Repository

Kou, Jisheng

2015-03-01

In this paper, we consider multi-component dynamic two-phase interface models, which are formulated by the Cahn-Hilliard system with Peng-Robinson equation of state and various boundary conditions. These models can be derived from the minimum problems of Helmholtz free energy or grand potential in the realistic thermodynamic systems. The resulted Cahn-Hilliard systems with various boundary conditions are fully coupled and strongly nonlinear. A linear transformation is introduced to decouple the relations between different components, and as a result, the models are simplified. From this, we further propose a semi-implicit unconditionally stable time discretization scheme, which allows us to solve the Cahn-Hilliard system by a decoupled way, and thus, our method can significantly reduce the computational cost and memory requirements. The mixed finite element methods are employed for the spatial discretization, and the approximate errors are also analyzed for both space and time. Numerical examples are tested to demonstrate the efficiency of our proposed methods. © 2015 Elsevier B.V.

The exact solution to the one-dimensional Poisson–Boltzmann equation with asymmetric boundary conditions

DEFF Research Database (Denmark)

Johannessen, Kim

2014-01-01

The exact solution to the one-dimensional Poisson–Boltzmann equation with asymmetric boundary conditions can be expressed in terms of the Jacobi elliptic functions. The boundary conditions determine the modulus of the Jacobi elliptic functions. The boundary conditions can not be solved analytically...

Effects of Collisionality on the Nonlinear Characteristics of Boundary Turbulence and Blob/hole Transport in Tokamak Plasmas

Energy Technology Data Exchange (ETDEWEB)

Li, J.; Yasuaki, K., E-mail: lijq@energy.kyoto-u.ac.jp [Kyoto University, Kyoto (Japan); Cheng, J.; Longwen, Y.; Jiaqi, D. [Southwestern Institute of Physics, Chengdu (China)

2012-09-15

Full text: Blob/hole dynamics near tokamak separatrix is of striking importance in determining the boundary transport. Based on simulations using an extended 2-region (edge/SOL) fluid model, we found that blob/hole dynamics are sensitively influenced by the plasma collisionality, i.e., ion-electron and ion-neutral collisions. Namely, the holes are enhanced in highly collisional edge whereas the blobs are weakened at the SOL, causing larger particle convection. These blob/hole dynamics are closely correlated with potential dipoles. The trends are experimentally evidenced on the HL-2A tokamak. Moreover, as the neutral-ion collision increases, the blobs at the SOL tend to develop into streamers propagating outwards with reduced amplitude while the holes inwards are suppressed, showing a key role in nonlinear structure regulation and resultant transport suppression. Results suggest that adjusting the plasma collisionality by fueling, e.g., gas puffing, could serve as a method to nonlinearly select turbulent structures, i.e., blobs, holes or streamers, to access the control of boundary transport. (author)

Effects of boundary conditions on thermomechanical calculations: Spent fuel test - climax

International Nuclear Information System (INIS)

Butkovich, T.R.

1982-10-01

The effects of varying certain boundary conditions on the results of finite-element calculations were studied in relation to the Spent Fuel Test - Climax. The study employed a thermomechanical model with the ADINA structural analysis. Nodal temperature histories were generated with the compatible ADINAT heat flow codes. The boundary conditions studied included: (1) The effect of boundary loading on three progressively larger meshes. (2) Plane strain vs plane stress conditions. (3) The effect of isothermal boundaries on a small mesh and on a significantly larger mesh. The results showed that different mesh sizes had an insignificant effect on isothermal boundaries up to 5 y, while on the smallest and largest mesh, the maximum temperature difference in the mesh was 0 C. In the corresponding ADINA calculation, these different mesh sizes produce insignificant changes in the stress field and displacements in the region of interest near the heat sources and excavations. On the other hand, plane stress produces horizontal and vertical stress differences approx. 9% higher than does plane strain

Ducks in space: from nonlinear absolute instability to noise-sustained structures in a pattern-forming system

Science.gov (United States)

Avitabile, D.; Desroches, M.; Knobloch, E.; Krupa, M.

2017-11-01

A subcritical pattern-forming system with nonlinear advection in a bounded domain is recast as a slow-fast system in space and studied using a combination of geometric singular perturbation theory and numerical continuation. Two types of solutions describing the possible location of stationary fronts are identified, whose origin is traced to the onset of convective and absolute instability when the system is unbounded. The former are present only for non-zero upstream boundary conditions and provide a quantitative understanding of noise-sustained structures in systems of this type. The latter correspond to the onset of a global mode and are present even with zero upstream boundary conditions. The role of canard trajectories in the nonlinear transition between these states is clarified and the stability properties of the resulting spatial structures are determined. Front location in the convective regime is highly sensitive to the upstream boundary condition, and its dependence on this boundary condition is studied using a combination of numerical continuation and Monte Carlo simulations of the partial differential equation. Statistical properties of the system subjected to random or stochastic boundary conditions at the inlet are interpreted using the deterministic slow-fast spatial dynamical system.

A nonlinear plate control without linearization

Directory of Open Access Journals (Sweden)

Yildirim Kenan

2017-03-01

Full Text Available In this paper, an optimal vibration control problem for a nonlinear plate is considered. In order to obtain the optimal control function, wellposedness and controllability of the nonlinear system is investigated. The performance index functional of the system, to be minimized by minimum level of control, is chosen as the sum of the quadratic 10 functional of the displacement. The velocity of the plate and quadratic functional of the control function is added to the performance index functional as a penalty term. By using a maximum principle, the nonlinear control problem is transformed to solving a system of partial differential equations including state and adjoint variables linked by initial-boundary-terminal conditions. Hence, it is shown that optimal control of the nonlinear systems can be obtained without linearization of the nonlinear term and optimal control function can be obtained analytically for nonlinear systems without linearization.

Influence of different boundary conditions on analysis of SSI

International Nuclear Information System (INIS)

Wang Jiachun

2005-01-01

In the discussions of structural response to earthquakes, it has been assumed that the foundation medium is very stiff and that the seismic motions applied at the structure support points are the same as the free-field earthquake motions at those locations; in other words, the effects of soil structure interaction (SSI) have been neglected. However, its effects can be significant when the structure supported on a soft soil. Structures on the ground are affected by ground motion when there is seismic loading. The inability of the foundation to resist to deformation of soil would cause huge damages on the structures. The different codes and boundary conditions affect on analysis results of SSI. A comparison of the reactor buildings response as predicted by CLASSI and FLUSH shows substantial differences. To absorb, rather than reflect, the outwardly radiated energy, transmitting boundary conditions and soil structure interface should be taken into consideration in analysis of SSI. The paper discusses influence of several different boundary conditions on analysis of SSI. (author)

Boundary controllability for a nonlinear beam equation

Directory of Open Access Journals (Sweden)

Xiao-Min Cao

2015-09-01

Full Text Available This article concerns a nonlinear system modeling the bending vibrations of a nonlinear beam of length $L>0$. First, we derive the existence of long time solutions near an equilibrium. Then we prove that the nonlinear beam is locally exact controllable around the equilibrium in $H^4(0,L$ and with control functions in $H^2(0,T$. The approach we used are open mapping theorem, local controllability established by linearization, and the induction.

Effect of boundary conditions on thermohydraulic behavior of clay buffer used in nuclear waste repository

International Nuclear Information System (INIS)

Arul Peter, A.; Murugesan, K.; Mamidi, Ganesh; Sharma, Umesh Kumar; Sharma, D. Akanshu; Arora, Puneet

2010-01-01

The use of nuclear energy is increasing dramatically in the world due to the fast depletion of fossil fuels, and hence the nuclear waste disposal and its short and long-term effects are of considerable importance. One of the options considered for nuclear waste disposal is underground nuclear waste repository facility. In this underground nuclear waste disposal system the waste filled canisters are placed in the rock surrounded by an engineered clay barrier and the whole system is buried in the geological formation, which serves as the natural or geological barrier. The important characteristic of the clay barrier is that it should not open up for radiation though it is continuously subjected to heat loading from the canisters. The heat and moisture transport mechanisms through the clay barrier plays an important role in deciding its mechanical strength. Clay behaves as an unsaturated porous material when it is used as a buffer material in nuclear waste facility. The governing equations for heat and moisture transfer through unsaturated porous media are coupled and nonlinear and hence they have to be solved using numerical solution technique. This paper reports the results of a numerical study on heat and moisture transport through a buffer layer made of clay as used in nuclear waste repository. Galerkin's weighted residual finite element method has been employed for the solution of the non-linear coupled governing equations used to represent the heat and moisture transport through unsaturated clay material. A detailed computational procedure has been established for the solution of the non-linear governing equations using Newton-Raphson technique. Initially the code has been validated with available experimental results. Then numerical simulation results were obtained for heat and moisture variations within the buffer material for Dirichlet temperature boundary conditions in the range, 50 deg C 2 2 , with an aim to simulate the boundary conditions which the clay

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

International Nuclear Information System (INIS)

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

1992-11-01

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

Effects of Uncertainties in Electric Field Boundary Conditions for Ring Current Simulations

Science.gov (United States)

Chen, Margaret W.; O'Brien, T. Paul; Lemon, Colby L.; Guild, Timothy B.

2018-01-01

Physics-based simulation results can vary widely depending on the applied boundary conditions. As a first step toward assessing the effect of boundary conditions on ring current simulations, we analyze the uncertainty of cross-polar cap potentials (CPCP) on electric field boundary conditions applied to the Rice Convection Model-Equilibrium (RCM-E). The empirical Weimer model of CPCP is chosen as the reference model and Defense Meteorological Satellite Program CPCP measurements as the reference data. Using temporal correlations from a statistical analysis of the "errors" between the reference model and data, we construct a Monte Carlo CPCP discrete time series model that can be generalized to other model boundary conditions. RCM-E simulations using electric field boundary conditions from the reference model and from 20 randomly generated Monte Carlo discrete time series of CPCP are performed for two large storms. During the 10 August 2000 storm main phase, the proton density at 10 RE at midnight was observed to be low (Dst index is bounded by the simulated Dst values. In contrast, the simulated Dst values during the recovery phases of the 10 August 2000 and 31 August 2005 storms tend to underestimate systematically the observed late Dst recovery. This suggests a need to improve the accuracy of particle loss calculations in the RCM-E model. Application of this technique can aid modelers to make efficient choices on either investing more effort on improving specification of boundary conditions or on improving descriptions of physical processes.

Three-party quantum teleportation using thermal states in Heisenberg XX model with open boundary condition

International Nuclear Information System (INIS)

Bhan, Jaemi; Kwon, Younghun

2007-01-01

Recently Yeo showed that thermal states in Heisenberg XX model with periodic boundary condition could be used for three-party quantum teleportation. However it is hard to implement the periodic boundary condition in spin chain. So instead of imposing the periodic boundary condition, we consider open boundary condition in Heisenberg XX model and investigate the possibility of using thermal states in Heisenberg XX model with open boundary condition. Using this way, we find the best fidelity conditions to three known protocols in three-party quantum teleportation. It turns out that the best fidelity in every protocol would be 23

Unbounded Perturbations of Nonlinear Second-Order Difference Equations at Resonance

Directory of Open Access Journals (Sweden)

Ma Ruyun

2007-01-01

Full Text Available We study the existence of solutions of nonlinear discrete boundary value problems , , , where is the first eigenvalue of the linear problem , , , satisfies some “asymptotic nonuniform” resonance conditions, and for .

New boundary conditions for 3D RF modelling

International Nuclear Information System (INIS)

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

1990-01-01

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

Antishocks in the ASEP with open boundaries conditioned on low current

International Nuclear Information System (INIS)

Belitsky, V; Schütz, G M

2013-01-01

We study the time evolution of the ASEP on a finite lattice with L sites and open boundaries, conditioned on an atypically low current up to a finite time t. By an exact computation, we show that for a one-parameter family of boundary densities and a special value of the conditioned current, an initial product measure with an antishock at site k evolves into a convex combination of such antishocks at sites k′. The weights p(k′, t|k, 0) are shown to be the transition probabilities of simple biased random walk with reflecting boundaries. We compute explicitly these transition rates. Our result implies that the antishock remains microscopically stable under the locally conditioned dynamics. (paper)

''Free-space'' boundary conditions for the time-dependent wave equation

International Nuclear Information System (INIS)

Lindman, E.L.

1975-01-01

Boundary conditions for the discrete wave equation which act like an infinite region of free space in contact with the computational region can be constructed using projection operators. Propagating and evanescent waves coming from within the computational region generate no reflected waves as they cross the boundary. At the same time arbitrary waves may be launched into the computational region. Well known projection operators for one-dimensional waves may be used for this purpose in one dimension. Extensions of these operators to higher dimensions along with numerically efficient approximations to them are described for higher-dimensional problems. The separation of waves into ingoing and outgoing waves inherent in these boundary conditions greatly facilitates diagnostics

An Irrotational Flow Field That Approximates Flat Plate Boundary Conditions

OpenAIRE

Ruffa, Anthony A.

2004-01-01

An irrotational solution is derived for the steady-state Navier-Stokes equations that approximately satisfies the boundary conditions for flow over a finite flat plate. The nature of the flow differs substantially from boundary layer flow, with severe numerical difficulties in some regions.

Schrödinger functional boundary conditions and improvement for N > 3

DEFF Research Database (Denmark)

Hietanen, A.; Karavirta, T.; Vilaseca, P.

2014-01-01

The standard method to calculate non-perturbatively the evolution of the running coupling of a SU(N ) gauge theory is based on the Schrodinger functional (SF). In this paper we construct a family of boundary fields for general values of N which enter the standard definition of the SF coupling. We...... provide spatial boundary conditions for fermions in several representations which reduce the condition number of the squared Dirac operator. In addition, we calculate the improvement coefficients for N > 3 needed to remove boundary cutoff effects from the gauge action. After this, residual cutoff effects...

Entropy generation due to double diffusive convective flow of Casson fluids over nonlinearity stretching sheets with slip conditions

Directory of Open Access Journals (Sweden)

Sameh E. Ahmed

2017-12-01

Full Text Available The present paper deals with the effects of slip boundary conditions and chemical reaction on the heat and mass transfer by mixed convective boundary layer flow of a non-Newtonian fluid over a nonlinear stretching sheet. The Casson fluid model is used to characterize the non-Newtonian fluid behavior. First order chemical reactions are considered. Similar solutions are used to convert the partial differential equations governing the problem to ordinary differential equations. The velocity, temperature and concentration profiles are obtained, numerically, using the MATLAB function bvp4c and those are used to compute the entropy generation number. The effect of increasing values of the Casson parameter is found to suppress the velocity field and temperature distribution. But the concentration is enhanced with the increasing of Casson parameter. The viscous dissipation, temperature and concentration irreversibility are determined and discussed in details.

Solution to random differential equations with boundary conditions

Directory of Open Access Journals (Sweden)

Fairouz Tchier

2017-04-01

Full Text Available We study a family of random differential equations with boundary conditions. Using a random fixed point theorem, we prove an existence theorem that yields a unique random solution.

A highly precise frequency-based method for estimating the tension of an inclined cable with unknown boundary conditions

Science.gov (United States)

Ma, Lin

2017-11-01

This paper develops a method for precisely determining the tension of an inclined cable with unknown boundary conditions. First, the nonlinear motion equation of an inclined cable is derived, and a numerical model of the motion of the cable is proposed using the finite difference method. The proposed numerical model includes the sag-extensibility, flexural stiffness, inclination angle and rotational stiffness at two ends of the cable. Second, the influence of the dynamic parameters of the cable on its frequencies is discussed in detail, and a method for precisely determining the tension of an inclined cable is proposed based on the derivatives of the eigenvalues of the matrices. Finally, a multiparameter identification method is developed that can simultaneously identify multiple parameters, including the rotational stiffness at two ends. This scheme is applicable to inclined cables with varying sag, varying flexural stiffness and unknown boundary conditions. Numerical examples indicate that the method provides good precision. Because the parameters of cables other than tension (e.g., the flexural stiffness and rotational stiffness at the ends) are not accurately known in practical engineering, the multiparameter identification method could further improve the accuracy of cable tension measurements.

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)

Bosonization relations as bag boundary conditions

International Nuclear Information System (INIS)

Nadkarni, S.; Nielsen, H.B.; Zahed, I.

1984-10-01

The more sophisticated bag models of hadrons become, the less precisely they seem to determine the bag radius. Idealizing this situation leads to the concept of exact bag models - ''Cheshire Cat'' models, CCM'S - where the physics is completely insensitive to changes in the bag radius. CCM's are constructed explitly in 1+1-dimensions, where exact bosonization relations are known. In the formalism of bag models, these relations appear as boundary conditions which ensure that the shifting of the bag wall has no physical effect. Other notable features of 1+1-dimensional CCM's are: (i) Fermion number, though classically confined, can escape the bag via a vector current anomaly at the surface. (ii) Essentially the same boundary action works for a variety of models and its symmetries determine those of the external boson fields. Remarkably enough, this 1+1-dimensional boundary action has precisely the same form as the one used in 3+1-dimensional chiral bag models, lending support to the belief that the latter are indeed approximateCCM's. These 1+1-dimensional results are expected to provide useful guidelines in the attempt to, at least approximately, besonize 3+1-dimensional QCD. (orig.)

Tear film dynamics with evaporation, wetting, and time-dependent flux boundary condition on an eye-shaped domain

Science.gov (United States)

Li, Longfei; Braun, R. J.; Maki, K. L.; Henshaw, W. D.; King-Smith, P. E.

2014-01-01

We study tear film dynamics with evaporation on a wettable eye-shaped ocular surface using a lubrication model. The mathematical model has a time-dependent flux boundary condition that models the cycles of tear fluid supply and drainage; it mimics blinks on a stationary eye-shaped domain. We generate computational grids and solve the nonlinear governing equations using the OVERTURE computational framework. In vivo experimental results using fluorescent imaging are used to visualize the influx and redistribution of tears for an open eye. Results from the numerical simulations are compared with the experiment. The model captures the flow around the meniscus and other dynamic features of human tear film observed in vivo. PMID:24926191

Boundary conditions for open quantum systems driven far from equilibrium

Science.gov (United States)

Frensley, William R.

1990-07-01

This is a study of simple kinetic models of open systems, in the sense of systems that can exchange conserved particles with their environment. The system is assumed to be one dimensional and situated between two particle reservoirs. Such a system is readily driven far from equilibrium if the chemical potentials of the reservoirs differ appreciably. The openness of the system modifies the spatial boundary conditions on the single-particle Liouville-von Neumann equation, leading to a non-Hermitian Liouville operator. If the open-system boundary conditions are time reversible, exponentially growing (unphysical) solutions are introduced into the time dependence of the density matrix. This problem is avoided by applying time-irreversible boundary conditions to the Wigner distribution function. These boundary conditions model the external environment as ideal particle reservoirs with properties analogous to those of a blackbody. This time-irreversible model may be numerically evaluated in a discrete approximation and has been applied to the study of a resonant-tunneling semiconductor diode. The physical and mathematical properties of the irreversible kinetic model, in both its discrete and its continuum formulations, are examined in detail. The model demonstrates the distinction in kinetic theory between commutator superoperators, which may become non-Hermitian to describe irreversible behavior, and anticommutator superoperators, which remain Hermitian and are used to evaluate physical observables.

Threshold condition for nonlinear tearing modes in tokamaks

Energy Technology Data Exchange (ETDEWEB)

Zabiego, M.F. [Association Euratom-CEA, Centre d`Etudes de Cadarache, 13 - Saint-Paul-lez-Durance (France). Dept. de Recherches sur la Fusion Controlee; Callen, J.D. [Wisconsin Univ., Madison, WI (United States). Dept. of Nuclear Engineering and Engineering Physics

1996-04-01

Low-mode-number tearing mode nonlinear evolution is analyzed emphasizing the need for a threshold condition, to account for observations in tokamaks. The discussion is illustrated by two models recently introduced in the literature. Introducing a threshold condition in the tearing mode stability analysis is found to reveal some bifurcation points and thus domains of intrinsic stability in the island dynamics operational space. (author). 19 refs.

Threshold condition for nonlinear tearing modes in tokamaks

International Nuclear Information System (INIS)

Zabiego, M.F.; Callen, J.D.

1996-04-01

Low-mode-number tearing mode nonlinear evolution is analyzed emphasizing the need for a threshold condition, to account for observations in tokamaks. The discussion is illustrated by two models recently introduced in the literature. Introducing a threshold condition in the tearing mode stability analysis is found to reveal some bifurcation points and thus domains of intrinsic stability in the island dynamics operational space. (author)

Imperfection Sensitivity of Nonlinear Vibration of Curved Single-Walled Carbon Nanotubes Based on Nonlocal Timoshenko Beam Theory

Directory of Open Access Journals (Sweden)

Iman Eshraghi

2016-09-01

Full Text Available Imperfection sensitivity of large amplitude vibration of curved single-walled carbon nanotubes (SWCNTs is considered in this study. The SWCNT is modeled as a Timoshenko nano-beam and its curved shape is included as an initial geometric imperfection term in the displacement field. Geometric nonlinearities of von Kármán type and nonlocal elasticity theory of Eringen are employed to derive governing equations of motion. Spatial discretization of governing equations and associated boundary conditions is performed using differential quadrature (DQ method and the corresponding nonlinear eigenvalue problem is iteratively solved. Effects of amplitude and location of the geometric imperfection, and the nonlocal small-scale parameter on the nonlinear frequency for various boundary conditions are investigated. The results show that the geometric imperfection and non-locality play a significant role in the nonlinear vibration characteristics of curved SWCNTs.

Propagation of Quasi-plane Nonlinear Waves in Tubes

OpenAIRE

P. Koníček; M. Bednařík; M. Červenka

2002-01-01

This paper deals with possibilities of using the generalized Burgers equation and the KZK equation to describe nonlinear waves in circular ducts. A new method for calculating of diffraction effects taking into account boundary layer effects is described. The results of numerical solutions of the model equations are compared. Finally, the limits of validity of the used model equations are discussed with respect to boundary conditions and the radius of the circular duct. The limits of applicabi...

A methodology for on-line calculation of temperature and thermal stress under non-linear boundary conditions

International Nuclear Information System (INIS)

Botto, D.; Zucca, S.; Gola, M.M.

2003-01-01

In the literature many works have been written dealing with the task of on-line calculation of temperature and thermal stress for machine components and structures, in order to evaluate fatigue damage accumulation and estimate residual life. One of the most widespread methodologies is the Green's function technique (GFT), by which machine parameters such as fluid temperatures, pressures and flow rates are converted into metal temperature transients and thermal stresses. However, since the GFT is based upon the linear superposition principle, it cannot be directly used in the case of varying heat transfer coefficients. In the present work, a different methodology is proposed, based upon CMS for temperature transient calculation and upon the GFT for the related thermal stress evaluation. This new approach allows variable heat transfer coefficients to be accounted for. The methodology is applied for two different case studies, taken from the literature: a thick pipe and a nozzle connected to a spherical head, both subjected to multiple convective boundary conditions

Uniqueness theorems for differential pencils with eigenparameter boundary conditions and transmission conditions

Science.gov (United States)

Yang, Chuan-Fu

Inverse spectral problems are considered for differential pencils with boundary conditions depending polynomially on the spectral parameter and with a finite number of transmission conditions. We give formulations of the associated inverse problems such as Titchmarsh-Weyl theorem, Hochstadt-Lieberman theorem and Mochizuki-Trooshin theorem, and prove corresponding uniqueness theorems. The obtained results are generalizations of the similar results for the classical Sturm-Liouville operator on a finite interval.

Simple Navier’s slip boundary condition for the non-Newtonian Lattice Boltzmann fluid dynamics solver

DEFF Research Database (Denmark)

Svec, Oldrich; Skoček, Jan

2013-01-01

The ability of the Lattice Boltzmann method, as the fluid dynamics solver, to properly simulate macroscopic Navier’s slip boundary condition is investigated. An approximate equation relating the Lattice Boltzmann variable slip boundary condition with the macroscopic Navier’s slip boundary condition...

Linear and nonlinear development of controlled disturbances in the supersonic boundary layer on a swept wing at Mach 2.5

International Nuclear Information System (INIS)

Kolosov, G L; Kosinov, A D

2016-01-01

Experimental data on the linear and nonlinear wave train development in 3D supersonic boundary layer over a 45° swept-wing at Mach number 2.5 are presented. Travelling artificial disturbances were introduced in the boundary layer by periodical glow discharge at frequencies 10 and 20 kHz. The spatial-temporal and spectral-wave characteristics of the wave train of unstable disturbances in the linear region are obtained. It is shown that the additional peaks in β '-spectra arise for both subharmonic and fundamental frequencies. The experiments indicate the presence of subharmonic resonance mechanism in 3D boundary layer at Mach number 2.5. (paper)

Nonlinear Schroedinger equation with U(p,q) isotopical group

International Nuclear Information System (INIS)

Makhankov, V.G.; Pashaev, O.K.

1981-01-01

The properties of the nonlinear Schroedinger equation (NLS) with U(1,1) isogroup are considered in detail. This example illustrates the essential difference between the system and the well-known ''vector'' NLS, i.e. the large set of allowed boundary conditions on the fields that leads to a rich set of solutions of the system. Four types of boundary conditions and related soliton solutions are considered. The Bohr-Sommerfeld quantization allows to interpret them in terms of ''drops'' and ''bubbles'' as bound states of a large number of constituent bosons subject to the thermodynamical relations for gas mixtures. The U(1,1) system under the vanishing boundary conditions may be considered as continuous analog of the Hubbard model and therefore the paper is concluded by studying the inverse scattering equations for this case [ru

BOUNDARY CONDITIONS IN GAP GEOMETRY

Energy Technology Data Exchange (ETDEWEB)

Rothenstein, W.; Helholtz, J.

1963-11-15

The procedure for calculnting the monoenergetic angular flux density in lattice cells including voids between fuel and moderator is discussed. Boundary conditions describThe thermal energy of a nuclear reactor may be conserved by using as the reactor coolant a hydrocarbon fraction boiling within the range 220 to 650 deg C (preferably 340 to 550 deg C) and containing not more than 5% of extraneous materials having neutron cross sections of > 10 barns. The hot coolant may either be cracked outside of the reactor or used to heat another petroleum hydrocarbon which is to be converted. (D.L.C.)

Nonlinear Ion Harmonics in the Paul Trap with Added Octopole Field: Theoretical Characterization and New Insight into Nonlinear Resonance Effect.

Science.gov (United States)

Xiong, Caiqiao; Zhou, Xiaoyu; Zhang, Ning; Zhan, Lingpeng; Chen, Yongtai; Nie, Zongxiu

2016-02-01

The nonlinear harmonics within the ion motion are the fingerprint of the nonlinear fields. They are exclusively introduced by these nonlinear fields and are responsible to some specific nonlinear effects such as nonlinear resonance effect. In this article, the ion motion in the quadrupole field with a weak superimposed octopole component, described by the nonlinear Mathieu equation (NME), was studied by using the analytical harmonic balance (HB) method. Good accuracy of the HB method, which was comparable with that of the numerical fourth-order Runge-Kutta (4th RK), was achieved in the entire first stability region, except for the points at the stability boundary (i.e., β = 1) and at the nonlinear resonance condition (i.e., β = 0.5). Using the HB method, the nonlinear 3β harmonic series introduced by the octopole component and the resultant nonlinear resonance effect were characterized. At nonlinear resonance, obvious resonant peaks were observed in the nonlinear 3β series of ion motion, but were not found in the natural harmonics. In addition, both resonant excitation and absorption peaks could be observed, simultaneously. These are two unique features of the nonlinear resonance, distinguishing it from the normal resonance. Finally, an approximation equation was given to describe the corresponding working parameter, q nr , at nonlinear resonance. This equation can help avoid the sensitivity degradation due to the operation of ion traps at the nonlinear resonance condition.

Atmospheric-radiation boundary conditions for high-frequency waves in time-distance helioseismology

Science.gov (United States)

Fournier, D.; Leguèbe, M.; Hanson, C. S.; Gizon, L.; Barucq, H.; Chabassier, J.; Duruflé, M.

2017-12-01

The temporal covariance between seismic waves measured at two locations on the solar surface is the fundamental observable in time-distance helioseismology. Above the acoustic cut-off frequency ( 5.3 mHz), waves are not trapped in the solar interior and the covariance function can be used to probe the upper atmosphere. We wish to implement appropriate radiative boundary conditions for computing the propagation of high-frequency waves in the solar atmosphere. We consider recently developed and published radiative boundary conditions for atmospheres in which sound-speed is constant and density decreases exponentially with radius. We compute the cross-covariance function using a finite element method in spherical geometry and in the frequency domain. The ratio between first- and second-skip amplitudes in the time-distance diagram is used as a diagnostic to compare boundary conditions and to compare with observations. We find that a boundary condition applied 500 km above the photosphere and derived under the approximation of small angles of incidence accurately reproduces the "infinite atmosphere" solution for high-frequency waves. When the radiative boundary condition is applied 2 Mm above the photosphere, we find that the choice of atmospheric model affects the time-distance diagram. In particular, the time-distance diagram exhibits double-ridge structure when using a Vernazza Avrett Loeser atmospheric model.

Mixed problems for linear symmetric hyperbolic systems with characteristic boundary conditions

International Nuclear Information System (INIS)

Secchi, P.

1994-01-01

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

Flow boundary conditions for chain-end adsorbing polymer blends.

Science.gov (United States)

Zhou, Xin; Andrienko, Denis; Delle Site, Luigi; Kremer, Kurt

2005-09-08

Using the phenol-terminated polycarbonate blend as an example, we demonstrate that the hydrodynamic boundary conditions for a flow of an adsorbing polymer melt are extremely sensitive to the structure of the epitaxial layer. Under shear, the adsorbed parts (chain ends) of the polymer melt move along the equipotential lines of the surface potential whereas the adsorbed additives serve as the surface defects. In response to the increase of the number of the adsorbed additives the surface layer becomes thinner and solidifies. This results in a gradual transition from the slip to the no-slip boundary condition for the melt flow, with a nonmonotonic dependence of the slip length on the surface concentration of the adsorbed ends.

Quantum osp-invariant non-linear Schroedinger equation

International Nuclear Information System (INIS)

Kulish, P.P.

1985-04-01

The generalizations of the non-linear Schroedinger equation (NS) associated with the orthosymplectic superalgebras are formulated. The simplest osp(1/2)-NS model is solved by the quantum inverse scattering method on a finite interval under periodic boundary conditions as well as on the wholeline in the case of a finite number of excitations. (author)

Scalar field dynamics in a BTZ background with generic boundary conditions

Energy Technology Data Exchange (ETDEWEB)

Garbarz, Alan; La Madrid, Joan [UBA y IFIBA, CONICET, Departamento de Fisica, FCEyN, Buenos Aires (Argentina); Leston, Mauricio [Pabellon IAFE-CONICET, Instituto de Astronomia y Fisica del Espacio, Buenos Aires (Argentina)

2017-11-15

We revisit the dynamics of a massive scalar field in a Banados, Teitelboim, and Zanelli background taking into account the lack of global hyperbolicity of the spacetime. We approach this issue using the strategy of Ishibashi and Wald which finds a unique smooth solution as the causal evolution of initial data, each possible evolution corresponding to a positive self-adjoint extension of certain operator in a Hilbert space on the initial surface. Moreover, solutions obtained this way are the most general ones satisfying a few physically sensible requirements. This procedure is intimately related to the choice of boundary conditions and the existence of bound states. We find that the scalar field dynamics in the (effective) mass window -3/4 ≤ m{sub e}{sup 2}l{sup 2} < 0 can be well defined within a one-parametric family of distinct boundary conditions (-3/4 being the conformally coupled case), while for m{sub e}{sup 2}l{sup 2} ≥ 0 the boundary condition is unique (only one self-adjoint extension is possible). It is argued that there is no sensible evolution possible for m{sub e}{sup 2}l{sup 2} < -1, and also it is shown that in the range m{sub e}{sup 2}l{sup 2} element of [-1, -3/4) there is a U(1) family of allowed boundary conditions, however, the positivity of the self-adjoint extensions is only motivated but not proven. We focus mainly on describing the dynamics of such evolutions given the initial data and all possible boundary conditions, and in particular we show the energy is always positive and conserved. (orig.)

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

Directory of Open Access Journals (Sweden)

Yurii M. Streliaiev

2016-06-01

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

Enhancement of single mode operation in coaxial optical waveguide using DB boundary conditions

Science.gov (United States)

Lohia, Pooja; Prajapati, Y.; Saini, J. P.; Rai, B. S.

2014-11-01

In this study, a competent numerical strategy to compute the dispersion of optical waveguides is presented and propagation of electromagnetic waves in a coaxial optical waveguide with DB boundary conditions is instigated. For this intend, cylindrical coordinates are here being used to derive the DB boundary conditions and to obtain field components for the modes. The propagation constant for the waveguide to be studied is determined by solving the Bessel and the modified Bessel functions. The cutoff frequencies for various lower order modes have been calculated and their dispersion characteristics are plotted correspondingly. The behavior of the coaxial optical waveguide under DB boundary conditions is shown to be significantly different from that of coaxial optical waveguide and conventional optical waveguide under traditional or tangential boundary conditions. Finally, the effect of waveguide dimensions on the mode cutoff frequencies and fabrication issues are also addressed.

Hydrogeological boundary settings in SR 97. Uncertainties in regional boundary settings and transfer of boundary conditions to site-scale models

International Nuclear Information System (INIS)

Follin, S.

1999-06-01

The SR 97 project presents a performance assessment (PA) of the overall safety of a hypothetical deep repository at three sites in Sweden arbitrarily named Aberg, Beberg and Ceberg. One component of this PA assesses the uncertainties in the hydrogeological modelling. This study focuses on uncertainties in boundary settings (size of model domain and boundary conditions) in the regional and site-scale hydrogeological modelling of the three sites used to simulating the possible transport of radionuclides from the emplacement waste packages through the host rock to the accessible environment. Model uncertainties associated with, for instance, parameter heterogeneity and structural interpretations are addressed in other studies. This study concludes that the regional modelling of the SR 97 project addresses uncertainties in the choice of boundary conditions and size of model domain differently at each site, although the overall handling is acceptable and in accordance with common modelling practice. For example, the treatment of uncertainties with regard to the ongoing post-glacial flushing of the Baltic Shield is creditably addressed although not exhaustive from a modelling point of view. A significant contribution of the performed modelling is the study of nested numerical models, i.e., the numerical interplay between regional and site-scale numerical models. In the site-scale modelling great efforts are made to address problems associated with (i) the telescopic mesh refinement (TMR) technique with regard to the stochastic continuum approach, and (ii) the transfer of boundary conditions between variable-density flow systems and flow systems that are constrained to treat uniform density flow. This study concludes that the efforts made to handle these problems are acceptable with regards to the objectives of the SR 97 project

Electron temperature gradient mode instability and stationary vortices with elliptic and circular boundary conditions in non-Maxwellian plasmas

Energy Technology Data Exchange (ETDEWEB)

Haque, Q. [Theoretical Physics Division, PINSTECH, P.O. Nilore, Islamabad (Pakistan); Zakir, U. [Department of Physics, University of Peshawar, Khyber Pakhtun Khwa 25000 (Pakistan); Department of Physics, University of Malakand, Khyber Pakhtun Khwa 18800 (Pakistan); Qamar, A. [Department of Physics, University of Peshawar, Khyber Pakhtun Khwa 25000 (Pakistan)

2015-12-15

Linear and nonlinear dynamics of electron temperature gradient mode along with parallel electron dynamics is investigated by considering hydrodynamic electrons and non-Maxwellian ions. It is noticed that the growth rate of η{sub e}-mode driven linear instability decreases by increasing the value of spectral index and increases by reducing the ion/electron temperature ratio along the magnetic field lines. The eigen mode dispersion relation is also found in the ballooning mode limit. Stationary solutions in the form of dipolar vortices are obtained for both circular and elliptic boundary conditions. It is shown that the dynamics of both circular and elliptic vortices changes with the inclusion of inhomogeneity and non-Maxwellian effects.

Electron temperature gradient mode instability and stationary vortices with elliptic and circular boundary conditions in non-Maxwellian plasmas

Science.gov (United States)

Haque, Q.; Zakir, U.; Qamar, A.

2015-12-01

Linear and nonlinear dynamics of electron temperature gradient mode along with parallel electron dynamics is investigated by considering hydrodynamic electrons and non-Maxwellian ions. It is noticed that the growth rate of ηe-mode driven linear instability decreases by increasing the value of spectral index and increases by reducing the ion/electron temperature ratio along the magnetic field lines. The eigen mode dispersion relation is also found in the ballooning mode limit. Stationary solutions in the form of dipolar vortices are obtained for both circular and elliptic boundary conditions. It is shown that the dynamics of both circular and elliptic vortices changes with the inclusion of inhomogeneity and non-Maxwellian effects.

The Impact of Dam-Reservoir-Foundation Interaction on Nonlinear Response of Concrete Gravity Dams

International Nuclear Information System (INIS)

Amini, Ali Reza; Motamedi, Mohammad Hossein; Ghaemian, Mohsen

2008-01-01

To study the impact of dam-reservoir-foundation interaction on nonlinear response of concrete gravity dams, a two-dimensional finite element model of a concrete gravity dam including the dam body, a part of its foundation and a part of the reservoir was made. In addition, the proper boundary conditions were used in both reservoir and foundation in order to absorb the energy of outgoing waves at the far end boundaries. Using the finite element method and smeared crack approach, some different seismic nonlinear analyses were done and finally, we came to a conclusion that the consideration of dam-reservoir-foundation interaction in nonlinear analysis of concrete dams is of great importance, because from the performance point of view, this interaction significantly improves the nonlinear response of concrete dams

Boundary conditions and dualities: vector fields in AdS/CFT

International Nuclear Information System (INIS)

Marolf, Donald; Ross, Simo F.

2006-01-01

In AdS, scalar fields with masses slightly above the Breitenlohner-Freedman bound admit a variety of possible boundary conditions which are reflected in the Lagrangian of the dual field theory. Generic small changes in the AdS boundary conditions correspond to deformations of the dual field theory by multi-trace operators. Here we extend this discussion to the case of vector gauge fields in the bulk spacetime using the results of Ishibashi and Wald [hep-th/0402184]. As in the context of scalar fields, general boundary conditions for vector fields involve multi-trace deformations which lead to renormalization-group flows. Such flows originate in ultra-violet CFTs which give new gauge/gravity dualities. At least for AdS 4 /CFT 3 , the dual of the bulk photon appears to be a propagating gauge field instead of the usual R-charge current. Applying similar reasoning to tensor fields suggests the existence of a duality between string theory on AdS 4 and a quantum gravity theory in three dimensions

Nonlinear effects of high temperature on buckling of structural elements

International Nuclear Information System (INIS)

Iyengar, N.G.R.

1975-01-01

Structural elements used in nuclear reactors are subjected to high temperatures. Since with increase in temperature there is a gradual fall in the elastic modulus and the stress-strain relationship is nonlinear at these operating load levels, a realistic estimate of the buckling load should include this nonlinearity. In this paper the buckling loads for uniform columns with circular and rectangular cross-sections and different boundary conditions under high temperature environment are estimated. The stress-strain relationship for the material has been assumed to follow inverse Ramberg-Osgood law. In view of the fact that no closed form solutions are possible, approximate methods like perturbation and Galerkin techniques are used. Further, the solution for general value for 'm' is quite involved. Results have been obtained with values for 'm' as 3 and 5. Studies reveal that the influence of material nonlinearity on the buckling load is of the softening type, and it increases with increase in the value of 'm'. The nonlinear effects are more for clamped boundaries than for simply supported boundaries. For the first mode analysis both the methods are powerful. It is, however, felt that for higher modes the Galerkin method might be better in view of its simplicity. This investigation may be considered as a step towards a more general solution

Transducer placement for robustness to variations in boundary conditions for active structural acoustic control

Science.gov (United States)

Sprofera, Joseph D.; Clark, Robert L.; Cabell, Randolph H.; Gibbs, Gary P.

2005-05-01

Turbulent boundary layer (TBL) noise is considered a primary contribution to the interior noise present in commercial airliners. There are numerous investigations of interior noise control devoted to aircraft panels; however, practical realization is a potential challenge since physical boundary conditions are uncertain at best. In most prior studies, pinned or clamped boundary conditions were assumed; however, realistic panels likely display a range of boundary conditions between these two limits. Uncertainty in boundary conditions is a challenge for control system designers, both in terms of the compensator implemented and the location of transducers required to achieve the desired control. The impact of model uncertainties, specifically uncertain boundaries, on the selection of transducer locations for structural acoustic control is considered herein. The final goal of this work is the design of an aircraft panel structure that can reduce TBL noise transmission through the use of a completely adaptive, single-input, single-output control system. The feasibility of this goal is demonstrated through the creation of a detailed analytical solution, followed by the implementation of a test model in a transmission loss apparatus. Successfully realizing a control system robust to variations in boundary conditions can lead to the design and implementation of practical adaptive structures that could be used to control the transmission of sound to the interior of aircraft. Results from this research effort indicate it is possible to optimize the design of actuator and sensor location and aperture, minimizing the impact of boundary conditions on the desired structural acoustic control.

Homogenized boundary conditions and resonance effects in Faraday cages

OpenAIRE

Hewett, DP; Hewitt, IJ

2016-01-01

We present a mathematical study of two-dimensional electrostatic and electromagnetic shielding by a cage of conducting wires (the so-called `Faraday cage e ect'). Taking the limit as the number of wires in the cage tends to in nity we use the asymptotic method of multiple scales to derive continuum models for the shielding, involving homogenized boundary conditions on an e ective cage boundary. We show how the resulting models depend on key cage parameters such as the...

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

Science.gov (United States)

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

2018-05-01

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

Existence and Nonexistence of Positive Solutions for Coupled Riemann-Liouville Fractional Boundary Value Problems

Directory of Open Access Journals (Sweden)

Johnny Henderson

2016-01-01

Full Text Available We investigate the existence and nonexistence of positive solutions for a system of nonlinear Riemann-Liouville fractional differential equations with two parameters, subject to coupled integral boundary conditions.

Boundary conditions for quasiclassical Green's function for superfluid Fermi systems

International Nuclear Information System (INIS)

Nagai, K.; Hara, J.

1988-01-01

The authors show that the quasiclassical Green's Function for Fermi liquids can be constructed from the solutions of the Bogoliubov-de Gennes equation within the Andreev approximation and derive self-consistent relations to be satisfied by the quasiclassical Green's function at the surfaces. The so-called normalization condition for the quasiclassical Green's function is obtained from this self-consistent relation. They consider a specularly reflecting wall, a randomly rippled wall, and a proximity boundary as model surfaces. Their boundary condition for the randomly rippled wall is different from that derived by Buchholtz and Rainer and Buchholtz

The neutron transport with general boundary conditions (II)

International Nuclear Information System (INIS)

Boulanouar, Mohamed

2012-01-01

This Note deals with the one-dimensional transport operator, on an unbounded domain, endowed with general boundary conditions. We show the generation of a strongly continuous semigroup and we study its spectral properties. In particular, we prove the existence of a leading eigenvalue. (author)

Mixing in thermally stratified nonlinear spin-up with uniform boundary fluxes

International Nuclear Information System (INIS)

Baghdasarian, Meline; Pacheco-Vega, Arturo; Pacheco, J. Rafael; Verzicco, Roberto

2014-01-01

Studies of stratified spin-up experiments in enclosed cylinders have reported the presence of small pockets of well-mixed fluids but quantitative measurements of the mixedness of the fluid has been lacking. Previous numerical simulations have not addressed these measurements. Here we present numerical simulations that explain how the combined effect of spin-up and thermal boundary conditions enhances or hinders mixing of a fluid in a cylinder. The energy of the system is characterized by splitting the potential energy into diabatic and adiabatic components, and measurements of efficiency of mixing are based on both, the ratio of dissipation of available potential energy to forcing and variance of temperature. The numerical simulations of the Navier–Stokes equations for the problem with different sets of thermal boundary conditions at the horizontal walls helped shed some light on the physical mechanisms of mixing, for which a clear explanation was absent

Mixing in thermally stratified nonlinear spin-up with uniform boundary fluxes

Energy Technology Data Exchange (ETDEWEB)

Baghdasarian, Meline; Pacheco-Vega, Arturo [Department of Mechanical Engineering, California State University, Los Angeles, Los Angeles, California 90032 (United States); Pacheco, J. Rafael, E-mail: rpacheco@asu.edu [SAP Americas Inc., Scottsdale, Arizona 85251 (United States); School of Mathematical and Statistical Sciences, Arizona State University, Tempe, Arizona 85287 (United States); Environmental Fluid Dynamics Laboratories, Department of Civil Engineering and Geological Sciences, The University of Notre Dame, South Bend, Indiana 46556 (United States); Verzicco, Roberto [Dipartimento di Ingegneria Meccanica, Universita di Roma “Tor Vergata”, Via del Politecnico 1, 00133, Roma (Italy); PoF, University of Twente, 7500 AE Enschede (Netherlands)

2014-09-15

Studies of stratified spin-up experiments in enclosed cylinders have reported the presence of small pockets of well-mixed fluids but quantitative measurements of the mixedness of the fluid has been lacking. Previous numerical simulations have not addressed these measurements. Here we present numerical simulations that explain how the combined effect of spin-up and thermal boundary conditions enhances or hinders mixing of a fluid in a cylinder. The energy of the system is characterized by splitting the potential energy into diabatic and adiabatic components, and measurements of efficiency of mixing are based on both, the ratio of dissipation of available potential energy to forcing and variance of temperature. The numerical simulations of the Navier–Stokes equations for the problem with different sets of thermal boundary conditions at the horizontal walls helped shed some light on the physical mechanisms of mixing, for which a clear explanation was absent.

Revisiting Johnson and Jackson boundary conditions for granular flows

Energy Technology Data Exchange (ETDEWEB)

Li, Tingwen; Benyahia, Sofiane

2012-07-01

In this article, we revisit Johnson and Jackson boundary conditions for granular flows. The oblique collision between a particle and a flat wall is analyzed by adopting the classic rigid-body theory and a more realistic semianalytical model. Based on the kinetic granular theory, the input parameter for the partial-slip boundary conditions, specularity coefficient, which is not measurable in experiments, is then interpreted as a function of the particle-wall restitution coefficient, the frictional coefficient, and the normalized slip velocity at the wall. An analytical expression for the specularity coefficient is suggested for a flat, frictional surface with a low frictional coefficient. The procedure for determining the specularity coefficient for a more general problem is outlined, and a working approximation is provided.

Airborne observations of newly formed boundary layer aerosol particles under cloudy conditions

Directory of Open Access Journals (Sweden)

B. Altstädter

2018-06-01

Full Text Available This study describes the appearance of ultrafine boundary layer aerosol particles under classical non-favourable conditions at the research site of TROPOS (Leibniz Institute for Tropospheric Research. Airborne measurements of meteorological and aerosol properties of the atmospheric boundary layer (ABL were repeatedly performed with the unmanned aerial system ALADINA (Application of Light-weight Aircraft for Detecting IN-situ Aerosol during three seasons between October 2013 and July 2015. More than 100 measurement flights were conducted on 23 different days with a total flight duration of 53 h. In 26 % of the cases, maxima of ultrafine particles were observed close to the inversion layer at altitudes between 400 and 600 m and the particles were rapidly mixed vertically and mainly transported downwards during short time intervals of cloud gaps. This study focuses on two measurement days affected by low-level stratocumulus clouds, but different wind directions (NE, SW and minimal concentrations (< 4.6 µg m−3 of SO2, as a common indicator for precursor gases at ground. Taken from vertical profiles, the onset of clouds led to a non-linearity of humidity that resulted in an increased turbulence at the local-scale and caused fast nucleation e.g., but in relation to rapid dilution of surrounding air, seen in sporadic clusters of ground data, so that ultrafine particles disappeared in the verticality. The typical banana shape of new particle formation (NPF and growth was not seen at ground and thus these days might not have been classified as NPF event days by pure surface studies.

Adaptive Neural Control of Nonaffine Nonlinear Systems without Differential Condition for Nonaffine Function

Directory of Open Access Journals (Sweden)

Chaojiao Sun

2016-01-01

Full Text Available An adaptive neural control scheme is proposed for nonaffine nonlinear system without using the implicit function theorem or mean value theorem. The differential conditions on nonaffine nonlinear functions are removed. The control-gain function is modeled with the nonaffine function probably being indifferentiable. Furthermore, only a semibounded condition for nonaffine nonlinear function is required in the proposed method, and the basic idea of invariant set theory is then constructively introduced to cope with the difficulty in the control design for nonaffine nonlinear systems. It is rigorously proved that all the closed-loop signals are bounded and the tracking error converges to a small residual set asymptotically. Finally, simulation examples are provided to demonstrate the effectiveness of the designed method.

Non-linear analysis of wave progagation using transform methods and plates and shells using integral equations

Science.gov (United States)

Pipkins, Daniel Scott

Two diverse topics of relevance in modern computational mechanics are treated. The first involves the modeling of linear and non-linear wave propagation in flexible, lattice structures. The technique used combines the Laplace Transform with the Finite Element Method (FEM). The procedure is to transform the governing differential equations and boundary conditions into the transform domain where the FEM formulation is carried out. For linear problems, the transformed differential equations can be solved exactly, hence the method is exact. As a result, each member of the lattice structure is modeled using only one element. In the non-linear problem, the method is no longer exact. The approximation introduced is a spatial discretization of the transformed non-linear terms. The non-linear terms are represented in the transform domain by making use of the complex convolution theorem. A weak formulation of the resulting transformed non-linear equations yields a set of element level matrix equations. The trial and test functions used in the weak formulation correspond to the exact solution of the linear part of the transformed governing differential equation. Numerical results are presented for both linear and non-linear systems. The linear systems modeled are longitudinal and torsional rods and Bernoulli-Euler and Timoshenko beams. For non-linear systems, a viscoelastic rod and Von Karman type beam are modeled. The second topic is the analysis of plates and shallow shells under-going finite deflections by the Field/Boundary Element Method. Numerical results are presented for two plate problems. The first is the bifurcation problem associated with a square plate having free boundaries which is loaded by four, self equilibrating corner forces. The results are compared to two existing numerical solutions of the problem which differ substantially. non-linear model are compared to those

Three Boundary Conditions for Computing the Fixed-Point Property in Binary Mixture Data.

Directory of Open Access Journals (Sweden)

Leendert van Maanen

Full Text Available The notion of "mixtures" has become pervasive in behavioral and cognitive sciences, due to the success of dual-process theories of cognition. However, providing support for such dual-process theories is not trivial, as it crucially requires properties in the data that are specific to mixture of cognitive processes. In theory, one such property could be the fixed-point property of binary mixture data, applied-for instance- to response times. In that case, the fixed-point property entails that response time distributions obtained in an experiment in which the mixture proportion is manipulated would have a common density point. In the current article, we discuss the application of the fixed-point property and identify three boundary conditions under which the fixed-point property will not be interpretable. In Boundary condition 1, a finding in support of the fixed-point will be mute because of a lack of difference between conditions. Boundary condition 2 refers to the case in which the extreme conditions are so different that a mixture may display bimodality. In this case, a mixture hypothesis is clearly supported, yet the fixed-point may not be found. In Boundary condition 3 the fixed-point may also not be present, yet a mixture might still exist but is occluded due to additional changes in behavior. Finding the fixed-property provides strong support for a dual-process account, yet the boundary conditions that we identify should be considered before making inferences about underlying psychological processes.

The Hardy inequality with boundary or intermediate conditions

Czech Academy of Sciences Publication Activity Database

Kufner, Alois

2017-01-01

Roč. 8, č. 2 (2017), s. 105-109 ISSN 2077-9879 Institutional support: RVO:67985840 Keywords : Hardy's inequality * boundary conditions Subject RIV: BA - General Mathematics OBOR OECD: Pure mathematics http://www.mathnet.ru/ php /archive.phtml?wshow=paper&jrnid=emj&paperid=259&option_lang=eng

The Hardy inequality with boundary or intermediate conditions

Czech Academy of Sciences Publication Activity Database

Kufner, Alois

2017-01-01

Roč. 8, č. 2 (2017), s. 105-109 ISSN 2077-9879 Institutional support: RVO:67985840 Keywords : Hardy's inequality * boundary conditions Subject RIV: BA - General Mathematics OBOR OECD: Pure mathematics http://www.mathnet.ru/php/archive.phtml?wshow=paper&jrnid=emj&paperid=259&option_lang=eng

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.

Boundary-value problems with free boundaries for elliptic systems of equations

CERN Document Server

Monakhov, V N

1983-01-01

This book is concerned with certain classes of nonlinear problems for elliptic systems of partial differential equations: boundary-value problems with free boundaries. The first part has to do with the general theory of boundary-value problems for analytic functions and its applications to hydrodynamics. The second presents the theory of quasiconformal mappings, along with the theory of boundary-value problems for elliptic systems of equations and applications of it to problems in the mechanics of continuous media with free boundaries: problems in subsonic gas dynamics, filtration theory, and problems in elastico-plasticity.

Framing the impact of external magnetic field on bioconvection of a nanofluid flow containing gyrotactic microorganisms with convective boundary conditions

Directory of Open Access Journals (Sweden)

Tanmoy Chakraborty

2018-03-01

Full Text Available The intention of this study is to examine the combined impacts of magnetic field and convective boundary state on bioconvection of a nanofluid flow along an expanding sheet co-existed with gyrotactic microorganisms. The fundamental partial differential equations are reduced to a set of nonlinear ordinary differential equations taking a guide of some appropriate similarity transformations. The numerical fallouts are calculated by considering the bvp4c function of Matlab. The impacts of magnetic field parameter, surface convection parameter, Eckert number and Peclet number on non-dimensional velocity, nanoparticle concentration, temperature and density of self-moving microorganisms are interpreted through graphs and charts. The fluid velocity near the surface and the Nusselt number is lessen with magnetic field. Surface convection parameter enhances the self-moving microorganism flux but a reverse result is noticed for Peclet number. Also, the contrast between the present results with formerly visited outcomes is in excellent harmony. Keywords: Nanofluid, Bioconvection, Gyrotactic microorganisms, Magnetic field, Convective boundary condition

Analysis of activation energy in Couette-Poiseuille flow of nanofluid in the presence of chemical reaction and convective boundary conditions

Science.gov (United States)

Zeeshan, A.; Shehzad, N.; Ellahi, R.

2018-03-01

The motivation of the current article is to explore the energy activation in MHD radiative Couette-Poiseuille flow nanofluid in horizontal channel with convective boundary conditions. The mathematical model of Buongiorno [1] effectively describes the current flow analysis. Additionally, the impact of chemical reaction is also taken in account. The governing flow equations are simplified with the help of boundary layer approximations. Non-linear coupled equations for momentum, energy and mass transfer are tackled with analytical (HAM) technique. The influence of dimensionless convergence parameter like Brownian motion parameter, radiation parameter, buoyancy ratio parameter, dimensionless activation energy, thermophoresis parameter, temperature difference parameter, dimensionless reaction rate, Schmidt number, Brinkman number, Biot number and convection diffusion parameter on velocity, temperature and concentration profiles are discussed graphically and in tabular form. From the results, it is elaborate that the nanoparticle concentration is directly proportional to the chemical reaction with activation energy and the performance of Brownian motion on nanoparticle concentration gives reverse pattern to that of thermophoresis parameter.

A One-Dimensional Wave Equation with White Noise Boundary Condition

International Nuclear Information System (INIS)

Kim, Jong Uhn

2006-01-01

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

Dirac perturbations on Schwarzschild-anti-de Sitter spacetimes: Generic boundary conditions and new quasinormal modes

Science.gov (United States)

Wang, Mengjie; Herdeiro, Carlos; Jing, Jiliang

2017-11-01

We study Dirac quasinormal modes of Schwarzschild-anti-de Sitter (Schwarzschild-AdS) black holes, following the generic principle for allowed boundary conditions proposed in [M. Wang, C. Herdeiro, and M. O. P. Sampaio, Phys. Rev. D 92, 124006 (2015)., 10.1103/PhysRevD.92.124006]. After deriving the equations of motion for Dirac fields on the aforementioned background, we impose vanishing energy flux boundary conditions to solve these equations. We find a set of two Robin boundary conditions are allowed. These two boundary conditions are used to calculate Dirac normal modes on empty AdS and quasinormal modes on Schwarzschild-AdS black holes. In the former case, we recover the known normal modes of empty AdS; in the latter case, the two sets of Robin boundary conditions lead to two different branches of quasinormal modes. The impact on these modes of the black hole size, the angular momentum quantum number and the overtone number are discussed. Our results show that vanishing energy flux boundary conditions are a robust principle, applicable not only to bosonic fields but also to fermionic fields.

Quantum communication through a spin ring with twisted boundary conditions

International Nuclear Information System (INIS)

Bose, S.; Jin, B.-Q.; Korepin, V.E.

2005-01-01

We investigate quantum communication between the sites of a spin ring with twisted boundary conditions. Such boundary conditions can be achieved by a magnetic flux through the ring. We find that a nonzero twist can improve communication through finite odd-numbered rings and enable high-fidelity multiparty quantum communication through spin rings (working near perfectly for rings of five and seven spins). We show that in certain cases, the twist results in the complete blockage of quantum-information flow to a certain site of the ring. This effect can be exploited to interface and entangle a flux qubit and a spin qubit without embedding the latter in a magnetic field

The effects of external conditions in turbulent boundary layers

Science.gov (United States)

Brzek, Brian G.

The effects of multiple external conditions on turbulent boundary layers were studied in detail. These external conditions include: surface roughness, upstream turbulence intensity, and pressure gradient. Furthermore, the combined effects of these conditions show the complicated nature of many realistic flow conditions. It was found that the effects of surface roughness are difficult to generalize, given the importance of so many parameters. These parameters include: roughness geometry, roughness regime, roughness height to boundary layer thickness, (k/delta), roughness parameter, ( k+), Reynolds number, and roughness function (Delta B+). A further complication, is the difficulty in computing the wall shear stress, tauw/rho. For the sand grain type roughness, the mean velocity and Reynolds stresses were studied in inner and outer variables, as well as, boundary layer parameters, anisotropy tensor, production term, and viscous stress and form drag contributions. To explore the effects of roughness and Reynolds number dependence in the boundary layer, a new experiment was carefully designed to properly capture the x-dependence of the single-point statistics. It was found that roughness destroys the viscous layer near the wall, thus, reducing the contribution of the viscous stress in the wall region. As a result, the contribution in the skin friction due to form drag increases, while the viscous stress decreases. This yields Reynolds number invariance in the skin friction, near-wall roughness parameters, and inner velocity profiles as k + increases into the fully rough regime. However, in the transitionally rough regime, (i.e., 5 component shows the largest influence of roughness, where the high peak near the wall was decreased and became nearly flat for the fully rough regime profiles. In addition, the Reynolds stresses in outer variables show self-similarity for fixed experimental conditions. However, as the roughness parameter, k +, increases, all Reynolds stress

The existence of positive solutions for nonlinear boundary system with $p$-Laplacian operator based on sign-changing nonlinearities

Directory of Open Access Journals (Sweden)

Fuyi Xu

2010-12-01

(\\phi_{p_1}(u''+a_1(tf(u,v=0, 01, i=1,2$. We obtain some sufficient conditions for the existence of two positive solutions or infinitely many positive solutions by using a fixed-point theorem in cones. Especially, the nonlinear terms $f,g $ are allowed to change sign. The conclusions essentially extend and improve the known results.

Symbolic computation of analytic approximate solutions for nonlinear fractional differential equations

Science.gov (United States)

Lin, Yezhi; Liu, Yinping; Li, Zhibin

2013-01-01

The Adomian decomposition method (ADM) is one of the most effective methods to construct analytic approximate solutions for nonlinear differential equations. In this paper, based on the new definition of the Adomian polynomials, Rach (2008) [22], the Adomian decomposition method and the Padé approximants technique, a new algorithm is proposed to construct analytic approximate solutions for nonlinear fractional differential equations with initial or boundary conditions. Furthermore, a MAPLE software package is developed to implement this new algorithm, which is user-friendly and efficient. One only needs to input the system equation, initial or boundary conditions and several necessary parameters, then our package will automatically deliver the analytic approximate solutions within a few seconds. Several different types of examples are given to illustrate the scope and demonstrate the validity of our package, especially for non-smooth initial value problems. Our package provides a helpful and easy-to-use tool in science and engineering simulations. Program summaryProgram title: ADMP Catalogue identifier: AENE_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AENE_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 12011 No. of bytes in distributed program, including test data, etc.: 575551 Distribution format: tar.gz Programming language: MAPLE R15. Computer: PCs. Operating system: Windows XP/7. RAM: 2 Gbytes Classification: 4.3. Nature of problem: Constructing analytic approximate solutions of nonlinear fractional differential equations with initial or boundary conditions. Non-smooth initial value problems can be solved by this program. Solution method: Based on the new definition of the Adomian polynomials [1], the Adomian decomposition method and the Pad

On MHD nonlinear stretching flow of Powell–Eyring nanomaterial

Directory of Open Access Journals (Sweden)

Tasawar Hayat

Full Text Available This communication addresses the magnetohydrodynamic (MHD flow of Powell–Eyring nanomaterial bounded by a nonlinear stretching sheet. Novel features regarding thermophoresis and Brownian motion are taken into consideration. Powell–Eyring fluid is electrically conducted subject to non-uniform applied magnetic field. Assumptions of small magnetic Reynolds number and boundary layer approximation are employed in the mathematical development. Zero nanoparticles mass flux condition at the sheet is selected. Adequate transformation yield nonlinear ordinary differential systems. The developed nonlinear systems have been computed through the homotopic approach. Effects of different pertinent parameters on velocity, temperature and concentration fields are studied and analyzed. Further numerical data of skin friction and heat transfer rate is also tabulated and interpreted. Keywords: Powell–Eyring fluid, Magnetohydrodynamics, Nanomaterial, Nonlinear stretching surface

Eigenstates of a particle in an array of hexagons with periodic boundary condition

Directory of Open Access Journals (Sweden)

A Nemati

2013-10-01

Full Text Available In this paper the problem of a particle in an array of hexagons with periodic boundary condition is solved. Using the projection operators, we categorize eigenfunctions corresponding to each of the irreducible representations of the symmetry group . Based on these results, the Dirichlet and Neumann boundary conditions are discussed.

Molecular Dynamics with Helical Periodic Boundary Conditions

Czech Academy of Sciences Publication Activity Database

Kessler, Jiří; Bouř, Petr

2014-01-01

Roč. 35, č. 21 (2014), s. 1552-1559 ISSN 0192-8651 R&D Projects: GA ČR GAP208/11/0105; GA MŠk(CZ) LH11033 Grant - others:GA AV ČR(CZ) M200551205; GA MŠk(CZ) LM2010005 Institutional support: RVO:61388963 Keywords : periodic boundary conditions * helical symmetry * molecular dynamics * protein structure * amyloid fibrils Subject RIV: CF - Physical ; Theoretical Chemistry Impact factor: 3.589, year: 2014

Boundary conditions in conformal and integrable theories

CERN Document Server

Petkova, V B

2000-01-01

The study of boundary conditions in rational conformal field theories is not only physically important. It also reveals a lot on the structure of the theory ``in the bulk''. The same graphs classify both the torus and the cylinder partition functions and provide data on their hidden ``quantum symmetry''. The Ocneanu triangular cells -- the 3j-symbols of these symmetries, admit various interpretations and make a link between different problems.

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

Boundary terms and junction conditions for the DGP π-Lagrangian and galileon

International Nuclear Information System (INIS)

Dyer, Ethan; Hinterbichler, Kurt

2009-01-01

In the decoupling limit of DGP, π describes the brane-bending degree of freedom. It obeys second order equations of motion, yet it is governed by a higher derivative Lagrangian. We show that, analogously to the Einstein-Hilbert action for GR, the π-Lagrangian requires Gibbons-Hawking-York type boundary terms to render the variational principle well-posed. These terms are important if there are other boundaries present besides the DGP brane, such as in higher dimensional cascading DGP models. We derive the necessary boundary terms in two ways. First, we derive them directly from the brane-localized π-Lagrangian by demanding well-posedness of the action. Second, we calculate them directly from the bulk, taking into account the Gibbons-Hawking-York terms in the bulk Einstein-Hilbert action. As an application, we use the new boundary terms to derive Israel junction conditions for π across a sheet-like source. In addition, we calculate boundary terms and junction conditions for the galileons which generalize the DGP π-Lagrangian, showing that the boundary term for the n-th order galileon is the (n-1)-th order galileon.

On the Existence of a Free Boundary for a Class of Reaction-Diffusion Systems.

Science.gov (United States)

1982-02-01

I. Diaz. "Soluciones con soporte compacto para alguno. problemas semilineales". Collect. Math. 30 (1979), 141-179. -26- [121 J. I. Diaz. Tecnica de ...supersoluciones locales para problemas estacionarios no lineales: applicacion al estudio de flujoe subsonicos. Memory of the Real Academia de Ciencias...nonlinearity, nonlinear boundary conditions, dead core set, chemical reactions Work Unit Number I - Applied Analysis (1) Seccion de Matematicas

First-principle proof of the modified collision boundary conditions for the hard-sphere system

International Nuclear Information System (INIS)

Tessarotto, Massimo; Cremaschini, Claudio

2014-01-01

A fundamental issue lying at the foundation of classical statistical mechanics is the determination of the collision boundary conditions that characterize the dynamical evolution of multi-particle probability density functions (PDF) and are applicable to systems of hard-spheres undergoing multiple elastic collisions. In this paper it is proved that, when the deterministic N-body PDF is included in the class of admissible solutions of the Liouville equation, the customary form of collision boundary conditions adopted in previous literature becomes physically inconsistent and must actually be replaced by suitably modified collision boundary conditions.

Slarti: A boundary condition editor for a coupled climate model

Science.gov (United States)

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

2006-12-01

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

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.

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.

Matrix factorisations for rational boundary conditions by defect fusion

International Nuclear Information System (INIS)

Behr, Nicolas; Fredenhagen, Stefan

2015-01-01

A large class of two-dimensional N=(2,2) superconformal field theories can be understood as IR fixed-points of Landau-Ginzburg models. In particular, there are rational conformal field theories that also have a Landau-Ginzburg description. To understand better the relation between the structures in the rational conformal field theory and in the Landau-Ginzburg theory, we investigate how rational B-type boundary conditions are realised as matrix factorisations in the SU(3)/U(2) Grassmannian Kazama-Suzuki model. As a tool to generate the matrix factorisations we make use of a particular interface between the Kazama-Suzuki model and products of minimal models, whose fusion can be realised as a simple functor on ring modules. This allows us to formulate a proposal for all matrix factorisations corresponding to rational boundary conditions in the SU(3)/U(2) model.

Matrix factorisations for rational boundary conditions by defect fusion

Energy Technology Data Exchange (ETDEWEB)

Behr, Nicolas [Department of Mathematics, Heriot-Watt University,Riccarton, Edinburgh, EH14 4AS (United Kingdom); Maxwell Institute for Mathematical Sciences,Edinburgh (United Kingdom); Fredenhagen, Stefan [Max-Planck-Institut für Gravitationsphysik, Albert-Einstein-Institut,D-14424 Golm (Germany)

2015-05-11

A large class of two-dimensional N=(2,2) superconformal field theories can be understood as IR fixed-points of Landau-Ginzburg models. In particular, there are rational conformal field theories that also have a Landau-Ginzburg description. To understand better the relation between the structures in the rational conformal field theory and in the Landau-Ginzburg theory, we investigate how rational B-type boundary conditions are realised as matrix factorisations in the SU(3)/U(2) Grassmannian Kazama-Suzuki model. As a tool to generate the matrix factorisations we make use of a particular interface between the Kazama-Suzuki model and products of minimal models, whose fusion can be realised as a simple functor on ring modules. This allows us to formulate a proposal for all matrix factorisations corresponding to rational boundary conditions in the SU(3)/U(2) model.

Nonlinear Dielectric Properties of Yeast Cells Cultured in Different Environmental Conditions

Science.gov (United States)

Kawanishi, Gomon; Fukuda, Naoki; Muraji, Masafumi

The harmonics of the electric current through yeast suspensions, the nonlinear dielectric properties of yeast cells, have particular patterns according to the biological activity of the cells and the measurement of these patterns is a technique for determining the activity of living cells. The concentration of glucose and oxygen in yeast culture medium influences the manifestation of fermentation or respiration of yeast cells. Measurements were made with yeast cells (Saccharomyces cerevisiae) cultured aerobically and anaerobically in sufficient glucose concentration, aerobic fermentation and anaerobic fermentation, and aerobically in limited glucose concentration, respiration. The results showed that the harmonics were barely apparent for yeast cells in aerobic fermentation and respiratory; however, cells in the anaerobic fermentation displayed substantial third and fifth harmonics. We can say that environmental condition affects the yeast cells' nonlinear properties, from another viewpoint, the measurements of the nonlinear properties are available to determine the activity of yeast cells adjusted to the conditions of their cultivation.

Comparison of mass transport using average and transient rainfall boundary conditions

International Nuclear Information System (INIS)

Duguid, J.O.; Reeves, M.

1976-01-01

A general two-dimensional model for simulation of saturated-unsaturated transport of radionuclides in ground water has been developed and is currently being tested. The model is being applied to study the transport of radionuclides from a waste-disposal site where field investigations are currently under way to obtain the necessary model parameters. A comparison of the amount of tritium transported is made using both average and transient rainfall boundary conditions. The simulations indicate that there is no substantial difference in the transport for the two conditions tested. However, the values of dispersivity used in the unsaturated zone caused more transport above the water table than has been observed under actual conditions. This deficiency should be corrected and further comparisons should be made before average rainfall boundary conditions are used for long-term transport simulations

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

Directory of Open Access Journals (Sweden)

Idris Addou

2000-01-01

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

Natural frequency analysis of fluid conveying pipeline with different boundary conditions

International Nuclear Information System (INIS)

Huang Yimin; Liu Yongshou; Li Baohui; Li Yanjiang; Yue Zhufeng

2010-01-01

In this study, the natural frequency of fluid-structure interaction in pipeline conveying fluid is investigated by eliminated element-Galerkin method, and the natural frequency equations with different boundary conditions are obtained. Furthermore, the expressions of first natural frequency are simplified in the case of different boundary conditions. Taking into account the Coriolis force, as an example, the natural frequency of a straight pipe simply supported at both ends is studied. In a given boundary condition, the four components (mass, stiffness, length and flow velocity) which relate to the natural frequency of pipeline conveying fluid are studied in detail and the results indicate that the effect of Coriolis force on natural frequency is inappreciable. Then the relationship between natural frequency of the pipeline conveying fluid and that of Euler beam is analyzed. Finally, a dimensionless flow velocity and limit values are presented, and they can be used to estimate the effect of the flow velocity on natural frequency. All the conclusions are well suited to nuclear plant and other industry fields.

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

Impact of quadratic non-linearity on the dynamics of periodic solutions of a wave equation

International Nuclear Information System (INIS)

Kolesov, Andrei Yu; Rozov, Nikolai Kh

2002-01-01

For the non-linear telegraph equation with homogeneous Dirichlet or Neumann conditions at the end-points of a finite interval the question of the existence and the stability of time-periodic solutions bifurcating from the zero equilibrium state is considered. The dynamics of these solutions under a change of the diffusion coefficient (that is, the coefficient of the second derivative with respect to the space variable) is investigated. For the Dirichlet boundary conditions it is shown that this dynamics substantially depends on the presence - or the absence - of quadratic terms in the non-linearity. More precisely, it is shown that a quadratic non-linearity results in the occurrence, under an unbounded decrease of diffusion, of an infinite sequence of bifurcations of each periodic solution. En route, the related issue of the limits of applicability of Yu.S. Kolesov's method of quasinormal forms to the construction of self-oscillations in singularly perturbed hyperbolic boundary value problems is studied

Boundary Conditions, Data Assimilation, and Predictability in Coastal Ocean Models

National Research Council Canada - National Science Library

Samelson, Roger M; Allen, John S; Egbert, Gary D; Kindle, John C; Snyder, Chris

2007-01-01

...: The specific objectives of this research are to determine the impact on coastal ocean circulation models of open ocean boundary conditions from Global Ocean Data Assimilation Experiment (GODAE...

Strongly nonlinear theory of rapid solidification near absolute stability

Science.gov (United States)

Kowal, Katarzyna N.; Altieri, Anthony L.; Davis, Stephen H.

2017-10-01

We investigate the nonlinear evolution of the morphological deformation of a solid-liquid interface of a binary melt under rapid solidification conditions near two absolute stability limits. The first of these involves the complete stabilization of the system to cellular instabilities as a result of large enough surface energy. We derive nonlinear evolution equations in several limits in this scenario and investigate the effect of interfacial disequilibrium on the nonlinear deformations that arise. In contrast to the morphological stability problem in equilibrium, in which only cellular instabilities appear and only one absolute stability boundary exists, in disequilibrium the system is prone to oscillatory instabilities and a second absolute stability boundary involving attachment kinetics arises. Large enough attachment kinetics stabilize the oscillatory instabilities. We derive a nonlinear evolution equation to describe the nonlinear development of the solid-liquid interface near this oscillatory absolute stability limit. We find that strong asymmetries develop with time. For uniform oscillations, the evolution equation for the interface reduces to the simple form f''+(βf')2+f =0 , where β is the disequilibrium parameter. Lastly, we investigate a distinguished limit near both absolute stability limits in which the system is prone to both cellular and oscillatory instabilities and derive a nonlinear evolution equation that captures the nonlinear deformations in this limit. Common to all these scenarios is the emergence of larger asymmetries in the resulting shapes of the solid-liquid interface with greater departures from equilibrium and larger morphological numbers. The disturbances additionally sharpen near the oscillatory absolute stability boundary, where the interface becomes deep-rooted. The oscillations are time-periodic only for small-enough initial amplitudes and their frequency depends on a single combination of physical parameters, including the

On non-linear boundary value problems and parametrisation at multiple nodes

Czech Academy of Sciences Publication Activity Database

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

2016-01-01

Roč. 2016, Č. 80 (2016), s. 1-18 ISSN 1417-3875 Institutional support: RVO:67985840 Keywords : non-local boundary conditions * parametrisation * successive approximations * interval division Subject RIV: BA - General Mathematics Impact factor: 0.926, year: 2016 http://www.math.u-szeged.hu/ejqtde/periodica.html?periodica=1¶mtipus_ertek=publication¶m_ertek=5302

Rate-Independent Processes with Linear Growth Energies and Time-Dependent Boundary Conditions

Czech Academy of Sciences Publication Activity Database

Kružík, Martin; Zimmer, J.

2012-01-01

Roč. 5, č. 3 (2012), s. 591-604 ISSN 1937-1632 R&D Projects: GA AV ČR IAA100750802 Grant - others:GA ČR(CZ) GAP201/10/0357 Institutional research plan: CEZ:AV0Z10750506 Keywords : concentrations * oscillations * time - dependent boundary conditions * rate-independent evolution Subject RIV: BA - General Mathematics http://library.utia.cas.cz/separaty/2011/MTR/kruzik-rate-independent processes with linear growth energies and time - dependent boundary conditions.pdf

Off-wall boundary conditions for turbulent flows obtained from buffer-layer minimal flow units

Science.gov (United States)

Garcia-Mayoral, Ricardo; Pierce, Brian; Wallace, James

2012-11-01

There is strong evidence that the transport processes in the buffer region of wall-bounded turbulence are common across various flow configurations, even in the embryonic turbulence in transition (Park et al., Phys. Fl. 24). We use this premise to develop off-wall boundary conditions for turbulent simulations. Boundary conditions are constructed from DNS databases using periodic minimal flow units and reduced order modeling. The DNS data was taken from a channel at Reτ = 400 and a zero-pressure gradient transitional boundary layer (Sayadi et al., submitted to J . FluidMech .) . Both types of boundary conditions were first tested on a DNS of the core of the channel flow with the aim of extending their application to LES and to spatially evolving flows. 2012 CTR Summer Program.

Nonlinear surface waves at ferrite-metamaterial waveguide structure

Science.gov (United States)

Hissi, Nour El Houda; Mokhtari, Bouchra; Eddeqaqi, Noureddine Cherkaoui; Shabat, Mohammed Musa; Atangana, Jacques

2016-09-01

A new ferrite slab made of a metamaterial (MTM), surrounded by a nonlinear cover cladding and a ferrite substrate, was shown to support unusual types of electromagnetic surface waves. We impose the boundary conditions to derive the dispersion relation and others necessary to formulate the proposed structure. We analyse the dispersion properties of the nonlinear surface waves and we calculate the associated propagation index and the film-cover interface nonlinearity. In the calculation, several sets of the permeability of the MTM are considered. Results show that the waves behaviour depends on the values of the permeability of the MTM, the thickness of the waveguide and the film-cover interface nonlinearity. It is also shown that the use of the singular solutions to the electric field equation allows to identify several new properties of surface waves which do not exist in conventional waveguide.

Construction of the Nuclear Effective Interaction from Energy Eigenstates and Boundary Conditions

Science.gov (United States)

McElvain, Kenneth; Haxton, Wick

2017-01-01

The original Harmonic Oscillator Based Effective Theory (HOBET) work by Haxton and Luu reduced H = T +VNN , with VNN a realistic potential, to Heff in a small basis defined by projection operator P while correctly including all scattering by H through an excluded space Q. Scattering by T is analytically included to all orders, leaving the ET expansion focused on the short range VNN. Results do not depend on the size P as the effect of scattering through Q is fully included, also distinguishing HOBET from other methods. In this talk we abandon VNN and determine the LECs of the ET expansion from energy levels and boundary conditions. In the infinite volume continuum case every energy is an eigenvalue of H with an associated scattering state. In the LQCD context boundary conditions are periodic. In either case the ET LECs can be determined from energy, boundary condition pairs. We show that the Cartesian HO ET LECs can be expressed in terms of the spherical ones, giving a spherical, infinite volume ET, bypassing the use of Luscher's method. The approach cleanly isolates operator mixing induced by the finite box, sequestering effects that vanish in the continuum limit in a Green's function constrained to match the boundary conditions. Supported by the DOE under contracts DE-SC00046548 and DE-AC02-98CH10886.

Global existence and nonexistence for the viscoelastic wave equation with nonlinear boundary damping-source interaction

KAUST Repository

Said-Houari, Belkacem

2012-09-01

The goal of this work is to study a model of the viscoelastic wave equation with nonlinear boundary/interior sources and a nonlinear interior damping. First, applying the Faedo-Galerkin approximations combined with the compactness method to obtain existence of regular global solutions to an auxiliary problem with globally Lipschitz source terms and with initial data in the potential well. It is important to emphasize that it is not possible to consider density arguments to pass from regular to weak solutions if one considers regular solutions of our problem where the source terms are locally Lipschitz functions. To overcome this difficulty, we use an approximation method involving truncated sources and adapting the ideas in [13] to show that the existence of weak solutions can still be obtained for our problem. Second, we show that under some restrictions on the initial data and if the interior source dominates the interior damping term, then the solution ceases to exist and blows up in finite time provided that the initial data are large enough.

Global existence and nonexistence for the viscoelastic wave equation with nonlinear boundary damping-source interaction

KAUST Repository

Said-Houari, Belkacem; Nascimento, Flá vio A Falcã o

2012-01-01

The goal of this work is to study a model of the viscoelastic wave equation with nonlinear boundary/interior sources and a nonlinear interior damping. First, applying the Faedo-Galerkin approximations combined with the compactness method to obtain existence of regular global solutions to an auxiliary problem with globally Lipschitz source terms and with initial data in the potential well. It is important to emphasize that it is not possible to consider density arguments to pass from regular to weak solutions if one considers regular solutions of our problem where the source terms are locally Lipschitz functions. To overcome this difficulty, we use an approximation method involving truncated sources and adapting the ideas in [13] to show that the existence of weak solutions can still be obtained for our problem. Second, we show that under some restrictions on the initial data and if the interior source dominates the interior damping term, then the solution ceases to exist and blows up in finite time provided that the initial data are large enough.

A fast Cauchy-Riemann solver. [differential equation solution for boundary conditions by finite difference approximation

Science.gov (United States)

Ghil, M.; Balgovind, R.

1979-01-01

The inhomogeneous Cauchy-Riemann equations in a rectangle are discretized by a finite difference approximation. Several different boundary conditions are treated explicitly, leading to algorithms which have overall second-order accuracy. All boundary conditions with either u or v prescribed along a side of the rectangle can be treated by similar methods. The algorithms presented here have nearly minimal time and storage requirements and seem suitable for development into a general-purpose direct Cauchy-Riemann solver for arbitrary boundary conditions.

Sublayer of Prandtl Boundary Layers

Science.gov (United States)

Grenier, Emmanuel; Nguyen, Toan T.

2018-03-01

The aim of this paper is to investigate the stability of Prandtl boundary layers in the vanishing viscosity limit {ν \\to 0} . In Grenier (Commun Pure Appl Math 53(9):1067-1091, 2000), one of the authors proved that there exists no asymptotic expansion involving one of Prandtl's boundary layer, with thickness of order {√{ν}} , which describes the inviscid limit of Navier-Stokes equations. The instability gives rise to a viscous boundary sublayer whose thickness is of order {ν^{3/4}} . In this paper, we point out how the stability of the classical Prandtl's layer is linked to the stability of this sublayer. In particular, we prove that the two layers cannot both be nonlinearly stable in L^∞. That is, either the Prandtl's layer or the boundary sublayer is nonlinearly unstable in the sup norm.

Boundary conditions estimation of scoop dosimeter for primary sorting during earthworks

International Nuclear Information System (INIS)

Batij, V.G.; Pravdivyj, A.A.; Stoyanov, A.I.

2005-01-01

Simple method of radioactive waste first separation using collimated dosimeter, which is placed on boom of power shovel, is proposed, and separation process mathematic modeling for boundary conditions definition of sorting under 'Ukryttya' object high gamma background condition are carry out

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

Coupling nonlinear Stokes and Darcy flow using mortar finite elements

KAUST Repository

Ervin, Vincent J.

2011-11-01

We study a system composed of a nonlinear Stokes flow in one subdomain coupled with a nonlinear porous medium flow in another subdomain. Special attention is paid to the mathematical consequence of the shear-dependent fluid viscosity for the Stokes flow and the velocity-dependent effective viscosity for the Darcy flow. Motivated by the physical setting, we consider the case where only flow rates are specified on the inflow and outflow boundaries in both subdomains. We recast the coupled Stokes-Darcy system as a reduced matching problem on the interface using a mortar space approach. We prove a number of properties of the nonlinear interface operator associated with the reduced problem, which directly yield the existence, uniqueness and regularity of a variational solution to the system. We further propose and analyze a numerical algorithm based on mortar finite elements for the interface problem and conforming finite elements for the subdomain problems. Optimal a priori error estimates are established for the interface and subdomain problems, and a number of compatibility conditions for the finite element spaces used are discussed. Numerical simulations are presented to illustrate the algorithm and to compare two treatments of the defective boundary conditions. © 2010 Published by Elsevier B.V. on behalf of IMACS.

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

CERN Document Server

Bleher, P M

2005-01-01

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

A universal nonlinear relation among boundary states in closed string field theory

International Nuclear Information System (INIS)

Kishimoto, Isao; Matsuo, Yutaka; Watanabe, Eitoku

2004-01-01

We show that the boundary states satisfy a nonlinear relation (the idempotency equation) with respect to the star product of closed string field theory. This relation is universal in the sense that various D-branes, including the infinitesimally deformed ones, satisfy the same equation, including the coefficient. This paper generalizes our analysis [hep-th/0306189] in the following senses. (1) We present a background-independent formulation based on conformal field theory. It illuminates the geometric nature of the relation and allows us to more systematically analyze the variations around the D-brane background. (2) We show that the Witten-type star product satisfies a similar relation but with a more divergent coefficient. (3) We determine the coefficient of the relation analytically. The result shows that the α parameter can be formally factored out, and the relation becomes universal. We present a conjecture on vacuum theory based on this computation. (author)

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

Science.gov (United States)

Zhang, Shuhai; Oskay, Caglar

2015-04-01

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

A new technique for observationally derived boundary conditions for space weather

Science.gov (United States)

Pagano, Paolo; Mackay, Duncan Hendry; Yeates, Anthony Robinson

2018-04-01

Context. In recent years, space weather research has focused on developing modelling techniques to predict the arrival time and properties of coronal mass ejections (CMEs) at the Earth. The aim of this paper is to propose a new modelling technique suitable for the next generation of Space Weather predictive tools that is both efficient and accurate. The aim of the new approach is to provide interplanetary space weather forecasting models with accurate time dependent boundary conditions of erupting magnetic flux ropes in the upper solar corona. Methods: To produce boundary conditions, we couple two different modelling techniques, MHD simulations and a quasi-static non-potential evolution model. Both are applied on a spatial domain that covers the entire solar surface, although they extend over a different radial distance. The non-potential model uses a time series of observed synoptic magnetograms to drive the non-potential quasi-static evolution of the coronal magnetic field. This allows us to follow the formation and loss of equilibrium of magnetic flux ropes. Following this a MHD simulation captures the dynamic evolution of the erupting flux rope, when it is ejected into interplanetary space. Results.The present paper focuses on the MHD simulations that follow the ejection of magnetic flux ropes to 4 R⊙. We first propose a technique for specifying the pre-eruptive plasma properties in the corona. Next, time dependent MHD simulations describe the ejection of two magnetic flux ropes, that produce time dependent boundary conditions for the magnetic field and plasma at 4 R⊙ that in future may be applied to interplanetary space weather prediction models. Conclusions: In the present paper, we show that the dual use of quasi-static non-potential magnetic field simulations and full time dependent MHD simulations can produce realistic inhomogeneous boundary conditions for space weather forecasting tools. Before a fully operational model can be produced there are a

Evolution of passive movement in advective environments: General boundary condition

Science.gov (United States)

Zhou, Peng; Zhao, Xiao-Qiang

2018-03-01

In a previous work [16], Lou et al. studied a Lotka-Volterra competition-diffusion-advection system, where two species are supposed to differ only in their advection rates and the environment is assumed to be spatially homogeneous and closed (no-flux boundary condition), and showed that weaker advective movements are more beneficial for species to win the competition. In this paper, we aim to extend this result to a more general situation, where the environmental heterogeneity is taken into account and the boundary condition at the downstream end becomes very flexible including the standard Dirichlet, Neumann and Robin type conditions as special cases. Our main approaches are to exclude the existence of co-existence (positive) steady state and to provide a clear picture on the stability of semi-trivial steady states, where we introduced new ideas and techniques to overcome the emerging difficulties. Based on these two aspects and the theory of abstract competitive systems, we achieve a complete understanding on the global dynamics.

New Conjugacy Conditions and Related Nonlinear Conjugate Gradient Methods

International Nuclear Information System (INIS)

Dai, Y.-H.; Liao, L.-Z.

2001-01-01

Conjugate gradient methods are a class of important methods for unconstrained optimization, especially when the dimension is large. This paper proposes a new conjugacy condition, which considers an inexact line search scheme but reduces to the old one if the line search is exact. Based on the new conjugacy condition, two nonlinear conjugate gradient methods are constructed. Convergence analysis for the two methods is provided. Our numerical results show that one of the methods is very efficient for the given test problems

Design of HIFU Transducers for Generating Specified Nonlinear Ultrasound Fields.

Science.gov (United States)

Rosnitskiy, Pavel B; Yuldashev, Petr V; Sapozhnikov, Oleg A; Maxwell, Adam D; Kreider, Wayne; Bailey, Michael R; Khokhlova, Vera A

2017-02-01

Various clinical applications of high-intensity focused ultrasound have different requirements for the pressure levels and degree of nonlinear waveform distortion at the focus. The goal of this paper is to determine transducer design parameters that produce either a specified shock amplitude in the focal waveform or specified peak pressures while still maintaining quasi-linear conditions at the focus. Multiparametric nonlinear modeling based on the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation with an equivalent source boundary condition was employed. Peak pressures, shock amplitudes at the focus, and corresponding source outputs were determined for different transducer geometries and levels of nonlinear distortion. The results are presented in terms of the parameters of an equivalent single-element spherically shaped transducer. The accuracy of the method and its applicability to cases of strongly focused transducers were validated by comparing the KZK modeling data with measurements and nonlinear full diffraction simulations for a single-element source and arrays with 7 and 256 elements. The results provide look-up data for evaluating nonlinear distortions at the focus of existing therapeutic systems as well as for guiding the design of new transducers that generate specified nonlinear fields.

Effect of boundary conditions on downstream vorticity from counter-rotating swirlers

Directory of Open Access Journals (Sweden)

Weiye Huo

2015-02-01

Full Text Available Particle image velocimetry (PIV is utilized to measure the non-reacting flow field in a reflow combustor with multiple and single swirlers. The velocity field, vortex structure and total vorticity levels are experimentally obtained using two different boundary conditions, representing a single confined swirler and multiple swirlers in an annular combustor. The influence of the boundary conditions on the flow field at several locations downstream of the swirlers is experimentally investigated, showing that the central vortex in the multi-swirler case is more concentrated than in the single-swirler case. The vorticity of the central vortex and average cross-sectional vorticity are relatively low at the swirler outlet in both cases. Both of these statistics gradually increase to the maximum values near 20 mm downstream of the swirler outlet, and subsequently decrease. It is also found that the central vortex in the multi-swirler case is consistently greater than the single-swirler case. These results demonstrate the critical influence of boundary conditions on flow characteristic of swirling flow, providing insight into the difference of the experiments on test-bed combustor and the full-scale annular combustors.

On solving wave equations on fixed bounded intervals involving Robin boundary conditions with time-dependent coefficients

Science.gov (United States)

van Horssen, Wim T.; Wang, Yandong; Cao, Guohua

2018-06-01

In this paper, it is shown how characteristic coordinates, or equivalently how the well-known formula of d'Alembert, can be used to solve initial-boundary value problems for wave equations on fixed, bounded intervals involving Robin type of boundary conditions with time-dependent coefficients. A Robin boundary condition is a condition that specifies a linear combination of the dependent variable and its first order space-derivative on a boundary of the interval. Analytical methods, such as the method of separation of variables (SOV) or the Laplace transform method, are not applicable to those types of problems. The obtained analytical results by applying the proposed method, are in complete agreement with those obtained by using the numerical, finite difference method. For problems with time-independent coefficients in the Robin boundary condition(s), the results of the proposed method also completely agree with those as for instance obtained by the method of separation of variables, or by the finite difference method.

Aeroelastic Limit-Cycle Oscillations resulting from Aerodynamic Non-Linearities

NARCIS (Netherlands)

van Rooij, A.C.L.M.

2017-01-01

Aerodynamic non-linearities, such as shock waves, boundary layer separation or boundary layer transition, may cause an amplitude limitation of the oscillations induced by the fluid flow around a structure. These aeroelastic limit-cycle oscillations (LCOs) resulting from aerodynamic non-linearities

Global existence and uniform stabilization of a generalized dissipative Klein-Gordon equation type with boundary damping

International Nuclear Information System (INIS)

Zhang Zaiyun; Miao Xiujin; Chen Yuezhong; Liu Zhenhai

2011-01-01

In this paper, we prove the existence, uniqueness, and uniform stability of strong and weak solutions of the nonlinear generalized Klein-Gordon equation (1.1) 1 (see Sec. I) in bounded domains with nonlinear damped boundary conditions given by (1.1) 3 (see Sec. I) with some restrictions on function f(u), h(∇u), g(u t ), and b(x), we prove the existence and uniqueness by means of nonlinear semigroup method and obtain the uniform stabilization by using the multiplier technique.

Growth and wall-transpiration control of nonlinear unsteady Görtler vortices forced by free-stream vortical disturbances

Science.gov (United States)

Marensi, Elena; Ricco, Pierre

2017-11-01

The generation, nonlinear evolution, and wall-transpiration control of unsteady Görtler vortices in an incompressible boundary layer over a concave plate is studied theoretically and numerically. Görtler rolls are initiated and driven by free-stream vortical perturbations of which only the low-frequency components are considered because they penetrate the most into the boundary layer. The formation and development of the disturbances are governed by the nonlinear unsteady boundary-region equations with the centrifugal force included. These equations are subject to appropriate initial and outer boundary conditions, which account for the influence of the upstream and free-stream forcing in a rigorous and mutually consistent manner. Numerical solutions show that the stabilizing effect on nonlinearity, which also occurs in flat-plate boundary layers, is significantly enhanced in the presence of centrifugal forces. Sufficiently downstream, the nonlinear vortices excited at different free-stream turbulence intensities Tu saturate at the same level, proving that the initial amplitude of the forcing becomes unimportant. At low Tu, the disturbance exhibits a quasi-exponential growth with the growth rate being intensified for more curved plates and for lower frequencies. At higher Tu, in the typical range of turbomachinery applications, the Görtler vortices do not undergo a modal stage as nonlinearity saturates rapidly, and the wall curvature does not affect the boundary-layer response. Good quantitative agreement with data from direct numerical simulations and experiments is obtained. Steady spanwise-uniform and spanwise-modulated zero-mass-flow-rate wall transpiration is shown to attenuate the growth of the Görtler vortices significantly. A novel modified version of the Fukagata-Iwamoto-Kasagi identity, used for the first time to study a transitional flow, reveals which terms in the streamwise momentum balance are mostly affected by the wall transpiration, thus

Role of Boundary Conditions in Monte Carlo Simulation of MEMS Devices

Science.gov (United States)

Nance, Robert P.; Hash, David B.; Hassan, H. A.

1997-01-01

A study is made of the issues surrounding prediction of microchannel flows using the direct simulation Monte Carlo method. This investigation includes the introduction and use of new inflow and outflow boundary conditions suitable for subsonic flows. A series of test simulations for a moderate-size microchannel indicates that a high degree of grid under-resolution in the streamwise direction may be tolerated without loss of accuracy. In addition, the results demonstrate the importance of physically correct boundary conditions, as well as possibilities for reducing the time associated with the transient phase of a simulation. These results imply that simulations of longer ducts may be more feasible than previously envisioned.

The effect of external boundary conditions on condensation heat transfer in rotating heat pipes

Science.gov (United States)

Daniels, T. C.; Williams, R. J.

1979-01-01

Experimental evidence shows the importance of external boundary conditions on the overall performance of a rotating heat pipe condenser. Data are presented for the boundary conditions of constant heat flux and constant wall temperature for rotating heat pipes containing either pure vapor or a mixture of vapor and noncondensable gas as working fluid.

Existence of weak solutions to a nonlinear reaction-diffusion system with singular sources

Directory of Open Access Journals (Sweden)

Ida de Bonis

2017-09-01

Full Text Available We discuss the existence of a class of weak solutions to a nonlinear parabolic system of reaction-diffusion type endowed with singular production terms by reaction. The singularity is due to a potential occurrence of quenching localized to the domain boundary. The kind of quenching we have in mind is due to a twofold contribution: (i the choice of boundary conditions, modeling in our case the contact with an infinite reservoir filled with ready-to-react chemicals and (ii the use of a particular nonlinear, non-Lipschitz structure of the reaction kinetics. Our working techniques use fine energy estimates for approximating non-singular problems and uniform control on the set where singularities are localizing.

Models for Predicting Boundary Conditions in L-Mode Tokamak Plasma

International Nuclear Information System (INIS)

Siriwitpreecha, A.; Onjun, T.; Suwanna, S.; Poolyarat, N.; Picha, R.

2009-07-01

Full text: The models for predicting temperature and density of ions and electrons at boundary conditions in L-mode tokamak plasma are developed using an empirical approach and optimized against the experimental data obtained from the latest public version of the International Pedestal Database (version 3.2). It is assumed that the temperature and density at boundary of L-mode plasma are functions of engineering parameters such as plasma current, toroidal magnetic field, total heating power, line averaged density, hydrogenic particle mass (A H ), major radius, minor radius, and elongation at the separatrix. Multiple regression analysis is carried out for these parameters with 86 data points in L-mode from Aug (61) and JT60U (25). The RMSE of temperature and density at boundary of L-mode plasma are found to be 24.41% and 18.81%, respectively. These boundary models are implemented in BALDUR code, which will be used to simulate the L-mode plasma in the tokamak

Self-adjoint elliptic operators with boundary conditions on not closed hypersurfaces

Science.gov (United States)

Mantile, Andrea; Posilicano, Andrea; Sini, Mourad

2016-07-01

The theory of self-adjoint extensions of symmetric operators is used to construct self-adjoint realizations of a second-order elliptic differential operator on Rn with linear boundary conditions on (a relatively open part of) a compact hypersurface. Our approach allows to obtain Kreĭn-like resolvent formulae where the reference operator coincides with the ;free; operator with domain H2 (Rn); this provides an useful tool for the scattering problem from a hypersurface. Concrete examples of this construction are developed in connection with the standard boundary conditions, Dirichlet, Neumann, Robin, δ and δ‧-type, assigned either on a (n - 1) dimensional compact boundary Γ = ∂ Ω or on a relatively open part Σ ⊂ Γ. Schatten-von Neumann estimates for the difference of the powers of resolvents of the free and the perturbed operators are also proven; these give existence and completeness of the wave operators of the associated scattering systems.

Coupling the Gaussian Free Fields with Free and with Zero Boundary Conditions via Common Level Lines

Science.gov (United States)

Qian, Wei; Werner, Wendelin

2018-06-01

We point out a new simple way to couple the Gaussian Free Field (GFF) with free boundary conditions in a two-dimensional domain with the GFF with zero boundary conditions in the same domain: Starting from the latter, one just has to sample at random all the signs of the height gaps on its boundary-touching zero-level lines (these signs are alternating for the zero-boundary GFF) in order to obtain a free boundary GFF. Constructions and couplings of the free boundary GFF and its level lines via soups of reflected Brownian loops and their clusters are also discussed. Such considerations show for instance that in a domain with an axis of symmetry, if one looks at the overlay of a single usual Conformal Loop Ensemble CLE3 with its own symmetric image, one obtains the CLE4-type collection of level lines of a GFF with mixed zero/free boundary conditions in the half-domain.

Involving the Navier-Stokes equations in the derivation of boundary conditions for the lattice Boltzmann method.

Science.gov (United States)

Verschaeve, Joris C G

2011-06-13

By means of the continuity equation of the incompressible Navier-Stokes equations, additional physical arguments for the derivation of a formulation of the no-slip boundary condition for the lattice Boltzmann method for straight walls at rest are obtained. This leads to a boundary condition that is second-order accurate with respect to the grid spacing and conserves mass. In addition, the boundary condition is stable for relaxation frequencies close to two.

Probability distribution of magnetization in the one-dimensional Ising model: effects of boundary conditions

Energy Technology Data Exchange (ETDEWEB)

Antal, T [Physics Department, Simon Fraser University, Burnaby, BC V5A 1S6 (Canada); Droz, M [Departement de Physique Theorique, Universite de Geneve, CH 1211 Geneva 4 (Switzerland); Racz, Z [Institute for Theoretical Physics, Eoetvoes University, 1117 Budapest, Pazmany setany 1/a (Hungary)

2004-02-06

Finite-size scaling functions are investigated both for the mean-square magnetization fluctuations and for the probability distribution of the magnetization in the one-dimensional Ising model. The scaling functions are evaluated in the limit of the temperature going to zero (T {yields} 0), the size of the system going to infinity (N {yields} {infinity}) while N[1 - tanh(J/k{sub B}T)] is kept finite (J being the nearest neighbour coupling). Exact calculations using various boundary conditions (periodic, antiperiodic, free, block) demonstrate explicitly how the scaling functions depend on the boundary conditions. We also show that the block (small part of a large system) magnetization distribution results are identical to those obtained for free boundary conditions.

Threshold condition for nonlinear tearing modes in tokamaks

International Nuclear Information System (INIS)

Zabiego, M.F.; Callen, J.D.

1996-03-01

Low-mode-number tearing, mode nonlinear evolution is analyzed emphasizing the need for a threshold condition, to account for observations in tokamaks. The discussion is illustrated by two models recently introduced in the literature. The models can be compared with the available data and/or serve as a basis for planning some experiments in order to either test theory (by means of beta-limit scaling laws, as proposed in this paper) or attempt to control undesirable tearing modes. Introducing a threshold condition in the tearing mode stability analysis is found to reveal some bifurcation points and thus domains of intrinsic stability in the island dynamics operational space

Generalized nematohydrodynamic boundary conditions with application to bistable twisted nematic liquid-crystal displays

KAUST Repository

Fang, Angbo

2008-12-08

Parallel to the highly successful Ericksen-Leslie hydrodynamic theory for the bulk behavior of nematic liquid crystals (NLCs), we derive a set of coupled hydrodynamic boundary conditions to describe the NLC dynamics near NLC-solid interfaces. In our boundary conditions, translational flux (flow slippage) and rotational flux (surface director relaxation) are coupled according to the Onsager variational principle of least energy dissipation. The application of our boundary conditions to the truly bistable π -twist NLC cell reveals a complete picture of the dynamic switching processes. It is found that the thus far overlooked translation-rotation dissipative coupling at solid surfaces can accelerate surface director relaxation and enhance the flow rate. This can be utilized to improve the performance of electro-optical nematic devices by lowering the required switching voltages and reducing the switching times. © 2008 The American Physical Society.

Impact of helical boundary conditions on nonlinear 3D magnetohydrodynamic simulations of reversed-field pinch

International Nuclear Information System (INIS)

Veranda, M; Bonfiglio, D; Cappello, S; Chacón, L; Escande, D F

2013-01-01

Helical self-organized reversed-field pinch (RFP) regimes emerge both numerically—in 3D visco-resistive magnetohydrodynamic (MHD) simulations—and experimentally, as in the RFX-mod device at high current (I P above 1 MA). These states, called quasi-single helicity (QSH) states, are characterized by the action of a MHD mode that impresses a quasi-helical symmetry to the system, thus allowing a high degree of magnetic chaos healing. This is in contrast with the multiple helicity (MH) states, where magnetic fluctuations create a chaotic magnetic field degrading the confinement properties of the RFP. This paper reports an extensive numerical study performed in the frame of 3D visco-resistive MHD which considers the effect of helical magnetic boundary conditions, i.e. of a finite value of the radial magnetic field at the edge (magnetic perturbation, MP). We show that the system can be driven to a selected QSH state starting from both spontaneous QSH and MH regimes. In particular, a high enough MP can force a QSH helical self-organization with a helicity different from the spontaneous one. Moreover, MH states can be turned into QSH states with a selected helicity. A threshold in the amplitude of MP is observed above which is able to influence the system. Analysis of the magnetic topology of these simulations indicates that the dominant helical mode is able to temporarily sustain conserved magnetic structures in the core of the plasma. The region occupied by conserved magnetic surfaces increases reducing secondary modes' amplitude to experimental-like values. (paper)

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

KAUST Repository

Leise, Tanya L.; Walton, Jay R.; Gorb, Yuliya

2009-01-01

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

A quenched study of the Schroedinger functional with chirally rotated boundary conditions. Non-preturbative tuning

Energy Technology Data Exchange (ETDEWEB)

Lopez, J. Gonzalez [Humboldt-Universitaet, Berlin (Germany). Inst. fuer Physik; Deutsches Elektronen-Synchrotron (DESY), Zeuthen (Germany); Jansen, K. [Deutsches Elektronen-Synchrotron (DESY), Zeuthen (Germany); Renner, D.B. [Jefferson Lab, Newport News, VA (United States); Shindler, A. [Humboldt-Universitaet, Berlin (Germany). Inst. fuer Physik

2012-08-23

The use of chirally rotated boundary conditions provides a formulation of the Schroedinger functional that is compatible with automatic O(a) improvement of Wilson fermions up to O(a) boundary contributions. The elimination of bulk O(a) effects requires the non-perturbative tuning of the critical mass and one additional boundary counterterm. We present the results of such a tuning in a quenched setup for several values of the renormalized gauge coupling, from perturbative to nonperturbative regimes, and for a range of lattice spacings. We also check that the correct boundary conditions and symmetries are restored in the continuum limit. (orig.)

A quenched study of the Schroedinger functional with chirally rotated boundary conditions. Non-preturbative tuning

International Nuclear Information System (INIS)

Lopez, J. Gonzalez; Jansen, K.; Renner, D.B.; Shindler, A.

2012-01-01

The use of chirally rotated boundary conditions provides a formulation of the Schroedinger functional that is compatible with automatic O(a) improvement of Wilson fermions up to O(a) boundary contributions. The elimination of bulk O(a) effects requires the non-perturbative tuning of the critical mass and one additional boundary counterterm. We present the results of such a tuning in a quenched setup for several values of the renormalized gauge coupling, from perturbative to nonperturbative regimes, and for a range of lattice spacings. We also check that the correct boundary conditions and symmetries are restored in the continuum limit. (orig.)

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

International Nuclear Information System (INIS)

Vasileva, D.P.